Identification of a combined hypoxia and lactate metabolism prognostic signature in lung adenocarcinoma

Abstract

Background In the tumor microenvironment (TME), a bidirectional relationship exists between hypoxia and lactate metabolism, with each component exerting a reciprocal influence on the other, forming an inextricable link. The aim of the present investigation was to develop a prognostic model by amalgamating genes associated with hypoxia and lactate metabolism. This model is intended to serve as a tool for predicting patient outcomes, including survival rates, the status of the immune microenvironment, and responsiveness to therapy in patients with lung adenocarcinoma (LUAD). Methods Transcriptomic sequencing data and patient clinical information specific to LUAD were obtained from comprehensive repositories of The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). A compendium of genes implicated in hypoxia and lactate metabolism was assembled from an array of accessible datasets. Univariate and multivariate Cox regression analyses were employed. Additional investigative procedures, including tumor mutational load (TMB), microsatellite instability (MSI), functional enrichment assessments and the ESTIMATE, CIBERSORT, and TIDE algorithms, were used to evaluate drug sensitivity and predict the efficacy of immune-based therapies. Results A novel prognostic signature comprising five lactate and hypoxia-related genes (LHRGs), PKFP, SLC2A1, BCAN, CDKN3, and ANLN, was established. This model demonstrated that LUAD patients with elevated LHRG-related risk scores exhibited significantly reduced survival rates. Both univariate and multivariate Cox analyses confirmed that the risk score was a robust prognostic indicator of overall survival. Immunophenotyping revealed increased infiltration of memory CD4 + T cells, dendritic cells and NK cells in patients classified within the high-risk category compared to their low-risk counterparts. Higher probability of mutations in lung adenocarcinoma driver genes in high-risk groups, and the MSI was associated with the risk-score. Functional enrichment analyses indicated a predominance of cell cycle-related pathways in the high-risk group, whereas metabolic pathways were more prevalent in the low-risk group. Moreover, drug sensitivity analyses revealed increased sensitivity to a variety of drugs in the high-risk group, especially inhibitors of the PI3K-AKT, EGFR, and ELK pathways. Conclusions This prognostic model integrates lactate metabolism and hypoxia parameters, offering predictive insights regarding survival, immune cell infiltration and functionality, as well as therapeutic responsiveness in LUAD patients. This model may facilitate personalized treatment strategies, tailoring interventions to the unique molecular profile of each patient’s disease.