Fractionation of non-polyadenylated and ribosomal-free RNAs from mammalian cells

Approaches and challenges in genome-wide circular RNA identification and quantification

Numerous circular RNAs (circRNAs) produced from back-splicing of exon(s) have been recently revealed on a genome-wide scale across species. Although generally expressed at a low level, some relatively abundant circRNAs can play regulatory roles in various biological processes, prompting continuous profiling of circRNA in broader conditions. Over the past decade, distinct strategies have been applied in both transcriptome enrichment and bioinformatic tools for detecting and quantifying circRNAs. Understanding the scope and limitations of these strategies is crucial for the subsequent annotation and characterization of circRNAs, especially those with functional potential. Here, we provide an overview of different transcriptome enrichment, deep sequencing and computational approaches for genome-wide circRNA identification, and discuss strategies for accurate quantification and characterization of circRNA.

Identification of circular RNAs of Cannabis sativa L. potentially involved in the biosynthesis of cannabinoids

We identified circRNAs in the Cannabis sativa L. genome and examined their association with 28 cannabinoids in three tissues of C. sativa. Nine circRNAs are potentially involved in the biosynthesis of six cannabinoids. Cannabis sativa L. has been widely used in the production of medicine, textiles, and food for over 2500 years. The main bioactive compounds in C. sativa are cannabinoids, which have multiple important pharmacological actions. Circular RNAs (circRNAs) play essential roles in growth and development, stress resistance, and the biosynthesis of secondary metabolites. However, the circRNAs in C. sativa remain unknown. In this study, to explore the role of circRNAs in cannabinoid biosynthesis, we performed RNA-Seq and metabolomics analysis on the leaves, roots, and stems of C. sativa. We identified 741 overlapping circRNAs by three tools, of which 717, 16, and 8 circRNAs were derived from exonic, intronic, and intergenic, respectively. Functional enrichment analysis indicated that the parental genes (PGs) of circRNAs were enriched in many processes related to biological stress responses. We found that most of the circRNAs showed tissue-specific expression and 65 circRNAs were significantly correlated with their PGs (P < 0.05, |r|≥ 0.5). We also determined 28 cannabinoids by High-performance liquid chromatography-ESI-triple quadrupole-linear ion trap mass spectrometry. Ten circRNAs, including ciR0159, ciR0212, ciR0153, ciR0149, ciR0016, ciR0044, ciR0022, ciR0381, ciR0006, and ciR0025 were found to be associated with six cannabinoids by weighted gene co-expression network analysis. Twenty-nine of 53 candidate circRNAs, including 9 cannabinoids related were validated successfully using PCR amplification and Sanger sequencing. Taken together, all these results would help to enhance our acknowledge of the regulation of circRNAs, and lay the foundation for breeding new C. sativa cultivars with high cannabinoids through manipulating circRNAs.

CIRCexplorer pipelines for circRNA annotation and quantification from non-polyadenylated RNA-seq datasets

Covalently closed circular RNAs (circRNAs) produced by back-splicing of exon(s) are co-expressed with their cognate linear RNAs from the same gene loci. Most circRNAs are fully overlapped with their cognate linear RNAs in sequences except the back-spliced junction (BSJ) site, thus challenging the computational detection, experimental validation and hence functional evaluation of circRNAs. Nevertheless, specific bioinformatic pipelines were developed to identify fragments mapped to circRNA-featured BSJ sites, and circRNAs were pervasively identified from non-polyadenylated RNA-seq datasets in different cell lines/tissues and across species. Precise identification and quantification of circRNAs provide a basis to further understand their functions. Here, we describe detailed computational steps to annotate and quantify circRNAs using a series of CIRCexplorer pipelines.

Fast and furious: insights of back splicing regulation during nascent RNA synthesis

Alternative splicing of eukaryotic precursor (messenger) RNAs in the nucleus not only increases transcriptomic complexity, but also expands proteomic and functional diversity. In addition to basic types of alternative splicing, recent transcriptome-wide analyses have also suggested other new types of non-canonical splicing, such as back splicing and recursive splicing, and their widespread expression across species Increasing lines of evidence have suggested mechanisms for back splicing, including insights from analyses of nascent RNA sequencing. In this review, we discuss our current understanding of back splicing regulation, and highlight its distinct characteristics in processing during nascent RNA synthesis by taking advantage of metabolic tagging nascent RNA sequencing. Features of recursive splicing are also discussed in the perspective of nascent RNA sequencing.

Mapping circular RNA structures in living cells by SHAPE-MaP

Circular RNAs are produced from back-splicing of exons of precursor mRNAs (pre-mRNAs). The sequences of exons in circular RNAs are identical to their linear cognate mRNAs, but the circular format may confer constraints on their folding and conformation, leading to potentially different functions from their linear RNA cognates. Here, we describe experimental and computational steps that optimize the selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) to probe circular RNA secondary structure at single-nucleotide resolution in living cells.

Characterization of Circular RNAs

Accumulated lines of evidence have revealed that a large number of circular RNAs are produced in transcriptomes from fruit fly to mouse and human. Unlike linear RNAs shaped with 5' cap and 3' tail, circular RNAs are characterized by covalently closed loop structures without open terminals, thus required specific treatments for their identification and validation. Here, we describe a detailed pipeline for the characterization of circular RNAs. It has been successfully applied to the study of circular intronic RNAs (ciRNAs) derived from intron lariats and circular RNAs (circRNAs) produced from back spliced exons in human.

Distinct Processing of lncRNAs Contributes to Non-conserved Functions in Stem Cells

Long noncoding RNAs (lncRNAs) evolve more rapidly than mRNAs. Whether conserved lncRNAs undergo conserved processing, localization, and function remains unexplored. We report differing subcellular localization of lncRNAs in human and mouse embryonic stem cells (ESCs). A significantly higher fraction of lncRNAs is localized in the cytoplasm of hESCs than in mESCs. This turns out to be important for hESC pluripotency. FAST is a positionally conserved lncRNA but is not conserved in its processing and localization. In hESCs, cytoplasm-localized hFAST binds to the WD40 domain of the E3 ubiquitin ligase β-TrCP and blocks its interaction with phosphorylated β-catenin to prevent degradation, leading to activated WNT signaling, required for pluripotency. In contrast, mFast is nuclear retained in mESCs, and its processing is suppressed by the splicing factor PPIE, which is highly expressed in mESCs but not hESCs. These findings reveal that lncRNA processing and localization are previously under-appreciated contributors to the rapid evolution of function.

Altered expression of circular RNA in human dental pulp cells during odontogenic differentiation

The alterations in expression and function of circular RNA (circRNA) in human dental pulp cells (hDPCs) during odontogenic differentiation were investigated. To induce odontogenic differentiation, hDPCs (passage 3) were cultured for 14 days in odontogenic induction medium. circRNA high-throughput sequencing was performed using Illumina HiSeqseq™ 2000. Subsequently, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were used to evaluate the bio-functions of the identified circRNAs. To validate the results of circRNA sequencing, reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was performed for two selected differentially expressed circRNAs. The RNA sequencing results revealed that 1,314 and 1,780 circRNAs were upregulated and downregulated, respectively, during odontogenic induction. Their predicted target miRNAs and genes are involved in several biological functions and signaling pathways, including the mitogen-associated protein kinase signaling pathway. The RT-qPCR results of the two selected circRNAs (hsa_circ_0015260 and hsa_circ_0006984) were consistent with the expression trend obtained using high-throughput sequencing. The results of the present study add to the current understanding of the regulatory mechanisms underlying hDPCs differentiation.

DNA damage and transcriptional regulation in iPSC-derived neurons from Ataxia Telangiectasia patients

Ataxia Telangiectasia (A-T) is neurodegenerative syndrome caused by inherited mutations inactivating the ATM kinase, a master regulator of the DNA damage response (DDR). What makes neurons vulnerable to ATM loss remains unclear. In this study we assessed on human iPSC-derived neurons whether the abnormal accumulation of DNA-Topoisomerase 1 adducts (Top1ccs) found in A-T impairs transcription elongation, thus favoring neurodegeneration. Furthermore, whether neuronal activity-induced immediate early genes (IEGs), a process involving the formation of DNA breaks, is affected by ATM deficiency. We found that Top1cc trapping by CPT induces an ATM-dependent DDR as well as an ATM-independent induction of IEGs and repression especially of long genes. As revealed by nascent RNA sequencing, transcriptional elongation and recovery were found to proceed with the same rate, irrespective of gene length and ATM status. Neuronal activity induced by glutamate receptors stimulation, or membrane depolarization with KCl, triggered a DDR and expression of IEGs, the latter independent of ATM. In unperturbed A-T neurons a set of genes (FN1, DCN, RASGRF1, FZD1, EOMES, SHH, NR2E1) implicated in the development, maintenance and physiology of central nervous system was specifically downregulated, underscoring their potential involvement in the neurodegenerative process in A-T patients.

Coordinated circRNA Biogenesis and Function with NF90/NF110 in Viral Infection

Circular RNAs (circRNAs) generated via back-splicing are enhanced by flanking complementary sequences. Expression levels of circRNAs vary under different conditions, suggesting participation of protein factors in their biogenesis. Using genome-wide siRNA screening that targets all human unique genes and an efficient circRNA expression reporter, we identify double-stranded RNA-binding domain containing immune factors NF90/NF110 as key regulators in circRNA biogenesis. NF90/NF110 promote circRNA production in the nucleus by associating with intronic RNA pairs juxtaposing the circRNA-forming exon(s); they also interact with mature circRNAs in the cytoplasm. Upon viral infection, circRNA expression is decreased, in part owing to the nuclear export of NF90/NF110 to the cytoplasm. Meanwhile, NF90/NF110 released from circRNP complexes bind to viral mRNAs as part of their functions in antiviral immune response. Our results therefore implicate a coordinated regulation of circRNA biogenesis and function by NF90/NF110 in viral infection.

Diverse alternative back-splicing and alternative splicing landscape of circular RNAs

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.

Characterization of Circular RNAs

Accumulated lines of evidence reveal that a large number of circular RNAs are produced in transcriptomes from fruit fly to mouse and human. Unlike linear RNAs shaped with 5' cap and 3' tail, circular RNAs are characterized by covalently closed loop structures without open terminals, thus requiring specific treatments for their identification and validation. Here, we describe a detailed pipeline for the characterization of circular RNAs. It has been successfully applied to the study of circular intronic RNAs derived from intron lariats (ciRNAs) and circular RNAs produced from back spliced exons (circRNAs) in human.