Structure-seq2 probing of RNA structure upon amino acid starvation reveals both known and novel RNA switches in Bacillus subtilis

RNA Structure: Past, Future, and Gene Therapy Applications

First believed to be a simple intermediary between the information encoded in deoxyribonucleic acid and that functionally displayed in proteins, ribonucleic acid (RNA) is now known to have many functions through its abundance and intricate, ubiquitous, diverse, and dynamic structure. About 70-90% of the human genome is transcribed into protein-coding and noncoding RNAs as main determinants along with regulatory sequences of cellular to populational biological diversity. From the nucleotide sequence or primary structure, through Watson-Crick pairing self-folding or secondary structure, to compaction via longer distance Watson-Crick and non-Watson-Crick interactions or tertiary structure, and interactions with RNA or other biopolymers or quaternary structure, or with metabolites and biomolecules or quinary structure, RNA structure plays a critical role in RNA's lifecycle from transcription to decay and many cellular processes. In contrast to the success of 3-dimensional protein structure prediction using AlphaFold, RNA tertiary and beyond structures prediction remains challenging. However, approaches involving machine learning and artificial intelligence, sequencing of RNA and its modifications, and structural analyses at the single-cell and intact tissue levels, among others, provide an optimistic outlook for the continued development and refinement of RNA-based applications. Here, we highlight those in gene therapy.

In vivo structure probing of RNA in Archaea: novel insights into the ribosome structure of Methanosarcina acetivorans

Structure probing combined with next-generation sequencing (NGS) has provided novel insights into RNA structure-function relationships. To date, such studies have focused largely on bacteria and eukaryotes, with little attention given to the third domain of life, archaea. Furthermore, functional RNAs have not been extensively studied in archaea, leaving open questions about RNA structure and function within this domain of life. With archaeal species being diverse and having many similarities to both bacteria and eukaryotes, the archaea domain has the potential to be an evolutionary bridge. In this study, we introduce a method for probing RNA structure in vivo in the archaea domain of life. We investigated the structure of ribosomal RNA (rRNA) from Methanosarcina acetivorans, a well-studied anaerobic archaeal species, grown with either methanol or acetate. After probing the RNA in vivo with dimethyl sulfate (DMS), Structure-seq2 libraries were generated, sequenced, and analyzed. We mapped the reactivity of DMS onto the secondary structure of the ribosome, which we determined independently with comparative analysis, and confirmed the accuracy of DMS probing in M. acetivorans Accessibility of the rRNA to DMS in the two carbon sources was found to be quite similar, although some differences were found. Overall, this study establishes the Structure-seq2 pipeline in the archaea domain of life and informs about ribosomal structure within M. acetivorans.

Transcriptome-wide probing reveals RNA thermometers that regulate translation of glycerol permease genes in Bacillus subtilis

RNA structure regulates bacterial gene expression by several distinct mechanisms via environmental and cellular stimuli, one of which is temperature. While some genome-wide studies have focused on heat shock treatments and the subsequent transcriptomic changes, soil bacteria are less likely to experience such rapid and extreme temperature changes. Though RNA thermometers (RNATs) have been found in 5' untranslated leader regions (5' UTRs) of heat shock and virulence-associated genes, this RNA-controlled mechanism could regulate other genes as well. Using Structure-seq2 and the chemical probe dimethyl sulfate (DMS) at four growth temperatures ranging from 23°C to 42°C, we captured a dynamic response of the Bacillus subtilis transcriptome to temperature. Our transcriptome-wide results show RNA structural changes across all four temperatures and reveal nonmonotonic reactivity trends with increasing temperature. Then, focusing on subregions likely to contain regulatory RNAs, we examined 5' UTRs to identify large, local reactivity changes. This approach led to the discovery of RNATs that control the expression of glpF (glycerol permease) and glpT (glycerol-3-phosphate permease); expression of both genes increased with increased temperature. Results with mutant RNATs indicate that both genes are controlled at the translational level. Increased import of glycerols at high temperatures could provide thermoprotection to proteins.

Structure-seq of tRNAs and other short RNAs in droplets and in vivo

There is a multitude of small (<100nt) RNAs that serve diverse functional roles in biology. Key amongst these is transfer RNA (tRNA), which is among the most ancient RNAs and is part of the translational apparatus in every domain of life. Transfer RNAs are also the most heavily modified class of RNAs. They are essential and their misregulation, due to mutated sequences or loss of modification, can lead to disease. Because of the severe phenotypes associated with mitochondrial tRNA defects in particular, the desire to deliver repaired tRNAs via droplets such as lipid nanoparticles or other compartments is an active area of research. Here we describe how to use our tRNA Structure-seq method to study tRNAs and other small RNAs in two different biologically relevant contexts, peptide-rich droplets and in vivo.

Rock, scissors, paper: How RNA structure informs function

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

Comprehensive transcription terminator atlas for Bacillus subtilis

The transcriptome-wide contributions of Rho-dependent and intrinsic (Rho-independent) transcription termination mechanisms in bacteria are unclear. By sequencing released transcripts in a wild-type strain and strains containing deficiencies in NusA, NusG and/or Rho (10 strains), we produced an atlas of terminators for the model Gram-positive bacterium Bacillus subtilis. We found that NusA and NusG stimulate 77% and 19% of all intrinsic terminators, respectively, and that both proteins participate in Rho-dependent termination. We also show that Rho stimulates termination at 10% of the intrinsic terminators in vivo. We recapitulated Rho-stimulated intrinsic termination at 5 terminators in vitro and found that Rho requires the KOW domain of NusG to stimulate this process at one of these terminators. Computational analyses of our atlas using RNAstructure, MEME suite and DiffLogo, combined with in vitro transcription experiments, revealed that Rho stimulates intrinsic terminators with weak hairpins and/or U-rich tracts by remodelling the RNA upstream of the intrinsic terminator to prevent the formation of RNA structures that could otherwise compete with the terminator hairpin. We also identified 56 putative examples of 'hybrid Rho-dependent termination', wherein classical Rho-dependent termination occurs after readthrough of a Rho-stimulated intrinsic terminator.

The Second Class of Tetrahydrofolate (THF-II) Riboswitches Recognizes the Tetrahydrofolic Acid Ligand via Local Conformation Changes

Riboswitches are regulatory noncoding RNAs found in bacteria, fungi and plants, that modulate gene expressions through structural changes in response to ligand binding. Understanding how ligands interact with riboswitches in solution can shed light on the molecular mechanisms of this ancient regulators. Previous studies showed that riboswitches undergo global conformation changes in response to ligand binding to relay information. Here, we report conformation switching models of the recently discovered tetrahydrofolic acid-responsive second class of tetrahydrofolate (THF-II) riboswitches in response to ligand binding. Using a combination of selective 2′-hydroxyl acylation, analyzed by primer extension (SHAPE) assay, 3D modeling and small-angle X-ray scattering (SAXS), we found that the ligand specifically recognizes and reshapes the THF-II riboswitch loop regions, but does not affect the stability of the P3 helix. Our results show that the THF-II riboswitch undergoes only local conformation changes in response to ligand binding, rearranging the Loop1-P3-Loop2 region and rotating Loop1 from a ~120° angle to a ~75° angle. This distinct conformation changes suggest a unique regulatory mechanism of the THF-II riboswitch, previously unseen in other riboswitches. Our findings may contribute to the fields of RNA sensors and drug design.

Nearest neighbor rules for RNA helix folding thermodynamics: improved end effects

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5′GGUC3′ paired to 3′CUGG5′ separately from other 5′GU3′/3′UG5′ stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37°C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model.

Analysis of mRNA Decay Intermediates in Bacillus subtilis 3' Exoribonuclease and RNA Helicase Mutant Strains

The Bacillus subtilis genome encodes four 3' exoribonucleases: polynucleotide phosphorylase (PNPase), RNase R, RNase PH, and YhaM. Previous work showed that PNPase, encoded by the pnpA gene, is the major 3' exonuclease involved in mRNA turnover; in a pnpA deletion strain, numerous mRNA decay intermediates accumulate. Whether B. subtilis mRNA decay occurs in the context of a degradosome complex is controversial. In this study, global mapping of mRNA decay intermediate 3' ends within coding sequences was performed in strains that were either deleted for or had an inactivating point mutation in the pnpA gene. The patterns of 3'-end accumulation in these strains were highly similar, which may have implications for the role of a degradosome in mRNA decay. A comparison with mapped 3' ends in a strain lacking CshA, the major RNA helicase, indicated that many mRNAs require both PNPase and CshA for efficient decay. Transcriptome sequencing (RNA-seq) analysis of strains lacking RNase R suggested that this enzyme did not play a major role in mRNA turnover in the wild-type strain. Strains were constructed that contained only one of the four known 3' exoribonucleases. When RNase R was the only 3' exonuclease present, it was able to degrade a model mRNA efficiently, showing processive decay even through a strong stem-loop structure that inhibits PNPase processivity. Strains containing only RNase PH or only YhaM were also insensitive to this RNA secondary structure, suggesting the existence of another, as-yet-unidentified, 3' exoribonuclease. IMPORTANCE The ability to rapidly change bacterial gene expression programs in response to environmental conditions is highly dependent on the efficient turnover of mRNA. While much is known about the regulation of gene expression at the transcriptional and translational levels, much less is known about the intermediate step of mRNA decay. Here, we mapped the 3' ends of mRNA decay intermediates in strains that were missing the major 3' exoribonuclease PNPase or the RNA helicase CshA. We also assessed the roles of three other B. subtilis 3' exonucleases in the mRNA decay process. The data confirm the primary role of PNPase in mRNA turnover and suggest the involvement of one or more unknown RNases.

Perspective: a stirring role for metabolism in cells

Based on recent findings indicating that metabolism might be governed by a limit on the rate at which cells can dissipate Gibbs energy, in this Perspective, we propose a new mechanism of how metabolic activity could globally regulate biomolecular processes in a cell. Specifically, we postulate that Gibbs energy released in metabolic reactions is used to perform work, allowing enzymes to self-propel or to break free from supramolecular structures. This catalysis-induced enzyme movement will result in increased intracellular motion, which in turn can compromise biomolecular functions. Once the increased intracellular motion has a detrimental effect on regulatory mechanisms, this will establish a feedback mechanism on metabolic activity, and result in the observed thermodynamic limit. While this proposed explanation for the identified upper rate limit on cellular Gibbs energy dissipation rate awaits experimental validation, it offers an intriguing perspective of how metabolic activity can globally affect biomolecular functions and will hopefully spark new research.

Small RNA Profiling in Mycobacterium Provides Insights Into Stress Adaptability

Mycobacteria encounter a number of environmental changes during infection and respond using different mechanisms. Small RNA (sRNA) is a post-transcriptionally regulatory system for gene functions and has been investigated in many other bacteria. This study used Mycobacterium tuberculosis and Mycobacterium bovis Bacillus Calmette-Guérin (BCG) infection models and sequenced whole bacterial RNAs before and after host cell infection. A comparison of differentially expressed sRNAs using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) and target prediction was carried out. Six pathogenically relevant stress conditions, growth rate, and morphology were used to screen and identify sRNAs. From these data, a subset of sRNAs was differentially expressed in multiple infection groups and stress conditions. Many were found associated with lipid metabolism. Among them, ncBCG427 was significantly downregulated when BCG entered into macrophages and was associated with increased biofilm formation. The reduction of virulence possibility depends on regulating lipid metabolism.

NusG is an intrinsic transcription termination factor that stimulates motility and coordinates gene expression with NusA

NusA and NusG are transcription factors that stimulate RNA polymerase pausing in Bacillus subtilis. While NusA was known to function as an intrinsic termination factor in B. subtilis, the role of NusG in this process was unknown. To examine the individual and combinatorial roles that NusA and NusG play in intrinsic termination, Term-seq was conducted in wild type, NusA depletion, ΔnusG, and NusA depletion ΔnusG strains. We determined that NusG functions as an intrinsic termination factor that works alone and cooperatively with NusA to facilitate termination at 88% of the 1400 identified intrinsic terminators. Our results indicate that NusG stimulates a sequence-specific pause that assists in the completion of suboptimal terminator hairpins with weak terminal A-U and G-U base pairs at the bottom of the stem. Loss of NusA and NusG leads to global misregulation of gene expression and loss of NusG results in flagella and swimming motility defects.