Key circRNAs Research Articles

Function of hsa_circ_0006646 as a competing endogenous RNA to promote progression in gastric cancer by regulating the miR-665-HMGB1 axis.

Mounting evidences indicate that circular RNAs (circRNAs) are a novel class of non-coding RNAs and play vital roles in the tumorigenesis and aggressiveness including gastric cancer (GC). Nevertheless, the precise functions and underlying mechanisms of circRNAs in GC remain largely unknown. The Gene Expression Omnibus (GEO) data set GSE163416 was analyzed to screen the key circRNAs in GC. hsa_circ_0006646 was chosen for further study. GC tissues and matched adjacent normal gastric mucosal epithelial tissues were obtained from the Fourth Hospital of Hebei Medical University. The expressions of hsa_circ_0006646 was detected using quantitative real-time polymerase chain reaction (qRT-PCR). hsa_circ_0006646 was knocked down to identify its effects on GC cells. Bioinformatics algorithms were analyzed to predict the microRNA (miRNAs) potentially sponged by hsa_circ_0006646 and its target genes. Fluorescence in situ hybridization (FISH) was conducted to determine the subcellular location of hsa_circ_0006646 and the predicted miRNA. Then, qRT-PCR, luciferase reporter assay, radioimmunoprecipitation assay, Western blotting, and miRNA rescue experiments were used to confirm the hsa_circ_0006646-related regulatory axis in GC. Cell Counting Kit-8 (CCK-8), colony formation, wound healing, and Transwell experiments were performed to determine the effect of the hsa_circ_0006646-related regulatory axis on GC cells' malignant behaviors in vitro. The xenograft tumor mouse model was established to evaluate the effect of hsa_circ_0006646 in vivo. hsa_circ_0006646 exhibited a high expression in GC tissues as compared to corresponding adjacent normal gastric mucosal epithelial tissues and its high expression was positively correlated with TNM stage, lymph node invasion and poor prognosis (P<0.05). Knockdown of hsa_circ_0006646 suppressed the proliferation, colony formation, migration, and invasion in GC cells (all P<0.05). hsa_circ_0006646 upregulated high mobility group box 1 (HMGB1) by sponging miR-665 in GC cells (P<0.05). The hsa_circ_0006646-miR-665-HMGB1 axis promoted malignant behaviors and epithelial-mesenchymal transition (EMT) in GC cells by activating the Wnt/β-catenin pathway (P<0.05). The existence of hsa_circ_0006646-miR-665-HMGB1 axis was confirmed in GC specimens (P<0.05). Consequently, down-regulated hsa_circ_0006646 inhibited the progression and EMT of GC cells in vivo (P<0.05). For the first time, we demonstrated that hsa_circ_0006646-miR-665-HMGB1 axis exerted its tumor-promoting effects in GC, which suggested that hsa_circ_0006646 could be potentially targeted for GC treatment.

CircSnd1 promotes atherosclerosis progression through the miR-485-3p/Olr1 signaling pathway

BackgroundCircular RNAs (circRNAs) participate in the development of atherosclerotic cardiovascular disease. Identifying and verifying the key competing endogenous RNA (ceRNA) network related to atherosclerosis (AS) is significant for understanding the development of AS. The aim of this study was to investigate the circRNA-miRNA‒mRNA network, identify a key circRNA and explore its role in the development of atherosclerosis. MethodsDifferentially expressed mRNAs (DEMs) and circRNAs (DECs) in the AS model were obtained from datasets in the Gene Expression Omnibus (GEO) database. R software and Cytoscape software were used to construct and visualize the ceRNA network. The dual-luciferase reporter experiment and the RNA pull-down experiment were used to verify the selected ceRNA axis. siRNA targeting circRNA, miRNA mimic, miRNA inhibitor, or gene overexpression plasmid was used for in vitro functional studies. ELISA and western blotting were used to detect inflammation and lipid transport-related proteins. Furthermore, an AS mouse model was established and treated with recombinant adeno-associated viral vectors to further verify the influence of the selected ceRNA axis on the occurrence and/or development of AS. ResultsA total of 497 DEMs were enriched in 25 pathways, based on which the circ_0082139 (circSnd1)/miR-485-3p/Olr1 axis was selected. In vitro, the interaction among the three molecules of this axis was validated and it was found to affect inflammation and lipid transport, which were characterized by the significant change of inflammatory factors (Il-6, Il-8, Tnf-α, Mcp-1, Vcam-1, and Icam-1), and lipid transport-related genes, including Abca1, Abcg1, Ldlr, Hdlbp, Lp-pla2, and Srebp-1c. Through animal experiments, we further verified that the circSnd1/miR-485-3p/Olr1 axis regulated these molecules and participated in the formation and/or development of AS in vivo. ConclusionsThe circSnd1/miR-485-3p/Olr1 axis participates in the formation and development of atherosclerosis by regulating inflammation and lipid transport.

Genome-wide analysis of circular RNA-mediated ceRNA regulation in porcine skeletal muscle development

BackgroundAs a diverse and abundant class of endogenous RNAs, circular RNAs (circRNAs) participate in various biological processes including cell proliferation and apoptosis. Nevertheless, few researchers have investigated the role of circRNAs in muscle development in cultivated pigs.ResultsIn this study, we used RNA-seq to construct circRNA expression profiles in skeletal muscle of Jinfen White pigs at the age of 1, 90, and 180 days. Among the 16,990 identified circRNAs, 584 circRNAs were differentially expressed. Moreover, the enrichment analysis of DE circRNA host genes showed that they were mainly involved in muscle contraction, muscle organ development and muscle system processes, as well as AMPK and cAMP-related signal pathways. We also constructed a circRNA–miRNA–mRNA co-expression network to find key circRNAs which many involved in the regulation of porcine skeletal muscle development through the competitive endogenous RNA (ceRNA) mechanism. It is noteworthy that circ_0018595/miR-1343/PGM1 axis may play a regulatory role in the development of porcine skeletal muscle.ConclusionsThis study identified the circRNAs and present the circRNA expression profile in the development of pigs, revealed that DE circRNA host genes participate in different cell fates and enriched the porcine ceRNA network. Thus, this work will become a valuable resource for further in-depth study of the regulatory mechanism of circRNA in the development of porcine skeletal muscle.

Construction of circRNA-miRNA-mRNA Network Reveal Functional circRNAs and Key Genes in Acute Myeloid Leukemia.

CircRNA is closely correlated with a wide variety of processes of acute myeloid leukemia (AML), whereas the novel circRNAs, their molecular mechanism and the specific function they played in AML should be explored in depth. The microarray chip data of AML patients and normal samples in the Gene Expression Omnibus (GEO) database were selected to differentially expressed (DE) circRNA, miRNA, and mRNA genes. The miRNA gene was the intersection of the circRNA target gene predicted using CSCD and the miRNA gene screened from AML patients, while the mRNA gene was the intersection of the target gene mRNA of miRNA predicted using miRanda and miRTarBase software and the mRNA gene screened from AML patients. The hub mRNAs related to survival were further screened through Cox proportional hazard regression. CircRNA/miRNA/mRNA interaction network was constructed by using Cytoscape software.10 circRNAs and 6 miRNAs in bone marrow mononuclear cells (BMMNCs) of AML patients (n=43) and healthy controls (n=35) were determined by RT-qPCR. Correlations between them were analyzed by Pearson correlation coefficient. 10 circRNAs, 6 miRNAs, and 33 mRNAs were identified. Subsequently, the network of circRNAs, miRNAs, and hub genes was built using Cystoscope. Four key circRNAs, seven hub genes and their regulatory pathways were identified. The result of RT-qPCRs showed that hsa_circ_0009581 and hsa_circ_0005273 were significantly upregulated in AML patients while hsa_circ_0000497 and hsa_circ_0001947 were significantly downregulated. Hsa-miR-150-5p was significantly downregulated; hsa-miR-454-3p was upregulated in AML patients. Hsa_circ_0009581 and hsa-miR-150-5p; hsa_ circ_0001947 and hsa-miR-454-3p were inversely correlated using Pearson's correlation coefficient. This study suggests that differentially expressed circRNAs take on a critical significance to AML development and may be the effective therapeutic targets. We suppose that hsa_circ_0009581 promotes leukemia development through hsa-miR-150-5p and hsa_circ_0001947 through hsa-miR-454-3p. hsa_circ_0001947 and hsa_circ_0009581 may provide new directions in the pathogenesis of AML.

Bioinformatic analysis confirms differences in circular RNA expression profiles of cumulus cells between patients with ovarian and peritoneal endometriosis-associated infertility.

Endometriosis has a detrimental effect on oocyte quality, and ovarian endometriosis (OEM) and peritoneal endometriosis (PEM) may have different effects on female fertility. Therefore, we conducted a study to explore the circular RNA (circRNA) expression profiles of cumulus cells (CCs) in patients with OEM (n = 3), PEM (n = 3), and tubal factor infertility (TFI, n = 3) using high-throughput sequencing techniques and attempted to identify common and unique circRNAs in the OEM and PEM groups. The CIRCexplorer2 program was used to identify circRNAs. Seven candidate circRNAs were validated in 30 samples using quantitative real-time polymerase chain reaction (qRT-PCR). Finally, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to annotate the function of circRNA-targeted genes, which were verified by sequencing results and constructed circRNA-miRNA-mRNA networks. A total of 11833 circRNAs were identified in nine samples. The numbers of differentially expressed circRNAs between the OEM and TFI groups, PEM and TFI groups, and OEM and PEM groups were 130, 71, and 191, respectively. After taking intersections, 11 circRNAs were considered common circRNAs in the OEM and PEM groups; 39 circRNAs in the OEM group and 17 circRNAs in the PEM group were identified as unique key circRNAs. During qRT-PCR validation, hsa_circ_0003638 was significantly upregulated in the PEM group compared to that in the OEM and TFI groups. Functional analysis of circRNA-targeted genes revealed that apoptosis, PI3K-AKT, and p53 signaling pathways were enriched in the PEM-TFI comparison groups, whereas the functions of target genes involved in the JAK-STAT and TGF-β signaling pathways were enriched in the PEM-OEM comparison groups. Our findings confirmed differences in circRNA expression profiles of CCs between patients with OEM and PEM infertility and provide new insights into the different effects of various endometriosis phenotypes on oocytes.

Systematic identification of potential key microRNAs and circRNAs in the dorsal root ganglia of mice with sciatic nerve injury.

Neuropathic pain (NeP) is a pathological condition arising from a lesion or disease affecting the somatosensory system. Accumulating evidence has shown that circular RNAs (circRNAs) exert critical functions in neurodegenerative diseases by sponging microRNAs (miRNAs). However, the functions and regulatory mechanisms of circRNAs as competitive endogenous RNAs (ceRNAs) in NeP remain to be determined. The sequencing dataset GSE96051 was obtained from the public Gene Expression Omnibus (GEO) database. First, we conducted a comparison of gene expression profiles in the L3/L4 dorsal root ganglion (DRG) of sciatic nerve transection (SNT) mice (n = 5) and uninjured mice (Control) (n = 4) to define the differentially expressed genes (DEGs). Then, critical hub genes were screened by exploring protein-protein interaction (PPI) networks with Cytoscape software, and the miRNAs bound to them were predicted and selected and then validated by qRT-PCR. Furthermore, key circRNAs were predicted and filtered, and the network of circRNA-miRNA-mRNA in NeP was constructed. A total of 421 DEGs were identified, including 332 upregulated genes and 89 downregulated genes. Ten hub genes, including IL6, Jun, Cd44, Timp1, and Csf1, were identified. Two miRNAs, mmu-miR-181a-5p and mmu-miR-223-3p, were preliminarily verified as key regulators of NeP development. In addition, circARHGAP5 and circLPHN3 were identified as key circRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis demonstrated that these differentially expressed mRNAs and targeting miRNAs were involved in signal transduction, positive regulation of receptor-mediated endocytosis and regulation of neuronal synaptic plasticity. These findings have useful implications for the exploration of new mechanisms and therapeutic targets for NeP. These newly identified miRNAs and circRNAs in networks reveal potential diagnostic or therapeutic targets for NeP.

METTL3-mediated m6A modification of has_circ_0007905 promotes age-related cataract progression through miR-6749-3p/EIF4EBP1.

Many cases of blindness are caused by age-related cataracts (ARCs). N6-methyladenosine (m6A)-modified circRNA widely participates in disease progression. However, the role of m6A modification of circRNA in ARC is unclear. We mined and elucidated the functions and mechanisms of key circRNAs with m6A modification involved in ARC progression. The GSE153722 dataset was used to mine m6A-mediated key circRNA. Loss-of-function assays and rescue assays were used to explore the effect and mechanism of circRNA on ARC cell proliferation and apoptosis. Has_circ_0007905 was a hypermethylated and upregulated expression in the ARC group relative to the control group both in vivo and in vitro. Silencing of has_circ_0007905 promoted proliferation and inhibited the apoptosis of HLE-B3 cells. METTL3 was upregulated in HLE-B3 cells after ARC modeling and had four binding sites with has_circ_0007905 and a mediated m6A modification of has_circ_0007905. Proliferation was significantly inhibited and apoptosis of HLE-B3 cells was facilitated by METTL3 overexpression, whereas these effects were prevented by has_circ_0007905 silencing. Silencing of has_circ_0007905 led to an alteration in the transcriptome landscape. Differentially expressed genes were mainly involved in immune-related processes and pathways. EIF4EBP1 overexpression promoted apoptosis and suppressed proliferation, and also significantly reversed effects of has_circ_0007905 silencing. Moreover, miR-6749-3p significantly decreased the luciferase activities of wild type plasmids with both of has_circ_0007905 and EIF4EBP1. MiR-6749-3p inhibitor blocked elevation in proliferation and reduced EIF4EBP1 expression and apoptosis conferred by has_circ_0007905 silencing. We reveal for the first time that the commitment of ARC progression is guided by METTL3/has_circ_0007905/miR-6749-3p/EIF4EBP1 axis, and the results provide new insights into ARC pathology.

BMSC-Derived Exosomal CircHIPK3 Promotes Osteogenic Differentiation of MC3T3-E1 Cells via Mitophagy

Exosome-based therapy is emerging as a promising strategy to promote bone regeneration due to exosomal bioactive cargos, among which circular RNA (circRNA) has recently been recognized as the key effector. The role of exosomal circRNA derived from bone marrow mesenchymal stem cells (BMSCs) has not been well-defined. The present study aimed to clarify the regulatory function and molecular mechanism of BMSC-derived exosomal circRNA in osteogenesis. Exosomes derived from bone marrow mesenchymal stem cells (BMSC-Exos) were isolated and identified. BMSC-Exos' pro-osteogenic effect on MC3T3-E1 cells was validated by alkaline phosphatase (ALP) activity and Alizarin Red staining. Through bioinformatic analysis and molecular experiments, circHIPK3 was selected and verified as the key circRNA of BMSC-Exos to promote osteoblast differentiation of MC3T3-E1 cells. Mechanistically, circHIPK3 acted as an miR-29a-5p sponge and functioned in mitophagy via targeting miR-29a-5p and PINK1. Additionally, we showed that the mitophagy level of MC3T3-E1 cells were mediated by BMSC-Exos, which promoted the osteogenic differentiation. Collectively, our results revealed an important role for BMSC-derived exosomal circHIPK3 in osteogenesis. These findings provide a potentially effective therapeutic strategy for bone regeneration.

Key circRNAs from goat: discovery, integrated regulatory network and their putative roles in the differentiation of intramuscular adipocytes

BackgroundThe procession of preadipocytes differentiation into mature adipocytes involves multiple cellular and signal transduction pathways. Recently. a seirces of noncoding RNAs (ncRNAs), including circular RNAs (circRNAs) were proved to play important roles in regulating differentiation of adipocytes.ResultIn this study, we aimed to identificate the potential circRNAs in the early and late stages of goat intramuscular adipocytes differentiation. Using bioinformatics methods to predict their biological functions and map the circRNA-miRNA interaction network. Over 104 million clean reads in goat intramuscular preadipocytes and adipocytes were mapped, of which16 circRNAs were differentially expressed (DE-circRNAs). Furthermore, we used real-time fluorescent quantitative PCR (qRT-PCR) technology to randomly detect the expression levels of 8 circRNAs among the DE-circRNAs, and our result verifies the accuracy of the RNA-seq data. From the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the DE-circRNAs, two circRNAs, circ_0005870 and circ_0000946, were found in Focal adhesion and PI3K-Akt signaling pathway. Then we draw the circRNA-miRNA interaction network and obtained the miRNAs that possibly interact with circ_0005870 and circ_0000946. Using TargetScan, miRTarBase and miR-TCDS online databases, we further obtained the mRNAs that may interact with the miRNAs, and generated the final circRNA-miRNA-mRNA interaction network. Combined with the following GO (Gene Ontology) and KEGG enrichment analysis, we obtained 5 key mRNAs related to adipocyte differentiation in our interaction network, which are FOXO3(forkhead box O3), PPP2CA (protein phosphatase 2 catalytic subunit alpha), EEIF4E (eukaryotic translation initiation factor 4), CDK6 (cyclin dependent kinase 6) and ACVR1 (activin A receptor type 1).ConclusionsBy using Illumina HiSeq and online databases, we generated the final circRNA-miRNA-mRNA interaction network that have valuable functions in adipocyte differentiation. Our work serves as a valuable genomic resource for in-depth exploration of the molecular mechanism of ncRNAs interaction network regulating adipocyte differentiation.

Key circRNAs, lncRNAs, and mRNAs of ShenQi Compound in Protecting Vascular Endothelial Cells From High Glucose Identified by Whole Transcriptome Sequencing.

Vascular endothelial cells, which make up the inner wall of blood arteries, are susceptible to damage from oxidative stress and apoptosis caused by hyperglycemia. According to certain reports, noncoding RNAs are involved in controlling oxidative stress and apoptosis. ShenQi Compound (SQC), a traditional herbal remedy, has been successfully treating diabetic vascular disease in China for more than 20 years. Although it is well established that SQC protects the vascular endothelium, the molecular mechanism remains unknown. Goto-Kakizaki rats, spontaneous type II diabetes rats, that consistently consume a high-fat diet were chosen as model animals. Six groups (control group, model group, metformin group, and 7.2 g/kg/d SQC group, 14.4 g/kg/d SQC group, and 28.8 g/kg/d SQC group) were included in this work, 15 rats each group. The approach of administration was gavage, and the same volume (5.0 mL/kg/d) was given in each group, once a day, 12 weeks. The thoracic aortas were removed after the rats were sacrificed. Oxidative reduction profile in thoracic aorta, histopathological observation of thoracic aorta, endothelial cell apoptosis in thoracic aorta, whole transcriptome sequencing, bioinformatic analyses, and qRT-PCR were conducted. As a result, SQC prevented the oxidative stress and apoptosis induced by a high glucose concentration. Under hyperglycemia condition, noncoding RNAs, including 1 downregulated novel circRNA (circRNA.3121), 3 downregulated lncRNAs (Skil.cSep08, Shawso.aSep08-unspliced, and MSTRG.164.2), and 1 upregulated mRNA (Pcdh17), were clearly reverse regulated by SQC. SQC plays a role in protecting vascular endothelial cells from high glucose mainly by mediating ncRNA to inhibit cell apoptosis and oxidative stress.

Profiling of differentially expressed circRNAs and functional prediction in peripheral blood mononuclear cells from patients with rheumatoid arthritis

Background Rheumatoid arthritis (RA) is a chronic autoimmune disease associated with an increased risk of death, but its underlying mechanisms are not fully understood. Circular RNAs (circRNAs) have recently been implicated in various biological processes. The aim of this study was to investigate the key circRNAs related to RA. Methods A microarray assay was used to identify the differentially expressed circRNAs (DEcircRNAs) in peripheral blood mononuclear cells (PBMCs) from patients with RA compared to patients with osteoarthritis (OA) and healthy controls. Then, quantitative real-time PCR was applied to verify the DEcircRNAs, and correlations between the levels of DEcircRNAs and laboratory indices were analysed. We also performed extensive bioinformatic analyses including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genome (KEGG) pathway and potential circRNA–miRNA–mRNA network analyses to predict the function of these DEcircRNAs. Results A total of 35,342 and 6146 DEcircRNAs were detected in RA patients compared to controls and OA patients, respectively. Nine out of the DEcircRNAs in RA were validated by real-time PCR. There were correlations between the levels of DEcircRNAs and some of the laboratory indices. GO analyses revealed that these DEcircRNAs in RA were closely related to cellular protein metabolic processes, gene expression, the immune system, cell cycle, posttranslational protein modification and collagen formation. Functional annotation of host genes of these DEcircRNAs was implicated in several significantly enriched pathways, including metabolic pathways, ECM–receptor interaction, the PI3K–Akt signalling pathway, the AMPK signalling pathway, leukocyte transendothelial migration, platelet activation and the cAMP signalling pathway, which might be responsible for the pathophysiology of RA. Conclusions The findings of this study may help to elucidate the role of circRNAs in the specific mechanism underlying RA. Key messages Microarray assays showed that a total of 35,342 and 6146 DEcircRNAs were detected in RA patients compared to controls and OA patients, respectively. Nine out of the DEcircRNAs in RA were validated by real-time PCR, and the levels of the DEcircRNAs were related to some of the laboratory indices. GO analyses revealed that the DEcircRNAs in RA were closely related to cellular protein metabolic processes, gene expression, the immune system, etc. Functional annotation of host genes of the DEcircRNAs in RA was implicated in several significantly enriched pathways, including metabolic pathways, ECM–receptor interaction, the PI3K–Akt signalling pathway, etc.

Identification of pathological-related and diagnostic potential circular RNAs in Stanford type A aortic dissection.

Stanford type A aortic dissection (TAAD) is one of the lethal macrovascular diseases caused by the invasion of blood into the media layer of ascending aortic wall. Inflammation, smooth muscle dysfunction, and extracellular matrix (ECM) degradation were regarded as the major pathology in affected tissue. However, the expression pattern and its regulation especially through circular RNAs (circRNAs) as an overall characteristic of TAAD molecular pathology remain unclear. We employed CIRCexplorer2 to identify circRNAs based on the RNA sequencing (RNA-seq) data of human ascending aortic tissues to systematically assess the role of circRNA in the massive alterations of gene expression in TAAD aortas. The key circRNAs were determined by LASSO model and functionally annotated by competing endogenous RNAs (ceRNA) network and co-analysis with mRNA profile. The expression level and diagnostic capability of the 4 key circRNAs in peripheral serum were confirmed by real-time polymerase chain reaction (RT-PCR). The 4 key circRNAs, namely circPTGR1 (chr9:114341075-114348445[-]), circNOX4 (chr11:89069012-89106660[-]), circAMN1 (chr12:31854796-31862359[-]) and circUSP3 (chr15:63845913-63855207[+]), demonstrated a high power to discriminate between TAAD and control tissues, suggesting that these molecules stand for a major difference between the tissues at gene regulation level. Functionally, the ceRNA network of circRNA-miRNA-mRNA predicted by the online databases, combining gene set enrichment analysis (GSEA) and cell component prediction, revealed that the identified circRNAs covered all the aspects of primary TAAD pathology, centralized with increasing inflammatory factors and cells, and ECM destruction and loss of vascular inherent cells along with the circRNAs. Importantly, we validated the high concentration and diagnostic capability of the 4 key circRNAs in the peripheral serum in TAAD patients. This study reinforces the vital status of circRNAs in TAAD and the possibility of serving as promising diagnostic biomarkers.

Evolutionary Landscape of Tea Circular RNAs and Its Contribution to Chilling Tolerance of Tea Plant.

Chilling stress threatens the yield and distribution pattern of global crops, including the tea plant (Camellia sinensis), one of the most important cash crops around the world. Circular RNA (circRNA) plays roles in regulating plant growth and biotic/abiotic stress responses. Understanding the evolutionary characteristics of circRNA and its feedbacks to chilling stress in the tea plant will help to elucidate the vital roles of circRNAs. In the current report, we systematically identified 2702 high-confidence circRNAs under chilling stress in the tea plant, and interestingly found that the generation of tea plant circRNAs was associated with the length of their flanking introns. Repetitive sequences annotation and DNA methylation analysis revealed that the longer flanking introns of circRNAs present more repetitive sequences and higher methylation levels, which suggested that repeat-elements-mediated DNA methylation might promote the circRNAs biogenesis in the tea plant. We further detected 250 differentially expressed circRNAs under chilling stress, which were functionally enriched in GO terms related to cold/stress responses. Constructing a circRNA-miRNA-mRNA interaction network discovered 139 differentially expressed circRNAs harboring potential miRNA binding sites, which further identified 14 circRNAs that might contribute to tea plant chilling responses. We further characterized a key circRNA, CSS-circFAB1, which was significantly induced under chilling stress. FISH and silencing experiments revealed that CSS-circFAB1 was potentially involved in chilling tolerance of the tea plant. Our study emphasizes the importance of circRNA and its preliminary role against low-temperature stress, providing new insights for tea plant cold tolerance breeding.

Analyses of lncRNAs, circRNAs, and the Interactions between ncRNAs and mRNAs in Goat Submandibular Glands Reveal Their Potential Function in Immune Regulation.

As part of one of the main ruminants, goat salivary glands hardly secrete digestive enzymes, but play an important role in immunity. The immune function of goat salivary glands significantly changes with age, while the expression profile and specific function of non-coding RNA during this process are unknown. In this study, transcriptome sequencing was performed on submandibular gland (SMG) tissues of 1-month-old, 12-month-old, and 24-month-old goats, revealing the expression patterns of lncRNA and circRNA at different ages. A total of 369 lncRNAs and 1699 circRNAs were found to be differentially expressed. Functional enrichment analyses showed that the lncRNA regulated target mRNAs and circRNA host genes were significantly enriched in immune-related GO terms and pathways. CeRNA network analysis showed that the key differentially expressed circRNAs and lncRNAs mainly regulate the key immune-related genes ITGB2, LCP2, PTPRC, SYK, and ZAP70 through competitive binding with miR-141-x, miR-29-y, and chi-miR-29b-3p, thereby affecting the natural killer cell-mediated cytotoxicity pathway, the T cell receptor signaling pathway, and other immune-related pathways. It should be noted that the expression of key circRNAs, lncRNAs, and key immune-related genes in goat SMGs decreased significantly with the growth of the goat. This is the first reporting of lncRNAs, circRNAs, and ceRNA network regulation in goat SMGs. Our study contributes to the knowledge of changes in the expression of non-coding RNAs during SMG development in goats and provides new insights into the relationship between non-coding RNAs and salivary gland immune function in goats.

Increased transcription of hsa_circ_0000644 upon RUNX family transcription factor 3 downregulation participates in the malignant development of bladder cancer

BackgroundStudies are ongoing to examine the versatile functions of circular RNAs (circRNAs) in human diseases. This research investigates the effects of hsa_circ_0000644 (circ_644) and its related molecules on the malignant behavior of bladder cancer (BCa) cells. MethodsAbundant bioinformatics analyses were performed to screen the key circRNA and its related molecules in BCa. Tumor tissues and the para-tumorous tissues were collected from 58 patients with BCa. Expression of RUNX family transcription factor 3 (RUNX3), circ_644, microRNA-143-3p (miR-143-3p), and musashi RNA binding protein 2 (MSI2) in BCa tissues or cells was determined. Molecular interactions were confirmed by chromatin immunoprecipitation, RNA pull-down, and luciferase assays. Gain and loss-of function assays were performed using two BCa cell lines (T24 and HT1376). ResultsCirc_644 was highly expressed whereas RUNX3, which could suppress circ_644 transcription, was lowly expressed in BCa tissues and cells. Upregulation of RUNX3 suppressed proliferation, colony formation, migration and invasion, and tumorigenicity of BCa cells and induced cell cycle arrest. However, the tumor-suppressive effects of RUNX3 were blocked by circ_644 upregulation. Circ_644 served as a sponge for miR-143-3p, and miR-143-3p bound to MSI2 mRNA. The rescue experiments showed that miR-143-3p inhibition or MSI2 overexpression restored the malignant behaviors of BCa cells induced by circ_644 knockdown or RUNX3 overexpression. ConclusionThis study demonstrates that transcriptional activation of circ_644 upon RUNX3 downregulation drives the malignant development of BCa through the miR-143-3p/MSI2 axis. RUNX3 restoration or specific inhibition of circ_644 or MSI2 may help block BCa progression.

Excavation and characterization of key circRNAs for milk fat percentage in Holstein cattle.

Milk fat percentage is one of the significant indicators governing the price and quality of milk and is regulated by a variety of non-coding RNAs. We used RNA sequencing (RNA-seq) techniques and bioinformatics approaches to explore potential candidate circular RNAs (circRNAs) regulating milk fat metabolism. After analysis, compared with low milk fat percentage (LMF) cows, 309 circRNAs were significantly differentially expressed in high milk fat percentage (HMF) cows. Functional enrichment and pathway analysis revealed that the main functions of the parental genes of differentially expressed circRNAs (DE-circRNAs) were related to lipid metabolism. We selected four circRNAs (Novel_circ_0000856, Novel_circ_0011157, novel_circ_0011944, and Novel_circ_0018279) derived from parental genes related to lipid metabolism as key candidate DE-circRNAs. Their head-to-tail splicing was demonstrated by linear RNase R digestion experiments and Sanger sequencing. However, the tissue expression profiles showed that only Novel_circ_0000856, Novel_circ_0011157, and Novel_circ_0011944 were expressed with high abundance in breast tissue. Based on the subcellular localization found that Novel_circ_0000856, Novel_circ_0011157, and Novel_circ_0011944 mainly function as competitive endogenous RNAs (ceRNAs) in the cytoplasm. Therefore, we constructed their ceRNA regulatory networks, and the five hub target genes (CSF1, TET2, VDR, CD34, and MECP2) in ceRNAs were obtained by CytoHubba and MCODE plugins in Cytoscape, as well as tissue expression profiles analysis of target genes. These genes play a key role as important target genes in lipid metabolism, energy metabolism, and cellular autophagy. The Novel_circ_0000856, Novel_circ_0011157, and Novel_circ_0011944 regulate the expression of hub target genes through interaction with miRNAs and constitute key regulatory networks that may be involved in milk fat metabolism. The circRNAs obtained in this study may act as miRNA sponges and thus influence mammary gland development and lipid metabolism in cows, which improves our understanding of the role of circRNAs in cow lactation.

Genome-wide characterization of circRNA expression profile in overexpression of RIP2 chicken macrophages associated with avian pathogenic E.coli infection

ABSTRACT Avian pathogenic E. coli (APEC) can cause localized and systemic diseases in poultry, threatening human health via meat or egg contamination and resulting in considerable economic losses to the poultry industry globally. Increasing evidence shows circRNAs were widely involved in various biological processes. However, the role of circRNAs in the host response against APEC infection, especially correlated with the regulation of RIP2, remains unclear. Herein, the RNAseq technology was used to identify the circRNA expression profiles in the overexpression of RIP2 macrophages with or without APEC infection. A total of 256 and 287 differentially expressed (DE) circRNAs were identified in the overexpression of RIP2 group (oeRIP2) vs. the wild-type group (WT) and oeRIP2 + APEC vs. APEC, respectively, whose parental genes were involved in MAPK signalling pathway, Wnt signalling pathway, focal adhesion, tight junction, and VEGF signalling pathways. Specifically, the key circRNAs, such as 5:814443-825127, 10:18922360-18928461, 2:8746306-8750639, and 2:124177751-124184063 might play a critical role in APEC infection and the regulation of RIP2. As a whole, these findings will facilitate understanding the molecular mechanism underlying circRNAs, especially related to the regulation of the RIP2 gene. Meanwhile, the study may offer new ideas to improve host immune and inflammatory response against APEC infection.

Bioinformatics-based study reveals that AP2M1 is regulated by the circRNA-miRNA-mRNA interaction network and affects Alzheimer's disease.

Alzheimer's disease (AD) is a progressive neurological disease that worsens with time. The hallmark illnesses include extracellular senile plaques caused by β-amyloid protein deposition, neurofibrillary tangles caused by tau protein hyperphosphorylation, and neuronal loss accompanying glial cell hyperplasia. Noncoding RNAs are substantially implicated in related pathophysiology, according to mounting data. However, the function of these ncRNAs is mainly unclear. Circular RNAs (circRNAs) include many miRNA-binding sites (miRNA response elements, MREs), which operate as miRNA sponges or competing endogenous RNAs (ceRNAs). The purpose of this study was to look at the role of circular RNAs (circRNAs) and microRNAs (miRNAs) in Alzheimer's disease (AD) as possible biomarkers. The Gene Expression Omnibus (GEO) database was used to obtain an expression profile of Alzheimer's disease patients (GSE5281, GSE122603, GSE97760, GSE150693, GSE1297, and GSE161435). Through preliminary data deletion, 163 genes with significant differences, 156 miRNAs with significant differences, and 153 circRNAs with significant differences were identified. Then, 10 key genes, led by MAPT and AP2M1, were identified by the mediation center algorithm, 34 miRNAs with obvious prognosis were identified by the cox regression model, and 16 key circRNAs were selected by the database. To develop competitive endogenous RNA (ceRNA) networks, hub circRNAs and mRNAs were used. Finally, GO analysis and clinical data verification of key genes were carried out. We discovered that a down-regulated circRNA (has_circ_002048) caused the increased expression of numerous miRNAs, which further inhibited the expression of a critical mRNA (AP2M1), leading to Alzheimer's disease pathology. The findings of this work contribute to a better understanding of the circRNA-miRNA-mRNA regulating processes in Alzheimer's disease. Furthermore, the ncRNAs found here might become novel biomarkers and potential targets for the development of Alzheimer's drugs.

Deregulated Expression of Circular RNAs Is Associated with Immune Evasion and Leukemia Relapse after Allogeneic Hematopoietic Stem Cell Transplantation.

Circular RNAs (circRNAs) are a novel class of epigenetic regulators that participate in leukemogenesis. However, their roles in leukemia relapse after transplantation remain unclear. We defined the circRNAs profile of the bone-marrow-enriched CD34+ cells from ten acute myeloid leukemia (AML) patients after transplantation (five relapse [RE] and five continuous complete remission [CR]) and four healthy controls (HCs) by RNA-seq. Differentially expressed circRNAs were validated using real-time quantitative polymerase chain reaction (RT-qPCR) in an independent cohort of six AML patients with pairwise samples at diagnosis and at relapse and six controls. The bioinformatics analysis revealed a distinct circRNAs profile in relapse patients compared with controls (CR or HCs), while there was no significant difference between CR and HCs. Functional enrichment analysis demonstrated that mRNAs co-expressed with identified circRNAs were primarily involved in immune-related pathways, including the T cell receptor signaling pathway and lymphocyte differentiation. Moreover, we performed a protein-protein interaction network based on the immune-related genes and annotated 20 hub genes. The abnormal expression of hub genes was responsible for impairing T cell co-stimulation and activation, thus contributing to the immune escape of relapse blasts. We further constructed competing endogenous RNAs (ceRNA) regulatory networks based on immune-related genes and identified 10 key circRNAs that are associated with immune evasion. Six candidate circRNAs and their associated miRNA/mRNAs in the ceRNA network were randomly selected to be validated in another set by RT-qPCR. CircRNAs dysregulation may be involved in the immune evasion of relapse blasts and is associated with AML relapse. Our results identify several promising biomarkers and might provide novel insights into the biology of AML relapse post-transplantation.

The comprehensive detection of miRNA and circRNA in the regulation of intramuscular and subcutaneous adipose tissue of Laiwu pig

circRNAs, as miRNA sponges, participate in many important biological processes. However, it remains unclear whether circRNAs can regulate lipid metabolism. This study aimed to explore the competing endogenouse RNA (ceRNA) regulatory network that affects the difference between intramuscular fat (IMF) and subcutaneous fat (SCF) deposition, and to screen key circRNAs and their regulatory genes. In this experiment, we identified 265 differentially expressed circRNAs, of which 187 up-regulated circRNA and 78 down-regulated circRNA in IMF. Subsequently, we annotated the function of DEcircRNA's host genes, and found that DEcircRNA's host genes were mainly involved in GO terms (including cellular response to fatty acids, lysophosphatidic acid acyltransferase activity, R-SMAD binding, etc.) and signaling pathways (fatty acid biosynthesis, Citrate cycle, TGF- β Signal pathway) related to adipogenesis, differentiation and lipid metabolism. By constructing a circRNA-miRNA network, we screened out DEcircRNA that can competitively bind to more miRNAs as key circRNAs (circRNA_06424 and circRNA_08840). Through the functional annotation of indirect target genes and protein network analysis, we found that circRNA_06424 affects the expression of PPARD, MMP9, UBA7 and other indirect target genes by competitively binding to miRNAs such as ssc-miR-339-5p, ssc-miR-744 and ssc-miR-328, and participates in PPAR signaling pathway, Wnt signaling pathway, unsaturated fatty acid and other signaling pathways, resulting in the difference of fat deposition between IMF and SCF. This study provide a theoretical basis for further research investigating the differences of lipid metabolism in different adipose tissues, providing potential therapeutic targets for ectopic fat deposition and lipid metabolism diseases.