Construction of a Potentially Functional circRNA-miRNA-mRNA Network in Intervertebral Disc Degeneration by Bioinformatics Analysis.
Ontology highlight
ABSTRACT: Background
The competing endogenous RNA- (ceRNA-) mediated regulatory mechanisms are known to play a pivotal role in intervertebral disc degeneration (IDD). Our research intended to establish a ceRNA regulatory network related to IDD through bioinformatics analyses.Methods
The expression profiles of circRNA, miRNA, and mRNA were obtained from the public Gene Expression Omnibus (GEO) datasets. Then, we use sequence-based bioinformatics methods to select differentially expressed mRNAs (DEmRNAs), microRNAs (DEmiRNAs), or circRNAs (DEcircRNAs) related to IDD. We used ChEA3 to verify the targets of transcription factors (TFs). Then, we used DAVID to annotate the DEmRNAs. Finally, we constructed a potentially circRNA-miRNA-mRNA network related to IDD by predicting in the database (ENCORI, TargetScan, miRecords, miRmap, and circBank).Results
We identified 31 common DEmRNAs by Venn analysis, of which MMP2 was regarded as the key hub genes. Simultaneously, miR-423-5p and miR-185-5p were predicted as the upstream molecules of MMP2. Furthermore, a total of six DEcircRNAs were predicted as the upstream circRNAs of miR-423-5p and miR-185-5p. Then, a potential circRNA-miRNA-mRNA network related to IDD was constructed by bioinformatics analysis.Conclusion
A comprehensive ceRNA regulatory network was constructed, which was found to be significant in IDD progression.
SUBMITTER: Huo Z
PROVIDER: S-EPMC8357516 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Accumulating evidence has indicated that noncoding RNAs are involved in intervertebral disc degeneration (IDD); however, the competing endogenous RNA (ceRNA)‑mediated regulatory mechanisms in IDD remain rarely reported. The present study aimed to comprehensively investigate the alterations in expression levels of circular RNA (circRNA), long noncoding RNA (lncRNA), microRNA (miRNA/miR) and mRNA in the nucleus pulposus (NP) of patients with IDD. In addition, crucial lncRNA/circRNA‑miRNA‑mRNA ceRNA interaction axes were screened using the GSE67567 microarray dataset obtained from the Gene Expression Omnibus database. After data preprocessing, differentially expressed circRNAs (DECs), lncRNAs (DELs), miRNAs (DEMs) or genes (DEGs) between IDD and normal controls were identified using the Linear Models for Microarray data method. A protein‑protein interaction (PPI) network was constructed for DEGs based on protein databases, followed by module analysis. The ceRNA network was constructed based on the interaction between miRNAs and mRNAs, and lncRNAs/circRNAs and miRNAs. The underlying functions of mRNAs were predicted using the Database for Annotation, Visualization and Integrated Discovery database. The present study identified 636 DECs, 115 DELs, 84 DEMs and 1,040 DEGs between patients with IDD and control individuals. PPI network analysis demonstrated that Fos proto‑oncogene, AP‑1 transcription factor subunit (FOS), mitogen‑activated protein kinase 1 (MAPK1), hypoxia inducible factor 1 subunit α (HIF1A) and transforming growth factor β1 (TGFB1) were hub genes and enriched in modules. Metastasis‑associated lung adenocarcinoma transcript 1 (MALAT1)/hsa_circRNA_102348‑hsa‑miR‑185‑5p‑TGFB1/FOS, MALAT1‑hsa‑miR‑155‑5p‑HIF1A, hsa_circRNA_102399‑hsa‑miR‑302a‑3p‑HIF1A, MALAT1‑hsa‑miR‑519d‑3p‑MAPK1 and hsa_circRNA_100086‑hsa‑miR‑509‑3p‑MAPK1 ceRNA axes were obtained by constructing the ceRNA networks. In conclusion, these identified ceRNA interaction axes may be crucial targets for the treatment of IDD.
| S-EPMC6755180 | biostudies-literature
Project description:ObjectiveAlterations in the structure and function of intervertebral discs by multifaceted chronic processes can result in intervertebral disc degeneration (IDD). The mechanisms involved in IDD are still unknown.MethodsWe investigated the possible mechanisms underlying IDD using a bioinformatics analysis of publicly available microarray expression datasets and built a circular RNA-microRNA-mRNA (circRNA-miRNA-mRNA) network based on the results. Datasets GSE67566 and GSE116726 were downloaded from the Gene Expression Omnibus (GEO) and analyzed using the limma package in R. The CircInteractome database was used to detect miRNAs related to circRNA, and TargetScan, miRDB, and miRTarBase were used to predict target mRNAs. Key target genes were annotated using Gene Ontology terms.ResultsThe circRNA hsa-circ-0040039 was found to have the top log fold-change score. Analysis using Metascape showed that the associated genes were enriched mainly in the cell cycle. The Cytoscape plugin MCODE predicted that two members of the RAS oncogene family-RAB1A and RAB1B-and multiple coagulation factor deficiency (MCFD2) may play key roles in IDD.ConclusionOur results suggested that hsa-circ-0040039 and the related network may be potential biomarkers for IDD.
| S-EPMC8733709 | biostudies-literature
Project description:Background: Acute myocardial infarction (AMI) is one of the main fatal diseases of cardiovascular diseases. Circular RNA (circRNA) is a non-coding RNA (ncRNA), which plays a role in cardiovascular disease as a competitive endogenous RNA (ceRNA). However, their role in AMI has not been fully clarified. This study aims to explore the mechanism of circRNA-related ceRNA network in AMI, and to identify the corresponding immune infiltration characteristics. Materials and Methods: The circRNA (GSE160717), miRNA (GSE24548), and mRNA (GSE60993) microarray datasets of AMI were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed circRNAs (DEcircRNAs), miRNAs (DEmiRNAs), and mRNAs (DEmRNAs) were identified by the "limma" package. After integrating the circRNA, miRNA and mRNA interaction, we constructed a circRNA-miRNA-mRNA network. The "clusterProfiler" package and String database were used for functional enrichment analysis and protein-protein interaction (PPI) analysis, respectively. After that, we constructed a circRNA-miRNA-hub gene network and validated the circRNAs and mRNAs using an independent dataset (GSE61144) as well as qRT-PCR. Finally, we used CIBERSORTx database to analyze the immune infiltration characteristics of AMI and the correlation between hub genes and immune cells. Results: Using the "limma" package of the R, 83 DEcircRNAs, 54 DEmiRNAs, and 754 DEmRNAs were identified in the microarray datasets of AMI. Among 83 DEcircRNAs, there are 55 exonic DEcircRNAs. Then, a circRNA-miRNA-mRNA network consists of 21 DEcircRNAs, 11 DEmiRNAs, and 106 DEmRNAs were predicted by the database. After that, 10 hub genes from the PPI network were identified. Then, a new circRNA-miRNA-hub gene network consists of 14 DEcircRNAs, 7 DEmiRNAs, and 9 DEmRNAs was constructed. After that, three key circRNAs (hsa_circ_0009018, hsa_circ_0030569 and hsa_circ_0031017) and three hub genes (BCL6, PTGS2 and PTEN) were identified from the network by qRT-PCR. Finally, immune infiltration analysis showed that hub genes were significantly positively correlated with up-regulated immune cells (neutrophils, macrophages and plasma cells) in AMI. Conclusion: Our study constructed a circRNA-related ceRNA networks in AMI, consists of hsa_circ_0031017/hsa-miR-142-5p/PTEN axis, hsa_circ_0030569/hsa-miR-545/PTGS2 axis and hsa_circ_0009018/hsa-miR-139-3p/BCL6 axis. These three hub genes were significantly positively correlated with up-regulated immune cells (neutrophils, macrophages and plasma cells) in AMI. It helps improve understanding of AMI mechanism and provides future potential therapeutic targets.
| S-EPMC9002054 | biostudies-literature
Project description:BACKGROUND As all we know, gastric cancer (GC) is a highly aggressive disease. Recently, circular RNA (circRNA) was found to play a vital role in regulation of GC. Some circRNAs could regulate messenger RNA (mRNA) expression by functioning as a microRNA (miRNA) sponge. Nevertheless, the circRNA-miRNA-mRNA regulatory network involved GC rarely has been explored and researched. MATERIAL AND METHODS All the differentially expressed circRNAs, miRNAs, and mRNA were derived from Gene Expression Omnibus (GEO) microarray data (GSE78092, GSE89143, GSE93415, and GSE54129). GC level 3 miRNA-sequencing data and clinical information were downloaded from The Cancer Genome Atlas (TCGA) database. Furthermore, a circRNA-miRNA-mRNA regulatory network was constructed by Cytoscape (version 3.6.1). Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway revealed the functions and signaling pathways associated with these target genes. Hub genes of protein-protein interaction (PPI) network were identified by STRING database and cytoHubba. RESULTS The regulatory network consists of 3 circRNAs, 22 miRNAs, and 128 mRNAs. Only 3 miRNAs of the network were consistent with the expression of TCGA and were associated with some clinical features. The results of the functional analysis of 128 mRNAs showed that GO analysis and KEGG pathways of inclusion criteria were 49 and 24, respectively. PPI network and Cytoscape showed that the top 10 hub genes were MYC, CTGF, TGFBR2, TGFBR1, SERPINE1, KRAS, ZEB1, THBS1, CDK6, and TNS1; 4 of which were verified by GEPIA based on TCGA. Highly expressed SERPINE1 had a poor OS (over survival) and DFS (disease-free survival), and TGFBR1 expression increased along with the increase of clinical stages. CONCLUSIONS This study looked at a circRNA-miRNA-mRNA regulatory network associated with GC and explored the potential functions of mRNA in the network, then identified a new molecular marker for prediction, prognosis, and therapeutic targets for clinical patients.
| S-EPMC6880644 | biostudies-literature
Project description:BackgroundLow back pain (LBP) is the main factor of global disease burden. Intervertebral disc degeneration (IVDD) has long been known as the leading reason of LBP. Increasing studies have verified that circular RNAs (circRNAs)-microRNAs (miRNAs)-mRNAs network is widely involved in the pathological processes of IVDD. However, no study was made to demonstrate the circRNAs-mediated ferroptosis, oxidative stress, extracellular matrix metabolism, and immune response in IVDD.MethodsWe collected 3 normal and 3 degenerative nucleus pulposus tissues to conduct RNA-sequencing to identify the key circRNAs and miRNAs in IVDD. Bioinformatics analysis was then conducted to construct circRNAs-miRNAs-mRNAs interaction network associated with ferroptosis, oxidative stress, extracellular matrix metabolism, and immune response. We also performed animal experiments to validate the therapeutic effects of key circRNAs in IVDD.ResultsWe found that circ_0015435 was most obviously upregulated and circ_0071922 was most obviously downregulated in IVDD using RNA-sequencing. Then we observed that hsa-miR-15a-5p was the key downstream of circ_0071922, and hsa-miR-15a-5p was the top upregulated miRNA in IVDD. Bioinformatics analysis was conducted to predict that 56 immunity-related genes, 29 ferroptosis-related genes, 23 oxidative stress-related genes and 8 ECM-related genes are the targets mRNAs of hsa-miR-15a-5p. Then we constructed a ceRNA network encompassing 24 circRNAs, 6 miRNAs, and 101 mRNAs. Additionally, we demonstrated that overexpression of circ_0071922 can alleviate IVDD progression in a rat model.ConclusionsThe findings of this study suggested that circ_0071922-miR-15a-5p-mRNA signaling network might affect IVDD by modulating the nucleus pulposus cells ferroptosis, oxidative stress, ECM metabolism, and immune response, which is an effective therapeutic targets of IVDD.
| S-EPMC10782064 | biostudies-literature
Project description:BackgroundCircular RNAs (circRNAs) were known to be related to the pathogenesis of many diseases through competing endogenous RNA (ceRNA) regulatory mechanisms. However, the function of circRNA in coronary artery disease (CAD) remains unclear. In this study, we aim to construct a circRNA-related competing endogenous RNA (ceRNA) network in CAD.MethodsThe gene expression profiles of CAD were obtained from Gene Expression Omnibus datasets. Bioinformatics analysis was performed to construct a ceRNA regulatory network, from which the hub genes involved were identified through protein-protein interaction (PPI) networks leading to the identification of the circRNA-miRNA-hub gene subnetwork. In addition, function enrichment analysis was performed to detect the potential biological mechanism in which circRNA might be involved.ResultsA total of 115 DEcircRNAs (differentially expressed circRNAs), 17 DEmiRNAs (differentially expressed microRNAs), and 790 DEmRNAs (differentially expressed mRNAs) were identified between CAD and control groups from microarray datasets. Functional enrichment analysis showed that DEmRNAs were significantly involved in inflammation-related pathways and ubiquitin-protein ligase binding. After identifying 20 DEcircRNA-DEmiRNA pairs and 561 DEmiRNA-DEmRNA pairs, we obtained a circRNA-miRNA-mRNA regulatory network. PPI network analysis showed that eight hub genes were closely related to CAD, leading to the identification of a circRNA-miRNA-hub gene subnetwork consisting of nine circRNAs (hsa_circ_0020275, hsa_circ_0020387, hsa_circ_0020417, hsa_circ_0045512, hsa_circ_0047336, hsa_circ_0069094, hsa_circ_0071326, hsa_circ_0071330, and hsa_circ_0085340), four miRNAs (hsa-miR-136-5p, hsa-miR-376c-3p, hsa-miR-411-5p, and hsa-miR-654-5p), and eight mRNAs (MKRN1, UBE2H, UBE2W, UBE2D1, UBE2F, BE2J1, ZNRF1, and SIAH2). In addition, we discovered these hub genes were enriched in the ubiquitin-mediated proteolysis pathway, suggesting circRNAs may be involved in the pathogenesis of CAD through this pathway.ConclusionsThis study may deepen our understanding of the potential role of circRNA-miRNA-mRNA regulation network in CAD and suggest novel diagnostic biomarkers and therapeutic targets for CAD.
| S-EPMC8866025 | biostudies-literature
Project description:CircRNAs are a group of endogenous small noncoding RNAs that are involved in multiple diseases including cancers. At present, the functions of circRNAs in breast cancer need to be further explored. In this study, 3 pairs of breast cancer and paracancer tissues with axillary lymph node metastasis were collected for circRNA high-throughput sequencing. We have identified 17,966 distinct circRNA candidates. Significant differential expressions were found in 136 circRNAs in breast cancer tissues relative to the matched paracancer tissues. We aslo identified differentially expressed 156 miRNAs and 1105 mRNAs in breast cancer tissues and normal breast tissues from public databases. Then we constructed a regulatory ceRNA network. 12 mRNAs were associated with prognosis of breast cancer. We also constructed a circRNAs-mediated subnetwork which might be related to prognosis of breast cancer. This article provides a better understanding of circRNAs-mediated ceRNA regulatory network by which circRNAs compete for endogenous RNAs in breast cancer.
| S-EPMC9588024 | biostudies-literature
Project description:Atrial fibrillation (AF) is a chronic and progressive disease, with advancing age, the morbidity of which will increase exponentially. Circular ribonucleic acids (RNAs; circRNAs) have gained a growing attention in the development of AF in recent years. The purpose of this study is to explore the mechanism of circRNA regulation in AF, in particular, the intricate interactions among circRNA, microRNA (miRNA), and messenger RNA (mRNA). Three datasets (GSE129409, GSE68475, and GSE79768) were obtained from the Gene Expression Omnibus database to screen differentially expressed (DE) circRNAs, DE miRNAs, and DE mRNAs in AF, respectively. Based on circRNA-miRNA pairs and miRNA-mRNA pairs, a competing endogenous RNAs (ceRNAs) network was built. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis of DE mRNAs in the network were performed and protein-protein interaction (PPI) networks were established to identify hub genes. Finally, a circRNA-miRNA-hub gene subnetwork was constructed. A total of 103 DE circRNAs, 16 DE miRNAs, and 110 DE mRNAs were screened in AF. Next, ceRNAs network in AF was constructed with 3 upregulated circRNAs, 2 downregulated circRNAs, 2 upregulated miRNAs, 2 downregulated miRNAs, 17 upregulated mRNAs, and 24 downregulated mRNAs. Thirty GO terms and 6 KEGG pathways were obtained. Besides, 6 hub genes (C-X-C chemokine receptor type 4 [CXCR4], C-X-C chemokine receptor type 2 [CXCR2], C-X-C motif chemokine 11 [CXCL11], neuromedin-U, B1 bradykinin receptor, and complement C3) were screened from constructing a PPI network. Finally, a circRNA-miRNA-hub gene subnetwork with 10 regulatory axes was constructed to describe the interactions among the differential circRNAs, miRNA, and hub genes. We speculated that hsa_circRNA_0056281/hsa_circRNA_0006665 -hsa-miR-613-CXCR4/CXCR2/CXCL11 regulatory axes and hsa_circRNA_0003638-hsa-miR-1207-3p-CXCR4 regulatory axis may be associated with the pathogenesis of AF.
| S-EPMC9410607 | biostudies-literature
Project description:Background: Osteoarthritis (OA) is a degenerative joint disease that seriously affects the quality of people. Unfortunately, the pathogenesis of OA has not been fully known. Therefore, this study aimed to construct a ceRNA regulatory network related to OA to explore the pathogenesis of OA. Methods: Differentially expressed circRNAs (DEcircRNAs), microRNAs (DEmiRNAs), and mRNAs (DEmRNAs) were obtained from the Gene Expression Omnibus microarray data (GSE175959, GSE105027, and GSE169077). The miRNA response elements and target mRNAs were identified using bioinformatics approaches. Additionally, a circRNA-miRNA-mRNA network was established using Cytoscape version 3.8.0. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of mRNAs in the network were conducted to explore the possible mechanisms underlying OA development. Protein-protein interaction (PPI) analysis was performed to determine the hub genes. Based on the hub genes, a sub network was constructed using Cytoscape 3.8.0 version. Finally, connectivity map (CMap) and drug-gene interaction database (DGIdb) analyses were performed to identify the potential therapeutic targets for OA. Results: Altogether, five DEcircRNAs, 89 DEmiRNAs, and 345 DEmRNAs were identified. Moreover, a circRNA-miRNA-mRNA network was established using three circRNAs, seven miRNAs, and 37 mRNAs. GO and KEGG analyses demonstrated that the mRNAs in the network could be related to the occurrence and development of OA. PPI analysis was performed and six key genes, namely serpin family H member 1 [SERPINH1], collagen type VIII alpha 2 chain [COL8A2], collagen type XV alpha 1 chain [COL15A1], collagen type VI alpha 3 chain [COL6A3], collagen type V alpha 1 chain [COL5A1], and collagen type XI alpha 1 chain [COL11A1], were identified. Furthermore, a circRNA-miRNA-hub gene subnetwork was established in accordance with two circRNAs (hsa_circ_0075320 and hsa_circ_0051428), two miRNAs (hsa-miR-6124 and hsa-miR-1207-5p), and six hub genes (COL11A1, SERPINH1, COL6A3, COL5A1, COL8A2, and COL15A1). Finally, three chemicals (noscapine, diazepam, and TG100-115) based on CMap analysis and two drugs (collagenase Clostridium histolyticum and ocriplasmin) based on DGIdb were discovered as potential treatment options for OA. Conclusion: This study presents novel perspectives on the pathogenesis and treatment of OA based on circRNA-related competitive endogenous RNA regulatory networks.
| S-EPMC9523487 | biostudies-literature
Project description:Circular RNA expression profiling of human nucleus pulposus derived from patients with IDD in comparison with those derived from cadaveric disc as normal control. We have identified the expression profiles of miRNAs (GSE63492), lncRNAs, mRNAs (GSE56081) in IDD using 5 normal discs as control and 5 IDD discs. Accumulating evidence indicates that circRNAs are key regulators of gene expression by interacting with miRNAs. circRNA is a novel type of RNA that, unlike linear RNA, forms a covalently closed continuous loop, and is highly represented in the eukaryotic transcriptome. Two-condition experiment: control nucleus pulposus vs. degenerative nucleus pulposus. Biological replicates: 5 control, 5 degenerated, independently harvested (the same samples as GSE56081 and GSE63492). Four replicates per array.
2016-08-04 | E-GEOD-67566 | biostudies-arrayexpress