Integrative single-cell transcriptomic analyses reveal the cellular ontological and functional heterogeneities of primary and metastatic liver tumors

Abstract Background The global cellular landscape of the tumor microenvironment (TME) combining primary and metastatic liver tumors has not been comprehensively characterized. Methods Based on the scRNA-seq and spatial transcriptomic data of non-tumor liver tissues (NTs), primary liver tumors (PTs) and metastatic liver tumors (MTs), we performed the tissue preference, trajectory reconstruction, transcription factor activity inference, cell–cell interaction and cellular deconvolution analyses to construct a comprehensive cellular landscape of liver tumors. Results Our analyses depicted the heterogeneous cellular ecosystems in NTs, PTs and MTs. The activated memory B cells and effector T cells were shown to gradually shift to inhibitory B cells, regulatory or exhausted T cells in liver tumors, especially in MTs. Among them, we characterized a unique group of TCF7+ CD8+ memory T cells specifically enriched in MTs that could differentiate into exhausted T cells likely driven by the p38 MAPK signaling. With regard to myeloid cells, the liver-resident macrophages and inflammatory monocyte/macrophages were markedly replaced by tumor-associated macrophages (TAMs), with TREM2+ and UBE2C+ TAMs enriched in PTs, while SPP1+ and WDR45B+ TAMs in MTs. We further showed that the newly identified WDR45B+ TAMs exhibit an M2-like polarization and are associated with adverse prognosis in patients with liver metastases. Additionally, we addressed that endothelial cells display higher immune tolerance and angiogenesis capacity, and provided evidence for the source of the mesenchymal transformation of fibroblasts in tumors. Finally, the malignant hepatocytes and fibroblasts were prioritized as the pivotal cell populations in shaping the microenvironments of PTs and MTs, respectively. Notably, validation analyses by using spatial or bulk transcriptomic data in clinical cohorts concordantly emphasized the clinical significance of these findings. Conclusions This study defines the ontological and functional heterogeneities in cellular ecosystems of primary and metastatic liver tumors, providing a foundation for future investigation of the underlying cellular mechanisms.