scholarly journals Construction and Analysis of Survival-Associated Competing Endogenous RNA Network in Lung Adenocarcinoma

2021 ◽  
Vol 2021 ◽  
pp. 1-17
Author(s):  
Lixian Chen ◽  
Zhonglu Ren ◽  
Yongming Cai

Increasing evidence has shown that noncoding RNAs play significant roles in the initiation, progression, and metastasis of tumours via participating in competing endogenous RNA (ceRNA) networks. However, the survival-associated ceRNA in lung adenocarcinoma (LUAD) remains poorly understood. In this study, we aimed to investigate the regulatory mechanisms underlying ceRNA in LUAD to identify novel prognostic factors. mRNA, lncRNA, and miRNA sequencing data obtained from the GDC data portal were utilized to identify differentially expressed (DE) RNAs. Survival-related RNAs were recognized using univariate Kaplan-Meier survival analysis. We performed functional enrichment analysis of survival-related mRNAs using the clusterProfiler package of R and STRING. lncRNA-miRNA and miRNA-mRNA interactions were predicted based on miRcode, Starbase, and miRanda. Subsequently, the survival-associated ceRNA network was constructed for LUAD. Multivariate Cox regression analysis was used to identify prognostic factors. Finally, we acquired 15 DE miRNAs, 49 DE lncRNAs, and 843 DE mRNAs associated with significant overall survival. Functional enrichment analysis indicated that survival-related DE mRNAs were enriched in cell cycle. The survival-associated lncRNA-miRNA-mRNA ceRNA network was constructed using five miRNAs, 49 mRNAs, and 21 lncRNAs. Furthermore, seven hub RNAs (LINC01936, miR-20a-5p, miR-31-5p, TNS1, TGFBR2, SMAD7, and NEDD4L) were identified based on the ceRNA network. LINC01936 and miR-31-5p were found to be significant using the multifactorial Cox regression model. In conclusion, we successfully constructed a survival-related lncRNA-miRNA-mRNA ceRNA regulatory network in LUAD and identified seven hub RNAs, which provide novel insights into the regulatory molecular mechanisms associated with survival of LUAD, and identified two independent prognostic predictors for LUAD.

2020 ◽  
Author(s):  
Yang Wang ◽  
Chengping Hu

Abstract Background: Long non-coding RNAs (lncRNAs) have been reported to play essential roles in tumorigenesis and cancers prognosis, and they can be a potential cancer prognostic markers. However, in lung adenocarcinoma(LUAD), how lncRNA signatures predict the survival of patients is poorly understood. Our study aims to explore lncRNA signatures and prognostic function in LUAD.Methods: The expression and prognosis data of lncRNAs in LUAD patients was collected from the Cancer Genome Atlas (TCGA) data. All analyses were performed using the R package (version 3.6.2). Metascape, STRING and Cytoscape were used for enrichment analysis and function prediction of the lncRNA co-expressed protein-coding genes.Results: We have collected lncRNA expression data in 466 LUAD tumors, and a six-lncRNA signature(RP11-79H23.3, RP11-309M7.1, CTD-2357A8.3, RP11-108P20.4, U47924.29, LHFPL3-AS2) has been shown to be significantly related to LUAD patients’ overall survival. According to the lncRNA signatures, the high-risk and low-risk groups were divided in LUAD patients with different survival rates. Further multivariable cox regression analysis showed that the prognostic value of this signature was independent of clinical factors. The potential functional roles and hub co-expressed protein-coding genes in the six prognostic lncRNAs are shown in the functional enrichment analysis.Conclusions: These results showed that these six lncRNAs could be independent predicted prognostic biomarkers in LUAD patients.


2020 ◽  
Vol 29 (3) ◽  
pp. 399-416
Author(s):  
Dan Yang ◽  
Yang He ◽  
Bo Wu ◽  
Ruxi Liu ◽  
Nan Wang ◽  
...  

BACKGROUND: Lung adenocarcinoma (LUAD) is the most common histological subtype of lung cancer worldwide. Until now, the molecular mechanisms underlying LUAD progression have not been fully explained. This study aimed to construct a competing endogenous RNA (ceRNA) network to predict the progression in LUAD. METHODS: Differentially expressed lncRNAs (DELs), miRNAs (DEMs), and mRNAs (DEGs) were identified from The Cancer Genome Atlas (TCGA) database with a |log2FC|> 1.0 and a false discovery rate (FDR) < 0.05. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), protein-protein interaction (PPI) network, and survival analyses were performed to analyse these DEGs involved in the ceRNA network. Subsequently, the drug-gene interaction database (DGIdb) was utilized to select candidate LUAD drugs interacting with significant DEGs. Then, lasso-penalized Cox regression and multivariate Cox regression models were used to construct the risk score system. Finally, based on the correlations between DELs and DEGs involved in the risk score system, the final ceRNA network was identified. Meanwhile, the GEPIA2 database and immunohistochemical (IHC) results were utilized to validate the expression levels of selected DEGs. RESULTS: A total of 340 DELs, 29 DEMs, and 218 DEGs were selected to construct the initial ceRNA network. Functional enrichment analyses indicated that 218 DEGs were associated with the KEGG pathway terms “microRNAs in cancer”, “pathways in cancer”, “cell cycle”, “HTLV-1 infection”, and the “PI3K-Akt signalling pathway”. K-M survival analysis of all differentially expressed genes involved in the ceRNA network identified 24 DELs, 4 DEMs, and 29 DEGs, all of which were significantly correlated with LUAD progression (P< 0.05). Furthermore, 15 LUAD drugs interacting with 29 significant DEGs were selected. After lasso-penalized Cox regression and multivariate Cox regression modelling, PRKCE, DLC1, LATS2, and DPY19L1 were incorporated into the risk score system, and the results suggested that LUAD patients who had the high-risk score always suffered from a poorer overall survival. Additionally, the correlation coefficients between these 4 DEGs and their corresponding DELs involved in the ceRNA network suggested that there were 2 significant DEL-DEG pairs, NAV2-AS2 – PRKCE (r= 0.430, P< 0.001) and NAV2-AS2 – LATS2 (r= 0.338, P< 0.001). And NAV2-AS2 – mir-31 – PRKCE and NAV2-SA2 – mir-31 – LATS2 were finally identified as ceRNA network involved in the progression of LUAD. CONCLUSIONS: The lncRNA-miRNA-mRNA ceRNA network plays an essential role in predicting the progression of LUAD. These results may improve our understanding and provide novel mechanistic insights to explore prognosis and therapeutic drugs for LUAD patients.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Wenjie Chen ◽  
Wen Li ◽  
Zhenkun Liu ◽  
Guangzhi Ma ◽  
Yunfu Deng ◽  
...  

AbstractTo identify the prognostic biomarker of the competitive endogenous RNA (ceRNA) and explore the tumor infiltrating immune cells (TIICs) which might be the potential prognostic factors in lung adenocarcinoma. In addition, we also try to explain the crosstalk between the ceRNA and TIICs to explore the molecular mechanisms involved in lung adenocarcinoma. The transcriptome data of lung adenocarcinoma were obtained from The Cancer Genome Atlas (TCGA) database, and the hypergeometric correlation of the differently expressed miRNA-lncRNA and miRNA-mRNA were analyzed based on the starBase. In addition, the Kaplan–Meier survival and Cox regression model analysis were used to identify the prognostic ceRNA network and TIICs. Correlation analysis was performed to analysis the correlation between the ceRNA network and TIICs. In the differently expressed RNAs between tumor and normal tissue, a total of 190 miRNAs, 224 lncRNAs and 3024 mRNAs were detected, and the constructed ceRNA network contained 5 lncRNAs, 92 mRNAs and 10 miRNAs. Then, six prognostic RNAs (FKBP3, GPI, LOXL2, IL22RA1, GPR37, and has-miR-148a-3p) were viewed as the key members for constructing the prognostic prediction model in the ceRNA network, and three kinds of TIICs (Monocytes, Macrophages M1, activated mast cells) were identified to be significantly related with the prognosis in lung adenocarcinoma. Correlation analysis suggested that the FKBP3 was associated with Monocytes and Macrophages M1, and the GPI was obviously related with Monocytes and Macrophages M1. Besides, the LOXL2 was associated with Monocytes and Activated mast cells, and the IL22RA1 was significantly associated with Monocytes and Macrophages M1, while the GPR37 and Macrophages M1 was closely related. The constructed ceRNA network and identified Monocytes, Macrophages M1 and activated Mast cells are all prognostic factors for lung adenocarcinoma. Moreover, the crosstalk between the ceRNA network and TIICs might be a potential molecular mechanism involved.


2020 ◽  
Author(s):  
Gaochen Lan ◽  
Xiaoling Yu ◽  
Yanna Zhao ◽  
Jinjian Lan ◽  
Wan Li ◽  
...  

Abstract Background: Breast cancer is the most common malignant disease among women. At present, more and more attention has been paid to long non-coding RNAs (lncRNAs) in the field of breast cancer research. We aimed to investigate the expression profiles of lncRNAs and construct a prognostic lncRNA for predicting the overall survival (OS) of breast cancer.Methods: The expression profiles of lncRNAs and clinical data with breast cancer were obtained from The Cancer Genome Atlas (TCGA). Differentially expressed lncRNAs were screened out by R package (limma). The survival probability was estimated by the Kaplan‑Meier Test. The Cox Regression Model was performed for univariate and multivariate analysis. The risk score (RS) was established on the basis of the lncRNAs’ expression level (exp) multiplied regression coefficient (β) from the multivariate cox regression analysis with the following formula: RS=exp a1 * β a1 + exp a2 * β a2 +……+ exp an * β an. Functional enrichment analysis was performed by Metascape.Results: A total of 3404 differentially expressed lncRNAs were identified. Among them, CYTOR, MIR4458HG and MAPT-AS1 were significantly associated with the survival of breast cancer. Finally, The RS could predict OS of breast cancer (RS=exp CYTOR * β CYTOR + exp MIR4458HG * β MIR4458HG + exp MAPT-AS1 * β MAPT-AS1). Moreover, it was confirmed that the three-lncRNA signature could be an independent prognostic biomarker for breast cancer (HR=3.040, P=0.000).Conclusions: This study established a three-lncRNA signature, which might be a novel prognostic biomarker for breast cancer.


Author(s):  
Peng Wang ◽  
Xin Li ◽  
Yue Gao ◽  
Qiuyan Guo ◽  
Shangwei Ning ◽  
...  

Abstract LnCeVar (http://www.bio-bigdata.net/LnCeVar/) is a comprehensive database that aims to provide genomic variations that disturb lncRNA-associated competing endogenous RNA (ceRNA) network regulation curated from the published literature and high-throughput data sets. LnCeVar curated 119 501 variation–ceRNA events from thousands of samples and cell lines, including: (i) more than 2000 experimentally supported circulating, drug-resistant and prognosis-related lncRNA biomarkers; (ii) 11 418 somatic mutation–ceRNA events from TCGA and COSMIC; (iii) 112 674 CNV–ceRNA events from TCGA; (iv) 67 066 SNP–ceRNA events from the 1000 Genomes Project. LnCeVar provides a user-friendly searching and browsing interface. In addition, as an important supplement of the database, several flexible tools have been developed to aid retrieval and analysis of the data. The LnCeVar–BLAST interface is a convenient way for users to search ceRNAs by interesting sequences. LnCeVar–Function is a tool for performing functional enrichment analysis. LnCeVar–Hallmark identifies dysregulated cancer hallmarks of variation–ceRNA events. LnCeVar–Survival performs COX regression analyses and produces survival curves for variation–ceRNA events. LnCeVar–Network identifies and creates a visualization of dysregulated variation–ceRNA networks. Collectively, LnCeVar will serve as an important resource for investigating the functions and mechanisms of personalized genomic variations that disturb ceRNA network regulation in human diseases.


2020 ◽  
Vol 2020 ◽  
pp. 1-16 ◽  
Author(s):  
Yun Zhong ◽  
Zhe Liu ◽  
Dangchi Li ◽  
Qinyuan Liao ◽  
Jingao Li

Background. An increasing number of studies have indicated that the abnormal expression of certain long noncoding RNAs (lncRNAs) is linked to the overall survival (OS) of patients with myeloma. Methods. Gene expression data of myeloma patients were downloaded from the Gene Expression Omnibus (GEO) database (GSE4581 and GSE57317). Cox regression analysis, Kaplan-Meier, and receiver operating characteristic (ROC) analysis were performed to construct and validate the prediction model. Single sample gene set enrichment (ssGSEA) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used to predict the function of a specified lncRNA. Results. In this study, a seven-lncRNA signature was identified and used to construct a risk score system for myeloma prognosis. This system was used to stratify patients with different survival rates in the training set into high-risk and low-risk groups. Test set, the entire test set, the external validation set, and the myeloma subtype achieved the authentication of the results. In addition, functional enrichment analysis indicated that 7 prognostic lncRNAs may be involved in the tumorigenesis of myeloma through cancer-related pathways and biological processes. The results of the immune score showed that IF_I was negatively correlated with the risk score. Compared with the published gene signature, the 7-lncRNA model has a higher C-index (above 0.8). Conclusion. In summary, our data provide evidence that seven lncRNAs could be used as independent biomarkers to predict the prognosis of myeloma, which also indicated that these 7 lncRNAs may be involved in the progression of myeloma.


BMC Cancer ◽  
2020 ◽  
Vol 20 (1) ◽  
Author(s):  
Fuqiang Zu ◽  
Peng Liu ◽  
Huaitao Wang ◽  
Ting Zhu ◽  
Jian Sun ◽  
...  

Abstract Background It is well acknowledged that cancer-related pathways play pivotal roles in the progression of pancreatic cancer (PC). Employing Integrated analysis, we aim to identify the pathway-related ceRNA network associated with PC progression. Methods We divided eight GEO datasets into three groups according to their platform, and combined TCGA and GTEx databases as a group. Additionally, we screened out the differentially expressed genes (DEGs) and performed functional enrichment analysis in each group, and recognized the top hub genes in the most enriched pathway. Furthermore, the upstream of miRNAs and lncRNAs were predicted and validated according to their expression and prognostic roles. Finally, the co-expression analysis was applied to identify a pathway-related ceRNA network in the progression of PC. Results A total of 51 significant pathways that common enriched in all groups were spotted. Enrichment analysis indicated that pathway in cancer was greatly linked with tumor formation and progression. Next, the top 20 hug genes in this pathway were recognized, and stepwise prediction and validation from mRNA to lncRNA, including 11 hub genes, 4 key miRNAs, and 2 key lncRNAs, were applied to identify a meaningful ceRNA network according to ceRNA rules. Ultimately, we identified the PVT1/miR-20b/CCND1 axis as a promising pathway-related ceRNA axis in the progression of PC. Conclusion Overall, we elucidate the pathway-related ceRNA regulatory network of PVT1/miR-20b/CCND1 in the progression of PC, which can be considered as therapeutic targets and encouraging prognostic biomarkers for PC.


2020 ◽  
Author(s):  
Fuqiang Zu ◽  
Peng Liu ◽  
Huaitao Wang ◽  
Ting Zhu ◽  
Jian Sun ◽  
...  

Abstract Background: It is well acknowledged that cancer-related pathways play pivotal roles in the progression of pancreatic cancer (PC). Employing Integrated analysis, we aim to identify the pathway-related ceRNA network associated with PC progression.Methods: We divided eight GEO datasets into three groups according to their platform, and combined TCGA and GTEx databases as a group. Additionally, we screened out the differentially expressed genes (DEGs) and performed functional enrichment analysis in each group, and recognized the top hub genes in the most enriched pathway. Furthermore, the upstream of miRNAs and lncRNAs were predicted and validated according to their expression and prognostic roles. Finally, the co-expression analysis was applied to identify a pathway-related ceRNA network in the progression of PC.Results: A total of 51 significant pathways that common enriched in all groups were spotted. Enrichment analysis indicated that pathway in cancer was greatly linked with tumor formation and progression. Next, the top 20 hug genes in this pathway were recognized, and stepwise prediction and validation from mRNA to lncRNA, including 11 hub genes, 4 key miRNAs, and 2 key lncRNAs, were applied to identify a meaningful ceRNA network according to ceRNA rules. Ultimately, we identified the PVT1/miR-20b-/CCND1 axis as a promising pathway-related ceRNA axis in the progression of PC.Conclusion: Overall, we elucidate the pathway-related ceRNA regulatory network of PVT1/miR-20b-/CCND1 in the progression of PC, which can be considered as therapeutic targets and encouraging prognostic biomarkers for PC.


2022 ◽  
Vol 2022 ◽  
pp. 1-16
Author(s):  
Jin Zhou ◽  
Zheming Liu ◽  
Huibo Zhang ◽  
Tianyu Lei ◽  
Jiahui Liu ◽  
...  

Purpose. Recent researches showed the vital role of BACH1 in promoting the metastasis of lung cancer. We aimed to explore the value of BACH1 in predicting the overall survival (OS) of early-stage (stages I-II) lung adenocarcinoma. Patients and Methods. Lung adenocarcinoma cases were screened from the Cancer Genome Atlas (TCGA) database. Functional enrichment analysis was performed to obtain the biological mechanisms of BACH1. Gene set enrichment analysis (GSEA) was performed to identify the difference of biological pathways between high- and low-BACH1 groups. Univariate and multivariate COX regression analysis had been used to screen prognostic factors, which were used to establish the BACH1 expression-based prognostic model in the TCGA dataset. The C-index and time-dependent AUC curve were used to evaluate predictive power of the model. External validation of prognostic value was performed in two independent datasets from Gene Expression Omnibus (GEO). Decision analysis curve was finally used to evaluate clinical usefulness of the BACH1-based model beyond pathologic stage alone. Results. BACH1 was an independent prognostic factor for lung adenocarcinoma. High-expression BACH1 cases had worse OS. BACH1-based prognostic model showed an ideal C-index and t -AUC and validated by two GEO datasets, independently. More importantly, the BACH1-based model indicated positive clinical applicability by DCA curves. Conclusion. Our research confirmed that BACH1 was an important predictor of prognosis in early-stage lung adenocarcinoma. The higher the expression of BACH1, the worse OS of the patients.


2021 ◽  
Author(s):  
Xiaoyu Ji ◽  
Guangdi Chu ◽  
Jinwen Jiao ◽  
Teng Lv ◽  
Yulong Chen ◽  
...  

Abstract Objective: Cervical cancer (CC) is one of the most common types of malignant female cancer, and its incidence and mortality are not optimistic. Protein panels can be a powerful prognostic factor for many types of cancer. The purpose of our study was to investigate a proteomic panel to predict survival of patients with common CC. Methods and results: The protein expression and clinicopathological data of CC were downloaded from The Cancer Proteome Atlas (TCPA) and The Cancer Genome Atlas (TCGA) database, respectively. We selected the prognosis-related proteins (PRPs) by univariate Cox regression analysis and found that the results of functional enrichment analysis were mainly related to apoptosis. We used Kaplan–Meier(K-M) analysis and multivariable Cox regression analysis further to screen PRPs to establish a prognostic model, including BCL2, SMAD3, and 4EBP1-pT70. The signature was verified to be independent predictors of OS by Cox regression analysis and the Area Under Curves. Nomogram and subgroup classification were established based on the signature to verify its clinical application. Furthermore, we looked for the co-expressed proteins of three-protein panel as potential prognostic proteins.Conclusion: A proteomic signature independently predicted OS of CC patients, and the predictive ability was better than the clinicopathological characteristics. This signature can help improve prediction for clinical outcome and provides new targets for CC treatment.


Sign in / Sign up

Export Citation Format

Share Document