Circular RNA

Circular RNA (or circRNA) is a type of RNA which, unlike the better known linear RNA, forms a covalently closed continuous loop, i.e., in circular RNA the 3' and 5' ends normally present in an RNA molecule have been joined together. This feature confers numerous properties to circular RNAs, many of which have only recently been identified.

Many circular RNAs arise from otherwise protein-coding genes, but circular RNAs produced in the cell have not been shown to code for proteins. They have therefore been categorized as noncoding RNA. Some circular RNAs have recently shown potential as gene regulators. Like many other alternative noncoding isoforms, the biological function of most circular RNAs are unclear.

Because circular RNAs do not have 5' or 3' ends, they are resistant to exonuclease-mediated degradation and are presumably more stable than most linear RNAs in cells.^[1]

RNA splicing

In contrast to genes in bacteria, eukaryotic genes are split by non-coding sequences known as introns. In eukaryotes, as a gene is transcribed from DNA into a messenger RNA (mRNA) transcript, intervening introns are removed, leaving only exons in the mature mRNA, which can subsequently be translated to produce the protein product.^[2] The spliceosome,^[2] a protein-RNA complex located in the nucleus, catalyzes splicing in the following manner:

1. The spliceosome recognizes an intron, which is flanked by specific sequences at its 5’ and 3’ ends, known as a donor splice site (or 5’ splice site) and an acceptor splice site (or 3’ splice site), respectively.
2. The 5’ splice site sequence is then subjected to a nucleophilic attack by a downstream sequence called the branch point, resulting in a circular structure called a lariat.
3. The free 5’ exon then attacks the 3’ splice site, joining the two exons and releasing the intron lariat. The intron lariat is subsequently debranched and quickly degraded.^[2]

Pre-mRNA to mRNA splicing

Alternative splicing

Alternative splicing is a phenomenon through which one RNA transcript can yield different protein products based on which segments are considered "introns" and which are considered "exons" during each splicing event.^[2] Although not specific for humans, it is a partial explanation for the fact that human and other, much more simple species (such as nematodes) have a similar number of genes (in the range of 20 - 25 thousand).^[3] One of the most striking examples of alternative splicing is in the Drosophila DSCAM gene. This single gene can give rise to approximately 30 thousand distinct alternatively spliced isoforms.^[4]

Non-canonical splicing

Exon scrambling

Exon scrambling, also called exon shuffling, describes an event in which exons are spliced in a "non-canonical" (atypical) order. There are three ways in which exon scrambling can occur:

1. Tandem exon duplication in the genome, which often occurs in cancers
2. Trans-splicing (in which two RNA transcripts fuse), which results in a linear transcript that contains exons that, for example, may be derived from genes encoded on two different chromosomes. Trans-splicing is very common in C. elegans
3. A splice donor site being joined to a splice acceptor site further upstream in the primary transcript, yielding a circular transcript.^[5]

The notion that circularized transcripts are byproducts from imperfect splicing is supported by the low abundance and the lack of sequence conservation of most circRNAs,^[6] but has recently been challenged.^[5][7][8]

Characteristics of circular RNA

Genome-wide identification of circRNAs

Several research groups have characterized circular RNAs by sorting through vast collections of RNA sequencing data.^[5][6][8][9] Since circRNAs are non-polyadenylated, poly(A)-selected RNA-seq data cannot be used for circRNA discovery. To identify potential circular RNA isoforms, one typically searches for sequencing reads showing a junction between two “scrambled” exons.

A brief summary of each group’s results follows:

Salzman et al. 2012^[5]

• Originally wanted to identify cancer-specific exon scrambling events
• Ended up finding scrambled exons in a large number of both normal and cancer cells
• Scrambled exon isoforms: ~10% of total transcript isoforms in leukocytes
• Identified 2,748 scrambled isoforms in HeLa and H9 embryonic stem cells
• About 1 in 50 expressed genes produced scrambled transcript isoforms at least 10% of the time
• Some tests for circularity: (1) Treated samples with RNase R, an enzyme which degrades linear RNAs but not circular RNAs, and (2) Tested for the presence of poly-A tails (shouldn’t be present in a circular molecule)
• Conclusion: 98% of scrambled isoforms represented circRNAs

Jeck et al. 2013^[8]

• Treated human fibroblast RNA with RNase R to enrich for circular RNAs
• Used three "stringency" categories (low, medium, high) to classify circular transcripts based on their levels of abundance
• Including the "low" category, ~1 in 8 expressed genes produced detectable levels of circRNA
• Significantly higher than Salzman’s number (above)
• May be due to greater sequencing depth

Memczak et al. 2013^[9]

• Developed a computational method to detect circRNAs
• de novo detected circRNAs in humans, mouse and C. elegans and extensively validated them
• Found that circRNAs are often expressed tissue/developmental stage specific
• Described that circRNAs can act as antagonists of miRNAs as exemplified by the circRNA CDR1as (see below)

Guo et al. 2014^[6]

• Identified and quantified human circRNAs from ENCODE Ribozero RNA-seq data
• Most circRNAs are minor splice isoforms and are expressed only in a few cell types
• 7,112 human circRNAs have circular fractions of at least 10%
• circRNAs are no more conserved than their linear controls
• Ribosome profiling indicates that circRNAs are not translated
• Except for CDR1as, very few circRNAs have the potential to act as microRNA sponges
• Concluded that most circRNAs are inconsequential side-products of imperfect splicing

Zhang et al. 2014^[10]

• Developed CIRCexplorer to identify thousands of circRNAs in human with p(A)－ w/o RNase R RNA-seq data
• The vast majority of highly expressed exonic circular RNAs identified were processed from exons located in the middle of RefSeq genes, suggesting that the circular RNA formation is generally coupled to RNA splicing.
• Most circular RNAs contain multiple exons, most commonly two to three exons.
• Exons from circular RNAs with only one circularized exon were much longer than those from circular RNAs with multiple circularised exons, indicating that processing may prefer a certain length to maximize exon(s) circularization.
• Alu elements that could form IRAlus, either convergent or divergent, are juxtaposed across flanking introns of circRNAs in a parallel way with similar distances to adjacent exons.
• Circular RNA formation is promoted by IRAlus or by other nonrepetitive but complementary sequences.
• Exon circularization efficiency is affected by competition of RNA pairing.
• Alternative RNA pairing and their competition lead to alternative circularization.
• Both exon circularization and its regulation are evolutionarily dynamic.

Length of circRNAs

A recent study of human circRNAs revealed that these molecules are usually composed of 1–5 exons.^[9] Each of these exons can be up to 3x longer than the average expressed exon,^[8] suggesting that exon length may play a role in deciding which exons to circularize. 85% of circularized exons overlap with exons that code for protein,^[9] although the circular RNAs themselves do not appear to be translated. During circRNA formation, exon 2 is often the upstream "acceptor" exon.^[5]

Introns surrounding exons that are selected to be circularized are, on average, up to 3x longer than those not flanking pre-circle exons,^[5][8] although it is not yet clear why this is the case. Compared to regions not resulting in circles, these introns are much more likely to contain complementary inverted Alu repeats; Alu is the most common transposon in the genome.^[8] By the Alu repeats base pairing to one another, it has been proposed that this may enable the splice sites to find each other, thus facilitating circularization.^[7][8]

Introns within the circRNAs are retained at a relatively high frequency (~25%),^[6] thus adding extra sequence to the mature circRNAs.

Location of circRNAs in the cell

In the cell, circRNAs are predominantly found in the cytoplasm, where the number of circular RNA transcripts derived from a gene can be up to ten times greater than the number of associated linear RNAs generated from that locus. It is unclear how circular RNAs exit the nucleus through a relatively small nuclear pore. Because the nuclear envelope breaks down during mitosis, one hypothesis is that the molecules exit the nucleus during this phase of the cell cycle.^[8] However, certain circRNAs, such as CiRS-7/CDR1as, are expressed in neuronal tissues,^[9][11] where mitotic division is not prevalent.

Human cell nucleus diagram

CircRNAs are stable compared to linear RNAs

CircRNAs lack polyadenylated tail and therefore are predicted to be less prone to degradation by exonucleases. In 2015, Enuka et al. measured the half-lives of 60 circRNAs and their linear counterparts expressed from the same host gene and revealed that the median half-life of circRNAs of mammary cells (18.8–23.7 h) is at least 2.5 times longer than the median half-life of their linear counterparts (4.0–7.4 h).^[1] Generally, the lifetime of RNA molecules defines their response time.^[12] Accordingly, it was reported that mammary circRNAs respond slowly to stimulation by growth factors.^[1]

Plausible functions of circular RNA

Evolutionary conservation of circularization mechanisms and signals

CircRNAs have been identified in various species across the domains of life. In 2011, Danan et al. sequenced RNA from Archaea. After digesting total RNA with RNase R, they were able to identify circular species, indicating that circRNAs are not specific to eukaryotes.^[13] However, these archaeal circular species are probably not made via splicing, suggesting that other mechanisms to generate circular RNA likely exist.

Three domains of life

In a closer evolutionary connection, a comparison of RNA from mouse testes vs. RNA from a human cell found 69 orthologous circRNAs. For example, both humans and mice encode the HIPK2 and HIPK3 genes, two paralogous kinases which produce a large amount of circRNA from one particular exon in both species.^[8] Evolutionary conservation reinforces the likelihood of a relevant and significant role for RNA circularization.

CDR1as/CiRS-7 as a miR-7 sponge

microRNAs (miRNAs) are small (~21nt) non-coding RNAs that repress translation of messenger RNAs involved in a large, diverse set of biological processes.^[14] They directly base-pair to target messenger RNAs (mRNAs), and can trigger cleavage of the mRNA depending on the degree of complementarity.

MicroRNAs are grouped in “seed families.” Family members share nucleotides #2–7, called the seed region.^[15] Argonaute proteins are the “effector proteins” which help miRNAs carry out their job, while microRNA sponges are RNAs that “sponge up” miRNAs of a particular family, thereby serving as competitive inhibitors that suppress the ability of the miRNA to bind its mRNA targets, thanks to the presence of multiple binding sites that recognize a specific seed region.^[15] Certain circular RNAs have many miRNA binding sites, which yielded a clue that they may function in sponging. Two recent papers confirmed this hypothesis by investigating a circular sponge called CDR1as/CiRS-7 in detail.

CDR1as/CiRS-7 is encoded in the genome antisense to the human CDR1 (gene) locus (hence the name CDR1as),^[9] and targets miR-7 (hence the name CiRS-7 – Circular RNA Sponge for miR-7).^[11] It has over 60 miR-7 binding sites, far more than any known linear miRNA sponge.^[9][11]

AGO2 is miR-7’s associated Argonaute protein (see above). Though CDR1as/CiRS-7 can be cleaved by miR-671 and its associated Argonaute protein,^[11] it cannot be cleaved by miR-7 and AGO2. MicroRNA cleavage activity depends on complementarity beyond the 12th nucleotide position; none of CiRS-7’s binding sites meet this requirement.

An experiment with zebrafish, which do not have the CDR1 locus in their genome, provides evidence for CiRS-7’s sponge activity. During development, miR-7 is strongly expressed in the zebrafish brain. To silence miR-7 expression in zebrafish, Memczak and colleagues took advantage of a tool called morpholino, which can base pair and sequester target molecules.^[16] Interestingly, morpholino treatment had the same severe effect on midbrain development as ectopically expressing CiRS-7 in zebrafish brains using injected plasmids. This indicates a significant interaction between CiRS-7 and miR-7 in vivo.^[9]

Another notable circular miRNA sponge is SRY. SRY, which is highly expressed in murine testes, functions as a miR-138 sponge.^[11][17] In the genome, SRY is flanked by long inverted repeats (IRs) over 15.5 kilobases (kb) in length. When one or both of the IRs are deleted, circularization does not occur. It was this finding that introduced the idea of inverted repeats enabling circularization.^[18]

Because circular RNA sponges are characterized by high expression levels, stability, and a large number of miRNA binding sites, they are likely to be more effective sponges than those that are linear.^[7]

Other possible functions for circRNAs

Though recent attention has been focused on circRNA’s “sponge” functions, scientists are considering several other functional possibilities as well. For example, some areas of the mouse adult hippocampus show expression of CiRS-7 but not miR-7, suggesting that CiRS-7 may have roles that are independent of interacting with the miRNA.^[9]

Potential roles include the following:

• Binding to RNA-binding proteins (RBPs) and RNAs besides miRNAs to form RNA-protein complexes.^[7] These complexes could regulate RBP & RNA interactions with, for example, the canonical linear transcript of the gene.^[5]
• Protein production
  • Chen and Sarnow 1995 showed that a synthetic circRNA that contained an IRES (internal ribosome entry site) produced a protein product ‘‘in vitro’’, whereas that without an IRES did not. Although the tested circRNA was a purely artificial construct, Chen and Sarnow stated in their paper that they would be interested to see whether circles naturally contain IRES elements.^[19]
  • Jeck et al. 2013: Tested natural circRNAs that contained a translation "start codon." However, none of these molecules bound to ribosomes, suggesting that many circRNAs may not be translated ‘‘in vivo’’.^[8]
• Transporting miRNAs inside the cell. The fact that CiRS-7 can be sliced by miR-671 might indicate the existence of a system to release a “load” of miRNAs at the appropriate time.^[20]
• Regulating mRNA in the cell through limited base pairing. It is formally possible that miR-7 moderates CiRS-7’s regulatory activity instead of the other way around!^[9][20]

Circular intronic long non-coding RNAs (ciRNAs)

Usually, intronic lariats (see above) are debranched and rapidly degraded. However, a debranching failure can lead to the formation of circular intronic long non-coding RNAs, also known as ciRNAs.^[21] CiRNA formation, rather than being a random process, seems to depend on the presence of specific elements near the 5’ splice site and the branchpoint site (see above).

CiRNAs are distinct from circRNAs in that they are prominently found in the nucleus rather than the cytoplasm. In addition, these molecules contain few (if any) miRNA binding sites. Instead of acting as sponges, ciRNAs seem to function in regulating the expression of their parent genes. For example, a relatively abundant ciRNA called ci-ankrd52 positively regulates Pol II transcription. Many ciRNAs remain at their “sites of synthesis” in the nucleus. However, ciRNA may have roles other than simply regulating their parent genes, as ciRNAs do localize to additional sites in the nucleus other than their “sites of synthesis”.^[21]

Circular RNA and disease

As with most topics in molecular biology, it is important to consider how circular RNA can be used as a tool to help mankind. Given its (1) abundance, (2) evolutionary conservation, and (3) potential regulatory role, it is worthwhile to look into how circular RNA can be used to study pathogenesis and devise therapeutic interventions. For example:

• Circular ANRIL (cANRIL) is the circular form of ANRIL, a long non-coding RNA (ncRNA). Expression of cANRIL is correlated with risk for atherosclerosis, a disease in which the arteries become hard. It has been proposed that cANRIL can modify INK4/ARF expression, which, in turn, increases risk for atherosclerosis.^[22] Further study of cANRIL expression could potentially be used to prevent or treat atherosclerosis.
• miR-7 plays an important regulatory role in several cancers and in Parkinson's disease, which is a degenerative neuronal illness.^[11] Perhaps CiRS-7’s sponge activity could help in countering miR-7 activity. If circular sponge activity can indeed help in countering harmful miRNA activity, scientists will need to figure out the best way to introduce sponge expression, perhaps via a transgene, which is a synthetic gene that is transferred between organisms. It is also important to consider how transgenes can be expressed only in specific tissues, or expressed only when induced.^[15]

Viroids as circular RNAs

Viroids are mostly plant pathogens, which consist of short stretches (a few hundred nucleobases) of highly complementary, circular, single-stranded, and non-coding RNAs without a protein coat. Compared with other infectious plant pathogens, viroids are extremely small in size, ranging from 246 to 467 nucleobases; they thus consist of fewer than 10,000 atoms. In comparison, the genome of the smallest known viruses capable of causing an infection by themselves are around 2,000 nucleobases long.^[23]

Databases

Recently were created various circRNA databases for the identification and study of this kind of RNAs.

References

External links

• CIRCpedia database^[1]
• circRNABase database ^[2]
• circBase.org database – recently published comprehensive database of circular RNAs^[3]
• Circ2Traits database

1. ↑ "CIRCpedia". www.picb.ac.cn. Retrieved 2016-08-09. 
2. ↑ Li, JH; Liu, S; Zhou, H; Qu, LH; Yang, JH (January 2014). "starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data.". Nucleic Acids Research. 42 (Database issue): D92–7. doi:10.1093/nar/gkt1248. PMID 24297251. 
3. ↑ Glazar, Petar (2014). "circBase: a database for circular RNAs.". RNA. 20 (11): 1666–70. doi:10.1261/rna.043687.113. PMID 25234927.