Accumulating evidences indicate significant alterations in the aerobic glycolysis in clear cell renal cell carcinoma (ccRCC). We aim to develop and validate a glycolysis-related genes signature for predicting the clinical outcomes of patients with ccRCC.
mRNA expression profiling of ccRCC was obtained from The Cancer Genome Atlas database. Univariate Cox regression analysis and lasso Cox regression model were performed to identify and construct the prognostic gene signature. The protein expression levels of the core genes were obtained from the Human Protein Atlas database. We used four external independent data sets to verify the predictive power of the model for prognosis, tyrosine kinase inhibitor (TKI) therapy, and immunotherapy responses, respectively. Finally, we explored the potential mechanism of this signature through gene set enrichment analysis (GSEA).
Through the GSEA, glycolysis-related gene sets were significantly different between ccRCC tissues and normal tissues. Next, we identified and constructed a seven-mRNA signature (GALM, TGFA, RBCK1, CD44, HK3, KIF20A, and IDUA), which was significantly correlated with worse survival outcome and was an independent prognostic indicator for ccRCC patients. Furthermore, the expression levels of hub genes were validated based on the Human Protein Atlas databases. More importantly, the model can predict patients’ response to TKI therapy and immunotherapy. These findings were successfully validated in the external independent ccRCC cohorts. The mechanism exploration showed that the model may influence the prognosis by influencing tumor proliferation, base mismatch repair system and immune status of patients.
Our study has built up a robust glycolysis-based molecular signature that predicts the prognosis and TKI therapy and immunotherapy responses of patients with ccRCC with high accuracy, which might provide important guidance for clinical assessment. Also, clinical investigations in large ccRCC cohorts are greatly needed to validate our findings.