Circular RNAs (circRNAs) derived from back-spliced exons have been widely identified as being co-expressed with their linear counterparts. A single gene locus can produce multiple circRNAs through alternative back-splice site selection and/or alternative splice site selection; however, a detailed map of alternative back-splicing/splicing in circRNAs is lacking. Here, with the upgraded CIRCexplorer2 pipeline, we systematically annotated different types of alternative back-splicing and alternative splicing events in circRNAs from various cell lines. Compared with their linear cognate RNAs, circRNAs exhibited distinct patterns of alternative back-splicing and alternative splicing. Alternative back-splice site selection was correlated with the competition of putative RNA pairs across introns that bracket alternative back-splice sites. In addition, all four basic types of alternative splicing that have been identified in the (linear) mRNA process were found within circRNAs, and many exons were predominantly spliced in circRNAs. Unexpectedly, thousands of previously unannotated exons were detected in circRNAs from the examined cell lines. Although these novel exons had similar splice site strength, they were much less conserved than known exons in sequences. Finally, both alternative back-splicing and circRNA-predominant alternative splicing were highly diverse among the examined cell lines. All of the identified alternative back-splicing and alternative splicing in circRNAs are available in the CIRCpedia database (http://www.picb.ac.cn/rnomics/ circpedia). Collectively, the annotation of alternative back-splicing and alternative splicing in circRNAs provides a valuable resource for depicting the complexity of circRNA biogenesis and for studying the potential functions of circRNAs in different cells.

SUMMARY Many eukaryotic genes generate linear mRNAs and circular RNAs, but it is largely unknown how the ratio of linear to circular RNA is controlled or modulated. Using RNAi screening in Drosophila cells, we identify many core spliceosome and transcription termination factors that control the RNA outputs of reporter and endogenous genes. When spliceosome components were depleted or inhibited pharmacologically, the steady-state levels of circular RNAs increased while expression of their associated linear mRNAs concomitantly decreased. Upon inhibiting RNA polymerase II termination via depletion of the cleavage/polyadenylation machinery, circular RNA levels were similarly increased. This is because readthrough transcripts now extend into downstream genes and are subjected to backsplicing. In total, these results demonstrate that inhibition or slowing of canonical pre-mRNA processing events shifts the steady-state output of protein-coding genes towards circular RNAs. This is in part because nascent RNAs become directed into alternative pathways that lead to circular RNA production.

We identify a type of polycistronic transcript-derived long noncoding RNAs (lncRNAs) that are 5' small nucleolar RNA (snoRNA) capped and 3' polyadenylated (SPAs). SPA processing is associated with nascent mRNA 3' processing and kinetic competition between XRN2 trimming and Pol II elongation. Following cleavage/polyadenylation of its upstream gene, the downstream uncapped pre-SPA is trimmed by XRN2 until this exonuclease reaches the co-transcriptionally assembled snoRNP. This snoRNP complex prevents further degradation, generates a snoRNA 5' end, and allows continuous Pol II elongation. The imprinted 15q11-q13 encodes two SPAs that are deleted in Prader-Willi syndrome (PWS) patients. These lncRNAs form a nuclear accumulation that is enriched in RNA binding proteins (RBPs) including TDP43, RBFOX2, and hnRNP M. Generation of a human PWS cellular model by depleting these lncRNAs results in altered patterns of RBPs binding and alternative splicing. Together, these results expand the diversity of lncRNAs and provide additional insights into PWS pathogenesis.

No abstract