RNA motif search with data-driven element ordering

Robo-Therm, a pipeline to RNA thermometer discovery and validation

RNA thermometers are highly structured noncoding RNAs located in the 5'-untranslated regions (UTRs) of genes that regulate expression by undergoing conformational changes in response to temperature. The discovery of RNA thermometers through bioinformatics is difficult because there is little sequence conservation among their structural elements. Thus, the abundance of these thermosensitive regulatory structures remains unclear. Herein, to advance the discovery and validation of RNA thermometers, we developed Robo-Therm, a pipeline that combines an adaptive and user-friendly in silico motif search with a well-established reporter system. Through our application of Robo-Therm, we discovered two novel RNA thermometers in bacterial and bacteriophage genomes found in the human gut. One of these thermometers is present in the 5'-UTR of a gene that codes for σ 70 RNA polymerase subunit in the bacteria Mediterraneibacter gnavus and Bacteroides pectinophilus, and in the bacteriophage Caudoviricetes, which infects B. pectinophilus The other thermometer is in the 5'-UTR of a tetracycline resistance gene (tetR) in the intestinal bacteria Escherichia coli and Shigella flexneri Our Robo-Therm pipeline can be applied to discover multiple RNA thermometers across various genomes.

Identification of HDV-like theta ribozymes involved in tRNA-based recoding of gut bacteriophages

Trillions of microorganisms, collectively known as the microbiome, inhabit our bodies with the gut microbiome being of particular interest in biomedical research. Bacteriophages, the dominant virome constituents, can utilize suppressor tRNAs to switch to alternative genetic codes (e.g., the UAG stop-codon is reassigned to glutamine) while infecting hosts with the standard bacterial code. However, what triggers this switch and how the bacteriophage manipulates its host is poorly understood. Here, we report the discovery of a subgroup of minimal hepatitis delta virus (HDV)-like ribozymes - theta ribozymes - potentially involved in the code switch leading to the expression of recoded lysis and structural phage genes. We demonstrate their HDV-like self-scission behavior in vitro and find them in an unreported context often located with their cleavage site adjacent to tRNAs, indicating a role in viral tRNA maturation and/or regulation. Every fifth associated tRNA is a suppressor tRNA, further strengthening our hypothesis. The vast abundance of tRNA-associated theta ribozymes - we provide 1753 unique examples - highlights the importance of small ribozymes as an alternative to large enzymes that usually process tRNA 3'-ends. Our discovery expands the short list of biological functions of small HDV-like ribozymes and introduces a previously unknown player likely involved in the code switch of certain recoded gut bacteriophages.

M. Passalacqua L, Abdelsayed MM. Characterization of a FourU RNA Thermometer in the 5' Untranslated Region of Autolysin Gene blyA in the Bacillus subtilis 168 Prophage SPβ

RNA thermometers are noncoding RNA structures located in the 5' untranslated regions (UTRs) of genes that regulate gene expression through temperature-dependent conformational changes. The fourU class of RNA thermometers contains a specific motif in which four consecutive uracil nucleotides are predicted to base pair with the Shine-Dalgarno (SD) sequence in a stem. We employed a bioinformatic search to discover a fourU RNA thermometer in the 5'-UTR of the blyA gene of the Bacillus subtilis phage SPβc2, a bacteriophage that infects B. subtilis 168. blyA encodes an autolysin enzyme, N-acetylmuramoyl-l-alanine amidase, which is involved in the lytic life cycle of the SPβ prophage. We have biochemically validated the predicted RNA thermometer in the 5'-UTR of the blyA gene. Our study suggests that RNA thermometers may play an underappreciated yet critical role in the lytic life cycle of bacteriophages.

High-throughput methods in aptamer discovery and analysis

Aptamers are small, functional nucleic acids that bind a variety of targets, often with high specificity and affinity. Genomic aptamers constitute the ligand-binding domains of riboswitches, whereas synthetic aptamers find applications as diagnostic and therapeutic tools, and as ligand-binding domains of regulatory RNAs in synthetic biology. Discovery and characterization of aptamers has been limited by a lack of high-throughput approaches that uncover the target-binding domains and the biochemical properties of individual sequences. With the advent of high-throughput sequencing, large-scale analysis of in vitro selected populations of aptamers (and catalytic nucleic acids, such as ribozymes and DNAzmes) became possible. In recent years the development of new experimental approaches and software tools has led to significant streamlining of the selection-pool analysis. This article provides an overview of post-selection data analysis and describes high-throughput methods that facilitate rapid discovery and biochemical characterization of aptamers.

Kinetic Parameters of trans Scission by Extended HDV-like Ribozymes and the Prospect for the Discovery of Genomic trans-Cleaving RNAs

Hepatitis delta virus (HDV)-like ribozymes are self-cleaving catalytic RNAs with a widespread distribution in nature and biological roles ranging from self-scission during rolling-circle replication in viroids to co-transcriptional processing of eukaryotic retrotransposons, among others. The ribozymes fold into a double pseudoknot with a common catalytic core motif and highly variable peripheral domains. Like other self-cleaving ribozymes, HDV-like ribozymes can be converted into trans-acting catalytic RNAs by bisecting the self-cleaving variants at non-essential loops. Here we explore the trans-cleaving activity of ribozymes derived from the largest examples of the ribozymes (drz-Agam-2 family), which contain an extended domain between the substrate strand and the rest of the RNA. When this peripheral domain is bisected at its distal end, the substrate strand is recognized through two helices, rather than just one 7 bp helix common among the HDV ribozymes, resulting in stronger binding and increased sequence specificity. Kinetic characterization of the extended trans-cleaving ribozyme revealed an efficient trans-cleaving system with a surprisingly high K_M', supporting a model that includes a recently proposed activation barrier related to the assembly of the catalytically competent ribozyme. The ribozymes also exhibit a very long k_off for the products (∼2 weeks), resulting in a trade-off between sequence specificity and turnover. Finally, structure-based searches for the catalytic cores of these ribozymes in the genome of the mosquito Anopheles gambiae, combined with sequence searches for their putative substrates, revealed two potential ribozyme-substrate pairs that may represent examples of natural trans-cleaving ribozymes.