A genome-wide circular RNA transcriptome in Rat

Abstract Background Circular RNAs are a novel class of non-coding RNAs that backsplice from 5’ donor site and 3’ acceptor site to form a circular structure. A number of circRNAs have been discovered in model organisms including human, mouse, Drosophila, among other organisms. There are a few candidate-based studies on circular RNAs in rat, a well studied model organism. The availability of a recent dataset of transcriptomes encompassing 11 tissues, 4 developmental stages and 2 genders motivated us to explore the landscape of circular RNAs in the organism. Methodology In order to understand the difference among different pipelines, we have used the same bodymap RNA sequencing dataset. A number of pipelines have been published to identify the backsplice junctions for the discovery of circRNAs but studies comparing these tools have suggested that a combination of tools would be a better approach to identify high-confidence circular RNAs. We employed 5 different combinations of tools including tophat_CIRCexplorer2, segemehl_CIRCexplorer2, star_CIRCexplorer, Bowtie2_findcirc and Bowtie2_findcirc (noHisat2) to identify circular RNAs from the dataset. Results Our analysis identified a number of tissue-specific, developmental stage specific and gender specific circular RNAs. We further independently validated 16 circRNA junctions out of 24 selected candidates in 5 tissue samples. We additionally estimated the quantitative expression of 5 circRNA candidates using real-time PCR and our analysis suggests 3 candidates as tissue-enriched Conclusion This study is one of the most comprehensive studies that provides a circular RNA transcriptome as well as to understand the difference among different computational pipelines in Rat.