PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Molecular Sequence Annotation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70Linked to original sources

Complete nucleotide sequence and genome analysis of bacteriophage BFK20--a lytic phage of the industrial producer Brevibacterium flavum.

The entire double-stranded DNA genome of bacteriophage BFK20, a lytic phage of the Brevibacterium flavum CCM 251--industrial producer of L-lysine--was sequenced and analyzed. It consists of 42,968 base pairs with an overall molar G + C content of 56.2%. Fifty-five potential open reading frames were identified and annotated using various bioinformatics tools. Clusters of functionally related putative genes were defined (structural, lytic, replication and regulatory). To verify the annotation of structural proteins, they were resolved by 2D gel electrophoresis and were submitted to N-terminal amino acid sequencing. Structural proteins identified included the portal and major and minor tail proteins. Based on the overall genome sequence comparison, similarities with other known bacteriophage genomes include primarily bacteriophages from Mycobacterium spp. and some regions of Corynebacterium spp. genomes--possible prophages. Our results support the theory that phage genomes are mosaics with respect to each other.

Bacteriophages↗

Deciphering the ghost proteome in ovarian cancer cells by deep proteogenomic characterization.

Proteogenomics is becoming a powerful tool in personalized medicine by linking genomics, transcriptomics and mass spectrometry (MS)-based proteomics. Due to increasing evidence of alternative open reading frame-encoded proteins (AltProts), proteogenomics has a high potential to unravel the characteristics, variants, expression levels of the alternative proteome, in addition to already annotated proteins (RefProts). To obtain a broader view of the proteome of ovarian cancer cells compared to ovarian epithelial cells, cell-specific total RNA-sequencing profiles and customized protein databases were generated. In total, 128 RefProts and 30 AltProts were identified exclusively in SKOV-3 and PEO-4 cells. Among them, an AltProt variant of IP_715944, translated from DHX8, was found mutated (p.Leu44Pro). We show high variation in protein expression levels of RefProts and AltProts in different subcellular compartments. The presence of 117 RefProt and two AltProt variants was described, along with their possible implications in the different physiological/pathological characteristics. To identify the possible involvement of AltProts in cellular processes, cross-linking-MS (XL-MS) was performed in each cell line to identify AltProt-RefProt interactions. This approach revealed an interaction between POLD3 and the AltProt IP_183088, which after molecular docking, was placed between POLD3-POLD2 binding sites, highlighting its possibility of the involvement in DNA replication and repair.

Humans↗

An integrated gene annotation and transcriptional profiling approach towards the full gene content of the Drosophila genome.

BACKGROUND: While the genome sequences for a variety of organisms are now available, the precise number of the genes encoded is still a matter of debate. For the human genome several stringent annotation approaches have resulted in the same number of potential genes, but a careful comparison revealed only limited overlap. This indicates that only the combination of different computational prediction methods and experimental evaluation of such in silico data will provide more complete genome annotations. In order to get a more complete gene content of the Drosophila melanogaster genome, we based our new D. melanogaster whole-transcriptome microarray, the Heidelberg FlyArray, on the combination of the Berkeley Drosophila Genome Project (BDGP) annotation and a novel ab initio gene prediction of lower stringency using the Fgenesh software. RESULTS: Here we provide evidence for the transcription of approximately 2,600 additional genes predicted by Fgenesh. Validation of the developmental profiling data by RT-PCR and in situ hybridization indicates a lower limit of 2,000 novel annotations, thus substantially raising the number of genes that make a fly. CONCLUSIONS: The successful design and application of this novel Drosophila microarray on the basis of our integrated in silico/wet biology approach confirms our expectation that in silico approaches alone will always tend to be incomplete. The identification of at least 2,000 novel genes highlights the importance of gathering experimental evidence to discover all genes within a genome. Moreover, as such an approach is independent of homology criteria, it will allow the discovery of novel genes unrelated to known protein families or those that have not been strictly conserved between species.

Animals↗

Applications of generalized pair hidden Markov models to alignment and gene finding problems.

Hidden Markov models (HMMs) have been successfully applied to a variety of problems in molecular biology, ranging from alignment problems to gene finding and annotation. Alignment problems can be solved with pair HMMs, while gene finding programs rely on generalized HMMs in order to model exon lengths. In this paper, we introduce the generalized pair HMM (GPHMM), which is an extension of both pair and generalized HMMs. We show how GPHMMs, in conjunction with approximate alignments, can be used for cross-species gene finding and describe applications to DNA-cDNA and DNA-protein alignment. GPHMMs provide a unifying and probabilistically sound theory for modeling these problems.

Algorithms↗

Heuristic approach to deriving models for gene finding.

Computer methods of accurate gene finding in DNA sequences require models of protein coding and non-coding regions derived either from experimentally validated training sets or from large amounts of anonymous DNA sequence. Here we propose a new, heuristic method producing fairly accurate inhomogeneous Markov models of protein coding regions. The new method needs such a small amount of DNA sequence data that the model can be built 'on the fly' by a web server for any DNA sequence >400 nt. Tests on 10 complete bacterial genomes performed with the GeneMark.hmm program demonstrated the ability of the new models to detect 93.1% of annotated genes on average, while models built by traditional training predict an average of 93.9% of genes. Models built by the heuristic approach could be used to find genes in small fragments of anonymous prokaryotic genomes and in genomes of organelles, viruses, phages and plasmids, as well as in highly inhomogeneous genomes where adjustment of models to local DNA composition is needed. The heuristic method also gives an insight into the mechanism of codon usage pattern evolution.

Codon↗

MHCDB--database of the human MHC.

Genetic and physical data relating to the human major histocompatibility complex (MHC) were compiled and analyzed, using a genome analysis program. The current contents of the database include: 1) location of over 100 genes and other markers; 2) location of over 250 YAC and cosmid clones; 3) 150 kilobases of genomic DNA sequence including full annotation (exon/intron boundaries, repeats, promoters, etc.); 4) cDNA sequences of currently-known class I and class II alleles; and 5) accompanying descriptive data--references, comments, laboratory addresses, and so on. The database allows rapid access, retrieval, and display of all these data, which should make it a useful tool for the study of the human MHC. MHCDB is publicly available and will be updated as new data become available.

Alleles↗

Characteristics of the Lotus japonicus gene repertoire deduced from large-scale expressed sequence tag (EST) analysis.

To perform a comprehensive analysis of genes expressed in a model legume, Lotus japonicus, a total of 74472 3'-end expressed sequence tags (EST) were generated from cDNA libraries produced from six different organs. Clustering of sequences was performed with an identity criterion of 95% for 50 bases, and a total of 20457 non-redundant sequences, 8503 contigs and 11954 singletons were generated. EST sequence coverage was analyzed by using the annotated L. japonicus genomic sequence and 1093 of the 1889 predicted protein-encoding genes (57.9%) were hit by the EST sequence(s). Gene content was compared to several plant species. Among the 8503 contigs, 471 were identified as sequences conserved only in leguminous species and these included several disease resistance-related genes. This suggested that in legumes, these genes may have evolved specifically to resist pathogen attack. The rate of gene sequence divergence was assessed by comparing similarity level and functional category based on the Gene Ontology (GO) annotation of Arabidopsis genes. This revealed that genes encoding ribosomal proteins, as well as those related to translation, photosynthesis, and cellular structure were more abundantly represented in the highly conserved class, and that genes encoding transcription factors and receptor protein kinases were abundantly represented in the less conserved class. To make the sequence information and the cDNA clones available to the research community, a Web database with useful services was created at http://www.kazusa.or.jp/en/plant/lotus/EST/.

DNA, Complementary↗

Complete genome sequencing of Anaplasma marginale reveals that the surface is skewed to two superfamilies of outer membrane proteins.

The rickettsia Anaplasma marginale is the most prevalent tick-borne livestock pathogen worldwide and is a severe constraint to animal health. A. marginale establishes lifelong persistence in infected ruminants and these animals serve as a reservoir for ticks to acquire and transmit the pathogen. Within the mammalian host, A. marginale generates antigenic variants by changing a surface coat composed of numerous proteins. By sequencing and annotating the complete 1,197,687-bp genome of the St. Maries strain of A. marginale, we show that this surface coat is dominated by two families containing immunodominant proteins: the msp2 superfamily and the msp1 superfamily. Of the 949 annotated coding sequences, just 62 are predicted to be outer membrane proteins, and of these, 49 belong to one of these two superfamilies. The genome contains unusual functional pseudogenes that belong to the msp2 superfamily and play an integral role in surface coat antigenic variation, and are thus distinctly different from pseudogenes described as byproducts of reductive evolution in other Rickettsiales.

Anaplasma marginale↗

Sequence of the genome of Salmonella bacteriophage P22.

The sequence of the nonredundant region of the Salmonella enterica serovar Typhimurium temperate, serotype-converting bacteriophage P22 has been completed. The genome is 41,724 bp with an overall moles percent GC content of 47.1%. Numerous examples of potential integration host factor and C1-binding sites were identified in the sequence. In addition, five potential rho-independent terminators were discovered. Sixty-five genes were identified and annotated. While many of these had been described previously, we have added several new ones, including the genes involved in serotype conversion and late control. Two of the serotype conversion gene products show considerable sequence relatedness to GtrA and -B from Shigella phages SfII, SfV, and SfX. We have cloned the serotype-converting cassette (gtrABC) and demonstrated that it results in Salmonella serovar Typhimurium LT2 cells which express antigen O1. Many of the putative proteins show sequence relatedness to proteins from a great variety of other phages, supporting the hypothesis that this phage has evolved through the recombinational exchange of genetic information with other viruses.

Bacteriophage P22↗

Functional annotation by identification of local surface similarities: a novel tool for structural genomics.

BACKGROUND: Protein function is often dependent on subsets of solvent-exposed residues that may exist in a similar three-dimensional configuration in non homologous proteins thus having different order and/or spacing in the sequence. Hence, functional annotation by means of sequence or fold similarity is not adequate for such cases. RESULTS: We describe a method for the function-related annotation of protein structures by means of the detection of local structural similarity with a library of annotated functional sites. An automatic procedure was used to annotate the function of local surface regions. Next, we employed a sequence-independent algorithm to compare exhaustively these functional patches with a larger collection of protein surface cavities. After tuning and validating the algorithm on a dataset of well annotated structures, we applied it to a list of protein structures that are classified as being of unknown function in the Protein Data Bank. By this strategy, we were able to provide functional clues to proteins that do not show any significant sequence or global structural similarity with proteins in the current databases. CONCLUSION: This method is able to spot structural similarities associated to function-related similarities, independently on sequence or fold resemblance, therefore is a valuable tool for the functional analysis of uncharacterized proteins. Results are available at http://cbm.bio.uniroma2.it/surface/structuralGenomics.html.

Algorithms↗

Defining a holoprosencephaly locus on human chromosome 14q13 and characterization of potential candidate genes.

Holoprosencephaly (HPE) is the most common developmental field defect in patterning of the human prosencephalon and associated craniofacial structures. The genetics is complex, with 12 loci defined on 11 chromosomes. We defined a locus for HPE (HPE8) on human chromosome 14q13 between markers D14S49 and AFM205XG5, by mapping deletion intervals of affected subjects with proximal chromosome 14q interstitial cytogenetic deletions. A 35-BAC contig was built by chromosome walking. By annotation of the 2.82-Mb minimal critical region, we identified 28 possible genes. Seven genes were expressed in human fetal brain: NPAS3, SNX6, C14ORF11, C14ORF10, PAX9, NKX2.1, and C14ORF19, the last an apparent gene fragment. Molecular embryology, animal modeling, and human mutation studies were reported elsewhere for PAX9 and NKX2.1. We focused on three genes, SNX6, NPAS3, and C14ORF11, as potential candidates for HPE. Genomic structure, human expression patterns, protein cellular localization, and embryonic expression patterns of orthologous murine genes were determined, showing that the three genes have properties similar to those of known HPE genes.

Animals↗

Genomic exploration of the hemiascomycetous yeasts: 3. Methods and strategies used for sequence analysis and annotation.

The primary analysis of the sequences for our Hemiascomycete random sequence tag (RST) project was performed using a combination of classical methods for sequence comparison and contig assembly, and of specifically written scripts and computer visualization routines. Comparisons were performed first against DNA and protein sequences from Saccharomyces cerevisiae, then against protein sequences from other completely sequenced organisms and, finally, against protein sequences from all other organisms. Blast alignments were individually inspected to help recognize genes within our random genomic sequences despite the fact that only parts of them were available. For each yeast species, validated alignments were used to infer the proper genetic code, to determine codon usage preferences and to calculate their degree of sequence divergence with S. cerevisiae. The quality of each genomic library was monitored from contig analysis of the DNA sequences. Annotated sequences were submitted to the EMBL database, and the general annotation tables produced served as a basis for our comparative description of the evolution, redundancy and function of the Hemiascomycete genomes described in other articles of this issue.

Amino Acid Sequence↗

A tool to assist the study of specific features at protein binding sites.

The Protein Data Bank contains a large amount of proteins that have been solved with small ligands bound to them. This constitutes a rich source of information for the study of the specific requirements of protein sites to bind small molecules with a favorable free energy. The specific atomic composition and three-dimensional geometric restraints of protein binding sites for different ligands could be easily obtained from there. The development of accurate binding site descriptors in proteins constitutes a valuable tool to assist in the large-scale prediction and annotation of protein function in whole genomes. In this work, an integrated database containing some processed and calculated protein/ligand information is described. It is expected that this database will constitute a useful tool for people working in the prediction of protein function from its structure. The database is accessible from the Internet through a web server located at: http://protein.bio.puc.cl

Binding Sites↗

On the inference of parsimonious indel evolutionary scenarios.

Given a multiple alignment of orthologous DNA sequences and a phylogenetic tree for these sequences, we investigate the problem of reconstructing a most parsimonious scenario of insertions and deletions capable of explaining the gaps observed in the alignment. This problem, called the Indel Parsimony Problem, is a crucial component of the problem of ancestral genome reconstruction, and its solution provides valuable information to many genome functional annotation approaches. We first show that the problem is NP-complete. Second, we provide an algorithm, based on the fractional relaxation of an integer linear programming formulation. The algorithm is fast in practice, and the solutions it produces are, in most cases, provably optimal. We describe a divide-and-conquer approach that makes it possible to solve very large instances on a simple desktop machine, while retaining guaranteed optimality. Our algorithms are tested and shown efficient and accurate on a set of 1.8 Mb mammalian orthologous sequences in the CFTR region.

Algorithms↗

inGeno--an integrated genome and ortholog viewer for improved genome to genome comparisons.

BACKGROUND: Systematic genome comparisons are an important tool to reveal gene functions, pathogenic features, metabolic pathways and genome evolution in the era of post-genomics. Furthermore, such comparisons provide important clues for vaccines and drug development. Existing genome comparison software often lacks accurate information on orthologs, the function of similar genes identified and genome-wide reports and lists on specific functions. All these features and further analyses are provided here in the context of a modular software tool "inGeno" written in Java with Biojava subroutines. RESULTS: InGeno provides a user-friendly interactive visualization platform for sequence comparisons (comprehensive reciprocal protein--protein comparisons) between complete genome sequences and all associated annotations and features. The comparison data can be acquired from several different sequence analysis programs in flexible formats. Automatic dot-plot analysis includes output reduction, filtering, ortholog testing and linear regression, followed by smart clustering (local collinear blocks; LCBs) to reveal similar genome regions. Further, the system provides genome alignment and visualization editor, collinear relationships and strain-specific islands. Specific annotations and functions are parsed, recognized, clustered, logically concatenated and visualized and summarized in reports. CONCLUSION: As shown in this study, inGeno can be applied to study and compare in particular prokaryotic genomes against each other (gram positive and negative as well as close and more distantly related species) and has been proven to be sensitive and accurate. This modular software is user-friendly and easily accommodates new routines to meet specific user-defined requirements.

Base Sequence↗

Identification of transcripts related to high egg production in the chicken hypothalamus and pituitary gland.

To identify transcripts related to high egg production expressed specifically in the hypothalamus and pituitary gland of the chicken, two subtracted cDNA libraries were constructed. Two divergently selected strains of Taiwan Country Chickens (TCCs), B (sire line) and L2 (dam line) were used; they had originated from a single population and were further subjected (since 1982) to selection for egg production to 40 wk of age and body weight/comb size, respectively. A total of 324 and 370 clones were identified from the L2-B (L2-subtract-B) and the B-L2 subtracted cDNA libraries, respectively. After sequencing and annotation, 175 and 136 transcripts that represented 53 known and 65 unknown non-redundant sequences were characterized in the L2-B subtracted cDNA library. Quantitative reverse-transcription (RT)-PCR was used to screen the mRNA expression levels of 32 randomly selected transcripts in another 78 laying hens from five different strains. These strains included the two original strains (B and L2) used to construct the subtracted cDNA libraries and an additional three commercial strains, i.e., Black- and Red-feather TCCs and Single-Comb White Leghorn (WL) layer. The mRNA expression levels of 16 transcripts were significantly higher in the L2 than in the B strain, whereas the mRNA expression levels of nine transcripts, BDH, NCAM1, PCDHA@, PGDS, PLAG1, PRL, SAR1A, SCG2 and STMN2, were significantly higher in two high egg production strains, L2 and Single-Comb WL; this indicated their usefulness as molecular markers of high egg production.

Animals↗

Plant Gene and Alternatively Spliced Variant Annotator. A plant genome annotation pipeline for rice gene and alternatively spliced variant identification with cross-species expressed sequence tag conservation from seven plant species.

The completion of the rice (Oryza sativa) genome draft has brought unprecedented opportunities for genomic studies of the world's most important food crop. Previous rice gene annotations have relied mainly on ab initio methods, which usually yield a high rate of false-positive predictions and give only limited information regarding alternative splicing in rice genes. Comparative approaches based on expressed sequence tags (ESTs) can compensate for the drawbacks of ab initio methods because they can simultaneously identify experimental data-supported genes and alternatively spliced transcripts. Furthermore, cross-species EST information can be used to not only offset the insufficiency of same-species ESTs but also derive evolutionary implications. In this study, we used ESTs from seven plant species, rice, wheat (Triticum aestivum), maize (Zea mays), barley (Hordeum vulgare), sorghum (Sorghum bicolor), soybean (Glycine max), and Arabidopsis (Arabidopsis thaliana), to annotate the rice genome. We developed a plant genome annotation pipeline, Plant Gene and Alternatively Spliced Variant Annotator (PGAA). Using this approach, we identified 852 genes (931 isoforms) not annotated in other widely used databases (i.e. the Institute for Genomic Research, National Center for Biotechnology Information, and Rice Annotation Project) and found 87% of them supported by both rice and nonrice EST evidence. PGAA also identified more than 44,000 alternatively spliced events, of which approximately 20% are not observed in the other three annotations. These novel annotations represent rich opportunities for rice genome research, because the functions of most of our annotated genes are currently unknown. Also, in the PGAA annotation, the isoforms with non-rice-EST-supported exons are significantly enriched in transporter activity but significantly underrepresented in transcription regulator activity. We have also identified potential lineage-specific and conserved isoforms, which are important markers in evolutionary studies. The data and the Web-based interface, RiceViewer, are available for public access at http://RiceViewer.genomics.sinica.edu.tw/.

Base Sequence↗

Identification, cloning, and expression of bacteriophage T5 dnk gene encoding a broad specificity deoxyribonucleoside monophosphate kinase (EC 2.7.4.13).

The nucleotide sequence corresponding to 13-19.5% of the bacteriophage T5 genome in early region C was determined (GenBank AY 140897). One of the five major single-stranded interruptions (nicks) of bacteriophage T5 DNA was identified at 18.5%. The sequenced region was annotated and the putative functions of some open reading frames were proposed by comparison with databases. The dnk gene, encoding a deoxyribonucleoside monophosphate kinase, was identified using a previously defined N-terminal amino acid sequence. The gene was cloned and expressed in Escherichia coli, the enzyme was purified to homogeneity with high yield using two alternative methods, and the recombinant deoxyribonucleoside monophosphate kinase was found to have the same activity and specificity as the native enzyme.

Amino Acid Sequence↗