SpliceV: analysis and publication quality printing of linear and circular RNA splicing, expression and regulation

BackgroundIn eukaryotes, most genes code for multiple transcript isoforms that are generated through the complex and tightly regulated process of RNA splicing. Despite arising from identical precursor transcripts, alternatively spliced RNAs can have dramatically different functions. Transcriptome complexity is elevated further by the production of circular RNAs (circRNAs), another class of mature RNA that results from the splicing of a downstream splice donor to an upstream splice acceptor. While there has been a rapid expansion of circRNA catalogs in the last few years through the utilization of next generation sequencing approaches, our understanding of the mechanisms and regulation of circular RNA biogenesis, the impact that circRNA generation has on parental transcript processing, and the functions carried out by circular RNAs remains limited.ResultsHere, we present a visualization and analysis tool, SpliceV, that rapidly plots all relevant forward- and back-splice data, with exon and single nucleotide level coverage information from RNA-seq experiments in a publication quality format. SpliceV also integrates analysis features that assist investigations into splicing regulation and transcript functions through the display of predicted RNA binding protein sites and the configuration of repetitive elements along the primary transcript.ConclusionsSpliceV is an easy-to-use splicing visualization tool, compatible with both Python 2.7 and 3+, and distributed under the GNU Public License. The source code is freely available for download at https://github.com/flemingtonlab/SpliceV and can be installed from PyPI using pip.

[1]  N. Rajewsky,et al.  Translation of CircRNAs , 2017, Molecular cell.

[2]  Christoph Dieterich,et al.  Analysis of intron sequences reveals hallmarks of circular RNA biogenesis in animals. , 2015, Cell reports.

[3]  Steven J. M. Jones,et al.  Comprehensive molecular characterization of gastric adenocarcinoma , 2014, Nature.

[4]  A. Schedl,et al.  Two Splice Variants of the Wilms' Tumor 1 Gene Have Distinct Functions during Sex Determination and Nephron Formation , 2001, Cell.

[5]  J. Kjems,et al.  Natural RNA circles function as efficient microRNA sponges , 2013, Nature.

[6]  J. Valcárcel,et al.  Distinct binding specificities and functions of higher eukaryotic polypyrimidine tract-binding proteins. , 1995, Science.

[7]  A. Weiner,et al.  A compensatory base change in human U2 snRNA can suppress a branch site mutation. , 1989, Genes & development.

[8]  D. Solnick Alternative splicing caused by RNA secondary structure , 1985, Cell.

[9]  Michael K. Slevin,et al.  Circular RNAs are abundant, conserved, and associated with ALU repeats. , 2013, RNA.

[10]  Roy Parker,et al.  Circular RNAs Co-Precipitate with Extracellular Vesicles: A Possible Mechanism for circRNA Clearance , 2016, PloS one.

[11]  Brendan J. Frey,et al.  A compendium of RNA-binding motifs for decoding gene regulation , 2013, Nature.

[12]  Charles Gawad,et al.  Circular RNAs Are the Predominant Transcript Isoform from Hundreds of Human Genes in Diverse Cell Types , 2012, PloS one.

[13]  Michael G. Rosenfeld,et al.  Alternative RNA processing in calcitonin gene expression generates mRNAs encoding different polypeptide products , 1982, Nature.

[14]  Sebastian D. Mackowiak,et al.  Circular RNAs are a large class of animal RNAs with regulatory potency , 2013, Nature.

[15]  B. Frey,et al.  Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing , 2008, Nature Genetics.

[16]  C. Mayr,et al.  Widespread Shortening of 3′UTRs by Alternative Cleavage and Polyadenylation Activates Oncogenes in Cancer Cells , 2009, Cell.

[17]  N. Rajewsky,et al.  Circ-ZNF609 Is a Circular RNA that Can Be Translated and Functions in Myogenesis , 2017, Molecular cell.

[18]  T. Curran,et al.  Binding of the Wilms' tumor locus zinc finger protein to the EGR-1 consensus sequence. , 1990, Science.

[19]  Li Yang,et al.  The Biogenesis of Nascent Circular RNAs. , 2016, Cell reports.

[20]  F. Zhao,et al.  CIRI: an efficient and unbiased algorithm for de novo circular RNA identification , 2015, Genome Biology.

[21]  S. Srivastav,et al.  Global Bidirectional Transcription of the Epstein-Barr Virus Genome during Reactivation , 2013, Journal of Virology.

[22]  Gaël Varoquaux,et al.  The NumPy Array: A Structure for Efficient Numerical Computation , 2011, Computing in Science & Engineering.

[23]  Tom Maniatis,et al.  A role for exon sequences and splice-site proximity in splice-site selection , 1986, Cell.

[24]  C. Will,et al.  The Spliceosome: Design Principles of a Dynamic RNP Machine , 2009, Cell.

[25]  Eric T. Wang,et al.  Alternative Isoform Regulation in Human Tissue Transcriptomes , 2008, Nature.

[26]  L. Swanson,et al.  Production of a novel neuropeptide encoded by the calcitonin gene via tissue-specific RNA processing , 1983, Nature.

[27]  M. Waterfield,et al.  Glial growth factors are alternatively spliced erbB2 ligands expressed in the nervous system , 1993, Nature.

[28]  Alain Bergeron,et al.  Widespread and Functional RNA Circularization in Localized Prostate Cancer , 2019, Cell.

[29]  Adrian R. Krainer,et al.  The essential pre-mRNA splicing factor SF2 influences 5′ splice site selection by activating proximal sites , 1990, Cell.

[30]  A. Levine,et al.  The p53 tumour suppressor gene , 1991, Nature.

[31]  Allison P. Heath,et al.  Toward a Shared Vision for Cancer Genomic Data. , 2016, The New England journal of medicine.

[32]  Tim Schneider,et al.  Exon circularization requires canonical splice signals. , 2015, Cell reports.

[33]  P. Brown,et al.  Circular RNA Is Expressed across the Eukaryotic Tree of Life , 2014, PloS one.

[34]  Jinhua Wang,et al.  ESEfinder: a web resource to identify exonic splicing enhancers , 2003, Nucleic Acids Res..

[35]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[36]  Ling-Ling Chen,et al.  Complementary Sequence-Mediated Exon Circularization , 2014, Cell.

[37]  C. Thompson,et al.  bcl-x, a bcl-2-related gene that functions as a dominant regulator of apoptotic cell death , 1993, Cell.

[38]  Petar Glažar,et al.  Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed. , 2015, Molecular cell.

[39]  G. Shan,et al.  Exon-intron circular RNAs regulate transcription in the nucleus , 2015, Nature Structural &Molecular Biology.

[40]  Sol Shenker,et al.  Genome-wide analysis of drosophila circular RNAs reveals their structural and sequence properties and age-dependent neural accumulation. , 2014, Cell reports.

[41]  Ying Li,et al.  An intron with a constitutive transport element is retained in a Tap messenger RNA , 2006, Nature.

[42]  D. Bentley Coupling mRNA processing with transcription in time and space , 2014, Nature Reviews Genetics.

[43]  Walter N. Moss,et al.  The Epstein Barr virus circRNAome , 2018, PLoS pathogens.