Abstract
Abnormal glycolysis is one of the hallmarks of cancer and plays an important role in its development. This study was devoted to identify glycolysis related genes as prognostic biomarkers for non-small cell lung cancer (NSCLC). The mRNA expression profile and clinical follow-up data were obtained using The Cancer Genome Atlas (TCGA) database. The validation set was obtained by bootstrap method of random repeated sampling. A total of 200 glycolysis-related genes were obtained from Gene Set Enrichment Analysis (GSEA) and 46 genes were significantly associated with overall survival (OS). Five genes (PKP2, LDHA, HMMR, COL5A1 and B3GNT3) were eventually identified to calculate risk score of NSCLC patients. The univariate and multivariate Cox regression analysis indicated that the risk score was an independent prognostic factor (training set: HR=2.126, 95% CI [1.605, 2.815], p<0.001; validation set: HR=2.298, 95%CI [1.450, 3.640], p<0.001). Patients assigned to the high-risk group were associated with poor OS compared with patients in the low-risk group (training set: P=7.946e-06; validation set: P=6.368e-07). Receiver operating characteristic (ROC) curve and stratification analysis also demonstrated the potential prognostic performance. In conclusion, we constructed a novel glycolysis related risk signature which might contribute to predicting the prognosis of NSCLC.