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.     

Comparative sequence analysis and patterns of covariation in RNA secondary structures

A novel method of RNA secondary structure prediction based on a comparison of nucleotide sequences is described. This method correctly predicts nearly all evolutionarily conserved secondary structures of five different RNAs: tRNA, 5S rRNA, bacterial ribonuclease P (RNase P) RNA, eukaryotic small subunit rRNA, and the 3' untranslated region (UTR) of the Drosophila bicoid (bcd) mRNA. Furthermore, covariations occurring in the helices of these conserved RNA structures are analyzed. Two physical parameters are found to be important determinants of the evolution of compensatory mutations: the length of a helix and the distance between base-pairing nucleotides. For the helices of bcd 3' UTR mRNA and RNase P RNA, a positive correlation between the rate of compensatory evolution and helix length is found. The analysis of Drosophila bcd 3' UTR mRNA further revealed that the rate of compensatory evolution decreases with the physical distance between base-pairing residues. This result is in qualitative agreement with Kimura's model of compensatory fitness interactions, which assumes that mutations occurring in RNA helices are individually deleterious but become neutral in appropriate combinations.     

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.     

A phylogeny of anisopterous dragonflies (Insecta,Odonata) using mtRNA genes and mixed nucleotide/doublet models

The application of mixed nucleotide/doublet substitution models has recently received attention in RNA‐based phylogenetics. Within a Bayesian approach, it was shown that mixed models outperformed analyses relying on simple nucleotide models. We analysed an mt RNA data set of dragonflies representing all major lineages of Anisoptera plus outgroups, using a mixed model in a Bayesian and parsimony (MP) approach. We used a published mt 16S rRNA secondary consensus structure model and inferred consensus models for the mt 12S rRNA and tRNA valine. Secondary structure information was used to set data partitions for paired and unpaired sites on which doublet or nucleotide models were applied, respectively. Several different doublet models are currently available of which we chose the most appropriate one by a Bayes factor test. The MP reconstructions relied on recoded data for paired sites in order to account for character covariance and an application of the ratchet strategy to find most parsimonious trees. Bayesian and parsimony reconstructions are partly differently resolved, indicating sensitivity of the reconstructions to model specification. Our analyses depict a tree in which the damselfly family Lestidae is sister group to a monophyletic clade Epiophlebia + Anisoptera, contradicting recent morphological and molecular work. In Bayesian analyses, we found a deep split between Libelluloidea and a clade ‘Aeshnoidea’ within Anisoptera largely congruent with Tillyard’s early ideas of anisopteran evolution, which had been based on evidently plesiomorphic character states. However, parsimony analysis did not support a clade ‘Aeshnoidea’, but instead, placed Gomphidae as sister taxon to Libelluloidea. Monophyly of Libelluloidea is only modestly supported, and many inter‐family relationships within Libelluloidea do not receive substantial support in Bayesian and parsimony analyses. We checked whether high Bayesian node support was inflated owing to either: (i) wrong secondary consensus structures; (ii) under‐sampling of the MCMC process, thereby missing other local maxima; or (iii) unrealistic prior assumptions on topologies or branch lengths. We found that different consensus structure models exert strong influence on the reconstruction, which demonstrates the importance of taxon‐specific realistic secondary structure models in RNA phylogenetics.     

Phylogeny of Sibling Species of Simulium venustrum and S. verecundum (Diptera: Simuliidae) Based on Sequences of the Mitochondrial 165 rRNA Gene

The phylogeny of four sibling species of Simulium venustum (CC, CC3, CC4, and AC(gB)) and two sibling species of S. verecundum (AA and ACD = S. rostratum) was reconstructed using nucleotide sequences of the mitochondrial large subunit rRNA gene. Separate phylogenetic analyses were performed by dividing the sequence data into: (1) helices of a computer-generated secondary structure of the rDNA; (2) loops of the computer-generated structure; (3) helices of a consensus secondary structure (deduced by comparing the computer-generated structure of black flies with the corresponding structures proposed for the fruit fly (Drosophila yakuba) and the mosquito (Aedes albopictus)); (4) loops of the consensus structure; (5) both helices and loops of the consensus secondary structure; and (6) the entire sequence regardless of secondary structure, including 11 variable sites in regions where the prediction of secondary structure was not possible. We found that different data sets led to different phylogenetic conclusions. The phylogenies based on data sets 4 and 6 were consistent with nonmolecular evidence, while the phylogenies based on other data sets were not. Our study suggests S. decorum, a morphospecies, might have shared a common ancestor with sibling species of S. venustum.     

A consensus prediction of the secondary structure for the 6-phospho-β-D-galactosidase superfamily

Two separate unrefined models for the secondary structure of two subfamilies of the 6-phospho-β-D -galactosidase superfamily were independently constructed by examining patterns of variation and conservation within homologous protein sequences, assigning surface, interior, parsing, and active site residues to positions in the alignment, and identifying periodicities in these. A consensus model for the secondary structure of the entire superfamily was then built. The prediction tests the limits of an unrefined prediction made using this approach in a large protein with substantial functional and sequence divergence within the family. The protein belongs to the (α–β class), with the core β strands aligned parallel. The supersecondary structural elements that are readily identified in this model is a parallel β sheet built by strands C, D, and E, with helices 2 and 3 connecting strands (C + D) and (D + E), respectively, and an analogous α–β unit (strand G and helix 7) toward the end of the sequence. The resemblance of the supersecondary model to the tertiary structure formed by 8-fold α–β barrel proteins is almost certainly not coincidental. © 1995 Wiley-Liss, Inc.     

The complete nucleotide sequence of mouse 28S rRNA gene. Implications for the process of size increase of the large subunit rRNA in higher eukaryotes.

We have determined the complete nucleotide sequence (4712 nucleotides) of the mouse 28S rRNA gene. Comparison with all other homologs indicates that the potential for major variations in size during the evolution has been restricted to a unique set of a few sites within a largely conserved secondary structure core. The D (divergent) domains, responsible for the large increase in size of the molecule from procaryotes to higher eukaryotes, represent half the mouse 28S rRNA length. They show a clear potential to form self-contained secondary structures. Their high GC content in vertebrates is correlated with the folding of very long stable stems. Their comparison with the two other vertebrates, xenopus and rat, reveals an history of repeated insertions and deletions. During the evolution of vertebrates, insertion or deletion of new sequence tracts preferentially takes place in the subareas of D domains where the more recently fixed insertions/deletions were located in the ancestor sequence. These D domains appear closely related to the transcribed spacers of rRNA precursor but a sizable fraction displays a much slower rate of sequence variation.     

RNA secondary structure and compensatory evolution

The classic concept of epistatic fitness interactions between genes has been extended to study interactions within gene regions, especially between nucleotides that are important in maintaining pre-mRNA/mRNA secondary structures. It is shown that the majority of linkage disequilibria found within the Drosophila Adh gene are likely to be caused by epistatic selection operating on RNA secondary structures. A recently proposed method of RNA secondary structure prediction based on DNA sequence comparisons is reviewed and applied to several types of RNAs, including tRNA, rRNA, and mRNA. The patterns of covariation in these RNAs are analyzed based on Kimura's compensatory evolution model. The results suggest that this model describes the substitution process in the pairing regions (helices) of RNA secondary structures well when the helices are evolutionarily conserved and thermodynamically stable, but fails in some other cases. Epistatic selection maintaining pre-mRNA/mRNA secondary structures is compared to weak selective forces that determine features such as base composition and synonymous codon usage. The relationships among these forces and their relative strengths are addressed. Finally, our mutagenesis experiments using the Drosophila Adh locus are reviewed. These experiments analyze long-range compensatory interactions between the 5' and 3' ends of Adh mRNA, the different constraints on secondary structures in introns and exons, and the possible role of secondary structures in RNA splicing.     

Epistasis Creates Invariant Sites and Modulates the Rate of Molecular Evolution

Invariant sites are a common feature of amino acid sequence evolution. The presence of invariant sites is frequently attributed to the need to preserve function through site-specific conservation of amino acid residues. Amino acid substitution models without a provision for invariant sites often fit the data significantly worse than those that allow for an excess of invariant sites beyond those predicted by models that only incorporate rate variation among sites (e.g., a Gamma distribution). An alternative is epistasis between sites to preserve residue interactions that can create invariant sites. Through computer-simulated sequence evolution, we evaluated the relative effects of site-specific preferences and site-site couplings in the generation of invariant sites and the modulation of the rate of molecular evolution. In an analysis of ten major families of protein domains with diverse sequence and functional properties, we find that the negative selection imposed by epistasis creates many more invariant sites than site-specific residue preferences alone. Further, epistasis plays an increasingly larger role in creating invariant sites over longer evolutionary periods. Epistasis also dictates rates of domain evolution over time by exerting significant additional purifying selection to preserve site couplings. These patterns illuminate the mechanistic role of epistasis in the processes underlying observed site invariance and evolutionary rates.     

A new sequence data set of SSU rRNA gene for Scleractinia and its phylogenetic and ecological applications

Scleractinian corals (i.e. hard corals) play a fundamental role in building and maintaining coral reefs, one of the most diverse ecosystems on Earth. Nevertheless, their phylogenies remain largely unresolved and little is known about dispersal and survival of their planktonic larval phase. The small subunit ribosomal RNA (SSU rRNA) is a commonly used gene for DNA barcoding in several metazoans, and small variable regions of SSU rRNA are widely adopted as barcode marker to investigate marine plankton community structure worldwide. Here, we provide a large sequence data set of the complete SSU rRNA gene from 298 specimens, representing all known extant reef coral families and a total of 106 genera. The secondary structure was extremely conserved within the order with few exceptions due to insertions or deletions occurring in the variable regions. Remarkable differences in SSU rRNA length and base composition were detected between and within acroporids (Acropora, Montipora, Isopora and Alveopora) compared to other corals. The V4 and V9 regions seem to be promising barcode loci because variation at commonly used barcode primer binding sites was extremely low, while their levels of divergence allowed families and genera to be distinguished. A time‐calibrated phylogeny of Scleractinia is provided, and mutation rate heterogeneity is demonstrated across main lineages. The use of this data set as a valuable reference for investigating aspects of ecology, biology, molecular taxonomy and evolution of scleractinian corals is discussed.     

Structural analysis of three prokaryotic 5S rRNA species and selected 5S rRNA--ribosomal-protein complexes by means of Pb(II)-induced hydrolysis.

Lead ions have been applied to the structural analysis of 5S rRNA from Thermus thermophilus, Bacillus stearothermophilus and Escherichia coli. Based on the distribution of Pb(II)-induced cleavages, some minor modifications of the consensus secondary structure model of 5S rRNA are proposed. They include the possible base pairing between nucleotides at positions 11 and 109, as well as changes in secondary interactions within the helix B region. The 'prokaryotic arm' region is completely resistant to hydrolysis in the three RNA species, suggesting that it is a relatively stable, highly ordered structure. Hydrolysis of E. coli 5S rRNA complexed with ribosomal protein L18 shows, besides the shielding effect of the bound protein, a highly enhanced cleavage between A108 and A109. It supports the concept that the major L18-induced conformational change involves the junction of helices A, B and D.     

Tracing the evolution of RNA structure in ribosomes

The elucidation of ribosomal structure has shown that the function of ribosomes is fundamentally confined to dynamic interactions established between the RNA components of the ribosomal ensemble. These findings now enable a detailed analysis of the evolution of ribosomal RNA (rRNA) structure. The origin and diversification of rRNA was studied here using phylogenetic tools directly at the structural level. A rooted universal tree was reconstructed from the combined secondary structures of large (LSU) and small (SSU) subunit rRNA using cladistic methods and considerations in statistical mechanics. The evolution of the complete repertoire of structural ribosomal characters was formally traced lineage-by-lineage in the tree, showing a tendency towards molecular simplification and a homogeneous reduction of ribosomal structural change with time. Character tracing revealed patterns of evolution in inter-subunit bridge contacts and tRNA-binding sites that were consistent with the proposed coupling of tRNA translocation and subunit movement. These patterns support the concerted evolution of tRNA-binding sites in the two subunits and the ancestral nature and common origin of certain structural ribosomal features, such as the peptidyl (P) site, the functional relay of the penultimate stem helix of SSU rRNA, and other structures participating in ribosomal dynamics. Overall results provide a rare insight into the evolution of ribosomal structure.     

The origin and evolution of variable-region helices in V4 and V7 of the small-subunit ribosomal RNA of branchiopod crustaceans

We sequenced the V4 and V7 regions of the small-subunit ribosomal RNA (SSU rRNA) from 38 species of branchiopod crustaceans (e.g., Artemia, Daphnia, Triops) representing all eight extant orders. Ancestral large-bodied taxa in the orders Anostraca, Notostraca, Laevicaudata, and Spinicaudata (limnadiids and cyzicids) possess the typical secondary structure in these regions, whereas the spinicaudatan Cyclestheria and all of the cladocerans (Anomopoda, Ctenopoda, Onychopoda, and Haplopoda) possess three unique helices. Although the lengths and primary sequences of the distal ends of these helices are extremely variable, their locations, secondary structures, and primary sequences at the proximal end are conserved, indicating that they are homologous. This evidence supports the classical view that Cladocera is a monophyletic group and that the cyclestheriids are transitional between spinicaudatans and cladocerans. The single origin and persistence since the Permian of the unique cladoceran helices suggests that births and deaths of variable region helices have been rare. The broad range of sequence divergences observed among the cladoceran helices permitted us to make inferences about their evolution. Although their proximal ends are very GC-biased, there is a significant negative correlation between length and GC content due to an increasing proportion of U at their distal ends. Slippage-like processes occurring at unpaired nucleotides or bulges, which are very U-biased, are associated with both helix origin and runaway length expansion. The overall GC contents and lengths of V4 and V7 are highly correlated. More surprisingly, the lengths of these SSU rRNA variable regions are also highly correlated with the length of the large-subunit rRNA expansion segment, D2, indicating that mechanisms affecting length variation do so both across single genes and across genes in the rRNA gene family.     

Secondary structure and phylogenetic analysis of the internal transcribed spacers 1 and 2 of bush crickets (Orthoptera: Tettigoniidae: Barbitistini)

We inferred secondary structure models of the internal transcribed spacers (ITS) 1 and 2 of bush crickets using a combined comparative and thermodynamic approach. The inferred secondary structure models were used to account for interdependency of interacting nucleotides in a phylogenetic analysis of the bush cricket genus Poecilimon. Our analysis indicates that the two previously reported conformational structures (i.e., hairpin and ring) of ITS2 are likely to fold in bush crickets as well and that both predicted structures are similar to those proposed for other eukaryotes. Comparing predicted ITS1 secondary structure models proved to be difficult because of substantial variation in their nucleotide sequence length. Our study revealed that the phylogenetic signal of ITS1 and ITS2 is largely congruent with that preserved in the mitochondrial genes 16S rRNA, tRNA‐Val and 12S rRNA. The phylogenetic signal in both the nuclear and the mitochondrial genome question the monophyly of the genus Poecilimon: species of the genera Poecilimonella, Parapoecilimon, Polysarcus and Phonochorion consistently cluster within Poecilimon.     

Evolution of Hypervariable Regions, V4 and V7, of Insect 18S rRNA and Their Phylogenetic Implications

We compared primary and secondary structures of V4 (helices E23-2 to E23-5) and V7 (helix 43) regions of 18S rRNAs in insects and the other three major arthropod groups (crustaceans, myriapods, and chelicerates) known so far. We found that the lengths of primary sequences and the shapes of secondary structures of these two hypervariable regions of insect 18S rRNA even at infraclass levels are phylogenetically informative and reflect major steps in insect evolution. The long sequence insertion and bifurcated shape of helices E23-2 to E23-5 in the V4 region are unique synapomorphic characters for winged insects (Pterygota). The long sequence insertion and expanded stem length of helix 43 in the V7 region are synapomorphic characters for holometabolous insects which conduct complete metamorphosis. The strongly conserved secondary structures suggest the possibility that these hypervariable regions may be related with certain important cellular functions unknown thus far. The comparison with insect fossil records revealed that the pterygote synapomorphy (V4) and the holometabolous synapomorphy (V7) were established prior to the acquisition of insect wings (flight system) and prior to the development of complete metamorphosis, respectively. These synapomorphies have been also relatively stable over at least 300 Myr and 280 Myr, respectively as well. It implies that the expansion events of the V4 and V7 regions have not occurred simultaneously but independently at different periods during the insect evolution. Then this suggests that V4 and V7 regions are not functionally correlated as recently suggested by Crease and Coulbourn.     

Secondary structure models of the nuclear internal transcribed spacer regions and 5.8S rRNA in Calciodinelloideae (Peridiniaceae) and other dinoflagellates

Secondary structure models of the 5.8S rRNA and both internal transcribed spacers (ITS1 and ITS2) are proposed for Calciodinelloideae (Peridiniaceae) and are also plausible for other dinoflagellates. The secondary structure of the 5.8S rRNA corresponds to previously developed models, with two internal paired regions and at least one 5.8S rRNA–28S rRNA interaction. A general secondary structure model of ITS1 for Calciodinelloideae (and other dinoflagellates), consisting of an open multibranch loop with three major helices, is proposed. The homology of these paired regions with those found in other taxa, published in previous studies (e.g. yeast, green algae and Platyhelmithes) remains to be determined. Finally, a general secondary structure model of ITS2 for Calciodinelloideae (and other dinoflagellates) is reconstructed. Based on the 5.8S rRNA–28S rRNA interaction, it consists of a closed multibranch loop, with four major helices. At least helix III and IV have homology with paired regions found in other eukaryotic taxa (e.g. yeast, green algae and vertebrates). Since the secondary structures of both ITS regions are more conserved than the nucleotide sequences, their analysis helps in understanding molecular evolution and increases the number of structural characters. Thus, the structure models developed in this study may be generally useful for future phylogenetic analyses.     

Coevolution in RNA Molecules Driven by Selective Constraints: Evidence from 5S rRNA

Understanding intra-molecular coevolution helps to elucidate various structural and functional constraints acting on molecules and might have practical applications in predicting molecular structure and interactions. In this study, we used 5S rRNA as a template to investigate how selective constraints have shaped the RNA evolution. We have observed the nonrandom occurrence of paired differences along the phylogenetic trees, the high rate of compensatory evolution, and the high TIR scores (the ratio of the numbers of terminal to intermediate states), all of which indicate that significant positive selection has driven the evolution of 5S rRNA. We found three mechanisms of compensatory evolution: Watson-Crick interaction (the primary one), complex interactions between multiple sites within a stem, and interplay of stems and loops. Coevolutionary interactions between sites were observed to be highly dependent on the structural and functional environment in which they occurred. Coevolution occurred mostly in those sites closest to loops or bulges within structurally or functionally important helices, which may be under weaker selective constraints than other stem positions. Breaking these pairs would directly increase the size of the adjoining loop or bulge, causing a partial or total structural rearrangement. In conclusion, our results indicate that sequence coevolution is a direct result of maintaining optimal structural and functional integrity.     

Secondary structure of the central region of bacteriophage MS2 RNA. Conservation and biological significance

The RNA of the Escherichia coli RNA phages is highly structured with 75% of the nucleotides estimated to take part in base-pairing. We have used enzymatic and chemical sensitivity of nucleotides, phylogenetic sequence comparison and the phenotypes of constructed mutants to develop a secondary structure model for the central region (900 nucleotides) of the group I phage MS2. The RNA folds into a number of, mostly irregular, helices and is further condensed by several long-distance interactions. There is substantial conservation of helices between the related groups I and II, attesting to the relevance of discrete RNA folding. In general, the secondary structure is thought to be needed to prevent annealing of plus and minus strand and to confer protection against RNase. Superimposed, however, are features required to regulate translation and replication. The MS2 RNA section studied here contains three translational start sites, as well as the binding sites for the coat protein and the replicase enzyme. Considering the density of helices along the RNA, it is not unexpected to find that all these sites lie in helical regions. This fact, however, does not mean that these sites are recognized as secondary structure elements by their interaction partners. This holds true only for the coat protein binding site. The other four sites function in the unfolded state and the stability of the helix in which they are contained serves to negatively control their accessibility. Mutations that stabilize helices containing ribosomal binding sites reduce their efficiency and vice versa. Comparison of homologous helices in different phage RNAs indicates that base substitutions have occurred in such a way that the thermodynamic stability of the helix is maintained. The evolution of individual helices shows several distinct size-reduction patterns. We have observed codon deletions from loop areas and shortening of hairpins by base-pair deletions from either the bottom, the middle or the top of stem structures. Evidence for the coaxial stacking of some helical segments is discussed.