Analysis of Intron Sequences Reveals Hallmarks of Circular RNA Biogenesis in Animals

Combinatorial control of Drosophila circular RNA expression by intronic repeats, hnRNPs, and SR proteins

Kramer et al. discovered that circularization of the Drosophila laccase2 gene is regulated by both intronic repeats and trans-acting splicing factors. Unlike the previously characterized Muscleblind (Mbl) circular RNA, which requires the Mbl protein for its biogenesis, Laccase2 circular RNA levels are not controlled by Mbl or the Laccase2 gene product but rather by multiple hnRNP and SR proteins acting in a combinatorial manner.

Thousands of eukaryotic protein-coding genes are noncanonically spliced to produce circular RNAs. Bioinformatics has indicated that long introns generally flank exons that circularize in Drosophila, but the underlying mechanisms by which these circular RNAs are generated are largely unknown. Here, using extensive mutagenesis of expression plasmids and RNAi screening, we reveal that circularization of the Drosophila laccase2 gene is regulated by both intronic repeats and trans-acting splicing factors. Analogous to what has been observed in humans and mice, base-pairing between highly complementary transposable elements facilitates backsplicing. Long flanking repeats (∼400 nucleotides [nt]) promote circularization cotranscriptionally, whereas pre-mRNAs containing minimal repeats (<40 nt) generate circular RNAs predominately after 3′ end processing. Unlike the previously characterized Muscleblind (Mbl) circular RNA, which requires the Mbl protein for its biogenesis, we found that Laccase2 circular RNA levels are not controlled by Mbl or the Laccase2 gene product but rather by multiple hnRNP (heterogeneous nuclear ribonucleoprotein) and SR (serine–arginine) proteins acting in a combinatorial manner. hnRNP and SR proteins also regulate the expression of other Drosophila circular RNAs, including Plexin A (PlexA), suggesting a common strategy for regulating backsplicing. Furthermore, the laccase2 flanking introns support efficient circularization of diverse exons in Drosophila and human cells, providing a new tool for exploring the functional consequences of circular RNA expression across eukaryotes.

CircNet: a database of circular RNAs derived from transcriptome sequencing data

Circular RNAs (circRNAs) represent a new type of regulatory noncoding RNA that only recently has been identified and cataloged. Emerging evidence indicates that circRNAs exert a new layer of post-transcriptional regulation of gene expression. In this study, we utilized transcriptome sequencing datasets to systematically identify the expression of circRNAs (including known and newly identified ones by our pipeline) in 464 RNA-seq samples, and then constructed the CircNet database (http://circnet.mbc.nctu.edu.tw/) that provides the following resources: (i) novel circRNAs, (ii) integrated miRNA-target networks, (iii) expression profiles of circRNA isoforms, (iv) genomic annotations of circRNA isoforms (e.g. 282 948 exon positions), and (v) sequences of circRNA isoforms. The CircNet database is to our knowledge the first public database that provides tissue-specific circRNA expression profiles and circRNA–miRNA-gene regulatory networks. It not only extends the most up to date catalog of circRNAs but also provides a thorough expression analysis of both previously reported and novel circRNAs. Furthermore, it generates an integrated regulatory network that illustrates the regulation between circRNAs, miRNAs and genes.

NCLscan: accurate identification of non-co-linear transcripts (fusion, trans-splicing and circular RNA) with a good balance between sensitivity and precision

Analysis of RNA-seq data often detects numerous ‘non-co-linear’ (NCL) transcripts, which comprised sequence segments that are topologically inconsistent with their corresponding DNA sequences in the reference genome. However, detection of NCL transcripts involves two major challenges: removal of false positives arising from alignment artifacts and discrimination between different types of NCL transcripts (trans-spliced, circular or fusion transcripts). Here, we developed a new NCL-transcript-detecting method (‘NCLscan’), which utilized a stepwise alignment strategy to almost completely eliminate false calls (>98% precision) without sacrificing true positives, enabling NCLscan outperform 18 other publicly-available tools (including fusion- and circular-RNA-detecting tools) in terms of sensitivity and precision, regardless of the generation strategy of simulated dataset, type of intragenic or intergenic NCL event, read depth of coverage, read length or expression level of NCL transcript. With the high accuracy, NCLscan was applied to distinguishing between trans-spliced, circular and fusion transcripts on the basis of poly(A)- and nonpoly(A)-selected RNA-seq data. We showed that circular RNAs were expressed more ubiquitously, more abundantly and less cell type-specifically than trans-spliced and fusion transcripts. Our study thus describes a robust pipeline for the discovery of NCL transcripts, and sheds light on the fundamental biology of these non-canonical RNA events in human transcriptome.

Splicing noncoding RNAs from the inside out

Eukaryotic precursor-messenger RNAs (pre-mRNAs) undergo splicing to remove intragenic regions (introns) and ligate expressed regions (exons) together. Unlike exons in the mature messenger RNAs (mRNAs) that are used for translation, introns that are spliced out of pre-mRNAs were generally believed to lack function and to be degraded. However, recent studies have revealed that a large group of spliced introns can escape complete degradation and are processed to generate noncoding RNAs (ncRNAs), including different types of small RNAs, long-noncoding RNAs, and circular RNAs. Strikingly, exonic sequences can be also back-spliced from pre-mRNAs to form stable circular RNAs. Together, the findings that ncRNAs can be spliced out of mRNA precursors not only expand the ever-growing repertoire of ncRNAs that originate from different genomic regions, but also reveal the unexpected transcriptomic complexity and functional capacity of eukaryotic genomes.

Elevated RNA Editing Activity Is a Major Contributor to Transcriptomic Diversity in Tumors

Genomic mutations in key genes are known to drive tumorigenesis and have been the focus of much attention in recent years. However, genetic content also may change farther downstream. RNA editing alters the mRNA sequence from its genomic blueprint in a dynamic and flexible way. A few isolated cases of editing alterations in cancer have been reported previously. Here, we provide a transcriptome-wide characterization of RNA editing across hundreds of cancer samples from multiple cancer tissues, and we show that A-to-I editing and the enzymes mediating this modification are significantly altered, usually elevated, in most cancer types. Increased editing activity is found to be associated with patient survival. As is the case with somatic mutations in DNA, most of these newly introduced RNA mutations are likely passengers, but a few may serve as drivers that may be novel candidates for therapeutic and diagnostic purposes.

Widespread noncoding circular RNAs in plants

A large number of noncoding circular RNAs (circRNAs) with regulatory potency have been identified in animals, but little attention has been given to plant circRNAs. We performed genome-wide identification of circRNAs in Oryza sativa and Arabidopsis thaliana using publically available RNA-Seq data, analyzed and compared features of plant and animal circRNAs. circRNAs (12037 and 6012) were identified in Oryza sativa and Arabidopsis thaliana, respectively, with 56% (10/18) of the sampled rice exonic circRNAs validated experimentally. Parent genes of over 700 exonic circRNAs were orthologues between rice and Arabidopsis, suggesting conservation of circRNAs in plants. The introns flanking plant circRNAs were much longer than introns from linear genes, and possessed less repetitive elements and reverse complementary sequences than the flanking introns of animal circRNAs. Plant circRNAs showed diverse expression patterns, and 27 rice exonic circRNAs were found to be differentially expressed under phosphate-sufficient and -starvation conditions. A significantly positive correlation was observed for the expression profiles of some circRNAs and their parent genes. Our results demonstrated that circRNAs are widespread in plants, revealed the common and distinct features of circRNAs between plants and animals, and suggested that circRNAs could be a critical class of noncoding regulators in plants.

The Biological Functions of Non-coding RNAs: From a Line to a Circle

Non-coding RNAs have gained increasing attention, as their physiological and pathological functions are being gradually uncovered. MicroRNAs are the most well-studied ncRNAs, which play essential roles in translational repression and mRNA degradation. In contrast, long non-coding RNAs are distinguished from other small/short non-coding RNAs by length and regulate chromatin remodeling, gene transcription and posttranscriptional modifications. Recently, circular RNAs have emerged as endogenous, abundant, conserved and stable in mammalian cells. It has been demonstrated that circular RNAs can function as miRNA sponges. Other possible biological functions of circular RNAs are still under investigation. In this review, the biogenesis and biological functions of the three major types of ncRNAs, including miRNAs, lncRNAs and circRNAs, are overviewed. In addition, the role of ncRNAs in human diseases and potential clinical applications of ncRNAs are discussed.

Noncoding RNA control of cellular senescence

Senescent cells accumulate in normal tissues with advancing age and arise by long-term culture of primary cells. Senescence develops following exposure to a range of stress-causing agents and broadly influences the physiology and pathology of tissues, organs, and systems in the body. While many proteins are known to control senescence, numerous noncoding (nc)RNAs are also found to promote or repress the senescent phenotype. Here, we review the regulatory ncRNAs (primarily microRNAs and lncRNAs) identified to-date as key modulators of senescence. We highlight the major senescent pathways (p53/p21 and pRB/p16), as well as the senescence-associated secretory phenotype (SASP) and other senescence-associated events governed by ncRNAs, and discuss the importance of understanding comprehensively the ncRNAs implicated in cell senescence.

Gear Up in Circles

Two studies published in this issue of Molecular Cell (Rybak-Wolf et al., 2015) and in the April issue of Nature Neuroscience (You et al., 2015) independently report the upregulated expression of back-spliced circular RNAs (circRNAs) in brains and suggest that they have a potential to regulate synaptic function.

Biogenesis, identification, and function of exonic circular RNAs

Circular RNAs (circRNAs) arise during post-transcriptional processes, in which a single-stranded RNA molecule forms a circle through covalent binding. Previously, circRNA products were often regarded to be splicing intermediates, by-products, or products of aberrant splicing. But recently, rapid advances in high-throughput RNA sequencing (RNA-seq) for global investigation of nonco-linear (NCL) RNAs, which comprised sequence segments that are topologically inconsistent with the reference genome, leads to renewed interest in this type of NCL RNA (i.e., circRNA), especially exonic circRNAs (ecircRNAs). Although the biogenesis and function of ecircRNAs are mostly unknown, some ecircRNAs are abundant, highly expressed, or evolutionarily conserved. Some ecircRNAs have been shown to affect microRNA regulation, and probably play roles in regulating parental gene transcription, cell proliferation, and RNA-binding proteins, indicating their functional potential for development as diagnostic tools. To date, thousands of ecircRNAs have been identified in multiple tissues/cell types from diverse species, through analyses of RNA-seq data. However, the detection of ecircRNA candidates involves several major challenges, including discrimination between ecircRNAs and other types of NCL RNAs (e.g., trans-spliced RNAs and genetic rearrangements); removal of sequencing errors, alignment errors, and in vitro artifacts; and the reconciliation of heterogeneous results arising from the use of different bioinformatics methods or sequencing data generated under different treatments. Such challenges may severely hamper the understanding of ecircRNAs. Herein, we review the biogenesis, identification, properties, and function of ecircRNAs, and discuss some unanswered questions regarding ecircRNAs. We also evaluate the accuracy (in terms of sensitivity and precision) of some well-known circRNA-detecting methods.

What happens at or after transcription: Insights into circRNA biogenesis and function

Circular RNAs (circRNAs) are a large family of noncoding RNAs (ncRNAs) found in metazoans. Systematic studies of circRNAs have just begun. Here, we discuss circRNA biogenesis and functions with a focus on studies indicating great diversification of circRNAs. We highlight the recent identification of a special subtype of circRNAs, called EIciRNAs, and their role in transcriptional regulation. New insights on RNA-RNA interaction and other features associated with circRNA biology are also discussed.

Circular RNAs: Identification, biogenesis and function

Circular RNAs are a novel class of non-coding RNA characterized by the presence of a covalent bond linking the 3' and 5' ends generated by backsplicing. Circular RNAs are widely expressed in a tissue and developmental-stage specific pattern and a subset displays conservation across species. Functional circRNAs have been shown to act as cytoplasmic microRNA sponges and RNA-binding protein sequestering agents as well as nuclear transcriptional regulators, illustrating the relevance of circular RNAs as participants in the regulatory networks governing gene expression. Here, we review the features that characterize circular RNAs, discuss putative circular RNA biogenesis pathways as well as review the uncovered functions of circular RNAs. This article is part of a Special Issue entitled: Clues to long noncoding RNA taxonomy1, edited by Dr. Tetsuro Hirose and Dr. Shinichi Nakagawa.

Biochemical and Transcriptome-Wide Identification of A-to-I RNA Editing Sites by ICE-Seq

Inosine (I) is a modified adenosine (A) in RNA. In Metazoa, I is generated by hydrolytic deamination of A, catalyzed by adenosine deaminase acting RNA (ADAR) in a process called A-to-I RNA editing. A-to-I RNA editing affects various biological processes by modulating gene expression. In addition, dysregulation of A-to-I RNA editing results in pathological consequences. I on RNA strands is converted to guanosine (G) during cDNA synthesis by reverse transcription. Thus, the conventional method used to identify A-to-I RNA editing sites compares cDNA sequences with their corresponding genomic sequences. Combined with deep sequencing, this method has been applied to transcriptome-wide screening of A-to-I RNA editing sites. This approach, however, produces a large number of false positives mainly owing to mapping errors. To address this issue, we developed a biochemical method called inosine chemical erasing (ICE) to reliably identify genuine A-to-I RNA editing sites. In addition, we applied the ICE method combined with RNA-seq, referred to as ICE-seq, to identify transcriptome-wide A-to-I RNA editing sites. In this chapter, we describe the detailed protocol for ICE-seq, which can be applied to various sources and taxa.

The dynamic epitranscriptome: A to I editing modulates genetic information

Adenosine to inosine editing (A to I editing) is a cotranscriptional process that contributes to transcriptome complexity by deamination of adenosines to inosines. Initially, the impact of A to I editing has been described for coding targets in the nervous system. Here, A to I editing leads to recoding and changes of single amino acids since inosine is normally interpreted as guanosine by cellular machines. However, more recently, new roles for A to I editing have emerged: Editing was shown to influence splicing and is found massively in Alu elements. Moreover, A to I editing is required to modulate innate immunity. We summarize the multiple ways in which A to I editing generates transcriptome variability and highlight recent findings in the field.

Statistically based splicing detection reveals neural enrichment and tissue-specific induction of circular RNA during human fetal development

Background

The pervasive expression of circular RNA is a recently discovered feature of gene expression in highly diverged eukaryotes, but the functions of most circular RNAs are still unknown. Computational methods to discover and quantify circular RNA are essential. Moreover, discovering biological contexts where circular RNAs are regulated will shed light on potential functional roles they may play.

Results

We present a new algorithm that increases the sensitivity and specificity of circular RNA detection by discovering and quantifying circular and linear RNA splicing events at both annotated and un-annotated exon boundaries, including intergenic regions of the genome, with high statistical confidence. Unlike approaches that rely on read count and exon homology to determine confidence in prediction of circular RNA expression, our algorithm uses a statistical approach. Using our algorithm, we unveiled striking induction of general and tissue-specific circular RNAs, including in the heart and lung, during human fetal development. We discover regions of the human fetal brain, such as the frontal cortex, with marked enrichment for genes where circular RNA isoforms are dominant.

Conclusions

The vast majority of circular RNA production occurs at major spliceosome splice sites; however, we find the first examples of developmentally induced circular RNAs processed by the minor spliceosome, and an enriched propensity of minor spliceosome donors to splice into circular RNA at un-annotated, rather than annotated, exons. Together, these results suggest a potentially significant role for circular RNA in human development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0690-5) contains supplementary material, which is available to authorized users.

Circular RNA biogenesis can proceed through an exon-containing lariat precursor

Pervasive expression of circular RNA is a recently discovered feature of eukaryotic gene expression programs, yet its function remains largely unknown. The presumed biogenesis of these RNAs involves a non-canonical 'backsplicing' event. Recent studies in mammalian cell culture posit that backsplicing is facilitated by inverted repeats flanking the circularized exon(s). Although such sequence elements are common in mammals, they are rare in lower eukaryotes, making current models insufficient to describe circularization. Through systematic splice site mutagenesis and the identification of splicing intermediates, we show that circular RNA in Schizosaccharomyces pombe is generated through an exon-containing lariat precursor. Furthermore, we have performed high-throughput and comprehensive mutagenesis of a circle-forming exon, which enabled us to discover a systematic effect of exon length on RNA circularization. Our results uncover a mechanism for circular RNA biogenesis that may account for circularization in genes that lack noticeable flanking intronic secondary structure.

Circular RNA: A new star of noncoding RNAs

Circular RNAs (circRNAs) are a novel type of RNA that, unlike linear RNAs, form a covalently closed continuous loop and are highly represented in the eukaryotic transcriptome. Recent studies have discovered thousands of endogenous circRNAs in mammalian cells. CircRNAs are largely generated from exonic or intronic sequences, and reverse complementary sequences or RNA-binding proteins (RBPs) are necessary for circRNA biogenesis. The majority of circRNAs are conserved across species, are stable and resistant to RNase R, and often exhibit tissue/developmental-stage-specific expression. Recent research has revealed that circRNAs can function as microRNA (miRNA) sponges, regulators of splicing and transcription, and modifiers of parental gene expression. Emerging evidence indicates that circRNAs might play important roles in atherosclerotic vascular disease risk, neurological disorders, prion diseases and cancer; exhibit aberrant expression in colorectal cancer (CRC) and pancreatic ductal adenocarcinoma (PDAC); and serve as diagnostic or predictive biomarkers of some diseases. Similar to miRNAs and long noncoding RNAs (lncRNAs), circRNAs are becoming a new research hotspot in the field of RNA and could be widely involved in the processes of life. Herein, we review the formation and properties of circRNAs, their functions, and their potential significance in disease.

RNA editing of non-coding RNA and its role in gene regulation

It has for a long time been known that repetitive elements, particularly Alu sequences in human, are edited by the adenosine deaminases acting on RNA, ADAR, family. The functional interpretation of these events has been even more difficult than that of editing events in coding sequences, but today there is an emerging understanding of their downstream effects. A surprisingly large fraction of the human transcriptome contains inverted Alu repeats, often forming long double stranded structures in RNA transcripts, typically occurring in introns and UTRs of protein coding genes. Alu repeats are also common in other primates, and similar inverted repeats can frequently be found in non-primates, although the latter are less prone to duplex formation. In human, as many as 700,000 Alu elements have been identified as substrates for RNA editing, of which many are edited at several sites. In fact, recent advancements in transcriptome sequencing techniques and bioinformatics have revealed that the human editome comprises at least a hundred million adenosine to inosine (A-to-I) editing sites in Alu sequences. Although substantial additional efforts are required in order to map the editome, already present knowledge provides an excellent starting point for studying cis-regulation of editing. In this review, we will focus on editing of long stem loop structures in the human transcriptome and how it can effect gene expression.

Repetitive elements regulate circular RNA biogenesis

It was long assumed that eukaryotic precursor mRNAs (pre-mRNAs) are almost always spliced to generate a linear mRNA that is subsequently translated to produce a protein. However, it is now clear that thousands of protein-coding genes can be non-canonically spliced to produce circular noncoding RNAs, some of which are expressed at much higher levels than their associated linear mRNAs. How then does the splicing machinery decide whether to generate a linear mRNA or a circular RNA? Recent work has revealed that intronic repetitive elements, including sequences derived from transposons, are critical regulators of this decision. In most cases, circular RNA biogenesis appears to be initiated when complementary sequences from 2 different introns base pair to one another. This brings the splice sites from the intervening exon(s) into close proximity and facilitates the backsplicing event that generates the circular RNA. As many pre-mRNAs contain multiple intronic repeats, distinct circular transcripts can be produced depending on which repeats base pair to one another. Intronic repeats are thus critical regulatory sequences that control the functional output of their host genes, and potentially cause the functions of protein-coding genes to be highly divergent across species.

Widespread exon skipping triggers degradation by nuclear RNA surveillance in fission yeast

Exon skipping is considered a principal mechanism by which eukaryotic cells expand their transcriptome and proteome repertoires, creating different splice variants with distinct cellular functions. Here we analyze RNA-seq data from 116 transcriptomes in fission yeast (Schizosaccharomyces pombe), covering multiple physiological conditions as well as transcriptional and RNA processing mutants. We applied brute-force algorithms to detect all possible exon-skipping events, which were widespread but rare compared to normal splicing events. Exon-skipping events increased in cells deficient for the nuclear exosome or the 5′-3′ exonuclease Dhp1, and also at late stages of meiotic differentiation when nuclear-exosome transcripts decreased. The pervasive exon-skipping transcripts were stochastic, did not increase in specific physiological conditions, and were mostly present at less than one copy per cell, even in the absence of nuclear RNA surveillance and during late meiosis. These exon-skipping transcripts are therefore unlikely to be functional and may reflect splicing errors that are actively removed by nuclear RNA surveillance. The average splicing rate by exon skipping was ∼0.24% in wild type and ∼1.75% in nuclear exonuclease mutants. We also detected approximately 250 circular RNAs derived from single or multiple exons. These circular RNAs were rare and stochastic, although a few became stabilized during quiescence and in splicing mutants. Using an exhaustive search algorithm, we also uncovered thousands of previously unknown splice sites, indicating pervasive splicing; yet most of these splicing variants were cryptic and increased in nuclear degradation mutants. This study highlights widespread but low frequency alternative or aberrant splicing events that are targeted by nuclear RNA surveillance.

RNA circularization strategies in vivo and in vitro

In the plenitude of naturally occurring RNAs, circular RNAs (circRNAs) and their biological role were underestimated for years. However, circRNAs are ubiquitous in all domains of life, including eukaryotes, archaea, bacteria and viruses, where they can fulfill diverse biological functions. Some of those functions, as for example playing a role in the life cycle of viral and viroid genomes or in the maturation of tRNA genes, have been elucidated; other putative functions still remain elusive. Due to the resistance to exonucleases, circRNAs are promising tools for in vivo application as aptamers, trans-cleaving ribozymes or siRNAs. How are circRNAs generated in vivo and what approaches do exist to produce ring-shaped RNAs in vitro? In this review we illustrate the occurrence and mechanisms of RNA circularization in vivo, survey methods for the generation of circRNA in vitro and provide appropriate protocols.

Identification of the specific interactors of the human lariat RNA debranching enzyme 1 protein

In eukaryotes, pre-mRNA splicing is an essential step for gene expression. We have been analyzing post-splicing intron turnover steps in higher eukaryotes. Here, we report protein interaction between human Debranching enzyme 1 (hDbr1) and several factors found in the Intron Large (IL) complex, which is an intermediate complex of the intron degradation pathway. The hDbr1 protein specifically interacts with xeroderma pigmentosum, complementeation group A (XPA)-binding protein 2 (Xab2). We also attempted to identify specific interactors of hDbr1. Co-immunoprecipitation experiments followed by mass spectrometry analysis identified a novel protein as one of the specific interactors of hDbr1. This protein is well conserved among many species and shows the highest similarity to yeast Drn1, so it is designated as human Dbr1 associated ribonuclease 1 (hDrn1). hDrn1 directly interacts with hDbr1 through protein–protein interaction. Furthermore, hDrn1 shuttles between the nucleus and the cytoplasm, as hDbr1 protein does. These findings suggest that hDrn1 has roles in both the nucleus and the cytoplasm, which are highly likely to involve hDbr1.