Abstract
BACKGROUND: Esophageal squamous cell carcinoma (ESCC) is a highly lethal malignant tumor lacking effective treatments; 20% of ESCC patients develop liver metastasis with an extremely short survival time of ≈5 months. The tumor microenvironment (TME) plays a crucial role in tumor homeostasis, but the relationship between the ESCC TME and liver metastasis is still unknown. METHODS: To identify potential cell populations contributing to ESCC liver metastasis, single-cell RNA (scRNA) sequencing data were analyzed to identify the major cell populations within the TME. Each of the major cell populations was re-clustered to define detailed cell subsets. Thereafter, the gene set variation analysis (GSVA) score was calculated for the bulk RNA-seq data based on the gene signatures of each cell subset. The relationship between the GSVA score of each cellular subset and clinical outcome was further analyzed to identify the cellular subset associated with ESCC liver metastasis, which was validated by multiplex immunohistochemistry. RESULTS: C1QA/B(+) tumor-associated macrophages (TAMs) acted as the central regulator of the ESCC TME, closely associated with several key cell subsets. Several immune checkpoints, including CD40, CD47 and LGALS9, were all positively expressed in C1QA/B(+) macrophages, which may exert central regulatory control of immune evasion by ESCC via these immune checkpoints expressions. CONCLUSIONS: Our results comprehensively revealed the landscape of tumor-infiltrating immune cells associated with ESCC prognosis and metastasis, and suggest a novel strategy for developing immunotherapies for ESCC liver metastasis by targeting C1QA/B(+) TAMs.