Background
Multiple factors affect the survival time of patients with lung adenocarcinoma (LUAD). Specifically, the therapeutic effect of medicines and the disease recurrence probability differs among patients with the same stage of LUAD. Thus, effective prognostic predictors need to be identified.
Methods
Based on the tumor mutation burden (TMB) data obtained by TCGA, LUAD was divided into high and low groups, and the differentially expressed glycolysis-related genes between the two groups were screened out. Cox regression was used to obtain a prognostic model. A receiver operating characteristic (ROC) curve and calibration curve were generated to evaluate the nomogram that was constructed based on clinicopathological characteristics and the risk score. Two datasets (GSE68465 and GSE11969) from Gene Expression Omnibus (GEO) were used to verify the prognostic performance of the gene. Furthermore, differences in immune cell distribution, immune-related molecules and drug susceptibility were assessed for their relationship with the risk score.
Results
We confirmed a 5-gene signature (FKBP4, HMMR, B4GALT1, ERO1L, ENO1) capable of dividing patients into two risk groups. There was a significant difference in overall survival (OS) times between the high-risk group and the low-risk group (P = 1.085e-4), with the low-risk group having a better survival outcome. Through multivariate Cox analysis, the risk score was confirmed to be an independent prognostic factor (HR = 1.289, 95% CI = 1.202-1.383, P < 0.001), and the ROC curve and nomogram exhibited accurate prediction performance. Validation of the data obtained in the GEO database yielded similar results. Additionally, there were significant differences in cisplatin, paclitaxel, gemcitabine, docetaxel, gefitiniband erlotinib sensitivity between the low-risk and high-risk groups.
Conclusions
Our results reveal that glycolysis-related gene are feasible predictors of LUAD patient survival and response to therapeutics.