SEARCH · Search PubMed
Results for “Genome alignment”
Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Analysis of the transcriptional complexity of Arabidopsis thaliana by massively parallel signature sequencing.
Large-scale sequencing of short mRNA-derived tags can establish the qualitative and quantitative characteristics of a complex transcriptome. We sequenced 12,304,362 tags from five diverse libraries of Arabidopsis thaliana using massively parallel signature sequencing (MPSS). A total of 48,572 distinct signatures, each representing a different transcript, were expressed at significant levels. These signatures were compared to the annotation of the A. thaliana genomic sequence; in the five libraries, this comparison yielded between 17,353 and 18,361 genes with sense expression, and between 5,487 and 8,729 genes with antisense expression. An additional 6,691 MPSS signatures mapped to unannotated regions of the genome. Expression was demonstrated for 1,168 genes for which expression data were previously unknown. Alternative polyadenylation was observed for more than 25% of A. thaliana genes transcribed in these libraries. The MPSS expression data suggest that the A. thaliana transcriptome is complex and contains many as-yet uncharacterized variants of normal coding transcripts.
Herpesviruses throw a curve ball: new insights into microRNA biogenesis and evolution.
Explore the source record for details and available documents.
Identification of microRNAs of the herpesvirus family.
Epstein-Barr virus (EBV or HHV4), a member of the human herpesvirus (HHV) family, has recently been shown to encode microRNAs (miRNAs). In contrast to most eukaryotic miRNAs, these viral miRNAs do not have close homologs in other viral genomes or in the genome of the human host. To identify other miRNA genes in pathogenic viruses, we combined a new miRNA gene prediction method with small-RNA cloning from several virus-infected cell types. We cloned ten miRNAs in the Kaposi sarcoma-associated virus (KSHV or HHV8), nine miRNAs in the mouse gammaherpesvirus 68 (MHV68) and nine miRNAs in the human cytomegalovirus (HCMV or HHV5). These miRNA genes are expressed individually or in clusters from either polymerase (pol) II or pol III promoters, and share no substantial sequence homology with one another or with the known human miRNAs. Generally, we predicted miRNAs in several large DNA viruses, and we could neither predict nor experimentally identify miRNAs in the genomes of small RNA viruses or retroviruses.
resLens: genomic language models to enhance antibiotic resistance gene detection.
The rise of antibiotic resistance necessitates advanced tools to detect and analyze antibiotic resistance genes (ARGs). We present resLens, a family of genomic language models that leverage latent genomic representations to enhance ARG detection and analysis. Unlike alignment-based methods constrained by reference databases, resLens fine-tunes a pre-trained DNA language model on curated ARG datasets, achieving competitive or superior performance in classifying resistance genes across multiple evaluation scenarios, including when ARGs exhibit sequences and mechanisms of resistance dissimilar to those in reference datasets.
The origin of the eukaryotic cell: a genomic investigation.
We have collected a set of 347 proteins that are found in eukaryotic cells but have no significant homology to proteins in Archaea and Bacteria. We call these proteins eukaryotic signature proteins (ESPs). The dominant hypothesis for the formation of the eukaryotic cell is that it is a fusion of an archaeon with a bacterium. If this hypothesis is accepted then the three cellular domains, Eukarya, Archaea, and Bacteria, would collapse into two cellular domains. We have used the existence of this set of ESPs to test this hypothesis. The evidence of the ESPs implicates a third cell (chronocyte) in the formation of the eukaryotic cell. The chronocyte had a cytoskeleton that enabled it to engulf prokaryotic cells and a complex internal membrane system where lipids and proteins were synthesized. It also had a complex internal signaling system involving calcium ions, calmodulin, inositol phosphates, ubiquitin, cyclin, and GTP-binding proteins. The nucleus was formed when a number of archaea and bacteria were engulfed by a chronocyte. This formation of the nucleus would restore the three cellular domains as the Chronocyte was not a cell that belonged to the Archaea or to the Bacteria.
Cloning, characterization and prokaryotic expression of cytosolic malate dehydrogenase from Oryza sativa.
cDNA sequences of malate dehydrogenase (MDH) were cloned from various species, and MDH was identified to play an important role in cell energy metabolism. Here, we present the isolation and characterization of its homologue (OscMDH) in Oryza sativa. Comparison of the results to the genome details indicated that OscMDH consisted of seven exons. Sequence alignment showed that the deduced amino acid sequence of OscMDH shared a significant similarity with cMDH protein in Zea mays, as well as with other cMDH gene products. The different expression patterns of OscMDH in different tissues revealed that OscMDH mRNA was highly transcribed in either young panicle or immature seed, while its abundance was much low in green leaf and root. A nearly 56-kDa fusion protein generated by adding a Trx-His-tag at the N-terminal of OscMDH was induced by IPTG in Escherichia coli BL21 and an obvious MDH activity was detected in the protein by native PAGE analysis. All these results suggest that OscMDH encodes a cytosolic MDH in rice.
A comparison of sequence and structure protein domain families as a basis for structural genomics.
MOTIVATION: Protein families can be defined based on structure or sequence similarity. We wanted to compare two protein family databases, one based on structural and one on sequence similarity, to investigate to what extent they overlap, the similarity in definition of corresponding families, and to create a list of large protein families with unknown structure as a resource for structural genomics. We also wanted to increase the sensitivity of fold assignment by exploiting protein family HMMs. RESULTS: We compared Pfam, a protein family database based on sequence similarity, to Scop, which is based on structural similarity. We found that 70% of the Scop families exist in Pfam while 57% of the Pfam families exist in Scop. Most families that occur in both databases correspond well to each other, but in some cases they are different. Such cases highlight situations in which structure and sequence approaches differ significantly. The comparison enabled us to compile a list of the largest families that do not occur in Scop; these are suitable targets for structure prediction and determination, and may be useful to guide projects in structural genomics. It can be noted that 13 out of the 20 largest protein families without a known structure are likely transmembrane proteins. We also exploited Pfam to increase the sensitivity of detecting homologs of proteins with known structure, by comparing query sequences to Pfam HMMs that correspond to Scop families. For SWISSPROT+TREMBL, this yielded an increase in fold assignment from 31% to 42% compared to using FASTA only. This method assigned a structure to 22% of the proteins in Saccharomyces cerevisiae, 24% in Escherichia coli, and 16% in Methanococcus jannaschii.
Analysis of genomic sequences by Chaos Game Representation.
MOTIVATION: Chaos Game Representation (CGR) is an iterative mapping technique that processes sequences of units, such as nucleotides in a DNA sequence or amino acids in a protein, in order to find the coordinates for their position in a continuous space. This distribution of positions has two properties: it is unique, and the source sequence can be recovered from the coordinates such that distance between positions measures similarity between the corresponding sequences. The possibility of using the latter property to identify succession schemes have been entirely overlooked in previous studies which raises the possibility that CGR may be upgraded from a mere representation technique to a sequence modeling tool. RESULTS: The distribution of positions in the CGR plane were shown to be a generalization of Markov chain probability tables that accommodates non-integer orders. Therefore, Markov models are particular cases of CGR models rather than the reverse, as currently accepted. In addition, the CGR generalization has both practical (computational efficiency) and fundamental (scale independence) advantages. These results are illustrated by using Escherichia coli K-12 as a test data-set, in particular, the genes thrA, thrB and thrC of the threonine operon.
me-PCR: a refined ultrafast algorithm for identifying sequence-defined genomic elements.
We have adapted the originally described electronic PCR (e-PCR) algorithm to perform string searches more accurately and much more rapidly than previously possible. Our implementation [multithreaded e-PCR (me-PCR)] runs sufficiently fast to allow even desktop machines to query quickly large genomes with very large genomic element sets. In addition, me-PCR is multithreaded, interprets all IUPAC nucleotide symbols, allows searches with elements specified by long sequences (such as SNPs), accepts ranges in the expected PCR size input field, requires substantially less memory for analysis of large sequences and corrects a number of minor flaws causing misreporting of hits in exceptional cases. Thus, me-PCR provides increased annotation capabilities for complex genomes to non-expert laboratories.
GeneSyn: a tool for detecting conserved gene order across genomes.
UNLABELLED: GeneSyn is a software tool that allows automatic detection of conserved gene order from annotated genomes. AVAILABILITY: Available free of charge for Unix/Linux/Cygwin platforms at ftp://159.149.110.11/pub/GeneSyn_1.0/ SUPPLEMENTARY INFORMATION: ftp://159.149.110.11/pub/GeneSyn_1.0/
Common structural patterns in human genes.
MOTIVATION: A graphical representation of the exon-intron structure of various genes, such as that presented by the National Center for Biotechnology Information Map Viewer, suggests a digital waveform or pattern that varies either in amplitude or frequency. This observation suggests that different genes may have portions of their total exon-intron structure in common. The existence of common structural patterns across unrelated genes suggests the repeated insertion of transposable elements throughout the human genome and/or a common structural function. RESULTS: We compared the exon-intron size patterns of a number of human genes and discovered numerous conserved arrangements with similarity at a high degree of stringency (>99%) across the otherwise unrelated and diverse genomic landscape. In our experimental analyses, more than 200 patterns of length 2 or greater at 99% stringency were found among the 72 genes we compared.
A practical and robust sequence search strategy for structural genomics target selection.
MOTIVATION: Target selection strategies for structural genomic projects must be able to prioritize gene regions on the basis of significant sequence similarity with proteins that have already been structurally determined. With the rapid development of protein comparison software a robust prioritization scheme should be independent of the choice of algorithm and be able to incorporate different sequence similarity thresholds. RESULTS: A robust target selection strategy has been developed that can assign a priority level to all genes in any genome. Structural assignments to genome sequences are calculated at two thresholds and six levels (1-6) describe the prioritization of all whole genes and partial gene regions. This simple two-threshold approach can be implemented with any fold recognition or homology detection algorithms. The results for 10 genomes are presented using the SSEARCH and PSI-BLAST programs. AVAILABILITY: Programs are available on request from the authors.
Describing ancient horizontal gene transfers at the nucleotide and gene levels by comparative pathogenicity island genometrics.
MOTIVATION: Lateral gene transfer is a major mechanism contributing to bacterial genome dynamics and pathovar emergence via pathogenicity island (PAI) spreading. However, since few of these genomic exchanges are experimentally reproducible, it is difficult to establish evolutionary scenarios for the successive PAI transmissions between bacterial genera. Methods initially developed at the gene and/or nucleotide level for genomics, i.e. comparisons of concatenated sequences, ortholog frequency, gene order or dinucleotide usage, were combined and applied here to homologous PAIs: we call this approach comparative PAI genometrics. RESULTS: YAPI, a Yersinia PAI, and related islands were compared with measure evolutionary relationships between related modules. Through use of our genometric approach designed for tracking codon usage adaptation and gene phylogeny, an ancient inter-genus PAI transfer was oriented for the first time by characterizing the genomic environment in which the ancestral island emerged and its subsequent transfers to other bacterial genera.
Hypoxia Response Is Associated with Reduced HPV Activity and Tumor Microenvironment Remodeling in Cervical Cancer.
Human papillomavirus (HPV) significantly influences cervical cancer progression and treatment, yet its interactions with the tumor microenvironment remain incompletely understood. We performed single-cell and spatial transcriptomic sequencing on cervical cancer samples to explore these interactions. By aligning sequencing reads to a merged HPV16-human reference genome, we characterized HPV16 heterogeneity and its association with host states at the single-cell and single-gene levels. E5 transcriptional activity was negatively associated with the host interferon response, indicating a role in immune evasion. A hypoxic environment was correlated with the downregulation of E5 activity and elevated MHC-I expression, which may contribute to stronger interactions between hypoxic cancer cells and cytotoxic CD8⁺ T cells. Additionally, HPV16 integration in host cells was associated with increased fatty acid metabolism. These findings suggest that combining anti-angiogenic drugs and fatty acid metabolism inhibitors has the potential to improve cervical cancer treatment.
Glutamine:phenylpyruvate aminotransferase from an extremely thermophilic bacterium, Thermus thermophilus HB8.
A subfamily I aminotransferase gene homologue containing an open reading frame encoding 381 amino acid residues (Mr=42,271) has been identified in the process of the genome project of an extremely thermophilic bacterium, Thermus thermophilus HB8. Alignment of the predicted amino acid sequence using FASTA shows that this protein is a member of aminotransferase subfamily Igamma. The protein shows around 40% identity with both T. thermophilus aspartate aminotransferase [EC 2.6.1.1] and mammalian glutamine:phenylpyruvate aminotransferase [EC 2.6.1.64]. The recombinant protein expressed in Escherichia coli is a homodimer with a subunit molecular weight of 42,000, has one pyridoxal 5'-phosphate per subunit, and is highly active toward glutamine, methionine, aromatic amino acids, and corresponding keto acids, but has no preference for alanine and dicarboxylic amino acids. These substrate specificities are similar to those described for mammalian glutamine: phenylpyruvate aminotransferase. This is the first enzyme reported so far that has the glutamine aminotransferase activity in non-eukaryotic cells. As the presence of aromatic amino acid:2-oxoglutarate aminotransferase [EC 2.6.1.57] has not been reported in T. thermophilus, this enzyme is expected to catalyze the last transamination step of phenylalanine and tyrosine biosynthesis. It may also be involved in the methionine regeneration pathway associated with polyamine biosynthesis. The enzyme shows a strikingly high pKa value (9.3) of the coenzyme Schiff base in comparison with other subfamily I aminotransferases. The origin of this unique pKa value and the substrate specificity is discussed based on the previous crystallographic data of T. thermophilus and E. coli aspartate aminotransferases.
The DNA replication origins of herpes simplex virus type 1 strain Angelotti.
The nucleotide sequences of the origins of DNA replication (ori) of the S- and L-component (oriS, oriL) of the herpes simplex virus type 1 (HSV-1) standard genome were determined from HSV-1 strain Angelotti (ANG). In contrast to other HSV-1 strains, the ANG oriS sequence revealed an insertion of an TA-dinucleotide in an otherwise very conserved but imperfect palindromic sequence of 47 bp. The oriL sequence of the standard ANG genome was found to be identical to that of an ANG class II defective genome which exhibits a duplication of a 144 bp palindrome. A model is presented to explain the origination of the amplified ANG oriL sequences from the parental genome with a single copy of oriL via illegitimate recombination. Alignment of the ori sequences of HSV, adeno- and papovaviruses unveiled that the HSV ori region can be subdivided into two distinct sites of homology to the DNA initiation signals of papova- and adenoviruses, suggesting that the HSV origins of replication comprise elements for DNA replication by both, cellular and virus-encoded DNA polymerases.
Compilation of tRNA sequences and sequences of tRNA genes.
Maintained at the Universitat Bayreuth, Bayreuth, Germany, the Compilation of tRNA Sequences and Sequences of tRNA Genes is accessible at the URL http://www.tRNA.uni-bayreuth.de with mirror site located at the Institute of Protein Research, Pushchino, Russia (http://alpha.protres.ru/trnadbase). The compilation is a searchable, periodically updated database of currently available tRNA sequences. The present version of the database contains a new Genomic tRNA Compilation including the sequences of tRNA genes from genomic sequences published up to July 2003. It consists of about 5800 tRNA gene sequences from 111 organisms covering archaea, bacteria, higher and lower eukarya. The former Compilation of tRNA Genes (up to the end of 1998) and the updated Compilation tRNA Sequences (561 entries) are also supported by the new software. The database can be explored by using multiple search criteria and sequence templates. The database provides a service that allows to obtain statistical information on the occurrences of certain bases at given positions of the tRNA sequences. This allows phylogenic studies and search for identity elements in respect to interactions of tRNAs with various enzymes.