Circular RNAs are abundant, conserved, and associated with ALU repeats

Retrozymes are a unique family of non-autonomous retrotransposons with hammerhead ribozymes that propagate in plants through circular RNAs

Background

Catalytic RNAs, or ribozymes, are regarded as fossils of a prebiotic RNA world that have remained in the genomes of modern organisms. The simplest ribozymes are the small self-cleaving RNAs, like the hammerhead ribozyme, which have been historically considered biological oddities restricted to some RNA pathogens. Recent data, however, indicate that small self-cleaving ribozymes are widespread in genomes, although their functions are still unknown.

Results

We reveal that hammerhead ribozyme sequences in plant genomes form part of a new family of small non-autonomous retrotransposons with hammerhead ribozymes, referred to as retrozymes. These elements contain two long terminal repeats of approximately 350 bp, each harbouring a hammerhead ribozyme that delimitates a variable region of 600–1000 bp with no coding capacity. Retrozymes are actively transcribed, which gives rise to heterogeneous linear and circular RNAs that accumulate differentially depending on the tissue or developmental stage of the plant. Genomic and transcriptomic retrozyme sequences are highly heterogeneous and share almost no sequence homology among species except the hammerhead ribozyme motif and two small conserved domains typical of Ty3-gypsy long terminal repeat retrotransposons. Moreover, we detected the presence of RNAs of both retrozyme polarities, which suggests events of independent RNA-RNA rolling-circle replication and evolution, similarly to that of infectious circular RNAs like viroids and viral satellite RNAs.

Conclusions

Our work reveals that circular RNAs with hammerhead ribozymes are frequently occurring molecules in plant and, most likely, metazoan transcriptomes, which explains the ubiquity of these genomic ribozymes and suggests a feasible source for the emergence of circular RNA plant pathogens.

Electronic supplementary material

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

LncRNA, a new component of expanding RNA-protein regulatory network important for animal sperm development

Spermatogenesis is one of the fundamental processes of sexual reproduction, present in almost all metazoan animals. Like many other reproductive traits, developmental features and traits of spermatogenesis are under strong selective pressure to change, both at morphological and underlying molecular levels. Yet evidence suggests that some fundamental features of spermatogenesis may be ancient and conserved among metazoan species. Identifying the underlying conserved molecular mechanisms could reveal core components of metazoan spermatogenic machinery and provide novel insight into causes of human infertility. Conserved RNA-binding proteins and their interacting RNA network emerge to be a common theme important for animal sperm development. We review research on the recent addition to the RNA family - Long non-coding RNA (lncRNA) and its roles in spermatogenesis in the context of the expanding RNA-protein network.

Altered expression pattern of circular RNAs in primary and metastatic sites of epithelial ovarian carcinoma

Recently, a class of endogenous species of RNA called circular RNA (circRNA) has been shown to regulate gene expression in mammals and their role in cellular function is just beginning to be understood. To investigate the role of circRNAs in ovarian cancer, we performed paired-end RNA sequencing of primary sites, peritoneal and lymph node metastases from three patients with stage IIIC ovarian cancer. We developed an in-house computational pipeline to identify and characterize the circRNA expression from paired-end RNA-Seq libraries. This pipeline revealed thousands of circular isoforms in Epithelial Ovarian Carcinoma (EOC). These circRNAs are enriched for potentially effective miRNA seed matches. A significantly larger number of circRNAs are differentially expressed between tumor sites than mRNAs. Circular and linear expression exhibits an inverse trend for many cancer related pathways and signaling pathways like NFkB, PI3k/AKT and TGF-β typically activated for mRNA in metastases are inhibited for circRNA expression. Further, circRNAs show a more robust expression pattern across patients than mRNA forms indicating their suitability as biomarkers in highly heterogeneous cancer transcriptomes. The consistency of circular RNA expression may offer new candidates for cancer treatment and prognosis.

Reprogramming, Circular Reasoning and Self versus Non-self: One-Stop Shopping with RNA Editing

Transcription of genetic information from archival DNA into RNA molecule working copies is vital for proper cellular function and is highly accurate. In turn, RNAs serve structural, enzymatic, and regulatory roles, as well as being informational templates for the ribosomal translation of proteins. Following RNA synthesis, maturing of RNA molecules occurs through various RNA processing events. One component of the collection of processes involving RNA species, broadly defined as RNA metabolism, is the RNA-editing pathway and is found in all animals. Acting specifically on RNA substrates with double-stranded character, RNA editing has been shown to regulate a plethora of genomic outputs, including gene recoding, RNA splicing, biogenesis and targeting actions of microRNAs and small interfering RNAs, and global gene expression. Recent evidence suggests that RNA modifications mediated via RNA editing influence the biogenesis of circular RNAs and safeguard against aberrant innate immune responses generated to endogenous RNA sources. These novel roles have the potential to contribute new insights into molecular mechanisms underlying pathogenesis mediated by mishandling of double-stranded RNA. Here, we discuss recent advances in the field, which highlight novel roles associated with the RNA-editing process and emphasize their importance during cellular RNA metabolism. In addition, we highlight the relevance of these newly discovered roles in the context of neurological disorders and the more general concept of innate recognition of self versus non-self.

Circular RNAs as potential biomarkers for cancer diagnosis and therapy

Circular RNAs (circRNAs) are a naturally occurring type of universal and diverse endogenous noncoding RNAs which unlike linear RNAs, have covalently linked ends. They are usually stable, abundant, conserved RNA molecules and often exhibit tissue/developmental-stage specific expression. Functional circRNAs have been identified to act as microRNA sponges and RNA-binding protein (RBP) sequestering agents as well as transcriptional regulators. These multiple functional roles elicit a great potential for circRNAs in biological applications. Emerging evidence shows that circRNAs play important roles in several diseases, particularly in cancer where they act through regulating protein expression of the pivotal genes that are critical for carcinogenesis. The presence of abundant circRNAs in saliva, exosomes and clinical standard blood samples will make them potential diagnostic or predictive biomarkers for diseases, particularly for cancer development, progression and prognosis. Here, we review the current literature and provide evidence for the impact of circRNAs in cancers and their potential significance in cancer prognosis and clinical treatment.

Lessons from non-canonical splicing

Recent improvements in experimental and computational techniques that are used to study the transcriptome have enabled an unprecedented view of RNA processing, revealing many previously unknown non-canonical splicing events. This includes cryptic events located far from the currently annotated exons and unconventional splicing mechanisms that have important roles in regulating gene expression. These non-canonical splicing events are a major source of newly emerging transcripts during evolution, especially when they involve sequences derived from transposable elements. They are therefore under precise regulation and quality control, which minimizes their potential to disrupt gene expression. We explain how non-canonical splicing can lead to aberrant transcripts that cause many diseases, and also how it can be exploited for new therapeutic strategies.

Circular RNA expression in cutaneous squamous cell carcinoma

BACKGROUND

CircularRNAs (circRNAs) are a reinvented class of abundant, stable, and evolutionary conserved non-coding RNAs with pivotal impacts on the cellular regulatory network and epigenetics by sequestering microRNAs (miRNAs) like a sponge.

OBJECTIVE

Purpose of the present study was to investigate circRNA expression in cutaneous squamous cell carcinoma (cSCC).

METHODS

A total of six cSCC and six non-lesional skin (control) biopsies were harvested. Microarray based circRNA expression was determined in the cSCC (n=3) and compared with the non-lesional skin (n=3) from a group of 13,617 distinct human circRNAs found in the Arraystar circRNA Array V2.0 (Arraystar, Rockville, USA). Microarray data were validated by quantitative real-time reverse transcription polymerase chain reaction in a separate group (cSCC, n=3 and non-lesional skin, n=3). miRNA binding to miRNA response elements (MREs) sequence data were acquired bioinformatically. Further data mining was performed to identify circRNAs containing MRE sequences that interacted with previously described miRNAs playing a role in cSCC formation.

RESULTS

A total of 322 circRNAs (143 up- and 179 down-regulated; fold change ≥2 and p<0.05) were identified as differentially expressed in cSCC. Furthermore, we identified a total of 1603 MREs that were part of the differentially expressed circRNAs. Among those circRNAs, a complementary MRE sequence was identified in 23 miRNAs previously known to be cSCC relevant.

CONCLUSION

This study showed that circRNAs are differentially expressed in cSCC and play an important role in tumor formation by interfering with cSCC relevant miRNAs through miRNA sequence complementary MREs participating in epigenetic control.

Comprehensive Circular RNA Profiling Reveals That hsa_circ_0005075, a New Circular RNA Biomarker, Is Involved in Hepatocellular Crcinoma Development

There is increasing evidence that circular RNAs (circRNAs) are involved in cancer development; however, their role in hepatocellular carcinoma (HCC) remains unclear. Here, we aimed to determine the circRNA expression profile in HCC, and investigate relevant mechanisms for cancer progression. The global circRNA expression profile between HCC (n = 3) and adjacent normal liver (n = 3) tissue was significantly different. Three circRNAs (hsa_circ_0000520, hsa_circ_0005075, and hsa_circ_0066444) showed significantly different expression levels in HCC tissues, which were further validated in 60 matched tissue samples using real-time qRT-PCR. Only hsa_circ_0005075 exhibited significant difference in expression (P <0.001) between HCC and normal tissues. Hsa_circ_0005075 expression correlated with HCC tumor size (P = 0.042), and showed good diagnostic potential (AUROC = 0.94). Finally, we constructed a network of hsa_circ_0005075-targeted miRNA-gene interactions, including miR-23b-5p, miR-93-3p, miR-581, miR-23a-5p, and their corresponding mRNAs. Gene oncology analysis revealed that hsa_circ_0005075 could participate in cell adhesion during HCC development. In summary, we identified hsa_circ_0005075 as a potential HCC biomarker; however, further studies are required to confirm the role of this circRNA, and others, in HCC development.

Circular RNA expression in basal cell carcinoma

Aim: Circular RNAs (circRNAs), are nonprotein coding RNAs consisting of a circular loop with multiple miRNA, binding sites called miRNA response elements (MREs), functioning as miRNA sponges. This study was performed to identify differentially expressed circRNAs and their MREs in basal cell carcinoma (BCC). Materials & methods: Microarray circRNA expression profiles were acquired from BCC and control followed by qRT-PCR validation. Bioinformatical target prediction revealed multiple MREs. Sequence analysis was performed concerning MRE interaction potential with the BCC miRNome. Results: We identified 23 upregulated and 48 downregulated circRNAs with 354 miRNA response elements capable of sequestering miRNA target sequences of the BCC miRNome. Conclusion: The present study describes a variety of circRNAs that are potentially involved in the molecular pathogenesis of BCC.

Profiling and Validation of the Circular RNA Repertoire in Adult Murine Hearts

For several decades, cardiovascular disease has been the leading cause of death throughout all countries. There is a strong genetic component to many disease subtypes (e.g., cardiomyopathy) and we are just beginning to understand the relevant genetic factors. Several studies have related RNA splicing to cardiovascular disease and circular RNAs (circRNAs) are an emerging player. circRNAs, which originate through back-splicing events from primary transcripts, are resistant to exonucleases and typically not polyadenylated. Initial functional studies show clear phenotypic outcomes for selected circRNAs. We provide, for the first time, a comprehensive catalogue of RNase R-resistant circRNA species for the adult murine heart. This work combines state-of-the-art circle sequencing with our novel DCC software to explore the circRNA landscape of heart tissue. Overall, we identified 575 circRNA species that pass a beta-binomial test for enrichment (false discovery rate of 1%) in the exonuclease-treated sequencing sample. Several circRNAs can be directly attributed to host genes that have been previously described as associated with cardiovascular disease. Further studies of these candidate circRNAs may reveal disease-relevant properties or functions of specific circRNAs.

HIPK family kinases bind and regulate the function of the CCR4-NOT complex

Down-regulation of the HIPK interactor CNOT2 leads to reduced HIPK2 protein levels, identifying the CCR4-NOT complex as a new regulator of HIPK2 abundance. Functional assays reveal that HIPK2 and HIPK1 restrict CNOT2-dependent mRNA decay, thus extending the regulatory potential of these kinases to the level of posttranscriptional gene regulation.

The serine/threonine kinase HIPK2 functions as a regulator of developmental processes and as a signal integrator of a wide variety of stress signals, such as DNA damage, hypoxia, and reactive oxygen intermediates. Because the kinase is generated in a constitutively active form, its expression levels are restricted by a variety of different mechanisms. Here we identify the CCR4-NOT complex as a new regulator of HIPK2 abundance. Down-regulation or knockout of the CCR4-NOT complex member CNOT2 leads to reduced HIPK2 protein levels without affecting the expression level of HIPK1 or HIPK3. A fraction of all HIPK family members associates with the CCR4-NOT components CNOT2 and CNOT3. HIPKs also phosphorylate the CCR4-NOT complex, a feature that is shared with their yeast progenitor kinase, YAK1. Functional assays reveal that HIPK2 and HIPK1 restrict CNOT2-dependent mRNA decay. HIPKs are well known regulators of transcription, but the mutual regulation between CCR4-NOT and HIPKs extends the regulatory potential of these kinases by enabling posttranscriptional gene regulation.

The Inescapable Influence of Noncoding RNAs in Cancer

This report summarizes information presented at the 2015 Keystone Symposium on "MicroRNAs and Noncoding RNAs in Cancer." Nearly two decades after the discovery of the first miRNA, the role of noncoding RNAs in developmental processes and the mechanisms behind their dysregulation in cancer has been steadily elucidated. Excitingly, miRNAs have begun making their way into the clinic to combat diseases such as hepatitis C and various forms of cancer. Therefore, at this Keystone meeting, novel findings were presented that enhance our view on how small and long noncoding RNAs control developmental timing and oncogenic processes. Recurring themes included (i) how miRNAs can be differentially processed, degraded, and regulated by ribonucleoprotein complexes, (ii) how particular miRNA genetic networks that control developmental process, when disrupted, can result in cancer disease, (iii) the technologies available to therapeutically deliver RNA to combat diseases such as cancer, and (iv) the elucidation of the mechanism of actions for long noncoding RNAs, currently a poorly understood class of noncoding RNA. During the meeting, there was an emphasis on presenting unpublished findings, and the breadth of topics covered reflected how inescapable the influence of noncoding RNAs is in development and cancer.

The insights of Let-7 miRNAs in oncogenesis and stem cell potency

The ability of the classic tumour‐suppressive let‐7 family to inhibit carcinogenesis, tumour progression, recurrence and pluripotency of cancer stem cells has generated significant interest in the field of cancer research. Through suppressing and degrading downstream‐targeted mRNAs, let‐7 affected most aspects of cell biology. It is perplexing how let‐7 affects oncogenesis, as the large influx of new miRNAs and other kinds of non‐coding RNAs are continuously defined. In this review, we delineate the complex functions of let‐7 and discuss the future direction of let‐7 research.

Roles of Circular RNAs in Neurologic Disease

Circular RNAs (circRNAs) are a novel type of endogenous noncoding RNA receiving increasing attention. They have been shown to act as a natural microRNA sponges that repress the activity of corresponding miRNAs by binding with them, thus regulating target genes. Numerous studies have shown that miRNAs are involved in the pathogenesis of neurological diseases. Therefore, circRNAs may act as important regulatory factors in the occurrence and development processes of neurological disease.

Circular RNAs in Eukaryotic Cells

Circular RNAs (circRNAs) are now recognized as large species of transcripts in eukaryotic cells. From model organisms such as C. elegans, Drosophila, mice to human beings, thousands of circRNAs formed from back-splicing of exons have been identified. The known complexity of transcriptome has been greatly expanded upon the discovery of these RNAs. Studies about the biogenesis and physiological functions have yielded substantial knowledge for the circRNAs, and they are now more likely to be viewed as regulatory elements coded by the genome rather than unavoidable noise of gene expression. Certain human diseases may also relate to circRNAs. These circRNAs show diversifications in features such as sequence composition and cellular localization, and thus we propose that they may be divided into subtypes such as cytoplasmic circRNAs, nuclear circRNAs, and exon-intron circRNAs (EIciRNAs). Here we summarize and discuss knowns and unknowns for these RNAs, and we need to keep in mind that the whole field is still at the beginning of exciting explorations.

Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs

Circular RNAs (circRNAs) represent a class of widespread and diverse endogenous RNAs that may regulate gene expression in eukaryotes. However, the regulation and function of human circRNAs remain largely unknown. Here we generate ribosomal-depleted RNA sequencing data from six normal tissues and seven cancers, and detect at least 27,000 circRNA candidates. Many of these circRNAs are differently expressed between the normal and cancerous tissues. We further characterize one abundant circRNA derived from Exon2 of the HIPK3 gene, termed circHIPK3. The silencing of circHIPK3 but not HIPK3 mRNA significantly inhibits human cell growth. Via a luciferase screening assay, circHIPK3 is observed to sponge to 9 miRNAs with 18 potential binding sites. Specifically, we show that circHIPK3 directly binds to miR-124 and inhibits miR-124 activity. Our results provide evidence that circular RNA produced from precursor mRNA may have a regulatory role in human cells.

Endogenous microRNA sponges: evidence and controversy

The competitive endogenous RNA (ceRNA) hypothesis proposes that transcripts with shared microRNA (miRNA) binding sites compete for post-transcriptional control. This hypothesis has gained substantial attention as a unifying function for long non-coding RNAs, pseudogene transcripts and circular RNAs, as well as an alternative function for messenger RNAs. Empirical evidence supporting the hypothesis is accumulating but not without attracting scepticism. Recent studies that model transcriptome-wide binding-site abundance suggest that physiological changes in expression of most individual transcripts will not compromise miRNA activity. In this Review, we critically evaluate the evidence for and against the ceRNA hypothesis to assess the impact of endogenous miRNA-sponge interactions.

Circular RNA Expression: Its Potential Regulation and Function

In 2012, a new feature of eukaryotic gene expression emerged: ubiquitous expression of circular RNA (circRNA) from genes traditionally thought to express messenger or linear noncoding (nc)RNA only. CircRNAs are covalently closed, circular RNA molecules that typically comprise exonic sequences and are spliced at canonical splice sites. This feature of gene expression was first recognized in humans and mouse, but it quickly emerged that it was common across essentially all eukaryotes studied by molecular biologists. CircRNA abundance, and even which alternatively spliced circRNA isoforms are expressed, varies by cell type and can exceed the abundance of the traditional linear mRNA or ncRNA transcript. CircRNAs are enriched in the brain and increase in abundance during fetal development. Together, these features raise fundamental questions regarding the regulation of circRNA in cis and in trans, and its function.

The Emerging Function and Mechanism of ceRNAs in Cancer

Complex diseases, such as cancer, are often associated with aberrant gene expression at both the transcriptional and post-transcriptional level. Over the past several years, competing endogenous RNAs (ceRNAs) have emerged as an important class of post-transcriptional regulators that alter gene expression through a miRNA-mediated mechanism. Recent studies in both solid tumors and hematopoietic malignancies showed that ceRNAs have significant roles in cancer pathogenesis by altering the expression of key tumorigenic or tumor-suppressive genes. Characterizing the identity, function, and mechanism of the ceRNAs will not only further our fundamental understanding of RNA-mediated cancer pathogenesis, but may also shed light on the development of new RNA-based therapeutic strategies for treating cancer.

Circular RNAs as a new field in gene regulation and their implications in translational research

Circular RNAs are a class of long noncoding RNA that were recently rediscovered as diverse, highly abundant, conserved and naturally occurring RNAs in eukaryotes. They are characterized by their 5' and 3' covalently joined ends. Some studies have attributed functions for circular RNAs, such as miRNAs sponges and transcriptional regulators, indicating that they may be largely biomarkers of both physiological and pathological processes. Circular RNAs have the potential to play important roles in transcription and post-transcription, giving rise to a whole complexity level to gene expression regulation. In this review, we discuss the biogenesis of circular RNAs, their properties and functions as well as different methods for their identification and their role in some diseases.

A Method for Expressing and Imaging Abundant, Stable, Circular RNAs In Vivo Using tRNA Splicing

Recent improvements in high-throughput sequencing technologies underscore the pervasiveness of circular RNA (circRNA) expression in animal cells. CircRNAs are distinct from their linear counterparts because they lack the 5' caps and 3' tails that typically help determine the cellular fate of a transcript. However, due to the lack of free ends, circRNAs are impervious to exonucleases and thus can evade normal RNA turnover mechanisms. Most circRNAs are derived from protein-coding pre-mRNAs, via a mechanism called "back-splicing." Existing methods of circRNA expression thus typically involve genes that have been engineered to contain sequence elements that promote back-splicing. We recently uncovered an anciently conserved mechanism of RNA circularization in metazoans that involves splicing of tRNA introns. This splicing mechanism is completely independent from that of pre-mRNAs. In this chapter, we detail an orthogonal method that involves splicing of intron-containing tRNAs in order to produce circRNAs in vivo. We utilize fluorescence-based RNA reporters to characterize the expression, localization, and stability of these so-called tRNA intronic circular RNAs. Because tRNA biogenesis is essential for all cellular life, this method provides a means to express ultrastable, high-copy, circRNA effectors in a wide variety of metazoan cell types.

Emerging functions of the Quaking RNA-binding proteins and link to human diseases

RNA-binding proteins (RBPs) are essential players in RNA metabolism including key cellular processes from pre-mRNA splicing to mRNA translation. The K homology-type QUAKING RBP is emerging as a vital factor for oligodendrocytes, monocytes/macrophages, endothelial cell, and myocyte function. Interestingly, the qkI gene has now been identified as the culprit gene for a patient with intellectual disabilities and is translocated in a pediatric ganglioglioma as a fusion protein with MYB. In this review, we will focus on the emerging discoveries of the QKI proteins as well as highlight the recent advances in understanding the role of QKI in human disease pathology including myelin disorders, schizophrenia and cancer. WIREs RNA 2016, 7:399-412. doi: 10.1002/wrna.1344 For further resources related to this article, please visit the WIREs website.

Circular RNA Related to the Chondrocyte ECM Regulates MMP13 Expression by Functioning as a MiR-136 'Sponge' in Human Cartilage Degradation

Circular RNAs (circRNAs) are involved in the development of various diseases, but there is little knowledge of circRNAs in osteoarthritis (OA). The aim of study was to identify circRNA expression in articular cartilage and to explore the function of chondrocyte extracellular matrix (ECM)-related circRNAs (circRNA-CER) in cartilage. To identify circRNAs that are specifically expressed in cartilage, we compared the expression of circRNAs in OA cartilage with that in normal cartilage. Bioinformatics was employed to predict the interaction of circRNAs and mRNAs in cartilage. Loss-of-function and rescue experiments for circRNA-CER were performed in vitro. A total of 71 circRNAs were differentially expressed in OA and normal cartilage. CircRNA-CER expression increased with interleukin-1 and tumor necrosis factor levels in chondrocytes. Silencing of circRNA-CER using small interfering RNA suppressed MMP13 expression and increased ECM formation. CircRNA-CER could compete for miR-136 with MMP13. Our results demonstrated that circRNA-CER regulated MMP13 expression by functioning as a competing endogenous RNA (ceRNA) and participated in the process of chondrocyte ECM degradation. We propose that circRNA-CER could be used as a potential target in OA therapy.

Circular RNA ITCH has inhibitory effect on ESCC by suppressing the Wnt/β-catenin pathway

Circular RNAs with exonic sequences represent a special form of non-coding RNAs, discovered by analyzing a handful of transcribed genes. It has been observed that circular RNAs function as microRNA sponges. In the present study, we investigated whether the expression of circular RNAs is altered during the development of esophageal squamous cell carcinoma (ESCC). Using a TaqMan-based reverse transcriptase polymerase chain reaction assay, the relationship between cir-ITCH and ESCC was analyzed in a total of 684 ESCC and paired adjacent non-tumor tissue samples from eastern and southern China. We found that cir-ITCH expression was usually low in ESCC compared to the peritumoral tissue. The functional relevance of cir-ITCH was further examined by biochemical assays. As sponge of miR-7, miR-17, and miR-214, cir-ITCH might increase the level of ITCH. ITCH hyper expression promotes ubiquitination and degradation of phosphorylated Dvl2, thereby inhibiting the Wnt/β-catenin pathway. These results indicate that cir-ITCH may have an inhibitory effect on ESCC by regulating the Wnt pathway.

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.

The biogenesis and emerging roles of circular RNAs

Circular RNAs (circRNAs) are produced from precursor mRNA (pre-mRNA) back-splicing of thousands of genes in eukaryotes. Although circRNAs are generally expressed at low levels, recent findings have shed new light on their cell type-specific and tissue-specific expression and on the regulation of their biogenesis. Furthermore, the data indicate that circRNAs shape gene expression by titrating microRNAs, regulating transcription and interfering with splicing, thus effectively expanding the diversity and complexity of eukaryotic transcriptomes.

Noncoding RNAs: Regulators of the Mammalian Transcription Machinery

Transcription by RNA polymerase II (Pol II) is required to produce mRNAs and some noncoding RNAs (ncRNAs) within mammalian cells. This coordinated process is precisely regulated by multiple factors, including many recently discovered ncRNAs. In this perspective, we will discuss newly identified ncRNAs that facilitate DNA looping, regulate transcription factor binding, mediate promoter-proximal pausing of Pol II, and/or interact with Pol II to modulate transcription. Moreover, we will discuss new roles for ncRNAs, as well as a novel Pol II RNA-dependent RNA polymerase activity that regulates an ncRNA inhibitor of transcription. As the multifaceted nature of ncRNAs continues to be revealed, we believe that many more ncRNA species and functions will be discovered.

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.

Circular RNA: a novel biomarker for progressive laryngeal cancer

Circular RNAs (circRNAs), a class of endogenous RNAs, are characterized by covalently closed continuous loop without 5' to 3' polarity and polyadenylated tail. Recent studies indicated that circRNAs might play an important role in cancer. However, the function of circRNA in laryngeal squamous cell cancer tissues (LSCC) is still unknown. In this study, we investigated the expression of circRNAs in 4 paired LSCC tissues and adjacent non-tumor tissues by microarray analysis. Results showed significant upregulation (n = 302) of or downregulation (n = 396) of 698 circRNAs in LSCC tissues. We further detected hsa_circRNA_100855 as the most upregulated circRNA and hsa_circRNA_104912 as the most downregulated circRNA using qRT-PCR methods. Results showed that hsa_circRNA_100855 level was significantly higher in LSCC than in the corresponding adjacent non-neoplastic tissues. Patients with T3-4 stage, neck nodal metastasis or advanced clinical stage had higher hsa_circRNA_100855 expression. The hsa_circRNA_104912 level was significantly lower in LSCC than in corresponding adjacent non-neoplastic tissues. Patients with T3-4 stage, neck nodal metastasis, poor differentiation or advanced clinical stage had a lower hsa_circRNA_104912 expression. Overall, our data suggest that circRNAs play an important role in the tumorigenesis of LSCC and may serve as novel and stable biomarkers for the diagnosis and progress of LSCC.

Circular RNA profile in gliomas revealed by identification tool UROBORUS

Recent evidence suggests that many endogenous circular RNAs (circRNAs) may play roles in biological processes. However, the expression patterns and functions of circRNAs in human diseases are not well understood. Computationally identifying circRNAs from total RNA-seq data is a primary step in studying their expression pattern and biological roles. In this work, we have developed a computational pipeline named UROBORUS to detect circRNAs in total RNA-seq data. By applying UROBORUS to RNA-seq data from 46 gliomas and normal brain samples, we detected thousands of circRNAs supported by at least two read counts, followed by successful experimental validation on 24 circRNAs from the randomly selected 27 circRNAs. UROBORUS is an efficient tool that can detect circRNAs with low expression levels in total RNA-seq without RNase R treatment. The circRNAs expression profiling revealed more than 476 circular RNAs differentially expressed in control brain tissues and gliomas. Together with parental gene expression, we found that circRNA and its parental gene have diversified expression patterns in gliomas and control brain tissues. This study establishes an efficient and sensitive approach for predicting circRNAs using total RNA-seq data. The UROBORUS pipeline can be accessed freely for non-commercial purposes at http://uroborus.openbioinformatics.org/.

Retrotransposons as regulators of gene expression

Transposable elements (TEs) are both a boon and a bane to eukaryotic organisms, depending on where they integrate into the genome and how their sequences function once integrated. We focus on two types of TEs: long interspersed elements (LINEs) and short interspersed elements (SINEs). LINEs and SINEs are retrotransposons; that is, they transpose via an RNA intermediate. We discuss how LINEs and SINEs have expanded in eukaryotic genomes and contribute to genome evolution. An emerging body of evidence indicates that LINEs and SINEs function to regulate gene expression by affecting chromatin structure, gene transcription, pre-mRNA processing, or aspects of mRNA metabolism. We also describe how adenosine-to-inosine editing influences SINE function and how ongoing retrotransposition is countered by the body's defense mechanisms.

Changing expression profiles of lncRNAs, mRNAs, circRNAs and miRNAs during osteoclastogenesis

Bone is a dynamic organ continuously undergoing shaping, repairing and remodeling. The homeostasis of bone is maintained by the balance between osteoblastic bone formation and osteoclastic bone resorption. Osteoclasts (OCs) are specialized multinucleated cells derived from hematopoietic stem cells (HSCs) or monocytes/macrophage progenitor cells. There are different stages during osteoclastogenesis, and one of the most important steps to form functional osteoclasts is realized by cell-cell fusion. In our study, microarray was performed to detect the expression profiles of lncRNA, mRNA, circRNA and miRNA at different stages during osteoclastogenesis of RAW264.7 cells. Often changed RNAs were selected and clustered among the four groups with Venn analysis. The results revealed that expressions of 518 lncRNAs, 207 mRNAs, 24 circRNAs and 37 miRNAs were often altered at each stage during OC differentiation. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) biological pathway analysis were performed to predict the functions of differentially expressed lncRNAs and co-expressed potential targeting genes. Co-expression networks of lncRNA-mRNA and circRNA-miRNA were constructed based on the correlation analysis between the differentially expressed RNAs. The present study provided a systematic perspective on the potential function of non-coding RNAs (ncRNAs) during osteoclastogenesis.

Circular RNAs Co-Precipitate with Extracellular Vesicles: A Possible Mechanism for circRNA Clearance

Backspliced circular RNAs (circRNAs) are prevalent in many eukaryotic systems and are spliced from thousands of different genes. Where examined, circRNAs are often highly stable and the mechanisms by which cells degrade and/or clear circRNAs from the cells are unknown. Here we investigated the possibility that cells can eliminate circRNAs into extracellular space, possibly within released vesicles such as exosomes and microvesicles. From three different cell lines and examining multiple circRNAs, we show that extracellular vesicle (EVs) preparations recovered from cell culture conditioned media contain established circRNAs. Moreover, these circRNAs are enriched over their linear counterparts within EV preparations when compared to the producing cells. This supports the idea that expulsion from cells into extracellular space, as by EVs release, can be a mechanism by which cells clear circRNAs. Moreover, since EVs can be taken up by other cells, excreted circRNAs could contribute to cell to cell communication.

Computational identification of circular RNAs based on conformational and thermodynamic properties in the flanking introns

Circular RNAs (circRNAs) were found more than 30 years ago, but have been treated as molecular flukes in a long time. Combining deep sequencing studies with bioinformatics technique, thousands of endogenous circRNAs have been found in mammalian cells, and some researchers have proved that several circRNAs act as competing endogenous RNAs (ceRNAs) to regulate gene expression. However, the mechanism by which the precursor mRNA to be transformed into a circular RNA or a linear mRNA is largely unknown. In this paper, we attempted to bioinformatically identify shared genomic features that might further elucidate the mechanism of formation and proposed a SVM-based model to distinguish circRNAs from non-circularized, expressed exons. Firstly, conformational and thermodynamic dinucleotide properties in the flanking introns were extracted as potential features. Secondly, two feature selection methods were applied to gain the optimal feature subset. Our 10-fold cross-validation results showed that the model can be used to distinguish circRNAs from non-circularized, expressed exons with an Sn of 0.884, Sp of 0.900, ACC of 0.892, MCC of 0.784, respectively. The identification results suggest that conformational and thermodynamic properties in the flanking introns are closely related to the formation of circRNAs. Datasets and the tool involved in this paper are all available at https://sourceforge.net/projects/predicircrnatool/files/.

Circular RNA expression alterations are involved in OGD/R-induced neuron injury

Cerebral ischemia-reperfusion injury (IRI) is a common clinical pathological process, and it is a key step in causing further ischemic organ damage. The mechanism of cerebral IRI is still not fully understood, leading to a lack of effective treatment. It has been demonstrated that circular RNAs (circRNAs) can act as miRNA sponges and play an important role in regulating gene expression through a circRNA-miRNA-gene pathway. The specific role of circRNAs in the pathogenesis of cerebral IRI, however, is still unclear. Thus, in the present study, we investigated circRNA expression differences in HT22 cells with oxygen-glucose deprivation/reoxygenation (OGD/R) versus normal controls. The results from circRNA microarrays revealed that 15 circRNAs were significantly altered in the OGD/R model (p < 0.05) compared with the control group. Among them, 3 were significantly up-regulated, and the other 12 were down-regulated. Furthermore, the up-regulated expression of mmu-circRNA-015947 was verified using quantitative real-time polymerase chain reaction (qRT-PCR). Bioinformatics analysis revealed that up-regulated expression of mmu-circRNA-015947 could interact with miRNAs (mmu-miR-188-3p, mmu-miR-329-5p, mmu-miR-3057-3p, mmu-miR-5098 and mmu-miR-683) and thereby enhance target gene expression. KEGG pathway analysis predicted that mmu-circRNA-015947 may participate in apoptosis-related, metabolism-related and immune-related pathways, which are known to be involved in the pathogenesis of IRI. This research suggests that the overlapping expression of mmu-circRNA-015947 might be involved in the process of cerebral IRI and presents a novel molecular target for clinical therapy.

Sequencing of lariat termini in S. cerevisiae reveals 5' splice sites, branch points, and novel splicing events

Pre-mRNA splicing is a central step in the shaping of the eukaryotic transcriptome and in the regulation of gene expression. Yet, due to a focus on fully processed mRNA, common approaches for defining pre-mRNA splicing genome-wide are suboptimal-especially with respect to defining the branch point sequence, a key cis-element that initiates the chemistry of splicing. Here, we report a complementary intron-centered approach designed to more efficiently, simply, and directly define splicing events genome-wide. Specifically, we developed a method distinguished by deep sequencing of lariat intron termini (LIT-seq). In a test of LIT-seq using the budding yeast Saccharomyces cerevisiae, we not only successfully captured the majority of annotated, expressed splicing events but also uncovered 45 novel splicing events, establishing the sensitivity of LIT-seq. Moreover, our libraries were highly enriched with reads that reported on splice sites; by a simple and direct inspection of sequencing reads, we empirically defined both 5' splice sites and branch sites, as well as their consensus sequences, with nucleotide resolution. Additionally, our study revealed that the 3' termini of lariat introns are subject to nontemplated addition of adenosines, characteristic of signals sensed by 3' to 5' RNA turnover machinery. Collectively, this work defines a novel, genome-wide approach for analyzing splicing with unprecedented depth, specificity, and resolution.

A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223

AIMS

Sustained cardiac hypertrophy accompanied by maladaptive cardiac remodelling represents an early event in the clinical course leading to heart failure. Maladaptive hypertrophy is considered to be a therapeutic target for heart failure. However, the molecular mechanisms that regulate cardiac hypertrophy are largely unknown.

METHODS AND RESULTS

Here we show that a circular RNA (circRNA), which we term heart-related circRNA (HRCR), acts as an endogenous miR-223 sponge to inhibit cardiac hypertrophy and heart failure. miR-223 transgenic mice developed cardiac hypertrophy and heart failure, whereas miR-223-deficient mice were protected from hypertrophic stimuli, indicating that miR-223 acts as a positive regulator of cardiac hypertrophy. We identified ARC as a miR-223 downstream target to mediate the function of miR-223 in cardiac hypertrophy. Apoptosis repressor with CARD domain transgenic mice showed reduced hypertrophic responses. Further, we found that a circRNA HRCR functions as an endogenous miR-223 sponge to sequester and inhibit miR-223 activity, which resulted in the increase of ARC expression. Heart-related circRNA directly bound to miR-223 in cytoplasm and enforced expression of HRCR in cardiomyocytes and in mice both exhibited attenuated hypertrophic responses.

CONCLUSIONS

These findings disclose a novel regulatory pathway that is composed of HRCR, miR-223, and ARC. Modulation of their levels provides an attractive therapeutic target for the treatment of cardiac hypertrophy and heart failure.

RNA Editing: A Contributor to Neuronal Dynamics in the Mammalian Brain

Post-transcriptional RNA modification by adenosine to inosine (A-to-I) editing expands the functional output of many important neuronally expressed genes. The mechanism provides flexibility in the proteome by expanding the variety of isoforms, and is a requisite for neuronal function. Indeed, targets for editing include key mediators of synaptic transmission with an overall significant effect on neuronal signaling. In addition, editing influences splice-site choice and miRNA targeting capacity, and thereby regulates neuronal gene expression. Editing efficiency at most of these sites increases during neuronal differentiation and brain maturation in a spatiotemporal manner. This editing-induced dynamics in the transcriptome is essential for normal brain development, and we are only beginning to understand its role in neuronal function. In this review we discuss the impact of RNA editing in the brain, with special emphasis on the physiological consequences for neuronal development and plasticity.

PTESFinder: a computational method to identify post-transcriptional exon shuffling (PTES) events

BACKGROUND

Transcripts, which have been subject to Post-transcriptional exon shuffling (PTES), have an exon order inconsistent with the underlying genomic sequence. These have been identified in a wide variety of tissues and cell types from many eukaryotes, and are now known to be mostly circular, cytoplasmic, and non-coding. Although there is no uniformly ascribed function, several have been shown to be involved in gene regulation. Accurate identification of these transcripts can, however, be difficult due to artefacts from a wide variety of sources.

RESULTS

Here, we present a computational method, PTESFinder, to identify these transcripts from high throughput RNAseq data. Uniquely, it systematically excludes potential artefacts emanating from pseudogenes, segmental duplications, and template switching, and outputs both PTES and canonical exon junction counts to facilitate comparative analyses. In comparison with four existing methods, PTESFinder achieves highest specificity and comparable sensitivity at a variety of read depths. PTESFinder also identifies between 13 % and 41.6 % more structures, compared to publicly available methods recently used to identify human circular RNAs.

CONCLUSIONS

With high sensitivity and specificity, user-adjustable filters that target known sources of false positives, and tailored output to facilitate comparison of transcript levels, PTESFinder will facilitate the discovery and analysis of these poorly understood transcripts.

Regulation of circRNA biogenesis

Unlike linear RNAs terminated with 5' caps and 3' tails, circular RNAs are characterized by covalently closed loop structures with neither 5' to 3' polarity nor polyadenylated tail. This intrinsic characteristic has led to the general under-estimation of the existence of circular RNAs in previous polyadenylated transcriptome analyses. With the advent of specific biochemical and computational approaches, a large number of circular RNAs from back-spliced exons (circRNAs) have been identified in various cell lines and across different species. Recent studies have uncovered that back-splicing requires canonical spliceosomal machinery and can be facilitated by both complementary sequences and specific protein factors. In this review, we highlight our current understanding of the regulation of circRNA biogenesis, including both the competition between splicing and back-splicing and the previously under-appreciated alternative circularization.

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.

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.

Deciphering the cross-talking of human competitive endogenous RNAs in K562 chronic myelogenous leukemia cell line

Chronic myelogenous leukemia (CML) is a myeloproliferative disorder characterized by increased proliferation or abnormal accumulation of the granulocytic cell line without the depletion of their capacity to differentiate.

RNA structure replaces the need for U2AF2 in splicing

RNA secondary structure plays an integral role in catalytic, ribosomal, small nuclear, micro, and transfer RNAs. Discovering a prevalent role for secondary structure in pre-mRNAs has proven more elusive. By utilizing a variety of computational and biochemical approaches, we present evidence for a class of nuclear introns that relies upon secondary structure for correct splicing. These introns are defined by simple repeat expansions of complementary AC and GT dimers that co-occur at opposite boundaries of an intron to form a bridging structure that enforces correct splice site pairing. Remarkably, this class of introns does not require U2AF2, a core component of the spliceosome, for its processing. Phylogenetic analysis suggests that this mechanism was present in the ancestral vertebrate lineage prior to the divergence of tetrapods from teleosts. While largely lost from land dwelling vertebrates, this class of introns is found in 10% of all zebrafish genes.

Long noncoding RNAs: Re-writing dogmas of RNA processing and stability

Most of the human genome is transcribed, yielding a complex network of transcripts that includes tens of thousands of long noncoding RNAs. Many of these transcripts have a 5' cap and a poly(A) tail, yet some of the most abundant long noncoding RNAs are processed in unexpected ways and lack these canonical structures. Here, I highlight the mechanisms by which several of these well-characterized noncoding RNAs are generated, stabilized, and function. The MALAT1 and MEN β (NEAT1_2) long noncoding RNAs each accumulate to high levels in the nucleus, where they play critical roles in cancer progression and the formation of nuclear paraspeckles, respectively. Nevertheless, MALAT1 and MEN β are not polyadenylated as the tRNA biogenesis machinery generates their mature 3' ends. In place of a poly(A) tail, these transcripts are stabilized by highly conserved triple helical structures. Sno-lncRNAs likewise lack poly(A) tails and instead have snoRNA structures at their 5' and 3' ends. Recent work has additionally identified a number of abundant circular RNAs generated by the pre-mRNA splicing machinery that are resistant to degradation by exonucleases. As these various transcripts use non-canonical strategies to ensure their stability, it is becoming increasingly clear that long noncoding RNAs may often be regulated by unique post-transcriptional control mechanisms. This article is part of a Special Issue entitled: Clues to long noncoding RNA taxonomy1, edited by Dr. Tetsuro Hirose and Dr. Shinichi Nakagawa.