Statistically based splicing detection reveals neural enrichment and tissue-specific induction of circular RNA during human fetal development

BackgroundThe pervasive expression of circular RNA is a recently discovered feature of gene expression in highly diverged eukaryotes, but the functions of most circular RNAs are still unknown. Computational methods to discover and quantify circular RNA are essential. Moreover, discovering biological contexts where circular RNAs are regulated will shed light on potential functional roles they may play.ResultsWe present a new algorithm that increases the sensitivity and specificity of circular RNA detection by discovering and quantifying circular and linear RNA splicing events at both annotated and un-annotated exon boundaries, including intergenic regions of the genome, with high statistical confidence. Unlike approaches that rely on read count and exon homology to determine confidence in prediction of circular RNA expression, our algorithm uses a statistical approach. Using our algorithm, we unveiled striking induction of general and tissue-specific circular RNAs, including in the heart and lung, during human fetal development. We discover regions of the human fetal brain, such as the frontal cortex, with marked enrichment for genes where circular RNA isoforms are dominant.ConclusionsThe vast majority of circular RNA production occurs at major spliceosome splice sites; however, we find the first examples of developmentally induced circular RNAs processed by the minor spliceosome, and an enriched propensity of minor spliceosome donors to splice into circular RNA at un-annotated, rather than annotated, exons. Together, these results suggest a potentially significant role for circular RNA in human development.

[1]  P. McCullagh,et al.  Generalized Linear Models, 2nd Edn. , 1990 .

[2]  J. H. Schuenemeyer,et al.  Generalized Linear Models (2nd ed.) , 1992 .

[3]  Peter Goodfellow,et al.  Circular transcripts of the testis-determining gene Sry in adult mouse testis , 1993, Cell.

[4]  J. Lytton,et al.  A Circularized Sodium-Calcium Exchanger Exon 2 Transcript* , 1999, The Journal of Biological Chemistry.

[5]  Hee-Sup Shin,et al.  The Na+-Ca2+ Exchanger Is Essential for Embryonic Heart Development in Mice , 2000, Molecules and cells.

[6]  Eric R. Ziegel,et al.  Generalized Linear Models , 2002, Technometrics.

[7]  J. Steitz,et al.  The divergent U12-type spliceosome is required for pre-mRNA splicing and is essential for development in Drosophila. , 2002, Molecular cell.

[8]  Lila R Collins,et al.  Cardiomyocytes derived from human embryonic stem cells in pro-survival factors enhance function of infarcted rat hearts , 2007, Nature Biotechnology.

[9]  J. Castle,et al.  A postnatal switch of CELF and MBNL proteins reprograms alternative splicing in the developing heart , 2008, Proceedings of the National Academy of Sciences.

[10]  Michael K. Slevin,et al.  The F-box protein Cdc4/Fbxw7 is a novel regulator of neural crest development in Xenopus laevis , 2010, Neural Development.

[11]  H. Zoghbi,et al.  Inactivation of hnRNP K by Expanded Intronic AUUCU Repeat Induces Apoptosis Via Translocation of PKCδ to Mitochondria in Spinocerebellar Ataxia 10 , 2010, PLoS genetics.

[12]  Olga K Afanasiev,et al.  Endogenous Wnt/β-Catenin Signaling Is Required for Cardiac Differentiation in Human Embryonic Stem Cells , 2010, PloS one.

[13]  T. Cooper,et al.  MicroRNAs coordinate an alternative splicing network during mouse postnatal heart development. , 2010, Genes & development.

[14]  M. Fishbein,et al.  Myocardial function with reduced expression of the sodium-calcium exchanger. , 2010, Journal of cardiac failure.

[15]  James J. Cai,et al.  Widespread establishment and regulatory impact of Alu exons in human genes , 2011, Proceedings of the National Academy of Sciences.

[16]  Akinobu Matsumoto,et al.  Fbxw7-dependent Degradation of Notch Is Required for Control of “Stemness” and Neuronal-Glial Differentiation in Neural Stem Cells* , 2011, The Journal of Biological Chemistry.

[17]  Istvan Mody,et al.  The splicing regulator Rbfox1 (A2BP1) controls neuronal excitation in the mammalian brain , 2011, Nature Genetics.

[18]  Christopher B. Burge,et al.  Alternative Splicing of RNA Triplets Is Often Regulated and Accelerates Proteome Evolution , 2012, PLoS biology.

[19]  C. Murry,et al.  Targeted Genomic Integration of a Selectable Floxed Dual Fluorescence Reporter in Human Embryonic Stem Cells , 2012, PloS one.

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

[21]  Sean P. Palecek,et al.  Robust cardiomyocyte differentiation from human pluripotent stem cells via temporal modulation of canonical Wnt signaling , 2012, Proceedings of the National Academy of Sciences.

[22]  Julia Salzman,et al.  Cell-Type Specific Features of Circular RNA Expression , 2013, PLoS genetics.

[23]  L. Pachter,et al.  Streaming fragment assignment for real-time analysis of sequencing experiments , 2012, Nature Methods.

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

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

[26]  B. Hadland,et al.  Transmembrane protein 88: a Wnt regulatory protein that specifies cardiomyocyte development , 2013, Development.

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

[28]  L. Stanton,et al.  The long noncoding RNA RMST interacts with SOX2 to regulate neurogenesis. , 2013, Molecular cell.

[29]  Andrea Tanzer,et al.  A multi-split mapping algorithm for circular RNA, splicing, trans-splicing and fusion detection , 2014, Genome Biology.

[30]  L. McIntyre,et al.  Consistent Errors in First Strand cDNA Due to Random Hexamer Mispriming , 2013, PloS one.

[31]  M. Frilander,et al.  The significant other: splicing by the minor spliceosome , 2012, Wiley interdisciplinary reviews. RNA.

[32]  N. Sharpless,et al.  Detecting and characterizing circular RNAs , 2014, Nature Biotechnology.

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

[34]  Eric T. Wang,et al.  Alternative splicing regulates vesicular trafficking genes in cardiomyocytes during postnatal heart development , 2014, Nature Communications.

[35]  S. Grimmond,et al.  Minor class splicing shapes the zebrafish transcriptome during development , 2014, Proceedings of the National Academy of Sciences.

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

[37]  D. Bartel,et al.  Expanded identification and characterization of mammalian circular RNAs , 2014, Genome Biology.

[38]  R. Padgett,et al.  Biochemical defects in minor spliceosome function in the developmental disorder MOPD I , 2014, RNA.

[39]  N. Rajewsky,et al.  circRNA biogenesis competes with pre-mRNA splicing. , 2014, Molecular cell.

[40]  Somnath Datta,et al.  Statistical analysis of next generation sequencing date , 2014 .

[41]  Charles E. Murry,et al.  Human Embryonic Stem Cell-Derived Cardiomyocytes Regenerate Non-Human Primate Hearts , 2014, Nature.

[42]  Michael B. Black,et al.  IVT-seq reveals extreme bias in RNA sequencing , 2014, Genome Biology.

[43]  Dongming Liang,et al.  Short intronic repeat sequences facilitate circular RNA production , 2014, Genes & development.

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

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

[46]  Trees-Juen Chuang,et al.  Is an observed non-co-linear RNA product spliced in trans, in cis or just in vitro? , 2014, Nucleic acids research.

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

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