Spatially resolved analysis of pancreatic cancer identifies therapy-associated remodeling of the tumor microenvironment

In combination with cell-intrinsic properties, interactions in the tumor microenvironment modulate therapeutic response. We leveraged single-cell spatial transcriptomics to dissect the remodeling of multicellular neighborhoods and cell–cell interactions in human pancreatic cancer associated with neoadjuvant chemotherapy and radiotherapy. We developed spatially constrained optimal transport interaction analysis (SCOTIA), an optimal transport model with a cost function that includes both spatial distance and ligand–receptor gene expression. Our results uncovered a marked change in ligand–receptor interactions between cancer-associated fibroblasts and malignant cells in response to treatment, which was supported by orthogonal datasets, including an ex vivo tumoroid coculture system. We identified enrichment in interleukin-6 family signaling that functionally confers resistance to chemotherapy. Overall, this study demonstrates that characterization of the tumor microenvironment using single-cell spatial transcriptomics allows for the identification of molecular interactions that may play a role in the emergence of therapeutic resistance and offers a spatially based analysis framework that can be broadly applied to other contexts. Spatial molecular imaging analysis of human pancreatic adenocarcinomas describes multicellular neighborhoods in the tumor microenvironment. Ligand–receptor analysis using optimal transport-based SCOTIA identifies interleukin-6 as a mediator of chemoresistance.