- Research highlight
- Open Access
The amazing world of bacterial structured RNAs
© BioMed Central Ltd 2010
- Published: 15 March 2010
The discovery of several new structured non-coding RNAs in bacterial and archaeal genomes and metagenomes raises burning questions about their biological and biochemical functions.
See related research article by Weinberg et al.: http://genomebiology.com/2010/11/3/R31
- Secondary Structure Prediction
- Poor Discriminative Power
- Correct Secondary Structure
- Secondary Structure Helix
The compact genomes of bacteria contain 10 to 15% non-coding DNA sequences, which are transcribed into non-coding RNAs. Several classes of non-coding RNAs are small, less than 80 to 150 nucleotides, and act as post-transcriptional regulators by targeting mRNAs. Another large class of non-coding RNAs act in cis by binding structured elements in the 5' untranslated regions of mRNAs. Perhaps the best known are called riboswitches; upon binding a metabolite, the fold of the transcript is modified and this influences either the termination of transcription or the initiation of translation .
Some longer non-coding RNAs have also been detected in recent years. For example, RNAIII present in several Gram-positive bacteria is 500 nucleotides long and contains structured regions framing an open-reading frame . However, two recent papers from Ron Breaker's group increase the number of large non-coding RNAs astonishingly [3, 4]. Several new smaller non-coding RNAs are also identified. Strikingly, most of the new non-coding RNAs are structurally very complex. The complexity of some of the larger ones seems similar to that of the large ribozymes, such as the self-splicing group I and group II introns. These observations show, once again, how little we know about the microbial world: a great proportion of these new non-coding RNAs were identified in metagenomes or in environmental DNA sequences.
The search for non-coding RNAs in genomes is far from trivial . Even for homologous and functionally well characterized RNA molecules, such as the ubitquitous RNaseP or the telomerase RNA, the search cannot be reliably automated because of the large and unpredictable variation in the length of the RNA transcript, with new insertions appearing in an otherwise globally similar secondary structure. On the other hand, the de novo search for the presence of non-coding RNAs within intergenic regions is plagued by false positives because of the poor discriminative power. Various computer tools have been produced for searching for potential non-coding RNAs in genomes by exploiting the thermodynamic stabilities of the helices formed [6, 7]. The tools are generally dedicated to searching for either cis-acting RNAs (such as riboswitches) or trans-acting RNAs (such as the RNAs binding by full or partial complementarity to another RNA, either non-coding or coding).
Computer tools have been around for some time for searching RNAs on the basis of a known element of secondary structure. It has also been established several years ago that secondary structure alone is not enough for predicting non-coding RNA . The computational pipeline followed by Weinberg and coworkers [3, 4] exploits the power of comparative sequence analysis and involves sophisticated automatic techniques combined with manual intervention. The central tool used by Weinberg and coworkers [3, 4] is CMfinder, which can derive RNA motifs and secondary structures from a set of unaligned RNA sequences . However, in order to appreciate what these programs attempt to do, it is worth recalling how complex the structures of non-coding RNAs can be.
What is meant with structural complexity? The first level of folding of the transcribed RNA is the fold-back hairpin capped by a loop. Such a simple single hairpin can have profound biological effects. In bacteria, insertion of selenocysteine (a version of cysteine containing selenium rather than sulfur) occurs because the stop codon to be read as a selenocysteine codon is followed by a small hairpin. Series of hairpins can form, which, upon binding a ligand (another RNA, a protein or a metabolite), will lead to a more complex fold or to cleavage of the RNA. Structural complexity starts to appear when hairpins branch off from a hairpin, forming a three-way or multi-way junction. Naturally, further branching off of hairpins can occur within an already branched off hairpin. Because hairpins in three-dimensional space form RNA helices, which are bulky, the available space they can occupy is restricted, leading to co-axial or parallel stacking of some of them and, consequently, intricate three-dimensional architectures.
Such non-Watson-Crick pairs and the modules they form are an integral part of the tertiary structure; consequently, they are not predicted by the usual secondary structure programs that consider only Watson-Crick pairs. Correct secondary structure predictions should leave the bases that are potentially involved in non-Watson-Crick interactions as unpaired and single-stranded. Incorrect secondary structure predictions tend to predict that the bases that, in the native fold, would be forming non-Watson-Crick pairs are, instead, involved in secondary structure helices; this mis-prediction prevents the correct identification of structural elements key for the tertiary structure. Consequently, secondary structure predictions that allow for the possibility that single-stranded regions can form a known and recurrent RNA module have a higher probability of being functionally correct. Furthermore, given that such RNA modules are key elements of the tertiary structure, their presence indicates a potentially highly structured molecule.
Some striking cases are present in some secondary structures proposed for the newly reported RNAs. For example, the GOLLD (stands for Giant, Ornate, Lake- and Lactobacillales-Derived) RNA  contains two G-bulged modules, one internal loop within a hairpin, and a second loop that forms a complex junction comprising four helices. In a very unusual example, the two strands forming the G-bulged modules exchange in the sequences (69% of the observed sequences start with 5'-AAA...AGUA-3' and 18% 5'-AGUA...AAA-3'; the remaining 5% adopt a simpler purine-rich module). Another RNA, dct-1, has a cluster of four G-bulged modules positioned around a three-way junction . Interestingly, dct-1 is observed only in Dictyoglomus thermophilum, an extreme thermophile.
As discussed by Weinberg and colleagues [3, 4], several of the new RNAs could not have been discovered in the genomes of cultured bacteria known so far because such genomes do not contain the reported RNAs (except for some of the most recently sequenced genomes). Thus, the large collection of new RNAs are most probably just the tip of the iceberg, and an incredible number of still-to-be-discovered non-coding RNAs may be present in environmental sequences. The naming of the RNAs will continue to reflect the harvest of the sequences (for example, whalefall-1, Ocean-5, Soil-1 or Rhodopirellula-1).
The two recent papers [3, 4] are extremely rich in information content, with large and complete supplementary material. They present many more RNAs, some of which are new riboswitches, with several containing various structural elements, such as interactions between loops or between a loop and a single-stranded region. Here, we highlight one particular aspect of the work. In the future, much more biochemical work, tedious and time-consuming, will be necessary to characterize the functions of the non-coding RNAs, to see whether they interact with a metabolite, another RNA or a protein and participate in regulatory networks, to identify those RNAs with catalytic power and to assess how widespread they are and why they were so elusive up to now.
- Barrick JE, Breaker RR: The distributions, mechanisms, and structures of metabolite-binding riboswitches. Genome Biol. 2007, 8: R239-10.1186/gb-2007-8-11-r239.PubMedPubMed CentralView ArticleGoogle Scholar
- Benito Y, Kolb FA, Romby P, Lina G, Etienne J, Vandenesch F: Probing the structure of RNAIII, the Staphylococcus aureus agr regulatory RNA, and identification of the RNA domain involved in repression of protein A expression. RNA. 2000, 6: 668-679. 10.1017/S1355838200992550.PubMedPubMed CentralView ArticleGoogle Scholar
- Weinberg Z, Perreault J, Meyer MM, Breaker RR: Exceptional structured noncoding RNAs revealed by bacterial metagenome analysis. Nature. 2009, 462: 656-659. 10.1038/nature08586.PubMedPubMed CentralView ArticleGoogle Scholar
- Weinberg Z, Wang JX, Bogue J, Yang J, Corbino KA, Moy RH, Breaker RM: Comparative genomics reveals 104 candidate structured RNAs from bacteria, archaea and their metagenomes. Genome Biol. 2010, 11: R31-PubMedPubMed CentralView ArticleGoogle Scholar
- Menzel P, Gorodkin J, Stadler PF: The tedious task of finding homologous noncoding RNA genes. RNA. 2009, 15: 2075-2082. 10.1261/rna.1556009.PubMedPubMed CentralView ArticleGoogle Scholar
- Gorodkin J, Hofacker IL, Torarinsson E, Yao Z, Havgaard JH, Ruzzo WL: De novo prediction of structured RNAs from genomic sequences. Trends Biotechnol. 2009, 28: 9-19. 10.1016/j.tibtech.2009.09.006.PubMedView ArticleGoogle Scholar
- Yao Z, Weinberg Z, Ruzzo WL: CMfinder - -a covariance model based RNA motif finding algorithm. Bioinformatics. 2006, 22: 445-452. 10.1093/bioinformatics/btk008.PubMedView ArticleGoogle Scholar
- Rivas E, Eddy SR: Secondary structure alone is generally not statistically significant for the detection of noncoding RNAs. Bioinformatics. 2000, 16: 583-605. 10.1093/bioinformatics/16.7.583.PubMedView ArticleGoogle Scholar
- Leontis NB, Stombaugh J, Westhof E: Motif prediction in ribosomal RNAs Lessons and prospects for automated motif prediction in homologous RNA molecules. Biochimie. 2002, 84: 961-973. 10.1016/S0300-9084(02)01463-3.PubMedView ArticleGoogle Scholar
- Schuwirth BS, Borovinskaya MA, Hau CW, Zhang W, Vila-Sanjurjo A, Holton JM, Cate JH: Structures of the bacterial ribosome at 3.5 A resolution. Science. 2005, 310: 827-834. 10.1126/science.1117230.PubMedView ArticleGoogle Scholar