Genome-wide discovery of structured noncoding RNAs in bacteria.

BACKGROUND: Structured noncoding RNAs (ncRNAs) play essential roles in many biological processes such as gene regulation, signaling, RNA processing, and protein synthesis. Among the most common groups of ncRNAs in bacteria are riboswitches. These cis-regulatory, metabolite-binding RNAs are present in many species where they regulate various metabolic and signaling pathways. Collectively, there are likely to be hundreds of novel riboswitch classes that remain hidden in the bacterial genomes that have already been sequenced, and potentially thousands of classes distributed among various other species in the biosphere. The vast majority of these undiscovered classes are proposed to be exceedingly rare, and so current bioinformatics search techniques are reaching their limits for differentiating between true riboswitch candidates and false positives. RESULTS: Herein, we exploit a computational search pipeline that can efficiently identify intergenic regions most likely to encode structured ncRNAs. Application of this method to five bacterial genomes yielded nearly 70 novel genetic elements including 30 novel candidate ncRNA motifs. Among the riboswitch candidates identified is an RNA motif involved in the regulation of thiamin biosynthesis. CONCLUSIONS: Analysis of other genomes will undoubtedly lead to the discovery of many additional novel structured ncRNAs, and provide insight into the range of riboswitches and other kinds of ncRNAs remaining to be discovered in bacteria and archaea.

Gene regulation by a glycine riboswitch singlet uses a finely tuned energetic landscape for helical switching.

Riboswitches contain structured aptamer domains that, upon ligand binding, facilitate helical switching in their downstream expression platforms to alter gene expression. To fully dissect how riboswitches function requires a better understanding of the energetic landscape for helical switching. Here, we report a sequencing-based high-throughput assay for monitoring in vitro transcription termination and use it to simultaneously characterize the functional effects of all 522 single point mutants of a glycine riboswitch type-1 singlet. Mutations throughout the riboswitch cause ligand-dependent defects, but only mutations within the terminator hairpin alter readthrough efficiencies in the absence of ligand. A comprehensive analysis of the expression platform reveals that ligand binding stabilizes the antiterminator by just 2-3 kcal/mol, indicating that the competing expression platform helices must be extremely close in energy to elicit a significant ligand-dependent response. These results demonstrate that gene regulation by this riboswitch is highly constrained by the energetics of ligand binding and conformational switching. These findings exemplify the energetic parameters of RNA conformational rearrangements driven by binding events.

Tandem riboswitches form a natural Boolean logic gate to control purine metabolism in bacteria.

Gene control systems sometimes interpret multiple signals to set the expression levels of the genes they regulate. In rare instances, ligand-binding riboswitch aptamers form tandem arrangements to approximate the function of specific two-input Boolean logic gates. Here, we report the discovery of riboswitch aptamers for phosphoribosyl pyrophosphate (PRPP) that naturally exist either in singlet arrangements, or occur in tandem with guanine aptamers. Tandem guanine-PRPP aptamers can bind the target ligands, either independently or in combination, to approximate the function expected for an IMPLY Boolean logic gate to regulate transcription of messenger RNAs for de novo purine biosynthesis in bacteria. The existence of sophisticated all-RNA regulatory systems that sense two ancient ribonucleotide derivatives to control synthesis of RNA molecules supports the hypothesis that RNA World organisms could have managed a complex metabolic state without the assistance of protein regulatory factors.

Challenges of ligand identification for the second wave of orphan riboswitch candidates.

Orphan riboswitch candidates are noncoding RNA motifs whose representatives are believed to function as genetic regulatory elements, but whose target ligands have yet to be identified. The study of certain orphans, particularly classes that have resisted experimental validation for many years, has led to the discovery of important biological pathways and processes once their ligands were identified. Previously, we highlighted details for four of the most common and intriguing orphan riboswitch candidates. This facilitated the validation of riboswitches for the signaling molecules c-di-AMP, ZTP, and ppGpp, the metal ion Mn2+, and the metabolites guanidine and PRPP. Such studies also yield useful linkages between the ligands sensed by the riboswitches and numerous biochemical pathways. In the current report, we describe the known characteristics of 30 distinct classes of orphan riboswitch candidates - some of which have remained unsolved for over a decade. We also discuss the prospects for uncovering novel biological insights via focused studies on these RNAs. Lastly, we make recommendations for experimental objectives along the path to finding ligands for these mysterious RNAs.

Riboswitch diversity and distribution.

Riboswitches are commonly used by bacteria to detect a variety of metabolites and ions to regulate gene expression. To date, nearly 40 different classes of riboswitches have been discovered, experimentally validated, and modeled at atomic resolution in complex with their cognate ligands. The research findings produced since the first riboswitch validation reports in 2002 reveal that these noncoding RNA domains exploit many different structural features to create binding pockets that are extremely selective for their target ligands. Some riboswitch classes are very common and are present in bacteria from nearly all lineages, whereas others are exceedingly rare and appear in only a few species whose DNA has been sequenced. Presented herein are the consensus sequences, structural models, and phylogenetic distributions for all validated riboswitch classes. Based on our findings, we predict that there are potentially many thousands of distinct bacterial riboswitch classes remaining to be discovered, but that the rarity of individual undiscovered classes will make it increasingly difficult to find additional examples of this RNA-based sensory and gene control mechanism.

The yjdF riboswitch candidate regulates gene expression by binding diverse azaaromatic compounds.

The yjdF motif RNA is an orphan riboswitch candidate that almost exclusively associates with the yjdF protein-coding gene in many bacteria. The function of the YjdF protein is unknown, which has made speculation regarding the natural ligand for this putative riboswitch unusually challenging. By using a structure-probing assay for ligand binding, we found that a surprisingly broad diversity of nitrogen-containing aromatic heterocycles, or "azaaromatics," trigger near-identical changes in the structures adopted by representative yjdF motif RNAs. Regions of the RNA that undergo ligand-induced structural modulation reside primarily in portions of the putative aptamer region that are highly conserved in nucleotide sequence, as is typical for riboswitches. Some azaaromatic molecules are bound by the RNA with nanomolar dissociation constants, and a subset of these ligands activate riboswitch-mediated gene expression in cells. Furthermore, genetic elements most commonly adjacent to the yjdF motif RNA or to the yjdF protein-coding region are homologous to protein regulators implicated in mitigating the toxic effects of diverse phenolic acids or polycyclic compounds. Although the precise type of natural ligand sensed by yjdF motif RNAs remains unknown, our findings suggest that this riboswitch class might serve as part of a genetic response system to toxic or signaling compounds with chemical structures similar to azaaromatics.

Control of bacterial exoelectrogenesis by c-AMP-GMP.

Major changes in bacterial physiology including biofilm and spore formation involve signaling by the cyclic dinucleotides c-di-GMP and c-di-AMP. Recently, another second messenger dinucleotide, c-AMP-GMP, was found to control chemotaxis and colonization by Vibrio cholerae. We have identified a superregulon of genes controlled by c-AMP-GMP in numerous Deltaproteobacteria, including Geobacter species that use extracellular insoluble metal oxides as terminal electron acceptors. This exoelectrogenic process has been studied for its possible utility in energy production and bioremediation. Many genes involved in adhesion, pilin formation, and others that are important for exoelectrogenesis are controlled by members of a variant riboswitch class that selectively bind c-AMP-GMP. These RNAs constitute, to our knowledge, the first known specific receptors for c-AMP-GMP and reveal that this molecule is used by many bacteria to control specialized physiological processes.

COBRA-LIKE2, a member of the glycosylphosphatidylinositol-anchored COBRA-LIKE family, plays a role in cellulose deposition in arabidopsis seed coat mucilage secretory cells.

Differentiation of the maternally derived seed coat epidermal cells into mucilage secretory cells is a common adaptation in angiosperms. Recent studies identified cellulose as an important component of seed mucilage in various species. Cellulose is deposited as a set of rays that radiate from the seed upon mucilage extrusion, serving to anchor the pectic component of seed mucilage to the seed surface. Using transcriptome data encompassing the course of seed development, we identified COBRA-LIKE2 (COBL2), a member of the glycosylphosphatidylinositol-anchored COBRA-LIKE gene family in Arabidopsis (Arabidopsis thaliana), as coexpressed with other genes involved in cellulose deposition in mucilage secretory cells. Disruption of the COBL2 gene results in substantial reduction in the rays of cellulose present in seed mucilage, along with an increased solubility of the pectic component of the mucilage. Light birefringence demonstrates a substantial decrease in crystalline cellulose deposition into the cellulosic rays of the cobl2 mutants. Moreover, crystalline cellulose deposition into the radial cell walls and the columella appears substantially compromised, as demonstrated by scanning electron microscopy and in situ quantification of light birefringence. Overall, the cobl2 mutants display about 40% reduction in whole-seed crystalline cellulose content compared with the wild type. These data establish that COBL2 plays a role in the deposition of crystalline cellulose into various secondary cell wall structures during seed coat epidermal cell differentiation.