Nanopore long-read sequencing of circRNAs

Global dysregulation of circular RNAs in frontal cortex and whole blood from DM1 and DM2

Myotonic dystrophy type 1 (DM1) and type 2 (DM2) are autosomal dominant neuromuscular disorders associated with expansions of microsatellites, respectively, in DMPK and CNBP. Their pathogenesis is linked to the global aberrant alternative splicing (AAS) of many genes and marks mostly muscular and neuronal tissues, while blood is the least affected. Recent data in DM1 skeletal muscles indicated that abnormalities in RNA metabolism also include global upregulation of circular RNAs (circRNAs). CircRNAs are a heterogeneous group considered splicing errors and by-products of canonical splicing. To elucidate whether circRNA dysregulation is an inherent feature of the myotonic environment, we perform their analysis in the frontal cortex and whole blood of DM1 and DM2 patients. We find a global elevation of circRNAs in both tissues, and its magnitude is neither correlated with the differences in their parental gene expression nor is associated with AAS published earlier. Aberrantly spliced cassette exons of linear transcripts affected in DM1 and DM2 are not among the circularized exons, which unique genomic features prerequisite back-splicing. However, the blueprint of the AAS of linear RNAs is found in a variety of circRNA isoforms. The heterogeneity of circRNAs also originates from the utilization of exonic and intronic cryptic donors/acceptors in back splice junctions, and intron-containing circRNAs are more characteristic of the blood. Overall, this study reveals circRNA dysregulation in various tissues from DM1 and DM2; however, their levels do not correlate with the AAS in linear RNAs, suggesting a potential independent regulatory mechanism underlying circRNA upregulation in myotonic dystrophy.

EasyCircR: Detection and reconstruction of circular RNAs post-transcriptional regulatory interaction networks

Circular RNAs (circRNAs) are regulatory RNAs that play a crucial role in various biological activities and have been identified as potential biomarkers for neurological disorders and cancer. CircRNAs have emerged as significant regulators of gene expression through different mechanisms, including regulation of transcription and splicing, modulation of translation, and post-translational modifications. Additionally, some circRNAs operate as microRNA (miRNA) sponges in the cytoplasm, boosting post-transcriptional expression of target genes by inhibiting miRNA activity. Although existing pipelines can reconstruct circRNAs, identify miRNAs sponged by them, retrieve cascade-regulated mRNAs, and represent the regulatory interactions as complex circRNA-miRNA-mRNA networks, none of the state-of-the-art approaches can discriminate the biological level at which the mRNAs involved in the interactions are regulated, avoiding considering potential target mRNAs not regulated at the post-transcriptional level. EasyCircR is a novel R package that combines circRNA detection and reconstruction with post-transcriptional gene expression analysis (exon-intron split analysis) and miRNA response element prediction. The package enables estimation and visualization of circRNA-miRNA-mRNA interactions through an intuitive Shiny application, leveraging the post-transcriptional regulatory nature of circRNA-miRNA relationship and excluding unrealistic regulatory interactions at the biological level. EasyCircR source code, Docker container and user guide are available at: https://github.com/InfOmics/EasyCircR.

Discovery of circular transcripts of the human BCL2-like 12 (BCL2L12) apoptosis-related gene, using targeted nanopore sequencing, provides new insights into circular RNA biology

Circular RNAs (circRNAs) constitute an RNA type formed by back-splicing. BCL2-like 12 (BCL2L12) is an apoptosis-related gene comprising 7 exons. In this study, we used targeted nanopore sequencing to identify circular BCL2L12 transcripts in human colorectal cancer cells and investigated the effect of circRNA silencing on mRNA expression of the parental gene. In brief, nanopore sequencing following nested PCR amplification of cDNAs of BCL2L12 circRNAs from 7 colorectal cancer cell lines unraveled 46 BCL2L12 circRNAs, most of which described for the first time. Interestingly, 40 novel circRNAs are likely to form via back-splicing between non-canonical back-splice sites residing in highly similar regions of the primary transcripts. All back-splice junctions were validated using next-generation sequencing (NGS) after circRNA enrichment. Surprisingly, 2 novel circRNAs also comprised a poly(A) tract after BCL2L12 exon 7; this poly(A) tract was back-spliced to exon 1, in both cases. Furthermore, the selective silencing of a BCL2L12 circRNA resulted in a subsequent decrease of BCL2L12 mRNA levels in HCT 116 cells, thus providing evidence of parental gene expression regulation by circRNAs. In conclusion, our study led to the discovery of many circular transcripts from a single human gene and provided new insights into circRNA biogenesis and mode of action.

Species-specific circular RNA circDS-1 enhances adaptive evolution in Talaromyces marneffei through regulation of dimorphic transition

Thermal adaptability is a crucial characteristic for mammalian pathogenic fungi that originally inhabit natural ecosystems. Thermally dimorphic fungi have evolved a unique ability to respond to host body temperature by shifting from mycelia to yeast. The high similarity of protein-coding genes between these fungi and their relatives suggests the indispensable but often overlooked roles of non-coding elements in fungal thermal adaptation. Here, we systematically delineated the landscape of full-length circRNAs in both mycelial and yeast conditions of Talaromyces marneffei, a typical thermally dimorphic fungus causing fatal Talaromycosis, by optimizing an integrative pipeline for circRNA detection utilizing next- and third-generation sequencing. We found T. marneffei circRNA demonstrated features such as shorter length, lower abundance, and circularization-biased splicing. We then identified and validated that circDS-1, independent of its parental gene, promotes the hyphae-to-yeast transition, maintains yeast morphology, and is involved in virulence regulation. Further analysis and experiments among Talaromyces confirmed that the generation of circDS-1 is driven by a T. marneffei-specific region in the flanking intron of circDS-1. Together, our findings not only provide fresh insights into the role of circRNA in fungal thermal adaptation but also reveal a novel molecular mechanism for the adaptive evolution of functional circRNAs derived from intronic mutations.

Past, present, and future strategies for detecting and quantifying circular RNA variants

Circular RNAs (circRNAs) are a family of covalently closed RNA transcripts ubiquitous across the eukaryotic kingdom. CircRNAs are generated by a class of alternative splicing called backsplicing, with the resultant circularization of a part of parental RNA producing the characteristic backsplice junction (BSJ). Because of the noncontiguous sequence of the BSJ with respect to the DNA genome, circRNAs remained hidden in plain sight through over a decade of RNA next-generation sequencing, yet over 3 million unique circRNA transcripts have been illuminated in the past decade alone. CircRNAs are expressed in a cell type-specific manner, are highly stable, with many examples of circRNAs being evolutionarily conserved and/or functional in specific contexts. However, circRNAs can be very lowly expressed and predictions of the circRNA context from BSJ-spanning reads alone can confound extrapolation of the exact sequence composition of the circRNA transcript. For these reasons, specific and ultrasensitive detection, combined with enrichment, bespoke bioinformatics pipelines and, more recently, long-read, highly processive sequencing is becoming critical for complete characterization of all circRNA variants. Concomitantly, the need for targeted detection and quantification of specific circRNAs has sparked numerous laboratory-based and commercial approaches to visualize circRNAs in cells and quantify them in biological samples, including biospecimens. This review focuses on advancements in the detection and quantification of circRNAs, with a particular focus on recent next-generation sequencing approaches to bolster detection of circRNA variants and accurately normalize between sequencing libraries.

Circular RNA Formation and Degradation Are Not Directed by Universal Pathways

Circular RNAs (circRNAs) are a class of unique transcripts characterized by a covalently closed loop structure, which differentiates them from conventional linear RNAs. The formation of circRNAs occurs co-transcriptionally and post-transcriptionally through a distinct type of splicing known as back-splicing, which involves the formation of a head-to-tail splice junction between a 5' splice donor and an upstream 3' splice acceptor. This process, along with exon skipping, intron retention, cryptic splice site utilization, and lariat-driven intron processing, results in the generation of three main types of circRNAs (exonic, intronic, and exonic-intronic) and their isoforms. The intricate biogenesis of circRNAs is regulated by the interplay of cis-regulatory elements and trans-acting factors, with intronic Alu repeats and RNA-binding proteins playing pivotal roles, at least in the formation of exonic circRNAs. Various hypotheses regarding pathways of circRNA turnover are forwarded, including endonucleolytic cleavage and exonuclease-mediated degradation; however, similarly to the inconclusive nature of circRNA biogenesis, the process of their degradation and the factors involved remain largely unclear. There is a knowledge gap regarding whether these processes are guided by universal pathways or whether each category of circRNAs requires special tools and particular mechanisms for their life cycles. Understanding these factors is pivotal for fully comprehending the biological significance of circRNAs. This review provides an overview of the various pathways involved in the biogenesis and degradation of different types of circRNAs and explores key factors that have beneficial or adverse effects on the formation and stability of these unique transcripts in higher eukaryotes.

Full-length nanopore sequencing of circular RNA landscape in peripheral blood cells following sequential BNT162b2 mRNA vaccination

Circular RNAs (circRNA) lack 5' or 3' ends; their unique covalently closed structures prevent RNA degradation by exonucleases. These characteristics provide circRNAs with high pharmaceutical stability and biostability relative to current standard-of-care linear mRNAs. CircRNA levels are reportedly associated with certain human diseases, making them novel disease biomarkers and a noncanonical class of therapeutic targets. In this study, the endogenous circRNAs underlying the response to BNT162b2 mRNA vaccination were evaluated. To this end, peripheral blood samples were subjected to full-length sequencing of circRNAs via nanopore sequencing and transcriptome sequencing. Fifteen samples, comprising pre-, first, and second vaccination cohorts, were obtained from five healthcare workers with no history of SARS-CoV-2 infection or previous vaccination. A total of 4706 circRNAs were detected; following full-length sequencing, 4217 novel circRNAs were identified as being specifically expressed during vaccination. These circRNAs were enriched in the binding motifs of stress granule assemblies and SARS-CoV-2 RNA binding proteins, namely poly(A) binding protein cytoplasmic 1 (PABPC1), pumilio RNA binding family member 1 (PUM1), and Y box binding protein 1 (YBX1). Moreover, 489 circRNAs were identified as previously reported miRNA sponges. The differentially expressed circRNAs putatively originated from plasma B cells compared to circRNAs reported in human blood single-cell RNA sequencing datasets. The pre- and post-vaccination differences observed in the circRNA expression landscape in response to the SARS-CoV-2 BNT162b2 mRNA vaccine.

Full-length circRNA sequencing method using low-input RNAs and profiling of circRNAs in MPTP-PD mice on a nanopore platform

Considering the importance of accurate information of full-length (FL) transcripts in functional analysis, researchers prefer to develop new sequencing methods based on third-generation sequencing (TGS) rather than short-read sequencing. Several FL circRNA sequencing strategies have been developed. However, the current methods are inapplicable to low-biomass samples, since a large amount of total RNAs are acquired for circRNA enrichment before library preparation. In this work, we developed an effective method to detect FL circRNAs from a nanogram level (1-100 ng) of total RNAs based on a nanopore platform. Additionally, prior to the library preparation process, we added a series of 24 nt barcodes for each sample to reduce the cost and operating time. Using this method, we profiled circRNA expression in the striatum, hippocampus and cerebral cortex of a Parkinson's disease (PD) mouse model. Over 6% of reads were effective for FL circRNA identification in most datasets. Notably, a reduction in the RNA initial input resulted in a lower correlation between replicates and the detection efficiency for longer circRNA, but the lowest input (1 ng) was able to detect numerous FL circRNAs. Next, we systematically identified over 263 934 circRNAs in PD and healthy mice using the lower-input FL sequencing method, some of which came from 50.52% of PD-associated genes. Moreover, significant changes were observed in the circRNA expression pattern at an isoform level, and high-confidence protein translation evidence was predicted. Overall, we developed an effective method to characterize FL circRNAs from low-input samples and provide a comprehensive insight into the biological function of circRNAs in PD at an isoform level.

Alterations of mRNAs and Non-coding RNAs Associated with Neuroinflammation in Alzheimer's Disease

Alzheimer's disease is a neurodegenerative pathology whose pathognomonic hallmarks are increased generation of β-amyloid (Aβ) peptide, production of hyperphosphorylated (pTau), and neuroinflammation. The last is an alteration closely related to the progression of AD and although it is present in multiple neurodegenerative diseases, the pathophysiological events that characterize neuroinflammatory processes vary depending on the disease. In this article, we focus on mRNA and non-coding RNA alterations as part of the pathophysiological events characteristic of neuroinflammation in AD and the influence of these alterations on the course of the disease through interaction with multiple RNAs related to the generation of Aβ, pTau, and neuroinflammation itself.

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.

Long- and short-read sequencing methods discover distinct circular RNA pools in Lotus japonicus

Circular RNAs (circRNAs) are covalently closed single-stranded RNAs, generated through a back-splicing process that links a downstream 5' site to an upstream 3' end. The only distinction in the sequence between circRNA and their linear cognate RNA is the back splice junction. Their low abundance and sequence similarity with their linear origin RNA have made the discovery and identification of circRNA challenging. We have identified almost 6000 novel circRNAs from Lotus japonicus leaf tissue using different enrichment, amplification, and sequencing methods as well as alternative bioinformatics pipelines. The different methodologies identified different pools of circRNA with little overlap. We validated circRNA identified by the different methods using reverse transcription polymerase chain reaction and characterized sequence variations using nanopore sequencing. We compared validated circRNA identified in L. japonicus to other plant species and showed conservation of high-confidence circRNA-expressing genes. This is the first identification of L. japonicus circRNA and provides a resource for further characterization of their function in gene regulation. CircRNAs identified in this study originated from genes involved in all biological functions of eukaryotic cells. The comparison of methodologies and technologies to sequence, identify, analyze, and validate circRNA from plant tissues will enable further research to characterize the function and biogenesis of circRNA in L. japonicus.

Nanopore-Mediated Sequencing of Circular RNA

Circular RNAs (circRNAs) constitute a group of RNAs defined by a covalent bond between the 5' and 3' end formed by a unique back-splicing event. Most circRNAs are composed of more than one exon, which are spliced together in a linear fashion. This protocol describes methods to sequence full-length circRNA across the back-splicing junction, allowing unambiguous characterization of circRNA-specific exon-intron structures by long-read sequencing (LRS). Two different sequencing approaches are provided: (1) Global circRNA sequencing (the circNick-LRS strategy) relying on circRNA enrichment from total RNA followed by total circRNA long-read sequencing, and (2) targeted circRNA sequencing (the circPanel-LRS strategy) where a preselected panel of circRNA are sequenced without prior circRNA enrichment. Both methods were originally described in Karim et al. (Rahimi et al., Nat Commun 12: 4825, 2021) where they were applied to characterize the exon-intron structure of >10.000 circRNAs in mouse and human brains.

Approaches and challenges in genome-wide circular RNA identification and quantification

Numerous circular RNAs (circRNAs) produced from back-splicing of exon(s) have been recently revealed on a genome-wide scale across species. Although generally expressed at a low level, some relatively abundant circRNAs can play regulatory roles in various biological processes, prompting continuous profiling of circRNA in broader conditions. Over the past decade, distinct strategies have been applied in both transcriptome enrichment and bioinformatic tools for detecting and quantifying circRNAs. Understanding the scope and limitations of these strategies is crucial for the subsequent annotation and characterization of circRNAs, especially those with functional potential. Here, we provide an overview of different transcriptome enrichment, deep sequencing and computational approaches for genome-wide circRNA identification, and discuss strategies for accurate quantification and characterization of circRNA.

Discovery and Comprehensive Characterization of Novel Circular RNAs of the Apoptosis-Related BOK Gene in Human Ovarian and Prostate Cancer Cells, Using Nanopore Sequencing

CircRNAs have become a novel scientific research hotspot, and an increasing number of studies have shed light on their involvement in malignant progression. Prompted by the apparent scientific gap in circRNAs from apoptosis-related genes, such as BOK, we focused on the identification of novel BOK circRNAs in human ovarian and prostate cancer cells. Total RNA was extracted from ovarian and prostate cancer cell lines and reversely transcribed using random hexamer primers. A series of PCR assays utilizing gene-specific divergent primers were carried out. Next, third-generation sequencing based on nanopore technology followed by extensive bioinformatics analysis led to the discovery of 23 novel circRNAs. These novel circRNAs consist of both exonic and intronic regions of the BOK gene. Interestingly, the exons that form the back-splice junction were truncated in most circRNAs, and multiple back-splice sites were found for each BOK exon. Moreover, several BOK circRNAs are predicted to sponge microRNAs with a key role in reproductive cancers, while the presence of putative open reading frames indicates their translational potential. Overall, this study suggests that distinct alternative splicing events lead to the production of novel BOK circRNAs, which could come into play in the molecular landscape and clinical investigation of ovarian and prostate cancer.

Roles and mechanisms of CircRNAs in ovarian cancer

Ovarian cancer (OC) is one of the female malignancies with nearly 45% 5-year survival rate. Circular RNAs (circRNAs), a kind of single-stranded non-coding RNAs, are generated from the back-splicing of cellular housekeeping noncoding RNAs and precursor messenger RNAs. Recent studies revealed that circRNAs have different biological function, including sponging miRNAs, encoding micropeptides, regulating stability of cytoplasmic mRNAs, affecting transcription and splicing, via interacting with DNA, RNA and proteins. Due to their stability, circRNAs have the potential of acting as biomarkers and treatment targets. In this review, we briefly illustrate the biogenesis mechanism and biological function of circRNAs in OC, and make a perspective of circRNAs drug targeting immune responses and signaling pathways in OC. This article can provide a systematic view into the current situation and future of circRNAs in OC.

Biogenesis and Regulatory Roles of Circular RNAs

Covalently closed, single-stranded circular RNAs can be produced from viral RNA genomes as well as from the processing of cellular housekeeping noncoding RNAs and precursor messenger RNAs. Recent transcriptomic studies have surprisingly uncovered that many protein-coding genes can be subjected to backsplicing, leading to widespread expression of a specific type of circular RNAs (circRNAs) in eukaryotic cells. Here, we discuss experimental strategies used to discover and characterize diverse circRNAs at both the genome and individual gene scales. We further highlight the current understanding of how circRNAs are generated and how the mature transcripts function. Some circRNAs act as noncoding RNAs to impact gene regulation by serving as decoys or competitors for microRNAs and proteins. Others form extensive networks of ribonucleoprotein complexes or encode functional peptides that are translated in response to certain cellular stresses. Overall, circRNAs have emerged as an important class of RNAmolecules in gene expression regulation that impact many physiological processes, including early development, immune responses, neurogenesis, and tumorigenesis.

Best practice standards for circular RNA research

Circular RNAs (circRNAs) are formed in all domains of life and via different mechanisms. There has been an explosion in the number of circRNA papers in recent years; however, as a relatively young field, circRNA biology has an urgent need for common experimental standards for isolating, analyzing, expressing and depleting circRNAs. Here we propose a set of guidelines for circRNA studies based on the authors' experience. This Perspective will specifically address the major class of circRNAs in Eukarya that are generated by a spliceosome-catalyzed back-splicing event. We hope that the implementation of best practice principles for circRNA research will help move the field forward and allow a better functional understanding of this fascinating group of RNAs.

The emerging roles and potential applications of circular RNAs in ovarian cancer: a comprehensive review

Ovarian cancer (OC) is among the most common human malignancies and the first cause of deaths among gynecologic cancers. Early diagnosis can help improving prognosis in those patients, and accordingly exploring novel molecular mechanisms may lead to find therapeutic targets. Circular RNAs (circRNAs) comprise a group of non-coding RNAs in multicellular organisms, which are identified with characteristic circular structure. CircRNAs have been found with substantial functions in regulating gene expression through interacting with RNA-binding proteins, targeting microRNAs, and transcriptional regulation. They have been found to be involved in regulating several critical processes such as cell growth, and death, organ development, signal transduction, and tumorigenesis. Accordingly, circRNAs have been implicated in a number of human diseases including malignancies. They are particularly reported to contribute to several hallmarks of cancer leading to cancer development and progression, although a number also are described with tumor-suppressor function. In OC, circRNAs are linked to regulation of cell growth, invasiveness, metastasis, angiogenesis, and chemoresistance. Notably, clinical studies also have shown potentials in diagnosis, prediction of prognosis, and therapeutic targets for OC. In this review, I have an overview to the putative mechanisms, and functions of circRNAs in regulating OC pathogenesis in addition to their clinical potentials.

PlantcircBase 7.0: Full-length transcripts and conservation of plant circRNAs

Circular RNA (circRNA) is a special type of non-coding RNA that participates in diverse biological processes in both animals and plants. Five years ago, we developed a comprehensive plant circRNA database (PlantcircBase), which has attracted much attention from the plant circRNA community. Here, we report an updated PlantcircBase (v.7.0), which contains 171,118 circRNAs from 21 plant species. Over 31,000 of the circRNAs have full-length sequences constructed based on analysis of 749 bulk RNA sequencing (RNA-seq) datasets downloaded from the public domain and Nanopore long-read sequencing results of rice RNAs newly generated in this study. A plant multiple conservation score (PMCS), based on the conservation of both sequence and expression profiles, was calculated for each circRNA to quantify and compare the conservation of all circRNAs. A new parameter, plant circRNA confidence level (PCCL), is introduced to measure the identity reliability of each circRNA based on experimental validation results and the number of references that support the circRNA. All this information and other details of circRNAs can be browsed, searched, and downloaded from PlantcircBase 7.0, which also provides online bioinformatics tools for visualization and sequence alignment. PlantcircBase 7.0 is publicly and freely accessible at http://ibi.zju.edu.cn/plantcircbase/.

Circular RNAs: Characterization, cellular roles, and applications

Most circular RNAs are produced from the back-splicing of exons of precursor mRNAs. Recent technological advances have in part overcome problems with their circular conformation and sequence overlap with linear cognate mRNAs, allowing a better understanding of their cellular roles. Depending on their localization and specific interactions with DNA, RNA, and proteins, circular RNAs can modulate transcription and splicing, regulate stability and translation of cytoplasmic mRNAs, interfere with signaling pathways, and serve as templates for translation in different biological and pathophysiological contexts. Emerging applications of RNA circles to interfere with cellular processes, modulate immune responses, and direct translation into proteins shed new light on biomedical research. In this review, we discuss approaches used in circular RNA studies and the current understanding of their regulatory roles and potential applications.