Single-cell sequencing reveals the immune microenvironment landscape related to anti-PD-1 resistance in metastatic colorectal cancer with high microsatellite instability

Abstract

Abstract Background The objective response rate of microsatellite instability-high (MSI-H) metastatic colorectal cancer (mCRC) patients with first-line anti-programmed cell death protein-1 (PD-1) monotherapy is only 40–45%. Single-cell RNA sequencing (scRNA-seq) enables unbiased analysis of the full variety of cells comprising the tumor microenvironment. Thus, we used scRNA-seq to assess differences among microenvironment components between therapy-resistant and therapy-sensitive groups in MSI-H/mismatch repair-deficient (dMMR) mCRC. Resistance-related cell types and genes identified by this analysis were subsequently verified in clinical samples and mouse models to further reveal the molecular mechanism of anti-PD-1 resistance in MSI-H or dMMR mCRC. Methods The response of primary and metastatic lesions to first-line anti-PD-1 monotherapy was evaluated by radiology. Cells from primary lesions of patients with MSI-H/dMMR mCRC were analyzed using scRNA-seq. To identify the marker genes in each cluster, distinct cell clusters were identified and subjected to subcluster analysis. Then, a protein‒protein interaction network was constructed to identify key genes. Immunohistochemistry and immunofluorescence were applied to verify key genes and cell marker molecules in clinical samples. Immunohistochemistry, quantitative real-time PCR, and western blotting were performed to examine the expression of IL-1β and MMP9. Moreover, quantitative analysis and sorting of myeloid-derived suppressor cells (MDSCs) and CD8+ T cells were performed using flow cytometry. Results Tumor responses in 23 patients with MSI-H/dMMR mCRC were evaluated by radiology. The objective response rate was 43.48%, and the disease control rate was 69.57%. ScRNA-seq analysis showed that, compared with the treatment-resistant group, the treatment-sensitive group accumulated more CD8+ T cells. Experiments with both clinical samples and mice indicated that infiltration of IL-1β-driven MDSCs and inactivation of CD8+ T cells contribute to anti-PD-1 resistance in MSI-H/dMMR CRC. Conclusions CD8+ T cells and IL-1β were identified as the cell type and gene, respectively, with the highest correlation with anti-PD-1 resistance. Infiltration of IL-1β-driven MDSCs was a significant factor in anti-PD-1 resistance in CRC. IL-1β antagonists are expected to be developed as a new treatment for anti-PD-1 inhibitor resistance.