Search PubMed⌕ Search

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.

At least 1,549 records · Page 86Linked to original sources
PubMed2006

Finding novel genes in bacterial communities isolated from the environment.

MOTIVATION: Novel sequencing techniques can give access to organisms that are difficult to cultivate using conventional methods. When applied to environmental samples, the data generated has some drawbacks, e.g. short length of assembled contigs, in-frame stop codons and frame shifts. Unfortunately, current gene finders cannot circumvent these difficulties. At the same time, the automated prediction of genes is a prerequisite for the increasing amount of genomic sequences to ensure progress in metagenomics. RESULTS: We introduce a novel gene finding algorithm that incorporates features overcoming the short length of the assembled contigs from environmental data, in-frame stop codons as well as frame shifts contained in bacterial sequences. The results show that by searching for sequence similarities in an environmental sample our algorithm is capable of detecting a high fraction of its gene content, depending on the species composition and the overall size of the sample. The method is valuable for hunting novel unknown genes that may be specific for the habitat where the sample is taken. Finally, we show that our algorithm can even exploit the limited information contained in the short reads generated by 454 technology for the prediction of protein coding genes. AVAILABILITY: The program is freely available upon request.

Algorithms↗
PubMed2004

Rapid sequence turnover at an intergenic locus in Drosophila.

Closely related species of Drosophila tend to have similar genome sizes. The strong imbalance in favor of small deletions relative to insertions implies that the unconstrained DNA in Drosophila is unlikely to be passively inherited from even closely related ancestors, and yet most DNA in Drosophila genomes is intergenic and potentially unconstrained. In an attempt to investigate the maintenance of this intergenic DNA, we studied the evolution of an intergenic locus on the fourth chromosome of the Drosophila melanogaster genome. This 1.2-kb locus is marked by two distinct, large insertion events: a nuclear transposition of a mitochondrial sequence and a transposition of a nonautonomous DNA transposon DNAREP1_DM. Because we could trace the evolutionary histories of these sequences, we were able to reconstruct the length evolution of this region in some detail. We sequenced this locus in all four species of the D. melanogaster species complex: D. melanogaster, D. simulans, D. sechellia, and D. mauritiana. Although this locus is similar in size in these four species, less than 10% of the sequence from the most recent common ancestor remains in D. melanogaster and all of its sister species. This region appears to have increased in size through several distinct insertions in the ancestor of the D. melanogaster species complex and has been shrinking since the split of these lineages. In addition, we found no evidence suggesting that the size of this locus has been maintained over evolutionary time; these results are consistent with the model of a dynamic equilibrium between persistent DNA loss through small deletions and more sporadic DNA gain through less frequent but longer insertions. The apparent stability of genome size in Drosophila may belie very rapid sequence turnover at intergenic loci.

Animals↗
PubMed2004

Evolution of prokaryotic DNA: intragenic and extragenic divergences observed with orthologs from three related species.

This study compared orthologous gene pairs from Escherichia coli K12, E. coli O157:H7 EDL933, Salmonella typhimurium LT2, and Yersinia pestis CO92 using only homologs of equal length, and differing nucleotides were counted and mapped. The data showed very clearly how the rates of divergence change with intragenic and extragenic position. The rate of synonymous mutation is lowest near start codons and near stop codons, and, a little surprisingly, the opposite seemed to be true for nonsynonymous substitutions. Analysis outside genes reveals that nucleotide divergences occur less frequently upstream of start codons than downstream of stop codons, and a remarkable drop in divergences is seen for two of the data sets around N = 9 nucleotides upstream of start codons; that is, the Shine-Dalgarno region changes at a lower rate. The explanation is likely to be the link between expressivity and sequence complementarity to the 3' end of 16S ribosomal rRNA. The latter is highly conserved across many bacterial and archaebacterial species.

Bacteria↗
PubMed2004

Selective expansion of the newly evolved genomic variants of retrotransposon 1731 in the Drosophila genomes.

The structural variants of the regulatory and coding regions of the LTR-retrotransposon 1731 are described. Two classes of genomic copies of retrotransposon 1731, with and without frameshifting strategy to express Gag and Pol proteins, were earlier revealed in the D. melanogaster genome. Copies without frameshifting are shown to be evolved from an ancient variant with frameshifting and are widespread in the genomes of the melanogaster complex species. Position of a rare codon responsible for ribosome pausing and efficient frameshifting is identified. Two structural variants of 1731 LTRs were detected in the melanogaster complex species: the predominant structural variant A1A2 of 1731 LTR in the D. melanogaster, D. simulans, and D. sechellia genomes contains duplicated and diverged copies of 28 bp in the U3 region, whereas A1 variant lacking this duplication is expanded in the D. mauritiana genome. Selective expansion of the A1A2 variant was detected in the independently established D. melanogaster cell cultures. A1A2 variant is expressed in embryos, cell culture, and testes, whereas A1 is expressed only in testes of D. melanogaster. Relief of expression of the A1A2 but not A1 variant in the ovaries as a result of mutation in the RNA interference (RNAi) spn-E gene is shown. Thus, expansion of the recently evolved genomic variants of the LTR retrotransposon 1731 possessing a new translation strategy, duplication in the U3 region, and extended profile of expression is revealed.

Animals↗
PubMed2005

Evolutionary implications of bacterial polyketide synthases.

Polyketide synthases (PKS) perform a stepwise biosynthesis of diverse carbon skeletons from simple activated carboxylic acid units. The products of the complex pathways possess a wide range of pharmaceutical properties, including antibiotic, antitumor, antifungal, and immunosuppressive activities. We have performed a comprehensive phylogenetic analysis of multimodular and iterative PKS of bacteria and fungi and of the distinct types of fatty acid synthases (FAS) from different groups of organisms based on the highly conserved ketoacyl synthase (KS) domains. Apart from enzymes that meet the classification standards we have included enzymes involved in the biosynthesis of mycolic acids, polyunsaturated fatty acids (PUFA), and glycolipids in bacteria. This study has revealed that PKS and FAS have passed through a long joint evolution process, in which modular PKS have a central position. They appear to have derived from bacterial FAS and primary iterative PKS and, in addition, share a common ancestor with animal FAS and secondary iterative PKS. Furthermore, we have carried out a phylogenomic analysis of all modular PKS that are encoded by the complete eubacterial genomes currently available in the database. The phylogenetic distribution of acyltransferase and KS domain sequences revealed that multiple gene duplications, gene losses, as well as horizontal gene transfer (HGT) have contributed to the evolution of PKS I in bacteria. The impact of these factors seems to vary considerably between the bacterial groups. Whereas in actinobacteria and cyanobacteria the majority of PKS I genes may have evolved from a common ancestor, several lines of evidence indicate that HGT has strongly contributed to the evolution of PKS I in proteobacteria. Discovery of new evolutionary links between PKS and FAS and between the different PKS pathways in bacteria may help us in understanding the selective advantage that has led to the evolution of multiple secondary metabolite biosyntheses within individual bacteria.

Amino Acid Sequence↗
PubMed1998

Combining diverse evidence for gene recognition in completely sequenced bacterial genomes.

Analysis of a newly sequenced bacterial genome starts with identification of protein-coding genes. Functional assignment of proteins requires the exact knowledge of protein N-termini. We present a new program ORPHEUS that identifies candidate genes and accurately predicts gene starts. The analysis starts with a database similarity search and identification of reliable gene fragments. The latter are used to derive statistical characteristics of protein-coding regions and ribosome-binding sites and to predict the complete set of genes in the analyzed genome. In a test on Bacillus subtilis and Escherichia coli genomes, the program correctly identified 93.3% (resp. 96.3%) of experimentally annotated genes longer than 100 codons described in the PIR-International database, and for these genes 96.3% (83.9%) of starts were predicted exactly. Furthermore, 98.9% (99.1%) of genes longer than 100 codons annotated in GenBank were found, and 92.9% (75.7%) of predicted starts coincided with the feature table description. Finally, for the complete gene complements of B.subtilis and E.coli , including genes shorter than 100 codons, gene prediction accuracy was 88.9 and 87.1%, respectively, with 94.2 and 76.7% starts coinciding with the existing annotation.

Algorithms↗
PubMed2000

Analysis of the Thermotoga maritima genome combining a variety of sequence similarity and genome context tools.

The proliferation of genome sequence data has led to the development of a number of tools and strategies that facilitate computational analysis. These methods include the identification of motif patterns, membership of the query sequences in family databases, metabolic pathway involvement and gene proximity. We re-examined the completely sequenced genome of Thermotoga maritima by employing the combined use of the above methods. By analyzing all 1877 proteins encoded in this genome, we identified 193 cases of conflicting annotations (10%), of which 164 are new function predictions and 29 are amendments of previously proposed assignments. These results suggest that the combined use of existing computational tools can resolve inconclusive sequence similarities and significantly improve the prediction of protein function from genome sequence.

Computational Biology↗
PubMed2001

Phylogenetic footprinting of transcription factor binding sites in proteobacterial genomes.

Toward the goal of identifying complete sets of transcription factor (TF)-binding sites in the genomes of several gamma proteobacteria, and hence describing their transcription regulatory networks, we present a phylogenetic footprinting method for identifying these sites. Probable transcription regulatory sites upstream of Escherichia coli genes were identified by cross-species comparison using an extended Gibbs sampling algorithm. Close examination of a study set of 184 genes with documented transcription regulatory sites revealed that when orthologous data were available from at least two other gamma proteobacterial species, 81% of our predictions corresponded with the documented sites, and 67% corresponded when data from only one other species were available. That the remaining predictions included bona fide TF-binding sites was proven by affinity purification of a putative transcription factor (YijC) bound to such a site upstream of the fabA gene. Predicted regulatory sites for 2097 E.coli genes are available at http://www.wadsworth.org/resnres/bioinfo/.

Base Sequence↗
PubMed2003

Efficient clustering of large EST data sets on parallel computers.

Clustering expressed sequence tags (ESTs) is a powerful strategy for gene identification, gene expression studies and identifying important genetic variations such as single nucleotide polymorphisms. To enable fast clustering of large-scale EST data, we developed PaCE (for Parallel Clustering of ESTs), a software program for EST clustering on parallel computers. In this paper, we report on the design and development of PaCE and its evaluation using Arabidopsis ESTs. The novel features of our approach include: (i) design of memory efficient algorithms to reduce the memory required to linear in the size of the input, (ii) a combination of algorithmic techniques to reduce the computational work without sacrificing the quality of clustering, and (iii) use of parallel processing to reduce run-time and facilitate clustering of larger data sets. Using a combination of these techniques, we report the clustering of 168 200 Arabidopsis ESTs in 15 min on an IBM xSeries cluster with 30 dual-processor nodes. We also clustered 327 632 rat ESTs in 47 min and 420 694 Triticum aestivum ESTs in 3 h and 15 min. We demonstrate the quality of our software using benchmark Arabidopsis EST data, and by comparing it with CAP3, a software widely used for EST assembly. Our software allows clustering of much larger EST data sets than is possible with current software. Because of its speed, it also facilitates multiple runs with different parameters, providing biologists a tool to better analyze EST sequence data. Using PaCE, we clustered EST data from 23 plant species and the results are available at the PlantGDB website.

Algorithms↗
PubMed2003

ETOPE: Evolutionary test of predicted exons.

Since a large number of computationally predicted exons are not supported by existing sequence (e.g. ESTs) or experimental (e.g. expression analysis) data they need to be validated by other methods. ETOPE is designed to test computational predictions by using signals that have not been included in any current computational prediction method. The test is based on the ratio of non-synonymous to synonymous substitution rates between sequences from different genomes. It has been previously shown, by empirical data and computer simulation, to be a powerful criterion for identifying protein-coding regions. The ETOPE is available at http://nekrut.uchicago.edu/etope/.

Amino Acid Substitution↗
PubMed2004

VirGen: a comprehensive viral genome resource.

VirGen is a comprehensive viral genome resource that organizes the 'sequence space' of viral genomes in a structured fashion. It has been developed with the objective of serving as an annotated and curated database comprising complete genome sequences of viruses, value-added derived data and data mining tools. The current release (v1.1) contains 559 complete genomes in addition to 287 putative genomes of viruses belonging to eight viral families for which the host range includes animals and plants. Viral genomes in VirGen are annotated using sequence-based Bioinformatics approaches. The genomic data is also curated to identify 'alternate names' of viral proteins, where available. VirGen archives the results of comparisons of genomes, proteomes and individual proteins within and between viral species. It is the first resource to provide phylogenetic trees of viral species computed using whole-genome sequence data. The module of predicted B-cell antigenic determinants in VirGen is an attempt to link the genome to its vaccinome. Comparative genome analysis data facilitate the study of genome organization and evolution of viruses, which would have implications in applied research to identify candidates for the design of vaccines and antiviral drugs. VirGen is a relational database and is available at http://bioinfo. ernet.in/virgen/virgen.html.

Antigens, Viral↗
PubMed2005

Ensembl 2005.

The Ensembl (http://www.ensembl.org/) project provides a comprehensive and integrated source of annotation of large genome sequences. Over the last year the number of genomes available from the Ensembl site has increased by 7 to 16, with the addition of the six vertebrate genomes of chimpanzee, dog, cow, chicken, tetraodon and frog and the insect genome of honeybee. The majority have been annotated automatically using the Ensembl gene build system, showing its flexibility to reliably annotate a wide variety of genomes. With the increased number of vertebrate genomes, the comparative analysis provided to users has been greatly improved, with new website interfaces allowing annotation of different genomes to be directly compared. The Ensembl software system is being increasingly widely reused in different projects showing the benefits of a completely open approach to software development and distribution.

Animals↗
PubMed2005

AMOD: a morpholino oligonucleotide selection tool.

AMOD is a web-based program that aids in the functional evaluation of nucleotide sequences through sequence characterization and antisense morpholino oligonucleotide (target site) selection. Submitted sequences are analyzed by translation initiation site prediction algorithms and sequence-to-sequence comparisons; results are used to characterize sequence features required for morpholino design. Within a defined subsequence, base composition and homodimerization values are computed for all putative morpholino oligonucleotides. Using these properties, morpholino candidates are selected and compared with genomic and transcriptome databases with the goal to identify target-specific enriched morpholinos. AMOD has been used at the University of Minnesota to design approximately 200 morpholinos for a functional genomics screen in zebrafish. The AMOD web server and a tutorial are freely available to both academic and commercial users at http://www.secretomes.umn.edu/AMOD/.

DNA Primers↗
PubMed2006

Data mining for proteins characteristic of clades.

A synapomorphy is a phylogenetic character that provides evidence of shared descent. Ideally a synapomorphy is ubiquitous within the clade of related organisms and nonexistent outside the clade, implying that it arose after divergence from other extant species and before the last common ancestor of the clade. With the recent proliferation of genetic sequence data, molecular synapomorphies have assumed great importance, yet there is no convenient means to search for them over entire genomes. We have developed a new program called Conserv, which can rapidly assemble orthologous sequences and rank them by various metrics, such as degree of conservation or divergence from out-group orthologs. We have used Conserv to conduct a largescale search for molecular synapomorphies for bacterial clades. The search discovered sequences unique to clades, such as Actinobacteria, Firmicutes and gamma-Proteobacteria, and shed light on several open questions, such as whether Symbiobacterium thermophilum belongs with Actinobacteria or Firmicutes. We conclude that Conserv can quickly marshall evidence relevant to evolutionary questions that would be much harder to assemble with other tools.

Amino Acid Sequence↗
PubMed2005

A simple, bead-based approach for multi-SNP molecular haplotyping.

Single nucleotide polymorphisms (SNPs) within a gene region have often been studied to determine their effect on phenotype. Although a single base pair change can produce a phenotypic change, phenotype is often influenced by the presence of multiple polymorphisms and their relative positions within a given region. For example, if multiple changes occur in a promoter region, how they influence gene expression will depend on their cis/trans configuration. As such, it is essential to consider the haplotype, or the alignment of multiple SNP alleles on each chromosome when attempting to associate genomic changes with phenotype. Unfortunately, no method of high-throughput molecular haplotyping of multiple SNPs currently exists. In response to this unmet need, we have developed an inexpensive, reliable bead-based capture-based haplotyping (CBH) assay to determine the phase, or haplotype, of multiple SNP alleles in a high-throughput manner. The CBH assay requires minimal setup and handling, requires no centrifugation steps and can be performed in <1 h. Data collection is performed via flow cytometry and the assay yields plus/minus results allowing for automated calling by a simple computer application. We will present data demonstrating the molecular haplotyping of 11 SNPs within exon 2 of the N-acetyltransferase-2 (NAT2) gene, which expresses an important drug-metabolizing enzyme. This assay has applications in diagnostic testing, promoter analysis, association studies and pharmacogenetic analysis.

Arylamine N-Acetyltransferase↗
PubMed2002

A comprehensive vertebrate phylogeny using vector representations of protein sequences from whole genomes.

We recently developed a method for producing comprehensive gene and species phylogenies from unaligned whole genome data using singular value decomposition (SVD) to analyze character string frequencies. This work provides an integrated gene and species phylogeny for 64 vertebrate mitochondrial genomes composed of 832 total proteins. In addition, to provide a theoretical basis for the method, we present a graphical interpretation of both the original frequency matrix and the SVD-derived matrix. These large matrices describe high-dimensional Euclidean spaces within which biomolecular sequences can be uniquely represented as vectors. In particular, the SVD-derived vector space describes each protein relative to a restricted set of newly defined, independent axes, each of which represents a novel form of conserved motif, termed a correlated peptide motif. A quantitative comparison of the relative orientations of protein vectors in this space provides accurate and straightforward estimates of sequence similarity, which can in turn be used to produce comprehensive gene trees. Alternatively, the vector representations of genes from individual species can be summed, allowing species trees to be produced.

Algorithms↗
PubMed1996

Molecular phylogeny and genome evolution in the Drosophila virilis species group: duplications of the alcohol dehydrogenase gene.

Drosophila virilis is a prominent reference species for comparison with Drosophila melanogaster in regard to patterns and mechanisms of molecular and genomic evolution. Sequences were determined for 11 Adh genes from 8 species of the D. virilis species group, including species from both the virilis phylad and the montana subphylad. The genome of D. virilis contains a 6-kb duplication that includes the entire Adh coding region. The pattern of sequence identity within the duplication strongly suggests a recent gene-conversion event bordered by 36-bp indels. As in other Drosophila, the amino-acid coding region of Adh is encoded by three exons interrupted by two short introns. The promoter region includes 16 blocks of sequence that are well conserved in D. virilis, Drosophila hydei, and D. melanogaster. The developmental profile of Adh transcription suggests a distal/proximal promoter switch analogous to that in D. melanogaster. Duplicate Adh genes were also found in Drosophila montana and Drosophila lacicola, which apparently originated independently of that in D. virilis. The Adh genes in all species of the D. virilis group have among the lowest codon bias of any Adh genes so far reported in the genus Drosophila. Taking the low codon bias into account, we estimate the time of divergence between the virilis and montana clades as 9.0 +/- 0.7 Mya and the approximate time of divergence of D. virilis from other members of the virilis phylad as 2.6 +/- 0.4 Mya. The region of the D. virilis genome containing Adh, as well as the chromosome as a whole, gives evidence of extensive rearrangements relative to the genome of D. melanogaster.

Alcohol Dehydrogenase↗
PubMed1998

Molecular epidemiology of rabbit haemorrhagic disease virus outbreaks in France during 1988 to 1995.

In order to evaluate genetic variation between rabbit haemorrhagic disease virus (RHDV) isolates and to derive phylogenetic relationships, 56 virus isolates collected from various parts of France over a 7 year period (1988 to 1995) were examined. Analyses were carried out by direct nucleotide sequencing of PCR fragments of three genomic regions encoding the capsid protein (VP60) (regions A and B) and a non-structural protein (region C). Multiple sequence alignments revealed maximum nucleotide divergence of 7.6, 9.4 and 8.7% for regions A, B and C, respectively, indicating a high level of conservation between isolates. Irrespective of the genomic region analysed, phylogenetic analyses carried out using various methods allowed the identification of three genogroups; distribution of isolates within these genogroups appears to be more related to the year of their collection than to their geographical origin. The possible evolution of RHDV is discussed.

Animals↗