The evolution of structural databases

Starting with the Protein Data Bank (PDB) as a common ancestor, the evolution of structural databases has been driven by the rapprochement of the structural world and the practical applications. The result is an impressive number of secondary structural databases that is welcomed by structural biologists and bioinformaticians but runs the risk of producing an embarrassment of riches among non-specialist users. Given that any profit depends on the number of customers, efficient interfaces between many structural data banks must be available to make their contents easily accessible. Increasing the information content of central structural repositories might be the best way to guide users through the many, sometimes overlapping databases.     

Protein sequence databases

A variety of protein sequence databases exist, ranging from simple sequence repositories, which store data with little or no manual intervention in the creation of the records, to expertly curated universal databases that cover all species and in which the original sequence data are enhanced by the manual addition of further information in each sequence record. As the focus of researchers moves from the genome to the proteins encoded by it, these databases will play an even more important role as central comprehensive resources of protein information. Several the leading protein sequence databases are discussed here, with special emphasis on the databases now provided by the Universal Protein Knowledgebase (UniProt) consortium.     

Short-term sequence evolution and vertical inheritance of the Naegleria twin-ribozyme group I intron

Background  

Ribosomal DNA of several species of the free-living Naegleria amoeba harbors an optional group I intron within the nuclear small subunit ribosomal RNA gene. The intron (Nae.S516) has a complex organization of two ribozyme domains (NaGIR1 and NaGIR2) and a homing endonuclease gene (NaHEG). NaGIR2 is responsible for intron excision, exon ligation, and full-length intron RNA circularization, reactions typical for nuclear group I intron ribozymes. NaGIR1, however, is essential for NaHEG expression by generating the 5' end of the homing endonuclease messenger RNA. Interestingly, this unusual class of ribozyme adds a lariat-cap at the mRNA.     

Patterns of intron sequence evolution in Drosophila are dependent upon length and GC content

Background

Introns comprise a large fraction of eukaryotic genomes, yet little is known about their functional significance. Regulatory elements have been mapped to some introns, though these are believed to account for only a small fraction of genome wide intronic DNA. No consistent patterns have emerged from studies that have investigated general levels of evolutionary constraint in introns.

Results

We examine the relationship between intron length and levels of evolutionary constraint by analyzing inter-specific divergence at 225 intron fragments in Drosophila melanogaster and Drosophila simulans, sampled from a broad distribution of intron lengths. We document a strongly negative correlation between intron length and divergence. Interestingly, we also find that divergence in introns is negatively correlated with GC content. This relationship does not account for the correlation between intron length and divergence, however, and may simply reflect local variation in mutational rates or biases.

Conclusion

Short introns make up only a small fraction of total intronic DNA in the genome. Our finding that long introns evolve more slowly than average implies that, while the majority of introns in the Drosophila genome may experience little or no selective constraint, most intronic DNA in the genome is likely to be evolving under considerable constraint. Our results suggest that functional elements may be ubiquitous within longer introns and that these introns may have a more general role in regulating gene expression than previously appreciated. Our finding that GC content and divergence are negatively correlated in introns has important implications for the interpretation of the correlation between divergence and levels of codon bias observed in Drosophila.     

Human sequence variation and mutation databases

This paper aims to give an overview of current resources onhuman sequence variations and give an idea about the directionin which these services are moving.      

Nucleic acid and protein sequence databases

Nucleic acid and protein sequences contain a wealth of informationof interest to molecular biologists. The advent of molecularsequence databases provides a unique opportunity for the computeranalysis of all available sequences. Sequence databases servetwo main functions: (i) to facilitate comparisons with newlydetermined sequences, and (ii) to act as a source of data forthe generation and testing of hypotheses concerning molecularsequence organisation and evolution. The large amounts of sequencedata now becoming available require that algorithms for databasesearching be fast and efficient and considerable progress isbeing made in this area.     

The role of pattern databases in sequence analysis

In the wake of the numerous now-fruitful genome projects, we are entering an era rich in biological data. The field of bioinformatics is poised to exploit this information in increasingly powerful ways, but the abundance and growing complexity both of the data and of the tools and resources required to analyse them are threatening to overwhelm us. Databases and their search tools are now an essential part of the research environment. However, the rate of sequence generation and the haphazard proliferation of databases have made it difficult to keep pace with developments. In an age of information overload, researchers want rapid, easy-to-use, reliable tools for functional characterisation of newly determined sequences. But what are those tools? How do we access them? Which should we use? This review focuses on a particular type of database that is increasingly used in the task of routine sequence analysis--the so-called pattern database. The paper aims to provide an overview of the current status of pattern databases in common use, outlining the methods behind them and giving pointers on their diagnostic strengths and weaknesses.     

Strategies for searching sequence databases

We provide a detailed overview of the choices inherent in performing a sequence database search, including the choice of algorithm, substitution matrix and gap model. Each of these choices has implications that can be described as restrictions on the underlying model of sequence evolution, the expected degree of divergence between the query sequence and the database sequences (if one uses an evolutionary based matrix), as well as the sensitivity and selectivity of the search. We conclude with a series of recommendations for researchers performing these searches based on our experience and literature studies.     

Geographic isolation and evolution of Mediterranean endemic Cyclamen: insights from chloroplast trnL (UAA) intron sequence variation

 In this study we construct a phylogenetic hypothesis for the relatedness among disjunct subspecies of Cyclamen repandum and its two allopatric congeners, C. creticum and C. balearicum in order to examine the evolutionary divergence of currently isolated populations across the western Mediterranean. The most parsimonious phylogenetic tree obtained from sequencing the cpDNA trnL (UAA) intron suggests a major phylogeographic divide in southern Greece between two clades. The first clade comprises samples of C. repandum subsp. peloponnesiacum (from the Peloponnese) and C. creticum (from Crete). The second comprises samples of C. repandum subsp. repandum (from Croatia, Italy, southern France, Corsica, Sardinia and Sicily), C. repandum subsp. rhodense (from Rhodes and Kos) and C. balearicum (from the Balearic Islands and southern France). These data suggest that C. creticum has evolved in allopatry from C. repandum subsp. peloponnesiacum and that C. balearicum and C. repandum ssp. rhodense have diverged from C. repandum subsp. repandum at its western and eastern distribution limits. At one small site on Corsica, a population of C. repandum may have introgressed with relictual populations of C. balearicum. These divergence patterns illustrate how a phylogenetic perspective can be used to better understand the evolution of endemism in the Mediterranean flora. Received February 19, 2001 Accepted August 22, 2001     

Patterns of intron sequence evolution in Drosophila are dependent upon length and GC content

Background  

Introns comprise a large fraction of eukaryotic genomes, yet little is known about their functional significance. Regulatory elements have been mapped to some introns, though these are believed to account for only a small fraction of genome wide intronic DNA. No consistent patterns have emerged from studies that have investigated general levels of evolutionary constraint in introns.     

Occurrence of Can-SINEs and intron sequence evolution supports robust phylogeny of pinniped carnivores and their terrestrial relatives

Investigating the dog genome we found 178 965 introns with a moderate length of 200–1000 bp. A screening of these sequences against 23 different repeat libraries to find insertions of short interspersed elements (SINEs) detected 45 276 SINEs. Virtually all of these SINEs (98%) belong to the tRNA-derived Can-SINE family. Can-SINEs arose about 55 million years ago before Carnivora split into two basal groups, the Caniformia (dog-like carnivores) and the Feliformia (cat-like carnivores). Genome comparisons of dog and cat recovered 506 putatively informative SINE loci for caniformian phylogeny. In this study we show how to use such genome information of model organisms to research the phylogeny of related non-model species of interest. Investigating a dataset including representatives of all major caniformian lineages, we analysed 24 randomly chosen loci for 22 taxa. All loci were amplifiable and revealed 17 parsimony-informative SINE insertions. The screening for informative SINE insertions yields a large amount of sequence information, in particular of introns, which contain reliable phylogenetic information as well. A phylogenetic analysis of intron- and SINE sequence data provided a statistically robust phylogeny which is congruent with the absence/presence pattern of our SINE markers. This phylogeny strongly supports a sistergroup relationship of Musteloidea and Pinnipedia. Within Pinnipedia, we see strong support from bootstrapping and the presence of a SINE insertion for a sistergroup relationship of the walrus with the Otariidae.     

Origin and evolution of the chloroplast trnK (matK) intron: a model for evolution of group II intron RNA structures

The trnK intron of plants encodes the matK open reading frame (ORF), which has been used extensively as a phylogenetic marker for classification of plants. Here we examined the evolution of the trnK intron itself as a model for group II intron evolution in plants. Representative trnK intron sequences were compiled from species spanning algae to angiosperms, and four introns were newly sequenced. Phylogenetic analyses showed that the matK ORFs belong to the ML (mitochondrial-like) subclass of group II intron ORFs, indicating that they were derived from a mobile group II intron of the class. RNA structures of the introns were folded and analyzed, which revealed progressive RNA structural deviations and degenerations throughout plant evolution. The data support a model in which plant organellar group II introns were derived from bacterial-like introns that had "standard" RNA structures and were competent for self-splicing and mobility and that subsequently the ribozyme structures degenerated to ultimately become dependent upon host-splicing factors. We propose that the patterns of RNA structure evolution seen for the trnK intron will apply to the other group II introns in plants.     

A group II intron in the Neurospora mitochondrial coI gene: nucleotide sequence and implications for splicing and molecular evolution.

The temperature-sensitive Neurospora nuclear mutant cyt18-1 is deficient in splicing many Group I mitochondrial introns when grown at its non-permissive temperature; however, splicing of intron 1 in the coI gene of the Adiopodoume (formerly called North Africa) strain is unaffected (R.A. Collins and A.M. Lambowitz, J. Mol. Biol. 184: 413-428, 1985). Here we show that coI intron 1 is a typical Group II intron, the only one identified to date in Neurospora. The differential effect of the cyt18-1 mutation suggests that splicing of certain introns could be regulated independently of others by nuclear-encoded proteins. The intron contains a long open reading frame (ORF) resembling that of the Neurospora Mauriceville mitochondrial plasmid. The intron and plasmid ORFs share unusual features of codon usage that suggest both evolved outside of the Neurospora mitochondrial genetic system.     

Compressing DNA sequence databases with coil

Background  

Publicly available DNA sequence databases such as GenBank are large, and are growing at an exponential rate. The sheer volume of data being dealt with presents serious storage and data communications problems. Currently, sequence data is usually kept in large flat files, which are then compressed using standard Lempel-Ziv (gzip) compression – an approach which rarely achieves good compression ratios. While much research has been done on compressing individual DNA sequences, surprisingly little has focused on the compression of entire databases of such sequences. In this study we introduce the sequence database compression software coil.