Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo

RNA has a dual role as an informational molecule and a direct effector of biological tasks. The latter function is enabled by RNA’s ability to adopt complex secondary and tertiary folds and thus has motivated extensive computational and experimental efforts for determining RNA structures. Existing approaches for evaluating RNA structure have been largely limited to in vitro systems, yet the thermodynamic forces which drive RNA folding in vitro may not be sufficient to predict stable RNA structures in vivo. Indeed, the presence of RNA-binding proteins and ATP-dependent helicases can influence which structures are present inside cells. Here we present an approach for globally monitoring RNA structure in native conditions in vivo with single-nucleotide precision. This method is based on in vivo modification with dimethyl sulphate (DMS), which reacts with unpaired adenine and cytosine residues, followed by deep sequencing to monitor modifications. Our data from yeast and mammalian cells are in excellent agreement with known messenger RNA structures and with the high-resolution crystal structure of the Saccharomyces cerevisiae ribosome. Comparison between in vivo and in vitro data reveals that in rapidly dividing cells there are vastly fewer structured mRNA regions in vivo than in vitro. Even thermostable RNA structures are often denatured in cells, highlighting the importance of cellular processes in regulating RNA structure. Indeed, analysis of mRNA structure under ATP-depleted conditions in yeast shows that energy-dependent processes strongly contribute to the predominantly unfolded state of mRNAs inside cells. Our studies broadly enable the functional analysis of physiological RNA structures and reveal that, in contrast to the Anfinsen view of protein folding whereby the structure formed is the most thermodynamically favourable, thermodynamics have an incomplete role in determining mRNA structure in vivo

The helicase Ded1p controls use of near-cognate translation initiation codons in 5' UTRs.

The conserved and essential DEAD-box RNA helicase Ded1p from yeast and its mammalian orthologue DDX3 are critical for the initiation of translation1. Mutations in DDX3 are linked to tumorigenesis2-4 and intellectual disability5, and the enzyme is targeted by a range of viruses6. How Ded1p and its orthologues engage RNAs during the initiation of translation is unknown. Here we show, by integrating transcriptome-wide analyses of translation, RNA structure and Ded1p-RNA binding, that the effects of Ded1p on the initiation of translation are connected to near-cognate initiation codons in 5' untranslated regions. Ded1p associates with the translation pre-initiation complex at the mRNA entry channel and repressing the activity of Ded1p leads to the accumulation of RNA structure in 5' untranslated regions, the initiation of translation from near-cognate start codons immediately upstream of these structures and decreased protein synthesis from the corresponding main open reading frames. The data reveal a program for the regulation of translation that links Ded1p, the activation of near-cognate start codons and mRNA structure. This program has a role in meiosis, in which a marked decrease in the levels of Ded1p is accompanied by the activation of the alternative translation initiation sites that are seen when the activity of Ded1p is repressed. Our observations indicate that Ded1p affects translation initiation by controlling the use of near-cognate initiation codons that are proximal to mRNA structure in 5' untranslated regions

Novel methods for the quantitative determination of RNA folding on a genome-wide scale and in a targeted manner

The coupling of structure-specific in vivo chemical modification to next-generation sequencing has revolutionized the study of RNA secondary structure in living cells. Nonetheless, current approaches are limited by the detection of structural information as reverse transcriptase truncation products, which creates biases, blurs signal from heterogeneous structures, and is restricted to abundant RNAs. Here we present dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq), which encodes DMS modifications as mismatches using a thermostable group II intron reverse transcriptase (TGIRT). DMS-MaPseq yields a high signal-to-noise ratio, can report multiple structural features for each molecule, and allows genome-wide studies as well as focused investigations of low abundance RNAs. We apply DMS-MaPseq to Drosophila melanogaster ovaries—the first experimental analysis of RNA structure in an animal tissue—and demonstrate its utility in the discovery of a functional RNA structure involved in the non-canonical GUG translation initiation of the human FXR2 mRNA. Additionally, we use DMS-MaPseq to compare the in vivo structure of messages in their pre-mRNA and mature forms. These applications illustrate DMS-MaPseq’s capacity to dramatically expand our ability to monitor RNA structure in vivo

MegZ03_index13_AGTCAA_1MISS_lowerM_wildtype_polyA_PLUSstrand_COUNTS

This is .WIG file (processed from .fastq files) for DMS-seq data in S. pombe cells (plus strand) (generated by Meghan Zubradt

Engineering ER-stress dependent non-conventional mRNA splicing

The endoplasmic reticulum (ER) protein folding capacity is balanced with the protein folding burden to prevent accumulation of un- or misfolded proteins. The ER membrane-resident kinase/RNase Ire1 maintains ER protein homeostasis through two fundamentally distinct processes. First, Ire1 can initiate a transcriptional response through a non-conventional mRNA splicing reaction to increase the ER folding capacity. Second, Ire1 can decrease the ER folding burden through selective mRNA decay. In Saccharomyces cerevisiae and Schizosaccharomyces pombe, the two Ire1 functions have been evolutionarily separated. Here, we show that the respective Ire1 orthologs have become specialized for their functional outputs by divergence of their RNase specificities. In addition, RNA structural features separate the splicing substrates from the decay substrates. Using these insights, we engineered an S. pombe Ire1 cleavage substrate into a splicing substrate, which confers S. pombe with both Ire1 functional outputs.ISSN:2050-084

DMS-MaPseq for genome-wide or targeted RNA structure probing in vivo

Coupling of structure-specific in vivo chemical modification to next-generation sequencing is transforming RNA secondary structure studies in living cells. The dominant strategy for detecting in vivo chemical modifications uses reverse transcriptase truncation products, which introduce biases and necessitate population-average assessments of RNA structure. Here we present dimethyl sulfate (DMS) mutational profiling with sequencing (DMS-MaPseq), which encodes DMS modifications as mismatches using a thermostable group II intron reverse transcriptase. DMS-MaPseq yields a high signal-to-noise ratio, can report multiple structural features per molecule, and allows both genome-wide studies and focused in vivo investigations of even low-abundance RNAs. We apply DMS-MaPseq for the first analysis of RNA structure within an animal tissue and to identify a functional structure involved in noncanonical translation initiation. Additionally, we use DMS-MaPseq to compare the in vivo structure of pre-mRNAs with their mature isoforms. These applications illustrate DMS-MaPseq's capacity to dramatically expand in vivo analysis of RNA structure