Abstract
In hepatocellular carcinoma (HCC), classical cancer stem cells (CSC) markers were shared by normal stem cells, targeting which may hinder hepatic regeneration and cause liver failure. Additionally, the spatial structure of CSC still remained elusive. To address these limitations, we undertook a comprehensive study combining single-cell data (56,022 cells from 20 samples) and spatial data (38,191 spots from eight samples) to obtain CSC signature and uncover its spatial structure. Utilizing the CytoTRACE algorithm, we discretely identified CSC, which displayed upregulated proliferation pathways regulated by HIF1A. A CSC signature of 107 genes was then developed using Weighted Gene Co-expression Network Analysis (WGCNA). Notably, HCC patients with high CSC levels exhibited an accumulation of SPP1+ macrophages (Macro_SPP1) expressing metalloproteinases (MMP9, MMP12, and MMP7) regulated by HIF1A, suggesting a hypoxic tumor region connecting Macro_SPP1 and CSC. Both CSC and Macro_SPP1 correlated with worse prognosis and undesirable immunotherapy response. Spatial analysis revealed the co-location of CSC and Macro_SPP1, with CD8 T cells excluded from the tumor region. The co-location area and non-tumor area of boundary exhibited a high level of hypoxia, with the HAVRC2 checkpoint highly expressed. Within the co-location area, the SPP1 signaling pathway was most active in cell-cell communication, with SPP1-CD44 and SPP1-ITGA/ITGB identified as the main ligand-receptor pairs. This study successfully constructed a CSC signature and demonstrated the co-location of CSC and Macro_SPP1 in a hypoxic region that exacerbates the tumor microenvironment in HCC.