Genome-wide analysis of drosophila circular RNAs reveals their structural and sequence properties and age-dependent neural accumulation

Breast Cancer; Discovery of Novel Diagnostic Biomarkers, Drug Resistance, and Therapeutic Implications

Breast cancer is the second most reported cancer in women with high mortality causing millions of cancer-related deaths annually. Early detection of breast cancer intensifies the struggle towards discovering, developing, and optimizing diagnostic biomarkers that can improve its prognosis and therapeutic outcomes. Breast cancer-associated biomarkers comprise macromolecules, such as nucleic acid (DNA/RNA), proteins, and intact cells. Advancements in molecular technologies have identified all types of biomarkers that are exclusively studied for diagnostic, prognostic, drug resistance, and therapeutic implications. Identifying biomarkers may solve the problem of drug resistance which is a challenging obstacle in breast cancer treatment. Dysregulation of non-coding RNAs including circular RNAs (circRNAs) and microRNAs (miRNAs) initiates and progresses breast cancer. The circulating multiple miRNA profiles promise better diagnostic and prognostic performance and sensitivity than individual miRNAs. The high stability and existence of circRNAs in body fluids make them a promising new diagnostic biomarker. Many therapeutic-based novels targeting agents have been identified, including ESR1 mutation (DNA mutations), Oligonucleotide analogs and antagonists (miRNA), poly (ADP-ribose) polymerase (PARP) in BRCA mutations, CDK4/6 (cell cycle regulating factor initiates tumor progression), Androgen receptor (a steroid hormone receptor), that have entered clinical validation procedure. In this review, we summarize the role of novel breast cancer diagnostic biomarkers, drug resistance, and therapeutic implications for breast cancer.

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.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Nrf2 improves hippocampal synaptic plasticity, learning and memory through the circ-Vps41/miR-26a-5p/CaMKIV regulatory network

Antioxidant response transcription factor nuclear factor erythroid-2-related factor 2 (Nrf2/Nfe2l2) is a neuroprotective agent in learning and memory impairment. This study provides a new perspective to explore the regulatory mechanisms of Nrf2. Here, we found that Nrf2 regulated circular RNA circ-Vps41 to increase hippocampal synaptic plasticity; Nrf2 bound the Vps41 promoter to activate transcription of the Vps41 gene and promote expression of circ-Vps41; circ-Vps41 positively correlated with Nrf2, synaptic plasticity, and learning and memory but negatively correlated with reactive oxygen species; and Nrf2 promoted CaMKIV expression by increasing levels of circ-Vps41, which can absorb miR-26a-5p that targets CaMKIV. Our findings revealed a new circRNA-based regulatory network regulated by Nrf2 and provided novel insights into the potential mechanism involved in the improvement of learning and memory impairment.

Circular RNA as a Potential Biomarker for Forensic Age Prediction

In forensic science, accurate estimation of the age of a victim or suspect can facilitate the investigators to narrow a search and aid in solving a crime. Aging is a complex process associated with various molecular regulations on DNA or RNA levels. Recent studies have shown that circular RNAs (circRNAs) upregulate globally during aging in multiple organisms such as mice and C.elegans because of their ability to resist degradation by exoribonucleases. In the current study, we attempted to investigate circRNAs’ potential capability of age prediction. Here, we identified more than 40,000 circRNAs in the blood of thirteen Chinese unrelated healthy individuals with ages of 20–62 years according to their circRNA-seq profiles. Three methods were applied to select age-related circRNA candidates including the false discovery rate, lasso regression, and support vector machine. The analysis uncovered a strong bias for circRNA upregulation during aging in human blood. A total of 28 circRNAs were chosen for further validation in 30 healthy unrelated subjects by RT-qPCR, and finally, 5 age-related circRNAs were chosen for final age prediction models using 100 samples of 19–73 years old. Several different algorithms including multivariate linear regression (MLR), regression tree, bagging regression, random forest regression (RFR), and support vector regression (SVR) were compared based on root mean square error (RMSE) and mean average error (MAE) values. Among five modeling methods, regression tree and RFR performed better than the others with MAE values of 8.767 years (S.rho = 0.6983) and 9.126 years (S.rho = 0.660), respectively. Sex effect analysis showed age prediction models significantly yielded smaller prediction MAE values for males than females (MAE = 6.133 years for males, while 10.923 years for females in the regression tree model). In the current study, we first used circRNAs as additional novel age-related biomarkers for developing forensic age estimation models. We propose that the use of circRNAs to obtain additional clues for forensic investigations and serve as aging indicators for age prediction would become a promising field of interest.

Circular RNA Cwc27 contributes to Alzheimer's disease pathogenesis by repressing Pur-α activity

Circular RNAs (circRNAs) have gained growing attention in participating in various biological processes and referring to multiply kinds of diseases. Although differentially expressed circRNA profiling in Alzheimer's disease (AD) has been established, little is known about the precise characteristic and functions of key circRNAs with direct relevance to AD in gene expression and disease-related cognition. Herein, we screened and identified circCwc27 as a novel circRNA implicated in AD. CircCwc27 was a neuronal-enriched circRNA that abundantly expressed in the brain and significantly upregulated in AD mice and patients. Knockdown of circCwc27 markedly improved AD-related pathological traits and ameliorated cognitive dysfunctions. Mechanistically, we excluded the miRNA decoy mechanism and focused on the important function of circRNA-RNA-binding protein (RBP) interaction in AD. CircCwc27 directly bound to purine-rich element-binding protein A (Pur-α), increased retention of cytoplasmic Pur-α, and suppressed Pur-α recruitment to the promoters of a cluster of AD genes, including amyloid precursor protein (APP), dopamine receptor D1 (Drd1), protein phosphatase 1, regulatory inhibitor subunit1B (Ppp1r1b), neurotrophic tyrosine kinase, receptor, type 1 (Ntrk1), and LIM homeobox 8 (Lhx8). Downregulation of circCwc27 enhanced the affinity of Pur-α binding to these promoters, leading to altered transcription of Pur-α targets. Moreover, Pur-α overexpression largely phenocopied circCwc27 knockdown in preventing Aβ deposition and cognitive decline. Together, our findings suggest significant functional consequences of a circRNA-protein interaction, that circCwc27, by associating with the regulatory protein Pur-α, may act as a crucial player in AD pathogenesis and represent a promising AD therapeutic target with clinical translational potential.

Circular RNA detection methods: A minireview

Circular RNA (circRNA), a novel type of covalently closed RNA, is implicated in several developmental and metabolic disease processes. CircRNAs exhibit tissue-specific expression, and are stable, abundant, and highly conserved, making them ideal biomarkers for diagnosis and prognosis. Accurate profiling of circRNA, however, is a prerequisite for their clinical application. Traditional methods such as northern blotting, RT-qPCR, and microarray analysis provide useful but limited information. To address these issues, a number of novel assays have recently emerged, such as droplet digital PCR (ddPCR), isothermal exponential amplification, and rolling cycle amplification, which increase the sensitivity and specificity of circRNA detection. Herein, we summarize the advantages and limitations of the new detection methods and discuss the challenges as well as future directions.

Circular RNAs exhibit limited evidence for translation, or translation regulation of the mRNA counterpart in terminal hematopoiesis

Each day, about 10¹² erythrocytes and platelets are released into the bloodstream. This substantial output from hematopoietic stem cells is tightly regulated by transcriptional and epigenetic factors. Whether and how circular RNAs (circRNAs) contribute to the differentiation and/or identity of hematopoietic cells is to date not known. We recently reported that erythrocytes and platelets contain the highest levels and numbers of circRNAs among hematopoietic cells. Here, we provide the first detailed analysis of circRNA expression during erythroid and megakaryoid differentiation. CircRNA expression not only significantly increased upon enucleation, but also had limited overlap between progenitor cells and mature cells, suggesting that circRNA expression stems from regulated processes rather than resulting from mere accumulation. To study circRNA function in hematopoiesis, we first compared the expression levels of circRNAs with the translation efficiency of their mRNA counterpart. We found that only one out of 2531 (0.04%) circRNAs associated with mRNA-translation regulation. Furthermore, irrespective of thousands of identified putative open reading frames, deep ribosome-footprinting sequencing, and mass spectrometry analysis provided little evidence for translation of endogenously expressed circRNAs. In conclusion, circRNAs alter their expression profile during terminal hematopoietic differentiation, yet their contribution to regulate cellular processes remains enigmatic.

Circular RNAs in stem cells: from basic research to clinical implications

Circular RNAs (circRNAs) are a special class of endogenous RNAs with a wide variety of pathophysiological functions via diverse mechanisms, including transcription, microRNA (miRNA) sponge, protein sponge/decoy, and translation. Stem cells are pluripotent cells with unique properties of self-renewal and differentiation. Dysregulated circRNAs identified in various stem cell types can affect stem cell self-renewal and differentiation potential by manipulating stemness. However, the emerging roles of circRNAs in stem cells remain largely unknown. This review summarizes the major functions and mechanisms of action of circRNAs in stem cell biology and disease progression. We also highlight circRNA-mediated common pathways in diverse stem cell types and discuss their diagnostic significance with respect to stem cell-based therapy.

Expression Profile Analysis to Identify Circular RNA Expression Signatures in the Prolificacy Trait of Yunshang Black Goat Pituitary in the Estrus Cycle

The pituitary gland is an important organ. It is a complex area of the brain involved in endocrine function and reproductive regulation. However, the function of the pituitary in goat reproduction is still unclear. Herein, RNA sequencing was used to explore the expression patterns of circle RNAs (circRNAs) in the pituitary of Yunshang black goats during the various estrus phases. Then the host genes of the circRNAs were predicted, and a competing endogenous RNA (ceRNA) network was constructed. The results showed a total of 6,705 circRNAs in the pituitary of Yunshang black goats, among which 388 differentially expressed (DE) circRNAs (214 were upregulated, while 174 were downregulated) were identified between high- and low-yield Yunshang black goats in the follicular phase (HF vs. LF); moreover, 361 DE circRNAs (136 were upregulated, while 225 were downregulated) were identified between high- and low-yield Yunshang black goats in the luteal phase (HL vs. LL). There were 65 DE circRNAs targeting 40 miRNAs in the HF vs. LF comparison and 46 DE circRNAs targeting 31 miRNAs in the HL vs. LL comparison. We identified chi_circ_0030920, chi_circ_0043017, chi_circ_0008353, chi_circ_0041580, and chi_circ_0016478 as the key circRNAs through functional enrichment analysis. The ceRNA network analysis showed that chi_circ_0031209 and chi_circ_0019448 might play an important role in reproduction by influencing the expression of prolactin receptor (PRLR) in high- and low-yielding goats during the luteal phase, whereas chi_circ_0014542 regulates the expression of WNT5A during the follicular phase. Our study provided the overall expression profiles of circRNAs in the goat pituitary during the estrus phase, which provides new insight into the mechanism of high-yield goats, which can be helpful to guide goat breeding.

Digital multiplexed analysis of circular RNAs in FFPE and fresh non‐small cell lung cancer specimens

Although many studies highlight the implication of circular RNAs (circRNAs) in carcinogenesis and tumor progression, their potential as cancer biomarkers has not yet been fully explored in the clinic due to the limitations of current quantification methods. Here, we report the use of the nCounter platform as a valid technology for the analysis of circRNA expression patterns in non‐small cell lung cancer (NSCLC) specimens. Under this context, our custom‐made circRNA panel was able to detect circRNA expression both in NSCLC cells and formalin‐fixed paraffin‐embedded (FFPE) tissues. CircFUT8 was overexpressed in NSCLC, contrasting with circEPB41L2, circBNC2, and circSOX13 downregulation even at the early stages of the disease. Machine learning (ML) approaches from different paradigms allowed discrimination of NSCLC from nontumor controls (NTCs) with an 8‐circRNA signature. An additional 4‐circRNA signature was able to classify early‐stage NSCLC samples from NTC, reaching a maximum area under the ROC curve (AUC) of 0.981. Our results not only present two circRNA signatures with diagnosis potential but also introduce nCounter processing following ML as a feasible protocol for the study and development of circRNA signatures for NSCLC.

Emergent Roles of Circular RNAs in Metabolism and Metabolic Disorders

Circular RNAs (circRNAs) are an emerging group of long non-coding RNAs (lncRNAs) and have attracted attention again according to the progress in high-throughput sequencing in recent years. circRNAs are genome transcripts produced from pre-messenger (m)RNA regions in a specific process called “back-splicing,” which forms covalently closed continuous loops. Due to their lack of a 5’ cap and 3’ poly-adenylated tails, circRNAs are remarkably more stable than linear RNAs. Functionally, circRNAs can endogenously sponge to microRNAs, interact with RNA-binding proteins (RBPs), or translate themselves. Moreover, circRNAs can be expressed in cell type- or tissue-specific expression patterns. Therefore, they are proposed to play essential roles in fine-tuning our body’s homeostasis by regulating transcription and translation processes. Indeed, there has been accumulating emergent evidence showing that dysregulation of circRNAs can lead to metabolic disorders. This study explored the current knowledge of circRNAs that regulate molecular processes associated with glucose and lipid homeostasis and related pathogeneses of metabolic disorders. We also suggest the potential role of circRNAs as disease biomarkers and therapeutic targets.

Genome-wide perturbations of Alu expression and Alu-associated post-transcriptional regulations distinguish oligodendroglioma from other gliomas

Alu is a primate-specific repeat element in the human genome and has been increasingly appreciated as a regulatory element in many biological processes. But the appreciation of Alu has been limited in tumorigenesis, especially for brain tumor. To investigate the relevance of Alu to the gliomagenesis, we studied Alu element-associated post-transcriptional processes and the RNA expression of the element by performing RNA-seq for a total of 41 pairs of neurotypical and diverse glioma brain tissues. We find that A-to-I editing and circular RNA levels, as well as Alu RNA expression, are decreased overall in gliomas, compared to normal tissue. Interestingly, grade 2 oligodendrogliomas are least affected in A-to-I editing and circular RNA levels among gliomas, whereas they have a higher proportion of down-regulated Alu subfamilies, compared to the other gliomas. These findings collectively imply a unique pattern of Alu-associated transcriptomes in grade 2 oligodendroglioma, providing an insight to gliomagenesis from the perspective of an evolutionary genetic element.

Competing endogenous RNA network mediated by circ_3205 in SARS-CoV-2 infected cells

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a new member of the Betacoronaviridae family, responsible for the recent pandemic outbreak of COVID-19. To start exploring the molecular events that follow host cell infection, we queried VirusCircBase and identified a circular RNA (circRNA) predicted to be synthesized by SARS-CoV-2, circ_3205, which we used to probe: (i) a training cohort comprised of two pools of cells from three nasopharyngeal swabs of SARS-CoV-2 infected (positive) or uninfected (negative, UCs) individuals; (ii) a validation cohort made up of 12 positive and 3 negative samples. The expression of circRNAs, miRNAs and miRNA targets was assayed through real-time PCR. CircRNA-miRNA interactions were predicted by TarpMiR, Analysis of Common Targets for circular RNAs (ACT), and STarMir tools. Enrichment of the biological processes and the list of predicted miRNA targets were retrieved from DIANA miRPath v3.0. Our results showed that the predicted SARS-CoV-2 circ_3205 was expressed only in positive samples and its amount positively correlated with that of SARS-CoV-2 Spike (S) mRNA and the viral load (r values = 0.80952 and 0.84867, Spearman's correlation test, respectively). Human (hsa) miR-298 was predicted to interact with circ_3205 by all three predictive tools. KCNMB4 and PRKCE were predicted as hsa-miR-298 targets. Interestingly, the function of both is correlated with blood coagulation and immune response. KCNMB4 and PRKCE mRNAs were upregulated in positive samples as compared to UCs (6 and 8.1-fold, p values = 0.049 and 0.02, Student's t test, respectively) and their expression positively correlated with that of circ_3205 (r values = 0.6 and 0.25, Spearman's correlation test, respectively). We propose that our results convincingly suggest that circ_3205 is a circRNA synthesized by SARS-CoV-2 upon host cell infection and that it may behave as a competitive endogenous RNA (ceRNA), sponging hsa-miR-298 and contributing to the upregulation of KCNMB4 and PRKCE mRNAs.

Knockout of circRNAs by base editing back-splice sites of circularized exons

Many circular RNAs (circRNAs) are produced from back-splicing of exons of precursor mRNAs and are generally co-expressed with cognate linear RNAs. Methods for circRNA-specific knockout are lacking, largely due to sequence overlaps between forms. Here, we use base editors (BEs) for circRNA depletion. By targeting splice sites involved in both back-splicing and canonical splicing, BEs can repress circular and linear RNAs. Targeting sites predominantly for circRNA biogenesis, BEs could efficiently repress the production of circular but not linear RNAs. As hundreds of exons are predominantly back-spliced to produce circRNAs, this provides an efficient method to deplete circRNAs for functional study.

Genome-wide identification and functional analysis of circRNAs in Trichophyton rubrum conidial and mycelial stages

Background

Circular RNAs (circRNAs) are a group of noncoding RNAs that participate in gene expression regulation in various pathways. The essential roles of circRNAs have been revealed in many species. However, knowledge of circRNAs in fungi is still not comprehensive.

Results

Trichophyton rubrum (T. rubrum) is considered a model organism of human pathogenic filamentous fungi and dermatophytes. In this study, we performed a genome-wide investigation of circRNAs in T. rubrum based on high-throughput sequencing and ultimately identified 4254 circRNAs. Most of these circRNAs were specific to the conidial or mycelial stage, revealing a developmental stage-specific expression pattern. In addition, 940 circRNAs were significantly differentially expressed between the conidial and mycelial stages. PCR experiments conducted on seven randomly selected differentially expressed (DE-) circRNAs confirmed the circularized structures and relative expression levels of these circRNAs. Based on their genome locations, most circRNAs originated from intergenic regions, unlike those in plants and animals. Furthermore, we constructed circRNA-miRNA-mRNA regulatory networks that included 661 DE-circRNAs targeting 140 miRNAs and further regulating 2753 mRNAs. The relative expression levels of two randomly selected circRNA-miRNA-mRNA axes were investigated by qRT-PCR, and the competing endogenous RNA (ceRNA) network theory was validated. Functional enrichment analysis of the target genes suggested that they were significantly involved in posttranscriptional processes and protein synthesis as well as some small-molecule metabolism processes. CircRNAs are relatively more conserved in closely related dermatophytes but rarely conserved in distantly related species. Tru_circ07138_001 is a highly conserved circRNA that was conserved in all ten dermatophytes analyzed in our study and three distantly related species. Its host gene TERG_07138 was also highly conserved in two of these distantly related species Gallus gallus and Caenorhabditis elegans. The specific role of this circRNA deserves further exploration.

Conclusions

Our study is the first to provide a global profile of circRNAs in T. rubrum as well as dermatophytes. These results could serve as valuable resources for research on circRNA regulatory mechanisms in fungi and reveal new insights for further investigation of the physical characteristics of these significant human fungal pathogens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08184-y.

Circular RNAs in Diabetic Nephropathy: Updates and Perspectives

Circular RNAs (circRNAs) are widespread endogenous transcripts lacking 5′-caps and 3′-polyadenylation tails. Their closed-loop structure confers exonuclease resistance and extreme stability. CircRNAs play essential roles in various diseases, including diabetes. Diabetic nephropathy (DN) is the leading cause of end-stage kidney disease and is one of the most common complications of diabetes. CircRNAs are key in DN and therefore important for understanding DN pathophysiology and developing new therapeutic strategies. In the present review, we briefly introduce the characteristics and functions of circRNAs and summarize recent discoveries on how circRNAs participate in DN. Based on these advances, we suggest future perspectives for studying circRNAs in DN to improve DN treatment and management.

Recent origination of circular RNAs in plants

Circular RNA (circRNA) is a kind of new regulatory RNA with diverse biological functions. Numerous circRNAs have been identified in many plant species; however, evolution of plant circRNAs remains largely unknown. In this study, we assembled full-length sequences of 6519 rice (Oryza sativa) circRNAs and analyzed their conservation in another 46 plant species based on comparison of sequences and expression patterns. We found that, at the genomic level, 8.7% of the 6519 circRNAs were conserved in dicotyledonous plants and 49.1% in Oryza genus. Meanwhile, 57.8% of parental protein-coding genes of the rice circRNAs originated recently after divergence of monocotyledonous plants, implying recent origin of the majority of rice circRNAs, a conclusion further supported by the results based on analysis of 4663 full-length circRNAs in Arabidopsis thaliana. Accordingly, we proposed three models to address the origination of different types of circRNAs. Taken together, the results obtained in this study provide new insights for the evolutionary dynamics of plant circRNAs and candidate circRNAs for further functional exploration.

Genetic tools for the stable overexpression of circular RNAs

Circular RNAs (circRNAs) are a class of non-coding RNAs featuring a covalently closed ring structure formed through backsplicing. circRNAs are broadly expressed and contribute to biological processes through a variety of functions. Standard gain-of-function and loss-of-function approaches to study gene functions have significant limitations when studying circRNAs. Overexpression studies in particular suffer from the lack of efficient genetic tools. While mammalian expression plasmids enable transient circRNA overexpression in cultured cells, most cell biological studies require long-term ectopic expression. Here we report the development and characterization of genetic tools enabling stable circRNA overexpression in vitro and in vivo. We demonstrated that circRNA expression constructs can be delivered to cultured cells via transposons, whereas lentiviral vectors have limited utility for the delivery of circRNA constructs due to viral RNA splicing in virus-producing cells. We further demonstrated ectopic circRNA expression in a hepatocellular carcinoma mouse model upon circRNA transposon delivery via hydrodynamic tail vein injection. Furthermore, we generated genetically engineered mice harbouring circRNA expression constructs. We demonstrated that this approach enables constitutive, global circRNA overexpression as well as inducible circRNA expression directed specifically to melanocytes in a melanoma mouse model. These tools expand the genetic toolkit available for the functional characterization of circRNAs.

Circ_0000745 promotes acute lymphoblastic leukemia progression through mediating miR-494-3p/NET1 axis

BACKGROUND

Circular RNAs (circRNAs) have shown important regulatory roles in tumorigenesis. However, the role and working mechanism of circ_0000745 in acute lymphoblastic leukemia (ALL) development remain largely unclear.

METHODS

The expression of circ_0000745, sperm antigen with calponin homology and coiled-coil domains 1 (SPECC1), microRNA-494-3p (miR-494-3p), and neuroepithelial cell transforming 1 (NET1) messenger RNA (mRNA) and protein was analyzed by reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and Western blot assay. Flow cytometry was performed to assess cell apoptosis and cell cycle progression. Extracellular acidification rate (ECAR) was assessed to analyze cell glycolysis. Cell viability was analyzed by 3-(4,5-Dimethylthiazol-2-yl)-2,5-Diphenyltetrazolium Bromide (MTT) assay. Ferroptosis was assessed through measuring the intracellular levels of iron and lipid reactive oxygen species (ROS). Dual-luciferase reporter assay and RNA immunoprecipitation (RIP) assay were conducted to validate the interaction between miR-494-3p and circ_0000745 or NET1.

RESULTS

Circ_0000745 expression was elevated in ALL patients and cell lines. Circ_0000745 knockdown restrained cell cycle progression and glycolysis and triggered cell apoptosis and ferroptosis. Circ_0000745 acted as a molecular sponge for miR-494-3p in ALL cells. miR-494-3p silencing partly diminished circ_0000745 knockdown-induced anti-tumor effects in ALL cells. NET1 was a target of miR-494-3p, and miR-494-3p overexpression-induced anti-tumor influences were partly counteracted by the accumulation of NET1 in ALL cells. Circ_0000745 can positively regulate NET1 expression by sponging miR-494-3p in ALL cells.

CONCLUSION

Circ_0000745 contributed to ALL development partly by binding to miR-494-3p to induce NET1 expression.0020.

The Cellular and Viral circRNAome Induced by Respiratory Syncytial Virus Infection

Circular RNAs (circRNAs) are a new class of noncoding RNAs that have gained increased attention. DNA virus infections have been reported to induce modifications in cellular circRNA transcriptomes and express viral circRNAs. However, the identification and expression of cellular and viral circRNAs are unknown in the context of respiratory syncytial virus (RSV), a human RNA virus with no effective treatments or vaccines. Here, we report a comprehensive identification of the cellular and viral circRNAs induced by RSV infection in A549 cells with high-throughput sequencing. In total, 53,719 cellular circRNAs and 2,280 differentially expressed cellular circRNAs were identified. Trend analysis further identified three significant expression pattern clusters, which were related to the antiviral immune response according to gene enrichment analysis. Subsequent results showed that not only RSV infection but also poly(I·C) treatment and another RNA virus infection induced the upregulation of the top 10 circRNAs from the focused cluster. The top 10 circRNAs generally inhibit RSV replication in turn. Moreover, 1,254 viral circRNAs were identified by the same circRNA sequencing. The induced expression of viral circRNAs by RSV infection was found not only in A549 cells but also in HEp-2 cells. Additionally, we profiled the general characteristics of both cellular and viral circRNAs such as back-splicing signals, etc. Collectively, RSV infection induced the differential expression of cellular circRNAs, some of which affected RSV infection, and RSV also expressed viral circRNAs. Our study reveals novel layers of host-RSV interactions and identifies cellular or viral circRNAs that may be novel therapeutic targets or biomarkers. IMPORTANCE Noncoding RNAs (ncRNAs) demonstrate substantial roles in cell-virus interactions. Circular RNAs (circRNAs) are a newly identified class of ncRNAs that have gained increased attention recently. DNA virus infections have been reported to induce modifications in cellular circRNA transcriptomes and express viral circRNAs. However, the identification and expression of cellular and viral circRNAs are unknown in the context of respiratory syncytial virus (RSV), a human RNA virus with no effective treatments or vaccines. Here, we report a comprehensive identification of the cellular and viral circRNAs induced by RSV infection by high-throughput sequencing. We revealed that RSV infection induces the differential expression of cellular circRNAs, some of which affected RSV infection, and that RSV also expresses viral circRNAs. Our study reveals novel layers of host-RSV interactions and identifies cellular or viral circRNAs that may be novel therapeutic targets or biomarkers.

Circ-Hdgfrp3 shuttles along neurites and is trapped in aggregates formed by ALS-associated mutant FUS

CircRNAs belong to a family of RNA molecules which are conserved in evolution, have tissue-specific expression, and are abundant in neuronal cells. Here, we define several features of circ-Hdgfrp3 and describe interesting alterations occurring in motor neurons (MNs) carrying ALS-associated FUS mutations. Through a highly sensitive in situ approach we describe that circ-Hdgfrp3 traffics along neurites, while upon oxidative stress it is retained in the perinuclear region. While in wild-type stressed MNs, circ-Hdgfrp3 localizes in stress granules (SGs), in MNs carrying mutant FUS, a higher proportion of circ-Hdgfrp3 was trapped into cytoplasmic aggregates. Upon stress removal, circ-Hdgfrp3 was easily freed from SGs whereas it was less efficiently released from FUS-aggregates. We found that the human circ-Hdgfrp3 counterpart was also similarly associated to mutant FUS-aggregates in stressed neuronal cells. Overall, the alteration of circ-Hdgfrp3 trafficking adds a further layer of complexity to the role of FUS-aggregates in ALS disease.

CircularLRRC7 is a Potential Tumor Suppressor Associated With miR-1281 and PDXP Expression in Glioblastoma

Circular RNAs (circRNAs) are usually enriched in neural tissues, yet about 80% circRNAs have lower expression in gliomas relative to normal brains, highlighting the importance of circRNAs as tumor suppressors. However, the clinical impact as well as the pathways regulated by the tumor-suppressive circRNAs remain largely unknown in glioblastoma (GBM). Through bioinformatic analysis followed by experimental validation, we found that hsa_circ_0114014 (circLRRC7) was dramatically down-regulated in GBM when compared with normal brain tissues (p < 0.0001). GBM patients with a lower circLRRC7 expression had poorer progression-free survival (PFS, p < 0.05) and overall survival (OS, p < 0.05). Analyses of the predicted target miRNAs of circLRRC7 in CSCD and CRI databases, in combination with the miRNA expression data in GBMs and normal brains from GSE database, revealed miR-1281 as a potential downstream target of circLRRC7. Subsequently, the target genes of hsa-mir-1281 were predicted by TargetScan, miRDB and miRNATAR databases. Intersection analysis and correlation test indicated that PDXP was a potential target of miR-1281. In summary, circLRRC7 may be a tumor suppressor that associated with miR-1281 and PDXP expression in GBM, which may provide novel therapeutic targets for GBM treatment.

Genome-Wide Identification of CircRNAs of Infective Larvae and Adult Worms of Parasitic Nematode, Haemonchus contortus

CircRNAs, a novel class of ncRNA family, are endogenous transcriptional products involved in various biological and physiological processes in plants and animals. However, almost no information is available for circRNAs of parasitic helminths. In the present study, the circRNAs repertoire was comprehensively explored in Haemonchus contortus, a blood-sucking parasitic nematode of ruminants. In total, 20073 circRNAs were identified and annotated from three key developmental stages/genders of H. contortus including the free-living infective third-stage larvae (L3, 18883), parasitic adult female (Af, 3491), and male worms (Am, 2550) via deep-sequencing technology and bioinformatic analysis. Among these identified circRNAs, 71% were derived from exonic regions of protein-coding genes. The number of circRNAs transcribed from the X chromosome (4704) was higher than that from Chromosome I-V (3143, 3273, 3041, 3030, 2882). The amount of highly expressed circRNAs in third-stage larvae was significantly more abundant than that in adult stage. 15948 and 16847 circRNAs were differentially expressed between Af and L3s and between Am and L3, respectively. Among them, 13409 circRNAs existed in both comparisons. Furthermore, 1119 circRNAs were differentially expressed between Af_and_Am. GO enrichment analysis indicated that source genes of circRNAs differentially expressed between Am and L3 as well as between Af and L3 were significantly enriched in many biological processes, primarily including signaling, signal transduction and cell communication terms. KEGG analysis revealed that parental genes of differentially expressed circRNAs were mainly related to metabolism (pyruvate metabolism, glycerophospholipid metabolism, and carbon metabolism), MAPK signaling pathway, and phosphatidylinositol signaling system. Moreover, many circRNAs contained one or more miRNA potential binding sites, suggesting that they could regulate gene expression at the post-transcriptional level. Furthermore, the correctness of head-to-tail back splicing site and alternative circularization events were verified by Sanger sequencing using both divergent and convergent primers. Finally, the reliability of RNA-Seq data and the resistance of circRNAs to RNase R digestion were confirmed by quantitative RT-PCR. Taken together, our findings provide a foundation for elucidating the regulatory mechanisms of circRNAs in H. contortus, which will advance the understanding of circRNAs in parasitic nematodes.

CircRNA75 and CircRNA72 Function as the Sponge of MicroRNA-200 to Suppress Coelomocyte Apoptosis Via Targeting Tollip in Apostichopus japonicus

Circular RNAs (circRNAs) act as essential regulators in many biological processes, especially in mammalian immune response. Nonetheless, the functions and mechanisms of circRNAs in the invertebrate immune system are largely unclarified. In our previous work, 261 differentially expressed circRNAs potentially related to the development of Apostichopus japonicus skin ulceration syndrome (SUS), which is a major problem restricting the sea cucumber breeding industry, were identified by genome-wide screening. In this study, via miRanda analysis, both circRNA75 and circrRNA72 were shown to share the miR-200 binding site, a key microRNA in the SUS. The two circRNAs were verified to be increased significantly in LPS-exposed primary coelomocytes, similar to the results of circRNA-seq in sea cucumber under Vibrio splendidus-challenged conditions. A dual-luciferase assay indicated that both circRNA75 and circRNA72 could bind miR-200 in vivo, in which circRNA75 had four binding sites of miR-200 and only one for circRNA72. Furthermore, we found that miR-200 could bind the 3’-UTR of Toll interacting protein (Tollip) to negatively mediate the expression of Tollip. Silencing Tollip increased primary coelomocyte apoptosis. Consistently, inference of circRNA75 and circRNA72 could also downregulate Tollip expression, thereby increasing the apoptosis of primary coelomocytes, which could be blocked by miR-200 inhibitor treatment. Moreover, the rate of si-circRNA75-downregulated Tollip expression was higher than that of si-circRNA72 under an equivalent amount. CircRNA75 and circRNA72 suppressed coelomocyte apoptosis by sponging miR-200 to promote Tollip expression. The ability of circRNA to adsorb miRNA might be positively related to the number of binding sites for miRNA.

Circular RNAs in physiology and non-immunological diseases

Circular RNAs (circRNAs) are covalently closed single-stranded RNAs. Four subclasses of circRNAs have been identified in animal cells, and they have unique features in their biogenesis, degradation, and transport. CircRNAs have diverse molecular functions in sponging miRNAs, regulating transcription, modulating RNA-binding proteins, and even encoding proteins. Some circRNAs are important regulators of various physiological processes to maintain homeostasis. Dysregulation of circRNAs is associated with human disorders, and individual circRNAs are crucial factors that contribute to major diseases including non-immunological diseases such as cancers, neurological disorders, cardiovascular disease, and metabolic disease. Debates on circRNAs have also been raised in recent years, and further studies would help to resolve these disputes and potentially lead to biomedical applications of circRNAs.

The potential of using blood circular RNA as liquid biopsy biomarker for human diseases

Circular RNA (circRNA) is a novel class of single-stranded RNAs with a closed loop structure. The majority of circRNAs are formed by a back-splicing process in pre-mRNA splicing. Their expression is dynamically regulated and shows spatiotemporal patterns among cell types, tissues and developmental stages. CircRNAs have important biological functions in many physiological processes, and their aberrant expression is implicated in many human diseases. Due to their high stability, circRNAs are becoming promising biomarkers in many human diseases, such as cardiovascular diseases, autoimmune diseases and human cancers. In this review, we focus on the translational potential of using human blood circRNAs as liquid biopsy biomarkers for human diseases. We highlight their abundant expression, essential biological functions and significant correlations to human diseases in various components of peripheral blood, including whole blood, blood cells and extracellular vesicles. In addition, we summarize the current knowledge of blood circRNA biomarkers for disease diagnosis or prognosis.

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.

CircARL8B Contributes to the Development of Breast Cancer Via Regulating miR-653-5p/HMGA2 Axis

Circular RNAs (circRNAs) act as essential regulators in breast cancer (BC) progression. In this paper, we aimed to investigate the functions of circARL8B in BC. The levels of circARL8B, ADP Ribosylation Factor Like GTPase 8B (ARL8B), miR-653-5p and high-mobility group AT-hook 2 (HMGA2) mRNA were examined by qRT-PCR. The stability of circARL8B was determined by RNase R assay and Actinomycin D assay. Cell viability and metastasis were evaluated by Cell Counting Kit-8 (CCK-8) assay and transwell assay, respectively. The levels of cellular phospholipids and triglycerides were measured using relevant kits. Protein levels were measured by western blot analysis. The association between miR-653-5p and circARL8B or HMGA2 was verified by dual-luciferase reporter assay. A murine xenograft model was established to explore the function of circARL8B in vivo. CircARL8B was increased in BC tissues and cells. CircARL8B silencing inhibited cell viability, migration, invasion and fatty acid metabolism in BC cells in vitro and blocked tumor growth in vivo. MiR-653-5p was identified as the target of circARL8B and miR-653-5p was negatively modulated by circARL8B. The suppressive role of circARL8B silencing in BC cell progression was abolished by miR-653-5p downregulation. Moreover, HMGA2 was the target gene of miR-653-5p. HMGA2 overexpression abrogated the effect of miR-653-5p on BC cell development. In addition, circARL8B knockdown might block PGE2/PI3K/AKT/GSK-3β/Wnt/β-catenin pathway. Silencing of circARL8B inhibited cell viability, migration, invasion and fatty acid metabolism via miR-653-5p/HMGA2 axis in BC.

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.

The design and synthesis of circular RNAs

Circular RNAs (circRNAs) are a novel class of RNAs distinguished by their single-stranded, covalently-closed topology. Although initially perceived as rare byproducts of aberrant splicing, circRNAs are now recognized as ubiquitously expressed and functionally significant. These discoveries have led to a growing need for ways to model circRNAs in living cells to advance our understanding of their biogenesis, regulation, and function, and to adopt them as new technologies for application within research and medicine. In this review, we provide an updated summary of approaches used to produce circRNAs in vitro and in vivo, the latter of which has grown considerably in recent years. Given increased interest in the unique functions carried out by individual circRNAs, we further dedicate a section on how to customize synthesized circRNAs for specific biological roles. We focus on the most common applications, including designing circRNAs for protein delivery, to target miRNAs and proteins, to act as fluorescent reporters, and to modulate cellular immunity.

Circular RNA: Biosynthesis in vitro

Circular RNA (circRNA) is a unique type of noncoding RNA molecule. Compared with traditional linear RNA, circRNA is a covalently closed circle produced by a process called backsplicing. CircRNA is abundant in many cells and has rich functions in cells, such as acting as miRNA sponge, protein sponge, protein scaffold, and mRNA regulator. With the continuous development of circRNA study, circRNA has also played an important role in medical applications, including circRNA vaccines and gene therapy. In this review, we illustrate the synthesis of circRNAs in vitro. We focus on biological ligation methods, such as enzymatic ligation from the bacteriophage T4 and ribozyme method. In addition, we summarize the current challenges in the design, synthesis, application, and production of circRNAs, and propose possible solutions in the future. CircRNA is expected to play an essential role in basic research and medical applications.

Identification of Potential circRNA-microRNA-mRNA Regulatory Network in Skeletal Muscle

Circular RNAs (circRNAs) are a newly discovered family of regulatory RNAs generated through backsplicing. Genome-wide profiling of circRNAs found that circRNAs are ubiquitously expressed and regulate gene expression by acting as a sponge for RNA-binding proteins (RBPs) and microRNAs (miRNAs). To identify circRNAs expressed in mouse skeletal muscle, we performed high-throughput RNA-sequencing of circRNA-enriched gastrocnemius muscle RNA samples, which identified more than 1,200 circRNAs. In addition, we have identified more than 14,000 and 15,000 circRNAs in aging human skeletal muscle tissue and satellite cells, respectively. A subset of abundant circRNAs was analyzed by RT-PCR, Sanger sequencing, and RNase R digestion assays to validate their expression in mouse skeletal muscle tissues. Analysis of the circRNA-miRNA-mRNA regulatory network revealed that conserved circNfix might associate with miR-204-5p, a suppressor of myocyte enhancer factor 2c (Mef2c) expression. To support the hypothesis that circNfix might regulate myogenesis by controlling Mef2c expression, silencing circNfix moderately reduced Mef2c mRNA expression and inhibited C2C12 differentiation. We propose that circNfix promotes MEF2C expression during muscle cell differentiation in part by acting as a sponge for miR-204-5p.

Genome-Wide Identification of Circular RNAs Potentially Involved in the Biosynthesis of Secondary Metabolites in Salvia miltiorrhiza

Circular RNAs (circRNAs) play various roles in cellular functions. However, no studies have been reported on the potential involvement of circRNAs in the biosynthesis of secondary metabolites in plants. Here, we performed a genome-wide discovery of circRNAs from root, stem and leaf samples of Salvia miltiorrhiza using RNA-Seq. We predicted a total of 2,476 circRNAs with at least two junction reads using circRNA_finder and CIRI, of which 2,096, 151 and 229 were exonic, intronic and intergenic circRNAs, respectively. Sequence similarity analysis showed that 294 out of 2,476 circRNAs were conserved amongst multiple plants. Of the 55 predicted circRNAs, 31 (56%) were validated successfully by PCR and Sanger sequencing using convergent and divergent primer pairs. Alternative circularisation analysis showed that most parental genes produced two circRNAs. Functional enrichment analyses of the parental genes showed that the primary metabolism pathways were significantly enriched, particularly the carbon metabolism. Differential expression analysis showed that the expression profiles of circRNAs were tissue-specific. Co-expression analysis showed 275 circRNAs, and their parental genes had significantly positive correlations. However, 14 had significantly negative correlations. Weighted gene co-expression network analysis showed that nine circRNAs were co-expressed with four modules of protein-coding genes. Next, we found 416 exonic circRNAs with miRNA-binding sites, suggesting possible interactions between circRNAs and miRNAs. Lastly, we found six validated circRNAs, namely, SMscf2473-46693-46978, SMscf3091-29256-29724, SMscf16-111773-112193, SMscf432-13232-13866, SMscf7007-10563-10888 and SMscf1730-1749-2013, which were originated from the genes involved in the biosynthesis of secondary metabolites. Their parental genes were acetyl-CoA C-acetyltransferase 1 (SmAACT1), 1-deoxy-d-xylulose-5-phosphate synthase 2 (SmDXS2), 4-hydroxy-3-methylbut-2-enyl diphosphate reductase 1 (SmHDR1), kaurene synthase-like 2 (SmKSL2), DWF4 and CYP88A3, respectively. In particular, the correlation coefficient of SMscf2473-46693-46978 and SmDXS2 gene was 0.86 (p = 0.003), indicating a potential interaction between this pair of circRNA and its parent gene. Our results provided the first comprehensive catalogue of circRNAs in S. miltiorrhiza and identified one circRNA that might play important roles in the biosynthesis of secondary metabolites.

Genome-wide identification and characterization of circular RNA m6A modification in pancreatic cancer

BACKGROUND

N⁶-methyladenosine (m⁶A) is the most abundant modification of RNA in eukaryotic cells and play critical roles in cancer. While most related studies focus on m⁶A modifications in linear RNA, transcriptome-wide profiling and exploration of m⁶A modification in circular RNAs in cancer is still lacking.

METHODS

For the detection of m⁶A modification in circRNAs, we developed a new bioinformatics tools called Circm6A and applied it to the m⁶A-seq data of 77 tissue samples from 58 individuals with pancreatic ductal adenocarcinoma (PDAC).

RESULTS

Circm6A performs better than the existing circRNA identification tools, which achieved highest F1 score among these tools in the detection of circRNAs with m⁶A modifications. By using Circm6A, we identified a total of 8807 m⁶A-circRNAs from our m⁶A-seq data. The m⁶A-circRNAs tend to be hypermethylated in PDAC tumor tissues compared with normal tissues. The hypermethylated m⁶A-circRNAs were associated with a significant gain of circRNA-mRNA coexpression network, leading to the dysregulation of many important cancer-related pathways. Moreover, we found the cues that hypermethylated m⁶A-circRNAs may promote the circularization and translation of circRNAs.

CONCLUSIONS

These comprehensive findings further bridged the knowledge gaps between m⁶A modification and circRNAs fields by depicting the m⁶A-circRNAs genomic landscape of PDAC patients and revealed the emerging roles played by m⁶A-circRNAs in pancreatic cancer. Circm6A is available at https://github.com/canceromics/circm6a .

The role of noncoding RNAs in Parkinson's disease: biomarkers and associations with pathogenic pathways

The discovery of various noncoding RNAs (ncRNAs) and their biological implications is a growing area in cell biology. Increasing evidence has revealed canonical and noncanonical functions of long and small ncRNAs, including microRNAs, long ncRNAs (lncRNAs), circular RNAs, PIWI-interacting RNAs, and tRNA-derived fragments. These ncRNAs have the ability to regulate gene expression and modify metabolic pathways. Thus, they may have important roles as diagnostic biomarkers or therapeutic targets in various diseases, including neurodegenerative disorders, especially Parkinson's disease. Recently, through diverse sequencing technologies and a wide variety of bioinformatic analytical tools, such as reverse transcriptase quantitative PCR, microarrays, next-generation sequencing and long-read sequencing, numerous ncRNAs have been shown to be associated with neurodegenerative disorders, including Parkinson's disease. In this review article, we will first introduce the biogenesis of different ncRNAs, including microRNAs, PIWI-interacting RNAs, circular RNAs, long noncoding RNAs, and tRNA-derived fragments. The pros and cons of the detection platforms of ncRNAs and the reproducibility of bioinformatic analytical tools will be discussed in the second part. Finally, the recent discovery of numerous PD-associated ncRNAs and their association with the diagnosis and pathophysiology of PD are reviewed, and microRNAs and long ncRNAs that are transported by exosomes in biofluids are particularly emphasized.

Internal Promoters and Their Effects on the Transcription of Operon Genes for Epothilone Production in Myxococcus xanthus

The biosynthetic genes for secondary metabolites are often clustered into giant operons with no transcription terminator before the end. The long transcripts are frangible and the transcription efficiency declines along with the process. Internal promoters might occur in operons to coordinate the transcription of individual genes, but their effects on the transcription of operon genes and the yield of metabolites have been less investigated. Epothilones are a kind of antitumor polyketides synthesized by seven multifunctional enzymes encoded by a 56-kb operon. In this study, we identified multiple internal promoters in the epothilone operon. We performed CRISPR-dCas9–mediated transcription activation of internal promoters, combined activation of different promoters, and activation in different epothilone-producing M. xanthus strains. We found that activation of internal promoters in the operon was able to promote the gene transcription, but the activation efficiency was distinct from the activation of separate promoters. The transcription of genes in the operon was influenced by not only the starting promoter but also internal promoters of the operon; internal promoters affected the transcription of the following and neighboring upstream/downstream genes. Multiple interferences between internal promoters thus changed the transcriptional profile of operon genes and the production of epothilones. Better activation efficiency for the gene transcription and the epothilone production was obtained in the low epothilone-producing strains. Our results highlight that interactions between promoters in the operon are critical for the gene transcription and the metabolite production efficiency.

Research Progress on Circular RNA in Glioma

The discovery of circular RNA (circRNA) greatly complements the traditional gene expression theory. CircRNA is a class of non-coding RNA with a stable cyclic structure. They are highly expressed, spatiotemporal-specific and conservative across species. Importantly, circRNA participates in the occurrence of many kinds of tumors and regulates the tumor development. Glioma is featured by limited therapy and grim prognosis. Cancer-associated circRNA compromises original function or creates new effects in glioma, thus contributing to oncogenesis. Therefore, this article reviews the biogenesis, metabolism, functions and properties of circRNA as a novel potential biomarker for gliomas. We elaborate the expression characteristics, interaction between circRNA and other molecules, aiming to identify new targets for early diagnosis and treatment of gliomas.

An autoregulation loop in fust-1 for circular RNA regulation in Caenorhabditis elegans

Many circular RNAs (circRNAs) are differentially expressed in different tissues or cell types, suggestive of specific factors that regulate their biogenesis. Here, taking advantage of available mutation strains of RNA-binding proteins (RBPs) in Caenorhabditis elegans, I performed a screening of circRNA regulation in 13 conserved RBPs. Among them, loss of FUST-1, the homolog of Fused in Sarcoma (FUS), caused downregulation of multiple circRNAs. By rescue experiments, I confirmed FUST-1 as a circRNA regulator. Through RNA sequencing using circRNA-enriched samples, circRNAs targets regulated by FUST-1 were identified globally, with hundreds of them significantly altered. Furthermore, I showed that FUST-1 regulates circRNA formation with only small to little effect on the cognate linear mRNAs. When recognizing circRNA pre-mRNAs, FUST-1 can affect both exon-skipping and circRNA in the same genes. Moreover, I identified an autoregulation loop in fust-1, where FUST-1, isoform a (FUST-1A) promotes the skipping of exon 5 of its own pre-mRNA, which produces FUST-1, isoform b (FUST-1B) with different N-terminal sequences. FUST-1A is the functional isoform in circRNA regulation. Although FUST-1B has the same functional domains as FUST-1A, it cannot regulate either exon-skipping or circRNA formation. This study provided an in vivo investigation of circRNA regulation, which will be helpful to understand the mechanisms that govern circRNA formation.

The role of circular RNAs in neuropathic pain

Neuropathic pain (NP) results from a lesion or disease of the nervous system and accompanied by chronic pain, leading a serious public health issue and economic burden. In terms of the NP mechanisms remaining poorly understood, circular RNAs (circRNAs), owing to their high stability and evolutionary conservation, are expected to be used as potential therapeutic targets or diagnostic biomarkers. In this review, we concisely discuss the characteristics and biological functions of circRNAs. As emerging evidence we reviewed, deregulation of circRNAs (e.g., circ-Filip1l, circHIPK3, ciRS-7, circRNA.2837, circ-Ankib1 and circAnks1a) were involved in NP development. It suggested that specific circRNAs modulated through sponging their target miRNAs and thus regulated the homologous downstream mRNAs and proteins in neuropathic pain. These findings provide a theoretical basis for circRNAs use as biomarkers and therapeutic targets in neuropathic pain.

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.

CircRNA_01477 influences axonal growth via regulating miR-3075/FosB/Stat3 axis

Circular RNAs (circRNAs) are important for the development and regeneration of the nervous system. We investigated the differential expression profiles of circRNA induced by spinal cord injury and reported that circRNA_01477 facilitates spinal astrocyte proliferation and migration after injury in rats. In this study, we further clarified the function and possible mechanism of action of circRNA_01477 in neurons. Fluorescence in situ hybridization assay revealed that circRNA_01477 is mainly localized in the neuronal cytoplasm. Knockdown of circRNA_01477 significantly increased axonal length. The circRNA_01477/microRNAs (miRNA)/messenger RNA (mRNA) interaction network was investigated using RNA sequencing. miRNA-3075 showed a remarkable increase after circRNA_01477 depletion, and either overexpression of miRNA-3075 or downregulation of its target gene FosB significantly promoted axonal growth. Luciferase reporter assay showed that miRNA-3075 could directly bind to the 3'UTR of FosB and negatively regulated FosB transcription. Dual silencing of circRNA_01477 and miR-3075 revealed that miR-3075 inhibition rescued the increased axon length caused by siCircRNA_01477. Finally, we verified that the Stat3 pathway was activated after FosB protein depletion in rat spinal neurons, while the NF-κB pathway was not altered. In summary, our study is the first to report that circRNA_01477 contributes to axon growth by functioning as miRNA sponge by regulating the miRNA-3075/FosB/Stat3 axis.

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.

Circular RNAs as Novel Diagnostic Biomarkers and Therapeutic Targets in Kidney Disease

Circular RNAs (circRNAs) are a novel type of non-coding RNAs that have aroused growing attention in this decade. They are widely expressed in eukaryotes and generally have high stability owing to their special closed-loop structure. Many circRNAs are abundant, evolutionarily conserved, and exhibit cell-type-specific and tissue-specific expression patterns. Mounting evidence suggests that circRNAs have regulatory potency for gene expression by acting as microRNA sponges, interacting with proteins, regulating transcription, or directly undergoing translation. Dysregulated expression of circRNAs were found in many pathological conditions and contribute to the pathogenesis and progression of various disorders, including renal diseases. Recent studies have revealed that circRNAs may serve as novel reliable biomarkers for the diagnosis and prognosis prediction of multiple kidney diseases, such as renal cell carcinoma (RCC), acute kidney injury (AKI), diabetic kidney disease (DKD), and other glomerular diseases. Furthermore, circRNAs expressed by intrinsic kidney cells are shown to play a substantial role in kidney injury, mostly reported in DKD and RCC. Herein, we review the biogenesis and biological functions of circRNAs, and summarize their roles as promising biomarkers and therapeutic targets in common kidney diseases.

Characterizing viral circRNAs and their application in identifying circRNAs in viruses

Circular RNAs (circRNAs) are non-coding RNAs with a special circular structure produced formed by the reverse splicing mechanism, which play an important role in a variety of biological activities. Viruses can encode circRNA, and viral circRNAs have been found in multiple single-stranded and double-stranded viruses. However, the characteristics and functions of viral circRNAs remain unknown. Sequence alignment showed that viral circRNAs are less conserved than circRNAs in animal, indicating that the viral circRNAs may evolve rapidly. Through the analysis of the sequence characteristics of viral circRNAs and circRNAs in animal, it was found that viral circRNAs and animals circRNAs are similar in nucleic acid composition, but have obvious differences in secondary structure and autocorrelation characteristics. Based on these characteristics of viral circRNAs, machine learning algorithms were employed to construct a prediction model to identify viral circRNA. Additionally, analysis of the interaction between viral circRNA and miRNAs showed that viral circRNA is expected to interact with 518 human miRNAs, and preliminary analysis of the role of viral circRNA. And it has been also found that viral circRNAs may be involved in many KEGG pathways related to nervous system and cancer. We curated an online server, and the data and code are available: http://server.malab.cn/viral-CircRNA/.

SRCP: a comprehensive pipeline for accurate annotation and quantification of circRNAs

Here we describe a new integrative approach for accurate annotation and quantification of circRNAs named Short Read circRNA Pipeline (SRCP). Our strategy involves two steps: annotation of validated circRNAs followed by a quantification step. We show that SRCP is more sensitive than other individual pipelines and allows for more comprehensive quantification of a larger number of differentially expressed circRNAs. To facilitate the use of SRCP, we generate a comprehensive collection of validated circRNAs in five different organisms, including humans. We then utilize our approach and identify a subset of circRNAs bound to the miRNA-effector protein AGO2 in human brain samples.

CSCD2: an integrated interactional database of cancer-specific circular RNAs

The significant function of circRNAs in cancer was recognized in recent work, so a well-organized resource is required for characterizing the interactions between circRNAs and other functional molecules (such as microRNA and RNA-binding protein) in cancer. We previously developed cancer-specific circRNA database (CSCD), a comprehensive database for cancer-specific circRNAs, which is widely used in circRNA research. Here, we updated CSCD to CSCD2 (http://geneyun.net/CSCD2 or http://gb.whu.edu.cn/CSCD2), which includes significantly more cancer-specific circRNAs identified from a large number of human cancer and normal tissues/cell lines. CSCD2 contains >1000 samples (825 tissues and 288 cell lines) and identifies a large number of circRNAs: 1 013 461 cancer-specific circRNAs, 1 533 704 circRNAs from only normal samples and 354 422 circRNAs from both cancer and normal samples. In addition, CSCD2 predicts potential miRNA–circRNA and RBP–circRNA interactions using binding motifs from >200 RBPs and 2000 microRNAs. Furthermore, the potential full-length and open reading frame sequence of these circRNAs were also predicted. Collectively, CSCD2 provides a significantly enhanced resource for exploring the function and regulation of circRNAs in cancer.

Circular RNA repertoires are associated with evolutionarily young transposable elements

Circular RNAs (circRNAs) are found across eukaryotes and can function in post-transcriptional gene regulation. Their biogenesis through a circle-forming backsplicing reaction is facilitated by reverse-complementary repetitive sequences promoting pre-mRNA folding. Orthologous genes from which circRNAs arise, overall contain more strongly conserved splice sites and exons than other genes, yet it remains unclear to what extent this conservation reflects purifying selection acting on the circRNAs themselves. Our analyses of circRNA repertoires from five species representing three mammalian lineages (marsupials, eutherians: rodents, primates) reveal that surprisingly few circRNAs arise from orthologous exonic loci across all species. Even the circRNAs from orthologous loci are associated with young, recently active and species-specific transposable elements, rather than with common, ancient transposon integration events. These observations suggest that many circRNAs emerged convergently during evolution - as a byproduct of splicing in orthologs prone to transposon insertion. Overall, our findings argue against widespread functional circRNA conservation.

A genome-wide circular RNA transcriptome in rat

Circular RNAs (circRNAs) are a novel class of noncoding RNAs that back-splice from 5ʹ donor site and 3ʹ acceptor sites to form a circular structure. A number of circRNAs have been discovered in model organisms including human, mouse, Drosophila, among other organisms. There are a few candidate-based studies on circRNAs in rat, a well-studied model organism as well. A number of pipelines have been published to identify the back splice junctions for the discovery of circRNAs but studies comparing these tools have suggested that a combination of tools would be a better approach to identify high-confidence circRNAs. The availability of a recent dataset of transcriptomes encompassing 11 tissues, 4 developmental stages, and 2 genders motivated us to explore the landscape of circRNAs in the organism in this context. In order to understand the difference among different pipelines, we employed five different combinations of tools to identify circular RNAs from the dataset. We compared the results of the different combination of tools/pipelines with respect to alignment, total number of circRNAs identified and read-coverage. In addition, we identified tissue-specific, development-stage specific and gender-specific circRNAs and further independently validated 16 circRNA junctions out of 24 selected candidates in 5 tissue samples and estimated the quantitative expression of five circRNA candidates using real-time polymerase chain reaction and our analysis suggests three candidates as tissue-enriched. This study is one of the most comprehensive studies which provides a map of circRNAs transcriptome as well as to understand the difference among different computational pipelines in rat.