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 739 records · Page 41Linked to original sources

The cyclope gene of Drosophila encodes a cytochrome c oxidase subunit VIc homolog.

Cytochrome c oxidase is the terminal enzyme of the mitochondrial electron transfer chain. In eukaryotes, the enzyme is composed of 3 mitochondrial DNA-encoded subunits and 7-10 (in mammals) nuclear DNA-encoded subunits. This enzyme has been extensively studied in mammals and yeast but, in Drosophila, very little is known and no mutant has been described so far. Here we report the genetic and molecular characterization of mutations in cyclope (cype) and the cloning of the gene encoding a cytochrome c oxidase subunit VIc homolog. cype is an essential gene whose mutations are lethal and show pleiotropic phenotypes. The 77-amino acid peptide encoded by cype is 46% identical and 59% similar to the human subunit (75 amino acids). The transcripts are expressed maternally and throughout development in localized regions. They are found predominantly in the central nervous system of the embryo; in the central region of imaginal discs; in the germarium, follicular, and nurse cells of the ovary; and in testis. A search in the Genome Annotation Database of Drosophila revealed the absence of subunit VIIb and the presence of 9 putative nuclear cytochrome c oxidase subunits with high identity scores when compared to the 10 human subunits.

Alleles↗

Transposable element-driven expansion of enhancer RNA repertoires underlies regulatory innovation and polyploid adaptation in cereal crops.

Cereal genomes have undergone repeated polyploidization and transposable element (TE) proliferation, collectively generating complex regulatory landscapes. However, the evolutionary trajectories and functional implications of these landscapes remain largely unexplored. Using chromatin-bound RNA sequencing across seven cereal species, we systematically mapped 45,952 regulatory element transcripts (RETs), including 32,867 distal RETs corresponding to enhancer RNAs (eRNAs). Our analysis revealed that 56% of lineage-specific eRNAs originated from TE expansions, indicating that TEs serve as major reservoirs of species-specific regulatory innovation in cereals. Notably, we identified remarkable conservation in defense-related functions, root-specific expression, and TE-derived origins of eRNAs across both ancient and recent evolutionary layers of Triticeae, suggesting recurrent recruitment of TE-derived, root-associated regulatory elements throughout Triticeae evolution. Furthermore, we found that young eRNA pairs in hexaploid wheat with high sequence similarity, many originating from RLG_famc8.3 and DTC_famc4.3, exhibited pronounced root specificity and coordinated expression, suggesting targeted amplification and refinement of successful ancestral regulatory strategies established after Triticeae divergence. To facilitate community access, we developed Cereal-eRNAdb (http://bioinfo.cemps.ac.cn/Cereal-eRNAdb/), a comprehensive database integrating 69,426 eRNAs with functional annotations across 296 samples. Our findings suggest that TE-mediated innovation of root-specific eRNAs may contribute to Triticeae adaptation and provide a foundational resource for exploiting regulatory variation in cereal crop breeding.

Enhancer RNAs↗

COCO-CL: hierarchical clustering of homology relations based on evolutionary correlations.

MOTIVATION: Determining orthology relations among genes across multiple genomes is an important problem in the post-genomic era. Identifying orthologous genes can not only help predict functional annotations for newly sequenced or poorly characterized genomes, but can also help predict new protein-protein interactions. Unfortunately, determining orthology relation through computational methods is not straightforward due to the presence of paralogs. Traditional approaches have relied on pairwise sequence comparisons to construct graphs, which were then partitioned into putative clusters of orthologous groups. These methods do not attempt to preserve the non-transitivity and hierarchic nature of the orthology relation. RESULTS: We propose a new method, COCO-CL, for hierarchical clustering of homology relations and identification of orthologous groups of genes. Unlike previous approaches, which are based on pairwise sequence comparisons, our method explores the correlation of evolutionary histories of individual genes in a more global context. COCO-CL can be used as a semi-independent method to delineate the orthology/paralogy relation for a refined set of homologous proteins obtained using a less-conservative clustering approach, or as a refiner that removes putative out-paralogs from clusters computed using a more inclusive approach. We analyze our clustering results manually, with support from literature and functional annotations. Since our orthology determination procedure does not employ a species tree to infer duplication events, it can be used in situations when the species tree is unknown or uncertain. CONTACT: jothi@mail.nih.gov, przytyck@mail.nih.gov SUPPLEMENTARY INFORMATION: Supplementary materials are available at Bioinformatics online.

Algorithms↗

Modeling a whole organ using proteomics: the avian bursa of Fabricius.

While advances in proteomics have improved proteome coverage and enhanced biological modeling, modeling function in multicellular organisms requires understanding how cells interact. Here we used the chicken bursa of Fabricius, a common experimental system for B cell function, to model organ function from proteomics data. The bursa has two major functional cell types: B cells and the supporting stromal cells. We used differential detergent fractionation-multidimensional protein identification technology (DDF-MudPIT) to identify 5198 proteins from all cellular compartments. Of these, 1753 were B cell specific, 1972 were stroma specific and 1473 were shared between the two. By modeling programmed cell death (PCD), cell differentiation and proliferation, and transcriptional activation, we have improved functional annotation of chicken proteins and placed chicken-specific death receptors into the PCD process using phylogenetics. We have identified 114 transcription factors (TFs); 42 of the bursal B cell TFs have not been reported before in any B cells. We have also improved the structural annotation of a newly sequenced genome by confirming the in vivo expression of 4006 "predicted", and 6623 ab initio, ORFs. Finally, we have developed a novel method for facilitating structural annotation, "expressed peptide sequence tags" (ePSTs) and demonstrate its utility by identifying 521 potential novel proteins from the chicken "unassigned chromosome".

Amino Acid Sequence↗

Identification of the mammalian Not gene via a phylogenomic approach.

Despite the great morphological diversity of early embryos, the underlying mechanisms of gastrulation are known to be broadly conserved in vertebrates. However, a number of genes characterized as fulfilling an essential function in this process in several model organisms display no clear ortholog in mammalian genomes. We have devised an in silico phylogenomic approach, based on exhaustive similarity searches in vertebrate genomes and subsequent bayesian phylogenetic analyses, to identify such missing genes, presumed to be highly divergent. This approach has been used to identify mammalian orthologs of Not, an homeodomain containing gene previously characterized in Xenopus, chick and zebrafish as playing a critical role in the formation of the notochord. This attempt led to the identification of a highly divergent mammalian Not-related gene in the mouse, human and rat. The results from phylogenetic reconstructions, synteny analyses, expression pattern analyses in wild-type and mutant mouse embryos, and overexpression experiments in Xenopus embryos converge to confirm these genes as representatives of the Not family in mammals. The identification of the mammalian Not gene delivers an important component for the understanding of the genetics underlying notochord formation in mammals and its evolution among vertebrates. The phylogenomic method used to retrieve this gene thus provides a tool, which can complement or validate genome annotations in situations when they are weakly supported.

Amino Acid Sequence↗

An evolutionary approach reveals a high protein-coding capacity of the human genome.

We developed a new evolutionary method for identifying exons from genomic sequences and found 19000 potential coding exons that are absent from all existing annotations of the human genome. Of these, 13700 satisfied very stringent criteria and can with confidence be considered as novel exons. Evidently, a large number of new human genes can be identified using evolutionary approaches.

Animals↗

Functional analysis of a group A streptococcal glycoside hydrolase Spy1600 from family 84 reveals it is a beta-N-acetylglucosaminidase and not a hyaluronidase.

Group A streptococcus (Streptococcus pyogenes) is the causative agent of severe invasive infections such as necrotizing fasciitis (the so-called 'flesh eating disease') and toxic-shock syndrome. Spy1600, a glycoside hydrolase from family 84 of the large superfamily of glycoside hydrolases, has been proposed to be a virulence factor. In the present study we show that Spy1600 has no activity toward galactosaminides or hyaluronan, but does remove beta-O-linked N-acetylglucosamine from mammalian glycoproteins--an observation consistent with the inclusion of eukaryotic O-glycoprotein 2-acetamido-2-deoxy-beta-D-glucopyranosidases within glycoside hydrolase family 84. Proton NMR studies, structure-reactivity studies for a series of fluorinated analogues and analysis of 1,2-dideoxy-2'-methyl-alpha-D-glucopyranoso-[2,1-d]-Delta2'-thiazoline as a competitive inhibitor reveals that Spy1600 uses a double-displacement mechanism involving substrate-assisted catalysis. Family 84 glycoside hydrolases are therefore comprised of both prokaryotic and eukaryotic beta-N-acetylglucosaminidases using a conserved catalytic mechanism involving substrate-assisted catalysis. Since these enzymes do not have detectable hyaluronidase activity, many family 84 glycoside hydrolases are most likely incorrectly annotated as hyaluronidases.

Acetylglucosamine↗

SynBrowse: a synteny browser for comparative sequence analysis.

MOTIVATION: The recent efforts of various sequence projects to sequence deeply into various phylogenies provide great resources for comparative sequence analysis. A generic and portable tool is essential for scientists to visualize and analyze sequence comparisons. RESULTS: We have developed SynBrowse, a synteny browser for visualizing and analyzing genome alignments both within and between species. It is intended to help scientists study macrosynteny, microsynteny and homologous genes between sequences. It can also aid with the identification of uncharacterized genes, putative regulatory elements and novel structural features of a species. SynBrowse is a GBrowse (the Generic Genome Browser) family software tool that runs on top of the open source BioPerl modules. It consists of two components: a web-based front end and a set of relational database back ends. Each database stores pre-computed alignments from a focus sequence to reference sequences in addition to the genome annotations of the focus sequence. The user interface lets end users select a key comparative alignment type and search for syntenic blocks between two sequences and zoom in to view the relationships among the corresponding genome annotations in detail. SynBrowse is portable with simple installation, flexible configuration, convenient data input and easy integration with other components of a model organism system. AVAILABILITY: The software is available at http://www.gmod.org CONTACT: vbrendel@iastate.edu

Algorithms↗

A hybrid gene team model and its application to genome analysis.

It is well-known that functionally related genes occur in a physically clustered form, especially operons in bacteria. By leveraging on this fact, there has recently been an interesting problem formulation known as gene team model, which searches for a set of genes that co-occur in a pair of closely related genomes. However, many gene teams, even experimentally verified operons, frequently scatter within other genomes. Thus, the gene team model should be refined to reflect this observation. In this paper, we generalized the gene team model, that looks for gene clusters in a physically clustered form, to multiple genome cases with relaxed constraints. We propose a novel hybrid pattern model that combines the set and the sequential pattern models. Our model searches for gene clusters with and/or without physical proximity constraint. This model is implemented and tested with 97 genomes (120 replicons). The result was analyzed to show the usefulness of our model. We also compared the result from our hybrid model to those from the traditional gene team model. We also show that predicted gene teams can be used for various genome analysis: operon prediction, phylogenetic analysis of organisms, contextual sequence analysis and genome annotation. Our program is fast enough to provide a service on the web at http://platcom.informatics.indiana.edu/platcom/. Users can select any combination of 97 genomes to predict gene teams.

Algorithms↗

Heterologous expression of proteins from Plasmodium falciparum: results from 1000 genes.

As part of a structural genomics initiative, 1000 open reading frames from Plasmodium falciparum, the causative agent of the most deadly form of malaria, were tested in an E. coli protein expression system. Three hundred and thirty-seven of these targets were observed to express, although typically the protein was insoluble. Sixty-three of the targets provided soluble protein in yields ranging from 0.9 to 406.6 mg from one liter of rich media. Higher molecular weight, greater protein disorder (segmental analysis, SEG), more basic isoelectric point (pI), and a lack of homology to E. coli proteins were all highly and independently correlated with difficulties in expression. Surprisingly, codon usage and the percentage of adenosines and thymidines (%AT) did not appear to play a significant role. Of those proteins which expressed, high pI and a hypothetical annotation were both strongly and independently correlated with insolubility. The overwhelmingly important role of pI in both expression and solubility appears to be a surprising and fundamental issue in the heterologous expression of P. falciparum proteins in E. coli. Twelve targets which did not express in E. coli from the native gene sequence were codon-optimized through whole gene synthesis, resulting in the (insoluble) expression of three of these proteins. Seventeen targets which were expressed insolubly in E. coli were moved into a baculovirus/Sf-21 system, resulting in the soluble expression of one protein at a high level and six others at a low level. A variety of factors conspire to make the heterologous expression of P. falciparum proteins challenging, and these observations lay the groundwork for a rational approach to prioritizing and, ultimately, eliminating these impediments.

Animals↗

Complexities in ETS-domain transcription factor function and regulation: lessons from the TCF (ternary complex factor) subfamily. The Colworth Medal Lecture.

The ETS-domain transcription factor family can be divided into a series of subfamilies. Elk-1 represents the founding member of the ternary complex factor (TCF) subfamily. By focusing on the TCF subfamily, we can demonstrate the complexities that exist in the function and regulation of ETS-domain transcription factors. This article focuses on Elk-1 in detail and summarizes the functions of other TCFs. The key themes covered include the domain structure of the TCFs, the mechanisms of complex formation with serum response factor, regulation of TCFs by mitogen-activated protein kinase cascades, and transcriptional regulatory properties of the TCFs. Finally, the emerging role of the TCFs in vivo is discussed. A picture is developing indicating that, while these proteins exhibit significant sequence and functional conservation, key differences in their structure and regulation are being identified which may relate to unique functions of these proteins in vivo.

Amino Acid Sequence↗

DNASTAR's Lasergene sequence analysis software.

Lasergene's eight modules provide tools that enable users to accomplish each step of sequence analysis, from trimming and assembly of sequence data, to gene discovery, annotation, gene product analysis, sequence similarity searches, sequence alignment, phylogenetic analysis, oligonucleotide primer design, cloning strategies, and publication of the results. The Lasergene software suite provides the functions and customization tools needed so that users can perform analyses the software writers never imagined.

Base Sequence↗

Molecular Characterization of Listeria monocytogenes Isolated from Retail Yak Meat in Nyingchi, Xizang, China.

Listeria monocytogenes is a Gram-positive zoonotic pathogen responsible for listeriosis, a severe foodborne disease with high mortality in humans and animals. This study aimed to investigate the molecular epidemiology and genomic characteristics of L. monocytogenes isolated from raw yak meat in Nyingchi, Xizang, China. A total of 231 yak-related samples were collected in Nyingchi, consisting of 214 retail raw yak meat samples, 14 farm environmental samples, and 3 nearby water source samples. L. monocytogenes isolates were identified and characterized using culture-based methods, PCR serotyping, and whole-genome sequencing (WGS). Bioinformatic analyses were performed for virulence, antimicrobial resistance, and functional gene annotation using KEGG and COG databases. The overall contamination rate of Lm was 13.08% (28/214) for retail raw yak meat samples, whereas no isolates were recovered from 14 farm environmental samples (0.00%, 0/14) and 3 nearby water source samples (0.00%, 0/3). The serotypes of isolates were 1/2a (9/28, 32.14%), 1/2b (7/28, 25.00%), and 1/2c (12/28, 42.86%). These 28 isolates exhibited varied antimicrobial resistance profiles, with universal resistance to trimethoprim-sulfamethoxazole, high resistance to erythromycin and clindamycin, and low resistance to vancomycin. MLST analysis revealed seven sequence types (STs): ST9 (12/28, 42.86%), ST619 (6/28, 21.43%), ST8 (6/28, 21.43%), ST7 (1/28, 3.57%), ST87 (1/28, 3.57%), ST121 (1/28, 3.57%), ST297 (1/28, 3.57%). ST619 isolates harbored multiple virulence genes, including those located on Listeria pathogenicity islands LIPI-1, LIPI-3, and LIPI-4, indicating high genomic potential for virulence. Representative isolate Y2 (ST619) possessed a 3,009,858 bp genome with 3036 coding genes, four genomic islands, and two prophages. Functional annotation revealed enrichment of genes involved in carbohydrate transport and metabolism and amino acid biosynthesis pathways. Our findings provide the first genomic insight into L. monocytogenes contamination in yak meat from Nyingchi, Xizang, China, highlighting the urgent need to strengthen food safety monitoring and hygiene management in this region.

Listeria monocytogenes↗

SMASHing regulatory sites in DNA by human-mouse sequence comparisons.

Regulatory sequence elements provide important clues to understanding and predicting gene expression. Although the binding sites for hundreds of transcription factors are known, there has been no systematic attempt to incorporate this information in the annotation of the human genome. Cross species sequence comparisons are critical to a meaningful annotation of regulatory elements since they generally reside in conserved non-coding regions. To take advantage of the recently completed drafts of the mouse and human genomes for annotating transcription factor binding sites, we developed SMASH, a computational pipeline that identifies thousands of orthologous human/ mouse proteins, maps them to genomic sequences, extracts and compares upstream regions and annotates putative regulatory elements in conserved, non-coding, upstream regions. Our current dataset consists of approximately 2,500 human/mouse gene pairs. Transcription start sites were estimated by mapping quasi-full length cDNA sequences. SMASH uses a novel probabilistic method to identify putative conserved binding sites that takes into account the competition between transcription factors for binding DNA. SMASH presents the results via a genome browser web interface which displays the predicted regulatory information together with the current annotations for the human genome. Our results are validated by comparison to previously published experimental data. SMASH results compare favorably to other existing computational approaches.

Algorithms↗

Syntenic relationships between Medicago truncatula and Arabidopsis reveal extensive divergence of genome organization.

Arabidopsis and Medicago truncatula represent sister clades within the dicot subclass Rosidae. We used genetic map-based and bacterial artificial chromosome sequence-based approaches to estimate the level of synteny between the genomes of these model plant species. Mapping of 82 tentative orthologous gene pairs reveals a lack of extended macrosynteny between the two genomes, although marker collinearity is frequently observed over small genetic intervals. Divergence estimates based on non-synonymous nucleotide substitutions suggest that a majority of the genes under analysis have experienced duplication in Arabidopsis subsequent to divergence of the two genomes, potentially confounding synteny analysis. Moreover, in cases of localized synteny, genetically linked loci in M. truncatula often share multiple points of synteny with Arabidopsis; this latter observation is consistent with the large number of segmental duplications that compose the Arabidopsis genome. More detailed analysis, based on complete sequencing and annotation of three M. truncatula bacterial artificial chromosome contigs suggests that the two genomes are related by networks of microsynteny that are often highly degenerate. In some cases, the erosion of microsynteny could be ascribed to the selective gene loss from duplicated loci, whereas in other cases, it is due to the absence of close homologs of M. truncatula genes in Arabidopsis.

Arabidopsis↗

Using hidden Markov models and observed evolution to annotate viral genomes.

MOTIVATION: ssRNA (single stranded) viral genomes are generally constrained in length and utilize overlapping reading frames to maximally exploit the coding potential within the genome length restrictions. This overlapping coding phenomenon leads to complex evolutionary constraints operating on the genome. In regions which code for more than one protein, silent mutations in one reading frame generally have a protein coding effect in another. To maximize coding flexibility in all reading frames, overlapping regions are often compositionally biased towards amino acids which are 6-fold degenerate with respect to the 64 codon alphabet. Previous methodologies have used this fact in an ad hoc manner to look for overlapping genes by motif matching. In this paper differentiated nucleotide compositional patterns in overlapping regions are incorporated into a probabilistic hidden Markov model (HMM) framework which is used to annotate ssRNA viral genomes. This work focuses on single sequence annotation and applies an HMM framework to ssRNA viral annotation. A description of how the HMM is parameterized, whilst annotating within a missing data framework is given. A Phylogenetic HMM (Phylo-HMM) extension, as applied to 14 aligned HIV2 sequences is also presented. This evolutionary extension serves as an illustration of the potential of the Phylo-HMM framework for ssRNA viral genomic annotation. RESULTS: The single sequence annotation procedure (SSA) is applied to 14 different strains of the HIV2 virus. Further results on alternative ssRNA viral genomes are presented to illustrate more generally the performance of the method. The results of the SSA method are encouraging however there is still room for improvement, and since there is overwhelming evidence to indicate that comparative methods can improve coding sequence (CDS) annotation, the SSA method is extended to a Phylo-HMM to incorporate evolutionary information. The Phylo-HMM extension is applied to the same set of 14 HIV2 sequences which are pre-aligned. The performance improvement that results from including the evolutionary information in the analysis is illustrated.

Algorithms↗

Identification of promoter regions in the human genome by using a retroviral plasmid library-based functional reporter gene assay.

Attempts to identify regulatory sequences in the human genome have involved experimental and computational methods such as cross-species sequence comparisons and the detection of transcription factor binding-site motifs in coexpressed genes. Although these strategies provide information on which genomic regions are likely to be involved in gene regulation, they do not give information on their functions. We have developed a functional selection for promoter regions in the human genome that uses a retroviral plasmid library-based system. This approach enriches for and detects promoter function of isolated DNA fragments in an in vitro cell culture assay. By using this method, we have discovered likely promoters of known and predicted genes, as well as many other putative promoter regions based on the presence of features such as CpG islands. Comparison of sequences of 858 plasmid clones selected by this assay with the human genome draft sequence indicates that a significantly higher percentage of sequences align to the 500-bp segment upstream of the transcription start sites of known genes than would be expected from random genomic sequences. We also observed enrichment for putative promoter regions of genes predicted in at least two annotation databases and for clones overlapping with CpG islands. Functional validation of randomly selected clones enriched by this method showed that a large fraction of these putative promoters can drive the expression of a reporter gene in transient transfection experiments. This method promises to be a useful genome-wide function-based approach that can complement existing methods to look for promoters.

3T3 Cells↗

Establishment of the epithelial-specific transcriptome of normal and malignant human breast cells based on MPSS and array expression data.

INTRODUCTION: Diverse microarray and sequencing technologies have been widely used to characterise the molecular changes in malignant epithelial cells in breast cancers. Such gene expression studies to identify markers and targets in tumour cells are, however, compromised by the cellular heterogeneity of solid breast tumours and by the lack of appropriate counterparts representing normal breast epithelial cells. METHODS: Malignant neoplastic epithelial cells from primary breast cancers and luminal and myoepithelial cells isolated from normal human breast tissue were isolated by immunomagnetic separation methods. Pools of RNA from highly enriched preparations of these cell types were subjected to expression profiling using massively parallel signature sequencing (MPSS) and four different genome wide microarray platforms. Functional related transcripts of the differential tumour epithelial transcriptome were used for gene set enrichment analysis to identify enrichment of luminal and myoepithelial type genes. Clinical pathological validation of a small number of genes was performed on tissue microarrays. RESULTS: MPSS identified 6,553 differentially expressed genes between the pool of normal luminal cells and that of primary tumours substantially enriched for epithelial cells, of which 98% were represented and 60% were confirmed by microarray profiling. Significant expression level changes between these two samples detected only by microarray technology were shown by 4,149 transcripts, resulting in a combined differential tumour epithelial transcriptome of 8,051 genes. Microarray gene signatures identified a comprehensive list of 907 and 955 transcripts whose expression differed between luminal epithelial cells and myoepithelial cells, respectively. Functional annotation and gene set enrichment analysis highlighted a group of genes related to skeletal development that were associated with the myoepithelial/basal cells and upregulated in the tumour sample. One of the most highly overexpressed genes in this category, that encoding periostin, was analysed immunohistochemically on breast cancer tissue microarrays and its expression in neoplastic cells correlated with poor outcome in a cohort of poor prognosis estrogen receptor-positive tumours. CONCLUSION: Using highly enriched cell populations in combination with multiplatform gene expression profiling studies, a comprehensive analysis of molecular changes between the normal and malignant breast tissue was established. This study provides a basis for the identification of novel and potentially important targets for diagnosis, prognosis and therapy in breast cancer.

Biomarkers, Tumor↗