Integrative analyses of bulk and single-cell RNA-seq reveals the correlation between SPP1+ macrophages and resistance to neoadjuvant chemoimmunotherapy in esophageal squamous cell carcinoma

Neoadjuvant chemoimmunotherapy (NACI) has significant implications for the treatment of esophageal cancer. However, its clinical efficacy varies considerably among patients, necessitating further investigation into the underlying mechanisms. The rapid advancement of single-cell RNA sequencing (scRNA-seq) technology facilitates the analysis of patient heterogeneity at the cellular level, particularly regarding treatment outcomes. In this study, we first analyzed scRNA-seq data of esophageal squamous cell carcinoma (ESCC) following NACI, obtained from the Gene Expression Omnibus (GEO) database. After performing dimensionality reduction, clustering, and annotation on the scRNA-seq data, we employed CellChat to investigate differences in cell-cell communication among samples from distinct efficacy groups. The results indicated that macrophages in the non-responder exhibited stronger cell communication intensity compared to those in responders, with SPP1 and GALECTIN signals showing the most significant differences between the two groups. This finding underscores the crucial role of macrophages in the efficacy of NACI. Subsequently, reclustering of macrophages revealed that Mac-SPP1 may be primarily responsible for treatment resistance, while Mac-C1QC appears to promote T cell activation. Finally, we conducted transcriptome sequencing on ESCC tissues obtained from 32 patients who underwent surgery following NACI. Utilizing CIBERSORT, CIBERSORTx, and WGCNA, we analyzed the heterogeneity of tumor microenvironment among different efficacy groups and validated the correlation between SPP1⁺ macrophages and resistance to NACI in ESCC using publicly available transcriptome sequencing datasets. These findings suggest that SPP1⁺ macrophages may represent a key factor contributing to resistance against NACI in ESCC.