Pairwise visual comparison of small RNA secondary structures with base pair probabilities

AbstractBackgroundPredicted RNA secondary structures are typically visualized using dot-plots for base pair binding probabilities and planar graphs for unique structures, such as the minimum free energy structure. These are however difficult to analyze simultaneously.ResultsThis work introduces a compact unified view of the most stable conformation of an RNA secondary structure and its base pair probabilities, which is called the Circular Secondary Structure Base Pairs Probabilities Plot (CS2BP2-Plot). Along with our design we provide access to a web server implementation of our solution that facilitates pairwise comparison of short RNA (and DNA) sequences up to 200 base pairs. The web server first calculates the minimum free energy secondary structure and the base pair probabilities for up to 10 RNA or DNA sequences using RNAfold and then provides a two panel comparative view that includes CS2BP2-Plots along with the traditional graph, planar and circular diagrams obtained with VARNA. The CS2BP2-Plots include highlighting of the nucleotide differences between two selected sequences using ClustalW local alignments. We also provide descriptive statistics, dot-bracket secondary structure representations and ClustalW local alignments for compared sequences.ConclusionsUsing circular diagrams and colour and weight-coded arcs, we demonstrate how a single image can replace the state-of-the-art dual representations (dot-plots and minimum free energy structures) for base-pair probabilities of RNA secondary structures while allowing efficient exploration and comparison of different RNA conformations via a web server front end. With that, we provide the community, especially the biologically oriented, with an intuitive tool for ncRNA visualization. Web-server:https://nrcmonsrv01.nrc.ca/cs2bp2plot

[1]  R. Micura,et al.  Bistable secondary structures of small RNAs and their structural probing by comparative imino proton NMR spectroscopy. , 2003, Journal of molecular biology.

[2]  Alexander Churkin,et al.  RNA dot plots: an image representation for RNA secondary structure analysis and manipulations , 2013, Wiley interdisciplinary reviews. RNA.

[3]  David H Mathews,et al.  RNA structure prediction: an overview of methods. , 2012, Methods in molecular biology.

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

[5]  Steven J. M. Jones,et al.  Circos: an information aesthetic for comparative genomics. , 2009, Genome research.

[6]  Tessa G Montague,et al.  Internal guide RNA interactions interfere with Cas9-mediated cleavage , 2016, Nature Communications.

[7]  Yu Zhang,et al.  RNA-TVcurve: a Web server for RNA secondary structure comparison based on a multi-scale similarity of its triple vector curve representation , 2017, BMC Bioinformatics.

[8]  D. Bartel,et al.  One sequence, two ribozymes: implications for the emergence of new ribozyme folds. , 2000, Science.

[9]  Katja Nowick,et al.  Temporal ordering of substitutions in RNA evolution: Uncovering the structural evolution of the Human Accelerated Region 1. , 2018, Journal of theoretical biology.

[10]  P. Stadler,et al.  Design of multistable RNA molecules. , 2001, RNA.

[11]  Eivind Tøstesen,et al.  RNA Folding Transitions and Cooperativity , 2001 .

[12]  Harald Schwalbe,et al.  NMR Studies of HAR1 RNA Secondary Structures Reveal Conformational Dynamics in the Human RNA , 2012, Chembiochem : a European journal of chemical biology.

[13]  J. Doudna,et al.  Insights into RNA structure and function from genome-wide studies , 2014, Nature Reviews Genetics.

[14]  Daniel Lai,et al.  R-chie: a web server and R package for visualizing RNA secondary structures , 2012, Nucleic acids research.

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

[16]  Jerrold R. Griggs,et al.  Algorithms for Loop Matchings , 1978 .

[17]  Robert Giegerich,et al.  Local similarity in RNA secondary structures , 2003, Computational Systems Bioinformatics. CSB2003. Proceedings of the 2003 IEEE Bioinformatics Conference. CSB2003.

[18]  Howard Y. Chang,et al.  Detecting riboSNitches with RNA folding algorithms: a genome-wide benchmark , 2015, Nucleic acids research.

[19]  D. Haussler,et al.  An RNA gene expressed during cortical development evolved rapidly in humans , 2006, Nature.

[20]  Eugenio Mattei,et al.  Web-Beagle: a web server for the alignment of RNA secondary structures , 2015, Nucleic Acids Res..

[21]  Eric Westhof,et al.  Distinctive structures between chimpanzee and human in a brain noncoding RNA. , 2008, RNA.

[22]  Anne Condon,et al.  RNAsoft: a suite of RNA secondary structure prediction and design software tools , 2003, Nucleic Acids Res..

[23]  Sebastian Will,et al.  Temperature-responsive in vitro RNA structurome of Yersinia pseudotuberculosis , 2016, Proceedings of the National Academy of Sciences.

[24]  A. Visel,et al.  Response to Comment on "Human-Specific Gain of Function in a Developmental Enhancer" , 2009, Science.

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

[26]  David Hoksza,et al.  TRAVeLer: a tool for template-based RNA secondary structure visualization , 2017, BMC Bioinformatics.

[27]  Michael Zuker,et al.  Mfold web server for nucleic acid folding and hybridization prediction , 2003, Nucleic Acids Res..

[28]  William Jannen,et al.  Visualizing RNA base-pairing probabilities with RNAbow diagrams. , 2013, RNA.

[29]  Eric Westhof,et al.  The Dynamic Landscapes of RNA Architecture , 2009, Cell.

[30]  Janet Kelso,et al.  Analysis of Human Accelerated DNA Regions Using Archaic Hominin Genomes , 2012, PloS one.

[31]  Alain Laederach,et al.  Disease-Associated Mutations That Alter the RNA Structural Ensemble , 2010, PLoS genetics.