Frameworks for large-scale RNA structure profiling in transcriptomes and disease

In addition to their role as intermediaries on the route to protein synthesis, RNA molecules have long been known to base-pair into complex structures that serve specific functions. Some structured RNAs play pathogenic roles, especially in viral illnesses and repeat-expansion disorders, and disease-associated RNA structures are potential therapeutic targets. SHAPE is a well-established chemical probing strategy to interrogate RNA flexibility and obtain high-quality structure models. The recent development of an unbiased experimental approach that allows SHAPE to characterize populations of diverse RNAs using massively parallel sequencing presented a challenging data analysis problem. In this work, I apply SHAPE to study the relevance of huntingtin mRNA structure to Huntington's disease and discover that a classical CAG hairpin is likely absent or short in healthy-length transcripts. The formation of this hairpin correlates with increasing repeat length, which is a predictor of disease severity. I develop a fully-automated data analysis pipeline allowing for the extension of the SHAPE strategy to larger scales using mutational profiling (MaP), an approach that was applied to identify highly-structured elements within an HIV-1 genome. I further provide a pilot analysis of a bacterial transcriptome MaP dataset obtained in a single experiment, demonstrate the nucleotide accuracy of MaP within this large sample, and apply alignment clustering to identify conserved motifs at the genomic scale. Together, these three projects highlight the power of SHAPE to identify specific RNA structures related to human disease and the value of robust experimental design and careful analysis in large-scale sequencing studies of RNA structure.Doctor of Philosoph

Role of Context in RNA Structure: Flanking Sequences Reconfigure CAG Motif Folding in Huntingtin Exon 1 Transcripts

The length of the CAG repeat region in the huntingtin messenger RNA is predictive of Huntington’s disease. Structural studies of CAG repeat-containing RNAs suggest that these sequences form simple hairpin structures; however, in the context of the full-length huntingtin mRNA, CAG repeats may form complex structures that could be targeted for therapeutic intervention. We examined the structures of transcripts spanning the first exon of the huntingtin mRNA with both healthy and disease-prone repeat lengths. In transcripts with 17 to 70 repeats, the CAG sequences base paired extensively with bases in the 5′ UTR and with a conserved region downstream of the CCG repeat region. In huntingtin transcripts with healthy numbers of repeats, the previously observed CAG hairpin was either absent or short. In contrast, in transcripts with disease-associated numbers of repeats, a CAG hairpin was present and extended from a three-helix junction. Our findings demonstrate the profound importance of sequence context in RNA folding and identify specific structural differences between healthy and disease-inducing huntingtin alleles that may be targets for therapeutic intervention

RNA SHAPE Analysis of Small RNAs and Riboswitches

We describe structural analysis of RNAs by SHAPE chemical probing. RNAs are treated with 1-methyl-7-nitroisatoic anhydride (1M7), a reagent that detects local nucleotides flexibility, and N-methylisatoic anhydride (NMIA) and 1-methyl-6-nitroisatoic anhydride (1M6), reagents which together detect higher-order and non-canonical interactions. Chemical adducts are detected as stops during reverse transcriptase-mediated primer extension. Probing information can be used to infer conformational changes and ligand binding, and to develop highly accurate models of RNA secondary structures

Selective 2′-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) for direct, versatile and accurate RNA structure analysis

SHAPE chemistries exploit small electrophilic reagents that react with the 2′-hydroxyl group to interrogate RNA structure at single-nucleotide resolution. Mutational profiling (MaP) identifies modified residues based on the ability of reverse transcriptase to misread a SHAPE-modified nucleotide and then counting the resulting mutations by massively parallel sequencing. The SHAPE-MaP approach measures the structure of large and transcriptome-wide systems as accurately as for simple model RNAs. This protocol describes the experimental steps, implemented over three days, required to perform SHAPE probing and construct multiplexed SHAPE-MaP libraries suitable for deep sequencing. These steps include RNA folding and SHAPE structure probing, mutational profiling by reverse transcription, library construction, and sequencing. Automated processing of MaP sequencing data is accomplished using two software packages. ShapeMapper converts raw sequencing files into mutational profiles, creates SHAPE reactivity plots, and provides useful troubleshooting information, often within an hour. SuperFold uses these data to model RNA secondary structures, identify regions with well-defined structures, and visualize probable and alternative helices, often in under a day. We illustrate these algorithms with the E. coli thiamine pyrophosphate riboswitch, E. coli 16S rRNA, and HIV-1 genomic RNAs. SHAPE-MaP can be used to make nucleotide-resolution biophysical measurements of individual RNA motifs, rare components of complex RNA ensembles, and entire transcriptomes. The straightforward MaP strategy greatly expands the number, length, and complexity of analyzable RNA structures

RNA motif discovery by SHAPE and mutational profiling (SHAPE-MaP)

Many biological processes are RNA-mediated, but higher-order structures for most RNAs are unknown, making it difficult to understand how RNA structure governs function. Here we describe SHAPE mutational profiling (SHAPE-MaP) that makes possible de novo and large-scale identification of RNA functional motifs. Sites of 2’-hydroxyl acylation by SHAPE are encoded as non-complementary nucleotides during cDNA synthesis, as measured by massively parallel sequencing. SHAPE-MaP-guided modeling identified greater than 90% of accepted base pairs in complex RNAs of known structure and was used to define a second-generation model for the HIV-1 RNA genome. The HIV-1 model contains all known structured motifs and previously unknown elements, including experimentally validated pseudoknots. SHAPE-MaP yields accurate and high-resolution secondary structure models, enables analysis of low abundance RNAs, disentangles sequence polymorphisms in single experiments, and will ultimately democratize RNA structure analysis

Presynaptic Inhibition of Gamma-Aminobutyric Acid Release in the Bed Nucleus of the Stria Terminalis by Kappa Opioid Receptor Signaling

The kappa opioid receptor (KOR) and its endogenous agonist, the neuropeptide dynorphin, are a critical component of the central stress system. Both dynorphin and KOR are expressed in the bed nucleus of the stria terminalis (BNST), a brain region associated with anxiety and stress. This suggests that KOR activation in this region may play a role in the regulation of emotional behaviors. To date, however, there has been no investigation of the ability of KOR to modulate synaptic transmission in the BNST

Direct detection of RNA modifications and structure using single-molecule nanopore sequencing

Summary: Modifications are present on many classes of RNA, including tRNA, rRNA, and mRNA. These modifications modulate diverse biological processes such as genetic recoding and mRNA export and folding. In addition, modifications can be introduced to RNA molecules using chemical probing strategies that reveal RNA structure and dynamics. Many methods exist to detect RNA modifications by short-read sequencing; however, limitations on read length inherent to short-read-based methods dissociate modifications from their native context, preventing single-molecule modification analysis. Here, we demonstrate direct RNA nanopore sequencing to detect endogenous and exogenous RNA modifications on long RNAs at the single-molecule level. We detect endogenous 2′-O-methyl and base modifications across E. coli and S. cerevisiae ribosomal RNAs as shifts in current signal and dwell times distally through interactions with the helicase motor protein. We further use the 2′-hydroxyl reactive SHAPE reagent acetylimidazole to probe RNA structure at the single-molecule level with readout by direct nanopore sequencing

Pervasive Regulatory Functions of mRNA Structure Revealed by High-Resolution SHAPE Probing

mRNAs can fold into complex structures that regulate gene expression. Resolving such structures de novo has remained challenging and has limited our understanding of the prevalence and functions of mRNA structure. We use SHAPE-MaP experiments in living E. coli cells to derive quantitative, nucleotide-resolution structure models for 194 endogenous transcripts encompassing approximately 400 genes. Individual mRNAs have exceptionally diverse architectures, and most contain well-defined structures. Active translation destabilizes mRNA structure in cells. Nevertheless, mRNA structure remains similar between in-cell and cell-free environments, indicating broad potential for structure-mediated gene regulation. We find that the translation efficiency of endogenous genes is regulated by unfolding kinetics of structures overlapping the ribosome binding site. We discover conserved structured elements in 35% of UTRs, several of which we validate as novel protein binding motifs. RNA structure regulates every gene studied here in a meaningful way, implying that most functional structures remain to be discovered