RNA Framework: an all-in-one toolkit for the analysis of RNA structures and post-transcriptional modifications

Abstract RNA is emerging as a key regulator of a plethora of biological processes. While its study has remained elusive for decades, the recent advent of high-throughput sequencing technologies provided the unique opportunity to develop novel techniques for the study of RNA structure and post-transcriptional modifications. Nonetheless, most of the required downstream bioinformatics analyses steps are not easily reproducible, thus making the application of these techniques a prerogative of few laboratories. Here we introduce RNA Framework, an all-in-one toolkit for the analysis of most NGS-based RNA structure probing and post-transcriptional modification mapping experiments. To prove the extreme versatility of RNA Framework, we applied it to both an in-house generated DMS-MaPseq dataset, and to a series of literature available experiments. Notably, when starting from publicly available datasets, our software easily allows replicating authors' findings. Collectively, RNA Framework provides the most complete and versatile toolkit to date for a rapid and streamlined analysis of the RNA epistructurome. RNA Framework is available for download at: http://www.rnaframework.com.

[1]  Yi Xing,et al.  m(6)A RNA modification controls cell fate transition in mammalian embryonic stem cells. , 2014, Cell stem cell.

[2]  M. Kupiec,et al.  Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq , 2012, Nature.

[3]  Timothy L. Bailey,et al.  Gene expression Advance Access publication May 4, 2011 DREME: motif discovery in transcription factor ChIP-seq data , 2011 .

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

[5]  Howard Y. Chang,et al.  Genome-wide measurement of RNA secondary structure in yeast , 2010, Nature.

[6]  Kevin Y. Yip,et al.  Improved prediction of RNA secondary structure by integrating the free energy model with restraints derived from experimental probing data , 2015, Nucleic acids research.

[7]  Terrence S. Furey,et al.  The UCSC Genome Browser Database , 2003, Nucleic Acids Res..

[8]  D. Mathews,et al.  Accurate SHAPE-directed RNA structure determination , 2009, Proceedings of the National Academy of Sciences.

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

[10]  Steven Busan,et al.  RNA motif discovery by SHAPE and mutational profiling (SHAPE-MaP) , 2014, Nature Methods.

[11]  Howard Y. Chang,et al.  Structural imprints in vivo decode RNA regulatory mechanisms , 2015, Nature.

[12]  Yin Tang,et al.  StructureFold2: Bringing chemical probing data into the computational fold of RNA structural analysis. , 2018, Methods.

[13]  J. Weissman,et al.  DMS-MaPseq for genome-wide or targeted RNA structure probing in vivo , 2016, Nature Methods.

[14]  Yufei Huang,et al.  DRME: Count-based differential RNA methylation analysis at small sample size scenario. , 2016, Analytical biochemistry.

[15]  Yufei Huang,et al.  HEPeak: an HMM-based exome peak-finding package for RNA epigenome sequencing data , 2015, BMC Genomics.

[16]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[17]  Maxwell R. Mumbach,et al.  Transcriptome-wide Mapping Reveals Widespread Dynamic-Regulated Pseudouridylation of ncRNA and mRNA , 2014, Cell.

[18]  D. Haussler,et al.  FragSeq: transcriptome-wide RNA structure probing using high-throughput sequencing , 2010, Nature Methods.

[19]  O. Elemento,et al.  Comprehensive Analysis of mRNA Methylation Reveals Enrichment in 3′ UTRs and near Stop Codons , 2012, Cell.

[20]  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.

[21]  Schraga Schwartz,et al.  Transcriptome-Wide Mapping of 5-methylcytidine RNA Modifications in Bacteria, Archaea, and Yeast Reveals m5C within Archaeal mRNAs , 2013, PLoS genetics.

[22]  Jan Gorodkin,et al.  Profiling of ribose methylations in RNA by high-throughput sequencing. , 2014, Angewandte Chemie.

[23]  W. Gilbert,et al.  Pseudouridine profiling reveals regulated mRNA pseudouridylation in yeast and human cells , 2014, Nature.

[24]  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.

[25]  Cole Trapnell,et al.  Ultrafast and memory-efficient alignment of short DNA sequences to the human genome , 2009, Genome Biology.

[26]  Yann Ponty,et al.  VARNA: Interactive drawing and editing of the RNA secondary structure , 2009, Bioinform..

[27]  Yuri Motorin,et al.  Illumina-based RiboMethSeq approach for mapping of 2′-O-Me residues in RNA , 2016, Nucleic acids research.

[28]  Cole Trapnell,et al.  Multiplexed RNA structure characterization with selective 2′-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq) , 2011, Proceedings of the National Academy of Sciences.

[29]  Xiaodong Cui,et al.  MeTDiff: A Novel Differential RNA Methylation Analysis for MeRIP-Seq Data , 2018, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[30]  Peter Clote,et al.  Integrating Chemical Footprinting Data into RNA Secondary Structure Prediction , 2012, PloS one.

[31]  Cole Trapnell,et al.  Modeling and automation of sequencing-based characterization of RNA structure , 2011, Proceedings of the National Academy of Sciences.

[32]  D. Mathews,et al.  Accurate SHAPE-directed RNA secondary structure modeling, including pseudoknots , 2013, Proceedings of the National Academy of Sciences.

[33]  Chengqi Yi,et al.  Transcriptome-wide mapping reveals reversible and dynamic N(1)-methyladenosine methylome. , 2016, Nature chemical biology.

[34]  J. Woolford,et al.  Mod-seq: high-throughput sequencing for chemical probing of RNA structure , 2014, RNA.

[35]  Yufei Huang,et al.  A hierarchical model for clustering m6A methylation peaks in MeRIP-seq data , 2016, BMC Genomics.

[36]  Aaron R. Quinlan,et al.  Bioinformatics Applications Note Genome Analysis Bedtools: a Flexible Suite of Utilities for Comparing Genomic Features , 2022 .

[37]  Robert D. Finn,et al.  Rfam 13.0: shifting to a genome-centric resource for non-coding RNA families , 2017, Nucleic Acids Res..

[38]  Manolis Kellis,et al.  Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo , 2013, Nature.

[39]  Kevin M Weeks,et al.  Selective 2′-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) for direct, versatile and accurate RNA structure analysis , 2015, Nature Protocols.

[40]  Francesco Neri,et al.  High-throughput single-base resolution mapping of RNA 2΄-O-methylated residues , 2016, Nucleic acids research.

[41]  Fei Wang,et al.  Transcriptome-wide distribution and function of RNA hydroxymethylcytosine , 2016, Science.

[42]  Jun Yu,et al.  MeRIP-PF: An Easy-to-use Pipeline for High-resolution Peak-finding in MeRIP-Seq Data , 2013, Genom. Proteom. Bioinform..

[43]  Nikolay V. Dokholyan,et al.  Single-molecule correlated chemical probing of RNA , 2014, Proceedings of the National Academy of Sciences.

[44]  S. Oliviero,et al.  In vivo probing of nascent RNA structures reveals principles of cotranscriptional folding , 2017, Nucleic acids research.

[45]  David Haussler,et al.  The UCSC Genome Browser database: 2018 update , 2017, Nucleic Acids Res..

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

[47]  K. Weeks,et al.  Accurate detection of chemical modifications in RNA by mutational profiling (MaP) with ShapeMapper 2 , 2018, RNA.

[48]  Helga Thorvaldsdóttir,et al.  Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration , 2012, Briefings Bioinform..

[49]  T. Preiss,et al.  Widespread occurrence of 5-methylcytosine in human coding and non-coding RNA , 2012, Nucleic acids research.

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

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

[52]  Danny Incarnato,et al.  The RNA Epistructurome: Uncovering RNA Function by Studying Structure and Post-Transcriptional Modifications. , 2017, Trends in biotechnology.

[53]  Peter F. Stadler,et al.  ViennaRNA Package 2.0 , 2011, Algorithms for Molecular Biology.

[54]  Isabelle Behm-Ansmant,et al.  Use of Specific Chemical Reagents for Detection of Modified Nucleotides in RNA , 2011, Journal of nucleic acids.

[55]  Gideon Rechavi,et al.  The dynamic N1-methyladenosine methylome in eukaryotic messenger RNA , 2016, Nature.

[56]  Howard Y. Chang,et al.  Comparison of SHAPE reagents for mapping RNA structures inside living cells , 2017, RNA.