Construction and validation of a prognostic model for lung adenocarcinoma based on endoplasmic reticulum stress-related genes.
Lung adenocarcinoma (LUAD) is one of the most universal types of cancer all over the world and its morbidity continues to rise year by year. Growing evidence has demonstrated that endoplasmic reticulum stress is highly activated in cancer cells and plays a key role in regulating the fate of cancer cells. However, the role and mechanism of endoplasmic reticulum stress in lung adenocarcinoma genesis and development remains unclear. In this research, we developed a prognostic model to predict the overall survival of patients with LUAD utilizing endoplasmic reticulum stress-related genes and screened out potential small molecular compounds, which could assist the clinician in making accurate decisions and better treat LUAD patients. Firstly, we downloaded 419 endoplasmic reticulum stress-related genes (ERSRGs) from Molecular Signatures Database (MSigDB). Secondly, we obtained information about the transcriptome profiling and corresponding clinical data of 59 normal samples and 535 lung adenocarcinoma samples from The Cancer Genome Atlas (TCGA) database. Next, we used the DESeq2 package to identify differentially expressed genes related to endoplasmic reticulum stress. We performed univariate Cox, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analysis to establish a prognostic model for LUAD patients based on ERSRGs. Then, we carried out univariate and multivariate independent prognostic analysis of endoplasmic reticulum stress-related gene (ERSRG) score and some clinical traits of lung adenocarcinoma. Additionally, we developed a clinically applicable nomogram for predicting survival for LUAD patients over one, three, and five years. Moreover, we carried out a drug sensitivity analysis to identify novel small molecule compounds for LUAD treatment. Finally, we examined the tumor microenvironment (TME) and immune cell infiltrating analysis to explore the interactions between immune and cancer cells. 142 differentially expressed ERSRGs were identified by using the DESeq2 package. A prognostic model was built based on 7 differentially expressed ERSRGs after performing univariate Cox regression, LASSO regression, and multivariate Cox regression analysis. According to the results of univariate and multivariate independent prognostic analysis, we found ERSRG score can be used as an independent prognostic maker. Using the Kaplan-Meier curves, we found low-risk patients had higher survival probability than high-risk patients in both training set and test set. A nomogram was drawn to predict 1-, 3-, and 5-year survival probability. The calibration curves explained good performance of the model for the prediction of survival. Phenformin, OSU-03012, GSK-650394 and KIN001-135 were identified as the drugs most likely to provide important information to clinicians about the treatment of LUAD patients. A prognostic prediction model was established based on 7 differentially expressed ERSRGs (PDX1, IGFBP1, DDIT4, PPP1R3G, CFTR, DERL3 and NUPR1), which could effectively predict the prognosis of LUAD patients and give a reference for clinical doctors to help LUAD patients to make better treatment tactics. Based on the 4 small molecule compounds (Phenformin, OSU-03012, GSK-650394 and KIN001-135) we discovered, targeting endoplasmic reticulum stress-related genes may also be a therapeutic approach for LUAD patients.
Li F
,Niu Y
,Zhao W
,Yan C
,Qi Y
... -
《Scientific Reports》
Development and validation of endoplasmic reticulum stress-related eight-gene signature for predicting the overall survival of lung adenocarcinoma.
The high case-fatality rate of patients with lung adenocarcinoma (LUAD) emphasizes the importance of identifying a robust and reliable prognostic signature for LUAD patients. Endoplasmic reticulum (ER) stress results from protein misfolding imbalance and has been shown to participate in the development of cancer. We aimed to develop and validation a reliable and robust ER stress-related prognostic signature to accurately predict prognosis for patients with LUAD.
The mRNA expressions data and the clinical information were downloaded from The Cancer Genome Atlas (TCGA) as training set. The data of external validation sets were downloaded from GEO database with the accession number GSE 30219, GSE 31210, GSE 50081 and GSE 37745. Univariate Cox regression analyses was performed to identify mRNAs associated with overall survival (OS) in LUAD. ER-associated genes were retrieved using GeneCards database. Next, we construct the best risk score model by the least absolute shrinkage and selection operator (LASSO) regression with tenfold cross-validation. Subsequently, predictive models and risk scores were developed in the TCGA training dataset. Cox proportional hazards regression models were used for univariate and multivariate analysis of risk score and clinicopathologic characteristics. As a validation set GSE30219, GSE31210 and (GSE50081+GSE37745) were used to validate the predictive performance of the model in TCGA. Finally, functional enrichment analysis, including the gene ontology (GO) enrichment analysis, the Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathways and gene set enrichment analysis (GSEA) were performed to further explore function and mechanisms.
A prognostic prediction model based on eight genes was developed in the TCGA training dataset. As expected, in validation sets, patients with higher risk scores were found to have worse prognosis. Time-dependent ROC curve analyses demonstrated that the risk score model was reliable. The nomograms showed excellent predictive ability. Multivariate Cox regression analyses indicated that the risk score was an independent prognostic factor for LUAD. Additionally, functional enrichment analysis showed that the relevant biomarkers were enriched in cell cycle and glycolysis related signaling pathways.
The 8-gene signature may enable improved the prediction of clinical events and decisions about management of LUAD.
Lin L
,Zhang W
《-》
Construction of an endoplasmic reticulum stress-related signature in lung adenocarcinoma by comprehensive bioinformatics analysis.
Lung Adenocarcinoma (LUAD) is a major component of lung cancer. Endoplasmic reticulum stress (ERS) has emerged as a new target for some tumor treatments.
The expression and clinical data of LUAD samples were downloaded from The Cancer Genome Atlas (TCGA) and The Gene Expression Omnibus (GEO) database, followed by acquiring ERS-related genes (ERSGs) from the GeneCards database. Differentially expressed endoplasmic reticulum stress-related genes (DE-ERSGs) were screened and used to construct a risk model by Cox regression analysis. Kaplan-Meier (K-M) curves and receiver operating characteristic (ROC) curves were plotted to determine the risk validity of the model. Moreover, enrichment analysis of differentially expressed genes (DEGs) between the high- and low- risk groups was conducted to investigate the functions related to the risk model. Furthermore, the differences in ERS status, vascular-related genes, tumor mutation burden (TMB), immunotherapy response, chemotherapy drug sensitivity and other indicators between the high- and low- risk groups were studied. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) was used to validate the mRNA expression levels of prognostic model genes.
A total of 81 DE-ERSGs were identified in the TCGA-LUAD dataset, and a risk model, including HSPD1, PCSK9, GRIA1, MAOB, COL1A1, and CAV1, was constructed by Cox regression analysis. K-M and ROC analyses showed that the high-risk group had a low survival, and the Area Under Curve (AUC) of ROC curves of 1-, 3- and 5-years overall survival was all greater than 0.6. In addition, functional enrichment analysis suggested that the risk model was related to collagen and extracellular matrix. Furthermore, differential analysis showed vascular-related genes FLT1, TMB, neoantigen, PD-L1 protein (CD274), Tumor Immune Dysfunction and Exclusion (TIDE), and T cell exclusion score were significantly different between the high- and low-risk groups. Finally, qRT-PCR results showed that the mRNA expression levels of 6 prognostic genes were consistent with the analysis.
A novel ERS-related risk model, including HSPD1, PCSK9, GRIA1, MAOB, COL1A1, and CAV1, was developed and validated, which provided a theoretical basis and reference value for ERS-related fields in the study and treatment of LUAD.
Wang Y
,Nie J
,Dai L
,Hu W
,Han S
,Zhang J
,Chen X
,Ma X
,Tian G
,Wu D
,Zhang Z
,Long J
,Fang J
... -
《BMC Pulmonary Medicine》