Integrated analysis revealing a novel stemness-metabolism-related gene signature for predicting prognosis and immunotherapy response in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a malignant lethal tumor and both cancer stem cells (CSCs) and metabolism reprogramming have been proven to play indispensable roles in HCC. This study aimed to reveal the connection between metabolism reprogramming and the stemness characteristics of HCC, established a new gene signature related to stemness and metabolism and utilized it to assess HCC prognosis and immunotherapy response. The clinical information and gene expression profiles (GEPs) of 478 HCC patients came from the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA). The one-class logistic regression (OCLR) algorithm was employed to calculate the messenger ribonucleic acid expression-based stemness index (mRNAsi), a new stemness index quantifying stemness features. Differentially expressed analyses were done between high- and low-mRNAsi groups and 74 differentially expressed metabolism-related genes (DEMRGs) were identified with the help of metabolism-related gene sets from Molecular Signatures Database (MSigDB). After integrated analysis, a risk score model based on the three most efficient prognostic DEMRGs, including Recombinant Phosphofructokinase Platelet (PFKP), phosphodiesterase 2A (PDE2A) and UDP-glucuronosyltransferase 1A5 (UGT1A5) was constructed and HCC patients were divided into high-risk and low-risk groups. Significant differences were found in pathway enrichment, immune cell infiltration patterns, and gene alterations between the two groups. High-risk group patients tended to have worse clinical outcomes and were more likely to respond to immunotherapy. A stemness-metabolism-related model composed of gender, age, the risk score model and tumor-node-metastasis (TNM) staging was generated and showed great discrimination and strong ability in predicting HCC prognosis and immunotherapy response.