Comparison of circular RNA prediction tools

EIciRNAs in focus: current understanding and future perspectives

Circular RNAs (circRNAs) are a unique class of covalently closed single-stranded RNA molecules that play diverse roles in normal physiology and pathology. Among the major types of circRNA, exon-intron circRNA (EIciRNA) distinguishes itself by its sequence composition and nuclear localization. Recent RNA-seq technologies and computational methods have facilitated the detection and characterization of EIciRNAs, with features like circRNA intron retention (CIR) and tissue-specificity being characterized. EIciRNAs have been identified to exert their functions via mechanisms such as regulating gene transcription, and the physiological relevance of EIciRNAs has been reported. Within this review, we present a summary of the current understanding of EIciRNAs, delving into their identification and molecular functions. Additionally, we emphasize factors regulating EIciRNA biogenesis and the physiological roles of EIciRNAs based on recent research. We also discuss the future challenges in EIciRNA exploration, underscoring the potential for novel functions and functional mechanisms of EIciRNAs for further investigation.

Dysregulated circular RNAs in rheumatoid arthritis: Cellular roles and clinical prospects

Rheumatoid arthritis (RA) is still a healthcare challenge, although current therapeutic strategies have substantially improved its clinical outcomes. The development of novel biomarkers and treatments can increase the likelihood of identification and disease remission in RA patients, especially for patients with seronegative RA and difficult-to-treat RA (D2T RA). Circular RNAs (circRNAs), a novel non-coding RNA species, have been reported to play crucial roles in various biological process of RA. The mechanistic functions of the dysregulated circRNAs in RA are primarily associated with miRNA sponging and regulating transcription. CircRNAs acting as miRNA sponges are further summarized by cell types, including fibroblast-like synoviocytes (FLSs), lymphocytes, macrophages, chondrocytes, and mesenchymal stem cells (MSCs)-/plasma-secreted exosomes. Besides, a description of dysregulated circRNAs in blood, synovial tissue and cartilage tissue suggests their diagnostic potential for RA. In addition, some directions for future research are provided to open the possibility that dysregulated cell- and tissue- specific circRNAs constituting a fresh reservoir of therapeutic targets, and biomarkers for diagnosis, predicting response to therapy, drug selection or patient stratification for RA.

Unlocking the potential of circular RNA vaccines: a bioinformatics and computational biology perspective

Bioinformatics has significantly advanced RNA-based therapeutics, particularly circular RNAs (circRNAs), which outperform mRNA vaccines, by offering superior stability, sustained expression, and enhanced immunogenicity due to their covalently closed structure. This review highlights how bioinformatics and computational biology optimise circRNA vaccine design, elucidates internal ribosome entry sites (IRES) selection, open reading frame (ORF) optimisation, codon usage, RNA secondary structure prediction, and delivery system development. While circRNA vaccines may not always surpass traditional vaccines in stability, their production efficiency and therapeutic efficacy can be enhanced through computational strategies. The discussion also addresses challenges and future prospects, emphasizing the need for innovative solutions to overcome current limitations and advance circRNA vaccine applications.

Intron-Derived Lariat RNAs Go Stable

During pre-mRNA splicing, introns are removed by the spliceosome, and the flanking exons are ligated to form mature mRNA, which is subsequently translated into protein. Traditionally, intronic RNAs have been regarded as "junk", presumed to be degraded for nucleotide turnover. Notably, after debranching, some linearized lariat RNAs can be further processed into snoRNAs, miRNAs, and other long non-coding RNAs. However, recent studies have shown that many intron-derived lariat RNAs can escape degradation and remain stable across various eukaryotic organisms, indicating they may play significant roles in cellular processes. Moreover, these naturally retained lariat RNAs are frequently observed in circular forms in vivo, suggesting that their linear tails are highly susceptible to degradation. This highlights lariat RNAs as an important source of circular RNAs. Furthermore, many lariat-derived circRNAs have been detected in the cytoplasm, implying active nuclear export and potential roles in cytoplasmic processes. In this review, we provide an overview of the life cycle of intron-derived lariat RNAs, focusing on their biogenesis, degradation, and retention. We also discuss the mechanisms that enable their resistance to degradation and the biological functions of stable lariat RNAs, shedding light on these seemingly "nonsense" yet inevitably produced non-coding intronic RNAs.

Circular RNAs in cancer stem cells: Insights into their roles and mechanisms (Review)

Cancer stem cells (CSCs) represent a small, yet pivotal subpopulation of tumor cells that play significant roles in tumor initiation, progression and therapeutic resistance. Circular RNAs (circRNAs) are a distinct class of RNAs characterized by their closed‑loop structures, lacking 5' to 3'ends. There is growing evidence that circRNAs are integral to the development and regulation of CSCs. Aberrant expression of circRNAs in CSCs can contribute to oncogenic properties and drug resistance. Specifically, oncogenic circRNAs modulate CSC behavior via key signaling pathways, thereby promoting CSC self‑renewal and maintenance, as well as tumor progression. This review summarizes the latest research on the functional roles and regulatory mechanisms of circRNAs in CSC behavior and discusses potential applications and challenges of targeting circRNAs in CSCs. Understanding the intricate interactions between circRNAs and CSCs may lead to novel therapeutic strategies that effectively combat treatment resistance and improve patient outcomes.

Camelina circRNA landscape: Implications for gene regulation and fatty acid metabolism

Circular RNAs (circRNAs) are closed-loop RNAs forming a covalent bond between their 3' and 5' ends, the back splice junction (BSJ), rendering them resistant to exonucleases and thus more stable compared to linear RNAs. Identification of circRNAs and distinction from their cognate linear RNA is only possible by sequencing the BSJ that is unique to the circRNA. CircRNAs are involved in the regulation of their cognate RNAs by increasing transcription rates, RNA stability, and alternative splicing. We have identified circRNAs from C. sativa that are associated with the regulation of germination, light response, and lipid metabolism. We sequenced light-grown and etiolated seedlings after 5 or 7 days post-germination and identified a total of 3447 circRNAs from 2763 genes. Most circRNAs originate from a single homeolog of the three subgenomes from allohexaploid camelina and correlate with higher ratios of alternative splicing of their cognate genes. A network analysis shows the interactions of select miRNA:circRNA:mRNAs for regulation of transcript stabilities where circRNA can act as a competing endogenous RNA. Several key lipid metabolism genes can generate circRNA, and we confirmed the presence of KASII circRNA as a true circRNA. CircRNA in camelina can be a novel target for breeding and engineering efforts.

Identifying insulin-responsive circRNAs in chicken pectoralis

BACKGROUND

Circular RNAs (circRNAs) are stable, covalently closed non-coding RNAs formed by reverse splicing of precursor mRNA. They play critical roles in various biological processes, including insulin secretion and metabolism. However, their function in avian skeletal muscle's response to insulin remains poorly understood. This study aimed to comprehensively identify insulin-responsive circRNAs and explore their temporal and breed-specific regulation in poultry.

RESULTS

Using strand-specific RNA sequencing (ssRNA-Seq) on the pectoralis muscles of both Arbor Acres (AA) broilers and Silky fowls following insulin administration (5 IU/kg.BW, PBS as control). We identified 2,027 muscle circRNAs. Insulin-responsive circRNAs were detected in Silky fowls (29) and broilers (45) at 120 min, and in broilers (20) at 15 min post-injection. These circRNAs are enriched in processes such as exocrine pancreas development, response to exogenous stimuli, and regulation of intracellular signal transduction, likely mediated through a circRNA-miRNA network. Fewer insulin-responsive circRNAs were shared between time points in broilers (1) or between breeds (3) at 120 min. We further characterized a conserved insulin-responsive circRNA (circINSR), formed by exon 2 of the Insulin Receptor (INSR). The circINSR showed a similar spatiotemporal expression pattern to INSR, but exhibited distinct changes post-insulin administration. In broilers, INSR expression was dynamically modulated, while circINSR was downregulated only at 15 min (P < 0.01). Conversely, glucose did not change muscle circINSR but increased INSR at 10 min (P < 0.01). Energy restriction significantly downregulated circINSR (P < 0.01) without affecting INSR levels, and pyruvate had no effect on either circINSR or INSR levels.

CONCLUSION

This study reveals the dynamic and breed-specific roles of circRNAs, particularly circINSR, in avian skeletal muscle's response to insulin. The distinct regulation of circINSR and INSR under various metabolic conditions suggests a complex regulatory mechanism. These findings provide novel insights into the molecular basis of insulin signaling in avian species and highlight the potential of circRNAs as biomarkers for metabolic regulation.

Use of synthetic circular RNA spike-ins (SynCRS) for normalization of circular RNA sequencing data

High-throughput RNA sequencing enables the quantification of transcript abundance and the identification of novel transcripts in biological samples. These include circular RNAs (circRNAs), a family of alternatively spliced RNA molecules that form a continuous loop. However, quantification and comparison of circRNAs between RNA sequencing libraries remain challenging due to confounding errors introduced during exonuclease digestion, library preparation and RNA sequencing itself. Here we describe a set of synthetic circRNA spike-ins-termed 'SynCRS'-that can be added directly to purified RNA samples before exonuclease digestion and library preparation. SynCRS, introduced either individually or in combinations of varying size and abundance, can be integrated into all next-generation sequencing workflows and, critically, facilitate the quantitative calibration of circRNA transcript abundance between samples, tissue types, species and laboratories. Our step-by-step protocol details the generation of SynCRS and guides users on the stoichiometry of SynCRS spike-in to RNA samples, followed by the bioinformatic steps required to facilitate quantitative comparisons of circRNAs between libraries. The laboratory steps to produce the SynCRS require an additional 3 d on top of the high throughput circRNA sequencing and bioinformatics. The protocol is suitable for users with basic experience in molecular biology and bioinformatics.

Recent advances in investigation of circRNA/lncRNA-miRNA-mRNA networks through RNA sequencing data analysis

Non-coding RNAs (ncRNAs) are RNA molecules that are transcribed from DNA but are not translated into proteins. Studies over the past decades have revealed that ncRNAs can be classified into small RNAs, long non-coding RNAs and circular RNAs by genomic size and structure. Accumulated evidences have eludicated the critical roles of these non-coding transcripts in regulating gene expression through transcription and translation, thereby shaping cellular function and disease pathogenesis. Notably, recent studies have investigated the function of ncRNAs as competitive endogenous RNAs (ceRNAs) that sequester miRNAs and modulate mRNAs expression. The ceRNAs network emerges as a pivotal regulatory function, with significant implications in various diseases such as cancer and neurodegenerative disease. Therefore, we highlighted multiple bioinformatics tools and databases that aim to predict ceRNAs interaction. Furthermore, we discussed limitations of using current technologies and potential improvement for ceRNAs network detection. Understanding of the dynamic interplay within ceRNAs may advance the biological comprehension, as well as providing potential targets for therapeutic intervention.

Transcriptome Analysis of Muscle Growth-Related circRNA in the Pacific Abalone Haliotis discus hanna

(1) Background: Animal growth is a complex process, involving the coordination of a wide variety of genes, non-coding RNAs, and pathways. Circular RNAs (circRNAs) belong to a novel class of functional non-coding RNAs (ncRNAs). They have a distinctive ring structure and are involved in various biological processes, including the proliferation, differentiation, and apoptosis of muscle cells. The Pacific abalone Haliotis discus hannai is an economically valuable mollusk species cultivated in China. However, the modulation of muscle growth by circRNAs in this species is poorly understood. (2) Methods: In this study, we analyzed the muscle transcriptomes of 6 H. discus hannai specimens: three small (S_HD) and three large (L_HD) groups via RNA-seq and bioinformatics technology. (3) Results: The results indicated the presence of 11,744 circRNAs in abalone adductor muscle. Furthermore, the L_HD group had 250 significantly differentially expressed circRNAs (106 upregulated and 144 downregulated) relative to the S_HD group. Moreover, the bioinformatics assessment revealed that circRNAs were related to lipid transporter activity, lipid biosynthetic process, fat digestion and absorption, the single-organism metabolic process, the thyroid hormone signaling pathway, and the hippo signaling pathway, which regulates growth. Seventeen key candidate circRNAs were identified, and a core functional circRNA-miRNA-mRNA network associated with abalone muscle growth was described. Gene expression was verified using qRT-PCR, confirming the accuracy of the RNA-seq data. (4) Conclusion: Overall, this investigation furnishes novel evidence for the potential muscle growth modulatory mechanisms in Pacific abalone. These high-quality circRNA data of abalone muscle provide a reference for functional studies on the abalone genome.

Exosomal circular RNAs in tumor microenvironment: An emphasis on signaling pathways and clinical opportunities

Exosomes can regulate the malignant progression of tumors by carrying a variety of genetic information and transmitting it to target cells. Recent studies indicate that exosomal circular RNAs (circRNAs) regulate multiple biological processes in carcinogenesis, such as tumor growth, metastasis, epithelial-mesenchymal transition, drug resistance, autophagy, metabolism, angiogenesis, and immune escape. In the tumor microenvironment (TME), exosomal circRNAs can be transferred among tumor cells, endothelial cells, cancer-associated fibroblasts, immune cells, and microbiota, affecting tumor initiation and progression. Due to the high stability and widespread presence of exosomal circRNAs, they hold promise as biomarkers for tumor diagnosis and prognosis prediction in blood and urine. In addition, designing nanoparticles targeting exosomal circRNAs and utilizing exosomal circRNAs derived from immune cells or stem cells provide new strategies for cancer therapy. In this review, we examined the crucial role of exosomal circRNAs in regulating tumor-related signaling pathways and summarized the transmission of exosomal circRNAs between various types of cells and their impact on the TME. Finally, our review highlights the potential of exosomal circRNAs as diagnostic and prognostic prediction biomarkers, as well as suggesting new strategies for clinical therapy.

Identification of pathological pathways centered on circRNA dysregulation in association with irreversible progression of Alzheimer's disease

BACKGROUND

Circular RNAs (circRNAs) are highly stable regulators, often accumulated in mammalian brains and thought to serve as "memory molecules" that govern the long process of aging. Mounting evidence demonstrated circRNA dysregulation in the brains of Alzheimer's disease (AD) patients. However, whether and how circRNA dysregulation underlies AD progression remains unexplored.

METHODS

We combined Poly(A)-tailing/RNase R digestion experimental approach with CARP, our published computational framework using pseudo-reference alignment for more sensitive and accurate circRNA detection to identify genome-wide circRNA dysregulation and their downstream pathways in the 5xFAD mouse cerebral cortex between 5 and 7 months of age, a critical window marks the transition from reversible to irreversible pathogenic progression. Dysregulated circRNAs and pathways associated with disease progression in 5xFAD cortex were systematically compared with circRNAs affected in postmortem subcortical areas of a large human AD cohort. A top-ranked circRNA conserved and commonly affected in AD patients and 5xFAD mice was depleted in cultured cells to examine AD-relevant molecular and cellular changes.

RESULTS

We discovered genome-wide circRNA alterations specifically in 5xFAD cortex associated with AD progression, many of which are commonly dysregulated in the subcortical areas of AD patients. Among these circRNAs, circGigyf2 is highly conserved and showed the highest net reduction specifically in the 7-month 5xFAD cortex. CircGIGYF2 level in AD patients' cortices negatively correlated with dementia severity. Mechanistically, we found multiple AD-affected splicing factors that are essential for circGigyf2 biogenesis. Functionally, we identified and experimentally validated the conserved roles of circGigyf2 in sponging AD-relevant miRNAs and AD-associated RNA binding proteins (RBPs), including the cleavage and polyadenylation factor 6 (CPSF6). Moreover, circGigyf2 downregulation in AD promoted silencing activities of its sponged miRNAs and enhanced polyadenylation site processing efficiency of CPSF6 targets. Furthermore, circGigyf2 depletion in a mouse neuronal cell line dysregulated circGigyf2-miRNA and circGigyf2-CPSF6 axes and potentiated apoptotic responses upon insults, which strongly support the causative roles of circGigyf2 deficiency in AD neurodegeneration.

CONCLUSIONS

Together, our results unveiled brain circRNAs associated with irreversible disease progression in an AD mouse model that is also affected in AD patients and identified novel molecular mechanisms underlying the dysregulation of conserved circRNA pathways contributing to AD pathogenesis.

Current Understandings and Open Hypotheses on Extracellular Circular RNAs

Circular RNAs (circRNAs) are closed RNA loops present in humans and other organisms. Various circRNAs have an essential role in diseases, including cancer. Cells can release circRNAs into the extracellular space of adjacent biofluids and can be present in extracellular vesicles. Due to their circular nature, extracellular circRNAs (excircRNAs) are more stable than their linear counterparts and are abundant in many biofluids, such as blood plasma and urine. circRNAs' link with disease suggests their extracellular counterparts have high biomarker potential. However, circRNAs and the extracellular space are challenging research domains, as they consist of complex biological systems plagued with nomenclature issues and a wide variety of protocols with different advantages and disadvantages. Here, we summarize what is known about excircRNAs, the current challenges in the field, and what is needed to improve extracellular circRNA research.

A Representation Learning Approach for Predicting circRNA Back-Splicing Event via Sequence-Interaction-Aware Dual Encoder

Circular RNAs (circRNAs) play a crucial role in gene regulation and association with diseases because of their unique closed continuous loop structure, which is more stable and conserved than ordinary linear RNAs. As fundamental work to clarify their functions, a large number of computational approaches for identifying circRNA formation have been proposed. However, these methods fail to fully utilize the important characteristics of back-splicing events, i.e., the positional information of the splice sites and the interaction features of its flanking sequences, for predicting circRNAs. To this end, we hereby propose a novel approach called SIDE for predicting circRNA back-splicing events using only raw RNA sequences. Technically, SIDE employs a dual encoder to capture global and interactive features of the RNA sequence, and then a decoder designed by the contrastive learning to fuse out discriminative features improving the prediction of circRNAs formation. Empirical results on three real-world datasets show the effectiveness of SIDE. Further analysis also reveals that the effectiveness of SIDE.

CircRNome-wide characterisation reveals the promoting role of circAATF in anti-PD-L1 immunotherapy of gallbladder carcinoma

Circular RNAs (circRNAs) have been shown to play important roles in tumour development and tumour immunology. However, genome-wide characterisation of circRNAs and their roles in the immunology and immunotherapy of gallbladder carcinoma (GBC) has been lacking. We present a comprehensive characterisation of the circRNA landscape in GBC, revealing GBC-specific circRNAs. Our analysis found that circRNAs are significantly enriched in cell proliferation and are involved in cancer-related hallmarks. In particular, circAATF was upregulated in GBC, which was positively correlated with AATF mRNA expression, and promoted GBC cell growth. Through integrating computational and experimental approaches, we revealed that circAATF is positively associated with the CD4+ T cell abundance and PD-L1 level, and enhances the clinical benefits of anti-PD-L1 immunotherapy for GBC. We further demonstrate that circAATF elevates the PD-L1 level by activating phosphorylated AKT and acting as a sponge for miR-142-5p. CircAATF is positively associated with CD4+ T cells and PD-L1 levels and shows potential to aid anti-PD-L1 immunotherapy for GBC. Our study provides insights into roles of circAATF in the tumour development and immunology of GBC and accelerates the development of therapeutic strategies for GBC immunotherapy. HIGHLIGHTS: We present a comprehensive characterisation of circRNA landscape in gallbladder carcinoma (GBC). CircAATF is positively associated with CD4+ T cell abundance and PD-L1 expression and is shown to promote PD-L1 treatment in mouse model. CircAATF can elevate PD-L1 level through phosphorylated AKT and linear AATF, which upregulates PD-L1 by acting as a sponge of miR-142-5p.

Whole transcriptome RNA sequencing provides novel insights into the molecular dynamics of ovarian development in mud crab, Scylla paramamosain after mating

Ovarian development in animals is a complicated biological process, requiring the simultaneous coordination among various genes and pathways. To understand the dynamic changes and molecular regulatory mechanisms of ovarian development in mud crab (Scylla paramamosain), both histological observation and whole transcriptome sequencing of ovarian tissues at different mating stages were implemented in this study. The histological results revealed that ovarian development was delayed in unmated females (60 days after courtship behavior but not mating), who exhibited an oocyte diameter of 56.38 ± 15.17 μm. Conversely, mated females exhibited accelerated the ovarian maturation process, with females reaching ovarian stage III (proliferative stage) 23 days after mating and attained an average oocyte diameter of 132.19 ± 15.07 μm. Thus, mating process is essential in promoting the rapid ovarian development in mud crab. Based on the whole transcriptome sequencing analysis, a total of 518 mRNAs, 1502 lncRNAs, 18 circRNAs and 151 miRNAs were identified to be differentially expressed between ovarian tissues at different mating stages. Notably, six differentially expressed genes (DEGs) associated with ovarian development were identified, including ovary development-related protein, red pigment concentrating hormone receptor, G2/mitotic-specific cyclin-B3-like, lutropin-chorio gonadotropic hormone receptor, renin receptor, and SoxB2. More importantly, both DEGs and targets of differentially expressed non-coding RNAs (DEncRNAs) were enriched in renin-angiotensin system, TGF-β signaling, cell adhesion molecules, MAPK signaling pathway, and ECM-receptor interaction, suggesting that these pathways may play significant roles in the ovarian development of mud crabs. Moreover, competition endogenous RNA (ceRNA) networks were constructed while mRNAs were differentially expressed between mating stages were involved in Gene Ontology (GO) biological processes such as developmental process, reproduction, and growth. These findings could provide solid foundations for the future development of female mud crab maturation enhancement strategy, and improve the understanding of the ovarian maturation process in crustaceans.

CircRNAs involved in the red light of effect on follicle selection in pigeons

Red light (RL) can enhance egg production in poultry. CircRNAs play a crucial role by serving as transcriptional regulators. However, their role in influencing follicle development in White King pigeons remains unexplored. In this study, 54 paired White King pigeons were chosen and divided into RL and white light (WL) groups, each with 3 subgroups. The egg production of paired pigeons in each replicate was recorded for 45 d, and the characteristics of follicle development were monitored during the laying interval (LI). The granulosa cell layer from follicles of the second-largest follicle (F2) was collected, and high-throughput sequencing was performed to elucidate the molecular mechanism of follicle development in pigeons. The study confirmed that RL enhances egg production in pigeons. Additionally, under RL, the F2 follicle was selected, while under WL, small follicles were kept on the third day (LI3). A total of 5,510 circRNAs were identified across all samples, revealing differentially expressed circRNAs (DECs) in various comparisons: 627 in RF1 vs. WF1, 900 in RF2 vs. WF2, 606 in RF1 vs. RF2, and 937 in WF1 vs. WF2. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that host genes of DECs were enriched in pathways like steroid hormone biosynthesis, oocyte meiosis, GnRH signaling pathway, and apoptosis pathway. Moreover, circRNA_5497, circRNA_2016, and circRNA_3328 were common DECs across 4 groups, sharing miRNA binding sites with follicle selection-associated genes. In conclusion, our findings suggest that RL promotes egg production by stimulating follicle selection during LI, offering insights into the regulatory mechanisms of circRNAs in follicle selection under RL. This knowledge can help enhance the reproductive performance of pigeons.

Regulatory role of N6-Methyladenosine on skeletal muscle development in Hu sheep

N6-Methyladenosine (m6A) RNA modification plays an essential role in many biological processes. To investigate the regulatory role of m6A on the skeletal muscle development in Hu sheep, this study took newborn Hu sheep (b_B Group) and six-month-old Hu sheep (s_B Group) as the objects. MeRIP-seq and RNA-Seq analysis techniques were used to detect differentially methylated genes (DMGs) and differentially expressed genes (DEGs) in the longissimus dorsi muscle of Hu sheep at different months of age. Then, conjoint analysis was further employed to screen for key genes involved in skeletal muscle development that are modified by m6A and expressed by mRNA. According to the results of the MeRIP-seq analysis, there were 285 m6A differentially methylated peaks (DMPs) in total between b_B Group and s_B Group, with 192 significant upregulated peaks and 93 significant downregulated peaks. GO and KEGG analysis revealed that DMGs are mainly enriched in actin-binding, cellular transport, and metabolic pathways. According to the results of the RNA-seq analysis, there were 4,349 DEGs in total between b_B Group and s_B Group, with 2010 upregulated genes and 2,339 downregulated genes. DEGs are found to be mainly enriched in the regulation of actin cytoskeleton tissue, AMPK and FoxO signaling pathways, etc. The conjoint analysis demonstrated that 283 genes were both modified by m6A and expressed by mRNA. Among them, three genes relevant to muscle growth (RGMB, MAPK8IP3, and RSPO3) were selected as candidates for quantitative validation, and the results were in line with the sequencing results. The results mentioned above all suggest that m6A plays a certain role in the skeletal muscle development in Hu sheep.

Comprehensive analysis of differentially expressed mRNAs, circRNAs, and miRNAs and their ceRNA network in the testis of cattle-yak, yak, and cattle

Cattle-yak is a hybrid offspring resulting from the crossbreeding of yak and cattle, and it exhibits substantial heterosis in production performance. However, male sterility in cattle-yak remains a concern. Reports suggest that noncoding RNAs are involved in the regulation of spermatogenesis. Therefore, in this study, we comprehensively compared testicular transcription profiles among cattle, yak, and cattle-yak. Numerous differentially expressed genes (DEGs), differentially expressed circRNAs (DECs), and differentially expressed miRNAs (DEMs) were identified in the intersection of two comparison groups, namely cattle versus cattle-yak and yak versus cattle-yak, with the number of DEGs, DECs, and DEMs being 4968, 360, and 59, respectively. The DEGs in cattle-yaks, cattle, and yaks were mainly associated with spermatogenesis, male gamete generation, and sexual reproduction. Concurrently, GO and KEGG analyses indicated that DEC host genes and DEM source genes were involved in the regulation of spermatogenesis. The construction of a potential competing endogenous RNA network revealed that some differentially expressed noncoding RNAs may be involved in regulating the expression of genes related to testicular spermatogenesis, including miR-423-5p, miR-449b, miR-34b/c, and miR-15b, as well as previously unreported miR-6123 and miR-1306, along with various miRNA-circRNA interaction pairs. This study serves as a valuable reference for further investigations into the mechanisms underlying male sterility in cattle-yaks.

Identification and characterization of Varicella Zoster Virus circular RNA in lytic infection

This study investigates the role of circular RNAs (circRNAs) in the context of Varicella-Zoster Virus (VZV) lytic infection. We employ two sequencing technologies, short-read sequencing and long-read sequencing, following RNase R treatment on VZV-infected neuroblastoma cells to identify and characterize both cellular and viral circRNAs. Our large scanning analysis identifies and subsequent experiments confirm 200 VZV circRNAs. Moreover, we discover numerous VZV latency-associated transcripts (VLTs)-like circRNAs (circVLTslytic), which contain multiple exons and different isoforms within the same back-splicing breakpoint. To understand the functional significance of these circVLTslytic, we utilize the Bacteria Artificial Chromosome system to disrupt the expression of viral circRNAs in genomic DNA location. We reveal that the sequence flanking circVLTs' 5' splice donor plays a pivotal role as a cis-acting element in the formation of circVLTslytic. The circVLTslytic is dispensable for VZV replication, but the mutation downstream of circVLTslytic exon 5 leads to increased acyclovir sensitivity in VZV infection models. This suggests that circVLTslytic may have a role in modulating the sensitivity to antiviral treatment. The findings shed new insight into the regulation of cellular and viral transcription during VZV lytic infection, emphasizing the intricate interplay between circRNAs and viral processes.

Computational approaches and challenges in the analysis of circRNA data

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

DeepCMI: a graph-based model for accurate prediction of circRNA-miRNA interactions with multiple information

Recently, the role of competing endogenous RNAs in regulating gene expression through the interaction of microRNAs has been closely associated with the expression of circular RNAs (circRNAs) in various biological processes such as reproduction and apoptosis. While the number of confirmed circRNA-miRNA interactions (CMIs) continues to increase, the conventional in vitro approaches for discovery are expensive, labor intensive, and time consuming. Therefore, there is an urgent need for effective prediction of potential CMIs through appropriate data modeling and prediction based on known information. In this study, we proposed a novel model, called DeepCMI, that utilizes multi-source information on circRNA/miRNA to predict potential CMIs. Comprehensive evaluations on the CMI-9905 and CMI-9589 datasets demonstrated that DeepCMI successfully infers potential CMIs. Specifically, DeepCMI achieved AUC values of 90.54% and 94.8% on the CMI-9905 and CMI-9589 datasets, respectively. These results suggest that DeepCMI is an effective model for predicting potential CMIs and has the potential to significantly reduce the need for downstream in vitro studies. To facilitate the use of our trained model and data, we have constructed a computational platform, which is available at http://120.77.11.78/DeepCMI/. The source code and datasets used in this work are available at https://github.com/LiYuechao1998/DeepCMI.

Current advances in circular RNA detection and investigation methods: Are we running in circles?

Circular RNAs (circRNAs), characterized by their closed-loop structure, have emerged as significant transcriptomic regulators, with roles spanning from microRNA sponging to modulation of gene expression and potential peptide coding. The discovery and functional analysis of circRNAs have been propelled by advancements in both experimental and bioinformatics tools, yet the field grapples with challenges related to their detection, isoform diversity, and accurate quantification. This review navigates through the evolution of circRNA research methodologies, from early detection techniques to current state-of-the-art approaches that offer comprehensive insights into circRNA biology. We examine the limitations of existing methods, particularly the difficulty in differentiating circRNA isoforms and distinguishing circRNAs from their linear counterparts. A critical evaluation of various bioinformatics tools and novel experimental strategies is presented, emphasizing the need for integrated approaches to enhance our understanding and interpretation of circRNA functions. Our insights underscore the dynamic and rapidly advancing nature of circRNA research, highlighting the ongoing development of analytical frameworks designed to address the complexity of circRNAs and facilitate the assessment of their clinical utility. As such, this comprehensive overview aims to catalyze further advancements in circRNA study, fostering a deeper understanding of their roles in cellular processes and potential implications in disease. This article is categorized under: RNA Methods > RNA Nanotechnology RNA Methods > RNA Analyses in Cells RNA Methods > RNA Analyses In Vitro and In Silico.

Systematic identification and characterization of exon-intron circRNAs

Exon-intron circRNAs (EIciRNAs) are a circRNA subclass with retained introns. Global features of EIciRNAs remain largely unexplored, mainly owing to the lack of bioinformatic tools. The regulation of intron retention (IR) in EIciRNAs and the associated functionality also require further investigation. We developed a framework, FEICP, which efficiently detected EIciRNAs from high-throughput sequencing (HTS) data. EIciRNAs are distinct from exonic circRNAs (EcircRNAs) in aspects such as with larger length, localization in the nucleus, high tissue specificity, and enrichment mostly in the brain. Deep learning analyses revealed that compared with regular introns, the retained introns of circRNAs (CIRs) are shorter in length, have weaker splice site strength, and have higher GC content. Compared with retained introns in linear RNAs (LIRs), CIRs are more likely to form secondary structures and show greater sequence conservation. CIRs are closer to the 5'-end, whereas LIRs are closer to the 3'-end of transcripts. EIciRNA-generating genes are more actively transcribed and associated with epigenetic marks of gene activation. Computational analyses and genome-wide CRISPR screening revealed that SRSF1 binds to CIRs and inhibits the biogenesis of most EIciRNAs. SRSF1 regulates the biogenesis of EIciLIMK1, which enhances the expression of LIMK1 in cis to boost neuronal differentiation, exemplifying EIciRNA physiological function. Overall, our study has developed the FEICP pipeline to identify EIciRNAs from HTS data, and reveals multiple features of CIRs and EIciRNAs. SRSF1 has been identified to regulate EIciRNA biogenesis. EIciRNAs and the regulation of EIciRNA biogenesis play critical roles in neuronal differentiation.

CircHIPK3 regulates fatty acid metabolism through miR-637/FASN axis to promote esophageal squamous cell carcinoma

The oncogenic role of circRNA in cancers including esophageal cancer (EC) has been well studied. However, whether and how circRNAs are involved in cancer cell metabolic processes remains largely unknown. Here, we reported that circRNA, circHIPK3, is highly expressed in ESCC cell lines and tissues. Knockdown of circHIPK3 significantly restrained cell proliferation, colony formation, migration, and invasion in vitro and inhibited tumor growth in vivo. Mechanistically, circHIPK3 was found to act as a ceRNA by sponging miR-637 to regulate FASN expression and fatty acid metabolism in ESCC cells. Anti-sense oligonucleotide (ASO) targeting circHIPK3 substantially inhibited ESCC both in vitro and in vivo. Therefore, these results uncover a modulatory axis constituting of circHIPK3/miR-637/FASN may be a potential biomarker and therapeutic target for ESCC in the clinic.

DCDA: CircRNA-Disease Association Prediction with Feed-Forward Neural Network and Deep Autoencoder

Circular RNA is a single-stranded RNA with a closed-loop structure. In recent years, academic research has revealed that circular RNAs play critical roles in biological processes and are related to human diseases. The discovery of potential circRNAs as disease biomarkers and drug targets is crucial since it can help diagnose diseases in the early stages and be used to treat people. However, in conventional experimental methods, conducting experiments to detect associations between circular RNAs and diseases is time-consuming and costly. To overcome this problem, various computational methodologies are proposed to extract essential features for both circular RNAs and diseases and predict the associations. Studies showed that computational methods successfully predicted performance and made it possible to detect possible highly related circular RNAs for diseases. This study proposes a deep learning-based circRNA-disease association predictor methodology called DCDA, which uses multiple data sources to create circRNA and disease features and reveal hidden feature codings of a circular RNA-disease pair with a deep autoencoder, then predict the relation score of the pair by a deep neural network. Fivefold cross-validation results on the benchmark dataset showed that our model outperforms state-of-the-art prediction methods in the literature with the AUC score of 0.9794.

Circular RNAs: Regulators of endothelial cell dysfunction in atherosclerosis

Endothelial cell (EC) dysfunction is associated with atherosclerosis. Circular RNAs (circRNAs) are covalently closed loops formed by back-splicing, are highly expressed in a tissue-specific or cell-specific manner, and regulate ECs mainly through miRNAs (mircoRNAs) or protein sponges. This review describes the regulatory mechanisms and physiological functions of circRNAs, as well as the differential expression of circRNAs in aberrant ECs. This review focuses on their roles in inflammation, proliferation, migration, angiogenesis, apoptosis, senescence, and autophagy in ECs from the perspective of signaling pathways, such as nuclear factor κB (NF-κB), nucleotide-binding domain, leucine-rich-repeat family, pyrin-domain-containing 3 (NLRP3)/caspase-1, Janus kinase/signal transducer and activator of transcription (JAK/STAT), and phosphoinositide-3 kinase/protein kinase B (PI3K/Akt). Finally, we address the issues and recent advances in circRNAs as well as circRNA-mediated regulation of ECs to improve our understanding of the molecular mechanisms underlying the progression of atherosclerosis and provide a reference for studies on circRNAs that regulate EC dysfunction and thus affect atherosclerosis.

cirCodAn: A GHMM-based tool for accurate prediction of coding regions in circRNA

Studies focusing on characterizing circRNAs with the potential to translate into peptides are quickly advancing. It is helping to elucidate the roles played by circRNAs in several biological processes, especially in the emergence and development of diseases. While various tools are accessible for predicting coding regions within linear sequences, none have demonstrated accurate open reading frame detection in circular sequences, such as circRNAs. Here, we present cirCodAn, a novel tool designed to predict coding regions in circRNAs. We evaluated the performance of cirCodAn using datasets of circRNAs with strong translation evidence and showed that cirCodAn outperformed the other tools available to perform a similar task. Our findings demonstrate the applicability of cirCodAn to identify coding regions in circRNAs, which reveals the potential of use of cirCodAn in future research focusing on elucidating the biological roles of circRNAs and their encoded proteins. cirCodAn is freely available at https://github.com/denilsonfbar/cirCodAn.

Virus-derived circular RNAs populate hepatitis C virus-infected cells

It is known that pre-mRNAs in eukaryotic cells can be processed to circular RNAs by a backsplicing mechanism. Circular RNAs have great stability and can sequester proteins or small RNAs to exert functions on cellular pathways. Because viruses often exploit host pathways, we explored whether the RNA genome of the cytoplasmic hepatitis C virus is processed to yield virus-derived circRNAs (vcircRNAs). Computational analyses of RNA-seq experiments predicted that the viral RNA genome is fragmented to generate hundreds of vcircRNAs. More than a dozen of them were experimentally verified by rolling-circle amplification. VcircRNAs that contained the viral internal ribosome entry site were found to be translated into proteins that displayed proviral functions. Furthermore, two highly abundant, nontranslated vcircRNAs were shown to enhance viral RNA abundance. These findings argue that novel vcircRNA molecules modulate viral amplification in cells infected by a cytoplasmic RNA virus.

Transcriptome-Wide N6-Methyladenosine Alternations in Pulmonary Arteries of Monocrotaline-Induced Pulmonary Arterial Hypertension in Rats and Novel Therapeutic Targets

N6-methyladenosine (m6A) is a post-transcriptional epigenetic change with transcriptional stability and functionality regulated by specific m6A-modifying enzymes. However, the significance of genes modified by m6A and enzymes specific to m6A regulation in the context of pulmonary arterial hypertension (PAH) remains largely unexplored. MeRIP-seq and RNA-seq were applied to explore variances in m6A and RNA expression within the pulmonary artery tissues of control and monocrotaline-induced PAH rats. Functional enrichments were analyzed using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. To screen candidate m6A-related genes, the STRING and Metascape databases were used to construct a protein-protein interaction network followed by a real-time PCR validation of their expression. The expression level of an m6A regulator was further investigated using immunohistochemical staining, immunofluorescence, and Western blot techniques. Additionally, proliferation assays were conducted on primary rat pulmonary artery smooth muscle cells (PASMCs). We identified forty-two differentially expressed genes that exhibited either hypermethylated or hypomethylated m6A. These genes are predominantly related to the extracellular matrix structure, MAPK, and PI3K/AKT pathways. A candidate gene, centromere protein F (CENPF), was detected with increased expression in the PAH group. Additionally, we first identified an m6A reader, leucine rich pentatricopeptide repeat containing (LRPPRC), which was downregulated in the PAH rat model. The in vitro downregulation of Lrpprc mediated by siRNA resulted in the enhanced proliferation and elevated expression of Cenpf mRNA in primary rat PASMCs. Our study revealed a modified transcriptome-wide m6A landscape and associated regulatory mechanisms in the pulmonary arteries of PAH rats, potentially offering a novel target for therapeutic strategies in the future.

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

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

State-of-the-Art Circular RNA Analytics Using the Circtools Software Suite

Circular RNAs (circRNAs) are types of RNA molecules that have been discovered relatively recently and have been found to be widely expressed in eukaryotic cells. Unlike canonical linear RNA molecules, circRNAs form a covalently closed continuous loop structure without a 5' or 3' end. They are generated by a process called back-splicing, in which a downstream splice donor site is joined to an upstream splice acceptor site. CircRNAs have been found to play important roles in various biological processes, including gene regulation, alternative splicing, and protein translation. They can act as sponges for microRNAs or RNA-binding proteins and can also encode peptides or proteins. Additionally, circRNAs have been implicated in several diseases, including cancer, neurological disorders, and cardiovascular diseases.This protocol provides all necessary steps to detect and analyze circRNAs in silico from RNA sequencing data using the circtools circRNA analytics software suite. The protocol starts from raw sequencing data with circRNA detection via back-splice events and includes statistical testing of circRNAs as well as primer design for follow-up wet lab experiments.

A systematic analysis of circRNAs in subnuclear compartments

CircRNAs are an important class of RNAs with diverse cellular functions in human physiology and disease. A thorough knowledge of circRNAs including their biogenesis and subcellular distribution is important to understand their roles in a wide variety of processes. However, the analysis of circRNAs from total RNA sequencing data remains challenging. Therefore, we developed Calcifer, a versatile workflow for circRNA annotation. Using Calcifer, we analysed APEX-Seq data to compare circRNA occurrence between whole cells, nucleus and subnuclear compartments. We generally find that circRNAs show higher abundance in whole cells compared to nuclear samples, consistent with their accumulation in the cytoplasm. The notable exception is the single-exon circRNA circCANX(9), which is unexpectedly enriched in the nucleus. In addition, we observe that circFIRRE prevails over the linear lncRNA FIRRE in both the cytoplasm and the nucleus. Zooming in on the subnuclear compartments, we show that circRNAs are strongly depleted from nuclear speckles, indicating that excess splicing factors in this compartment counteract back-splicing. Our results thereby provide valuable insights into the subnuclear distribution of circRNAs. Regarding circRNA function, we surprisingly find that the majority of all detected circRNAs possess complete open reading frames with potential for cap-independent translation. Overall, we show that Calcifer is an easy-to-use, versatile and sustainable workflow for the annotation of circRNAs which expands the repertoire of circRNA tools and allows to gain new insights into circRNA distribution and function.

Circular RNAs: Biogenesis, Functions, and Role in Myocardial Hypertrophy

Circular RNAs (circRNAs) are a large class of endogenous single-stranded covalently closed RNA molecules. High-throughput RNA sequencing and bioinformatic algorithms have identified thousands of eukaryotic circRNAs characterized by high stability and tissue-specific expression pattern. Recent studies have shown that circRNAs play an important role in the regulation of physiological processes in the norm and in various diseases, including cardiovascular disorders. The review presents current concepts of circRNA biogenesis, structural features, and biological functions, describes the methods of circRNA analysis, and summarizes the results of studies on the role of circRNAs in the pathogenesis of hypertrophic cardiomyopathy, the most common inherited heart disease.

HNRNPD regulates the biogenesis of circRNAs and the ratio of mRNAs to circRNAs for a set of genes

Exonic circular RNAs (ecircRNAs) in animal cells are generated by backsplicing, and the biogenesis of ecircRNAs is regulated by an array of RNA binding proteins (RBPs). HNRNPD is a heterogeneous nuclear ribonucleoprotein family member with both cytoplasmic and nuclear roles, and whether HNRNPD regulates the biogenesis of circRNAs remains unknown. In this study, we examine the role of HNRNPD in the biogenesis of ecircRNAs. The levels of ecircRNAs are primarily increased upon depletion of HNRNPD. HNRNPD preferentially binds to motifs enriched with A and U nucleotides, and the flanking introns of ecircRNAs tend to have more numbers and higher intensity of HNRNPD binding sites. The levels of mRNAs are generally not significantly altered in HNRNPD knockout cells. For a small set of genes, the circRNA:mRNA ratio is substantially affected, and the mRNA levels of some of these genes demonstrate a significant decrease in HNRNPD knockout cells. CDK1 is identified as a key gene modulated by HNRNPD in the context of circRNA biogenesis. HNRNPD suppresses the biogenesis of circCDK1 and favours the generation of CDK1 mRNA, and the CDK1 protein is a critical regulator of the cell cycle and apoptosis. HNRNPD can participate in cellular physiology, including the cell cycle and apoptosis, and plays roles in clear cell renal cell carcinoma (ccRCC) by modulating circRNA biogenesis and the mRNA levels of key genes, such as CDK1.

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

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

Identification and analysis of circRNAs in the prefrontal cortices of wild boar and domestic pig

Domestication caused significant differences in morphology and behavior between wild and domestic pigs. However, the regulatory role of circRNA in this event is unclear. Here, we analyzed circRNA expression patterns in the prefrontal cortices of wild boar and domestic pigs to determine the potential role of circRNAs in domestication. We identified a total of 11,375 circRNAs and found that 349 and 354 circRNAs were up-regulated in wild boar and Rongchang pig, respectively. Functional enrichment analysis showed that host genes of significantly highly-expressed circRNAs in wild boar were significantly enriched in neural synapse-related categories and the categories of 'regulation of defense response (p = 0.028)' and 'neural retina development (p = 4.32 × 10-3)'. Host genes of significantly highly-expressed circRNAs in Rongchang pig were specifically involved in 'chordate embryonic development (p = 2.38 × 10-4)'. Additionally, we constructed circRNA-miRNA-mRNA regulatory axes in wild boar and Rongchang pig and found more regulatory axes in wild boar that potentially regulate synaptic activities. We identified multiple circRNAs that may be related to domesticated characteristics, such as ssc_circ_6179 (ssc_circ_6179-ssc-miR-9847-HRH3, related to aggression) and ssc_circ_3027 (ssc_circ_3027-ssc-miR-4334-5p-HCRTR1, related to attention). This study provides a resource for further investigation of the molecular basis of pig domestication.

Comprehensive Analysis of CircRNA Expression Profiles in Multiple Tissues of Pigs

Circular RNAs (circRNAs) are a class of non-coding RNAs with diverse functions, and previous studies have reported that circRNAs are involved in the growth and development of pigs. However, studies about porcine circRNAs over the past few years have focused on a limited number of tissues. Based on 215 publicly available RNA sequencing (RNA-seq) samples, we conducted a comprehensive analysis of circRNAs in nine pig tissues, namely, the gallbladder, heart, liver, longissimus dorsi, lung, ovary, pituitary, skeletal muscle, and spleen. Here, we identified a total of 82,528 circRNAs and discovered 3818 novel circRNAs that were not reported in the CircAtlas database. Moreover, we obtained 492 housekeeping circRNAs and 3489 tissue-specific circRNAs. The housekeeping circRNAs were enriched in signaling pathways regulating basic biological tissue activities, such as chromatin remodeling, nuclear-transcribed mRNA catabolic process, and protein methylation. The tissue-specific circRNAs were enriched in signaling pathways related to tissue-specific functions, such as muscle system process in skeletal muscle, cilium organization in pituitary, and cortical cytoskeleton in ovary. Through weighted gene co-expression network analysis, we identified 14 modules comprising 1377 hub circRNAs. Additionally, we explored circRNA-miRNA-mRNA networks to elucidate the interaction relationships between tissue-specific circRNAs and tissue-specific genes. Furthermore, our conservation analysis revealed that 19.29% of circRNAs in pigs shared homologous positions with their counterparts in humans. In summary, this extensive profiling of housekeeping, tissue-specific, and co-expressed circRNAs provides valuable insights into understanding the molecular mechanisms of pig transcriptional expression, ultimately deepening our understanding of genetic and biological processes.

An integrative pipeline for circular RNA quantitative trait locus discovery with application in human T cells

MOTIVATION

Molecular quantitative trait locus (QTL) mapping has proven to be a powerful approach for prioritizing genetic regulatory variants and causal genes identified by genome-wide association studies. Recently, this success has been extended to circular RNA (circRNA), a potential group of RNAs that can serve as markers for the diagnosis, prognosis, or therapeutic targets of various human diseases. However, a well-developed computational pipeline for circRNA QTL (circQTL) discovery is still lacking.

RESULTS

We introduce an integrative method for circQTL mapping and implement it as an automated pipeline based on Nextflow, named cscQTL. The proposed method has two main advantages. Firstly, cscQTL improves the specificity by systematically combining outputs of multiple circRNA calling algorithms to obtain highly confident circRNA annotations. Secondly, cscQTL improves the sensitivity by accurately quantifying circRNA expression with the help of pseudo references. Compared to the single method approach, cscQTL effectively identifies circQTLs with an increase of 20%-100% circQTLs detected and recovered all circQTLs that are highly supported by the single method approach. We apply cscQTL to a dataset of human T cells and discover genetic variants that control the expression of 55 circRNAs. By colocalization tests, we further identify circBACH2 and circYY1AP1 as potential candidates for immune disease regulation.

AVAILABILITY AND IMPLEMENTATION

cscQTL is freely available at: https://github.com/datngu/cscQTL and https://doi.org/10.5281/zenodo.7851982.

Dynamic characterization and interpretation for protein-RNA interactions across diverse cellular conditions using HDRNet

RNA-binding proteins play crucial roles in the regulation of gene expression, and understanding the interactions between RNAs and RBPs in distinct cellular conditions forms the basis for comprehending the underlying RNA function. However, current computational methods pose challenges to the cross-prediction of RNA-protein binding events across diverse cell lines and tissue contexts. Here, we develop HDRNet, an end-to-end deep learning-based framework to precisely predict dynamic RBP binding events under diverse cellular conditions. Our results demonstrate that HDRNet can accurately and efficiently identify binding sites, particularly for dynamic prediction, outperforming other state-of-the-art models on 261 linear RNA datasets from both eCLIP and CLIP-seq, supplemented with additional tissue data. Moreover, we conduct motif and interpretation analyses to provide fresh insights into the pathological mechanisms underlying RNA-RBP interactions from various perspectives. Our functional genomic analysis further explores the gene-human disease associations, uncovering previously uncharacterized observations for a broad range of genetic disorders.

Circular RNAs in the human brain are tailored to neuron identity and neuropsychiatric disease

Little is known about circular RNAs (circRNAs) in specific brain cells and human neuropsychiatric disease. Here, we systematically identify over 11,039 circRNAs expressed in vulnerable dopamine and pyramidal neurons laser-captured from 190 human brains and non-neuronal cells using ultra-deep, total RNA sequencing. 1526 and 3308 circRNAs are custom-tailored to the cell identity of dopamine and pyramidal neurons and enriched in synapse pathways. 29% of Parkinson's and 12% of Alzheimer's disease-associated genes produced validated circRNAs. circDNAJC6, which is transcribed from a juvenile-onset Parkinson's gene, is already dysregulated during prodromal, onset stages of common Parkinson's disease neuropathology. Globally, addiction-associated genes preferentially produce circRNAs in dopamine neurons, autism-associated genes in pyramidal neurons, and cancers in non-neuronal cells. This study shows that circular RNAs in the human brain are tailored to neuron identity and implicate circRNA-regulated synaptic specialization in neuropsychiatric diseases.

Going circular: history, present, and future of circRNAs in cancer

To date, thousands of highly abundant and conserved single-stranded RNA molecules shaped into ring structures (circRNAs) have been identified. CircRNAs are multifunctional molecules that have been shown to regulate gene expression transcriptionally and post-transcriptionally and exhibit distinct tissue- and development-specific expression patterns associated with a variety of normal and disease conditions, including cancer pathogenesis. Over the past years, due to their intrinsic stability and resistance to ribonucleases, particular attention has been drawn to their use as reliable diagnostic and prognostic biomarkers in cancer diagnosis, treatment, and prevention. However, there are some critical caveats to their utility in the clinic. Their circular shape limits their annotation and a complete functional elucidation is lacking. This makes their detection and biomedical application still challenging. Herein, we review the current knowledge of circRNA biogenesis and function, and of their involvement in tumorigenesis and potential utility in cancer-targeted therapy.

Large-scale benchmarking of circRNA detection tools reveals large differences in sensitivity but not in precision

The detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed using computational tools. Numerous such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools detected more than 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were validated using three orthogonal methods. Generally, tool-specific precision is high and similar (median of 98.8%, 96.3% and 95.5% for qPCR, RNase R and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant differentiators. Of note, precision values are lower when evaluating low-abundance circRNAs. We also show that the tools can be used complementarily to increase detection sensitivity. Finally, we offer recommendations for future circRNA detection and validation.

A review on circular RNAs and bacterial infections

Bacterial infections and related diseases have been a major burden on social public health and economic stability around the world. However, the effective diagnostic methods and therapeutic approaches to treat bacterial infections are still limited. As a group of non-coding RNA, circular RNAs (circRNAs) that were expressed specifically in host cells and played a key regulatory role have the potential to be of diagnostic and therapeutic value. In this review, we systematically summarize the role of circRNAs in common bacterial infections and their potential roles as diagnostic markers and therapeutic targets.

CircPrime: a web-based platform for design of specific circular RNA primers

BACKGROUND

Circular RNAs (circRNAs) are covalently closed-loop RNAs with critical regulatory roles in cells. Tens of thousands of circRNAs have been unveiled due to the recent advances in high throughput RNA sequencing technologies and bioinformatic tools development. At the same time, polymerase chain reaction (PCR) cross-validation for circRNAs predicted by bioinformatic tools remains an essential part of any circRNA study before publication.

RESULTS

Here, we present the CircPrime web-based platform, providing a user-friendly solution for DNA primer design and thermocycling conditions for circRNA identification with routine PCR methods.

CONCLUSIONS

User-friendly CircPrime web platform ( http://circprime.elgene.net/ ) works with outputs of the most popular bioinformatic predictors of circRNAs to design specific circular RNA primers. CircPrime works with circRNA coordinates and any reference genome from the National Center for Biotechnology Information database).

Computational approaches for circRNAs prediction and in silico characterization

Circular RNAs (circRNAs) are single-stranded and covalently closed non-coding RNA molecules originated from RNA splicing. Their functions include regulatory potential over other RNA species, such as microRNAs, messenger RNAs and RNA binding proteins. For circRNA identification, several algorithms are available and can be classified in two major types: pseudo-reference-based and split-alignment-based approaches. In general, the data generated from circRNA transcriptome initiatives is deposited on public specific databases, which provide a large amount of information on different species and functional annotations. In this review, we describe the main computational resources for the identification and characterization of circRNAs, covering the algorithms and predictive tools to evaluate its potential role in a particular transcriptomics project, including the public repositories containing relevant data and information for circRNAs, recapitulating their characteristics, reliability and amount of data reported.

Assessing the impacts of various factors on circular RNA reliability

Circular RNAs (circRNAs) are non-polyadenylated RNAs with a continuous loop structure characterized by a non-colinear back-splice junction (BSJ). Although millions of circRNA candidates have been identified, it remains a major challenge for determining circRNA reliability because of various types of false positives. Here, we systematically assess the impacts of numerous factors related to circRNA identification, conservation, biogenesis, and function on circRNA reliability by comparisons of circRNA expression from mock and the corresponding colinear/polyadenylated RNA-depleted datasets based on three different RNA treatment approaches. Eight important indicators of circRNA reliability are determined. The relative contribution to variability explained analyses reveal that the relative importance of these factors in affecting circRNA reliability in descending order is the conservation level of circRNA, full-length circular sequences, supporting BSJ read count, both BSJ donor and acceptor splice sites at the same colinear transcript isoforms, both BSJ donor and acceptor splice sites at the annotated exon boundaries, BSJs detected by multiple tools, supporting functional features, and both BSJ donor and acceptor splice sites undergoing alternative splicing. This study thus provides a useful guideline and an important resource for selecting high-confidence circRNAs for further investigations.

Genome-wide identification and characterization of circRNAs in wheat tiller

KEY MESSAGE

This study found that the intergenic circRNAs of wheat were more abundant than those of other plants. More importantly, a circRNA-mediated network associated with tillering was constructed for the first time. Circular RNAs (circRNAs) are a class of endogenous non-coding RNAs with covalently closed circular structures, which play an important role in transcriptional and post-transcriptional regulation. Tiller is an important agronomic trait that determines plant morphological architecture and affects spike number in wheat. However, no studies on the characteristics and functions of circRNAs involved in the regulation of wheat tiller. Here, we performed a genome-wide identification of circRNAs using ribosomal-depleted RNA-seq from wheat tiller of two pairs near-isogenic lines. A total of 686 circRNAs were identified and distributed on 21 chromosomes of wheat, of which 537 were novel circRNAs. Unlike other plants, the majority of these circRNAs (61.8%) were derived from intergenic regions. One circRNA-mediated network associated with tillering was constructed through weighted gene co-expression network analysis, including 323 circRNAs, 117 miRNAs, and 968 mRNAs. GO and pathway enrichment analysis of mRNAs suggested that these circRNAs are involved in cell cycle, ncRNA export from nucleus, developmental process, plant hormone signal transduction, MAPK signaling pathway, RNA degradation. Of these circRNAs, ten circRNAs are associated with known tillering/branching genes in rice or Arabidopsis thaliana, including OsCesA7, EBR1, DTE1, CRD1, LPA1, PAY1, LRK1, OsNR2, OsCCA1, OsBZR1. In summary, we present the first study of the identification and characterization of circRNAs in wheat tiller, and the results suggest these circRNAs associated with tillering could play an important role in wheat tiller formation and development.

Circular RNAs in the human brain are tailored to neuron identity and neuropsychiatric disease

Little is known about circular RNAs (circRNAs) in specific brain cells and human neuropsychiatric disease. Here, we systematically identified over 11,039 circRNAs expressed in vulnerable dopamine and pyramidal neurons laser-captured from 190 human brains and non-neuronal cells using ultra-deep, total RNA sequencing. 1,526 and 3,308 circRNAs were custom-tailored to the cell identity of dopamine and pyramidal neurons and enriched in synapse pathways. 88% of Parkinson's and 80% of Alzheimer's disease-associated genes produced circRNAs. circDNAJC6, produced from a juvenile-onset Parkinson's gene, was already dysregulated during prodromal, onset stages of common Parkinson's disease neuropathology. Globally, addiction-associated genes preferentially produced circRNAs in dopamine neurons, autism-associated genes in pyramidal neurons, and cancers in non-neuronal cells. This study shows that circular RNAs in the human brain are tailored to neuron identity and implicate circRNA- regulated synaptic specialization in neuropsychiatric diseases.