Dataset Information

Construction of circRNA-miRNA-mRNA network and identification of novel potential biomarkers for non-small cell lung cancer

ABSTRACT:

Background

The underlying circular RNAs (circRNAs)-related competitive endogenous RNA (ceRNA) mechanisms of pathogenesis and prognosis in non-small cell lung cancer (NSCLC) remain unclear.

Methods

Differentially expressed circRNAs (DECs) in two Gene Expression Omnibus datasets (GSE101684 and GSE112214) were identified by utilizing R package (Limma). Circinteractome and StarBase databases were used to predict circRNA associated-miRNAs and mRNAs, respectively. Then, protein–protein interaction (PPI) network of hub genes and ceRNA network were constructed by STRING and Cytoscape. Also, analyses of functional enrichment, genomic mutation and diagnostic ROC were performed. TIMER database was used to analyze the association between immune infiltration and target genes. Kaplan–Meier analysis, cox regression and the nomogram prediction model were used to evaluate the prognostic value of target genes. Finally, the expression of potential circRNAs and target genes was validated in cell lines and tissues by quantitative real-time PCR (qRT-PCR) and Human Protein Atlas (HPA) database.

Results

In this study, 15 DECs were identified between NSCLC tissues and adjacent-normal tissues in two GEO datasets. Following the qRT-PCR corroboration, 7 DECs (hsa_circ_0002017, hsa_circ_0069244, hsa_circ_026337, hsa_circ_0002346, hsa_circ_0007386, hsa_circ_0008234, hsa_circ_0006857) were dramatically downregulated in A549 and SK-MES-1 compared with HFL-1 cells. Then, 12 circRNA-sponged miRNAs were screened by Circinteractome and StarBase, especially, hsa-miR-767-3p and hsa-miR-767-5p were significantly up-regulated and relevant to the prognosis. Utilizing the miRDB and Cytoscape, 12 miRNA-target genes were found. Functional enrichment, genomic mutation and diagnostic analyses were also performed. Among them, FNBP1, AKT3, HERC1, COL4A1, TOLLIP, ARRB1, FZD4 and PIK3R1 were related to the immune infiltration via TIMER database. The expression of ARRB1, FNBP1, FZD4, and HERC1 was correlated with poor overall survival (OS) in NSCLC patients by cox regression and nomogram. Furthermore, the hub-mRNAs were validated in cell lines and tissues.

Conclusion

We constructed the circRNA-miRNA-mRNA network that might provide novel insights into the pathogenesis of NSCLC and reveal promising immune infiltration and prognostic biomarkers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02278-z.

SUBMITTER: Yang J

PROVIDER: S-EPMC8605517 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundOsteosarcoma is a common malignant primary bone tumor in adolescents and children. Numerous studies have shown that circRNAs were involved in the proliferation and invasion of various tumors. However, the role of circRNAs in osteosarcoma remains unclear. Here, we aimed to explore the regulatory network among circRNA-miRNA-mRNA in osteosarcoma.MethodsThe circRNA (GSE140256), microRNA (GSE28423), and mRNA (GSE99671) expression profiles of osteosarcoma were collected from the Gene Expression Omnibus (GEO) database. Differentially expressed circRNAs, miRNAs and mRNAs were identified. CircRNA-miRNA interactions and miRNA-mRNA interactions were determined by Circular RNA Interactome (CircInteractome) database and microRNA Data Integration Portal (mirDIP) database, respectively. Then, we constructed a regulatory network. Function enrichment analysis of miRNA and mRNA was performed by DIANA-miRPath v3.0 and Metascape database, respectively. mRNAs with significant prognostic value were identified based on expression profiles from The Cancer Genome Atlas (TCGA) database, and we constructed a subnetwork for them. To make the most of the network, we used the CLUE database to predict potential drugs for the treatment of osteosarcoma based on mRNA expression in the network. And we used the STITCH database to analyze and validate the interactions among these drugs and mRNAs, and to further screen for potential drugs.ResultsA total of 9 circRNAs, 19 miRNAs, 67 mRNAs, 54 pairs of circRNA-miRNA interactions and 110 pairs of miRNA-mRNA interactions were identified. A circRNA-miRNA-mRNA network was constructed. Function enrichment analysis indicated that these miRNAs and mRNAs in the network were involved in the process of tumorigenesis and immune response. Among these mRNAs, STC2 and RASGRP2 with significantly prognostic value were identified, and we constructed a subnetwork for them. Based on mRNA expression in the network, three potential drugs, quinacridine, thalidomide and zonisamide, were screened for the treatment of osteosarcoma. Among them, quinacridine and thalidomide have been proved to have anti-tumor effects in previous studies, while zonisamide has not been reported. And a corresponding drug-protein interaction network was constructed.ConclusionOverall, we constructed a circRNA-miRNA-mRNA regulatory network to investigate the possible mechanism in osteosarcoma, and predicted that quinacridine, thalidomide and zonisamide could be potential drugs for the treatment of osteosarcoma.

Project description:BackgroundCircular RNA (circRNA) plays an important role in the regulation of gene expression and the occurrence of human diseases. However, studies on the role of circRNA in acute myocardial infarction (AMI) are limited. This study was performed to explore novel circRNA-related regulatory networks in AMI, aiming to better understand the molecular mechanism of circRNAs involvement in AMI and provide basis for further scientific research and clinical decision-making.MethodsThe AMI-related microarray datasets GSE160717 (circRNA), GSE31568 (miRNA), GSE61741 (miRNA), and GSE24519 (mRNA) were obtained from the Gene Expression Omnibus (GEO) database. After differential expression analysis, the regulatory relationships between these DERNAs were identified by online databases circBank, circInteractome, miRDB, miRWalk, Targetscan, and then two circRNA-miRNA-mRNA regulatory networks were constructed. Differentially expressed genes (DEGs) in this network were selected followed by enrichment analysis and protein-protein interaction (PPI) analysis. Hub genes were identified using Cytohubba plug-in of Cytoscape software. Hub genes and hub gene-related miRNAs were used for receiver operating characteristic curve (ROC) analysis to identify potential biomarkers. The relative expression levels of these biomarkers were further assessed by GSE31568 (miRNA) and GSE66360 (mRNA). Finally, on the basis of the above analysis, myocardial hypoxia model was constructed to verify the expression of Hub genes and related circRNAs.ResultsA total of 83 DEcircRNAs, 109 CoDEmiRNAs and 1204 DEGs were significantly differentially expressed in these datasets. The up-regulated circRNAs and down-regulated circRNAs were used to construct a circRNA-miRNA-mRNA regulatory network respectively. These circRNA-related DEGs were mainly enriched in the terms of "FOXO signaling pathway," "T cell receptor signaling pathway," "MAPK signaling pathway," "Insulin resistance," "cAMP signaling pathway," and "mTOR signaling pathway." The top 10 hub genes ATP2B2, KCNA1, GRIN2A, SCN2B, GPM6A, CACNA1E, HDAC2, SRSF1, ANK2, and HNRNPA2B1 were identified from the PPI network. Hub genes GPM6A, SRSF1, ANK2 and hub gene-related circRNAs hsa_circ_0023461, hsa_circ_0004561, hsa_circ_0001147, hsa_circ_0004771, hsa_circ_0061276, and hsa_circ_0045519 were identified as potential biomarkers in AMI.ConclusionIn this study, the potential circRNAs associated with AMI were identified and two circRNA-miRNA-mRNA regulatory networks were constructed. This study explored the mechanism of circRNA involvement in AMI and provided new clues for the selection of new diagnostic markers and therapeutic targets for AMI.

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.

Project description:BackgroundAbdominal aortic aneurysm (AAA) is a progressive cardiovascular disease, which is a permanent and localized dilatation of the abdominal aorta with potentially fatal consequence of aortic rupture. Dysregulation of circRNAs is correlated with the development of various pathological events in cardiovascular diseases. However, the function of circRNAs in abdominal aortic aneurysm (AAA) is unknown and remains to be explored. This study is aimed at determining the regulatory mechanisms of circRNAs in AAAs. This study was aimed at exploring the underlying molecular mechanisms of abdominal aortic aneurysms based on the competing endogenous RNA (ceRNA) regulatory hypothesis of circRNA, miRNA, and mRNA.MethodsThe expression profiles of circRNAs (GSE144431), miRNAs (GSE62179), and mRNAs (GSE7084, GSE57691, and GSE47472) in human tissue sample from the aneurysm group and normal group were obtained from the Gene Expression Omnibus database, respectively. The circRNA-miRNA-mRNA network was constructed by using Cytoscape 3.7.2 software; then, the protein-protein interaction (PPI) network was constructed by using the STRING database, and the hub genes were identified by using the cytoHubba plug-in. The circRNA-miRNA-hub gene regulatory subnetwork was formed to understand the regulatory axis of hub genes in AAAs.ResultsThe present study identified 40 differentially expressed circRNAs (DECs) in the GSE144431, 90 differentially expressed miRNAs (DEmiRs) in the GSE62179, and 168 differentially expressed mRNAs (DEGs) with the same direction regulation (130 downregulated and 38 upregulated) in the GSE7084, GSE57691, and GSE47472 datasets identified regarding AAAs. The miRNA response elements (MREs) of three DECs were then predicted. Four overlapping miRNAs were obtained by intersecting the predicted miRNA and DEmiRs. Then, 17 overlapping mRNAs were obtained by intersecting the predicted target mRNAs of 4 miRNAs with 168 DEGs. Furthermore, the circRNA-miRNA-mRNA network was constructed through 3 circRNAs, 4 miRNAs, and 17 mRNAs, and three hub genes (SOD2, CCR7, and PGRMC1) were identified. Simultaneously, functional enrichment and pathway analysis were performed within genes in the circRNA-miRNA-mRNA network. Three of them (SOD2, CCR7, and PGRMC1) were suggested to be crucial based on functional enrichment, protein-protein interaction, and ceRNA network analysis. Furthermore, the expression of SOD2 and CCR7 may be regulated by hsa_circ_0011449/hsa_circ_0081968/hsa-let-7f-5p; the expression of PGRMC1 may be regulated by hsa_circ_0011449/hsa_circ_0081968-hsa-let-7f-5p/hsa-let-7e-5p.ConclusionIn conclusion, the ceRNA interaction axis we identified may be an important target for the treatment of abdominal aortic aneurysms. This study provided further understanding of the potential pathogenesis from the perspective of the circRNA-related competitive endogenous RNA network in AAAs.