Comprehensive analysis of the competing endogenous circRNA-lncRNA-miRNA-mRNA network and identification of a novel potential biomarker for hepatocellular carcinoma.
Ontology highlight
ABSTRACT: Background
The competing endogenous RNAs (ceRNAs) hypothesis has received increasing attention as a novel explanation for tumorigenesis and cancer progression. However, there is still a lack of comprehensive analysis of the circular RNA (circRNA)-long non-coding RNA (lncRNA)-miRNA-mRNA ceRNA network in hepatocellular carcinoma (HCC).Methods
RNA sequencing data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database were employed to identify Differentially Expressed mRNAs (DEmRNAs), DElncRNAs, and DEcircRNAs between HCC and normal tissues. Candidates were identified to construct networks through a comprehensive bioinformatics strategy. A prognostic mRNA signature was established based on data from TCGA database and validated using data from the GEO database. Then, the HCC prognostic circRNA-lncRNA-miRNA-mRNA ceRNA network was established. Finally, the expression and function of an unexplored hub gene, deoxythymidylate kinase (DTYMK), was explored through data mining. The results were examined using clinical samples and in vitro experiments.Results
We constructed a prognostic signature with seven target mRNAs by univariate, lasso and multivariate Cox regression analyses, which yielded 1, 3 and 5-year AUC values of 0.797, 0.733 and 0.721, respectively, indicating its sensitivity and specificity in the prognosis of HCC. Moreover, the prognostic signature could be validated in GSE14520. The prognostic ceRNA network of 21 circRNAs, 15 lncRNAs, 5 miRNAs, and 7 mRNAs was established according to the targeting relationship between 7 hub mRNAs and other RNAs. Our experiment results indicated that the depletion of DTYMK inhibited liver cancer cell proliferation and invasion.Conclusions
The network revealed in this study may help comprehensively elucidate the ceRNA mechanisms driving HCC, and provide novel candidate biomarkers for evaluating the prognosis of HCC.
SUBMITTER: Zhang L
PROVIDER: S-EPMC8266324 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundPulpitis is a common inflammatory disease that affects dental pulp. It is important to understand the molecular signals of inflammation and repair associated with this process. Increasing evidence has revealed that long noncoding RNAs (lncRNAs), via competitively sponging microRNAs (miRNAs), can act as competing endogenous RNAs (ceRNAs) to regulate inflammation and reparative responses. The aim of this study was to elucidate the potential roles of lncRNA, miRNA and messenger RNA (mRNA) ceRNA networks in pulpitis tissues compared to normal control tissues.MethodsThe oligo and limma packages were used to identify differentially expressed lncRNAs and mRNAs (DElncRNAs and DEmRNAs, respectively) based on expression profiles in two datasets, GSE92681 and GSE77459, from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were further analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Protein-protein interaction (PPI) networks and modules were established to screen hub genes using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and the Molecular Complex Detection (MCODE) plugin for Cytoscape, respectively. Furthermore, an lncRNA-miRNA-mRNA-hub genes regulatory network was constructed to investigate mechanisms related to the progression and prognosis of pulpitis. Then, quantitative real-time polymerase chain reaction (qRT-PCR) was applied to verify critical lncRNAs that may significantly affect the pathogenesis in inflamed and normal human dental pulp.ResultsA total of 644 upregulated and 264 downregulated differentially expressed genes (DEGs) in pulpitis samples were identified from the GSE77459 dataset, while 8 up- and 19 downregulated probes associated with lncRNA were identified from the GSE92681 dataset. Protein-protein interaction (PPI) based on STRING analysis revealed a network of DEGs containing 4,929 edges and 623 nodes. Upon combined analysis of the constructed PPI network and the MCODE results, 10 hub genes, including IL6, IL8, PTPRC, IL1B, TLR2, ITGAM, CCL2, PIK3CG, ICAM1, and PIK3CD, were detected in the network. Next, a ceRNA regulatory relationship consisting of one lncRNA (PVT1), one miRNA (hsa-miR-455-5p) and two mRNAs (SOCS3 and PLXNC1) was established. Then, we constructed the network in which the regulatory relationship between ceRNA and hub genes was summarized. Finally, our qRT-PCR results confirmed significantly higher levels of PVT1 transcript in inflamed pulp than in normal pulp tissues (p = 0.03).ConclusionOur study identified a novel lncRNA-mediated ceRNA regulatory mechanisms in the pathogenesis of pulpitis.
| S-EPMC6609876 | biostudies-literature
Project description:BackgroundIncreasing evidence has underscored the role of long non-coding RNAs (lncRNAs) acting as competing endogenous RNAs (ceRNAs) in the development and progression of tumors. Nevertheless, lncRNA biomarkers in lncRNA-related ceRNA network that can predict the prognosis of breast cancer (BC) are still lacking. The aim of our study was to identify potential lncRNA signatures capable of predicting overall survival (OS) of BC patients.MethodsThe RNA sequencing data and clinical characteristics of BC patients were obtained from the Cancer Genome Atlas database, and differentially expressed lncRNA (DElncRNAs), DEmRNAs, and DEmiRNAs were then identified between BC and normal breast tissue samples. Subsequently, the lncRNA-miRNA-mRNA ceRNA network of BC was established, and the gene oncology enrichment analyses for the DEmRNAs interacting with lncRNAs in the ceRNA network was implemented. Using univariate and multivariate Cox regression analyses, a four-lncRNA signature was developed and used for predicting the survival in BC patients. We applied receiver operating characteristic analysis to assess the performance of our model.ResultsA total of 1061 DElncRNAs, 2150 DEmRNAs, and 82 DEmiRNAs were identified between BC and normal breast tissue samples. A lncRNA-miRNA-mRNA ceRNA network of BC was established, which comprised of 8 DEmiRNAs, 48 DElncRNAs, and 10 DEmRNAs. Further gene oncology enrichment analyses revealed that the DEmRNAs interacting with lncRNAs in the ceRNA network participated in cell leading edge, protease binding, alpha-catenin binding, gamma-catenin binding, and adenylate cyclase binding. A univariate regression analysis of the DElncRNAs revealed 7 lncRNAs (ADAMTS9-AS1, AC061992.1, LINC00536, HOTAIR, AL391421.1, TLR8-AS1 and LINC00491) that were associated with OS of BC patients. A multivariate Cox regression analysis demonstrated that 4 of those lncRNAs (ADAMTS9-AS1, LINC00536, AL391421.1 and LINC00491) had significant prognostic value, and their cumulative risk score indicated that this 4-lncRNA signature independently predicted OS in BC patients. Furthermore, the area under the curve of the 4-lncRNA signature associated with 3-year survival was 0.696.ConclusionsThe current study provides novel insights into the lncRNA-related ceRNA network in BC and the 4 lncRNA biomarkers may be independent prognostic signatures in predicting the survival of BC patients.
| S-EPMC6161429 | biostudies-literature
Project description:Background: Hepatocellular carcinoma (HCC) accounts for the majority of liver cancer, with the incidence and mortality rates increasing every year. Despite the improvement of clinical management, substantial challenges remain due to its high recurrence rates and short survival period. This study aimed to identify potential diagnostic and prognostic biomarkers in HCC through bioinformatic analysis. Methods: Datasets from GEO and TCGA databases were used for the bioinformatic analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were carried out by WebGestalt website and clusterProfiler package of R. The STRING database and Cytoscape software were used to establish the protein-protein interaction (PPI) network. The GEPIA website was used to perform expression analyses of the genes. The miRDB, miRWalk, and TargetScan were employed to predict miRNAs and the expression levels of the predicted miRNAs were explored via OncomiR database. LncRNAs were predicted in the StarBase and LncBase while circRNA prediction was performed by the circBank. ROC curve analysis and Kaplan-Meier (KM) survival analysis were performed to evaluate the diagnostic and prognostic value of the gene expression, respectively. Results: A total of 327 upregulated and 422 downregulated overlapping DEGs were identified between HCC tissues and noncancerous liver tissues. The PPI network was constructed with 89 nodes and 178 edges and eight hub genes were selected to predict upstream miRNAs and ceRNAs. A lncRNA/circRNA-miRNA-mRNA network was successfully constructed based on the ceRNA hypothesis, including five lncRNAs (DLGAP1-AS1, GAS5, LINC00665, TYMSOS, and ZFAS1), six circRNAs (hsa_circ_0003209, hsa_circ_0008128, hsa_circ_0020396, hsa_circ_0030051, hsa_circ_0034049, and hsa_circ_0082333), eight miRNAs (hsa-miR-150-5p, hsa-miR-19b-3p, hsa-miR-23b-3p, hsa-miR-26a-5p, hsa-miR-651-5p, hsa-miR-10a-5p, hsa-miR-214-5p and hsa-miR-486-5p), and five mRNAs (CDC6, GINS1, MCM4, MCM6, and MCM7). The ceRNA network can promote HCC progression via cell cycle, DNA replication, and other pathways. Clinical diagnostic and survival analyses demonstrated that the ZFAS1/hsa-miR-150-5p/GINS1 ceRNA regulatory axis had a high diagnostic and prognostic value. Conclusion: These results revealed that cell cycle and DNA replication pathway could be potential pathways to participate in HCC development. The ceRNA network is expected to provide potential biomarkers and therapeutic targets for HCC management, especially the ZFAS1/hsa-miR-150-5p/GINS1 regulatory axis.
| S-EPMC8977626 | biostudies-literature
Project description:Growing evidence has illustrated critical roles of competing endogenous RNA (ceRNA) regulatory network in human cancers including hepatocellular carcinoma. In this study, we aimed to find promising diagnostic and prognostic biomarkers for patients with hepatocellular carcinoma. Three novel unfavorable prognosis-associated genes (CELSR3, GPSM2, and CHEK1) was first identified. We also demonstrated that these genes were significantly upregulated in hepatocellular carcinoma cell lines and tissues. Next, 154 potential miRNAs of CELSR3, GPSM2, and CHEK1 were predicted. CHEK1-hsa-mir-195-5p/hsa-mir-497-5p and GPSM2-hsa-mir-122-5p axes were defined as two key pathways in carcinogenesis of hepatocellular carcinoma by combination of in silico analysis and experimental validation. Subsequently, lncRNAs binding to hsa-mir-195-5p, hsa-mir-497-5p, and hsa-mir-122-5p were predicted via starBase and miRNet databases. After performing expression analysis and survival analysis for these predicted lncRNAs, we showed that nine lncRNAs (SNHG1, SNHG12, LINC00511, HCG18, FGD5-AS1, CERS6-AS1, NUTM2A-AS1, SNHG16, and ASB16-AS1) were markedly increased in hepatocellular carcinoma and their upregulation indicated poor prognosis. Moreover, a similar mRNA-miRNA-lncRNA analysis for six "known" genes (CLEC3B, DNASE1L3, PTTG1, KIF2C, XPO5, and UBE2S) was performed. Subsequently, a comprehensive mRNA-miRNA-lncRNA triple ceRNA network linked to prognosis of patients with hepatocellular carcinoma was established. Moreover, all RNAs in this network exhibited significantly diagnostic values for patients with hepatocellular carcinoma. In summary, the current study constructed a mRNA-miRNA-lncRNA ceRNA network associated with diagnosis and prognosis of hepatocellular carcinoma.
| S-EPMC7092636 | biostudies-literature
Project description:BackgroundCholangiocarcinoma (CCA) is the second most common malignant primary liver tumor and has shown an alarming increase in incidence over the last two decades. However, the mechanisms behind tumorigenesis and progression remain insufficient. The present study aimed to uncover the underlying regulatory mechanism on CCA and find novel biomarkers for the disease prognosis.MethodThe RNA-sequencing (RNA-seq) datasets of lncRNAs, miRNAs, and mRNAs in CCA as well as relevant clinical information were obtained from the Cancer Genome Atlas (TCGA) database. After pretreatment, differentially expressed RNAs (DERNAs) were identified and further interrogated for their correlations with clinical information. Prognostic RNAs were selected using univariate Cox regression. Then, a ceRNA network was constructed based on these RNAs.ResultsWe identified a total of five prognostic DEmiRNAs, 63 DElncRNAs, and 90 DEmRNAs between CCA and matched normal tissues. Integrating the relationship between the different types of RNAs, an lncRNA-miRNA-mRNA network was established and included 28 molecules and 47 interactions. Screened prognostic RNAs involved in the ceRNA network included 3 miRNAs (hsa-mir-1295b, hsa-mir-33b, and hsa-mir-6715a), 7 lncRNAs (ENSG00000271133, ENSG00000233834, ENSG00000276791, ENSG00000241155, COL18A1-AS1, ENSG00000274737, and ENSG00000235052), and 18 mRNAs (ANO9, FUT4, MLLT3, ABCA3, FSCN2, GRID2IP, NCK2, MACC1, SLC35E4, ST14, SH2D3A, MOB3B, ACTL10, RAB36, ATP1B3, MST1R, SEMA6A, and SEL1L3).ConclusionsOur study identified novel prognostic makers and predicted a previously unknown ceRNA regulatory network in CCA and may provide novel insight into a further understanding of lncRNA-mediated ceRNA regulatory mechanisms in CCA.
| S-EPMC6451803 | biostudies-literature
Project description:The competing endogenous RNA (ceRNA) axis has been shown to play a critical role in the pathogenesis of various viral infections. Generally, the ceRNA network involves long non-coding RNAs (lncRNAs) that act as sponges for miRNA to regulate mRNA expression. However, no information is available regarding the involvement of ceRNA networks in Enterovirus type 71 (EV71) infections. In the present study, data obtained from Gene Expression Omnibus (GEO) database was analyzed using various bioinformatics tools. EV71 infection in rhabdomyosarcoma (RD) cells was associated with differential expression of six lncRNAs, 28 miRNAs, and 349 mRNAs. Gene function enrichment analysis suggested induction of cytoplasmic vesicle process upon EV71 infection. The ceRNA networks were constructed, in which 20 hub genes were predicted by protein-protein interaction. To confirm the MALAT1/miR-194-5p/DUSP1 ceRNA regulatory axis in EV71 infection, real-time quantitative polymerase chain reaction (qRT-PCR) and luciferase reporter assay were performed. The results of the study also revealed the involvement of the MALAT1/miR-194-5p axis in apoptosis induced by EV71 infection, while no association with autophagy was observed. Thus, the present study provided novel insights into the pathogenic mechanism of EV71 infection.
| S-EPMC8340214 | biostudies-literature
Project description:AimsTo construct the lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network based on our microarray chip data for providing new insights into the pathogenesis of autoimmune hepatitis.MethodsThe ceRNA pairs were obtained by calculating the co-expression relationships among the differentially expressed lncRNAs (DELs), differentially expressed microRNAs (DEMis), and differentially expressed mRNAs (DEMs) with Pearson correlation analysis and hypergeometric distribution. The data of the differentially expressed genes were obtained from our previous studies in the concanavalin A-induced AIH mouse model. The biological functions of the ceRNA network were revealed by carrying out the GO and KEGG enrichment analysis. The expression of some differentially expressed genes constructed in the ceRNA pair was validated, and the correlation to liver injury was analyzed.ResultsThe mRNAs constructed in the ceRNA network were most significantly annotated in the GO terms of "inflammatory response" and enriched in "Cytokine-cytokine receptor interaction" and "MAPK signaling pathway". The differences in the expression of Gm38975, mmu-miR-125a-3p, and Map3k13 between the model group and control group were significant, and the expression of these genes at a transcriptional level was positively or negatively correlated to the activity of ALT and AST as well as the amount of MDA and NO.ConclusionOur work is the first in its kind to predict and illustrate the comprehensive lncRNA-miRNA-mRNA ceRNA network associated with the etiopathogenesis of AIH. This study indicates to lay the foundation for revealing the potential roles of ceRNAs in the occurrence of AIH and provide novel treatment targets for this disease.
| S-EPMC9790135 | biostudies-literature
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:BackgroundRecently, increasing evidence has uncovered the roles of mRNA-miRNA-lncRNA network in multiple human cancers. However, a systematic mRNA-miRNA-lncRNA network linked to pancreatic cancer prognosis is still absent.MethodsDifferentially expressed genes (DEGs) were first identified by mining GSE16515 and GSE15471 datasets. DAVID database was utilized to conduct functional enrichment analysis. Protein-protein interaction (PPI) network was built using STRING database, and hub genes were identified by Cytoscape plug-in CytoHubba. Upstream miRNAs and lncRNAs of mRNAs were predicted by miRTarBase and miRNet, respectively. Expression, survival and correlation analysis for genes, miRNAs and lncRNAs were performed via GEPIA, Kaplan-Meier plotter and starBase.Results734 and 180 upregulated and downregulated significant DEGs were identified, respectively. Functional enrichment analysis revealed that they were significantly enriched in focal adhesion, pathways in cancer and metabolic pathways. According to node degree, hub genes in the PPI networks were screened, such as TGFB1 and ALB. Among the top 20 hub genes, 7 upregulated genes and 2 downregulated hub genes had significant prognostic values in pancreatic cancer. 33 miRNAs were predicted to target the 9 key genes. But only high expression of 8 miRNAs indicated favorable prognosis in pancreatic cancer. Then, 90 lncRNAs were predicted to potentially bind to the 8 miRNAs. SCAMP1, HCP5, MAL2 and LINC00511 were finally identified as key lncRNAs. By combination of results from expression, survival and correlation analysis demonstrated that MMP9/ITGB1-miR-29b-3p-HCP5 competing endogenous RNA (ceRNA) sub-network was linked to prognosis of pancreatic cancer.ConclusionsIn a word, we established a novel mRNA-miRNA-lncRNA sub-network, among which each RNA may be utilized as a prognostic biomarker of pancreatic cancer.
| S-EPMC6535056 | biostudies-literature
Project description:BackgroundAccumulated studies have revealed that long non-coding RNAs (lncRNAs) play critical roles in human diseases by acting as competing endogenous RNAs (ceRNAs). However, functional roles and regulatory mechanisms of lncRNA-mediated ceRNA in atrial fibrillation (AF) remain unknown. In the present study, we aimed to construct the lncRNA-miRNA-mRNA network based on ceRNA theory in AF by using bioinformatic analyses of public datasets.MethodsMicroarray data sets of GSE115574 and GSE79768 from the Gene Expression Omnibus database were downloaded. Twenty-one AF right atrial appendage (RAA) samples and 22 sinus rhythm (SR) subjects RAA samples were selected for subsequent analyses. After merging all microarray data and adjusting for batch effect, differentially expressed genes were identified. Gene Ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were carried out. A ceRNA network was constructed.ResultA total of 8 lncRNAs and 43 mRNAs were significantly differentially expressed with fold change >1.5 (p < 0.05) in RAA samples of AF patients when compared with SR. GO and KEGG pathway analysis showed that cardiac muscle contraction pathway were involved in AF development. The ceRNA was predicted by co-expressing LOC101928304/ LRRC2 from the constructional network analysis, which was competitively combined with miR-490-3p. The expression of LOC101928304 and LRRC were up-regulated in myocardial tissue of patients with AF, while miR-490-3p was down-regulated.ConclusionWe constructed the LOC101928304/miR-490-3p/LRRC2 network based on ceRNA theory in AF in the bioinformatic analyses of public datasets. The ceRNA network found from this study may help improve our understanding of lncRNA-mediated ceRNA regulatory mechanisms in the pathogenesis of AF.
| S-EPMC8818759 | biostudies-literature