CIRCpedia v2: An Updated Database for Comprehensive Circular RNA Annotation and Expression Comparison

Circular RNAs (circRNAs) from back-splicing of exon(s) have been recently identified to be broadly expressed in eukaryotes, in tissue- and species-specific manners. Although functions of most circRNAs remain elusive, some circRNAs are shown to be functional in gene expression regulation and potentially relate to diseases. Due to their stability, circRNAs can also be used as biomarkers for diagnosis. Profiling circRNAs by integrating their expression among different samples thus provides molecular basis for further functional study of circRNAs and their potential application in clinic. Here, we report CIRCpedia v2, an updated database for comprehensive circRNA annotation from over 180 RNA-seq datasets across six different species. This atlas allows users to search, browse, and download circRNAs with expression features in various cell types/tissues, including disease samples. In addition, the updated database incorporates conservation analysis of circRNAs between humans and mice. Finally, the web interface also contains computational tools to compare circRNA expression among samples. CIRCpedia v2 is accessible at http://www.picb.ac.cn/rnomics/circpedia.     

CircPlant:An Integrated Tool for circRNA Detection and Functional Prediction in Plants

The recent discovery of circular RNAs (circRNAs) and characterization of their functional roles have opened a new avenue for understanding the biology of genomes. circRNAs have been implicated to play important roles in a variety of biological processes, but their precise functions remain largely elusive. Currently, a few approaches are available for novel circRNA prediction, but almost all these methods are intended for animal genomes. Considering that the major differences between the organization of plant and mammal genomes cannot be neglected, a plant-specific method is needed to enhance the validity of plant circRNA identification. In this study, we present CircPlant, an integrated tool for the exploration of plant circRNAs, potentially acting as competing endogenous RNAs (ceRNAs), and their potential functions. With the incorporation of several unique plant-specific criteria, CircPlant can accurately detect plant circRNAs from high-throughput RNA-seq data. Based on comparison tests on simulated and real RNA-seq datasets from Arabidopsis thaliana and Oryza sativa, we show that CircPlant outperforms all evaluated competing tools in both accuracy and efficiency. CircPlant is freely available at http://bis.zju.edu.cn/circplant.     

Identification of potential proteins translated from circular RNA splice variants

Circular RNAs (circRNAs) are covalently closed RNA molecules generated from precursor RNAs by the head-to-tail backsplicing of exons. Hundreds of studies demonstrated that circRNAs are ubiquitously expressed and regulate cellular events by modulating microRNA (miRNA) and RNA-binding protein (RBP) activities. A few circRNAs are also known to translate into functional polypeptides regulating cellular physiology. All these functions primarily depend on the full-length sequence of the circRNAs. CircRNA backsplice junction sequence is the key to identifying circRNAs and their full-length mature sequence. However, some multi-exonic circRNAs exist in different isoforms sharing identical backsplice junction sequences and are termed circRNA splice variants. Here, we analyzed the previously published HeLa cell RNA-seq datasets to identify circRNA splice variants using the de novo module of the CIRCexplorer2 circRNA annotation pipeline. A subset of circRNAs with splice variants was validated by the circRNA-rolling circle amplification (circRNA-RCA) method. Interestingly, several validated circRNAs were predicted to translate into proteins by the riboCIRC database. Furthermore, polyribosome fractionation followed by quantitative PCR confirmed the association of a subset of circRNAs with polyribosome supporting their protein-coding potential. Finally, bioinformatics analysis of proteins derived from splice variants of circCORO1C and circASPH suggested altered protein sequences and structures that could affect their physiological functions. Together, our study identified novel circRNA splice variants and their potential translation into protein isoforms which may regulate various physiological processes.     

Circular RNA expression profiles and features in human tissues: a study using RNA-seq data

Background

Circular RNA (circRNA) is one type of noncoding RNA that forms a covalently closed continuous loop. Similar to long noncoding RNA (lncRNA), circRNA can act as microRNA (miRNA) ‘sponges’ to regulate gene expression, and its abnormal expression is related to diseases such as atherosclerosis, nervous system disorders and cancer. So far, there have been no systematic studies on circRNA abundance and expression profiles in human adult and fetal tissues.

Results

We explored circRNA expression profiles using RNA-seq data for six adult and fetal normal tissues (colon, heart, kidney, liver, lung, and stomach) and four gland normal tissues (adrenal gland, mammary gland, pancreas, and thyroid gland). A total of 8120, 25,933 and 14,433 circRNAs were detected by at least two supporting junction reads in adult, fetal and gland tissues, respectively. Among them, 3092, 14,241 and 6879 circRNAs were novel when compared to the published results. In each adult tissue type, we found at least 1000 circRNAs, among which 36.97–50.04% were tissue-specific. We reported 33 circRNAs that were ubiquitously expressed in all the adult tissues we examined. To further explore the potential “housekeeping” function of these circRNAs, we constructed a circRNA-miRNA-mRNA regulatory network containing 17 circRNAs, 22 miRNAs and 90 mRNAs. Furthermore, we found that both the abundance and the relative expression level of circRNAs were higher in fetal tissue than adult tissue. The number of circRNAs in gland tissues, especially in mammary gland (9665 circRNA candidates), was higher than that of other adult tissues (1160–3777).

Conclusions

We systematically investigated circRNA expression in a variety of human adult and fetal tissues. Our observation of different expression level of circRNAs in adult and fetal tissues suggested that circRNAs might play their role in a tissue-specific and development-specific fashion. Analysis of circRNA-miRNA-mRNA network provided potential targets of circRNAs. High expression level of circRNAs in mammary gland might be attributed to the rich innervation.

     

NOVA2 regulates neural circRNA biogenesis

Circular RNAs (circRNAs) are highly expressed in the brain and their expression increases during neuronal differentiation. The factors regulating circRNAs in the developing mouse brain are unknown. NOVA1 and NOVA2 are neural-enriched RNA-binding proteins with well-characterized roles in alternative splicing. Profiling of circRNAs from RNA-seq data revealed that global circRNA levels were reduced in embryonic cortex of Nova2 but not Nova1 knockout mice. Analysis of isolated inhibitory and excitatory cortical neurons lacking NOVA2 revealed an even more dramatic reduction of circRNAs and establishes a widespread role for NOVA2 in enhancing circRNA biogenesis. To investigate the cis-elements controlling NOVA2-regulation of circRNA biogenesis, we generated a backsplicing reporter based on the Efnb2 gene. We found that NOVA2-mediated backsplicing of circEfnb2 was impaired when YCAY clusters located in flanking introns were mutagenized. CLIP (cross-linking and immunoprecipitation) and additional reporter analyses demonstrated the importance of NOVA2 binding sites located in both flanking introns of circRNA loci. NOVA2 is the first RNA-binding protein identified to globally promote circRNA biogenesis in the developing brain.     

IRESbase:A Comprehensive Database of Experimentally Validated Internal Ribosome Entry Sites

Internal ribosome entry sites (IRESs) are functional RNA elements that can directly recruit ribosomes to an internal position of the mRNA in a cap-independent manner to initiate translation. Recently, IRES elements have attracted much attention for their critical roles in various processes including translation initiation of a new type of RNA, circular RNA (circRNA), with no 5′ cap to support classical cap-dependent translation. Thus, an integrative data resource of IRES elements with experimental evidence will be useful for further studies. In this study, we present IRESbase, a comprehensive database of IRESs, by curating the experimentally validated functional minimal IRES elements from literature and annotating their host linear and circular RNAs. The current version of IRESbase contains 1328 IRESs, including 774 eukaryotic IRESs and 554 viral IRESs from 11 eukaryotic organisms and 198 viruses, respectively. As IRESbase collects only IRES of minimal length with functional evidence, the median length of IRESs in IRESbase is 174 nucleotides. By mapping IRESs to human circRNAs and long non-coding RNAs (lncRNAs), 2191 circRNAs and 168 lncRNAs were found to contain at least one entire or partial IRES sequence. IRESbase is available at http://reprod.njmu.edu.cn/cgi-bin/iresbase/index.php.     

Characteristics of circular RNA expression of pulmonary macrophages in mice with sepsis‐induced acute lung injury

Circular RNAs (circRNAs) make up a large class of non‐coding RNAs and play important roles in the pathology of a variety of diseases. However, their roles in pulmonary macrophage polarization after sepsisinduced lung injury is unknown. In this study, mice were divided into two groups: Sham control group and cecal ligation and puncture (CLP)‐induced ALI group. Macrophages were isolated from lung homogenates 24 hours after SCLP/CLP. We started with RNA‐seq of circRNA changes in macrophages and validated by RT‐PCR in the following experiments. A total of 4318 circRNAs were detected in the two groups. Of these, 11 and 126 circRNAs were found to be significantly upregulated and downregulated, respectively, compared to the control (p≤0.05, Fold Change ≥2). Differentially expressed circRNAs with a high foldchange (fold‐change >4, P<0.05) were selected for validation by qRT‐PCR, 10 of which were verified. Furthermore, the most differentially expressed circRNAs within all the comparisons were annotated in detail with circRNA/miRNA interaction information using miRNA target prediction software. The network of circRNA‐miRNA‐mRNA was illustrated by cytoscape software. Gene ontology analyses indicated the upregulated circRNAs were involved in the multiple biological functions such as regulation of mitochondrion distribution and Notch binding, while the down‐regulated circRNAs mainly involved in the biological process as histone H3K27 methylation. KEGG pathway analysis revealed TGF‐beta signaling pathway was related to the upregulated circRNAs. The present study provides a novel insight into the roles of circRNAs in pulmonary macrophage differentiation and polarization post septic lung injury.     

条纹斑竹鲨肝脏环形RNA的筛选、鉴定及功能分析北大核心CSCD

为筛选和验证条纹斑竹鲨肝脏中环形RNA(circRNA)及探究其在肝癌细胞HepG2中过表达对肝癌细胞增殖、迁移能力的影响,本研究主要进行了两项实验:对条纹斑竹鲨肝脏circRNA进行高通量测序和预测,随后设计正、反向引物验证其真实性;构建circRNA过表达载体,将其瞬时转染进肝癌细胞HepG2,进行CCK-8和划痕实验来评价其对肝癌细胞增殖和迁移能力的影响。结果显示:预测到有4558条circRNAs,并确认了14条circRNAs的真实性;qRT-PCR实验表明在肝癌细胞HepG2中能瞬时过表达circRNA 13-566、circRNA 4-475、circRNA 5-402、circRNA 294-177、circRNA 30-219;且CCK-8和划痕实验显示,这5条circRNAs过表达后,均能不同程度地抑制肝癌细胞的增殖和迁移能力,其中circRNA 4-475、circRNA 294-177作用尤为显著。上述结果为深入研究条纹斑竹鲨肝脏中circRNA及其在肝再生、肝癌治疗方面的功能提供了新思路和基础。     

意大利蜜蜂工蜂中肠的环状RNA及其调控网络分析

【目的】环状RNA(circRNA)在可变剪接、转录调控和来源基因的表达调控等方面具有重要功能。本研究旨在分析意大利蜜蜂Apis mellifera ligustica工蜂中肠circRNA的数量、种类、结构特征和作用,并通过构建和分析circRNA的调控网络探索circRNA的调控功能。【方法】在实验室条件下人工饲养意大利蜜蜂工蜂,利用circRNA-seq技术对意大利蜜蜂7和10日龄成年工蜂中肠样品进行深度测序。利用find_circ软件从质控后的数据中预测circRNA。通过BLAST比对GO和KEGG数据库,对circRNA的来源基因进行功能和代谢通路注释。利用TargetFinder软件预测circRNA靶向结合的miRNA及miRNA靶向结合的mRNA,通过Cytoscape v.3.2.1软件对circRNA-miRNA和circRNA-miRNA-mRNA调控网络进行构建及可视化。通过设计背靠背引物和线性扩增引物RT-PCR对预测出的circRNA进行验证。【结果】意大利蜜蜂工蜂中肠样品的测序平均得到136 463 071条clean reads,去除rRNA后各样品的anchor reads均在136 779 122条及以上。共预测出10 833个circRNA,长度主要介于15~1 000 nt;上述circRNA的类型丰富,其中已注释的外显子circRNA数量最多,分布在西方蜜蜂1号染色体的circRNA数量最多,其次为8号染色体。CircRNA的来源基因可注释到包括结合、细胞进程和细胞在内的45个GO条目,以及包括内吞作用、内质网蛋白加工及核糖体在内的121条KEGG代谢通路,表明circRNA在意大利蜜蜂工蜂中肠的生长、发育、新陈代谢和细胞生命活动等生物学过程中发挥重要作用。进一步构建circRNA-miRNA和circRNA-miRNA-mRNA调控网络,分析结果显示部分circRNA可能作为竞争性内源RNA吸附结合microRNA,从而调控基因的表达水平。最后,对随机选择的3个circRNA的RT-PCR结果验证了其真实存在。【结论】本研究对意大利蜜蜂工蜂中肠中的circRNA进行预测、分析及鉴定。研究结果提供了中肠circRNA的数量、种类、结构特征、作用和调控网络的信息,揭示了circRNA可能通过作用于来源基因和作为竞争性内源RNA在意大利蜜蜂工蜂中肠的生长发育和免疫防御中发挥作用,为深入研究circRNA在意大利蜜蜂中肠发育及胁迫响应过程中的功能奠定了基础。     

Differential expression of circular RNAs in bone marrow-derived exosomes from essential thrombocythemia patients

Circular RNAs (circRNA) are closely associated with the pathogenesis of various hematological diseases. However, little is known about the potential functions of circRNAs in essential thrombocythemia (ET) development. The circRNA profile alterations in the bone marrow of ET patients were mainly investigated in this study. The sizes of exosomes derived from human bone marrow tissues were validated by the nanoparticle tracking analysis (NTA) method. CD63 and TSG101 expressions in exosomes were analyzed by western blot analysis. The profiles and differential expression of circRNAs in bone-derived exosomes were characterized by high-throughput sequencing. Herein, circular structures and expression of circRNAs were verified by Sanger sequencing and real-time polymerase chain reaction, respectively. The circRNA-miRNA-mRNA networks were predicted using the Cytoscape software. And we detected the effect of circ_0014614 on the transformation of K562 cells into megakaryocytes. Exosomes derived from the bone marrow of ET patients and healthy volunteers showed a diameter between 70 and 140 nm and expressed high CD63 and TSG101. Meanwhile, the circRNA profiles were significantly altered in bone marrow–derived exosomes from ET patients, among which circDAP3, circASXL1, and circRUNX1 were significantly downregulated in ET patients, thus conferring a new insight into the role of circRNAs in the pathogenesis of ET. Besides this, circRNA-encoding genes and miRNA-mRNA networks targeted by this three circRNA were involved in various biological processes and signaling pathways. And circ_0014614 could inhibit K562 cells' differentiation into megakaryocytes. The predictions of the potential function of these three differentially expressed circRNAs along with their interaction with specific miRNAs could provide a basis for circRNA-based ET diagnosis and treatment.     

CircRNA expression profile and functional analysis in retinal ischemia-reperfusion injury

Retinal ischemia-reperfusion (I/R) is involved in the pathogenesis of many vision-threatening diseases. circRNAs act as key players in gene regulation and human diseases. However, the global circRNA expression profile in retinal I/R injury has not been fully uncovered. Herein, we established a murine model of retinal I/R injury and performed circRNA microarrays to identify I/R-related circRNAs. 1265 differentially expressed circRNAs were identified between I/R retinas and normal retinas. Notably, the detection of cWDR37 level in aqueous humor could discriminate glaucoma patients from cataract patients (AUC = 0.9367). cWdr37 silencing protected against hypoxic stress- or oxidative stress-induced retinal ganglion cell (RGC) injury. cWdr37 silencing alleviated IR-induced retinal neurodegeneration as shown by increased NeuN staining, reduced retinal reactive gliosis, and decreased retinal apoptosis. Collectively, this study provides a novel insight into the pathogenesis of retinal I/R injury. cWdr37 is a promising target for the diagnosis or treatment of I/R-related ocular diseases.     

Circular RNA circDtx1 regulates IRF3-mediated antiviral immune responses through suppression of miR-15a-5p-dependent TRIF downregulation in teleost fish

Circular RNAs (circRNAs) represent a class of widespread and diverse covalently closed circular endogenous RNAs that exert crucial functions in regulating gene expression in mammals. However, the function and regulation mechanism of circRNAs in lower vertebrates are still unknown. Here, we discovered a novel circRNA derived from Deltex E3 ubiquitin ligase 1 (Dtx1) gene, namely, circDtx1, which was related to the antiviral responses in teleost fish. Results indicated that circDtx1 played essential roles in host antiviral immunity and inhibition of SCRV replication. Our study also found a microRNA miR-15a-5p, which could inhibit antiviral immune response and promote viral replication by targeting TRIF. Moreover, we also found that the antiviral effect inhibited by miR-15a-5p could be reversed with the circDtx1. In mechanism, our data revealed that circDtx1 was a competing endogenous RNA (ceRNA) of TRIF by sponging miR-15a-5p, leading to activation of the NF-κB/IRF3 pathway, and then enhancing the innate antiviral responses. Our results indicated that circRNAs played a regulatory role in immune responses in teleost fish.