Circular RNA profile in gliomas revealed by identification tool UROBORUS

Recent evidence suggests that many endogenous circular RNAs (circRNAs) may play roles in biological processes. However, the expression patterns and functions of circRNAs in human diseases are not well understood. Computationally identifying circRNAs from total RNA-seq data is a primary step in studying their expression pattern and biological roles. In this work, we have developed a computational pipeline named UROBORUS to detect circRNAs in total RNA-seq data. By applying UROBORUS to RNA-seq data from 46 gliomas and normal brain samples, we detected thousands of circRNAs supported by at least two read counts, followed by successful experimental validation on 24 circRNAs from the randomly selected 27 circRNAs. UROBORUS is an efficient tool that can detect circRNAs with low expression levels in total RNA-seq without RNase R treatment. The circRNAs expression profiling revealed more than 476 circular RNAs differentially expressed in control brain tissues and gliomas. Together with parental gene expression, we found that circRNA and its parental gene have diversified expression patterns in gliomas and control brain tissues. This study establishes an efficient and sensitive approach for predicting circRNAs using total RNA-seq data. The UROBORUS pipeline can be accessed freely for non-commercial purposes at http://uroborus.openbioinformatics.org/.

[1]  D. Riesner,et al.  Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures. , 1976, Proceedings of the National Academy of Sciences of the United States of America.

[2]  C. Cocquerelle,et al.  Splicing with inverted order of exons occurs proximal to large introns. , 1992, The EMBO journal.

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

[4]  C. Cocquerelle,et al.  Mis‐splicing yields circular RNA molecules , 1993, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[5]  M. Garcia-Blanco,et al.  Exon circularization in mammalian nuclear extracts. , 1996, RNA.

[6]  H. Domdey,et al.  Inverse splicing of a discontinuous pre-mRNA intron generates a circular exon in a HeLa cell nuclear extract. , 1996, Nucleic acids research.

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

[8]  Lior Pachter,et al.  Sequence Analysis , 2020, Definitions.

[9]  Cole Trapnell,et al.  Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. , 2010, Nature biotechnology.

[10]  William R. Jeck,et al.  Expression of Linear and Novel Circular Forms of an INK4/ARF-Associated Non-Coding RNA Correlates with Atherosclerosis Risk , 2010, PLoS genetics.

[11]  S. Salzberg,et al.  TopHat-Fusion: an algorithm for discovery of novel fusion transcripts , 2011, Genome Biology.

[12]  Schraga Schwartz,et al.  Transcriptome-wide discovery of circular RNAs in Archaea , 2011, Nucleic acids research.

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

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

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

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

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

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

[19]  A. Yee,et al.  Versican V2 isoform enhances angiogenesis by regulating endothelial cell activities and fibronectin expression , 2013, FEBS letters.

[20]  Shanshan Zhu,et al.  Circular intronic long noncoding RNAs. , 2013, Molecular cell.

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

[22]  R. Parker,et al.  Circular RNAs: diversity of form and function , 2014, RNA.

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

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

[25]  Susan M. Chang,et al.  Where are we now? And where are we going? A report from the Accelerate Brain Cancer Cure (ABC2) low-grade glioma research workshop. , 2014, Neuro-Oncology.

[26]  Petar Glažar,et al.  circBase: a database for circular RNAs , 2014, RNA.

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

[28]  R. Zeillinger,et al.  Correlation of circular RNA abundance with proliferation – exemplified with colorectal and ovarian cancer, idiopathic lung fibrosis, and normal human tissues , 2015, Scientific Reports.

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

[30]  G. von Heijne,et al.  Tissue-based map of the human proteome , 2015, Science.

[31]  Andreas W. Schreiber,et al.  The RNA Binding Protein Quaking Regulates Formation of circRNAs , 2015, Cell.

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

[33]  L. Christophorou.  Science , 2018, Emerging Dynamics: Science, Energy, Society and Values.