Evolved RNA Secondary Structure and the Rooting of the Universal Tree of Life

The origin and diversification of RNA secondary structure were traced using cladistic methods. Structural components were coded as polarized and ordered multi-state characters, following a model of character state transformation outlined by considerations in statistical mechanics. Several classes of functional RNA were analyzed, including ribosomal RNA (rRNA). Considerable phylogenetic signal was present in their secondary structure. The intrinsically rooted phylogenies reconstructed from evolved RNA structure depicted those derived from nucleic acid sequence at all taxonomical levels, and grouped organisms in concordance with traditional classification, especially in the archaeal and eukaryal domains. Natural selection appears therefore to operate early in the information flow that originates in sequence and ends in an adapted phenotype. When examining the hierarchical classification of the living world, phylogenetic analysis of secondary structure of the small and large rRNA subunits reconstructed a universal tree of life that branched in three monophyletic groups corresponding to Eucarya, Archaea, and Bacteria, and was rooted in the eukaryotic branch. Ribosomal characters involved in the translational cycle could be easily traced and showed that transfer RNA (tRNA) binding domains in the large rRNA subunit evolved concurrently with the rest of the rRNA molecule. Results suggest it is equally parsimonious to consider that ancestral unicellular eukaryotes or prokaryotes gave rise to all extant life forms and provide a rare insight into the early evolution of nucleic acid and protein biosynthesis. Received: 13 September 2000 / Accepted: 27 August 2001     

Use of RNA Secondary Structure for Studying the Evolution of RNase P and RNase MRP

Secondary structure is evaluated for determining evolutionary relationships between catalytic RNA molecules that are so distantly related they are scarcely alignable. The ribonucleoproteins RNase P (P) and RNase MRP (MRP) have been suggested to be evolutionarily related because of similarities in both function and secondary structure. However, their RNA sequences cannot be aligned with any confidence, and this leads to uncertainty in any trees inferred from sequences. We report several approaches to using secondary structures for inferring evolutionary trees and emphasize quantitative tests to demonstrate that evolutionary information can be recovered. For P and MRP, three hypotheses for the relatedness are considered. The first is that MRP is derived from P in early eukaryotes. The next is that MRP is derived from P from an early endosymbiont. The third is that both P and MRP evolved in the RNA-world (and the need for MRP has since been lost in prokaryotes). Quantitative comparisons of the pRNA and mrpRNA secondary structures have found that the possibility of an organellar origin of MRP is unlikely. In addition, comparison of secondary structures support the identity of an RNase P–like sequence in the maize chloroplast genome. Overall, it is concluded that RNA secondary structure is useful for evaluating evolutionary relatedness, even with sequences that cannot be aligned with confidence. Received: 19 July 1999 / Accepted: 3 May 2000     

Primary and Secondary Structure of the Small-Subunit Ribosomal RNA of the Naked, Marine Amoeba Vannella anglica: Phylogenetic Implications

The primary and secondary structure of the small-subunit ribosomal RNA (ssrRNA) gene from the naked, marine amoeba, Vannella anglica (subclass Gymnamoebia), was determined. The ssrRNA is 1962 nucleotides in length, with a low G+C content of 37.1%. The ssrRNA is composed of several uncommon secondary structure features including helix E8-1, which may be a useful target for rRNA probes for the direct identification of isolates in mixed culture. Phylogenetic analysis of sequence data showed that V. anglica branched prior to the rapid diversification of the eukaryotes. It did not associate with the other naked, lobose amoebae represented by Acanthamoeba and Hartmannella, indicating that Vannella represents a separate amoeboid lineage and the subclass Gymnamoebia is polyphyletic. Received: 9 July 1998 / Accepted: 16 November 1998     

Selective Pressures on RNA Hairpins In Vivo and In Vitro

Comparison of the most stable potential hairpins in the sequences of natural ribozymes with those in the randomized sequences has revealed that the hairpin loop energies are lower than expected by chance. Although these hairpins are not necessarily parts of functional structures, there is a selective pressure to diminish the destabilizing free energies of the hairpin loops. In contrast, no significant bias is observed in the stacking values of the most stable stems. In the ribozymes isolated in vitro the loops of potential hairpins are closer to random values, which can result in less efficient folding rates. Furthermore, the effects of kinetic traps seem to be more significant in the folding pathways of the in vitro isolates due to a potential to form stable stacks incompatible with the functional folds. Similarly to natural ribozyme sequences, the untranslated regions of viral RNAs also form hairpins with relatively low loop free energies. These evolutionary trends suggest ways for efficient engineering of improved RNA constructs on the basis of analysis of in vitro isolates and approaches for the search of regions coding for functional RNA structures in large genome sequences. Received: 12 January 2001 / Accepted: 21 May 2001     

Analysis of the Primary Sequence and Secondary Structure of the Unusually Long SSU rRNA of the Soil Bug, Armadillidium vulgare

The complete nucleotide sequence of the SSU rRNA gene from the soil bug, Armadillidium vulgare (Crustacea, Isopoda), was determined. It is 3214 bp long, with a GC content of 56.3%. It is not only the longest SSU rRNA gene among Crustacea but also longer than any other SSU rRNA gene except that of the strepsipteran insect, Xenos vesparum (3316 bp). The unusually long sequence of this species is explained by the long sequences of variable regions V4 and V7, which make up more than half of the total length. RT-PCR analysis of these two regions showed that the long sequences also exist in the mature rRNA and sequence simplicity analysis revealed the presence of slippage motifs in these two regions. The putative secondary structure of the rRNA is typical for eukaryotes except for the length and shape variations of the V2, V4, V7, and V9 regions. Each of the V2, V4, and V7 regions was elongated, while the V9 region was shortened. In V2, two bulges, located between helix 8 and helix 9 and between helix 9 and helix 10, were elongated. In V4, stem E23-3 was dramatically expanded, with several small branched stems. In V7, stem 43 was branched and expanded. Comparisons with the unusually long SSU rRNAs of other organisms imply that the increase in total length of SSU rRNA is due mainly to expansion in the V4 and V7 regions. Received: 2 March 1999 / Accepted: 22 July 1999     

RNA Folding Argues Against a Hot-Start Origin of Life

Opinion is strongly divided on whether life arose on earth under hot or cold conditions, the hot-start and cold-start scenarios, respectively. The origin of life close to deep thermal vents appears as the majority opinion among biologists, but there is considerable biochemical evidence that high temperatures are incompatible with an RNA world. To be functional, RNA has to fold into a three-dimensional structure. We report both theoretical and experimental results on RNA folding and show that (as expected) hot conditions strongly reduce RNA folding. The theoretical results come from energy-minimization calculations of the average extent of folding of RNA, mainly from 0–90°C, for both random sequences and tRNA sequences. The experimental results are from circular-dichroism measurements of tRNA over a similar range of temperatures. The quantitative agreement between calculations and experiment is remarkable, even to the shape of the curves indicating the cooperative nature of RNA folding and unfolding. These results provide additional evidence for a lower temperature stage being necessary in the origin of life. Received: 1 March 2000 / Accepted: 14 June 2000     

Trematode and Monogenean rRNA ITS2 Secondary Structures Support a Four-Domain Model

The secondary structure of rRNA internal transcribed spacer 2 is important in the process of ribosomal biogenesis. Trematode ITS sequences are poorly conserved and difficult to align for phylogenetic comparisons above a family level. If a conserved secondary structure can be identified, it can be used to guide primary sequence alignments. ITS2 sequences from 39 species were compared. These species span four orders of trematodes (Echinostomiformes, Plagiorchiformes, Strigeiformes, and Paramphistomiformes) and one monogenean (Gyrodactyliformes). The sequences vary in length from 251 to 431 bases, with an average GC content of 48%. The monogenean sequence could not be aligned with confidence to the trematodes. Above the family level trematode sequences were alignable from the 5′ end for 139 bases. Secondary structure foldings predicted a four-domain model. Three folding patterns were required for the apex of domain B. The folding pattern of domains C and D varies for each family. The structures display a high GC content within stems. Bases A and U are favored in unpaired regions and variable sites cluster. This produces a mosaic of conserved and variable regions with a structural conformation resistant to change. Two conserved strings were identified, one in domain B and the other in domain C. The first site can be aligned to a processing site identified in yeast and rat. The second site has been found in plants, and structural location appears to be important. A phylogenetic tree of the trematode sequences, aligned with the aid of secondary structures, distinguishes the four recognized orders. Received: 21 November 1997 / Accepted: 9 February 1998     

Relics from the RNA World

An RNA world is widely accepted as a probable stage in the early evolution of life. Two implications are that proteins have gradually replaced RNA as the main biological catalysts and that RNA has not taken on any major de novo catalytic function after the evolution of protein synthesis, that is, there is an essentially irreversible series of steps RNA → RNP → protein. This transition, as expected from a consideration of catalytic perfection, is essentially complete for reactions when the substrates are small molecules. Based on these principles we derive criteria for identifying RNAs in modern organisms that are relics from the RNA world and then examine the function and phylogenetic distribution of RNA for such remnants of the RNA world. This allows an estimate of the minimum complexity of the last ribo-organism—the stage just preceding the advent of genetically encoded protein synthesis. Despite the constraints placed on its size by a low fidelity of replication (the Eigen limit), we conclude that the genome of this organism reached a considerable level of complexity that included several RNA-processing steps. It would include a large protoribosome with many smaller RNAs involved in its assembly, pre-tRNAs and tRNA processing, an ability for recombination of RNA, some RNA editing, an ability to copy to the end of each RNA strand, and some transport functions. It is harder to recognize specific metabolic reactions that must have existed but synthetic and bio-energetic functions would be necessary. Overall, this requires that such an organism maintained a multiple copy, double-stranded linear RNA genome capable of recombination and splicing. The genome was most likely fragmented, allowing each ``chromosome' to be replicated with minimum error, that is, within the Eigen limit. The model as developed serves as an outgroup to root the tree of life and is an alternative to using sequence data for inferring properties of the earliest cells. Received: 14 January 1997 / Accepted: 19 May 1997     

The Prebiotic Synthesis of Modified Purines and Their Potential Role in the RNA World

Modified purines are found in all organisms in the tRNA, rRNA, and even DNA, raising the possibility of an early role for these compounds in the evolution of life. These include N ⁶-methyladenine, 1-methyladenine, N ⁶,N ⁶-dimethyladenine, 1-methylhypoxanthine, 1-methylguanine, and N ²-methylguanine. We find that these bases as well as a number of nonbiological modified purines can be synthesized from adenine and guanine by the simple reaction of an amine or an amino group with adenine and guanine under the concentrated conditions of the drying-lagoon or drying-beach model of prebiotic synthesis with yields as high as 50%. These compounds are therefore as prebiotic as adenine and guanine and could have played an important role in the RNA world by providing additional functional groups in ribozymes, especially for the construction of hydrophobic binding pockets. Received: 7 August 1998 / Accepted: 31 December 1998     

Amino Acids as RNA Ligands: A Direct-RNA-Template Theory for the Code's Origin

Numerous RNA binding sites for specific amino acids are now known, coming predominantly from selection-amplification experiments. These sites are chemically discriminating despite being predominantly small, simple RNA structures: internal and bulge loops. Recent studies of sites for hydrophobic side chains suggest that there are other generalizable structural features which recur in hydrophobic RNA sites. Further, sites for hydrophobic side chains can contain codons for the bound amino acid, as has also long been known for the polar amino acid arginine. Such findings are comprehensively reviewed, and the implications for the origin of coded peptide synthesis are considered. An origins hypothesis which accommodates all the data, DRT (direct RNA templating), is formulated. Received: 22 December 1997 / Accepted: 13 February 1998     

Contribution of Mutation and RNA Recombination to the Evolution of a Plant Pathogenic RNA

The nucleotide sequence of 17 variants of the satellite RNA of cucumber mosaic virus (CMV-satRNA) isolated from field-infected tomato plants in the springs of 1989, 1990, and 1991 was determined. The sequence of each of the 17 satRNAs was unique and was between 334 and 340 nucleotides in length; 57 positions were polymorphic. There was much genetic divergence, ranging from 0.006 to 0.141 nucleotide substitutions per site for pairwise comparisons, and averaging 0.074 for any pair. When the polymorphic positions were analyzed relative to a secondary structure model proposed for CMV-satRNAs, it was found that there were significantly different numbers of changes in base-paired and non–base-paired positions, and that mutations that did not disrupt base pairing were preferred at the putatively paired sites. This supports the concept that the need to maintain a functional structure may limit genetic divergence of CMV-satRNA. Phylogenetic analyses showed that the 17 CMV-satRNA variants clustered into two subgroups, I and II, and evolutionary lines proceeding by the sequential accumulation of mutations were apparent. Three satRNA variants were outliers for these two phylogenetic groups. They were shown to be recombinants of subgroup I and II satRNAs by calculating phylogenies for different molecular regions and by using Sawyer's test for gene conversion. At least two recombination events were required to produce these three recombinant satRNAs. Thus, recombinants were found to be frequent (∼17%) in natural populations of CMV-satRNA, and recombination may make an important contribution to the generation of new variants. To our knowledge this is the first report of data allowing the frequency of recombinant isolates in natural populations of an RNA replicon to be estimated. Received: 14 May 1996 / Accepted: 17 July 1996     

Conserved Gene Clusters in Bacterial Genomes Provide Further Support for the Primacy of RNA

Five complete bacterial genome sequences have been released to the scientific community. These include four (eu)Bacteria, Haemophilus influenzae, Mycoplasma genitalium, M. pneumoniae, and Synechocystis PCC 6803, as well as one Archaeon, Methanococcus jannaschii. Features of organization shared by these genomes are likely to have arisen very early in the history of the bacteria and thus can be expected to provide further insight into the nature of early ancestors. Results of a genome comparison of these five organisms confirm earlier observations that gene order is remarkably unpreserved. There are, nevertheless, at least 16 clusters of two or more genes whose order remains the same among the four (eu)Bacteria and these are presumed to reflect conserved elements of coordinated gene expression that require gene proximity. Eight of these gene orders are essentially conserved in the Archaea as well. Many of these clusters are known to be regulated by RNA-level mechanisms in Escherichia coli, which supports the earlier suggestion that this type of regulation of gene expression may have arisen very early. We conclude that although the last common ancestor may have had a DNA genome, it likely was preceded by progenotes with an RNA genome. Received: 10 March 1996 / Accepted: 20 May 1997     

Synonymous Codon Choices in the Extremely GC-Poor Genome of Plasmodium falciparum: Compositional Constraints and Translational Selection

We have analyzed the patterns of synonymous codon preferences of the nuclear genes of Plasmodium falciparum, a unicellular parasite characterized by an extremely GC-poor genome. When all genes are considered, codon usage is strongly biased toward A and T in third codon positions, as expected, but multivariate statistical analysis detects a major trend among genes. At one end genes display codon choices determined mainly by the extreme genome composition of this parasite, and very probably their expression level is low. At the other end a few genes exhibit an increased relative usage of a particular subset of codons, many of which are C-ending. Since the majority of these few genes is putatively highly expressed, we postulate that the increased C-ending codons are translationally optimal. In conclusion, while codon usage of the majority of P. falciparum genes is determined mainly by compositional constraints, a small number of genes exhibit translational selection. Received: 10 November 1998 / Accepted: 28 January 1999     

Plant Mitochondrial RNA Editing

RNA editing affects messenger RNAs and transfer RNAs in plant mitochondria by site-specific exchange of cytidine and uridine bases in both seed and nonseed plants. Distribution of the phenomenon among bryophytes has been unclear since RNA editing has been detected in some but not all liverworts and mosses. A more detailed understanding of RNA editing in plants required extended data sets for taxa and sequences investigated. Toward this aim an internal region of the mitochondrial nad5 gene (1104 nt) was analyzed in a large collection of bryophytes and green algae (Charales). The genomic nad5 sequences predict editing in 30 mosses, 2 hornworts, and 7 simple thalloid and leafy liverworts (Jungermanniidae). No editing is, however, required in seven species of the complex thalloid liverworts (Marchantiidae) and the algae. RNA editing among the Jungermanniidae, on the other hand, reaches frequencies of up to 6% of codons being modified. Predictability of RNA editing from the genomic sequences was confirmed by cDNA analysis in the mosses Schistostega pennata and Rhodobryum roseum, the hornworts Anthoceros husnotii and A. punctatus, and the liverworts Metzgeria conjugata and Moerckia flotoviana. All C-to-U nucleotide exchanges predicted to reestablish conserved codons were confirmed. Editing in the hornworts includes the removal of genomic stop codons by frequent reverse U-to-C edits. Expectedly, no RNA editing events were identified by cDNA analysis in the marchantiid liverworts Ricciocarpos natans, Corsinia coriandra, and Lunularia cruciata. The findings are discussed in relation to models on the phylogeny of land plants. Received: 2 April 1998 / Accepted: 4 August 1998     

Relationships Between Genomic G+C Content,RNA Secondary Structures,and Optimal Growth Temperature in Prokaryotes

G:C pairs are more stable than A:T pairs because they have an additional hydrogen bond. This has led to many studies on the correlation between the guanine+cytosine (G+C) content of nucleic acids and temperature over the last 20 years. We collected the optimal growth temperatures (T_opt) and the G+C contents of genomic DNA; 23S, 16S, and 5S ribosomal RNAs; and transfer RNAs for 764 prokaryotic species. No correlation was found between genomic G+C content and T_opt, but there were striking correlations between the G+C content of ribosomal and transfer RNA stems and T_opt. Two explanations have been proposed—neutral evolution and selection pressure—for the approximate equalities of G and C (respectively, A and T) contents within each strand of DNA molecules. Our results do not support the notion that selection pressure induces complementary oligonucleotides in close proximity and therefore numerous secondary structures in prokaryotic DNA, as the genomic G+C content does not behave in the same way as that of folded RNA with respect to optimal growth temperature. Received: 25 September 1996 / Accepted: 21 January 1997     

The Vertebrate 7S K RNA Separates Hagfish (Myxine glutinosa) and Lamprey (Lampetra fluviatilis)

7S RNA sequences from the hagfish (Myxiniformes) and lamprey (Petromyzontiformes) were cloned and analyzed. In both species, 7S L RNA (also designated SRP RNA, since it represents the RNA constituent of the signal recognition particle) was clearly detectable. The sequence similarity between the two species was 86%, compared with about 75% similarity between either of these species and mammals. 7S K RNA was also cloned from the lamprey. The similarity between the 7S K RNA of the lamprey and that of mammals was 68%. Interestingly, several interspersed elements were found with nearly 100% similarity compared with mammals. In contrast to the lamprey, no 7S K RNA-related sequences were detectable among hagfish RNA, neither in northern blots nor with the PCR assay. In view of the significant conservation between the 7S K RNA of lamprey and that of mammals (human), this unexpected result clearly separates lamprey and hagfish. In addition, the lack of detectable 7S K RNA sequences in an outgroup, such as amphioxus, indicates that these results do not reflect an autapomorphy of hagfish. Therefore, our data provide additional support to the notion of a sister group relationship between Petromyzontiformes and gnathostomous vertebrates to the exclusion of Myxiniformes. Received: 24 September 1999 / Accepted: 9 February 2000