Abstract
Glycogen metabolism plays a key role in the development of hepatocellular carcinoma (HCC), but the function of glycogen metabolism genes in the tumor microenvironment (TME) is still to be elucidated. Single-cell RNA-seq data were obtained from ten HCC tumor samples totaling 64,545 cells, and 65 glycogen metabolism genes were analyzed by a nonnegative matrix factorization (NMF). The prognosis and immune response of new glycogen TME cell clusters were predicted by using HCC and immunotherapy cohorts from public databases. HCC single-cell analysis was divided into fibroblasts, NT T cells, macrophages, endothelial cells, and B cells, which were separately divided into new cell clusters by glycogen metabolism gene annotation. Pseudo-temporal trajectory analysis demonstrated the temporal differentiation trajectory of different glycogen subtype cell clusters. Cellular communication analysis revealed extensive interactions between endothelial cells with glycogen metabolizing TME cell-related subtypes and different glycogen subtype cell clusters. SCENIC analysis of transcription factors upstream of TME cell clusters with different glycogen metabolism. In addition, TME cell clusters of glycogen metabolism were found to be enriched in expression in CAF subtypes, CD8 depleted, M1, and M2 types. Bulk-seq analysis showed the prognostic significance of glycogen metabolism-mediated TME cell clusters in HCC, while a significant immune response was found in the immunotherapy cohort in patients treated with immune checkpoint blockade (ICB), especially for CAFs, T cells, and macrophages. In summary, our study reveals for the first time that glycogen metabolism mediates intercellular communication in the hepatocellular carcinoma microenvironment while elucidating the anti-tumor mechanisms and immune prognostic responses of different subtypes of cell clusters.