Glycolysis-related genes predict prognosis and indicate immune microenvironment features in gastric cancer

Gastric cancer (GC) is a major contributor to cancer-related mortality. Glycolysis plays a pivotal role in tumor microenvironment (TME) reprogramming. In this research, the functions of glycolysis-associated genes (GRGs) were evaluated to predict the outcome and reveal the characteristics of the immune microenvironment in individuals with stomach cancer.

The Cancer Genome Atlas (TCGA)-stomach adenocarcinoma (STAD) cohort provided gene expression and clinical data for gastric cancer (GC) patients, which were further authenticated using datasets sourced from the Gene Expression Omnibus (GEO). By referencing the Molecular Signatures Database (MSigDB), a total of 326 GRGs were pinpointed. The various subtypes of GC were outlined through consensus clustering, derived from the expression patterns of these GRGs. Utilizing multivariate Cox regression analysis, a multigene risk score model was formulated. Both the CIBERSORT and ESTIMATE algorithms played a pivotal role in assessing the immune microenvironment. To delve into the biological functions of the key genes, wound healing, transwell invasion, and MTT assays were conducted.

Based on the expression patterns of GRGs, patients were categorized into two distinct groups: the metabolic subtype, designated as cluster A, and the immune subtype, labeled as cluster B. Patients belonging to cluster B exhibited a poorer prognosis. A prognostic risk score model, formulated upon the expression levels of six key GRGs - ME1, PLOD2, NUP50, CXCR4, SLC35A3, and SRD35A3 - emerged as a viable tool for predicting patient outcomes. The downregulation of CXCR4 notably diminished the glycolytic capacity of gastric cancer (GC) cells, alongside their migratory, invasive, and proliferative capabilities. Intriguingly, despite the adverse prognostic implications associated with both the immune subtype (cluster B) and the high-risk cohort, these groups exhibited a favorable immune microenvironment coupled with elevated expression of immune checkpoint genes. Our investigations revealed a positive correlation between high CXCR4 expression and low ME1 expression with the infiltration of CD8⁺ T cells, as well as an enhanced responsiveness to treatment with an anti-PD-1 immune checkpoint inhibitor.

In this study, we discovered that the expression profiles of GRGs hold the potential to forecast the prognosis of gastric cancer (GC) patients, thereby possibly aiding in clinical treatment decision-making.