The accuracy of ribosomal RNA comparative structure models

The determination of the 16S and 23S rRNA secondary structure models was initiated shortly after the first complete 16S and 23S rRNA sequences were determined in the late 1970s. The structures that are common to all 16S rRNAs and all 23S rRNAs were determined using comparative methods from the analysis of thousands of rRNA sequences. Twenty-plus years later, the 16S and 23S rRNA comparative structure models have been evaluated against the recently determined high-resolution crystal structures of the 30S and 50S ribosomal subunits. Nearly all of the predicted covariation-based base pairs, including the regular base pairs and helices, and the irregular base pairs and tertiary interactions, were present in the 30S and 50S crystal structures.     

RNA secondary structure prediction based on free energy and phylogenetic analysis.

We describe a computational method for the prediction of RNA secondary structure that uses a combination of free energy and comparative sequence analysis strategies. Using a homology-based sequence alignment as a starting point, all favorable pairings with respect to the Turner energy function are identified. Each potentially paired region within a multiple sequence alignment is scored using a function that combines both predicted free energy and sequence covariation with optimized weightings. High scoring regions are ranked and sequentially incorporated to define a growing secondary structure. Using a single set of optimized parameters, it is possible to accurately predict the foldings of several test RNAs defined previously by extensive phylogenetic and experimental data (including tRNA, 5 S rRNA, SRP RNA, tmRNA, and 16 S rRNA). The algorithm correctly predicts approximately 80% of the secondary structure. A range of parameters have been tested to define the minimal sequence information content required to accurately predict secondary structure and to assess the importance of individual terms in the prediction scheme. This analysis indicates that prediction accuracy most strongly depends upon covariational information and only weakly on the energetic terms. However, relatively few sequences prove sufficient to provide the covariational information required for an accurate prediction. Secondary structures can be accurately defined by alignments with as few as five sequences and predictions improve only moderately with the inclusion of additional sequences.     

Nucleotide sequence of a crustacean 18S ribosomal RNA gene and secondary structure of eukaryotic small subunit ribosomal RNAs

The primary structure of the gene for 18 S rRNA of the crustacean Artemia salina was determined. The sequence has been aligned with 13 other small ribosomal subunit RNA sequences of eukaryotic, archaebacterial, eubacterial, chloroplastic and plant mitochondrial origin. Secondary structure models for these RNAs were derived on the basis of previously proposed models and additional comparative evidence found in the alignment. Although there is a general similarity in the secondary structure models for eukaryotes and prokaryotes, the evidence seems to indicate a different topology in a central area of the structures.     

An Empirical Analysis of mt 16S rRNA Covarion-Like Evolution in Insects: Site-Specific Rate Variation Is Clustered and Frequently Detected

The structural and functional analysis of rRNA molecules has attracted considerable scientific interest. Empirical studies have demonstrated that sequence variation is not directly translated into modifications of rRNA secondary structure. Obviously, the maintenance of secondary structure and sequence variation are in part governed by different selection regimes. The nature of those selection regimes still remains quite elusive. The analysis of individual bacterial models cannot adequately explore this topic. Therefore, we used primary sequence data and secondary structures of a mitochondrial 16S rRNA fragment of 558 insect species from 15 monophyletic groups to study patterns of sequence variation, and variation of secondary structure. Using simulation studies to establish significance levels of change, we found that despite conservation of secondary structure, the location of sequence variation within the conserved rRNA structure changes significantly between groups of insects. Despite our conservative estimation procedure we found significant site-specific rate changes at 56 sites out of 184. Additionally, site-specific rate variation is somewhat clustered in certain helices. Both results confirm what has been predicted from an application of non-stationary maximum likelihood models to rRNA sequences. Clearly, constraints on sequence variation evolve and leave footprints in the form of evolutionary plasticity in rRNA sequences. Here, we show that a better understanding of the evolution of rRNA sequences can be obtained by integrating both phylogenetic and structural information.     

Sequences of the 5S rRNAs of Azotobacter vinelandii, Pseudomonas aeruginosa and Pseudomonas fluorescens with some notes on 5S RNA secondary structure

Recently published alignments of available 5 S rRNA sequences have shown that a rigid base pairing pattern, pointing to the existence of a universal five-helix secondary structure for all 5 S RNAs, can be superimposed on such alignments. For a few species, the alignment and the base pairing pattern show distortions with respect to the large majority of sequences. Their 5 S RNAs may form exceptional secondary structures, or there may just be errors in the published sequences. We have examined such a case, Pseudomonas fluorescens, and found the sequence to be in error. The corrected sequence, as well as those of the related species Azotobacter vinelandii and Pseudomonas aeruginosa, fit perfectly in the 5 S RNA sequence alignment and in the five-helix secondary structure model. There exists comparative evidence for the frequent presence of non-standard base pairs at several points of the 5 S RNA secondary structure.     

Evolution of arthropod hemocyanins and insect storage proteins (hexamerins)

Crustacean and cheliceratan hemocyanins (oxygen-transport proteins) and insect hexamerins (storage proteins) are homologous gene products, although the latter do not bind oxygen and do not possess the copper- binding histidines present in the hemocyanins. An alignment of 19 amino acid sequences of hemocyanin subunits and insect hexamerins was made, based on the conservation of elements of secondary structure observed in X-ray structures of two hemocyanin subunits. The alignment was analyzed using parsimony and neighbor-joining methods. Results provide strong indications for grouping together the sequences of the 2 crustacean hemocyanin subunits, the 5 cheliceratan hemocyanin subunits, and the 12 insect hexamerins. Within the insect clade, four methionine- rich proteins, four arylphorins, and two juvenile hormone-suppressible proteins from Lepidoptera, as well as two dipteran proteins, form four separate groups. In the absence of an outgroup sequence, it is not possible to present information about the ancestral state from which these proteins are derived. Although this family of proteins clearly consists of homologous gene products, there remain striking differences in gene organization and site of biosynthesis of the proteins within the cell. Because studies on 18S and 12S rRNA sequences indicate a rather close relationship between insects and crustaceans, we propose that hemocyanin is the ancestral arthropod protein and that insect hexamerins lost their copper-binding capability after divergence of the insects from the crustaceans.      

Tensor decomposition reveals concurrent evolutionary convergences and divergences and correlations with structural motifs in ribosomal RNA

Evolutionary relationships among organisms are commonly described by using a hierarchy derived from comparisons of ribosomal RNA (rRNA) sequences. We propose that even on the level of a single rRNA molecule, an organism's evolution is composed of multiple pathways due to concurrent forces that act independently upon different rRNA degrees of freedom. Relationships among organisms are then compositions of coexisting pathway-dependent similarities and dissimilarities, which cannot be described by a single hierarchy. We computationally test this hypothesis in comparative analyses of 16S and 23S rRNA sequence alignments by using a tensor decomposition, i.e., a framework for modeling composite data. Each alignment is encoded in a cuboid, i.e., a third-order tensor, where nucleotides, positions and organisms, each represent a degree of freedom. A tensor mode-1 higher-order singular value decomposition (HOSVD) is formulated such that it separates each cuboid into combinations of patterns of nucleotide frequency variation across organisms and positions, i.e., "eigenpositions" and corresponding nucleotide-specific segments of "eigenorganisms," respectively, independent of a-priori knowledge of the taxonomic groups or rRNA structures. We find, in support of our hypothesis that, first, the significant eigenpositions reveal multiple similarities and dissimilarities among the taxonomic groups. Second, the corresponding eigenorganisms identify insertions or deletions of nucleotides exclusively conserved within the corresponding groups, that map out entire substructures and are enriched in adenosines, unpaired in the rRNA secondary structure, that participate in tertiary structure interactions. This demonstrates that structural motifs involved in rRNA folding and function are evolutionary degrees of freedom. Third, two previously unknown coexisting subgenic relationships between Microsporidia and Archaea are revealed in both the 16S and 23S rRNA alignments, a convergence and a divergence, conferred by insertions and deletions of these motifs, which cannot be described by a single hierarchy. This shows that mode-1 HOSVD modeling of rRNA alignments might be used to computationally predict evolutionary mechanisms.     

Evaluation of sequence alignments and oligonucleotide probes with respect to three-dimensional structure of ribosomal RNA using ARB software package

Background  

Availability of high-resolution RNA crystal structures for the 30S and 50S ribosomal subunits and the subsequent validation of comparative secondary structure models have prompted the biologists to use three-dimensional structure of ribosomal RNA (rRNA) for evaluating sequence alignments of rRNA genes. Furthermore, the secondary and tertiary structural features of rRNA are highly useful and successfully employed in designing rRNA targeted oligonucleotide probes intended for in situ hybridization experiments. RNA3D, a program to combine sequence alignment information with three-dimensional structure of rRNA was developed. Integration into ARB software package, which is used extensively by the scientific community for phylogenetic analysis and molecular probe designing, has substantially extended the functionality of ARB software suite with 3D environment.     

18S rRNA secondary structure and phylogenetic position of Peloridiidae (Insecta, hemiptera)

A secondary structure model for 18S rRNA of peloridiids, relict insects with a present-day circumantarctic distribution, is constructed using comparative sequence analysis, thermodynamic folding, a consensus method using 18S rRNA models of other taxa, and support of helices based on compensatory substitutions. Results show that probable in vivo configuration of 18S rRNA is not predictable using current free-energy models to fold the entire molecule concurrently. This suggests that refinements in free-energy minimization algorithms are needed. Molecular phylogenetic datasets were created using 18S rRNA nucleotide alignments produced by CLUSTAL and rigorous interpretation of homologous position based on certain secondary substructures. Phylogenetic analysis of a hemipteran data matrix of 18S rDNA sequences placed peloridiids sister to Heteroptera. Resolution of affiliations between the three main euhemipteran lineages was unresolved. The peloridiid 18S RNA model presented here provides the most accurate template to date for aligning homologous nucleotides of hemipteran taxa. Using folded 18S rRNA to infer homology of character as morpho-molecular structures or nucleotides and scoring particular sites or substructures is discussed.     

The European database on small subunit ribosomal RNA

The European database on SSU rRNA can be consulted via the World WideWeb at http://rrna.uia.ac.be/ssu/ and compiles all complete or nearly complete small subunit ribosomal RNA sequences. Sequences are provided in aligned format. The alignment takes into account the secondary structure information derived by comparative sequence analysis of thousands of sequences. Additional information such as literature references, taxonomy, secondary structure models and nucleotide variability maps, is also available.     

A secondary structural model of the 28S rRNA expansion segments D2 and D3 for Chalcidoid wasps (Hymenoptera: Chalcidoidea)

We analyze the secondary structure of two expansion segments (D2, D3) of the 28S ribosomal (rRNA)-encoding gene region from 527 chalcidoid wasp taxa (Hymenoptera: Chalcidoidea) representing 18 of the 19 extant families. The sequences are compared in a multiple sequence alignment, with secondary structure inferred primarily from the evidence of compensatory base changes in conserved helices of the rRNA molecules. This covariation analysis yielded 36 helices that are composed of base pairs exhibiting positional covariation. Several additional regions are also involved in hydrogen bonding, and they form highly variable base-pairing patterns across the alignment. These are identified as regions of expansion and contraction or regions of slipped-strand compensation. Additionally, 31 single-stranded locales are characterized as regions of ambiguous alignment based on the difficulty in assigning positional homology in the presence of multiple adjacent indels. Based on comparative analysis of these sequences, the largest genetic study on any hymenopteran group to date, we report an annotated secondary structural model for the D2, D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related apocritan taxa.     

Evaluation of the suitability of free-energy minimization using nearest-neighbor energy parameters for RNA secondary structure prediction

Background

A detailed understanding of an RNA's correct secondary and tertiary structure is crucial to understanding its function and mechanism in the cell. Free energy minimization with energy parameters based on the nearest-neighbor model and comparative analysis are the primary methods for predicting an RNA's secondary structure from its sequence. Version 3.1 of Mfold has been available since 1999. This version contains an expanded sequence dependence of energy parameters and the ability to incorporate coaxial stacking into free energy calculations. We test Mfold 3.1 by performing the largest and most phylogenetically diverse comparison of rRNA and tRNA structures predicted by comparative analysis and Mfold, and we use the results of our tests on 16S and 23S rRNA sequences to assess the improvement between Mfold 2.3 and Mfold 3.1.

Results

The average prediction accuracy for a 16S or 23S rRNA sequence with Mfold 3.1 is 41%, while the prediction accuracies for the majority of 16S and 23S rRNA structures tested are between 20% and 60%, with some having less than 20% prediction accuracy. The average prediction accuracy was 71% for 5S rRNA and 69% for tRNA. The majority of the 5S rRNA and tRNA sequences have prediction accuracies greater than 60%. The prediction accuracy of 16S rRNA base-pairs decreases exponentially as the number of nucleotides intervening between the 5' and 3' halves of the base-pair increases.

Conclusion

Our analysis indicates that the current set of nearest-neighbor energy parameters in conjunction with the Mfold folding algorithm are unable to consistently and reliably predict an RNA's correct secondary structure. For 16S or 23S rRNA structure prediction, Mfold 3.1 offers little improvement over Mfold 2.3. However, the nearest-neighbor energy parameters do work well for shorter RNA sequences such as tRNA or 5S rRNA, or for larger rRNAs when the contact distance between the base-pairs is less than 100 nucleotides.     

Secondary structure of mouse 28S rRNA and general model for the folding of the large rRNA in eukaryotes.

We present a secondary structure model for the entire sequence of mouse 28S rRNA (1) which is based on an extensive comparative analysis of the available eukaryotic sequences, i.e. yeast (2, 3), Physarum polycephalum (4), Xenopus laevis (5) and rat (6). It has been derived with close reference to the models previously proposed for yeast 26S rRNA (2) and for prokaryotic 23S rRNA (7-9). Examination of the recently published eukaryotic sequences confirms that all pro- and eukaryotic large rRNAs share a largely conserved secondary structure core, as already apparent from the previous analysis of yeast 26S rRNA (2). These new comparative data confirm most features of the yeast model (2). They also provide the basis for a few modifications and for new proposals which extend the boundaries of the common structural core (now representing about 85% of E. coli 23S rRNA length) and bring new insights for tracing the structural evolution, in higher eukaryotes, of the domains which have no prokaryotic equivalent and are inserted at specific locations within the common structural core of the large subunit rRNA.     

Aligned 18S and insect phylogeny

The nuclear small subunit rRNA (18S) has played a dominant role in the estimation of relationships among insect orders from molecular data. In previous studies, 18S sequences have been aligned by unadjusted automated approaches (computer alignments that are not manually readjusted), most recently with direct optimization (simultaneous alignment and tree building using a program called "POY"). Parsimony has been the principal optimality criterion. Given the problems associated with the alignment of rRNA, and the recent availability of the doublet model for the analysis of covarying sites using Bayesian MCMC analysis, a different approach is called for in the analysis of these data. In this paper, nucleotide sequence data from the 18S small subunit rRNA gene of insects are aligned manually with reference to secondary structure, and analyzed under Bayesian phylogenetic methods with both GTR+I+G and doublet models in MrBayes. A credible phylogeny of Insecta is recovered that is independent of the morphological data and (unlike many other analyses of 18S in insects) not contradictory to traditional ideas of insect ordinal relationships based on morphology. Hexapoda, including Collembola, are monophyletic. Paraneoptera are the sister taxon to a monophyletic Holometabola but weakly supported. Ephemeroptera are supported as the sister taxon of Neoptera, and this result is interpreted with respect to the evolution of direct sperm transfer and the evolution of flight. Many other relationships are well-supported but several taxa remain problematic, e.g., there is virtually no support for relationships among orthopteroid orders. A website is made available that provides aligned 18S data in formats that include structural symbols and Nexus formats.     

Seasonal dynamics of bacterioplankton community structure in a eutrophic lake as determined by 5S rRNA analysis.

Community structure of bacterioplankton was studied during the major growth season for phytoplankton (April to October) in the epilimnion of a temperate eutrophic lake (Lake Plusssee, northern Germany) by using comparative 5S rRNA analysis. Estimates of the relative abundances of single taxonomic groups were made on the basis of the amounts of single 5S rRNA bands obtained after high-resolution electrophoresis of RNA directly from the bacterioplankton. Full-sequence analysis of single environmental 5S rRNAs enabled the identification of single taxonomic groups of bacteria. Comparison of partial 5S rRNA sequences allowed the detection of changes of single taxa over time. Overall, the whole bacterioplankton community showed two to eight abundant (>4% of the total 5S rRNA) taxa. A distinctive seasonal succession was observed in the taxonomic structure of this pelagic community. A rather-stable community structure, with seven to eight different taxonomic units, was observed beginning in April during the spring phytoplankton bloom. A strong reduction in this diversity occurred at the beginning of the clear-water phase (early May), when only two to four abundant taxa were observed, with one taxon dominating (up to 72% of the total 5S rRNA). The community structure during summer stagnation (June and July) was characterized by frequent changes of different dominating taxa. During late summer, a dinoflagellate bloom (Ceratium hirudinella) occurred, with Comamonas acidovorans (beta-subclass of the class Proteobacteria) becoming the dominant bacterial species (average abundance of 43% of the total 5S rRNA). Finally, the seasonal dynamics of the community structure of bacterioplankton were compared with the abundances of other major groups of the aquatic food web, such as phyto- and zooplankton, revealing that strong grazing pressure by zooplankton can reduce microbial diversity substantially in pelagic environments.     

Homology modeling revealed more than 20,000 rRNA internal transcribed spacer 2 (ITS2) secondary structures

Structural genomics meets phylogenetics and vice versa: Knowing rRNA secondary structures is a prerequisite for constructing rRNA alignments for inferring phylogenies, and inferring phylogenies is a precondition to understand the evolution of such rRNA secondary structures. Here, both scientific worlds go together. The rRNA internal transcribed spacer 2 (ITS2) region is a widely used phylogenetic marker. Because of its high variability at the sequence level, correct alignments have to take into account structural information. In this study, we examine the extent of the conservation in structure. We present (1) the homology modeled secondary structure of more than 20,000 ITS2 covering about 14,000 species; (2) a computational approach for homology modeling of rRNA structures, which additionally can be applied to other RNA families; and (3) a database providing about 25,000 ITS2 sequences with their associated secondary structures, a refined ITS2 specific general time reversible (GTR) substitution model, and a scoring matrix, available at http://its2.bioapps.biozentrum.uni-wuerzburg.de.     

An overview of the secondary structure of the V4 region of eukaryotic small-subunit ribosomal RNA.

The V4 region of the small subunit (18S) ribosomal RNA was examined in 72 different sequences representing a broad sample eukaryotic diversity. This domain is the most variable region of the 18S rRNA molecule and ranges in length from ca. 230 to over 500 bases. Based upon comparative analysis, secondary structural models were constructed for all sequences and the resulting generalized model shows that most organisms possess seven helices for this region. The protists and two insects show from one to as many as four helices in addition to the above seven. In this report, we summarize secondary structure information presented elsewhere for the V4 region, describe the general features for helical and apical regions, and identify signature sequences useful in helix identification. Our model generally agrees with other current concepts; however, we propose modifications or alternative structures for the start of the V4 region, the large protist inserts, and the sector that may possibly contain a pseudoknot.     

Changes in bacterioplankton community structure and activity with depth in a eutrophic lake as revealed by 5S rRNA analysis

The community structure of bacterioplankton was studied at different depths (0 to 25 m) of a temperate eutrophic lake (Lake Plusssee in northern Germany) by using comparative 5S rRNA analysis. The relative amounts of taxonomic groups were estimated from 5S rRNA bands separated by high-resolution electrophoresis. Comparison of partial 5S rRNA sequences enabled detection of changes in single taxa over space and during seasons. Overall, the bacterioplankton community was dominated by 3 to 14 abundant (>4% of the total 5S rRNA) taxa. In general, the number of 5S rRNA bands (i.e., the number of bacterial taxa) decreased with depth. In the fall, when thermal stratification and chemical stratification were much more pronounced than they were in the spring, the correlation between the depth layers and the community structure was more pronounced. Therefore, in the fall each layer had its own community structure; i.e., there were different community structures in the epilimnion, the metalimnion, and the hypolimnion. Only three 5S rRNA bands were detected in the hypolimnion during the fall, and one band accounted for about 70% of the total 5S rRNA. The sequences of individual 5S rRNA bands from the spring and fall were different for all size classes analyzed except two bands, one of which was identified as Comamonas acidivorans. In the overall analysis of the depth profiles, the diversity in the epilimnion contrasted with the reduced diversity of the bacterioplankton communities in the hypolimnion, and large differences occurred in the composition of the epilimnion at different seasons except for generalists like C. acidivorans.     

The 23S ribosomal RNA higher-order structure of Pseudomonas cepacia and other prokaryotes

A 23S ribosomal RNA gene of Pseudomonas cepacia has been cloned and sequenced. A general higher-order structure model based on earlier published models has been derived from comparative analysis of 23S-like rRNAs of eubacteria, archaebacteria, organelles and eukaryotes. Differences between the previous models were carefully analyzed and controversial regions evaluated. Moderately large insertions and deletions have been found at new points in the secondary structure. The analysis of 50 published as well as unpublished 23S rRNA sequences provide additional proof for six of the seven previously suggested tertiary interactions within the 23S rRNA. P. cepacia is the first representative of the beta subgroup of the Proteobacteria phylum whose 23S rRNA has been sequenced. A tree reflecting evolutionary relationships of prokaryotes was constructed. The topology of this tree is in good agreement with the 16S rRNA tree.     

A universal model for the secondary structure of 5.8S ribosomal RNA molecules, their contact sites with 28S ribosomal RNAs, and their prokaryotic equivalent

The phylogenetic approach (ref. 1) has been utilized in construction of a universal 5.8S rRNA secondary structure model, in which about 65% of the residues exist in paired structures. Conserved nucleotides primarily occupy unpaired regions. Multiple compensating base changes are demonstrated to be present in each of the five postulated helices, thereby forming a major basis for their proof. The results of chemical and enzymatic probing of 5.8S rRNAs (ref. 13, 32) are fully consistent with, and support, our model. This model differs in several ways from recently proposed 5.8S rRNA models (ref. 3, 4), which are discussed. Each of the helices in our model has been extended to the corresponding bacterial, chloroplast and mitochondrial sequences, which are demonstrated to be positionally conserved by alignment with their eukaryotic counterparts. This extension is also made for the base paired 5.8S/28S contact points, and their prokaryotic and organelle counterparts. The demonstrated identity of secondary structure in these diverse molecules strongly suggests that they perform equivalent functions in prokaryotic and eukaryotic ribosomes.