The mechanism and detection of alternative splicing events in circular RNAs

Novel circular RNAs of the apoptosis-related BAX and BCL2L12 genes identified in a chronic lymphocytic leukemia cell line using nanopore sequencing

Circular RNAs (circRNAs), a novel RNA type generated by back-splicing, are key regulators of gene expression, with deregulated expression and established involvement in leukemia. The products of BCL2 and its homologs, including BAX and BCL2L12, are implicated in chronic lymphocytic leukemia (CLL). However, to the best of our knowledge, nothing is known about circRNAs produced by these two genes and their role in CLL. We sought to further elucidate the contribution of BAX and BCL2L12 in CLL by unraveling the identity, localization, and potential role of their circRNAs. Therefore, total RNA from the EHEB cell line and peripheral blood mononuclear cells (PBMCs) of CLL patients and non-leukemic blood donors was extracted and reverse-transcribed using random hexamers. Next, nested PCRs with divergent primers were performed and the purified PCR products were subjected to 3rd generation nanopore sequencing. Nested PCRs were also applied to first-strand cDNAs synthesized from total RNA extracts of PBMCs from CLL patients and non-leukemic blood donors. Lastly, a single-molecule resolution fluorescent in situ hybridization method called circFISH was used to visualize the circRNA distribution in EHEB cells. We discovered several novel circRNAs produced by BAX and BCL2L12, which were characterized by great exon structure diversity. In addition, intriguing findings regarding their formation emerged. Interestingly, visualization of the most abundant circRNAs showed distinct intracellular localization. Moreover, a complex BAX and BCL2L12 circRNA expression pattern was revealed in CLL patients and non-leukemic blood donors. Our data suggest a multifaceted role of BAX and BCL2L12 circRNAs in B-cell CLL.

JCcirc: circRNA full-length sequence assembly through integrated junction contigs

Recent studies have shed light on the potential of circular RNA (circRNA) as a biomarker for disease diagnosis and as a nucleic acid vaccine. The exploration of these functionalities requires correct circRNA full-length sequences; however, existing assembly tools can only correctly assemble some circRNAs, and their performance can be further improved. Here, we introduce a novel feature known as the junction contig (JC), which is an extension of the back-splice junction (BSJ). Leveraging the strengths of both BSJ and JC, we present a novel method called JCcirc (https://github.com/cbbzhang/JCcirc). It enables efficient reconstruction of all types of circRNA full-length sequences and their alternative isoforms using splice graphs and fragment coverage. Our findings demonstrate the superiority of JCcirc over existing methods on human simulation datasets, and its average F1 score surpasses CircAST by 0.40 and outperforms both CIRI-full and circRNAfull by 0.13. For circRNAs below 400 bp, 400-800 bp, 800 bp-1200 bp and above 1200 bp, the correct assembly rates are 0.13, 0.09, 0.04 and 0.03 higher, respectively, than those achieved by existing methods. Moreover, JCcirc also outperforms existing assembly tools on other five model species datasets and real sequencing datasets. These results show that JCcirc is a robust tool for accurately assembling circRNA full-length sequences, laying the foundation for the functional analysis of circRNAs.

The epigenetic landscape of oligodendrocyte lineage cells

The epigenetic landscape of oligodendrocyte lineage cells refers to the cell-specific modifications of DNA, chromatin, and RNA that define a unique gene expression pattern of functionally specialized cells. Here, we focus on the epigenetic changes occurring as progenitors differentiate into myelin-forming cells and respond to the local environment. First, modifications of DNA, RNA, nucleosomal histones, key principles of chromatin organization, topologically associating domains, and local remodeling will be reviewed. Then, the relationship between epigenetic modulators and RNA processing will be explored. Finally, the reciprocal relationship between the epigenome as a determinant of the mechanical properties of cell nuclei and the target of mechanotransduction will be discussed. The overall goal is to provide an interpretative key on how epigenetic changes may account for the heterogeneity of the transcriptional profiles identified in this lineage.

Circular RNAs and tuberculosis infection

Tuberculosis (TB) is a deadly infectious disease caused by Mycobacterium tuberculosis (Mtb) that affects the immune system chronically. Therefore, effective control and treatment of tuberculosis requires rapid and accurate diagnostic strategies. Tuberculosis has always been a global burden on health, social and economic systems due to the lack of standard curative and diagnostic (bio)markers. Accordingly, the management and monitoring of patients with active TB at the primary care level may be possible through new, rapid and cost-effective non-sputum-based diagnostic procedures. Biomarkers can help diagnose various diseases, including circular RNA (circRNA), which has recently been introduced as an endogenous, abundant and stable RNA in the cytoplasm with unique tissue specificity. There are frequent reports of circRNA involvement in many pathological and physiological processes in human beings. Recent studies have highlighted the presence of circRNAs in serum and their role as promising biomarkers in the diagnosis of the disease, potentially due to the continuous, stable, closed covalent circular structures and lack of easy degradation by nucleases. The purpose of this review article is to scrutinize the behavior of circulating plasma RNAs in relation to the pathogenesis and diagnosis of tuberculosis.

Competing Endogenous RNA (ceRNA) Networks and Splicing Switches in Cervical Cancer: HPV Oncogenesis, Clinical Significance and Therapeutic Opportunities

Cervical cancer (CC) is the primary cause of female cancer fatalities in low-middle-income countries (LMICs). Persistent infections from the human papillomavirus (HPV) can result in cervical cancer. However, numerous different factors influence the development and progression of cervical cancer. Transcriptomic knowledge of the mechanisms with which HPV causes cervical cancer pathogenesis is growing. Nonetheless, there is an existing gap hindering the development of therapeutic approaches and the improvement of patient outcomes. Alternative splicing allows for the production of numerous RNA transcripts and protein isoforms from a single gene, increasing the transcriptome and protein diversity in eukaryotes. Cancer cells exhibit astounding transcriptome modifications by expressing cancer-specific splicing isoforms. High-risk HPV uses cellular alternative splicing events to produce viral and host splice variants and proteins that drive cancer progression or contribute to distinct cancer hallmarks. Understanding how viruses utilize alternative splicing to drive pathogenesis and tumorigenesis is essential. Although research into the role of miRNAs in tumorigenesis is advancing, the function of other non-coding RNAs, including lncRNA and circRNA, has been understudied. Through their interaction with mRNA, non-coding RNAs form a network of competing endogenous RNAs (ceRNAs), which regulate gene expression and promote cervical cancer development and advancement. The dysregulated expression of non-coding RNAs is an understudied and tangled process that promotes cervical cancer development. This review will present the role of aberrant alternative splicing and immunosuppression events in HPV-mediated cervical tumorigenesis, and ceRNA network regulation in cervical cancer pathogenesis will also be discussed. Furthermore, the therapeutic potential of splicing disruptor drugs in cervical cancer will be deliberated.

Alternative splicing signatures of congenital heart disease and induced pluripotent stem cell-derived cardiomyocytes from congenital heart disease patients

Congenital heart disease (CHD) is the most serious congenital defect in newborns with higher mortality. Alternative splicing (AS) plays an essential role in numerous heart diseases. However, our understanding of the link between mRNA splicing and CHD in humans is limited. Here, we try to investigate the genome-wide AS events in CHD using bioinformatics methods. We collected available RNA-seq datasets of CHD-induced pluripotent stem cell-cardiomyocytes (iPSC-CMs) (including single ventricle disease [SVD] and tetralogy of Fallot [TOF]) and non-CHD from the Gene Expression Omnibus database. Then, we unprecedentedly performed AS profiles in CHD-iPSC-CMs and non-CHD-iPSC-CMs. The rMAPS was used to generate RNA-maps for the analysis of RNA-binding proteins' (RBPs) binding sites. We used StringTie to identify and quantify the transcripts from aligned RNA-Seq reads. A quantification matrix was generated with respect to different groups by extracting the transcripts per million values from StringTie outputs. Then, this matrix was used for correlation analysis between the expression level of RBP and AS level. Finally, we validated our AS results using RNA-seq data from CHD and non-CHD patient tissue samples. We identified CHD-related AS events using CHD-iPSC-CMs and CHD samples from patients. The results showed that functional enrichment of abnormal AS in SVD and TOF was transcription factor-related. Using rMAPS, RNA-binding proteins which regulated these AS were also determined, and RBP-AS regulatory network was constructed. Overall, we identified abnormal AS in CHD-iPSC-CMs and CHD samples from patients. We predicted AS regulators in SVD and TOF, respectively. At last, we concluded that AS played a key role in the pathogenesis of CHD.

CircRNA-Associated CeRNAs Regulatory Axes in Retinoblastoma: A Systematic Scoping Review

Retinoblastoma (RB) is one of the most common childhood cancers caused by RB gene mutations (tumor suppressor gene in various patients). A better understanding of molecular pathways and the development of new diagnostic approaches may lead to better treatment for RB patients. The number of studies on ceRNA axes is increasing, emphasizing the significance of these axes in RB. Circular RNAs (circRNAs) play a vital role in competing endogenous RNA (ceRNA) regulatory axes by sponging microRNAs and regulating gene expression. Because of the broadness of ceRNA interaction networks, they may assist in investigating treatment targets in RB. This study conducted a systematic scoping review to evaluate verified loops of ceRNA in RB, focusing on the ceRNA axis and its relationship to circRNAs. This scoping review was carried out using a six-step strategy and the Prisma guideline, and it involved systematically searching the publications of seven databases. Out of 363 records, sixteen articles were entirely consistent with the defined inclusion criteria and were summarized in the relevant table. The majority of the studies focused on the circRNAs circ_0000527, circ_0000034, and circTET1, with approximately two-fifths of the studies focusing on a single circRNA. Understanding the many features of this regulatory structure may help elucidate RB’s unknown causative factors and provide novel molecular potential therapeutic targets and medical fields.

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.

Evaluation of CircRNA Sequence Assembly Methods Using Long Reads

The functional study on circRNAs has been increasing in the past decade due to its important roles in micro RNA sponge, protein coding, the initiation, and progression of diseases. The study of circRNA functions depends on the full-length sequences of circRNA, and current sequence assembly methods based on short reads face challenges due to the existence of linear transcript. Long reads produced by long-read sequencing techniques such as Nanopore technology can cover full-length sequences of circRNA and therefore can be used to evaluate the correctness and completeness of circRNA full sequences assembled from short reads of the same sample. Using long reads of the same samples, one from human and the other from mouse, we have comprehensively evaluated the performance of several well-known circRNA sequence assembly algorithms based on short reads, including circseq_cup, CIRI_full, and CircAST. Based on the F1 score, the performance of CIRI-full was better in human datasets, whereas in mouse datasets CircAST was better. In general, each algorithm was developed to handle special situations or circumstances. Our results indicated that no single assembly algorithm generated better performance in all cases. Therefore, these assembly algorithms should be used together for reliable full-length circRNA sequence reconstruction. After analyzing the results, we have introduced a screening protocol that selects out exonic circRNAs with full-length sequences consisting of all exons between back splice sites as the final result. After screening, CIRI-full showed better performance for both human and mouse datasets. The average F1 score of CIRI-full over four circRNA identification algorithms increased from 0.4788 to 0.5069 in human datasets, and it increased from 0.2995 to 0.4223 in mouse datasets.

Circular RNA: A promising new star for the diagnosis and treatment of colorectal cancer

Background

Colorectal cancer (CRC) is one of the most common malignant tumors of the digestive tract. According to the research of circular RNAs in the CRC field, compared with linear RNAs, circular RNAs are a special type of noncoding RNA that are covalently closed circular structures, which have no 5' cap structure and 3' polyA tail and are not affected by RNA exonuclease and actinomycin D.

Biological functions

Notably, circular RNAs have a high degree of stability and potential effect on gene regulation. Meanwhile, circular RNAs are involved in the sponge action of microRNAs and mediate protein translation and direct binding, alternative splicing, and histone modification.

Relationships with CRC

Studies have shown that circular RNAs are related to the proliferation, invasion, recurrence, metastasis, ferroptosis, apoptosis, and chemotherapy resistance of CRC.

Conclusions

This article provides a brief review based on the source, structural characteristics, mechanisms, biological functions of circular RNAs, and the relationships between CRC.

From circRNAs to fusion circRNAs in hematological malignancies

Circular RNAs (circRNAs) represent a type of endogenous noncoding RNA generated by back-splicing events. Unlike the majority of RNAs, circRNAs are covalently closed, without a 5' end or a 3' poly(A) tail. A few circRNAs can be associated with polysomes, suggesting a protein-coding potential. CircRNAs are not degraded by RNA exonucleases or ribonuclease R and are enriched in exosomes. Recent developments in experimental methods coupled with evolving bioinformatic approaches have accelerated functional investigation of circRNAs, which exhibit a stable structure, a long half-life, and tumor specificity and can be extracted from body fluids and used as potential biological markers for tumors. Moreover, circRNAs may regulate the occurrence and development of cancers and contribute to drug resistance through a variety of molecular mechanisms. Despite the identification of a growing number of circRNAs, their effects in hematological cancers remain largely unknown. Recent studies indicate that circRNAs could also originate from fusion genes (fusion circRNAs, f-circRNAs) next to chromosomal translocations, which are considered the primary cause of various cancers, notably hematological malignancies. This Review will focus on circRNAs and f-circRNAs in hematological cancers.

Functions and mechanisms of circular RNAs in regulating stem cell differentiation

Stem cells are a class of undifferentiated cells with great self-renewal and differentiation capabilities that can differentiate into mature cells in specific tissue types. Stem cell differentiation plays critical roles in body homoeostasis, injury repair and tissue generation. The important functions of stem cell differentiation have resulted in numerous studies focusing on the complex molecular mechanisms and various signalling pathways controlling stem cell differentiation. Circular RNAs (circRNAs) are a novel class of noncoding RNAs with a covalently closed structure present in eukaryotes. Numerous studies have highlighted important biological functions of circRNAs, and they play multiple regulatory roles in various physiological and pathological processes. Importantly, multiple lines of evidence have shown the abnormal expression of numerous circRNAs during stem cell differentiation, and some play a role in regulating stem cell differentiation, highlighting the role of circRNAs as novel biomarkers of stem cell differentiation and novel targets for stem cell-based therapy. In this review, we systematically summarize and discuss recent advances in our understanding of the roles and underlying mechanisms of circRNAs in modulating stem cell differentiation, thus providing guidance for future studies to investigate stem cell differentiation and stem cell-based therapy.Abbreviations: CircRNAs: circular RNAs; ESCs: embryonic stem cells; ADSCs: adipose-derived mesenchymal stem cells; ecircRNAs: exonic circRNAs; EIciRNAs: exon-intron circRNAs; eiRNAs: circular intronic RNAs; tricRNAs: tRNA intronic circRNAs; pol II: polymerase II; snRNP: small nuclear ribonucleoprotein; m6A: N6-methyladenosine; AGO2: Argonaute 2; RBPs: RNA-binding proteins; MBNL: muscleblind-like protein 1; MSCs: mesenchymal stem cells; hiPSCs: human induced pluripotent stem cells; hiPSC-CMs: hiPSC-derived cardiomyocytes; hBMSCs: human bone marrow mesenchymal stem cells; hADSCs: human adipose-derived mesenchymal stem cells; hDPSCs: human dental pulp stem cells; RNA-seq: high-throughput RNA sequencing; HSCs: haematopoietic stem cells; NSCs: neural stem cells; EpSCs: epidermal stem cells; hESCs: human embryonic stem cells; mESCs: murine embryonic stem cells; MNs: motor neurons; SSUP: small subunit processome; BMSCs: bone marrow-derived mesenchymal stem cells; OGN: osteoglycin; GIOP: glucocorticoid‑induced osteoporosis; CDR1as: cerebellar degeneration-related protein 1 transcript; SONFH: steroid-induced osteogenesis of the femoral head; rBMSCs: rat bone marrow-derived mesenchymal stem cells; QUE: quercetin; AcvR1b: activin A receptor type 1B; BSP: bone sialoprotein; mADSCs: mouse ADSCs; PTBP1: polypyrimidine tract-binding protein; ER: endoplasmic reticulum; hUCMSCs: MSCs derived from human umbilical cord; MSMSCs: maxillary sinus membrane stem cells; SCAPs: stem cells from the apical papilla; MyoD: myogenic differentiation protein 1; MSTN: myostatin; MEF2C: myocyte enhancer factor 2C; BCLAF1: BCL2-associated transcription factor 1; EpSCs: epidermal stem cells; ISCs: intestinal stem cells; NSCs: neural stem cells; Lgr5+ ISCs: crypt base columnar cells; ILCs: innate lymphoid cells.