Defining 5S rRNA structure space: point mutation data can be used to predict the phenotype of multichange variants

A portion of the 5S ribosomal RNA (rRNA) structure space in the vicinity of the Vibrio proteolyticus 5S rRNA sequence is explored in detail with the intention of establishing principles that will allow a priori prediction of which sequences would be valid members of a particular RNA structure space. Four hundred and one sequence variants differing from the V. proteolyticus 5S rRNA wild-type sequence in 1-7 positions were characterized using an in vivo assay system. Most significantly, it was found that in general, the phenotypic effects of single changes were independent of the phenotypic effect of a second change. As a result, it was possible to use the new data in conjunction with results from prior studies of the same RNA to develop "truth tables" to predict which multiple change variants would be functional and which would be nonfunctional. The actual phenotype of 93.8% of the multichange variants studied was consistent with the predictions made using truth tables thereby providing for perhaps the first time an upper limit estimate of how frequent unexpected interactions are. It was also observed that single changes at positions involved in secondary structure were no more likely to be invalid than changes in other regions. In particular, internal changes in long standard stems were in fact almost always tolerated. Changes at positions that were hypervariable in the context of an alignment of related sequences were, as expected, usually found to be valid. However, the potential validity of changes that were idiosyncratic to a single lineage of related sequences when placed in the V. proteolyticus 5S rRNA context was unpredictable.     

Equally Parsimonious Pathways Through an RNA Sequence Space Are Not Equally Likely

An experimental system for determining the potential ability of sequences resembling 5S ribosomal RNA (rRNA) to perform as functional 5S rRNAs in vivo in the Escherichia coli cellular environment was devised previously. Presumably, the only 5S rRNA sequences that would have been fixed by ancestral populations are ones that were functionally valid, and hence the actual historical paths taken through RNA sequence space during 5S rRNA evolution would have most likely utilized valid sequences. Herein, we examine the potential validity of all sequence intermediates along alternative equally parsimonious trajectories through RNA sequence space which connect two pairs of sequences that had previously been shown to behave as valid 5S rRNAs in E. coli. The first trajectory requires a total of four changes. The 14 sequence intermediates provide 24 apparently equally parsimonious paths by which the transition could occur. The second trajectory involves three changes, six intermediate sequences, and six potentially equally parsimonious paths. In total, only eight of the 20 sequence intermediates were found to be clearly invalid. As a consequence of the position of these invalid intermediates in the sequence space, seven of the 30 possible paths consisted of exclusively valid sequences. In several cases, the apparent validity/invalidity of the intermediate sequences could not be anticipated on the basis of current knowledge of the 5S rRNA structure. This suggests that the interdependencies in RNA sequence space may be more complex than currently appreciated. If ancestral sequences predicted by parsimony are to be regarded as actual historical sequences, then the present results would suggest that they should also satisfy a validity requirement and that, in at least limited cases, this conjecture can be tested experimentally. Received: 27 August 1996 / Accepted: 14 April 1997     

Experimental investigation of an RNA sequence space

Modern rRNAs are the historic consequence of an ongoing evolutionary exploration of a sequence space. These extant sequences belong to a special subset of the sequence space that is comprised only of those primary sequences that can validly perform the biological function(s) required of the particular RNA. If it were possible to readily identify all such valid sequences, stochastic predictions could be made about the relative likelihood of various evolutionary pathways available to an RNA. Herein an experimental system which can assess whether a particular sequence is likely to have validity as a eubacterial 5S rRNA is described. A total of ten naturally occurring, and hence known to be valid, sequences and two point mutants of unknown validity were used to test the usefulness of the approach. Nine of the ten valid sequences tested positive whereas both mutants tested as clearly defective. The tenth valid sequence gave results that would be interpreted as reflecting a borderline status were the answer not known. These results demonstrate that it is possible to experimentally determine which sequences in local regions of the sequence space are potentially valid 5S rRNAs. This approach will allow direct study of the constraints governing RNA evolution and allow inquiry into how the last common ancestor of extant life apparently came to have very complex ribosomal RNAs that subsequently were very conserved.     

Human ribosomal RNA variants from a single individual and their expression in different tissues.

We have investigated the extent of sequence variation in human ribosomal RNA (rRNA) genes and the expression of specific rRNA gene variants in different tissues of an individual. Focusing on the fifth variable region (V5; nt 2065-2244) of the 28S rRNA gene, we find that sequence differences between rRNA genes of a single individual are characterized by differences in number of repeats of simple sequences at four specific sites. These data support and extend previous findings which show similar V5 sequence variation in rRNA genes from a group of individuals. We performed experiments to determine if there is differential gene expression within the rRNA multigene family. From the analysis of data of six variant V5 probes protected from RNase digestion by rRNAs isolated from different tissues of the individual, we conclude that each variant rRNA is present in a similar proportion in these tissues, whereas the actual contributions of variants differ, their relative proportion is maintained from tissue to tissue in an individual. We favor the explanation of a gene dosage effect over that of a regulated gene effect to account for this pattern of rRNA gene expression. In addition, computer generated secondary structure models of each V5 clone structure predict the same three helix structure with the regions of sequence variation contained in one stem-loop structure.     

Composite structure of the chloroplast 23 S ribosomal RNA genes of Chlamydomonas reinhardii: Evolutionary and functional implications

The chloroplast ribosomal unit of Chlamydomonas reinhardii displays two features which are not shared by other chloroplast ribosomal units. These include the presence of an intron in the 23 S ribosomal RNA gene and of two small genes coding for 3 S and 7 S rRNA in the spacer between the 16 S and 23 S rRNA genes (Rochaix & Malnoë, 1978). Sequencing of the 7 S and 3 S rRNAs as well as their genes and neighbouring regions has shown that: (1) the 7 S and 3 S rRNA genes are 282 and 47 base-pairs long, respectively, and are separated by a 23 base-pair A + T-rich spacer. (2) A sequence microheterogeneity exists within the 3 S RNA genes. (3) The sequences of the 7 S and 3 S rRNAs are homologous to the 5′ termini of prokaryotic and other chloroplast 23 S rRNAs, indicating that the C. reinhardii counterparts of 23 S rRNA have a composite structure. (4) The sequences of the 7 S and 3 S rRNAs are related to that of cytoplasmic 5.8 S rRNA, suggesting that these RNAs may perform similar functions in the ribosome. (5) Partial nucleotide sequence complementarity is observed between the 5′ ends of the 7 S and 3 S RNAs on one hand and the 23 S rRNA sequences which flank the ribosomal intron on the other. These data are compatible with the idea that these small rRNAs may play a role in the processing of the 23 S rRNA precursor.     

The 3''-terminal region of bacterial 23S ribosomal RNA: structure and homology with the 3''-terminal region of eukaryotic 28S rRNA and with chloroplast 4.5s rRNA.

The sequence of the 110 nucleotide fragment located at the 3'-end of E.coli, P.vulgaris and A.punctata 23S rRNAs has been determined. The homology between the E.coli and P.vulgaris fragments is 90%, whereas that between the E.coli and A.punctate fragments is only 60%. The three rRNA fragments have sequences compatible with a secondary structure consisting of two hairpins. Using chemical and enzymatic methods recently developed for the study of the secondary structure of RNA, we demonstrated that one of these hairpins and part of the other are actually present in the three 3'-terminal fragments in solution. This supports the existence of these two hairpins in the intact molecule. Indeed, results obtained upon limited digestion of intact 23S RNA with T1 RNase were in good agreement with the existence of these two hairpins. We observed that the primary structures of the 3'-terminal regions of yeast 26S rRNA and X.laevis 28S rRNA are both compatible with a secondary structure similar to that found at the 3'-end of bacterial 23S rRNAs. Furthermore, both tobacco and wheat chloroplast 4.5S rRNAs can also be folded in a similar way as the 3'-terminal region of bacterial 23S rRNA, the 3'-end of chloroplast 4.5S rRNAs being complementary to the 5'-end of chloroplast 23S rRNA. This strongly reinforces the hypothesis that chloroplast 4.5S rRNA originates from the 3'-end of bacterial 23S rRNA and suggests that this rRNA may be base-paired with the 5'-end of chloroplast 23S rRNA. Invariant oligonucleotides are present at identical positions in the homologous secondary structures of E.coli 23S, yeast 26S, X.laevis 28S and wheat and tobacco 4.5S rRNAs. Surprisingly, the sequences of these oligonucleotides are not all conserved in the 3'-terminal regions of A.punctata or even P.vulgaris 23S rRNAs. Results obtained upon mild methylation of E.coli 50S subunits with dimethylsulfate strongly suggest that these invariant oligonucleotides are involved in RNA tertiary structure or in RNA-protein interactions.     

Comparison of the nucleotide sequence of soybean 18S rRNA with the sequences of other small-subunit rRNAs

We present the sequence of the nuclear-encoded ribosomal small-subunit RNA from soybean. The soybean 18S rRNA sequence of 1807 nucleotides (nt) is contained in a gene family of approximately 800 closely related members per haploid genome. This sequence is compared with the ribosomal small-subunit RNAs of maize (1805 nt), yeast (1789 nt), Xenopus (1825 nt), rat (1869 nt), and Escherichia coli (1541 nt). Significant sequence homology is observed among the eukaryotic small-subunit rRNAs examined, and some sequence homology is observed between eukaryotic and prokaryotic small-subunit rRNAs. Conserved regions are found to be interspersed among highly diverged sequences. The significance of these comparisons is evaluated using computer simulation of a random sequence model. A tentative model of the secondary structure of soybean 18S rRNA is presented and discussed in the context of the functions of the various conserved regions within the sequence. On the basis of this model, the short base-paired sequences defining the four structural and functional domains of all 18S rRNAs are seen to be well conserved. The potential roles of other conserved soybean 18S rRNA sequences in protein synthesis are discussed.     

Nucleotide sequence of 5S ribosomal RNA from Aspergillus nidulans and Neurospora crassa

The nucleotide sequences of 5S rRNA molecules isolated from the cytosol and the mitochondria of the ascomycetes A. nidulans and N. crassa were determined by partial chemical cleavage of 3'-terminally labelled RNA. The sequence identity of the cytosolic and mitochondrial RNA preparations confirms the absence of mitochondrion-specific 5S rRNA in these fungi. The sequences of the two organisms differ in 35 positions, and each sequence differs from yeast 5S rRNA in 44 positions. Both molecules contain the sequence GCUC in place of GAAC or GAUY found in all other 5S rRNAs, indicating that this region is not universally involved in base-pairing to the invariant GTpsiC sequence of tRNAs.     

Rapid in vivo exploration of a 5S rRNA neutral network

A partial knockout compensation method to screen 5S ribosomal RNA sequence variants in vivo is described. The system utilizes an Escherichia coli strain in which five of eight genomic 5S rRNA genes were deleted in conjunction with a plasmid which is compensatory when carrying a functionally active 5S rRNA. The partial knockout strain is transformed with a population of potentially compensatory plasmids each carrying a randomly generated 5S rRNA gene variant. a The ability to compensate the slow growth rate of the knockout strain is used in conjunction with sequencing to rapidly identify variant 5S rRNAs that are functional as well as those that likely are not. The assay is validated by showing that the growth rate of 15 variants separately expressed in the partial knockout strain can be accurately correlated with in vivo assessments of the potential validity of the same variants. A region of 5S rRNA was mutagenized with this approach and nine novel variants were recovered and characterized. Unlike a complete knockout system, the method allows recovery of both deleterious and functional variants.. The method can be used to study variants of any 5S rRNA in the E. coli context including those of E. coli.     

Comparison of the nucleotide sequence of soybean 18S rRNA with the sequences of other small-subunit rRNAs

Summary We present the sequence of the nuclearencoded ribosomal small-subunit RNA from soybean. The soybean 18S rRNA sequence of 1807 nucleotides (nt) is contained in a gene family of approximately 800 closely related members per haploid genome. This sequence is compared with the ribosomal small-subunit RNAs of maize (1805 nt), yeast (1789 nt),Xenopus (1825 nt), rat (1869 nt), andEscherichia coli (1541 nt). Significant sequence homology is observed among the eukaryotic small-subunit rRNAs examined, and some sequence homology is observed between eukaryotic and prokaryotic small-subunit rRNAs. Conserved regions are found to be interspersed among highly diverged sequences. The significance of these comparisons is evaluated using computer simulation of a random sequence model. A tentative model of the secondary structure of soybean 18S rRNA is presented and discussed in the context of the functions of the various conserved regions within the sequence. On the basis of this model, the short basepaired sequences defining the four structural and functional domains of all 18S rRNAs are seen to be well conserved. The potential roles of other conserved soybean 18S rRNA sequences in protein synthesis are discussed.     

The secondary structure of the protein L1 binding region of ribosomal 23S RNA. Homologies with putative secondary structures of the L11 mRNA and of a region of mitochondrial 16S rRNA.

An heterologous complex was formed between E. coli protein L1 and P. vulgaris 23S RNA. We determined the primary structure of the RNA region which remained associated with protein L1 after RNase digestion of this complex. We also identified the loci of this RNA region which are highly susceptible to T1, S1 and Naja oxiana nuclease digestions respectively. By comparison of these results with those previously obtained with the homologous regions of E. coli and B. stearothermophilus 23S RNAs, we postulate a general structure for the protein L1 binding region of bacterial 23S RNA. Both mouse and human mit 16S rRNAs and Xenopus laevis and Tetrahymena 28S rRNAs contain a sequence similar to the E. coli 23s RNS region preceding the L1 binding site. The region of mit 16S rRNA which follows this sequence has a potential secondary structure bearing common features with the L1-associated region of bacterial 23S rRNA. The 5'-end region of the L11 mRNA also has several sequence potential secondary structures displaying striking homologies with the protein L1 binding region of 23S rRNA and this probably explains how protein L1 functions as a translational repressor. One of the L11 mRNA putative structures bears the features common to both the L1-associated region of bacterial 23S rRNA and the corresponding region of mit 16S rRNA.     

Characterization of a separate small domain derived from the 5' end of 23S rRNA of an alpha-proteobacterium.

We demonstrate the presence of a separate processed domain derived from the 5' end of 23S rRNA in ribosomes of Rhodopseudomonas palustris, a member of the alpha-++proteobacteria. Previous sequencing studies predicted intervening sequences (IVS) at homologous positions within the 23S rRNA genes of several alpha-proteobacteria, including R.palustris, and we find a processed 23S rRNA 5' domain in unfractionated RNA from several species. 5.8S rRNA from eukaryotic cytoplasmic large subunit ribosomes and the bacterial processed 23S rRNA 5' domain share homology, possess similar structures and are both derived by processing of large precursors. However, the internal transcribed spacer regions or IVSs separating them from the main large subunit rRNAs are evolutionarily unrelated. Consistent with the difference in sequence, we find that the site and mechanism of IVS processing also differs. Rhodopseudomonas palustris IVS-containing RNA precursors are cleaved in vitro by Escherichia coli RNase III or a similar activity present in R.palustris extracts at a processing site distinct from that found in eukaryotic systems and this results in only partial processing of the IVS. Surprisingly, in a reaction unlike characterized cases of eubacterial IVS processing, an RNA segment larger than the corresponding DNA insertion is removed which contains conserved sequences. These sequences, by analogy, serve to link the 23S rRNA 5' rRNA domains or 5.8S rRNAs to the main portion of other prokaryotic 23S rRNAs or to eukaryotic 28S rRNAs, respectively.     

A misfolded form of 5S rRNA is complexed with the Ro and La autoantigens.

In both vertebrate and invertebrate cells, the 60-kDa Ro autoantigen is bound to small cytoplasmic RNAs known as Y RNAs. In Xenopus oocytes, the 60-kDa Ro protein is also complexed with a class of 5S rRNA precursors that contain internal mutations. Because these 5S rRNA precursors are processed inefficiently and degraded eventually, the Ro protein may function in a quality control pathway for 5S rRNA biosynthesis. We have investigated the sequence and secondary structure determinants in the mutant 5S rRNAs that confer binding by the 60-kDa Ro protein. The mutant 5S rRNAs fold to form an alternative helix that is required for recognition by the 60-kDa Ro protein. Mutations that disrupt the alternative helix eliminate Ro protein binding, whereas compensatory changes that restore the helix are bound efficiently by the Ro protein. When the structure of the mutant RNA was probed using dimethylsulfate and oligonucleotide-directed RNase H cleavage, the results were consistent with the formation of the alternative structure. The La protein, which is also complexed with the mutant 5S rRNA precursors, protects similar sequences from nuclease digestion as does the 60-kDa Ro protein. Thus, the binding sites for these two proteins are either nearby on the RNA, or the two proteins may be complexed through protein-protein interactions. When the human Ro protein is expressed in the yeast Saccharomyces cerevisiae, the protein binds wild-type 5S rRNA precursors, suggesting that a population of wild-type precursors also folds into the alternative structure.