Assessing the impacts of various factors related to identification, conservation, biogenesis, and function on circular RNA reliability

Circular RNAs (circRNAs) are non-polyadenylated RNAs with a continuous loop structure characterized by a non-co-linear back-splice junction (BSJ). While dozens of computational tools have been developed and identified millions of circRNA candidates in diverse species, it remains a major challenge for determining circRNA reliability due to various types of false positives. Here, we systematically assess the impacts of numerous factors related to identification, conservation, biogenesis, and function on circRNA reliability by comparisons of circRNA expression from mock (total RNAs) and the corresponding co-linear/polyadenylated RNA-depleted datasets based on three different RNA treatment approaches. Eight important indicators of circRNA reliability are determined. The relative contribution to variability explained analyses further reveal that the relative importance of these factors in affecting circRNA reliability is conservation level of circRNA > full-length circular sequences > supporting BSJ read count > both BSJ donor and acceptor splice sites at the same co-linear transcript isoforms > both BSJ donor and acceptor splice sites at the annotated exon boundaries > BSJs detected by multiple tools > supporting functional features > both BSJ donor and acceptor splice sites undergoing alternative splicing. By extracting RT-independent circRNAs, circRNAs passing multiple experimental validations, and database-specific circRNAs, we showed the additive effects of these important factors in determining circRNA reliability. This study thus provides a useful guideline and an important resource for selecting high-confidence circRNAs for further investigations.