Single-cell RNA-sequencing atlas reveals an MDK-dependent immunosuppressive environment in ErbB pathway-mutated gallbladder cancer

Our previous genomic whole-exome sequencing (WES) data identified the key ErbB pathway mutations that play an essential role in regulating the malignancy of gallbladder cancer (GBC). Herein, we tested the hypothesis that individual cellular components of the tumor microenvironment (TME) in GBC function differentially to participate in ErbB pathway mutation-dependent tumor progression.We engaged single-cell RNA-sequencing to reveal transcriptomic heterogeneity and intercellular crosstalk from 13 human GBCs and adjacent normal tissues. In addition, we performed WES analysis to reveal the genomic variations related to tumor malignancy. A variety of bulk RNA-sequencing, immunohistochemical staining, immunofluorescence staining and functional experiments were employed to study the difference between tissues with or without ErbB pathway mutations.We identified 16 cell types from a total of 114,927 cells, in which epithelial cells, M2 macrophages, and regulatory T cells were predominant in tumors with ErbB pathway mutations. Furthermore, epithelial cell subtype 1, 2 and 3 were mainly found in adenocarcinoma and subtype 4 was present in adenosquamous carcinoma. The tumors with ErbB pathway mutations harbored larger populations of epithelial cell subtype 1 and 2, and expressed higher levels of secreted midkine (MDK) than tumors without ErbB pathway mutations. Increased MDK resulted in an interaction with its receptor LRP1, which is expressed by tumor-infiltrating macrophages, and promoted immunosuppressive macrophage differentiation. Moreover, the crosstalk between macrophage-secreted CXCL10 and its receptor CXCR3 on regulatory T cells was induced in GBC with ErbB pathway mutations. Elevated MDK was correlated with poor overall survival in patients with GBC.This study has provided valuable insights into transcriptomic heterogeneity and the global cellular network in the TME, which coordinately functions to promote the progression of GBC with ErbB pathway mutations; thus, unveiling novel cellular and molecular targets for cancer therapy.We employed single-cell RNA-sequencing and functional assays to uncover the transcriptomic heterogeneity and intercellular crosstalk present in gallbladder cancer. We found that ErbB pathway mutations reduced anti-cancer immunity and led to cancer development. ErbB pathway mutations resulted in immunosuppressive macrophage differentiation and regulatory T cell activation, explaining the reduced anti-cancer immunity and worse overall survival observed in patients with these mutations.