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

Project description:BackgroundPreeclampsia (PE), one of hypertension-related disorders of pregnancy, is a common cause of maternal death worldwide. This study aimed to identify a circRNA-miRNA-mRNA-associated ceRNA network and related pathways in PE.Material and methodsWe downloaded 3 microarray datasets from the Gene Expression Omnibus database, obtained 163 differentially expressed circRNAs (dif-circRNAs) (61 upregulated and 102 downregulated), 39 differentially expressed microRNAs (dif-miRNAs) (22 upregulated and 17 downregulated), and 271 differentially expressed mRNAs (dif-mRNAs) (168 upregulated and 103 downregulated) from placenta tissues of PE. Functional enrichment analysis and protein-protein interaction (PPI) network with module analysis of dif-mRNAs were performed. The regulatory relationship between dif-miRNAs and dif-mRNAs/circRNAs was predicted via related databases. A circRNA-miRNA-mRNA regulatory network was constructed.ResultsA total of 53 pairs were obtained, including 13 circRNAs (10 upregulated and 3 downregulated), 9 miRNAs (3 upregulated and 6 downregulated) and 31 mRNAs (22 upregulated and 9 downregulated). GNB5 and IL2RB were obtained. KEGG enrichment analysis showed that both of them were closely related with the PI3K-Akt signalling pathway. Therefore, ceRNAs might affect the PI3K-Akt signalling pathway via the upregulation of GNB5 by binding to miR-1248 in PE. Meanwhile, hsa_circ_0052661 might upregulate IL2RB by binding miR-4303 to play a role in PE in the same way.ConclusionGNB5 and IL2RB might be key genes involved in the PI3K-Akt signalling pathway in PE, and hsa_circ_0087208, hsa_circ_0035443, hsa_circ_0067557 and hsa_circ_0052661 might regulate these key genes in PE by binding miR-1248 or miR-4303.Key messagesThere is still a lack of predictive and diagnostic factors for preeclampsia, which is a common cause of maternal death worldwide.This study identified a novel circRNA-associated ceRNA network and related pathways in preeclampsia.GNB5 and IL2RB might be key genes in their related circRNA-associated ceRNA network, and probably take an important role in preeclampsia via PI3K-Akt signalling pathway, which made them to be potential markers of preeclampsia.