Sampled ensemble neutrality as a feature to classify potential structured RNAs.

BACKGROUND: Structured RNAs have many biological functions ranging from catalysis of chemical reactions to gene regulation. Yet, many homologous structured RNAs display most of their conservation at the secondary or tertiary structure level. As a result, strategies for structured RNA discovery rely heavily on identification of sequences sharing a common stable secondary structure. However, correctly distinguishing structured RNAs from surrounding genomic sequence remains challenging, especially during de novo discovery. RNA also has a long history as a computational model for evolution due to the direct link between genotype (sequence) and phenotype (structure). From these studies it is clear that evolved RNA structures, like protein structures, can be considered robust to point mutations. In this context, an RNA sequence is considered robust if its neutrality (extent to which single mutant neighbors maintain the same secondary structure) is greater than that expected for an artificial sequence with the same minimum free energy structure. RESULTS: In this work, we bring concepts from evolutionary biology to bear on the structured RNA de novo discovery process. We hypothesize that alignments corresponding to structured RNAs should consist of neutral sequences. We evaluate several measures of neutrality for their ability to distinguish between alignments of structured RNA sequences drawn from Rfam and various decoy alignments. We also introduce a new measure of RNA structural neutrality, the structure ensemble neutrality (SEN). SEN seeks to increase the biological relevance of existing neutrality measures in two ways. First, it uses information from an alignment of homologous sequences to identify a conserved biologically relevant structure for comparison. Second, it only counts base-pairs of the original structure that are absent in the comparison structure and does not penalize the formation of additional base-pairs. CONCLUSION: We find that several measures of neutrality are effective at separating structured RNAs from decoy sequences, including both shuffled alignments and flanking genomic sequence. Furthermore, as an independent feature classifier to identify structured RNAs, SEN yields comparable performance to current approaches that consider a variety of features including stability and sequence identity. Finally, SEN outperforms other measures of neutrality at detecting mutational robustness in bacterial regulatory RNA structures.

Discovery and validation of novel and distinct RNA regulators for ribosomal protein S15 in diverse bacterial phyla.

BACKGROUND: Autogenous cis-regulators of ribosomal protein synthesis play a critical role in maintaining the stoichiometry of ribosome components. Structured portions within an mRNA transcript typically interact with specific ribosomal proteins to prevent expression of the entire operon, thus balancing levels of ribosomal proteins across transcriptional units. Three distinct RNA structures from different bacterial phyla have demonstrated interactions with S15 to regulate gene expression; however, these RNAs are distributed across a small fraction of bacterial diversity. RESULTS: We used comparative genomics in combination with analysis of existing transcriptomic data to identify three novel putative RNA structures associated with the S15 coding region in microbial genomes. These structures are completely distinct from those previously published and encompass potential regulatory regions including ribosome-binding sites. To validate the biological relevance of our findings, we demonstrate that an example of the Alphaproteobacterial RNA from Rhizobium radiobacter specifically interacts with S15 in vitro, and allows in vivo regulation of gene expression in an E. coli reporter system. In addition, structural probing and nuclease protection assays confirm the predicted secondary structure and indicate nucleotides required for protein interaction. CONCLUSIONS: This work illustrates the importance of integrating comparative genomic and transcriptomic approaches during de novo ncRNA identification and reveals a diversity of distinct natural RNA regulators that support analogous biological functions. Furthermore, this work indicates that many additional uncharacterized RNA regulators likely exist within bacterial genomes and that the plasticity of RNA structure allows unique, and likely independently derived, solutions to the same biological problem.