Large-Scale Profiling of RBP-circRNA Interactions from Public CLIP-Seq Datasets

Computational approaches and challenges in the analysis of circRNA data

Circular RNAs (circRNA) are a class of non-coding RNA, forming a single-stranded covalently closed loop structure generated via back-splicing. Advancements in sequencing methods and technologies in conjunction with algorithmic developments of bioinformatics tools have enabled researchers to characterise the origin and function of circRNAs, with practical applications as a biomarker of diseases becoming increasingly relevant. Computational methods developed for circRNA analysis are predicated on detecting the chimeric back-splice junction of circRNAs whilst mitigating false-positive sequencing artefacts. In this review, we discuss in detail the computational strategies developed for circRNA identification, highlighting a selection of tool strengths, weaknesses and assumptions. In addition to circRNA identification tools, we describe methods for characterising the role of circRNAs within the competing endogenous RNA (ceRNA) network, their interactions with RNA-binding proteins, and publicly available databases for rich circRNA annotation.

FL-circAS: an integrative resource and analysis for full-length sequences and alternative splicing of circular RNAs with nanopore sequencing

Circular RNAs (circRNAs) are RNA molecules with a continuous loop structure characterized by back-splice junctions (BSJs). While analyses of short-read RNA sequencing have identified millions of BSJ events, it is inherently challenging to determine exact full-length sequences and alternatively spliced (AS) isoforms of circRNAs. Recent advances in nanopore long-read sequencing with circRNA enrichment bring an unprecedented opportunity for investigating the issues. Here, we developed FL-circAS (https://cosbi.ee.ncku.edu.tw/FL-circAS/), which collected such long-read sequencing data of 20 cell lines/tissues and thereby identified 884 636 BSJs with 1 853 692 full-length circRNA isoforms in human and 115 173 BSJs with 135 617 full-length circRNA isoforms in mouse. FL-circAS also provides multiple circRNA features. For circRNA expression, FL-circAS calculates expression levels for each circRNA isoform, cell line/tissue specificity at both the BSJ and isoform levels, and AS entropy for each BSJ across samples. For circRNA biogenesis, FL-circAS identifies reverse complementary sequences and RNA binding protein (RBP) binding sites residing in flanking sequences of BSJs. For functional patterns, FL-circAS identifies potential microRNA/RBP binding sites and several types of evidence for circRNA translation on each full-length circRNA isoform. FL-circAS provides user-friendly interfaces for browsing, searching, analyzing, and downloading data, serving as the first resource for discovering full-length circRNAs at the isoform level.

Circ_0001666 upregulation promotes intestinal epithelial cell fibrosis in pediatric Crohn's disease via the SRSF1/BMP7 axis

The epithelial-mesenchymal transition (EMT) is closely associated with Crohn's disease (CD) related intestinal fibrosis, a condition whose prevalence is increasing annually among children. Recently, the CD marker gene microarray screening revealed an upregulation of circ_0001666 in the colon tissues of CD patients, but its underlying mechanisms remain unclear. In this study, we explored the molecular mechanism of circ_0001666 in regulating EMT-mediated fibrosis in CD in vitro. The levels of circ_0001666 and EMT-associated proteins were assessed in CD clinical samples, and a CD cell model was established using TGF-β1 to induce human intestinal epithelial cells (HIECs). Additionally, the expression levels of genes and proteins related to EMT and fibrosis were analyzed by quantitative real-time PCR and western blot, cell migration, and invasion were assessed via wound healing assay and transwell, respectively, and RNA pull-down and RNA immunoprecipitation assays were performed to verify the relationship between SRSF1 and BMP7 or circ_0001666. Circ_0001666 was overexpressed in the intestinal mucosal tissues of CD patients and was positively correlated with EMT. Silencing circ_0001666 inhibited the migration, invasion, EMT, and fibrosis of HIECs induced by TGF-β1. Mechanistically, circ_0001666 regulated BMP7 expression by interacting with SRSF1. Furthermore, the effects of inhibiting circ_0001666 on HIECs could be partially reversed by overexpressing SRSF1 or silencing BMP7. Collectively, circ_0001666 regulates TGF-β1-induced HIEC migration, invasion, EMT, and fibrosis. Circ_0001666 also promoted EMT-mediated fibrosis by interacting with SRSF1 to accelerate BMP7 mRNA decay. These findings provide new insights into the pathogenesis of CD and suggest that circ_0001666 might be a potential therapeutic target for CD.

CRMSS: predicting circRNA-RBP binding sites based on multi-scale characterizing sequence and structure features

Circular RNAs (circRNAs) are reverse-spliced and covalently closed RNAs. Their interactions with RNA-binding proteins (RBPs) have multiple effects on the progress of many diseases. Some computational methods are proposed to identify RBP binding sites on circRNAs but suffer from insufficient accuracy, robustness and explanation. In this study, we first take the characteristics of both RNA and RBP into consideration. We propose a method for discriminating circRNA-RBP binding sites based on multi-scale characterizing sequence and structure features, called CRMSS. For circRNAs, we use sequence ${k}\hbox{-}{mer}$ embedding and the forming probabilities of local secondary structures as features. For RBPs, we combine sequence and structure frequencies of RNA-binding domain regions to generate features. We capture binding patterns with multi-scale residual blocks. With BiLSTM and attention mechanism, we obtain the contextual information of high-level representation for circRNA-RBP binding. To validate the effectiveness of CRMSS, we compare its predictive performance with other methods on 37 RBPs. Taking the properties of both circRNAs and RBPs into account, CRMSS achieves superior performance over state-of-the-art methods. In the case study, our model provides reliable predictions and correctly identifies experimentally verified circRNA-RBP pairs. The code of CRMSS is freely available at https://github.com/BioinformaticsCSU/CRMSS.

CRIT: Identifying RNA-binding protein regulator in circRNA life cycle via non-negative matrix factorization

Circular RNAs (circRNAs) are endogenous non-coding RNAs that regulate gene expression and participate in carcinogenesis. However, the RNA-binding proteins (RBPs) involved in circRNAs biogenesis and modulation remain largely unclear. We developed the circRNA regulator identification tool (CRIT), a non-negative matrix-factorization-based pipeline to identify regulating RBPs in cancers. CRIT uncovered 73 novel regulators across thousands of samples by effectively leveraging genomics data and functional annotations. We demonstrated that known RBPs involved in circRNA control are significantly enriched in these predictions. Analysis of circRNA-RBP interactions using two large cross-linking immunoprecipitation (CLIP) databases, we validated the consistency between CRIT prediction and the CLIP experiments. Furthermore, newly discovered RBPs are functionally connected with authentic circRNA regulators by various biological associations, such as physical interaction, similar binding motifs, common transcription factor modulation, and co-expression. When analyzing RNA sequencing (RNA-seq) datasets after short hairpin RNA (shRNA)/small interfering RNA (siRNA) knockdown, we found several novel RBPs that can affect global circRNA expression, which strengthens their role in the circRNA life cycle. The above evidence provided independent confirmation that CRIT is a useful tool to capture RBPs in circRNA processing. Finally, we show that authentic regulators are more likely the core splicing proteins and peripheral factors and usually harbor more alterations in the vast majority of cancers.

circRIP: an accurate tool for identifying circRNA-RBP interactions

Circular ribonucleic acids (RNAs) (circRNAs) are formed by covalently linking the downstream splice donor and the upstream splice acceptor. One of the most important functions of circRNAs is mainly exerted through binding RNA-binding proteins (RBPs). However, there is no efficient algorithm for identifying genome-wide circRNA-RBP interactions. Here, we developed a unique algorithm, circRIP, for identifying circRNA-RBP interactions from RNA immunoprecipitation sequencing (RIP-Seq) data. A simulation test demonstrated the sensitivity and specificity of circRIP. By applying circRIP, we identified 95 IGF2BP3-binding circRNAs based on the IGF2BP3 RIP-Seq dataset. We further identified 2823 and 1333 circRNAs binding to >100 RBPs in K562 and HepG2 cell lines, respectively, based on enhanced cross-linking immunoprecipitation (eCLIP) data, demonstrating the significance to survey the potential interactions between circRNAs and RBPs. In this study, we provide an accurate and sensitive tool, circRIP (https://github.com/bioinfolabwhu/circRIP), to systematically identify RBP and circRNA interactions from RIP-Seq and eCLIP data, which can significantly benefit the research community for the functional exploration of circRNAs.

Exosome-Associated circRNAs as Key Regulators of EMT in Cancer

Epithelial-to-mesenchymal transition (EMT) is a dynamic program of cell plasticity aberrantly reactivated in cancer. The crosstalk between tumor cells and the tumoral microenvironment (TME) has a pivotal importance for the induction of the EMT and the progression toward a malignant phenotype. Notably, exosomes are key mediators of this crosstalk as vehicles of specific molecular signals that include the class of circular RNAs (circRNAs). This review specifically focuses on the role of exosome-associated circRNAs as key regulators of EMT in cancer. The relevance of these molecules in regulating the intercellular communication in TME and tumor progression is highlighted. Moreover, the here-presented evidence indicates that exosome-associated circRNA modulation should be taken in account for cancer diagnostic and therapeutic approaches.

Evaluation of methods to detect circular RNAs from single-end RNA-sequencing data

Background

Circular RNA (circRNA), a class of RNA molecule with a loop structure, has recently attracted researchers due to its diverse biological functions and potential biomarkers of human diseases. Most of the current circRNA detection methods from RNA-sequencing (RNA-Seq) data utilize the mapping information of paired-end (PE) reads to eliminate false positives. However, much of the practical RNA-Seq data such as cross-linking immunoprecipitation sequencing (CLIP-Seq) data usually contain single-end (SE) reads. It is not clear how well these tools perform on SE RNA-Seq data.

Results

In this study, we present a systematic evaluation of six advanced RNA-based methods and two CLIP-Seq based methods for detecting circRNAs from SE RNA-Seq data. The performances of the methods are rigorously assessed based on precision, sensitivity, F1 score, and true discovery rate. We investigate the impacts of read length, false positive ratio, sequencing depth and PE mapping information on the performances of the methods using simulated SE RNA-Seq simulated datasets. The real datasets used in this study consist of four experimental RNA-Seq datasets with ≥100bp read length and 124 CLIP-Seq samples from 45 studies that contain mostly short-read (≤50bp) RNA-Seq data. The simulation study shows that the sensitivities of most of the methods can be improved by increasing either read length or sequencing depth, and that the levels of false positive rates significantly affect the precision of all methods. Furthermore, the PE mapping information can improve the method’s precision but can not always guarantee the increase of F1 score. Overall, no method is dominant for all SE RNA-Seq data. The RNA-based methods perform better for the long-read datasets but are worse for the short-read datasets. In contrast, the CLIP-Seq based methods outperform the RNA-Seq based methods for all the short-read samples. Combining the results of these methods can significantly improve precision in the CLIP-Seq data.

Conclusions

The results provide a systematic evaluation of circRNA detection methods on SE RNA-Seq data that would facilitate researchers’ strategies in circRNA analysis.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08329-7).

FUS driven circCNOT6L biogenesis in mouse and human spermatozoa supports zygote development

Circular RNA (circRNA) biogenesis requires a backsplicing reaction, promoted by inverted repeats in cis-flanking sequences and trans factors, such as RNA-binding proteins (RBPs). Among these, FUS plays a key role. During spermatogenesis and sperm maturation along the epididymis such a molecular mechanism has been poorly explored. With this in mind, we chose circCNOT6L as a study case and wild-type (WT) as well as cannabinoid receptor type-1 knock-out (Cb1^−/−) male mice as animal models to analyze backsplicing mechanisms. Our results suggest that spermatozoa (SPZ) have an endogenous skill to circularize mRNAs, choosing FUS as modulator of backsplicing and under CB1 stimulation. A physical interaction between FUS and CNOT6L as well as a cooperation among FUS, RNA Polymerase II (RNApol2) and Quaking (QKI) take place in SPZ. Finally, to gain insight into FUS involvement in circCNOT6L biogenesis, FUS expression was reduced through RNA interference approach. Paternal transmission of FUS and CNOT6L to oocytes during fertilization was then assessed by using murine unfertilized oocytes (NF), one-cell zygotes (F) and murine oocytes undergoing parthenogenetic activation (PA) to exclude a maternal contribution. The role of circCNOT6L as an active regulator of zygote transition toward the 2-cell-like state was suggested using the Embryonic Stem Cell (ESC) system. Intriguingly, human SPZ exactly mirror murine SPZ.

Methods to study circRNA-protein interactions

Circular RNAs (circRNAs) have been studied extensively in the last few years, uncovering functional roles in a diverse range of cell types and organisms. As shown for a few cases, these functions may be mediated by trans-acting factors, in particular RNA-binding proteins (RBPs). However, the specific interaction partners for most circRNAs remain unknown. This is mainly due to technical difficulties in their identification and in differentiating between interactors of circRNAs and their linear counterparts. Here we review the currently used methodology to systematically study circRNA-protein complexes (circRNPs), focusing either on a specific RNA or protein, both on the gene-specific or global level, and discuss advantages and challenges of the available approaches.

CircRNA-Mediated Regulation of Angiogenesis: A New Chapter in Cancer Biology

Angiogenesis is necessary for carcinoma progression and is regulated by a variety of pro- and anti-angiogenesis factors. CircRNAs are RNA molecules that do not have a 5'-cap or a 3'-polyA tail and are involved in a variety of biological functions. While circRNA-mediated regulation of tumor angiogenesis has received much attention, the detailed biological regulatory mechanism remains unclear. In this review, we investigated circRNAs in tumor angiogenesis from multiple perspectives, including its upstream and downstream factors. We believe that circRNAs have natural advantages and great potential for the diagnosis and treatment of tumors, which deserves further exploration.

isoCirc catalogs full-length circular RNA isoforms in human transcriptomes

Circular RNAs (circRNAs) have emerged as an important class of functional RNA molecules. Short-read RNA sequencing (RNA-seq) is a widely used strategy to identify circRNAs. However, an inherent limitation of short-read RNA-seq is that it does not experimentally determine the full-length sequences and exact exonic compositions of circRNAs. Here, we report isoCirc, a strategy for sequencing full-length circRNA isoforms, using rolling circle amplification followed by nanopore long-read sequencing. We describe an integrated computational pipeline to reliably characterize full-length circRNA isoforms using isoCirc data. Using isoCirc, we generate a comprehensive catalog of 107,147 full-length circRNA isoforms across 12 human tissues and one human cell line (HEK293), including 40,628 isoforms ≥500 nt in length. We identify widespread alternative splicing events within the internal part of circRNAs, including 720 retained intron events corresponding to a class of exon-intron circRNAs (EIciRNAs). Collectively, isoCirc and the companion dataset provide a useful strategy and resource for studying circRNAs in human transcriptomes.

The Biomarker and Therapeutic Potential of Circular Rnas in Schizophrenia

Circular RNAs (circRNAs) are endogenous, single-stranded, most frequently non-coding RNA (ncRNA) molecules that play a significant role in gene expression regulation. Circular RNAs can affect microRNA functionality, interact with RNA-binding proteins (RBPs), translate proteins by themselves, and directly or indirectly modulate gene expression during different cellular processes. The affected expression of circRNAs, as well as their targets, can trigger a cascade of events in the genetic regulatory network causing pathological conditions. Recent studies have shown that altered circular RNA expression patterns could be used as biomarkers in psychiatric diseases, including schizophrenia (SZ); moreover, circular RNAs together with other cell molecules could provide new insight into mechanisms of this disorder. In this review, we focus on the role of circular RNAs in the pathogenesis of SZ and analyze their biomarker and therapeutic potential in this disorder.

Circ09863 Regulates Unsaturated Fatty Acid Metabolism by Adsorbing miR-27a-3p in Bovine Mammary Epithelial Cells

Fatty acid composition plays a key role in regulating flavor and quality of milk. Therefore, in order to improve milk quality, it is particularly important to investigate regulatory mechanisms of milk fatty acid metabolism. Circular RNAs (circRNAs) regulate expression genes associated with several biological processes including fatty acid metabolism. In this study, high-throughput sequencing was used to detect differentially expressed genes in bovine mammary tissue at early lactation and peak lactation. Circ09863 profiles were influenced by the lactation stage. Functional studies in bovine mammary epithelial cells (BMECs) revealed that circ09863 promotes triglyceride (TAG) synthesis together with increased content of unsaturated fatty acids (C16:1 and C18:1). These results suggested that circ09863 is partly responsible for modulating fatty acid metabolism. Additionally, software prediction identified a miR-27a-3p binding site in the circ09863 sequence. Overexpression of miR-27a-3p in BMECs led to decreased TAG synthesis. However, overexpression of circ09863 (pcDNA-circ09863) in BMECs significantly reduced expression of miR-27a-3p and enhanced gene expression of fatty acid synthase (FASN), a target of miR-27a-3p. Overall, data suggest that circ09863 relieves the inhibitory effect of miR-27a-3p on FASN expression by binding miR-27a-3p and subsequently regulating TAG synthesis and fatty acid composition. Together, these mechanisms provide new research avenues and theoretical bases to improve milk quality.

Macrophage Long Non-Coding RNAs in Pathogenesis of Cardiovascular Disease

Chronic inflammation is inextricably linked to cardiovascular disease (CVD). Macrophages themselves play important roles in atherosclerosis, as well as acute and chronic heart failure. Although the role of macrophages in CVD pathophysiology is well-recognized, little is known regarding the precise mechanisms influencing their function in these contexts. Long non-coding RNAs (lncRNAs) have emerged as significant regulators of macrophage function; as such, there is rising interest in understanding how these nucleic acids influence macrophage signaling, cell fate decisions, and activity in health and disease. In this review, we summarize current knowledge regarding lncRNAs in directing various aspects of macrophage function in CVD. These include foam cell formation, Toll-like receptor (TLR) and NF-kβ signaling, and macrophage phenotype switching. This review will provide a comprehensive understanding concerning previous, ongoing, and future studies of lncRNAs in macrophage functions and their importance in CVD.