Common structures of the 5' non-coding RNA in enteroviruses and rhinoviruses. Thermodynamical stability and statistical significance

A total of 4051 suboptimal secondary structures are predicted by folding the 5' non-coding region of ten polioviruses, five human rhinoviruses and three coxsackieviruses using our new suboptimal folding algorithm for the prediction of both optimal and suboptimal RNA secondary structures. A comparative analysis of these RNA secondary structures reveals the conservation of common secondary structure that can be supported by phylogenetic data. The thermodynamic stability and statistical significance of these predicted, conserved helical elements are assessed and significant structure motifs in the 5' non-coding region are proposed. The possible roles of these structure motifs in the virus life cycle are discussed.     

Phylogenetic analysis and evolution of RNase P RNA in proteobacteria.

The secondary structures of the eubacterial RNase P RNAs are being elucidated by a phylogenetic comparative approach. Sequences of genes encoding RNase P RNA from each of the recognized subgroups (alpha, beta, gamma, and delta) of the proteobacteria have now been determined. These sequences allow the refinement, to nearly the base pair level, of the phylogenetic model for RNase P RNA secondary structure. Evolutionary change among the RNase P RNAs was found to occur primarily in four discrete structural domains that are peripheral to a highly conserved core structure. The new sequences were used to examine critically the proposed similarity (C. Guerrier-Takada, N. Lumelsky, and S. Altman, Science 246:1578-1584, 1989) between a portion of RNase P RNA and the "exit site" of the 23S rRNA of Escherichia coli. Phylogenetic comparisons indicate that these sequences are not homologous and that any similarity in the structures is, at best, tenuous.     

Conserved features of Y RNAs revealed by automated phylogenetic secondary structure analysis.

Y RNAs are small 'cytoplasmic' RNAs which are components of the Ro ribonucleoprotein (RNP) complex. The core of this complex, which is found in the cell nuclei of higher eukaryotes as well as the cytoplasm, is composed of a complex between the 60 kDa Ro protein and Y RNAs. Human cells contain four distinct Y RNAs (Y1, Y3, Y4 and Y5), while other eukaryotes contain a variable number of Y RNA homologues. When detected in a particular species, the Ro RNP has been present in every cell type within that particular organism. This characteristic, along with its high conservation among vertebrates, suggests an important function for Ro RNP in cellular metabolism; however, this function has not yet been definitively elucidated. In order to identify conserved features of Y RNA sequences and structures which may be directly involved in Ro RNP function, a phylogenetic comparative analysis of Y RNAs has been performed. Sequences of Y RNA homologues from five vertebrate species have been obtained and, together with previously published Y RNA sequences, used to predict Y RNA secondary structures. A novel RNA secondary structure comparison algorithm, the suboptimal RNA analysis program, has been developed and used in conjunction with available algorithms to find phylogenetically conserved secondary structure models for YI, Y3 and Y4 RNAs. Short, conserved sequences within the Y RNAs have been identified and are invariant among vertebrates, consistent with a direct role for Y RNAs in Ro function. A subset of these are located wholly or partially in looped regions in the Y3 and Y4 RNA predicted model structures, in accord with the possibility that these Y RNAs base pair with other cellular nucleic acids or are sites of interaction between the Ro RNP and other macromolecules.     

Prediction of common folding structures of homologous RNAs.

We have developed an algorithm and a computer program for simultaneously folding homologous RNA sequences. Given an alignment of M homologous sequences of length N, the program performs phylogenetic comparative analysis and predicts a common secondary structure conserved in the sequences. When the structure is not uniquely determined, it infers multiple structures which appear most plausible. This method is superior to energy minimization methods in the sense that it is not sensitive to point mutation of a sequence. It is also superior to usual phylogenetic comparative methods in that it does not require manual scrutiny for covariation or secondary structures. The most plausible 1-5 structures are produced in O(MN2 + N3) time and O(N2) space, which are the same requirements as those of widely used dynamic programs based on energy minimization for folding a single sequence. This is the first algorithm probably practical both in terms of time and space for finding secondary structures of homologous RNA sequences. The algorithm has been implemented in C on a Sun SparcStation, and has been verified by testing on tRNAs, 5S rRNAs, 16S rRNAs, TAR RNAs of human immunodeficiency virus type 1 (HIV-1), and RRE RNAs of HIV-1. We have also applied the program to cis-acting packaging sequences of HIV-1, for which no generally accepted structures yet exist, and propose potentially stable structures. Simulation of the program with random sequences with the same base composition and the same degree of similarity as the above sequences shows that structures common to homologous sequences are very unlikely to occur by chance in random sequences.     

The suboptimal structures find the optimal RNAs: homology search for bacterial non-coding RNAs using suboptimal RNA structures

Non-coding RNAs (ncRNAs) are regulatory molecules encoded in the intergenic or intragenic regions of the genome. In prokaryotes, biocomputational identification of homologs of known ncRNAs in other species often fails due to weakly evolutionarily conserved sequences, structures, synteny and genome localization, except in the case of evolutionarily closely related species. To eliminate results from weak conservation, we focused on RNA structure, which is the most conserved ncRNA property. Analysis of the structure of one of the few well-studied bacterial ncRNAs, 6S RNA, demonstrated that unlike optimal and consensus structures, suboptimal structures are capable of capturing RNA homology even in divergent bacterial species. A computational procedure for the identification of homologous ncRNAs using suboptimal structures was created. The suggested procedure was applied to strongly divergent bacterial species and was capable of identifying homologous ncRNAs.     

A conserved secondary structure for telomerase RNA.

The RNA moiety of the ribonucleoprotein enzyme telomerase contains the template for telomeric DNA synthesis. We present a secondary structure model for telomerase RNA, derived by a phylogenetic comparative analysis of telomerase RNAs from seven tetrahymenine ciliates. The telomerase RNA genes from Tetrahymena malaccensis, T. pyriformis, T. hyperangularis, T. pigmentosa, T. hegewishii, and Glaucoma chattoni were cloned, sequenced, and compared with the previously cloned RNA gene from T. thermophila and with each other. To define secondary structures of these RNAs, homologous complementary sequences were identified by the occurrence of covariation among putative base pairs. Although their primary sequences have diverged rapidly overall, a strikingly conserved secondary structure was identified for all these telomerase RNAs. Short regions of nucleotide conservation include a block of 22 totally conserved nucleotides that contains the telomeric templating region.     

Identification and analysis of ribonuclease P and MRP RNA in a broad range of eukaryotes

RNases P and MRP are ribonucleoprotein complexes involved in tRNA and rRNA processing, respectively. The RNA subunits of these two enzymes are structurally related to each other and play an essential role in the enzymatic reaction. Both of the RNAs have a highly conserved helical region, P4, which is important in the catalytic reaction. We have used a bioinformatics approach based on conserved elements to computationally analyze available genomic sequences of eukaryotic organisms and have identified a large number of novel nuclear RNase P and MRP RNA genes. For MRP RNA for instance, this investigation increases the number of known sequences by a factor of three. We present secondary structure models of many of the predicted RNAs. Although all sequences are able to fold into the consensus secondary structure of P and MRP RNAs, a striking variation in size is observed, ranging from a Nosema locustae MRP RNA of 160 nt to much larger RNAs, e.g. a Plasmodium knowlesi P RNA of 696 nt. The P and MRP RNA genes appear in tandem in some protists, further emphasizing the close evolutionary relationship of these RNAs.     

Modelling the secondary structures of slippage-prone hypervariable RNA regions: the example of the tiger beetle 18S rRNA variable region V4.

Variable regions within ribosomal RNAs frequently vary in length as a result of incorporating products of slippage. This makes constructing secondary structure models problematic because base homology is difficult or impossible to establish between species. Here, we model such a region by comparing the results of the MFOLD suboptimal folding algorithm for different species to identify conserved structures. Based on the reconstruction of base change on a phylogenetic tree of the species and comparison against null models of character change, we devise a statistical analysis to assess support of these structures from compensatory and semi-compensatory (i.e. G.C to G.U or A.U to G.U) mutations. As a model system we have used variable region V4 from cicindelid (tiger beetle) small subunit ribosomal RNAs (SSU rRNAs). This consists of a mixture of conserved and highly variable subregions and has been subject to extensive comparative analysis in the past. The model that results is similar to a previously described model of this variable region derived from a different set of species and contains a novel structure in the central, highly variable part. The method we describe may be useful in modelling other RNA regions that are subject to slippage.     

The 5S ribosomal RNA of Euglena gracilis cytoplasmic ribosomes is closely homologous to the 5S RNA of the trypanosomatid protozoa.

The complete nucleotide sequence of the major species of cytoplasmic 5S ribosomal RNA of Euglena gracilis has been determined. The sequence is: 5' GGCGUACGGCCAUACUACCGGGAAUACACCUGAACCCGUUCGAUUUCAGAAGUUAAGCCUGGUCAGGCCCAGUUAGUAC UGAGGUGGGCGACCACUUGGGAACACUGGGUGCUGUACGCUUOH3'. This sequence can be fitted to the secondary structural models recently proposed for eukaryotic 5S ribosomal RNAs (1,2). Several properties of the Euglena 5S RNA reveal a close phylogenetic relationship between this organism and the protozoa. Large stretches of nucleotide sequences in predominantly single-stranded regions of the RNA are homologous to that of the trypanosomatid protozoan Crithidia fasticulata. There is less homology when compared to the RNAs of the green alga Chlorella or to the RNAs of the higher plants. The sequence AGAAC near position 40 that is common to plant 5S RNAs is CGAUU in both Euglena and Crithidia. The Euglena 5S RNA has secondary structural features at positions 79-99 similar to that of the protozoa and different from that of the plants. The conclusions drawn from comparative studies of cytochrome c structures which indicate a close phylogenetic relatedness between Euglena and the trypanosomatid protozoa are supported by the comparative data with 5S ribosomal RNAs.     

Pattern analysis of RNA secondary structure similarity and consensus of minimal-energy folding

We describe an automated procedure to search for consensus structures or substructures in a set of homologous or related RNA molecules. The procedure is based on the calculation of optimal and sub-optimal secondary structures using thermodynamic rules for base-pairing by energy-minimization. A linear representation of the secondary structures of the related RNAs is used so that they can be compared and classified using standard alignment and clusterings programs. We illustrate the method by means of two sets of homologous small RNAs, U2 and U3, and a set of alpha-globin mRNAs and show that biologically interesting consensus structures are obtained.     

Phylogenetic comparative chemical footprint analysis of the interaction between ribonuclease P RNA and tRNA.

Ribonuclease P RNA is the catalytic moiety of the ribonucleoprotein enzyme that endonucleolytically cleaves precursor sequences from the 5' ends of pre-tRNAs. The bacterial RNase P RNA-tRNA complex was examined with a footprinting approach, utilizing chemical modification to determine RNase P RNA nucleotides that potentially contact tRNA. RNase P RNA was modified with dimethylsulfate or kethoxal in the presence or absence of tRNA, and sites of modification were detected by primer extension. Comparison of the results reveals RNase P bases that are protected from modification upon binding tRNA. Analyses were carried out with RNase P RNAs from three different bacteria: Escherichia coli, Chromatium vinosum and Bacillus subtilis. Discrete bases of these RNAs that lie within conserved, homologous portions of the secondary structures are similarly protected. One protection among all three RNAs was attributed to the precursor segment of pre-tRNA. Experiments using pre-tRNAs containing precursor segments of variable length demonstrate that a precursor segment of only 2-4 nucleotides is sufficient to confer this protection. Deletion of the 3'-terminal CCA sequence of tRNA correlates with loss of protection of a particular loop in the RNase P RNA secondary structure. Analysis of mutant tRNAs containing sequential 3'-terminal deletions suggests a relative orientation of the bound tRNA CCA to that loop.     

Efficient system of homologous RNA recombination in brome mosaic virus: sequence and structure requirements and accuracy of crossovers.

Brome mosaic virus (BMV), a tripartite positive-stranded RNA virus of plants engineered to support intersegment RNA recombination, was used for the determination of sequence and structural requirements of homologous crossovers. A 60-nucleotide (nt) sequence, common between wild-type RNA2 and mutant RNA3, supported efficient repair (90%) of a modified 3' noncoding region in the RNA3 segment by homologous recombination with wild-type RNA2 3' noncoding sequences. Deletions within this sequence in RNA3 demonstrated that a nucleotide identity as short as 15 nt can support efficient homologous recombination events, while shorter (5-nt) sequence identity resulted in reduced recombination frequency (5%) within this region. Three or more mismatches within a downstream portion of the common 60-nt RNA3 sequence affected both the incidence of recombination and the distribution of crossover sites, suggesting that besides the length, the extent of sequence identity between two recombining BMV RNAs is an important factor in homologous recombination. Site-directed mutagenesis of the common sequence in RNA3 did not reveal a clear correlation between the stability of predicted secondary structures and recombination activity. This indicates that homologous recombination does not require similar secondary structures between two recombining RNAs at the sites of crossovers. Nearly 20% of homologous recombinants were imprecise (aberrant), containing either nucleotide mismatches, small deletions, or small insertions within the region of crossovers. This implies that homologous RNA recombination is not as accurate as proposed previously. Our results provide experimental evidence that the requirements and thus the mechanism of homologous recombination in BMV differ from those of previously described heteroduplex-mediated nonhomologous recombination (P. D. Nagy and J. J. Bujarski, Proc. Natl. Acad. Sci. USA 90:6390-6394, 1993).     

Partial characterization of an RNA component that copurifies with Saccharomyces cerevisiae RNase P.

Saccharomyces cerevisiae cellular RNase P is composed of both protein and RNA components that are essential for activity. The isolated holoenzyme contains a highly structured RNA of 369 nucleotides that has extensive sequence similarities to the 286-nucleotide RNA associated with Schizosaccharomyces pombe RNase P but bears little resemblance to the analogous RNA sequences in procaryotes or S. cerevisiae mitochondria. Even so, the predicted secondary structure of S. cerevisiae RNA is strikingly similar to the bacterial phylogenetic consensus rather than to previously predicted structures of other eucaryotic RNase P RNAs.     

An algorithm for searching RNA motifs in genomic sequences

RNA molecules, which are found in all living cells, fold into characteristic structures that account for their diverse functional activities. Many of these RNA structures consist of a collection of fundamental RNA motifs. The various combinations of RNA basic components form different RNA classes and define their unique structural and functional properties. The availability of many genome sequences makes it possible to search computationally for functional RNAs. Biological experiments indicate that functional RNAs have characteristic RNA structural motifs represented by specific combinations of base pairings and conserved nucleotides in the loop regions. The searching for those well-ordered RNA structures and their homologues in genomic sequences is very helpful for the understanding of RNA-based gene regulation. In this paper, we consider the following problem: given an RNA sequence with a known secondary structure, efficiently determine candidate segments in genomic sequences that can potentially form RNA secondary structures similar to the given RNA secondary structure. Our new bottom-up approach searches all potential stem-loops similar to ones of the given RNA secondary structure first, and then based on located stem-loops, detects potential homologous structural RNAs in genomic sequences.     

Mutant prohead RNAs in the in vitro packaging of bacteriophage phi 29 DNA-gp3.

The 174-base prohead RNA encoded by bacteriophage phi 29 of Bacillus subtilis, essential for packaging of the DNA-gp3 (DNA-gene product 3) complex, was expressed efficiently from the cloned gene. Computer programs for RNA structure analysis were used to fold hypothetical RNA mutants and thus to target mutagenesis of the RNA for studies of structure and function. Five mutants of the RNA were then produced by oligonucleotide-directed mutagenesis that were altered in the primary sequence at selected sites; two of these mutants were predicted to be altered in secondary structure from a model established previously by a phylogenetic analysis. The binding of the 32P end-labeled mutant RNAs to RNA-free proheads was comparable with that of the wild-type RNA. However, the capability of the mutant RNAs to reconstitute RNA-free proheads for DNA-gp3 packaging in the defined in vitro system and for assembly of phage in RNA-free extracts was variable, depending upon the alteration. Changes of highly conserved bases that retained the predicted secondary structure of the RNA model were tolerated to a much greater extent than changes predicted to alter the RNA secondary structure.     

The nucleotide sequence of the chloroplast 5S ribosomal RNA from spinach.

Spinacia oleracia cholorplast 5S ribosomal RNA was end-labeled with [32P] and the complete nucleotide sequence was determined. The sequence is: pUAUUCUGGUGUCCUAGGCGUAGAGGAACCACACCAAUCCAUCCCGAACUUGGUGGUUAAACUCUACUGCGGUGACGAU ACUGUAGGGGAGGUCCUGCGGAAAAAUAGCUCGACGCCAGGAUGOH. This sequence can be fitted to the secondary structural model proposed for prokaryotic 5S ribosomal RNAs by Fox and Woese (1). However, the lengths of several single- and double-stranded regions differ from those common to prokaryotes. The spinach chloroplast 5S ribosomal RNA is homologous to the 5S ribosomal RNA of Lemna chloroplasts with the exception that the spinach RNA is longer by one nucleotide at the 3' end and has a purine base substitution at position 119. The sequence of spinach chloroplast 5S RNA is identical to the chloroplast 5S ribosomal RNA gene of tobacco. Thus the structures of the chloroplast 5S ribosomal RNAs from some of the higher plants appear to be almost totally conserved. This does not appear to be the case for the higher plant cytoplasmic 5S ribosomal RNAs.     

Phylogenetic comparative analysis and the secondary structure of ribonuclease P RNA--a review

The most incisive a priori approach to inferring the higher order structure of large RNAs has proven to be the use of phylogenetic comparisons. This article provides guidelines to the method, using as an illustration the elucidation of the secondary structure of the catalytic RNA subunit of ribonuclease P (RNase P). The resultant structure is compared to the possibilities that are predicted thermodynamically for the RNase P RNA sequences of nine eubacteria.     

The secondary structure of ribonuclease P RNA, the catalytic element of a ribonucleoprotein enzyme

Secondary structure models for the ribonuclease (RNAase) P RNAs of Bacillus subtilis and E. coli were derived by a phylogenetic comparative analysis of published sequences as well as four novel ones. The RNAase P RNA genes from Bacillus megaterium, Bacillus brevis, Bacillus stearothermophilus, and Pseudomonas fluorescens were cloned, sequenced, and compared with the other available sequences. Regions of pairing were identified by the occurrence of homologous complementary sequences that vary among the compared molecules. A common core of primary and secondary structure can be identified in all these RNAase P RNAs. The previously noted striking differences between the Bacillus and the enteric RNAase P RNAs arise not only from point mutations, but from the addition or deletion of structural domains. The primary and secondary structural features that are common to all of the RNAase P RNAs are likely to be the elements involved in the binding and cleavage of tRNA precursors, and in the interaction with the RNAase P protein.     

Phylogenetic analysis of the structure of RNase MRP RNA in yeasts

RNase MRP is a ribonucleoprotein enzyme involved in processing precursor rRNA in eukaryotes. To facilitate our structure-function analysis of RNase MRP from Saccharomyces cerevisiae, we have determined the likely secondary structure of the RNA component by a phylogenetic approach in which we sequenced all or part of the RNase MRP RNAs from 17 additional species of the Saccharomycetaceae family. The structure deduced from these sequences contains the helices previously suggested to be common to the RNA subunit of RNase MRP and the related RNA subunit of RNase P, an enzyme cleaving tRNA precursors. However, outside this common region, the structure of RNase MRP RNA determined here differs from a previously proposed universal structure for RNase MRPs. Chemical and enzymatic structure probing analyses were consistent with our revised secondary structure. Comparison of all known RNase MRP RNA sequences revealed three regions with highly conserved nucleotides. Two of these regions are part of a helix implicated in RNA catalysis in RNase P, suggesting that RNase MRP may cleave rRNA using a similar catalytic mechanism.     

Evolutionary Evidence for Alternative Structure in RNA Sequence Co-variation

Sequence conservation and co-variation of base pairs are hallmarks of structured RNAs. For certain RNAs (e.g. riboswitches), a single sequence must adopt at least two alternative secondary structures to effectively regulate the message. If alternative secondary structures are important to the function of an RNA, we expect to observe evolutionary co-variation supporting multiple conformations. We set out to characterize the evolutionary co-variation supporting alternative conformations in riboswitches to determine the extent to which alternative secondary structures are conserved. We found strong co-variation support for the terminator, P1, and anti-terminator stems in the purine riboswitch by extending alignments to include terminator sequences. When we performed Boltzmann suboptimal sampling on purine riboswitch sequences with terminators we found that these sequences appear to have evolved to favor specific alternative conformations. We extended our analysis of co-variation to classic alignments of group I/II introns, tRNA, and other classes of riboswitches. In a majority of these RNAs, we found evolutionary evidence for alternative conformations that are compatible with the Boltzmann suboptimal ensemble. Our analyses suggest that alternative conformations are selected for and thus likely play functional roles in even the most structured of RNAs.