Common 5S rRNA Variants Are Likely To Be Accepted in Many Sequence Contexts

Over evolutionary time RNA sequences which are successfully fixed in a population are selected from among those that satisfy the structural and chemical requirements imposed by the function of the RNA. These sequences together comprise the structure space of the RNA. In principle, a comprehensive understanding of RNA structure and function would make it possible to enumerate which specific RNA sequences belong to a particular structure space and which do not. We are using bacterial 5S rRNA as a model system to attempt to identify principles that can be used to predict which sequences do or do not belong to the 5S rRNA structure space. One promising idea is the very intuitive notion that frequently seen sequence changes in an aligned data set of naturally occurring 5S rRNAs would be widely accepted in many other 5S rRNA sequence contexts. To test this hypothesis, we first developed well-defined operational definitions for a Vibrio region of the 5S rRNA structure space and what is meant by a highly variable position. Fourteen sequence variants (10 point changes and 4 base-pair changes) were identified in this way, which, by the hypothesis, would be expected to incorporate successfully in any of the known sequences in the Vibrio region. All 14 of these changes were constructed and separately introduced into the Vibrio proteolyticus 5S rRNA sequence where they are not normally found. Each variant was evaluated for its ability to function as a valid 5S rRNA in an E. coli cellular context. It was found that 93% (13/14) of the variants tested are likely valid 5S rRNAs in this context. In addition, seven variants were constructed that, although present in the Vibrio region, did not meet the stringent criteria for a highly variable position. In this case, 86% (6/7) are likely valid. As a control we also examined seven variants that are seldom or never seen in the Vibrio region of 5S rRNA sequence space. In this case only two of seven were found to be potentially valid. The results demonstrate that changes that occur multiple times in a local region of RNA sequence space in fact usually will be accepted in any sequence context in that same local region.     

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     

A novel approach for monitoring genetically engineered microorganisms by using artificial, stable RNAs.

Further improvements in technology for efficient monitoring of genetically engineered microorganisms (GEMs) in the environment are needed. Technology for monitoring rRNA is well established but has not generally been applicable to GEMs because of the lack of unique rRNA target sequences. In the work described herein, it is demonstrated that a deletion mutant of a plasmid-borne Vibrio proteolyticus 5S rRNA gene continues to accumulate to high levels in Escherichia coli although it is no longer incorporated into 70S ribosomes. This deletion construct was subsequently modified by mutagenesis to create a unique recognition site for the restriction endonuclease BstEII, into which new sequences could be readily inserted. Finally, a novel 17-nucleotide identifier sequence from Pennisetum purpureum was embedded into the construct to create an RNA identification cassette. The artificial identifier RNA, expressed from this cassette in vivo, accumulated in E. coli to levels comparable to those of wild-type 5S rRNA without being seriously detrimental to cell survival in laboratory experiments and without entering the ribosomes. These results demonstrate that artificial, stable RNAs containing sequence segments remarkably different from those present in any known rRNA can be designed and that neither the deleted sequence segment nor ribosome incorporation is essential for accumulation of an RNA product.     

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.     

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.     

Probing the conformation of 26S rRNA in yeast 60S ribosomal subunits with kethoxal

The conformation and accessibility of 26S rRNA in yeast 60S ribosomal subunits were probed with kethoxal. Oligonucleotides originating from reactive sites were isolated by diagonal electrophoresis and sequenced. From over 70 oligonucleotide sequences, 26 kethoxal-reactive sites could be placed in the 26S rRNA sequence. These are in close agreement with a proposed secondary structure model for the RNA that is based on comparative sequence analysis. At least seven kethoxal-reactive sites in yeast 26S rRNA are in positions that are exactly homologous to reactive positions in E. coli 23S rRNA; each of these sites has previously been implicated in some aspect of ribosomal function.     

A computer method for finding common base paired helices in aligned sequences: application to the analysis of random sequences.

We describe a new computer program that identifies conserved secondary structures in aligned nucleotide sequences of related single-stranded RNAs. The program employs a series of hash tables to identify and sort common base paired helices that are located in identical positions in more than one sequence. The program gives information on the total number of base paired helices that are conserved between related sequences and provides detailed information about common helices that have a minimum of one or more compensating base changes. The program is useful in the analysis of large biological sequences. We have used it to examine the number and type of complementary segments (potential base paired helices) that can be found in common among related random sequences similar in base composition to 16S rRNA from Escherichia coli. Two types of random sequences were analyzed. One set consisted of sequences that were independent but they had the same mononucleotide composition as the 16S rRNA. The second set contained sequences that were 80% similar to one another. Different results were obtained in the analysis of these two types of random sequences. When 5 sequences that were 80% similar to one another were analyzed, significant numbers of potential helices with two or more independent base changes were observed. When 5 independent sequences were analyzed, no potential helices were found in common. The results of the analyses with random sequences were compared with the number and type of helices found in the phylogenetic model of the secondary structure of 16S ribosomal RNA. Many more helices are conserved among the ribosomal sequences than are found in common among similar random sequences. In addition, conserved helices in the 16S rRNAs are, on the average, longer than the complementary segments that are found in comparable random sequences. The significance of these results and their application in the analysis of long non-ribosomal nucleotide sequences is discussed.     

7 S RNA: A single site substrate for the RNA processing enzyme ribonuclease E of Escherichia coli

7 S RNA accumulates at non-permissive temperatures in an RNAase E strain containing the recombinant plasmid pJR3Δ which carries a single 5 S rRNA gene and expression sequences. 7 S RNA is a processing intermediate that contains the complete sequence of 5 S rRNA as well as a stem-and-loop structure encoded by the terminator of rrnD. 7 S RNA can be processed in vitro by RNAase E. Structural analysis of the products (5 S rRNA and the stem) of in vitro processing of 7 S RNA revealed that the cleavage site of RNAase E in 7 S RNA is 3 nucleotides downstream from the 3′ end of the mature 5 S rRNA. The cleavage generates 3′-hydroxyl and 5′-phosphate termini.     

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.     

Isolation and characterization of the 5S rRNA gene of Leptospira interrogans.

The gene encoding the 5S rRNA for Leptospira interrogans serovar canicola strain Moulton was isolated and sequenced. The 5S rRNA gene occurs as a single copy within the genome and encodes a 117-nucleotide-long RNA molecule. The 5S rRNA gene is flanked at both the 5' and 3' ends by regions of A + T-rich sequences, and the 5'-flanking region contains a promoter sequence. L. interrogans has a unique and remarkable organization of the 5S rRNA gene. The 5S rRNA molecule exhibits a strong similarity to typical eubacterial 5S rRNA in terms of overall secondary structure, while the primary sequence is conserved to a lesser degree. Restriction analysis of the 5S rRNA gene indicated that the DNA sequence including the 5S rRNA gene is highly conserved in the genomes of parasitic leptospires.     

Three new small nucleolar RNAs that are psoralen cross-linked in vivo to unique regions of pre-rRNA.

We have recently described three novel human small nucleolar RNA species with unique nucleotide sequences, which were named E1, E2, and E3. The present article describes specific psoralen photocross-linking in whole HeLa cells of E1, E2, and E3 RNAs to nucleolar pre-rRNA. These small RNAs were cross-linked to different sections of pre-rRNA. E1 RNA was cross-linked to two segments of nucleolar pre-rRNA; one was within residues 697 to 1163 of the 5' external transcribed spacer, and the other one was between nucleotides 664 and 1021 of the 18S rRNA sequence. E2 RNA was cross-linked to a region within residues 3282 to 3667 of the 28S rRNA sequence. E3 RNA was cross-linked to a sequence between positions 1021 and 1639 of the 18S rRNA sequence. Primer extension analysis located psoralen adducts in E1, E2, and E3 RNAs that were enriched in high-molecular-weight fractions of nucleolar RNA. Some of these psoralen adducts might be cross-links of E1, E2, and E3 RNAs to large nucleolar RNA. Antisense oligodeoxynucleotide-targeted RNase H digestion of nucleolar extracts revealed accessible segments in these three small RNAs. The accessible regions were within nucleotide positions 106 to 130 of E1 RNA, positions 24 to 48 and 42 to 66 of E2 RNA, and positions 7 to 16 and about 116 to 122 of E3 RNA. Some of the molecules of these small nucleolar RNAs sedimented as if associated with larger structures when both nondenatured RNA and a nucleolar extract were analyzed.     

Identifying constraints on the higher-order structure of RNA: continued development and application of comparative sequence analysis methods.

Comparative sequence analysis addresses the problem of RNA folding and RNA structural diversity, and is responsible for determining the folding of many RNA molecules, including 5S, 16S, and 23S rRNAs, tRNA, RNAse P RNA, and Group I and II introns. Initially this method was utilized to fold these sequences into their secondary structures. More recently, this method has revealed numerous tertiary correlations, elucidating novel RNA structural motifs, several of which have been experimentally tested and verified, substantiating the general application of this approach. As successful as the comparative methods have been in elucidating higher-order structure, it is clear that additional structure constraints remain to be found. Deciphering such constraints requires more sensitive and rigorous protocols, in addition to RNA sequence datasets that contain additional phylogenetic diversity and an overall increase in the number of sequences. Various RNA databases, including the tRNA and rRNA sequence datasets, continue to grow in number as well as diversity. Described herein is the development of more rigorous comparative analysis protocols. Our initial development and applications on different RNA datasets have been very encouraging. Such analyses on tRNA, 16S and 23S rRNA are substantiating previously proposed associations and are now beginning to reveal additional constraints on these molecules. A subset of these involve several positions that correlate simultaneously with one another, implying units larger than a basepair can be under a phylogenetic constraint.     

A proposed nucleotide sequence for the 5S ribosomal ribonucleic acid of rainbow trout (Salmo gairdneri).

Rainbow trout cell cultures have been exposed to 32P-labelled inorganic phosphate and the labelled RNA has been isolated. The 5S ribosomal ribonucleic acid (5S rRNA) was purified by polyacrylamide gel electrophoresis, then digested with RNase T1 or pancreatic RNase. The products of complete digestion were separated and their sequences determined. These analyses have allowed a sequence to be proposed which differs in eight positions from that of mammalian 5S rRNAs.     

裸子植物5S rRNA基因序列变异及二级结构特征

在高等植物中,5SrRNA基因一级结构是高度保守的,二级结构也相当一致。通过比较18种裸子植物5SrRNA基因序列和二级结构变异,发现55％的核苷酸位点是可变的,这种变异有68％发生在干区（双链区）,其中一些变异,如双链的互补性核苷酸替代,GU配对等能够维系5SrRNA二级结构的稳定性。环区相对保守,这与5SrRNA三级结构折叠或在转录翻译过程中蛋白质、RNA的结合相关。另外,首次报道了松属环E区核苷酸的变异性,这可能与其他区域的变异一样,是假基因造成的结果。5SrRNA基因信息可反映大分类群的系统进化关系,但由于基因长度短,信息量小,其在近缘种系统分类的应用受到限制。     

Effect of point mutations on 5.8S ribosomal ribonucleic acid secondary structure and the 5.8S--28S ribosomal ribonucleic acid junction

Naturally occurring differences in the nucleotide sequences of 5.8S ribosomal ribonucleic acids (rRNAs) from a variety of organisms have been used to study the role of specific nucleotides in the secondary structure and intermolecular interactions of this RNA. Significant differences in the electrophoretic mobilities of free 5.8S RNAs and the thermal stabilities of 5.8S--28S rRNA complexes were observed even in such closely related sequences as those of man, rat, turtle, and chicken. A single base transition from a guanylic acid residue in position 2 in mammalian 5.8S rRNA to an adenylic acid residue in turtle and chicken 5.8S rRNA results both in a more open molecular conformation and in a 5.8S--28S rRNA junction which is 3.5 degrees C more stable to thermal denaturation. Other changes such as the deletion of single nucleotides from either the 5' or the 3' terminals have no detectable effect on these features. The results support secondary structure models for free 5.8S rRNA in which the termini interact to various degrees and 5.8S--28S rRNA junctions in which both termini of the 5.8S molecule interact with the cognate high molecular weight RNA component.