Abstract
Cancer immunotherapy has shown to be a promising method in treating hepatocellular carcinoma (HCC), but suboptimal responses in patients are attributed to cellular and molecular heterogeneity. Iron metabolism-related genes (IRGs) are important in maintaining immune system homeostasis and have the potential to help develop new strategies for HCC treatment. Herein, we constructed and validated the iron-metabolism gene prognostic index (IPX) using univariate Cox proportional hazards regression and LASSO Cox regression analysis, successfully categorizing HCC patients into two groups with distinct survival risks. Then, we performed single-sample gene set enrichment analysis, weighted correlation network analysis, gene ontology enrichment analysis, cellular lineage analysis, and SCENIC analysis to reveal the key determinants underlying the ability of this model based on bulk and single-cell transcriptomic data. We identified several driver transcription factors specifically activated in specific malignant cell sub-populations to contribute to the adverse survival outcomes in the IPX-high subgroup. Within the tumor microenvironment (TME), T cells displayed significant diversity in their cellular characteristics and experienced changes in their developmental paths within distinct clusters identified by IPX. Interestingly, the proportion of Treg cells was increased in the high-risk group compared with the low-risk group. These results suggest that iron-metabolism could be involved in reshaping the TME, thereby disrupting the cell cycle of immune cells. This study utilized IRGs to construct a novel and reliable model, which can be used to assess the prognosis of patients with HCC and further clarify the molecular mechanisms of IRGs in HCC at single-cell resolution.