Structure-seq2: sensitive and accurate genome-wide profiling of RNA structure in vivo

Abstract RNA serves many functions in biology such as splicing, temperature sensing, and innate immunity. These functions are often determined by the structure of RNA. There is thus a pressing need to understand RNA structure and how it changes during diverse biological processes both in vivo and genome-wide. Here, we present Structure-seq2, which provides nucleotide-resolution RNA structural information in vivo and genome-wide. This optimized version of our original Structure-seq method increases sensitivity by at least 4-fold and improves data quality by minimizing formation of a deleterious by-product, reducing ligation bias, and improving read coverage. We also present a variation of Structure-seq2 in which a biotinylated nucleotide is incorporated during reverse transcription, which greatly facilitates the protocol by eliminating two PAGE purification steps. We benchmark Structure-seq2 on both mRNA and rRNA structure in rice (Oryza sativa). We demonstrate that Structure-seq2 can lead to new biological insights. Our Structure-seq2 datasets uncover hidden breaks in chloroplast rRNA and identify a previously unreported N1-methyladenosine (m1A) in a nuclear-encoded Oryza sativa rRNA. Overall, Structure-seq2 is a rapid, sensitive, and unbiased method to probe RNA in vivo and genome-wide that facilitates new insights into RNA biology.

[1]  David H. Mathews,et al.  RNAstructure: software for RNA secondary structure prediction and analysis , 2010, BMC Bioinformatics.

[2]  Cameron S. Osborne,et al.  Large Scale Loss of Data in Low-Diversity Illumina Sequencing Libraries Can Be Recovered by Deferred Cluster Calling , 2011, PloS one.

[3]  C. Leaver,et al.  Molecular integrity of chloroplast ribosomal ribonucleic acid. , 1971, The Biochemical journal.

[4]  Yin Tang,et al.  The RNA structurome: transcriptome-wide structure probing with next-generation sequencing. , 2015, Trends in biochemical sciences.

[5]  Y. Zhang,et al.  In vivo genome-wide profiling of RNA secondary structure reveals novel regulatory features , 2013, Nature.

[6]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[7]  Catherine Tran,et al.  Progress and challenges for chemical probing of RNA structure inside living cells. , 2015, Nature chemical biology.

[8]  S. Oliviero,et al.  Genome-wide profiling of mouse RNA secondary structures reveals key features of the mammalian transcriptome , 2014, Genome Biology.

[9]  Stuart Aitken,et al.  Snapshots of pre-rRNA structural flexibility reveal eukaryotic 40S assembly dynamics at nucleotide resolution , 2014, Nucleic acids research.

[10]  Lukasz Jan Kielpinski,et al.  Massive parallel-sequencing-based hydroxyl radical probing of RNA accessibility , 2014, Nucleic acids research.

[11]  Rhiju Das,et al.  Massively parallel RNA chemical mapping with a reduced bias MAP-seq protocol. , 2013, Methods in molecular biology.

[12]  Zhao Su,et al.  Genome-Wide Analysis of RNA Secondary Structure. , 2016, Annual review of genetics.

[13]  Sarah M Assmann,et al.  Genome-wide profiling of in vivo RNA structure at single-nucleotide resolution using structure-seq , 2015, Nature Protocols.

[14]  Daniel Boehringer,et al.  The complete structure of the chloroplast 70S ribosome in complex with translation factor pY , 2016, The EMBO journal.

[15]  R. Bock,et al.  The Conserved Endoribonuclease YbeY Is Required for Chloroplast Ribosomal RNA Processing in Arabidopsis1 , 2015, Plant Physiology.

[16]  A. Krogh,et al.  SHAPE Selection (SHAPES) enrich for RNA structure signal in SHAPE sequencing-based probing data. , 2015, RNA.

[17]  Kyle E. Watters,et al.  SHAPE-Seq 2.0: systematic optimization and extension of high-throughput chemical probing of RNA secondary structure with next generation sequencing , 2014, Nucleic acids research.

[18]  Yiliang Ding,et al.  A hybridization-based approach for quantitative and low-bias single-stranded DNA ligation. , 2013, Analytical biochemistry.

[19]  Nan Yu,et al.  The Comparative RNA Web (CRW) Site: an online database of comparative sequence and structure information for ribosomal, intron, and other RNAs , 2002, BMC Bioinformatics.

[20]  V. Ambros,et al.  An efficient and sensitive method for preparing cDNA libraries from scarce biological samples , 2014, Nucleic acids research.

[21]  D. Mathews,et al.  Bridging the gap between in vitro and in vivo RNA folding , 2016, Quarterly Reviews of Biophysics.

[22]  Andreas Hildebrandt,et al.  The reverse transcription signature of N-1-methyladenosine in RNA-Seq is sequence dependent , 2015, Nucleic acids research.

[23]  Yves Van de Peer,et al.  PLAZA 3.0: an access point for plant comparative genomics , 2014, Nucleic Acids Res..

[24]  Anton Nekrutenko,et al.  StructureFold: genome-wide RNA secondary structure mapping and reconstruction in vivo , 2015, Bioinform..

[25]  Maurille J. Fournier,et al.  The 3D rRNA modification maps database: with interactive tools for ribosome analysis , 2007, Nucleic Acids Res..

[26]  Kyle E. Watters,et al.  RNA systems biology: uniting functional discoveries and structural tools to understand global roles of RNAs. , 2016, Current opinion in biotechnology.

[27]  Andreas Hildebrandt,et al.  CoverageAnalyzer (CAn): A Tool for Inspection of Modification Signatures in RNA Sequencing Profiles , 2016, Biomolecules.

[28]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.