Abstract
Hepatocellular carcinoma (HCC) is a common malignancy with a poor prognosis, and its heterogeneity affects the response to clinical treatments. Glycolysis is highly associated with HCC therapy and prognosis. The present study aimed to identify a novel biomarker for HCC by exploring the heterogeneity of glycolysis in HCC. The intersection of both marker genes of glycolysis‑related cell clusters from single‑cell RNA sequencing analysis and mRNA data of liver HCC from The Cancer Genome Atlas were used to construct a prognostic model through Cox proportional hazard regression and the least absolute shrinkage and selection operator Cox regression. Data from the International Cancer Genome Consortium were used to validate the results of the analysis. Immune status analysis was then conducted. A significant gene in the prognostic model was identified as a potential biomarker and was verified through in vitro experiments. The results revealed that the glycolysis‑related prognostic model divided patients with HCC into high‑ and low‑risk groups. A nomogram combining the model and clinical features exhibited accurate predictive ability, with an area under the curve of 0.763 at 3 years. The high‑risk group exhibited a higher expression of checkpoint genes and lower tumor immune dysfunction and exclusion scores, suggesting that this group may be more likely to benefit from immunotherapy. The tumor tissues had a higher zinc finger protein (ZFP)41 mRNA and protein expression compared with the adjacent tissues. In vitro analyses revealed that ZFP41 played a crucial role in cell viability, proliferation, migration, invasion and glycolysis. On the whole, the present study demonstrates that the glycolysis‑related prognostic gene, ZFP41, is a potential prognostic biomarker and therapeutic target, and may play a crucial role in glycolysis and malignancy in HCC.