Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Molecular Sequence Annotation”

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.

At least 937 records · Page 52Linked to original sources

Comprehensive proteomic analysis of breast cancer cell membranes reveals unique proteins with potential roles in clinical cancer.

Proteins associated with cancer cell plasma membranes are rich in known drug and antibody targets as well as other proteins known to play key roles in the abnormal signal transduction processes required for carcinogenesis. We describe here a proteomics process that comprehensively annotates the protein content of breast tumor cell membranes and defines the clinical relevance of such proteins. Tumor-derived cell lines were used to ensure an enrichment for cancer cell-specific plasma membrane proteins because it is difficult to purify cancer cells and then obtain good membrane preparations from clinical material. Multiple cell lines with different molecular pathologies were used to represent the clinical heterogeneity of breast cancer. Peptide tandem mass spectra were searched against a comprehensive data base containing known and conceptual proteins derived from many public data bases including the draft human genome sequences. This plasma membrane-enriched proteome analysis created a data base of more than 500 breast cancer cell line proteins, 27% of which were of unknown function. The value of our approach is demonstrated by further detailed analyses of three previously uncharacterized proteins whose clinical relevance has been defined by their unique cancer expression profiles and the identification of protein-binding partners that elucidate potential functionality in cancer.

Amino Acid Sequence↗

The SWISS-MODEL Repository: new features and functionalities.

The SWISS-MODEL Repository is a database of annotated 3D protein structure models generated by the SWISS-MODEL homology-modelling pipeline. As of September 2005, the repository contained 675,000 models for 604,000 different protein sequences of the UniProt database. Regular updates ensure that the content of the repository reflects the current state of sequence and structure databases, integrating new or modified target sequences, and making use of new template structures. Each Repository entry consists of one or more 3D models accompanied by detailed information about the target protein and the model building process: functional annotation, a detailed template selection log, target-template alignment, summary of the model building and model quality assessment. The SWISS-MODEL Repository is freely accessible at http://swissmodel.expasy.org/repository/.

Computer Graphics↗

Discovery of eight novel divergent homologs expressed in cattle placenta.

Ten divergent homologs were identified using a subtractive bioinformatic analysis of 12,614 cattle placenta expressed sequence tags followed by comparative, evolutionary, and gene expression studies. Among the 10 divergent homologs, 8 have not been identified previously. These were named as follows: cattle cerebrum and skeletal muscle-specific transcript 1 (CSSMST1), cattle intestine-specific transcript 1 (CIST1), hepatitis A virus cellular receptor 1 amino-terminal domain-containing protein (HAVCRNDP), prolactin-related proteins 8, 9, and 11 (PRP8, PRP9, and PRP11, respectively) and secreted and transmembrane protein 1A and 1B (SECTM1A and SECTM1B, respectively). In addition, two previously known divergent genes were identified, trophoblast Kunitz domain protein 1 (TKDP1) and a new splice variant of TKDP4. Nucleotide substitution analysis provided evidence for positive selection in members of the PRP gene family, SECTM1A and SECTM1B. Gene expression profiles, motif predictions, and annotations of homologous sequences indicate immunological and reproductive functions of the divergent homologs. The genes identified in this study are thus of evolutionary and physiological importance and may have a role in placental adaptations.

Amino Acid Sequence↗

Phylogenomic analysis of the GIY-YIG nuclease superfamily.

BACKGROUND: The GIY-YIG domain was initially identified in homing endonucleases and later in other selfish mobile genetic elements (including restriction enzymes and non-LTR retrotransposons) and in enzymes involved in DNA repair and recombination. However, to date no systematic search for novel members of the GIY-YIG superfamily or comparative analysis of these enzymes has been reported. RESULTS: We carried out database searches to identify all members of known GIY-YIG nuclease families. Multiple sequence alignments together with predicted secondary structures of identified families were represented as Hidden Markov Models (HMM) and compared by the HHsearch method to the uncharacterized protein families gathered in the COG, KOG, and PFAM databases. This analysis allowed for extending the GIY-YIG superfamily to include members of COG3680 and a number of proteins not classified in COGs and to predict that these proteins may function as nucleases, potentially involved in DNA recombination and/or repair. Finally, all old and new members of the GIY-YIG superfamily were compared and analyzed to infer the phylogenetic tree. CONCLUSION: An evolutionary classification of the GIY-YIG superfamily is presented for the very first time, along with the structural annotation of all (sub)families. It provides a comprehensive picture of sequence-structure-function relationships in this superfamily of nucleases, which will help to design experiments to study the mechanism of action of known members (especially the uncharacterized ones) and will facilitate the prediction of function for the newly discovered ones.

Archaeal Proteins↗

IMGT databases, web resources and tools for immunoglobulin and T cell receptor sequence analysis, http://imgt.cines.fr.

IMGT, the international ImMunoGeneTics database((R)) (http://imgt.cines.fr), is a high-quality integrated information system specializing in immunoglobulins (IG), T cell receptors (TR) and major histocompatibility complex (MHC) of human and other vertebrates, created in 1989, by LIGM, at the Université Montpellier II, CNRS, Montpellier, France. IMGT provides a common access to standardized data which include nucleotide and protein sequences, oligonucleotide primers, gene maps, genetic polymorphisms, specificities, 2D and 3D structures. IMGT includes several databases (IMGT/LIGM-DB, IMGT/3Dstructure-DB, IMGT/HLA-DB), Web resources ('IMGT Marie-Paule page') and interactive tools (IMGT/V-QUEST, IMGT/JunctionAnalysis). IMGT expertly annotated data and tools described in this paper are particularly useful for the analysis of the IG and TR rearrangements in leukemia, lymphoma and myeloma, and in translocations involving the antigen receptor loci. IMGT is freely available at http://imgt.cines.fr.

Amino Acid Sequence↗

Quantitative analysis of nucleic acid three-dimensional structures.

A new computer program to annotate DNA and RNA three-dimensional structures, MC-Annotate, is introduced. The goals of annotation are to efficiently extract and manipulate structural information, to simplify further structural analyses and searches, and to objectively represent structural knowledge. The input of MC-Annotate is a PDB formatted DNA or RNA three-dimensional structure. The output of MC-Annotate is composed of a structural graph that contains the annotations, and a series of HTML documents, one for each nucleotide conformation and base-base interaction present in the input structure. The atomic coordinates of all nucleotides and the homogeneous transformation matrices of all base-base interactions are stored in the structural graph. Symbolic classifications of nucleotide conformations, using sugar puckering modes and nitrogen base orientations around the glycosyl bond, and base-base interactions, using stacking and hydrogen bonding information, are introduced. Peculiarity factors of nucleotide conformations and base-base interactions are defined to indicate their marginalities with all other examples. The peculiarity factors allow us to identify irregular regions and possible stereochemical errors in 3-D structures without interactive visualization. The annotations attached to each nucleotide conformation include its class, its torsion angles, a distribution of the root-mean-square deviations with examples of the same class, the list of examples of the same class, and its peculiarity value. The annotations attached to each base-base interaction include its class, a distribution of distances with examples of the same class, the list of examples of the same class, and its peculiarity value. The distance between two homogeneous transformation matrices is evaluated using a new metric that distinguishes between the rotation and the translation of a transformation matrix in the context of nitrogen bases. MC-Annotate was used to build databases of nucleotide conformations and base-base interactions. It was applied to the ribosomal RNA fragment that binds to protein L11, which annotations revealed peculiar nucleotide conformations and base-base interactions in the regions where the RNA contacts the protein. The question of whether the current database of RNA three-dimensional structures is complete is addressed.

Base Pairing↗

Improving genome annotations using phylogenetic profile anomaly detection.

MOTIVATION: A promising strategy for refining genome annotations is to detect features that conflict with known functional or evolutionary relationships between groups of genes. Previous work in this area has been focused on investigating the absence of 'housekeeping' genes or components of well-studied pathways. We have sought to develop a method for improving new annotations that can automatically synthesize and use the information available in a database of other annotated genomes. RESULTS: We show that a probabilistic model of phylogenetic profiles, trained from a database of curated genome annotations, can be used to reliably detect errors in new annotations. We use our method to identify 22 genes that were missed in previously published annotations of prokaryotic genomes. AVAILABILITY: The method was evaluated using MATLAB and open source software referenced in this work. Scripts and datasets are available from the authors upon request. CONTACT: tarjei@broad.mit.edu.

Algorithms↗

UTRdb and UTRsite: specialized databases of sequences and functional elements of 5' and 3' untranslated regions of eukaryotic mRNAs.

The 5' and 3' untranslated regions of eukaryotic mRNAs may play a crucial role in the regulation of gene expression controlling mRNA localization, stability and translational efficiency. For this reason we developed UTRdb, a specialized database of 5' and 3' untranslated sequences of eukaryotic mRNAs cleaned from redundancy. UTRdb entries are enriched with specialized information not present in the primary databases including the presence of nucleotide sequence patterns already demonstrated by experimental analysis to have some functional role. All these patterns have been collected in the UTRsite database so that it is possible to search any input sequence for the presence of annotated functional motifs. Furthermore, UTRdb entries have been annotated for the presence of repetitive elements. All internet resources implemented for retrieval and functional analysis of 5' and 3' untranslated regions of eukaryotic mRNAs are accessible at http://bigarea.area.ba.cnr.it:8000/EmbIT/UTRH ome/

3' Untranslated Regions↗

Re-annotating the Mycoplasma pneumoniae genome sequence: adding value, function and reading frames.

Four years after the original sequence submission, we have re-annotated the genome of Mycoplasma pneumoniae to incorporate novel data. The total number of ORFss has been increased from 677 to 688 (10 new proteins were predicted in intergenic regions, two further were newly identified by mass spectrometry and one protein ORF was dismissed) and the number of RNAs from 39 to 42 genes. For 19 of the now 35 tRNAs and for six other functional RNAs the exact genome positions were re-annotated and two new tRNA(Leu) and a small 200 nt RNA were identified. Sixteen protein reading frames were extended and eight shortened. For each ORF a consistent annotation vocabulary has been introduced. Annotation reasoning, annotation categories and comparisons to other published data on M.pneumoniae functional assignments are given. Experimental evidence includes 2-dimensional gel electrophoresis in combination with mass spectrometry as well as gene expression data from this study. Compared to the original annotation, we increased the number of proteins with predicted functional features from 349 to 458. The increase includes 36 new predictions and 73 protein assignments confirmed by the published literature. Furthermore, there are 23 reductions and 30 additions with respect to the previous annotation. mRNA expression data support transcription of 184 of the functionally unassigned reading frames.

Amino Acid Sequence↗

The relationship between protein sequences and their gene ontology functions.

BACKGROUND: One main research challenge in the post-genomic era is to understand the relationship between protein sequences and their biological functions. In recent years, several automated annotation systems have been developed for the functional assignment of uncharacterized proteins. The underlying assumption of these systems is that similar sequences imply similar biological functions. However, it has been noted that matching sequences do not always infer similar functions. RESULTS: In this paper, we present the correlation between protein sequences and protein functions for the yeast proteome in the context of gene ontology. A novel measure is introduced to define the overall similarity between two protein sequences. The effects of the level as well as the size of a gene ontology group on the degree of similarity were studied. The similarity distributions at different levels of gene ontology trees are presented. To evaluate the theoretical prediction power of similar sequences, we computed the posterior probability of correct predictions. CONCLUSION: The results indicate that protein pairs of similar biological functions tend to have higher sequence similarity, although the similarity distribution in each functional group is heterogeneous and varies from group to group. We conclude that sequence similarity can serve as a key measure in protein function prediction. However, the resulting annotations must be verified through other means. A method that combines a broader range of measures is more likely to provide more accurate prediction. Our study indicates that the posterior probability of a correct prediction could serve as one of the key measures.

Amino Acid Sequence↗

PhyloGibbs: a Gibbs sampling motif finder that incorporates phylogeny.

A central problem in the bioinformatics of gene regulation is to find the binding sites for regulatory proteins. One of the most promising approaches toward identifying these short and fuzzy sequence patterns is the comparative analysis of orthologous intergenic regions of related species. This analysis is complicated by various factors. First, one needs to take the phylogenetic relationship between the species into account in order to distinguish conservation that is due to the occurrence of functional sites from spurious conservation that is due to evolutionary proximity. Second, one has to deal with the complexities of multiple alignments of orthologous intergenic regions, and one has to consider the possibility that functional sites may occur outside of conserved segments. Here we present a new motif sampling algorithm, PhyloGibbs, that runs on arbitrary collections of multiple local sequence alignments of orthologous sequences. The algorithm searches over all ways in which an arbitrary number of binding sites for an arbitrary number of transcription factors (TFs) can be assigned to the multiple sequence alignments. These binding site configurations are scored by a Bayesian probabilistic model that treats aligned sequences by a model for the evolution of binding sites and "background" intergenic DNA. This model takes the phylogenetic relationship between the species in the alignment explicitly into account. The algorithm uses simulated annealing and Monte Carlo Markov-chain sampling to rigorously assign posterior probabilities to all the binding sites that it reports. In tests on synthetic data and real data from five Saccharomyces species our algorithm performs significantly better than four other motif-finding algorithms, including algorithms that also take phylogeny into account. Our results also show that, in contrast to the other algorithms, PhyloGibbs can make realistic estimates of the reliability of its predictions. Our tests suggest that, running on the five-species multiple alignment of a single gene's upstream region, PhyloGibbs on average recovers over 50% of all binding sites in S. cerevisiae at a specificity of about 50%, and 33% of all binding sites at a specificity of about 85%. We also tested PhyloGibbs on collections of multiple alignments of intergenic regions that were recently annotated, based on ChIP-on-chip data, to contain binding sites for the same TF. We compared PhyloGibbs's results with the previous analysis of these data using six other motif-finding algorithms. For 16 of 21 TFs for which all other motif-finding methods failed to find a significant motif, PhyloGibbs did recover a motif that matches the literature consensus. In 11 cases where there was disagreement in the results we compiled lists of known target genes from the literature, and found that running PhyloGibbs on their regulatory regions yielded a binding motif matching the literature consensus in all but one of the cases. Interestingly, these literature gene lists had little overlap with the targets annotated based on the ChIP-on-chip data. The PhyloGibbs code can be downloaded from http://www.biozentrum.unibas.ch/~nimwegen/cgi-bin/phylogibbs.cgi or http://www.imsc.res.in/~rsidd/phylogibbs. The full set of predicted sites from our tests on yeast are available at http://www.swissregulon.unibas.ch.

Algorithms↗

Species-specific transcriptomic changes upon respiratory syncytial virus infection in cotton rats.

The cotton rat (Sigmodon) is the gold standard pre-clinical small animal model for respiratory viral pathogens, especially for respiratory syncytial virus (RSV). However, without a reference genome or a published transcriptome, studies requiring gene expression analysis in cotton rats are severely limited. The aims of this study were to generate a comprehensive transcriptome from multiple tissues of two species of cotton rats that are commonly used as animal models (Sigmodon fulviventer and Sigmodon hispidus), and to compare and contrast gene expression changes and immune responses to RSV infection between the two species. Transcriptomes were assembled from lung, spleen, kidney, heart, and intestines for each species with a contig N50 > 1600. Annotation of contigs generated nearly 120,000 gene annotations for each species. The transcriptomes of S. fulviventer and S. hispidus were then used to assess immune response to RSV infection. We identified 238 unique genes that are significantly differentially expressed, including several genes implicated in RSV infection (e.g., Mx2, I27L2, LY6E, Viperin, Keratin 6A, ISG15, CXCL10, CXCL11, IRF9) as well as novel genes that have not previously described in RSV research (LG3BP, SYWC, ABEC1, IIGP1, CREB1). This study presents two comprehensive transcriptome references as resources for future gene expression analysis studies in the cotton rat model, as well as provides gene sequences for mechanistic characterization of molecular pathways. Overall, our results provide generalizable insights into the effect of host genetics on host-virus interactions, as well as identify new host therapeutic targets for RSV treatment and prevention.

Animals↗

Protein dispensability and rate of evolution.

If protein evolution is due in large part to slightly deleterious amino acid substitutions, then the rate of evolution should be greater in proteins that contribute less to individual fitness. The rationale for this prediction is that relatively dispensable proteins should be subject to weaker purifying selection, and should therefore accumulate mildly deleterious substitutions more rapidly. Although this argument was presented over twenty years ago, and is fundamental to many applications of evolutionary theory, the prediction has proved difficult to confirm. In fact, a recent study showed that essential mouse genes do not evolve more slowly than non-essential ones. Thus, although a variety of factors influencing the rate of protein evolution have been supported by extensive sequence analysis, the relationship between protein dispensability and evolutionary rate has remained unconfirmed. Here we use the results from a highly parallel growth assay of single gene deletions in yeast to assess protein dispensability, which we relate to evolutionary rate estimates that are based on comparisons of sequences drawn from twenty-one fully annotated genomes. Our analysis reveals a highly significant relationship between protein dispensability and evolutionary rate, and explains why this relationship is not detectable by categorical comparison of essential versus non-essential proteins. The relationship is highly conserved, so that protein dispensability in yeast is also predictive of evolutionary rate in a nematode worm.

Amino Acid Substitution↗

Gene expression profiling of bovine in vitro adipogenesis using a cDNA microarray.

The gene expression profile of bovine bone marrow stromal cells undergoing adipogenesis was established using a custom cDNA microarray. Cells that were treated with adipogenic stimulants and those that were not were collected at each of the six time points, and gene expression differences between the treated and untreated samples within each time point were compared using a microarray. Statistical analyses revealed that 158 genes showed a minimum fold change of 2 in at least one of the five post-differentiation time points. These genes are involved in various cellular pathways and functions, including lipogenesis, glycolysis, cytoskeleton remodelling, extracellular matrix, transcription as well as various signalling pathways such as insulin, calcium and wingless signalling. The experiment also identified 17 differentially expressed (DE) microarray elements with no assigned function. Quantitative real-time PCR was employed to validate eight DE genes, and the PCR data were found to reproduce the microarray data for these eight genes. Subsequent gene ontology annotation was able to provide a global overview of the molecular function of DE genes during adipogenesis. This analysis was able to indicate the importance of different gene categories at various stages of adipogenic conversion, thereby providing further insights into the molecular changes during bovine adipogenesis.

Adipogenesis↗

Genetic and molecular analysis of region 88E9;88F2 in Drosophila melanogaster, including the ear gene related to human factors involved in lineage-specific leukemias.

We identified and characterized the Drosophila gene ear (ENL/AF9-related), which is closely related to mammalian genes that have been implicated in the onset of acute lymphoblastic and myelogenous leukemias when their products are fused as chimeras with those of human HRX, a homolog of Drosophila trithorax. The ear gene product is present in all early embryonic cells, but becomes restricted to specific tissues in late embryogenesis. We mapped the ear gene to cytological region 88E11-13, near easter, and showed that it is deleted by Df(3R)ea(5022rx1), a small, cytologically invisible deletion. Annotation of the completed Drosophila genome sequence suggests that this region might contain as many as 26 genes, most of which, including ear, are not represented by mutant alleles. We carried out a large-scale noncomplementation screen using Df(3R)ea(5022rx1) and chemical (EMS) mutagenesis from which we identified seven novel multi-allele recessive lethal complementation groups in this region. An overlapping deficiency, Df(3R)Po(4), allowed us to map several of these groups to either the proximal or the distal regions of Df(3R)ea(5022rx1). One of these complementation groups likely corresponds to the ear gene as judged by map location, terminal phenotype, and reduction of EAR protein levels.

Animals↗

Molecular evolution and population genetics of duplicated accessory gland protein genes in Drosophila.

To investigate the potential importance of gene duplication in D. melanogaster accessory gland protein (Acp) gene evolution we carried out a computational analysis comparing annotated D. melanogaster Acp genes to the entire D. melanogaster genome. We found that two known Acp genes are actually members of small multigene families. Polymorphism and divergence data from these duplicated genes suggest that in at least four cases, protein divergence between D. melanogaster and D. simulans is a result of directional selection. One putative Acp revealed by our computational analysis shows evidence of a recent selective sweep in a non-African population (but not in an African population). These data support the idea that selection on reproduction-related genes may drive divergence of populations within species, and strengthen the conclusion that Acps may often be under directional selection in Drosophila.

Amino Acid Sequence↗

Orthologs, paralogs, and evolutionary genomics.

Orthologs and paralogs are two fundamentally different types of homologous genes that evolved, respectively, by vertical descent from a single ancestral gene and by duplication. Orthology and paralogy are key concepts of evolutionary genomics. A clear distinction between orthologs and paralogs is critical for the construction of a robust evolutionary classification of genes and reliable functional annotation of newly sequenced genomes. Genome comparisons show that orthologous relationships with genes from taxonomically distant species can be established for the majority of the genes from each sequenced genome. This review examines in depth the definitions and subtypes of orthologs and paralogs, outlines the principal methodological approaches employed for identification of orthology and paralogy, and considers evolutionary and functional implications of these concepts.

Animals↗

Mapping of immunoreactive antigens of Francisella tularensis live vaccine strain.

Francisella tularensis subsp. holarctica is the common causal agent of tularemia in Europe. Besides clinical signs, the diagnosis of the disease mostly depends on serological tests. To date, there is a lack of information about the F. tularensis antigens that induce antibody response. Therefore, we have started comprehensive mapping of immunoreactive antigens using the attenuated live vaccine strain of F. tularensis LVS originating from the European virulent strain. For this purpose, the immunoreactivity of sera collected from patients suffering from tularemia, together with the control sera of patients with Lyme disease and healthy blood donors, were examined by means of one-dimensional and two-dimensional immunoblotting. Furthermore, whole cell bacterial lysates, isolated integral membrane proteins and basic proteins were exploited as antigens. By this approach more than 80 different immunorelevant antigens were detected. Most of them came from whole cell bacterial lysate and integral membrane proteins. Conversely, only a negligible reaction was found in the case of basic proteins. Forty-five spots were further selected for mass spectrometric analyses and 22 of them were annotated. Among the spots that provided characteristic reactions with sera from patients with tularemia, 60 kDa and 10 kDa chaperonins that occurred in several charge and mass variants, predominated.

Animals↗