BackgroundThis study aimed to examine glycolysis/gluconeogenesis-related genes in hepatocellular carcinoma (HCC) and evaluate their potential roles in HCC progression and immunotherapy response.MethodsData analyzed in this study were collected from GSE14520, GSE76427, GSE174570, The Cancer Genome Atlas (TCGA), PXD006512, and GSE149614 datasets, metabolic pathways were collected from MSigDB database. Differentially expressed genes (DEGs) were identified between HCC and controls. Differentially expressed glycolysis/gluconeogenesis-related genes (candidate genes) were obtained and consensus clustering was performed based on the expression of candidate genes. Bioinformatics analysis was used to evaluate candidate genes and screen prognostic genes. Finally, the key results were tested in HCC patients.ResultsThirteen differentially expressed glycolysis/gluconeogenesis-related genes were validated in additional datasets. Consensus clustering analysis identified two distinct patient clusters (C1 and C2) with different prognoses and immune microenvironments. Immune score and tumor purity were significantly higher in C1 than in C2, and CD4+ memory activated T cell, Tfh, Tregs, and macrophage M0 were higher infiltrated in HCC and C1 group. The study also identified five intersecting DEGs from candidate genes in TCGA, GSE14520, and GSE141198 as prognostic genes, which had a protective role in HCC patient prognosis. Compared with the control group, the prognostic genes all showed decreased expression in HCC patients in RT-qPCR and Western blot analyses. Flow cytometry verified the abnormal infiltration level of immune cells in HCC patients.ConclusionResults showed that glycolysis/gluconeogenesis-related genes were associated with patient prognosis, immune microenvironment, and response to immunotherapy in HCC. It suggests that the model based on five prognostic genes may valuable for predicting the prognosis and immunotherapy response of HCC patients.