Accurate quantification of circular RNAs identifies extensive circular isoform switching events

The hitchhikers' guide to RNA sequencing and functional analysis

DNA and RNA sequencing technologies have revolutionized biology and biomedical sciences, sequencing full genomes and transcriptomes at very high speeds and reasonably low costs. RNA sequencing (RNA-Seq) enables transcript identification and quantification, but once sequencing has concluded researchers can be easily overwhelmed with questions such as how to go from raw data to differential expression (DE), pathway analysis and interpretation. Several pipelines and procedures have been developed to this effect. Even though there is no unique way to perform RNA-Seq analysis, it usually follows these steps: 1) raw reads quality check, 2) alignment of reads to a reference genome, 3) aligned reads' summarization according to an annotation file, 4) DE analysis and 5) gene set analysis and/or functional enrichment analysis. Each step requires researchers to make decisions, and the wide variety of options and resulting large volumes of data often lead to interpretation challenges. There also seems to be insufficient guidance on how best to obtain relevant information and derive actionable knowledge from transcription experiments. In this paper, we explain RNA-Seq steps in detail and outline differences and similarities of different popular options, as well as advantages and disadvantages. We also discuss non-coding RNA analysis, multi-omics, meta-transcriptomics and the use of artificial intelligence methods complementing the arsenal of tools available to researchers. Lastly, we perform a complete analysis from raw reads to DE and functional enrichment analysis, visually illustrating how results are not absolute truths and how algorithmic decisions can greatly impact results and interpretation.

Evolutionary Landscape of Tea Circular RNAs and Its Contribution to Chilling Tolerance of Tea Plant

Chilling stress threatens the yield and distribution pattern of global crops, including the tea plant (Camellia sinensis), one of the most important cash crops around the world. Circular RNA (circRNA) plays roles in regulating plant growth and biotic/abiotic stress responses. Understanding the evolutionary characteristics of circRNA and its feedbacks to chilling stress in the tea plant will help to elucidate the vital roles of circRNAs. In the current report, we systematically identified 2702 high-confidence circRNAs under chilling stress in the tea plant, and interestingly found that the generation of tea plant circRNAs was associated with the length of their flanking introns. Repetitive sequences annotation and DNA methylation analysis revealed that the longer flanking introns of circRNAs present more repetitive sequences and higher methylation levels, which suggested that repeat-elements-mediated DNA methylation might promote the circRNAs biogenesis in the tea plant. We further detected 250 differentially expressed circRNAs under chilling stress, which were functionally enriched in GO terms related to cold/stress responses. Constructing a circRNA-miRNA-mRNA interaction network discovered 139 differentially expressed circRNAs harboring potential miRNA binding sites, which further identified 14 circRNAs that might contribute to tea plant chilling responses. We further characterized a key circRNA, CSS-circFAB1, which was significantly induced under chilling stress. FISH and silencing experiments revealed that CSS-circFAB1 was potentially involved in chilling tolerance of the tea plant. Our study emphasizes the importance of circRNA and its preliminary role against low-temperature stress, providing new insights for tea plant cold tolerance breeding.

Systematic benchmarking of statistical methods to assess differential expression of circular RNAs

Circular RNAs (circRNAs) are covalently closed transcripts involved in critical regulatory axes, cancer pathways and disease mechanisms. CircRNA expression measured with RNA-seq has particular characteristics that might hamper the performance of standard biostatistical differential expression assessment methods (DEMs). We compared 38 DEM pipelines configured to fit circRNA expression data's statistical properties, including bulk RNA-seq, single-cell RNA-seq (scRNA-seq) and metagenomics DEMs. The DEMs performed poorly on data sets of typical size. Widely used DEMs, such as DESeq2, edgeR and Limma-Voom, gave scarce results, unreliable predictions or even contravened the expected behaviour with some parameter configurations. Limma-Voom achieved the most consistent performance throughout different benchmark data sets and, as well as SAMseq, reasonably balanced false discovery rate (FDR) and recall rate. Interestingly, a few scRNA-seq DEMs obtained results comparable with the best-performing bulk RNA-seq tools. Almost all DEMs' performance improved when increasing the number of replicates. CircRNA expression studies require careful design, choice of DEM and DEM configuration. This analysis can guide scientists in selecting the appropriate tools to investigate circRNA differential expression with RNA-seq experiments.

CYCLeR-a novel tool for the full isoform assembly and quantification of circRNAs

Splicing is one key mechanism determining the state of any eukaryotic cell. Apart from linear splice variants, circular splice variants (circRNAs) can arise via non-canonical splicing involving a back-splice junction (BSJ). Most existing methods only identify circRNAs via the corresponding BSJ, but do not aim to estimate their full sequence identity or to identify different, alternatively spliced circular isoforms arising from the same BSJ. We here present CYCLeR, the first computational method for identifying the full sequence identity of new and alternatively spliced circRNAs and their abundances while simultaneously co-estimating the abundances of known linear splicing isoforms. We show that CYCLeR significantly outperforms existing methods in terms of F score and quantification of transcripts in simulated data. In a in a comparative study with long-read data, we also show the advantages of CYCLeR compared to existing methods. When analysing Drosophila melanogaster data, CYCLeR uncovers biological patterns of circRNA expression that other methods fail to observe.

Expression Profiles and Functional Analysis of Plasma Exosomal Circular RNAs in Acute Myocardial Infarction

Acute myocardial infarction (AMI) is a common cardiovascular disease with high rates of morbidity and mortality globally. The dysregulation of circular RNAs (circRNAs) has been shown to be closely related to various pathological aspects of AMI. However, the function of exosomal circRNAs in AMI has yet to be investigated. The purpose of this study was to investigate the expression profiles of plasma exosomal circRNAs in AMI and explore their potential functionality. The expression profiles of plasma exosomal circRNAs in patients with AMI, stable coronary heart atherosclerotic disease (CAD), and healthy controls were obtained from a GEO expression dataset (GSE159657). We also analyzed bioinformatics functionality, potential pathways, and interaction networks related to the microRNAs associated with the differentially expressed circRNAs. A total of 253 exosomal circRNAs (184 up- and 69 down-regulated) and 182 exosomal circRNAs (94 up- and 88 down-regulated) were identified as being differentially expressed between the control group and the AMI and CAD patients, respectively. Compared with the CAD group, 231 different exosomal circRNAs (177 up- and 54 down-regulated) were identified in the AMI group. Functional analysis suggested that the parental genes of exosomal has_circ_0061776 were significantly enriched in the biological process of lysine degradation. Pathway interaction network analysis further indicated that exosomal has_circ_0061776 was associated with has-miR-133a, has-miR-214, has-miR-423, and has-miR-217 and may play a role in the pathogenesis of AMI through the MAPK signaling pathway. This study identified the differential expression and functionality of exosomal circRNAs in AMI and provided further understanding of the potential pathogenesis of an exosomal circRNA-related competing endogenous RNA (ceRNA) network in AMI.

Association of Circular RNA and Long Non-Coding RNA Dysregulation with the Clinical Response to Immune Checkpoint Blockade in Cutaneous Metastatic Melanoma

Cutaneous melanoma (CM) is the most lethal form of skin cancer if it becomes metastatic, where treatment options and survival chances decrease dramatically. Immunotherapy treatments based on the immunologic checkpoint inhibitors programmed death cell protein 1 (PD-1) and cytotoxic T-lymphocyte antigen 4 (CTLA-4) constituted a main breakthrough in the treatment of metastatic CM, particularly for the achievement of long-term benefits. Even though it is a very promising therapy, resistance to primary immune checkpoint blockade (ICB) arises in about 70% of CM patients treated with a CTLA-4 inhibitor, and 40–65% of CM patients administered with a PD-1-targeting treatment. Some long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) are implicated in triggering pro- and anti-tumorigenic responses to various cancer treatments. The relationship between lncRNAs, circRNAs and ICB immunotherapy has not been explored in cutaneous metastatic melanoma (CMM). The aim of this pilot study is to evaluate the potential role of circRNA and lncRNA expression as pre-treatment predictors of the clinical response to immunotherapy in CMM patients. RNA-seq from 12 formalin-fixed paraffin-embedded (FFPE) samples from the metastatic biopsies of CMM patients treated with nivolumab was used to identify response-associated transcripts. Our findings indicate that specific lncRNAs and circRNAs, probably acting as competitive endogenous RNAs (ceRNAs), are involved in the regulatory networks of the immune response against metastatic melanoma that these patients have under treatment with nivolumab. Moreover, we established a risk score that yields predictions of the overall survival (OS) and progression-free survival (PFS) of CMM patients with high accuracy. This proof-of-principle work provides a possible insight into the function of ceRNAs, contributing to efforts to decipher the complex molecular mechanisms of ICB cancer treatment response.

Identification of circular RNAs in cardiac hypertrophy and cardiac fibrosis

Cardiac hypertrophy initially serves as an adaptive response to physiological and pathological stimuli. Sustained hypertrophy progress to pathological cardiac hypertrophy, cardiac fibrosis and ultimately lead to heart failure, one of the leading medical causes of mortality worldwide. Intervention of pathological cardiac hypertrophy can effectively reduce the occurrence of heart failure. Abundant factors, such as adrenergic, angiotensin, and endothelin (ET-1) receptors, have been shown to participate in the regulation of pathological cardiac hypertrophy. Recently, an increasing number of studies have indicated that circRNA and circRNA-miRNA–mRNA network regulation is indispensable for the posttranscriptional regulation of mRNA in cardiac hypertrophy. In our study, the morphological, cardiac function and pathological changes during cardiac hypertrophy were investigated. RNA sequencing identified 93 circRNAs that were differentially expressed in the TAC_2w group, and 55 circRNAs in the TAC_4w group compared with the sham group. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses identified several significant pathways, including hypertrophic cardiomyopathy, extracellular matrix (ECM)-receptor interaction and focal adhesion. Coexpression analyses were performed for differentially expressed circRNAs and differentially expressed mRNAs. Based on gene set enrichment analysis (GSEA), 8 circRNAs (mmu-Nfkb1_0001, mmu-Smad4_0007, mmu-Hecw2_0009, mmu-Itgbl1_0002, mmu-Lrrc2_0005, mmu-Cpeb3_0007, mmu-Ryr2_0040, and mmu-Rtn4_0001) involved in cardiac hypertrophy and cardiac fibrosis were identified. We validated some key circRNAs by qPCR. The crucial coexpression of circRNA–mRNA and its interaction with miRNA showed the possible mechanism of circRNAs in the process of cardiac dysfunction. Our results may provide promising targets for the treatment of pathological cardiac hypertrophy and fibrosis.

Sex chromosome aneuploidies give rise to changes in the circular RNA profile: A circular transcriptome-wide study of Turner and Klinefelter syndrome across different tissues

Purpose: The landscape of circular RNAs (circRNAs), an important class of non-coding RNAs that regulate gene expression, has never been described in human disorders of sex chromosome aneuploidies. We profiled circRNAs in Turner syndrome females (45,X; TS) and Klinefelter syndrome males (47,XXY; KS) to investigate how circRNAs respond to a missing or an extra X chromosome.

Methods: Samples of blood, muscle and fat were collected from individuals with TS (n = 33) and KS (n = 22) and from male (n = 16) and female (n = 44) controls. CircRNAs were identified using a combination of circRNA identification pipelines (CIRI2, CIRCexplorer2 and circRNA_finder).

Results: Differential expression of circRNAs was observed throughout the genome in TS and KS, in all tissues. The host-genes from which several of these circRNAs were derived, were associated with known phenotypic traits. Furthermore, several differentially expressed circRNAs had the potential to capture micro RNAs that targeted protein-coding genes with altered expression in TS and KS.

Conclusion: Sex chromosome aneuploidies introduce changes in the circRNA transcriptome, demonstrating that the genomic changes in these syndromes are more complex than hitherto thought. CircRNAs may help explain some of the genomic and phenotypic traits observed in these syndromes.

SARS-CoV-2, SARS-CoV, and MERS-CoV encode circular RNAs of spliceosome-independent origin

Circular RNAs (circRNAs) are a newly recognized component of the transcriptome with critical roles in autoimmune diseases and viral pathogenesis. To address the importance of circRNA in RNA viral transcriptome, we systematically identified and characterized circRNAs encoded by the RNA genomes of betacoronaviruses using both bioinformatical and experimental approaches. We predicted 351, 224, and 2764 circRNAs derived from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), SARS-CoV, and Middle East respiratory syndrome coronavirus, respectively. We experimentally identified 75 potential SARS-CoV-2 circRNAs from RNA samples extracted from SARS-CoV-2-infected Vero E6 cells. A systematic comparison of viral and host circRNA features, including abundance, strand preference, length distribution, circular exon numbers, and breakpoint sequences, demonstrated that coronavirus-derived circRNAs had a spliceosome-independent origin. We further showed that back-splice junctions (BSJs) captured by inverse reverse-transcription polymerase chain reaction have different level of resistance to RNase R. Through northern blotting with a BSJ-spanning probe targeting N gene, we identified three RNase R-resistant bands that represent SARS-CoV-2 circRNAs that are detected cytoplasmic by single-molecule and amplified fluorescence in situ hybridization assays. Lastly, analyses of 169 sequenced BSJs showed that both back-splice and forward-splice junctions were flanked by homologous and reverse complementary sequences, including but not limited to the canonical transcriptional regulatory sequences. Our findings highlight circRNAs as an important component of the coronavirus transcriptome, offer important evaluation of bioinformatic tools in the analysis of circRNAs from an RNA genome, and shed light on the mechanism of discontinuous RNA synthesis.

Targeted Analysis of circRNA Expression in Patient Samples by Lexo-circSeq

Recently, circular RNAs (circRNAs) have been extensively studied in animals and plants. circRNAs are generated by backsplicing from the same linear transcripts that are canonically spliced to produce, for example, mature mRNAs. circRNAs exhibit tissue-specific expression and are potentially involved in many diseases, among them cardiovascular diseases. The comprehensive analysis of circRNA expression patterns across larger patient cohorts requires a streamlined and cost-effective workflow designed to meet small input requirements. In this article, we present Lexo-circSeq, a targeted RNA sequencing approach that can profile up to 110 circRNAs and their corresponding linear transcripts in one experiment. We established Lexo-circSeq employing total human heart RNA and show that our protocol can detect depletion of a specific circRNA in hiPSC-derived cardiomyocytes. Finally, Lexo-circSeq was applied to biopsies from patients diagnosed with dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM), respectively. Interestingly, our results indicate that circular-to-linear-ratios for circSLC8A1 and circRBM33 are deregulated in cardiomyopathy.

Usability evaluation of circRNA identification tools: Development of a heuristic-based framework and analysis

BACKGROUND AND OBJECTIVE

Circular RNAs (circRNAs) are endogenous molecules of non-coding RNA that form a covalently closed loop at the 3' and 5' ends. Recently, the role of these molecules in the regulation of gene expression and their involvement in several human pathologies has gained notoriety. The identification of circRNAs is highly dependent on computational methods for analyzing RNA sequencing data. However, bioinformatics software is known to be problematic in terms of usability. Evidence points out that tools for identifying circRNAs can have such problems, negatively impacting researchers in this field. Here we present a heuristic-based framework for evaluating the usability of command-line circRNA identification software.

METHODS

We used heuristics evaluation to comprehensively identify the usability issues in a sample of circRNA identification tools.

RESULTS

We identified 46 usability issues presented individually in four tools. Most of the issues had cosmetic or minor severity. These are unlikely to challenge experienced users but may cause inconvenience for novice users. We also identified severe issues with the potential to harm users regardless of their experience. The areas most affected were the documentation and the installability of the tools.

CONCLUSIONS

With the proposed framework, we formally describe, for the first time, the usability problems that can affect users in this area of circRNA research. We hope that our framework can help researchers evaluate their software's usability during development.

Exploring the cellular landscape of circular RNAs using full-length single-cell RNA sequencing

Previous studies have demonstrated the highly specific expression of circular RNAs (circRNAs) in different tissues and organisms, but the cellular architecture of circRNA has never been fully characterized. Here, we present a collection of 171 full-length single-cell RNA-seq datasets to explore the cellular landscape of circRNAs in human and mouse tissues. Through large-scale integrative analysis, we identify a total of 139,643 human and 214,747 mouse circRNAs in these scRNA-seq libraries. We validate the detected circRNAs with the integration of 11 bulk RNA-seq based resources, where 216,602 high-confidence circRNAs are uniquely detected in the single-cell cohort. We reveal the cell-type-specific expression pattern of circRNAs in brain samples, developing embryos, and breast tumors. We identify the uniquely expressed circRNAs in different cell types and validate their performance in tumor-infiltrating immune cell composition deconvolution. This study expands our knowledge of circRNA expression to the single-cell level and provides a useful resource for exploring circRNAs at this unprecedented resolution.

Studies of circular RNAs have often been limited to the tissue or organism level. Here, authors investigate the comprehensive expression landscape of circRNAs in human and mouse at single-cell resolution, revealing highly specific and dynamic changes of circRNAs during multiple biological processes.

Roles of circRNAs in the Tumorigenesis and Metastasis of HCC: A Mini Review

Circular RNAs (circRNAs) are a class of non-coding RNAs with loop structures that are stable and widely distributed in different tumor tissues. The development of high-throughput sequencing and in silico tools has enabled the discovery of numerous functional circRNAs. Hepatocellular carcinoma (HCC) is a malignant tumor, and the mechanism involved in its progression has remained unclear. In recent years, an increasing number of circRNAs have been identified in HCC, contributing to tumorigenesis and metastasis and with the potential role as biomarkers through competitive endogenous RNAs (ceRNAs) as miRNA sponges or by interacting with RNA binding proteins (RBPs). In this review, we summarize the regulatory roles of circRNAs in HCC development as well as the use of bioinformatics tools in the annotation and prioritization of circRNA and highlight the participation of exosomal circRNAs in HCC metastasis and drug resistance.

DEBKS: A Tool to Detect Differentially Expressed Circular RNAs

Circular RNAs (circRNAs) are involved in various biological processes and disease pathogenesis. However, only a small number of functional circRNAs have been identified among hundreds of thousands of circRNA species, partly because most current methods are based on circular junction counts and overlook the fact that a circRNA is formed from the host gene by back-splicing (BS). To distinguish the expression difference originating from BS or the host gene, we present differentially expressed back-splicing (DEBKS), a software program to streamline the discovery of differential BS events between two rRNA-depleted RNA sequencing (RNA-seq) sample groups. By applying to real and simulated data and employing RT-qPCR for validation, we demonstrate that DEBKS is efficient and accurate in detecting circRNAs with differential BS events between paired and unpaired sample groups. DEBKS is available at https://github.com/yangence/DEBKS as open-source software.

High-throughput analysis of ANRIL circRNA isoforms in human pancreatic islets

The antisense non-coding RNA in the INK locus (ANRIL) is a hotspot for genetic variants associated with cardiometabolic disease. We recently found increased ANRIL abundance in human pancreatic islets from donors with certain Type II Diabetes (T2D) risk-SNPs, including a T2D risk-SNP located within ANRIL exon 2 associated with beta cell proliferation. Recent studies have found that expression of circular species of ANRIL is linked to the regulation of cardiovascular phenotypes. Less is known about how the abundance of circular ANRIL may influence T2D phenotypes. Herein, we sequence circular RNA in pancreatic islets to characterize circular isoforms of ANRIL. We identify several consistently expressed circular ANRIL isoforms whose expression is correlated across dozens of individuals and characterize ANRIL splice sites that are commonly involved in back-splicing. We find that samples with the T2D risk allele in ANRIL exon 2 had higher ratios of circular to linear ANRIL compared to protective-allele carriers, and that higher circular:linear ANRIL was associated with decreased beta cell proliferation. Our study points to a combined involvement of both linear and circular ANRIL species in T2D phenotypes and opens the door for future studies of the molecular mechanisms by which ANRIL impacts cellular function in pancreatic islets.

Circular RNAs as novel biomarkers in triple-negative breast cancer: a systematic review

More effective prognostic and diagnostic tools are urgently required for early detecting and treating triple-negative breast cancer, which is the most acute type of breast cancer because of its lower survival rate, aggressiveness, and non-response to various common treatments. So, it remains the most harmful malignancy for women worldwide. Recently, circular RNAs, as a group of non-coding RNAs, with covalently closed loop and high stability have been discovered, which can modulate gene expression through competing with endogenous microRNA sponges. This finding provided further insight into novel approaches for controlling genes affected in many disorders and malignancies. This review concentrates on the dysregulated expression of circRNAs like their diagnostic and prognostic values in TNBC. This review aims to focus on the abnormal expression of circRNAs and their diagnostic and prognostic values in TNBC. We used PubMed, Embase, and Web of Science databases and ClinicalTrials.gov to systematically search for all relevant clinical studies. This review is based on articles published in databases up to April 2022 with the following keywords: "Circular RNA", "CircRNA", "Triple-Negative Breast Cancer" and "TNBC". We conducted a review of published CircRNA profiled-research articles to identify candidate CircRNA biomarkers for TNBC. The review is registered on JBI at https://jbi.global/systematic-review-register . Accumulating evidence has shown that several circRNAs are downregulated and some are upregulated in TNBC. The results of these studies confirm that circRNAs might be potential biomarkers with the diagnostic, prognostic, and therapeutic target value for TNBC. We also consider the connection between circRNAs and TNBC cell proliferation, apoptosis, metastasis, and chemotherapy resistance and sensitivity.

Urinary extracellular vesicles contain mature transcriptome enriched in circular and long noncoding RNAs with functional significance in prostate cancer

Long noncoding (lnc)RNAs modulate gene expression alongside presenting unexpected source of neoantigens. Despite their immense interest, their ability to be transferred and control adjacent cells is unknown. Extracellular Vesicles (EVs) offer a protective environment for nucleic acids, with pro and antitumourigenic functions by controlling the immune response. In contrast to extracellular nonvesicular RNA, few studies have addressed the full RNA content within human fluids' EVs and have compared them with their tissue of origin. Here, we performed Total RNA-Sequencing on six Formalin-Fixed-Paraffin-Embedded (FFPE) prostate cancer (PCa) tumour tissues and their paired urinary (u)EVs to provide the first whole transcriptome comparison from the same patients. UEVs contain simplified transcriptome with intron-free cytoplasmic transcripts and enriched lnc/circular (circ)RNAs, strikingly common to an independent 20 patients' urinary cohort. Our full cellular and EVs transcriptome comparison within three PCa cell lines identified a set of overlapping 14 uEV-circRNAs characterized as essential for prostate cell proliferation in vitro and 28 uEV-lncRNAs belonging to the cancer-related lncRNA census (CLC2). In addition, we found 15 uEV-lncRNAs, predicted to encode 768 high-affinity neoantigens, and for which three of the encoded-ORF produced detectable unmodified peptides by mass spectrometry. Our dual analysis of EVs-lnc/circRNAs both in urines' and in vitro's EVs provides a fundamental resource for future uEV-lnc/circRNAs phenotypic characterization involved in PCa.

Comprehensive comparison of two types of algorithm for circRNA detection from short-read RNA-Seq

MOTIVATION

Circular RNA is generally formed by the "back-splicing" process between the upstream splice acceptor and the downstream donor in/not in the regulation of the corresponding RNA-binding proteins or cis-elements. Therefore, more and more software packages have been developed and they are mostly based on the identification of the back-spliced junction reads. However, recent studies developed two software tools that can detect circRNA candidates by constructing k-mer table or/and de bruijn graph rather than reads mapping.

RESULTS

Here, we compared the precision, sensitivity and detection efficiency between software tools based on different algorithms. Eleven representative detection tools with two types of algorithm were selected for the overall pipeline analysis of RNA-seq datasets with/without RNase R treatment in two cell lines. Precision, sensitivity, AUC, F1 score and detection efficiency metrics were assessed to compare the prediction tools. Meanwhile, the sensitivity and distribution of highly expressed circRNAs before and after RNase R treatment were also revealed by their enrichment, unaffected and depleted candidate frequencies. Eventually, we found that compared to the k-mer based tools, CIRI2 and KNIFE based on reads mapping had relatively superior and more balanced detection performance regardless of the cell line or RNase R (-/+) datasets.

AVAILABILITY AND IMPLEMENTATION

All predicted results and source codes can be retrieved from https://github.com/luffy563/circRNA_tools_comparison.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of potentially functional circRNAs and prediction of circRNA-miRNA-mRNA regulatory network in periodontitis: Bridging the gap between bioinformatics and clinical needs

Background and Objective

Periodontitis is a multifactorial chronic inflammatory disease that can lead to the irreversible destruction of dental support tissues. As an epigenetic factor, the expression of circRNA is tissue‐dependent and disease‐dependent. This study aimed to identify novel periodontitis‐associated circRNAs and predict relevant circRNA‐periodontitis regulatory network by using recently developed bioinformatic tools and integrating sequencing profiling with clinical information for getting a better and more thorough image of periodontitis pathogenesis, from gene to clinic.

Material and Methods

High‐throughput sequencing and RT‐qPCR were conducted to identify differentially expressed circRNAs in gingival tissues from periodontitis patients. The relationship between upregulated circRNAs expression and probing depth (PD) was performed using Spearman's correlation analysis. Bioinformatic analyses including GO analysis, circRNA‐disease association prediction, and circRNA‐miRNA‐mRNA network prediction were performed to clarify potential regulatory functions of identified circRNAs in periodontitis. A receiver‐operating characteristic (ROC) curve was established to assess the diagnostic significance of identified circRNAs.

Results

High‐throughput sequencing identified 70 differentially expressed circRNAs (68 upregulated and 2 downregulated circRNAs) in human periodontitis (fold change >2.0 and p < .05). The top five upregulated circRNAs were validated by RT‐qPCR that had strong associations with multiple human diseases, including periodontitis. The upregulation of circRNAs were positively correlated with PD (R = .40–.69, p < .05, moderate). A circRNA‐miRNA‐mRNA network with the top five upregulated circRNAs, differentially expressed mRNAs, and overlapped predicted miRNAs indicated potential roles of circRNAs in immune response, cell apoptosis, migration, adhesion, and reaction to oxidative stress. The ROC curve showed that circRNAs had potential value in periodontitis diagnosis (AUC = 0.7321–0.8667, p < .05).

Conclusion

CircRNA‐disease associations were predicted by online bioinformatic tools. Positive correlation between upregulated circRNAs, circPTP4A2, chr22:23101560‐23135351+, circARHGEF28, circBARD1 and circRASA2, and PD suggested function of circRNAs in periodontitis. Network prediction further focused on downstream targets regulated by circRNAs during periodontitis pathogenesis.

Recent progress of circular RNAs in different types of human cancer: Technological landscape, clinical opportunities and challenges (Review)

Circular RNAs (circRNAs) are a novel class of endogenous non‑coding RNAs that have been recently regarded as functionally active. CircRNAs are remarkably stable and known to possess several biological functions such as microRNA sponging, regulating transcription and splicing and occasionally acting as polypeptide‑producing templates. CircRNAs show tissue‑specific expression and have been reported to be associated with the progression of several types of malignancies. Given the recent progress in genome sequencing and bioinformatics techniques, a rapid increment in the biological role of circRNAs has been observed. Concurrently, the patent search from different patent databases shows that the patent number of circRNA is increasing very quickly. These phenomena reveal a rapid development of the technological landscape. In the present review, the recent progress on circRNAs in various kinds of cancer has been investigated and their function as biomarkers or therapeutic targets and their technological landscape have been appreciated. A new insight into circRNAs structure and functional capabilities in cancer has been reviewed. Continually increasing knowledge on their critical role during cancer progression is projecting them as biomarkers or therapeutic targets for various kinds of cancer. Thus, recent updates on the functional role of circRNAs in terms of the technological landscape, clinical opportunities (biomarkers and therapeutic targets), and challenges in cancer have been illustrated.

Circular RNA Profiles in Viremia and ART Suppression Predict Competing circRNA–miRNA–mRNA Networks Exclusive to HIV-1 Viremic Patients

Since the onset of the HIV-1/AIDS epidemic in 1981, 75 million people have been infected with the virus, and the disease remains a public health crisis worldwide. Circular RNAs (circRNAs) are derived from excised exons and introns during backsplicing, a form of alternative splicing. The relevance of unconventional, non-capped, and non-poly(A) transcripts to transcriptomics studies remains to be routinely investigated. Knowledge gaps to be filled are the interface between host-encoded circRNAs and viral replication in chronically progressed patients and upon treatment with antiviral drugs. We implemented a bioinformatic pipeline and repurpose publicly archived RNA sequence reads from the blood of 19 HIV-1-positive patients that previously compared transcriptomes during viremia and viremia suppression by antiretroviral therapy (ART). The in silico analysis identified viremic patients’ circRNA that became undetectable after ART. The circRNAs originated from a subset of host genes enriched in the HDAC biological pathway. These circRNAs and parental mRNAs held in common a small collection of miRNA response elements (MREs), some of which were present in HIV-1 mRNAs. The function of the MRE-containing target mRNA enriched the RNA polymerase II GO pathway. To visualize the interplay between individual circRNA–miRNA–target mRNA, important for HIV-1 and potentially other diseases, an Interactive Circos tool was developed to efficiently parse the intricately competing endogenous network of circRNA–miRNA–mRNA interactions originating from seven circRNA singled out in viremic versus non-viremic patients. The combined downregulation of the identified circRNAs warrants investigation as a novel antiviral targeting strategy.

Circular RNAs as Potential Biomarkers in Breast Cancer

Due to the high heterogeneity and initially asymptomatic nature of breast cancer (BC), the management of this disease depends on imaging together with immunohistochemical and molecular evaluations. These tests allow early detection of BC and patient stratification as they guide clinicians in prognostication and treatment decision-making. Circular RNAs (circRNAs) represent a class of newly identified long non-coding RNAs. These molecules have been described as key regulators of breast carcinogenesis and progression. Moreover, circRNAs play a role in drug resistance and are associated with clinicopathological features in BC. Accumulating evidence reveals a clinical interest in deregulated circRNAs as diagnostic, prognostic and predictive biomarkers. Furthermore, due to their covalently closed structure, circRNAs are highly stable and easily detectable in body fluids, making them ideal candidates for use as non-invasive biomarkers. Herein, we provide an overview of the biogenesis and pleiotropic functions of circRNAs, and report on their clinical relevance in BC.

Integrative analysis reveals a lineage-specific circular RNA landscape for adipo-osteogenesis of human mesenchymal stem cells

Background

The balance between osteogenesis and adipogenesis of mesenchymal stem cells (MSCs) is critical to skeletal development and diseases. As a research hotspot, circular RNAs (circRNAs) have expanded our understanding of a hidden layer of the transcriptome. Yet, their roles during adipo-osteogenesis remain poorly described.

Methods

The identity of human MSCs derived from bone marrow and adipose were first determined by flow cytometry, cellular staining, and quantitative polymerase chain reaction (qPCR). Multi-strategic RNA-sequencing was performed using Poly A, RiboMinus and RiboMinus/RNase R methods. Integrative analysis was performed to identify lineage-specific expressed circRNAs. The structural and expressional characteristics were identified by Sanger sequencing and qPCR, respectively. The regulatory effects of adipogenesis-specific circ-CRLF1 were confirmed using siRNA transcfection and qPCR.

Results

We generated a whole transcriptome map during adipo-osteogenesis based on 10 Poly A, 20 RiboMinus and 20 RiboMinus/ RNase R datasets. A total of 31,326 circRNAs were identified and quantified from ~ 3.4 billion paired-end reads. Furthermore, the integrative analysis revealed that 1166 circRNA genes exhibited strong lineage-specific expression patterns. Their host genes were enriched in distinct biological functions, such as cell adhesion, cytokine signaling, and cell division. We randomly selected and validated the back-spliced junction sites and expression patterns of 12 lineage-specific circRNAs. Functional analysis indicated that circ-CRLF1 negatively regulated adipogenesis.

Conclusions

Our integrative analysis reveals an accurate and generally applicable lineage-specific circRNA landscape for adipo-osteogenesis of MSCs and provides a potential therapeutic target, circ-CRLF1, for the treatment of skeleton-related disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-022-02792-5.

circVAMP3 Drives CAPRIN1 Phase Separation and Inhibits Hepatocellular Carcinoma by Suppressing c-Myc Translation

Previous studies have identified the regulatory roles of circular RNAs (circRNAs) in human cancers. However, the molecular mechanisms of circRNAs in hepatocellular carcinoma (HCC) remain largely unknown. This study screens the expression profile of circRNAs in HCC and identifies circVAMP3 as a significantly downregulated circRNA in HCC tissues. HCC patients with low circVAMP3 expression present poor prognosis. circVAMP3 negatively regulates the proliferation and metastasis of HCC cells in vitro and in vivo by driving phase separation of CAPRIN1 and promoting stress granule formation in cells, which can downregulate the protein level of Myc proto-oncogene protein by inhibiting c-Myc translation. Furthermore, circVAMP3 is widely expressed in many human tissues and is downregulated in related cancers. Therefore, circVAMP3 is a potential prognostic indicator for HCC and may serve as a therapeutic target for HCC treatment.

Identification and Characterization of Rice Circular RNAs Responding to Xanthomonas oryzae pv. oryzae Invasion

Emerging roles of circular RNAs (circRNAs) in various biological processes have advanced our knowledge of transcriptional and posttranscriptional gene regulation. To date, no research has been conducted to explore their roles in the rice-Xanthomonas oryzae pv. oryzae interaction. Therefore, we identified 3,517 circRNAs from rice leaves infected with the highly virulent X. oryzae pv. oryzae strain PXO99^A by using rRNA depleted RNA sequencing technique coupled with the CIRI2 and CIRCexplorer2 pipeline. Characterization analyses showed that these circRNAs were distributed across the whole genome of rice, and most circRNAs arose from exons (85.13%), ranged from 200 to 1,000 bp, and were with a noncanonical GT/AG (including CT/AC equivalent) splicing signal. Functional annotation and enrichment analysis of the host genes that produced the differentially expressed circRNAs (DEcircRNAs) suggested that these identified circRNAs might play an important role in reprogramming rice responses to PXO99^A invasion, mainly by mediating photorespiration and chloroplast, peroxisome, and diterpenoid biosynthesis. Moreover, 31 DEcircRNAs were predicted to act as microRNA decoys in rice. The expression profile of four DEcircRNAs were validated by quantitative real-time PCR with divergent primers, and the back-splicing sites of seven DEcircRNAs were verified by PCR analysis and Sanger sequencing. Collectively, these results inferred a potential functional role of circRNAs in the regulation of rice immunity and provide novel clues about the molecular mechanisms of rice-PXO99^A interaction.

Accurate identification of circRNA landscape and complexity reveals their pivotal roles in human oligodendroglia differentiation

Background

Circular RNAs (circRNAs), a novel class of poorly conserved non-coding RNAs that regulate gene expression, are highly enriched in the human brain. Despite increasing discoveries of circRNA function in human neurons, the circRNA landscape and function in developing human oligodendroglia, the myelinating cells that govern neuronal conductance, remains unexplored. Meanwhile, improved experimental and computational tools for the accurate identification of circRNAs are needed.

Results

We adopt a published experimental approach for circRNA enrichment and develop CARP (CircRNA identification using A-tailing RNase R approach and Pseudo-reference alignment), a comprehensive 21-module computational framework for accurate circRNA identification and quantification. Using CARP, we identify developmentally programmed human oligodendroglia circRNA landscapes in the HOG oligodendroglioma cell line, distinct from neuronal circRNA landscapes. Numerous circRNAs display oligodendroglia-specific regulation upon differentiation, among which a subclass is regulated independently from their parental mRNAs. We find that circRNA flanking introns often contain cis-regulatory elements for RNA editing and are predicted to bind differentiation-regulated splicing factors. In addition, we discover novel oligodendroglia-specific circRNAs that are predicted to sponge microRNAs, which co-operatively promote oligodendroglia development. Furthermore, we identify circRNA clusters derived from differentiation-regulated alternative circularization events within the same gene, each containing a common circular exon, achieving additive sponging effects that promote human oligodendroglia differentiation.

Conclusions

Our results reveal dynamic regulation of human oligodendroglia circRNA landscapes during early differentiation and suggest critical roles of the circRNA-miRNA-mRNA axis in advancing human oligodendroglia development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02621-1.

Circular RNAs in Intervertebral Disc Degeneration: An Updated Review

Low back pain, a common medical condition, could result in severe disability and inflict huge economical and public health burden. Its pathogenesis is attributed to multiple etiological factors, including intervertebral disc degeneration (IDD). Emerging evidence suggests that circular RNAs (circRNAs), a major type of regulatory non-coding RNA, play critical roles in cellular processes that are pertinent to IDD development, including nucleus pulposus cell proliferation and apoptosis as well as extracellular matrix deposition. Increasing number of translational studies also indicated that circRNAs could serve as novel biomarkers for the diagnosis of IDD and/or predicting its clinical outcomes. Our review aims to discuss the recent progress in the functions and mechanisms of newly discovered IDD-related circRNAs.

Emerging Roles of Circular RNAs in Vascular Smooth Muscle Cell Dysfunction

Atherosclerosis is the major pathophysiological basis of cerebrovascular and cardiovascular diseases. Vascular smooth muscle cells (VSMCs) constitute the main structure of vasculature and play important roles in maintaining vascular tone and blood pressure. Many biological processes and cellular signaling events involved in atherosclerogenesis have been shown to converge on deregulating VSMC functions. However, the molecular mechanisms underlying dysfunctional VSMC in atherosclerosis are still poorly defined. Recent evidence revealed that circular RNAs (circRNAs) are closely related to diseases such as degenerative diseases, tumor, congenital diseases, endocrine diseases and cardiovascular diseases. Several studies demonstrated that circRNAs (e.g., circACTA2, Circ-SATB2, circDiaph3, circ_0020397, circTET3, circCCDC66) played critical roles in the regulation of VSMC proliferation, migration, invasion, and contractile-to-synthetic phenotype transformation by sponging microRNAs (e.g., miR-548f-5p, miR-939, miR-148a-5p, miR-138, miR-351-5p, miR-342-3p). This review describes recent progress in the profiling of circRNAs by transcriptome analysis in VSMCs and their molecular functions in regulating VSMC proliferation and migration.

circMRPS35 promotes malignant progression and cisplatin resistance in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the major causes of cancer-related death worldwide. Circular RNAs (circRNAs), a novel class of non-coding RNA, have been reported to be involved in the etiology of various malignancies. However, the underlying cellular mechanisms of circRNAs implicated in the pathogenesis of HCC remain unknown. In this study, we identified a functional RNA, hsa_circ_0000384 (circMRPS35), from public tumor databases using a set of computational analyses, and we further identified that circMRPS35 was highly expressed in 35 pairs of HCC from patients. Moreover, knockdown of the expression of circMRPS35 in Huh-7 and HCC-LM3 cells suppressed their proliferation, migration, invasion, clone formation, and cell cycle in vitro, and it suppressed tumor growth in vivo as well. Mechanically, circMRPS35 sponged microRNA-148a-3p (miR-148a), regulating the expression of Syntaxin 3 (STX3), which modulated the ubiquitination and degradation of phosphatase and tensin homolog (PTEN). Unexpectedly, we detected a peptide encoded by circMRPS35 (circMRPS35-168aa), which was significantly induced by chemotherapeutic drugs and promoted cisplatin resistance in HCC. These results demonstrated that circMRPS35 might be a novel mediator in HCC progress, and they raise the potential of a new biomarker for HCC diagnosis and prognosis, as well as a novel therapeutic target for HCC patients.

Whole Transcriptome Sequencing of Peripheral Blood Shows That Immunity/GnRH/PI3K-Akt Pathways Are Associated With Opioid Use Disorder

BACKGROUND

Opioid use disorder (OUD), which is most commonly exhibited as addiction, is a persistent chronic disease that places a burden on families and society. Various peripheral traits have been linked to OUD in the past, but research on this topic is insufficient.

METHODS

Seven male patients with OUD and 7 male healthy controls with matched demographic and clinical data were enrolled in this study. Peripheral blood RNA was used to construct an rRNA-removed library and a small RNA library. The peripheral transcriptomic differences between the two groups were investigated using RNA-seq. Differentially expressed messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs) and microRNAs (miRNAs) were identified by bioinformatics methods, and functional enrichment analysis with differentially expressed RNAs was performed to investigate the potential biological mechanisms of OUD.

RESULTS

A total of 229 mRNAs (115 upregulated, 114 downregulated), 416 lncRNAs (191 upregulated, 225 downregulated), 17 circRNAs (16 upregulated, 1 downregulated) and 74 miRNAs (42 upregulated, 32 downregulated) were differentially expressed between the OUD group and the healthy control group. Functional enrichment analysis with differentially expressed mRNAs showed that immunity, GnRH secretion, and PI3K-Akt signaling pathways were associated with OUD. Immunity-, JAK-STAT-, and insulin-related pathways were enriched in functional enrichment analysis of target genes predicted by differentially expressed miRNAs.

CONCLUSION

We identified hundreds of differentially expressed genes that were enriched in immunity, GnRH secretion and PI3K-Akt signaling pathways. Some genes with significant changes might be used as potential biomarkers for progression and treatment of OUD.

Non-coding RNA-associated competitive endogenous RNA regulatory networks: Novel diagnostic and therapeutic opportunities for hepatocellular carcinoma

Hepatocellular carcinoma (HCC), as the most prevalent liver malignancy, is annually diagnosed in more than half a million people worldwide. HCC is strongly associated with hepatitis B and C viral infections as well as alcohol abuse. Obesity and nonalcoholic fatty liver disease (NAFLD) also significantly enhance the risk of liver cancer. Despite recent improvements in therapeutic approaches, patients diagnosed in advanced stages show poor prognosis. Accumulating evidence provides support for the regulatory role of non-coding RNAs (ncRNAs) in cancer. There are a variety of reports indicating the regulatory role of microRNAs (miRNAs) in different stages of HCC. Long non-coding RNAs (LncRNAs) exert their effects by sponging miRNAs and controlling the expression of miRNA-targeted genes. Circular RNAs (circRNAs) perform their biological functions by acting as transcriptional regulators, miRNA sponges and protein templates. Diverse studies have illustrated that dysregulation of competing endogenous RNA networks (ceRNETs) is remarkably correlated with HCC-causing diseases such as chronic viral infections, nonalcoholic steatohepatitis and liver fibrosis/cirrhosis. The aim of the current article was to provide an overview of the role and molecular mechanisms underlying the function of ceRNETs that modulate the characteristics of HCC such as uncontrolled cell proliferation, resistance to cell death, metabolic reprogramming, immune escape, angiogenesis and metastasis. The current knowledge highlights the potential of these regulatory RNA molecules as novel diagnostic biomarkers and therapeutic targets in HCC.

Emerging roles of circRNAs in the pathological process of myocardial infarction

Myocardial infarction (MI) is defined as cardiomyocyte death in a clinical context consistent with ischemic insult. MI remains one of the leading causes of morbidity and mortality worldwide. Although there are a number of effective clinical methods for the diagnosis and treatment of MI, further investigation of novel biomarkers and molecular therapeutic targets is required. Circular RNAs (circRNAs), novel non-coding RNAs, have been reported to function mainly by acting as microRNA (miRNA) sponges or binding to RNA-binding proteins (RBPs). The circRNA-miRNA-mRNA (protein) regulatory pathway regulates gene expression and affects the pathological mechanisms of various diseases. Undoubtedly, a more comprehensive understanding of the relationship between MI and circRNA will lay the foundation for the development of circRNA-based diagnostic and therapeutic strategies for MI. Therefore, this review summarizes the pathophysiological process of MI and various approaches to measure circRNA levels in MI patients, tissues, and cells; highlights the significance of circRNAs in the regulation MI pathogenesis and development; and provides potential clinical insight for the diagnosis, prognosis, and treatment of MI.

Using Drosophila to uncover molecular and physiological functions of circRNAs

Circular RNAs (circRNAs) are a class of covalently closed RNA molecules generated by backsplicing. circRNAs are expressed in a tissue-specific manner, accumulate with age in neural tissues, and are highly stable. In many cases, circRNAs are generated at the expense of a linear transcript as back-splicing competes with linear splicing. Some circRNAs regulate gene expression in cis, and some circRNAs can be translated into protein. The advent of deep sequencing and new bioinformatic tools has allowed detection of thousands of circRNAs in eukaryotes. Studying the functions of circRNAs is done using a combination of molecular and genetic methods. The unique genetic tools that can be used in studies of Drosophila melanogaster are ideal for deciphering the functions of circRNAs in vivo. These tools include the GAL4-UAS system, which can be used to manipulate the levels of circRNAs with exquisite temporal and spatial control, and genetic interaction screening, which could be used to identify pathways regulated by circRNAs. Research performed in Drosophila has revealed circRNAs production mechanisms, details of their translation, and their physiological functions. Due to their short lifecycle and the existence of excellent neurodegeneration models, Drosophila can also be used to study the role of circRNAs in aging and age-related disorders. Here, we review molecular and genetic tools and methods for detecting, manipulating, and studying circRNAs in Drosophila.

Quantifying full-length circular RNAs in cancer

Circular RNAs (circRNAs) are abundantly expressed in cancer. Their resistance to exonucleases enables them to have potentially stable interactions with different types of biomolecules. Alternative splicing can create different circRNA isoforms that have different sequences and unequal interaction potentials. The study of circRNA function thus requires knowledge of complete circRNA sequences. Here we describe psirc, a method that can identify full-length circRNA isoforms and quantify their expression levels from RNA sequencing data. We confirm the effectiveness and computational efficiency of psirc using both simulated and actual experimental data. Applying psirc on transcriptome profiles from nasopharyngeal carcinoma and normal nasopharynx samples, we discover and validate circRNA isoforms differentially expressed between the two groups. Compared with the assumed circular isoforms derived from linear transcript annotations, some of the alternatively spliced circular isoforms have 100 times higher expression and contain substantially fewer microRNA response elements, showing the importance of quantifying full-length circRNA isoforms.

Comprehensive Analysis of Large-Scale Transcriptomes from Multiple Cancer Types

Various abnormalities of transcriptional regulation revealed by RNA sequencing (RNA-seq) have been reported in cancers. However, strategies to integrate multi-modal information from RNA-seq, which would help uncover more disease mechanisms, are still limited. Here, we present PipeOne, a cross-platform one-stop analysis workflow for large-scale transcriptome data. It was developed based on Nextflow, a reproducible workflow management system. PipeOne is composed of three modules, data processing and feature matrices construction, disease feature prioritization, and disease subtyping. It first integrates eight different tools to extract different information from RNA-seq data, and then used random forest algorithm to study and stratify patients according to evidences from multiple-modal information. Its application in five cancers (colon, liver, kidney, stomach, or thyroid; total samples n = 2024) identified various dysregulated key features (such as PVT1 expression and ABI3BP alternative splicing) and pathways (especially liver and kidney dysfunction) shared by multiple cancers. Furthermore, we demonstrated clinically-relevant patient subtypes in four of five cancers, with most subtypes characterized by distinct driver somatic mutations, such as TP53, TTN, BRAF, HRAS, MET, KMT2D, and KMT2C mutations. Importantly, these subtyping results were frequently contributed by dysregulated biological processes, such as ribosome biogenesis, RNA binding, and mitochondria functions. PipeOne is efficient and accurate in studying different cancer types to reveal the specificity and cross-cancer contributing factors of each cancer.It could be easily applied to other diseases and is available at GitHub.

A Concise Review on the Role of CircPVT1 in Tumorigenesis, Drug Sensitivity, and Cancer Prognosis

CircPVT1 (hsa_circ_0001821) is a cancer-related circular RNA (circRNA) that originated from a genomic locus on chromosome 8q24. This locus has been previously found to encode the oncogenic long non-coding RNA PVT1. Expression of this circRNA has been found to be upregulated in diverse neoplastic conditions. CircPVT1 acts as a sponge for miR-125a, miR-125b, miR-124-3p, miR-30a-5p, miR-205-5p, miR-423-5p, miR-526b, miR-137, miR-145-5p, miR-497, miR-30d/e, miR-455-5p, miR-29a-3p, miR-204-5p, miR-149, miR-106a-5p, miR-377, miR-3666, miR-203, and miR-199a-5p. Moreover, it can regulate the activities of PI3K/AKT, Wnt5a/Ror2, E2F2, and HIF-1α. Upregulation of circPVT1 has been correlated with decreased survival of patients with different cancer types. In the current review, we explain the oncogenic impact of circPVT1 in different tissues based on evidence from in vitro, in vivo, and clinical investigations.

Circling back to epigenetic regulation in osteosarcoma: Comment on 'Hsa_circ_0088212-mediated miR-520h/APOA1 axis inhibits the osteosarcoma progression' by ' Hao Peng'

•

CircRNAs are useful biomarkers and therapeutic targets in osteosarcoma.

•

The circ_0088212/miR-520h/APOA1 signaling axis highlights dysregulation of opposing signals that facilitate tumor growth.

Osteosarcoma is a genetically complex cancer, thus there are increasing efforts to identify biomarkers and regulators within the epigenome that can better predict patient outcomes and provide new therapeutic targets. In this commentary, we have evaluated the work of Peng and colleagues on the identification and function of the signaling axis circ_0088212/miR_520 h/APOA1 in osteosarcoma. We provide the context of how novel it is to demonstrate the anti-tumor functions of APOA1 and not just correlative gene expression with patient outcome. We further postulate why some studies involving circRNAs in osteosarcoma contradict one another. We conclude that the study performed by Peng and colleagues was performed with enough rigor to give confidence that circ_0088212 is a bona fide tumor suppressor in osteosarcoma.

Sensitive, reliable and robust circRNA detection from RNA-seq with CirComPara2

Circular RNAs (circRNAs) are a large class of covalently closed RNA molecules originating by a process called back-splicing. CircRNAs are emerging as functional RNAs involved in the regulation of biological processes as well as in disease and cancer mechanisms. Current computational methods for circRNA identification from RNA-seq experiments are characterized by low discovery rates and performance dependent on the analysed data set. We developed CirComPara2 (https://github.com/egaffo/CirComPara2), a new automated computational pipeline for circRNA discovery and quantification, which consistently achieves high recall rates without losing precision by combining multiple circRNA detection methods. In our benchmark analysis, CirComPara2 outperformed state-of-the-art circRNA discovery tools and proved to be a reliable and robust method for comprehensive transcriptome characterization.

A survey of circular RNAs in complex diseases: databases, tools and computational methods

Circular RNAs (circRNAs) are a category of novelty discovered competing endogenous non-coding RNAs that have been proved to implicate many human complex diseases. A large number of circRNAs have been confirmed to be involved in cancer progression and are expected to become promising biomarkers for tumor diagnosis and targeted therapy. Deciphering the underlying relationships between circRNAs and diseases may provide new insights for us to understand the pathogenesis of complex diseases and further characterize the biological functions of circRNAs. As traditional experimental methods are usually time-consuming and laborious, computational models have made significant progress in systematically exploring potential circRNA-disease associations, which not only creates new opportunities for investigating pathogenic mechanisms at the level of circRNAs, but also helps to significantly improve the efficiency of clinical trials. In this review, we first summarize the functions and characteristics of circRNAs and introduce some representative circRNAs related to tumorigenesis. Then, we mainly investigate the available databases and tools dedicated to circRNA and disease studies. Next, we present a comprehensive review of computational methods for predicting circRNA-disease associations and classify them into five categories, including network propagating-based, path-based, matrix factorization-based, deep learning-based and other machine learning methods. Finally, we further discuss the challenges and future researches in this field.

Circall: fast and accurate methodology for discovery of circular RNAs from paired-end RNA-sequencing data

BACKGROUND

Circular RNA (circRNA) is an emerging class of RNA molecules attracting researchers due to its potential for serving as markers for diagnosis, prognosis, or therapeutic targets of cancer, cardiovascular, and autoimmune diseases. Current methods for detection of circRNA from RNA sequencing (RNA-seq) focus mostly on improving mapping quality of reads supporting the back-splicing junction (BSJ) of a circRNA to eliminate false positives (FPs). We show that mapping information alone often cannot predict if a BSJ-supporting read is derived from a true circRNA or not, thus increasing the rate of FP circRNAs.

RESULTS

We have developed Circall, a novel circRNA detection method from RNA-seq. Circall controls the FPs using a robust multidimensional local false discovery rate method based on the length and expression of circRNAs. It is computationally highly efficient by using a quasi-mapping algorithm for fast and accurate RNA read alignments. We applied Circall on two simulated datasets and three experimental datasets of human cell-lines. The results show that Circall achieves high sensitivity and precision in the simulated data. In the experimental datasets it performs well against current leading methods. Circall is also substantially faster than the other methods, particularly for large datasets.

CONCLUSIONS

With those better performances in the detection of circRNAs and in computational time, Circall facilitates the analyses of circRNAs in large numbers of samples. Circall is implemented in C++ and R, and available for use at https://www.meb.ki.se/sites/biostatwiki/circall and https://github.com/datngu/Circall.

Widespread translational control regulates retinal development in mouse

Retinal development is tightly regulated to ensure the generation of appropriate cell types and the assembly of functional neuronal circuitry. Despite remarkable advances have been made in understanding regulation of gene expression during retinal development, how translational regulation guides retinogenesis is less understood. Here, we conduct a comprehensive translatome and transcriptome survey to the mouse retinogenesis from the embryonic to the adult stages. We discover thousands of genes that have dynamic changes at the translational level and pervasive translational regulation in a developmental stage-specific manner with specific biological functions. We further identify genes whose translational efficiencies are frequently controlled by changing usage in upstream open reading frame during retinal development. These genes are enriched for biological functions highly important to neurons, such as neuron projection organization and microtubule-based protein transport. Surprisingly, we discover hundreds of previously uncharacterized micropeptides, translated from putative long non-coding RNAs and circular RNAs. We validate their protein products in vitro and in vivo and demonstrate their potentials in regulating retinal development. Together, our study presents a rich and complex landscape of translational regulation and provides novel insights into their roles during retinogenesis.

Evolutionary dynamics of circular RNAs in primates

Many primate genes produce circular RNAs (circRNAs). However, the extent of circRNA conservation between closely related species remains unclear. By comparing tissue-specific transcriptomes across over 70 million years of primate evolution, we identify that within 3 million years circRNA expression profiles diverged such that they are more related to species identity than organ type. However, our analysis also revealed a subset of circRNAs with conserved neural expression across tens of millions of years of evolution. By comparing to species-specific circRNAs, we identified that the downstream intron of the conserved circRNAs display a dramatic lengthening during evolution due to the insertion of novel retrotransposons. Our work provides comparative analyses of the mechanisms promoting circRNAs to generate increased transcriptomic complexity in primates.

Dysregulation of circRNA expression in the peripheral blood of individuals with schizophrenia and bipolar disorder

Circular RNAs (circRNAs) are head-to-tail back-spliced RNA transcripts that have been linked to several biological processes and their perturbation is evident in human disease, including neurological disorders. There is also emerging research suggesting circRNA expression may also be altered in psychiatric and behavioural syndromes. Here, we provide a comprehensive analysis of circRNA expression in peripheral blood mononuclear cells (PBMCs) from 39 patients with schizophrenia and bipolar disorder as well as 20 healthy individuals using deep RNA-seq. We observed systematic alternative splicing leading to a complex and diverse profile of RNA transcripts including 8762 high confidence circRNAs. More specific scrutiny of the circular transcriptome in schizophrenia and bipolar disorder, compared to a non-psychiatric control group, revealed significant dysregulation of 55 circRNAs with a bias towards downregulation. These molecules were predicted to interact with a large number of miRNAs that target genes enriched in psychiatric disorders. Further replication and cross-validation to determine the specificity of these circRNAs across broader diagnostic groups and subgroups in psychiatry will enable their potential utility as biomarkers to be established. KEY MESSAGES: • We identified 8762 high confidence circRNAs with systematic alternative splicing in human PBMCs. • CircRNAs were dysregulated in schizophrenia and bipolar disorder, compared to a non-psychiatric control group. • The DE circRNAs were predicted to interact with miRNAs with target genes enriched in psychiatric disorders. • Some circRNAs have the potential to serve as biomarkers in psychiatry.

Comprehensive profiling of circular RNAs with nanopore sequencing and CIRI-long

Reconstructing the sequence of circular RNAs (circRNAs) from short RNA sequencing reads has proved challenging given the similarity of circRNAs and their corresponding linear messenger RNAs. Previous sequencing methods were unable to achieve high-throughput detection of full-length circRNAs. Here we describe a protocol for enrichment and full-length sequencing of circRNA isoforms using nanopore technology. Circular reverse transcription and size selection achieves a 20-fold higher enrichment of circRNAs from total RNA compared to previous methods. We developed an algorithm, called circRNA identifier using long-read sequencing data (CIRI-long), to reconstruct the sequence of circRNAs. The workflow was validated with simulated data and by comparison to Illumina sequencing as well as quantitative real-time RT-PCR. We used CIRI-long to analyze adult mouse brain samples and systematically profile circRNAs, including mitochondria-derived and transcriptional read-through circRNAs. We identified a new type of intronic self-ligated circRNA that exhibits special splicing and expression patterns. Our method takes advantage of nanopore long reads and enables unbiased reconstruction of full-length circRNA sequences.

Function of Circular RNAs in Fish and Their Potential Application as Biomarkers

Circular RNAs (circRNAs) are an emerging class of regulatory RNAs with a covalently closed-loop structure formed during pre-mRNA splicing. Recent advances in high-throughput RNA sequencing and circRNA-specific computational tools have driven the development of novel approaches to their identification and functional characterization. CircRNAs are stable, developmentally regulated, and show tissue- and cell-type-specific expression across different taxonomic groups. They play a crucial role in regulating various biological processes at post-transcriptional and translational levels. However, the involvement of circRNAs in fish immunity has only recently been recognized. There is also broad evidence in mammals that the timely expression of circRNAs in muscle plays an essential role in growth regulation but our understanding of their expression and function in teleosts is still very limited. Here, we discuss the available knowledge about circRNAs and their role in growth and immunity in vertebrates from a comparative perspective, with emphasis on cultured teleost fish. We expect that the interest in teleost circRNAs will increase substantially soon, and we propose that they may be used as biomarkers for selective breeding of farmed fish, thus contributing to the sustainability of the aquaculture sector.

Dissecting the multi-omics atlas of the exosomes released by human lung adenocarcinoma stem-like cells

Lung adenocarcinoma is heterogeneous and hierarchically organized, with a subpopulation of stem-like cells (CSCs) that reside at the apex of the hierarchy, in which exosomes act as important mediators by transporting specific molecules among different cell populations. Although there have been numerous studies on tumor exosomes, the constituents and functional properties of CSC-derived exosomes are still poorly characterized. Here we present a detail transcriptome and proteome atlas of the exosomes released by human lung adenocarcinoma stem-like cells (LSLCs). The transcriptome analysis indicates the specific patterns of exosomal constituents, including the fragmentation of transcripts and the low-level presence of circular RNAs, and identifies multiple exosomal-enriched mRNAs and lncRNAs. Integrative analysis of transcriptome and proteome data reveals the diverse functions of exosomal-enriched RNAs and proteins, many of which are associated with tumorigenesis. Importantly, several LSLC markers we identified are highly expressed in LSLC-derived exosomes and associate with poor survival, which may serve as promising liquid biopsy biomarkers for lung adenocarcinoma diagnosis. Our study provides a resource for the future elucidation of the functions of tumor-derived exosomes and their regulatory mechanisms in mediating lung cancer development.

AQUARIUM: accurate quantification of circular isoforms using model-based strategy

SUMMARY

Currently, most computational methods estimate the expression of circular RNAs (circRNAs) using the number of sequencing reads that support back-splicing junctions (BSJ) in RNA-seq data, which may introduce biased estimation of circRNA expression due to the uneven distribution of sequencing reads. To overcome this, we previously developed a model-based strategy for circRNA quantification, enabling consideration of sequencing reads from the entire transcript. Yet, the lack of exact transcript structure of circRNAs may limit its accuracy. Here, we proposed a substantially improved circRNA quantification tool, AQUARIUM, by introducing the full-length RNA structure of circular isoforms. We assessed its performance in circRNA quantification using both biological and simulated rRNA-depleted RNA-seq datasets, and demonstrated its superior performance at both BSJ and isoform level.

AVAILABILITY AND IMPLEMENTATION

AQUARIUM is freely available at https://github.com/wanjun-group-seu/AQUARIUM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

RNAdetector: a free user-friendly stand-alone and cloud-based system for RNA-Seq data analysis

Background

RNA-Seq is a well-established technology extensively used for transcriptome profiling, allowing the analysis of coding and non-coding RNA molecules. However, this technology produces a vast amount of data requiring sophisticated computational approaches for their analysis than other traditional technologies such as Real-Time PCR or microarrays, strongly discouraging non-expert users. For this reason, dozens of pipelines have been deployed for the analysis of RNA-Seq data. Although interesting, these present several limitations and their usage require a technical background, which may be uncommon in small research laboratories. Therefore, the application of these technologies in such contexts is still limited and causes a clear bottleneck in knowledge advancement.

Results

Motivated by these considerations, we have developed RNAdetector, a new free cross-platform and user-friendly RNA-Seq data analysis software that can be used locally or in cloud environments through an easy-to-use Graphical User Interface allowing the analysis of coding and non-coding RNAs from RNA-Seq datasets of any sequenced biological species.

Conclusions

RNAdetector is a new software that fills an essential gap between the needs of biomedical and research labs to process RNA-Seq data and their common lack of technical background in performing such analysis, which usually relies on outsourcing such steps to third party bioinformatics facilities or using expensive commercial software.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04211-7.

A Scalable Strand-Specific Protocol Enabling Full-Length Total RNA Sequencing From Single Cells

RNA sequencing (RNAseq) has been widely used to generate bulk gene expression measurements collected from pools of cells. Only relatively recently have single-cell RNAseq (scRNAseq) methods provided opportunities for gene expression analyses at the single-cell level, allowing researchers to study heterogeneous mixtures of cells at unprecedented resolution. Tumors tend to be composed of heterogeneous cellular mixtures and are frequently the subjects of such analyses. Extensive method developments have led to several protocols for scRNAseq but, owing to the small amounts of RNA in single cells, technical constraints have required compromises. For example, the majority of scRNAseq methods are limited to sequencing only the 3' or 5' termini of transcripts. Other protocols that facilitate full-length transcript profiling tend to capture only polyadenylated mRNAs and are generally limited to processing only 96 cells at a time. Here, we address these limitations and present a novel protocol that allows for the high-throughput sequencing of full-length, total RNA at single-cell resolution. We demonstrate that our method produced strand-specific sequencing data for both polyadenylated and non-polyadenylated transcripts, enabled the profiling of transcript regions beyond only transcript termini, and yielded data rich enough to allow identification of cell types from heterogeneous biological samples.