Back-splicing Sites Research Articles

Unraveling the modulatory manner and function of circRNAs in the Asian honey bee larval guts.

Circular RNAs (circRNAs) are a class of non-coding RNAs (ncRNAs) that can participate in biological processes such as gene expression, growth, and development. However, little has been explored about the function of circRNAs in the development of Apis cerana larval guts. By using our previously gained deep sequencing data from the guts of A. cerana worker larvae at 4-, 5-, and 6-day-old (Ac4, Ac5, and Ac6 groups), the expression pattern and regulatory role of circular RNAs (circRNAs) during the development process was comprehensively investigated, with a focus on differentially expressed circRNAs (DEcircRNAs) relevant to immunity pathways and developmental signaling pathways, followed by validation of the binding relationships among a key competing endogenous RNA (ceRNA) axis. Here, 224 (158) DEcircRNAs were detected in the Ac4 vs. Ac5 (Ac5 vs. Ac6) comparison group. It's suggested that 172 (123) parental genes of DEcircRNAs were involved in 26 (20) GO terms such as developmental process and metabolic process and 138 (136) KEGG pathways like Hippo and Wnt signaling pathways. Additionally, ceRNA network analysis indicated that 21 (11) DEcircRNAs could target seven (three) DEmiRNAs, further targeting 324 (198) DEmRNAs. These DEmRNAs can be annotated to 33 (26) GO terms and 168 (200) KEGG pathways, including 12 (16) cellular and humoral immune pathways (endocytosis, lysosome, Jak-STAT, etc.) and 10 (nine) developmental signaling pathways (Hippo, mTOR, Hedgehog, etc.). Interestingly, DEcircRNAs in these two comparison groups could target the same ace-miR-6001-y, forming complex sub-networks. The results of PCR and Sanger sequencing confirmed the back-splicing sites within four randomly selected DEcircRNAs. RT-qPCR detection of these four DEcircRNAs verified the reliability of the used transcriptome data. The results of dual-luciferase reporter assay verified the binding relationships between novel_circ_001627 and ace-miR-6001-y and between ace-miR-6001-y and apterous-like. Our data demonstrated that DEcircRNAs were likely to modulate the developmental process of the A. cerana worker larval guts via regulation of parental gene transcription and ceRNA network, and novel_circ_001627/ace-miR-6001-y/apterous-like was a potential regulatory axis in the larval gut development. Findings from this work offer a basis and a candidate ceRNA axis for illustrating the circRNA-modulated mechanisms underlying the A. cerana larval guts.

Dynamics and regulatory role of circRNAs in Asian honey bee larvae following fungal infection

Non-coding RNA (ncRNA) plays a vital part in the regulation of immune responses, growth, and development in plants and animals. Here, the identification, characteristic analysis, and molecular verification of circRNAs in Apis cerana cerana worker larval guts were conducted, followed by in-depth investigation of the expression pattern of larval circRNAs during Ascosphaera apis infection and exploration of the potential regulatory part of differentially expressed circRNAs (DEcircRNAs) in host immune responses. A total of 3178 circRNAs in the larval guts of A. c. cerana were identified, with a length distribution ranging from 15 to 96,007 nt. Additionally, 155, 95, and 86 DEcircRNAs were identified in the in the 4-, 5-, and 6-day-old larval guts following A. apis infection. These DEcircRNAs were predicted to target 29, 25, and 18 parental genes relevant to 12, 20, and 17 GO terms as well as 144, 114, and 61 KEGG pathways, including 5 cellular and 4 humoral immune pathways. Complex competing endogenous RNA (ceRNA) regulatory networks were detected as being formed among DEcircRNAs, DEmiRNAs, and DEmRNAs. The target DEmRNAs were engaged in 36, 47, and 47 GO terms as well as 331, 332, and 331 pathways, including 6 cellular and 6 humoral immune pathways. Further, 19 DEcircRNAs, 5 DEmiRNAs, and 3 mRNAs were included in the sub-networks relative to 3 antioxidant enzymes. Finally, back-splicing sites within 15 circRNAs and the difference in the 15 DEcircRNAs’ expression between uninoculated and A. apis–inoculated larval guts were confirmed based on molecular methods. These findings not only enrich our understanding of bee host–fungal pathogen interactions but also lay a foundation for illuminating the mechanism underlying the DEcircRNA-mediated immune defense of A. c. cerana larvae against A. apis invasion.Key points• The expression pattern of circRNAs was altered in the A. cerana worker larval guts following A. apis infection.• Back-splicing sites within 15 A. cerana circRNAs were verified using molecular approaches. DEcircRNAs potentially modulated immune responses and antioxidant enzymes in A. apis–challenged host guts.

Identification and characterization of circular RNAs involved in the fertility stability of cotton CMS-D2 restorer line under heat stress

BackgroundAs a vital type of noncoding RNAs, circular RNAs (circRNAs) play important roles in plant growth and development and stress response. However, little is known about the biological roles of circRNAs in regulating the stability of male fertility restoration for cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) cotton under high-temperature (HT) stress.ResultsIn this study, RNA-sequencing and bioinformatics analysis were performed on pollen grains of isonuclear alloplasmic near-isogenic restorer lines NH [N(Rf1rf1)] and SH [S(Rf1rf1)] with obvious differences in fertility stability under HT stress at two environments. A total of 967 circRNAs were identified, with 250 differentially expressed under HT stress. We confirmed the back-splicing sites of eight selected circRNAs using divergent primers and Sanger sequencing. Tissue-specific expression patterns of five differentially expressed circRNAs (DECs) were also verified by RT-PCR and qRT-PCR. Functional enrichment and metabolic pathway analysis revealed that the parental genes of DECs were significantly enriched in fertility-related biological processes such as pollen tube guidance and cell wall organization, as well as the Pentose and glucuronate interconversions, Steroid biosynthesis, and N-Glycan biosynthesis pathways. Moreover, we also constructed a putative circRNA-mediated competing endogenous RNA (ceRNA) network consisting of 21 DECs, eight predicted circRNA-binding miRNAs, and their corresponding 22 mRNA targets, especially the two ceRNA modules circRNA346-miR159a-MYB33 and circRNA484-miR319e-MYB33, which might play important biological roles in regulating pollen fertility stability of cotton CMS-D2 restorer line under HT stress.ConclusionsThrough systematic analysis of the abundance, characteristics and expression patterns of circRNAs, as well as the potential functions of their parent genes, our findings suggested that circRNAs and their mediated ceRNA networks acted vital biological roles in cotton pollen development, and might be also essential regulators for fertility stability of CMS-D2 restorer line under heat stress. This study will open a new door for further unlocking complex regulatory mechanisms underpinning the fertility restoration stability for CMS-D2 in cotton.

Base-Editor-Mediated circRNA Knockout by Targeting Predominantly Back-Splice Sites.

Back-splicing of eukaryotic exon(s) leads to the production of covalently closed circular RNAs (circRNAs). Generally, most circRNAs contain overlapping sequences to their cognate linear RNAs from the same gene loci, leading to difficulties in distinguishing them from each other. A recent study has shown that some circRNAs can be specifically depleted by using base editing systems to target their predominantly back-splice sites for circularization, suggesting an efficient approach for circRNA knockout (KO). Here, we describe the detailed protocol for applying base editors to disrupt back-splice sites of predominantly circularized exons for circRNA KO at the genomic DNA level.

The Functional Circular RNA Screening via RfxCas13d/BSJ-gRNA System.

Although discovered decades ago, functions of circular RNAs (circRNAs) produced from exon(s) back-splicing of pre-mRNAs have only been unveiled recently. As circRNAs share overlapping sequences with their cognate linear RNAs, except for the back-splicing junction sites, it is difficult to distinguish circRNAs from cognate mRNAs in functional studies. In this chapter, we describe a programmable method for the large-scale functional circRNA screening based on the RNA-guided, RNA-targeting CRISPR-Cas13 (RfxCas13d) system. This method can be applied both in vivo and in cell to explore highly expressed circRNAs that may influence cell growth, either under natural conditions or in response to environmental stimulation, without disturbing cognate linear mRNAs.

Novel Insights into the circRNA-Modulated Developmental Mechanism of Western Honey Bee Larval Guts.

Circular RNAs (circRNAs) are a class of novel non-coding RNAs (ncRNAs) that play essential roles in the development and growth of vertebrates through multiple manners. However, the mechanism by which circRNAs modulate the honey bee gut development is currently poorly understood. Utilizing the transcriptome data we obtained earlier, the highly expressed circRNAs in the Apis mellifera worker 4-, 5-, and 6-day-old larval guts were analyzed, which was followed by an in-depth investigation of the expression pattern of circRNAs during the process of larval guts development and the potential regulatory roles of differentially expressed circRNAs (DEcircRNAs). In total, 1728 expressed circRNAs were detected in the A. mellifera larval guts. Among the most highly expressed 10 circRNAs, seven (novel_circ_000069, novel_circ_000027, novel_circ_000438, etc.) were shared by the 4-, 5-, and 6-day-old larval guts. In addition, 21 (46) up-regulated and 22 (27) down-regulated circRNAs were, respectively, screened in the Am4 vs. Am5 (Am5 vs. Am6) comparison groups. Additionally, nine DEcircRNAs, such as novel_circ_000340, novel_circ_000758 and novel_circ_001116, were shared by these two comparison groups. These DEcircRNAs were predicted to be transcribed from 14 and 29 parental genes; these were respectively annotated to 15 and 22 GO terms such as biological regulation and catalytic activity as well as 16 and 21 KEGG pathways such as dorsoventral axis formation and apoptosis. Moreover, a complicated competing endogenous RNA (ceRNA) network was observed; novel_circ_000838 in the Am4 vs. Am5 comparison group potentially targeted ame-miR-6000a-3p, further targeting 518 mRNAs engaged in several developmental signaling pathways (e.g., TGF-beta, hedgehog, and wnt signaling pathway) and immune pathways (e.g., phagosome, lysosome, and MAPK signaling pathway). The results demonstrated that the novel_circ_000838-ame-miR-6000a-3p axis may plays a critical regulatory part in the larval gut development and immunity. Furthermore, back-splicing sites of six randomly selected DEcircRNAs were amplified and verified by PCR; an RT-qPCR assay of these six DEcircRNAs confirmed the reliability of the used high-throughput sequencing data. Our findings provide a novel insight into the honey bee gut development and pave a way for illustration of the circRNA-modulated developmental mechanisms underlying the A. mellifera worker larval guts.

Regulation of circular RNAs in patients with cardiogenic shock

Abstract Background The mortality rate of cardiogenic shock (CS) is around 50%. There is a need for biomarkers to predict the prognosis of CS. Due to their shape, circular RNAs (circRNAs) exhibit higher stability than other classes of RNA. Regulation of circRNAs has been observed in cardiovascular diseases such as myocardial infarction, but their regulation in CS is largely unknown. Monocytes can be classified into different populations based on their surface markers. Non-classical and intermediate monocytes are characterized by CD16 expression and have been associated with non-inflammatory function such as vascular healing. Methods and Results Blood samples were taken from patients with CS after myocardial infarction (MI) (n=4) on admission and 24 hours after revascularization and from healthy controls. Bioinformatic analysis by RNA-Sequencing revealed the expression of a total of 1819 circRNAs in the blood. 457 circRNAs showed altered expression in CS at the time of admission compared to controls, with 35 circRNAs upregulated and 422 circRNAs downregulated. Of these 457 circRNAs, 80 circRNAs were no longer regulated 24 hours after revascularization. This indicates dynamic regulation of circRNA during CS. Using divergently oriented primers spanning the backsplice junction of selected circRNA candidates, the presence of 17 of the most highly regulated or expressed circRNAs was confirmed by RT-PCR in peripheral blood mononuclear cells. We confirmed the circularity of the 17 circRNA candidates by RNAse-R digestion and oligo-dT selection. For one of the top candidates, circRNA-4, we confirmed upregulation in a larger group of patients with CS after MI (n=52) compared to healthy controls (n=38, 10.9-fold, p < 0.0001). Interestingly, patients with MI without CS (n=37) showed no increase in circRNA-4 compared to healthy controls (n=31, p=0.2), suggesting a potential CS-specific regulation of circRNA-4 independent of myocardial infarction. Single cell sequencing of CS patients (N=4) revealed an increase in CD16+ monocytes and enriched expression of the circRNA-4 host gene. Interestingly, knockdown of circRNA-4 in human monocytes using siRNA targeting the circRNA-4-specific backsplice site resulted in a reduction of the non-classical monocyte marker CD16, suggesting a role for circRNA-4 in monocyte subpopulation switching. Conclusion In conclusion, the data reveal a dynamic and significant regulation of circRNAs in the blood of patients with CS after MI. circRNA-4 showed a CS-specific upregulation that was not primarily dependent on MI. After knockdown of circRNA-4, a reduced expression of the surface marker CD16 was also shown in monocytes indicating the involvement in monocyte subpopulation switch toward non-classical monocytes. Thus, these data set the stage to further study circRNAs as potential new specific markers in acute diseases such as cardiogenic shock.

CircMSH3 is a potential biomarker for the diagnosis of colorectal cancer and affects the distant metastasis of colorectal cancer

Objectives To identify the most significantly differentially expressed circular RNAs (circRNAs) in colorectal cancer (CRC) tissues in terms of their expression levels and circularity, and to analyze the relationship between their expression levels and the clinical characteristics of patients. Methods circRNA RNA-seq technology was used to screen differentially expressed circRNAs in CRC. Sanger sequencing was used to identify circRNA back-splice junction sites. The relative expression levels of hsa_circ_0003761 (circMSH3) in CRC tissues and cell lines were detected by quantitative real-time fluorescence PCR technology. An RNA-protein pull-down assay was used to detect protein binding to circRNAs. Dual-luciferase reporter gene vectors were constructed to verify that circRNAs bind to microRNAs. Results Four hundred twenty circRNAs were found to be upregulated, and 616 circRNAs were downregulated. circMSH3 was derived from the MutS homolog 3 (MSH3) gene and was formed by a loop of exons 9, 10, 11, and 12. In 110 pairs of CRC and adjacent tissues, circMSH3 expression was 4.487-fold higher in CRC tissues. circMSH3 was also highly expressed in the HT-29 and LOVO CRC cell lines. The expression level of circMSH3 was associated with distant metastasis in CRC patients (P = 0.043); the area under the curve (AUC) of circMSH3 for CRC diagnosis was 0.75, with a sensitivity and specificity of 70.9% and 66.4%, respectively. circMSH3 could bind to a variety of proteins, mainly those involved in RNA transcription, splicing, cell cycle, and cell junctions. Furthermore, circMSH3 could bind to miR-1276, miR-942-5p, and miR-409-3p. Conclusion circMSH3 is a potential biomarker for the diagnosis of CRC and affects the distant metastasis of CRC. Multiple RNA-binding protein binds to circMSH3, and circMSH3 binds to miR-1276, miR-942-5p, and miR-409-3p, thereby affecting the expression of circMSH3.

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

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

Circular RNA ame_circ_000115 regulates expression of genes in larval gusts of Apis mellifera ligustica stressed by Ascosphaera apis

Circular RNAs (circRNAs) are a new class of non-coding RNAs, which have been confirmed to regulate insect gene expression and immune response through multiple manners such as competing endogenous RNA (ceRNA) regulatory network. Currently, function of circRNA in honey bee immune response remains unclear. In this study, PCR and Sanger sequencing were performed to validate the back splicing (BS) site of ame_circ_000115 (in short ac115). RT-qPCR was used to detect the expression profile of ac115 in larval guts of Apis mellifera ligustica stressed by Ascosphaera apis. Dual-luciferase reporter gene assay was conducted to verify the binding relationship between ac115 and ame-miR-13b. Interference of ac115 in larval guts was carried out by feeding specific siRNA, followed by determination of the effect of ac115 interference on expression of six genes relevant to host immune response. The results confirmed the existence of BS site within ac115. Compared with the un-inoculated group, the expression of ac115 in 4-day-old larval gut of the A. apis-inoculated group was up-regulated with extreme significance (P < 0.000 1), while that in 5- and 6-day-old larval guts were significantly up-regulated (P < 0.05). The brightness of specific band for ac115 in 4-, 5- and 6-day-old larval guts of the siRNA-circ_000115-fed group gradually became weak, whereas that of the siRNA-scrambl-fed group was pretty high without obvious variation. Compared with that of the siRNA-scramble-fed group, the expression of ac115 in 4-day-old larval gut of the siRNA-circ_000115-fed group was significantly down-regulated (P < 0.05), whereas that of the 5- and 6-day-old larval guts were down-regulated with extreme significance (P < 0.001). Ame-miR-13b was truly existed and expressed in A. m. ligustica larval guts, and there was true binding relationship between ac115 and ame-miR-13b. Compared with that of the siRNA-scramble-fed group, the expression of antimicrobial peptide genes hymenoptaecin and abaecin in 6-day-old larval gut of the siRNA-circ_000115-fed group was significantly up-regulated (P < 0.05), while that of ecdysone receptor (Ecr) was down-regulated with extreme significance (P < 0.01). These results indicate that ac115 is truly expressed in A. m. ligustica larval guts, BS site truly exists within ac115, and effective interference of ac115 in A. m. ligustica larval guts can be achieved via feeding siRNA. Moreover, ac115 potentially regulates Ecr expression through adsorption of ame-miR-13b and expression of hymenoptaecin and abaecin using a non-ceRNA manner, further participating in host stress-response.

Unveiling the circRNA-Mediated Immune Responses of Western Honey Bee Larvae to Ascosphaera apis Invasion.

Western honey bee (Apis mellifera), a eusocial insect with a superior economic and ecological value, is widely used in the beekeeping industry throughout the world. As a new class of non-coding RNAs (ncRNAs), circular RNAs (circRNAs) participate in the modulation of considerable biological processes, such as the immune response via diverse manners. Here, the identification, characteristic investigation, and molecular verification of circRNAs in the Apis mellifera ligustica larval guts were conducted, and the expression pattern of larval circRNAs during the Ascosphaera apis infection was analyzed, followed by the exploration of the potential regulatory part of differentially expressed circRNAs (DEcircRNAs) in host immune responses. A total of 2083 circRNAs in the larval guts of A. m. ligustcia were identified, with a length distribution ranging from 106 nt to 92,798 nt. Among these, exonic circRNAs were the most abundant type and LG1 was the most distributed chromosome. Additionally, 25, 14, and 30 up-regulated circRNAs as well as 26, 25, and 62 down-regulated ones were identified in the A. apis-inoculated 4-, 5-, and 6-day-old larval guts in comparison with the corresponding un-inoculated larval guts. These DEcircRNAs were predicted to target 35, 70, and 129 source genes, which were relative to 12, 23, and 20 GO terms as well as 11, 10, and 27 KEGG pathways, including 5 cellular and humoral immune pathways containing apoptosis, autophagy, endocytosis, MAPK, Toll, and Imd signaling pathways. Furthermore, complex competing endogenous RNA (ceRNA) regulatory networks were detected to be formed among DEcircRNAs, DEmiRNAs, and DEmRNAs. The Target DEmRNAs were engaged in 24, 20, and 25 functional terms as well as 62, 80, and 159 pathways, including several vital immune defense-associated pathways, namely the lysosome, endocytosis, phagosome, autophagy, apoptosis, MAPK, Jak-STAT, Toll, and Imd signaling pathways. Finally, back-splicing sites within 15 circRNAs and the difference in the 9 DEcircRNAs' expression between un-inoculated and A. apis-inoculated larval guts were confirmed utilizing molecular methods. These findings not only enrich our understanding of bee host-fungal pathogen interactions, but also lay a foundation for illuminating the mechanism underlying the DEcircRNA-mediated immune defense of A. m. ligustica larvae against A. apis invasion.

Tetrahedron supported CRISPR/Cas13a cleavage for electrochemical detection of circular RNA in bladder cancer

As a diagnostic biomarker, the detection of circular RNA (circRNA) is vital for the early screening of bladder cancer. Usually, the low abundance of circRNA in clinic samples results in the necessarily complicated extraction before detection. In this work, a tetrahedron supported CRISPR/Cas13a cleavage has been explored for direct electrochemical detection of circRNA in urine from bladder cancer. CRISPR/Cas13a system has been reasonably designed to recognize the characteristic back-splice junction site of circRNA. The activated CRISPR/Cas13a by circRNA can cleave uracil bases composed of DNA tetrahedron immobilized on the surface of gold electrode, resulting in the breakage of DNA tetrahedron and the release of electrochemical active molecule methylene blue. By virtue of highly catalytic efficiency of CRISPR/Cas13a and rigid structure of tetrahedron, the developed electrochemical biosensor can directly detect circRNA in 25 μL urine sample with the lowest detection limit of 0.089 fM and linear detection range from 10 fM to 50 nM in less than 10 min, so as to avoid complicated extraction process and benefit for on-site detection.

CircRNA-regulated immune responses of asian honey bee workers to microsporidian infection.

Nosema ceranae is a widespread fungal parasite for honey bees, causing bee nosemosis. Based on deep sequencing and bioinformatics, identification of circular RNAs (circRNAs) in Apis cerana workers’ midguts and circRNA-regulated immune response of host to N. ceranae invasion were conducted in this current work, followed by molecular verification of back-splicing sites and expression trends of circRNAs. Here, 10185 and 7405 circRNAs were identified in the midguts of workers at 7 days (AcT1) and 10 days (AcT2) post inoculation days post-inoculation with N. ceranae. PCR amplification result verified the back-splicing sites within three specific circRNAs (novel_circ_005123, novel_circ_007177, and novel_circ_015140) expressed in N. ceranae-inoculated midgut. In combination with transcriptome data from corresponding un-inoculated midguts (AcCK1 and AcCK2), 2266 circRNAs were found to be shared by the aforementioned four groups, whereas the numbers of specific ones were 2618, 1917, 5691, and 3723 respectively. Further, 83 52) differentially expressed circRNAs (DEcircRNAs) were identified in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group. Source genes of DEcircRNAs in workers’ midgut at seven dpi were involved in two cellular immune-related pathways such as endocytosis and ubiquitin mediated proteolysis. Additionally, competing endogenous RNA (ceRNA) network analysis showed that 23 13) DEcircRNAs in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group could target 18 14) miRNAs and further link to 1111 (1093) mRNAs. These target mRNAs were annotated to six cellular immunity pathways including endocytosis, lysosome, phagosome, ubiquitin mediated proteolysis, metabolism of xenobiotics by cytochrome P450, and insect hormone biosynthesis. Moreover, 284 164) internal ribosome entry site and 54 26) ORFs were identified from DEcircRNAs in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group; additionally, ORFs in DEcircRNAs in midgut at seven dpi with N. ceranae were associated with several cellular immune pathways including endocytosis and ubiquitin-mediated proteolysis. Ultimately, RT-qPCR results showed that the expression trends of six DEcircRNAs were consistent with those in transcriptome data. These results demonstrated that N. ceranae altered the expression pattern of circRNAs in A. c. cerana workers’ midguts, and DEcircRNAs were likely to regulate host cellular and humoral immune response to microsporidian infection. Our findings lay a foundation for clarifying the mechanism underlying host immune response to N. ceranae infection and provide a new insight into interaction between Asian honey bee and microsporidian.

Expression profiling and functional analysis of circular RNAs in vitro model of intermittent hypoxia-induced liver injury.

Intermittent hypoxia (IH) is a prominent feature of obstructive sleep apnea (OSA) which is increasingly recognized as a key risk factor for liver injury. Circular RNAs (circRNAs) has been suggested to act as a regulator of multiple biological processes. However, there is no study evaluating circRNAs alterations and potential role of circRNAs in OSA-related liver injury. The present study aimed to investigate circRNA expression profiles in vitro model of IH-induced liver injury, as well as potential functional characterization of the differentially expressed circRNAs (DE circRNAs). BRL-3A cells were exposed to IH or normoxia. Cell apoptosis and cell viability were evaluated using flow cytometry and cell counting kit-8, respectively. The expression profile of circRNAs was depicted by circRNA sequencing. The selected circRNAs were verified by quantitative real-time PCR (qRT-PCR). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) analyses were employed to predict DE circRNAs functions. The circRNA-miRNA-mRNA regulatory network was constructed. IH treatment caused cell injury in BRL-3A cells. 98 circRNAs were identified as being dysregulated in IH-treated BRL-3A cells. Among them, 58 were up-regulated and 40 were down-regulated. Go and KEGG analyses suggested that the DE circRNAs were predominantly enriched in the biological process such as positive regulation of NF−kappaB transcription factor activity and pathways such as circadian entrainment, Wnt signaling pathway, MAPK signaling pathway, and protein export. 3 up-regulated circRNAs and 3 down-regulated circRNAs with high number of back-splicing sites were chosen for qRT-PCR validation and were consistent with the sequencing data. CircRNA1056 and circRNA805 were predicted to interact with microRNAs that might thereby regulate downstream genes. The study characterized a profile of dysregulated circRNAs in IH-induced BRL-3A cell injury. DE circRNAs may play vital roles in the pathophysiology of IH-induced liver injury. Our findings provide preliminary support for further research in mechanisms and a new theory for the pathogenesis of OSA-related liver injury.

Prediction of Back-splicing sites for CircRNA formation based on convolutional neural networks

BackgroundCircular RNAs (CircRNAs) play critical roles in gene expression regulation and disease development. Understanding the regulation mechanism of CircRNAs formation can help reveal the role of CircRNAs in various biological processes mentioned above. Back-splicing is important for CircRNAs formation. Back-splicing sites prediction helps uncover the mysteries of CircRNAs formation. Several methods were proposed for back-splicing sites prediction or circRNA-realted prediction tasks. Model performance was constrained by poor feature learning and using ability.ResultsIn this study, CircCNN was proposed to predict pre-mRNA back-splicing sites. Convolution neural network and batch normalization are the main parts of CircCNN. Experimental results on three datasets show that CircCNN outperforms other baseline models. Moreover, PPM (Position Probability Matrix) features extract by CircCNN were converted as motifs. Further analysis reveals that some of motifs found by CircCNN match known motifs involved in gene expression regulation, the distribution of motif and special short sequence is important for pre-mRNA back-splicing.ConclusionsIn general, the findings in this study provide a new direction for exploring CircRNA-related gene expression regulatory mechanism and identifying potential targets for complex malignant diseases. The datasets and source code of this study are freely available at: https://github.com/szhh521/CircCNN.

Circular RNA circPOLR2A promotes clear cell renal cell carcinoma progression by facilitating the UBE3C-induced ubiquitination of PEBP1 and, thereby, activating the ERK signaling pathway

BackgroundIncreasing evidence has demonstrated that circular RNAs (circRNAs) are implicated in cancer progression. However, the aberrant expression and biological functions of circRNAs in clear cell renal cell carcinoma (cRCC) remain largely elusive.MethodDifferentially expressed circRNAs in cRCC were filtered via bioinformatics analysis. Aberrant circPOLR2A expression was validated in cRCC tissues and cell lines via qRT-PCR. Sanger sequencing was used to identify the backsplicing site of circPOLR2A. In vitro and in vivo functional experiments were performed to evaluate the role of circPOLR2A in cRCC malignancy. RNA pull-down, mass spectrometry, RIP, FISH and immunofluorescence assays were used to identify and validate the circPOLR2A-interacting proteins. Ubiquitination modification and interaction between proteins were detected via Co-IP and western blotting. The m6A modification in circPOLR2A was validated by the meRIP assay.ResultsBioinformatics analysis revealed that circPOLR2A was highly expressed in cRCC tissues and metastatic cRCC tissues. CircPOLR2A expression was associated with tumor size and TNM stage in cRCC patients. In vitro and in vivo functional assays revealed that circPOLR2A accelerated cRCC cell proliferation, migration, invasion and angiogenesis, while inhibiting apoptosis. Further mechanistic research suggested that circPOLR2A could interact with UBE3C and PEBP1 proteins, and that UBE3C could act as a specific ubiquitin E3 ligase for the PEBP1 protein. The UBE3C/circPOLR2A/PEBP1 protein-RNA ternary complex enhanced the UBE3C-mediated ubiquitination and degradation of the PEBP1 protein which could inactivate the ERK signaling pathway. Rescue experiments revealed that the PEBP1 protein was the functional downstream target of circPOLR2A. Furthermore, m6A modification in circPOLR2A was confirmed, and the m6A reader YTHDF2 could regulate circPOLR2A expression.ConclusionOur study demonstrated that circPOLR2A modulated the UBE3C-mediated ubiquitination and degradation of the PEBP1 protein, and further activated the ERK pathway during cRCC progression and metastasis. The m6A reader, YTHDF2, regulated circPOLR2A expression in cRCC. Hence, circPOLR2A could be a potential target for the diagnosis and treatment of cRCC.

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

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

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.

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 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.