SUMMARY Prognostically relevant RNA expression states exist in pancreatic ductal adenocarcinoma (PDAC), but our understanding of their drivers, stability, and relationship to therapeutic response is limited. To examine these attributes systematically, we profiled metastatic biopsies and matched organoid models at single-cell resolution. In vivo , we identify a new intermediate PDAC transcriptional cell state and uncover distinct site- and state-specific tumor microenvironments (TMEs). Benchmarking models against this reference map, we reveal strong culture-specific biases in cancer cell transcriptional state representation driven by altered TME signals. We restore expression state heterogeneity by adding back in vivo -relevant factors and show plasticity in culture models. Further, we prove that non-genetic modulation of cell state can strongly influence drug responses, uncovering state-specific vulnerabilities. This work provides a broadly applicable framework for aligning cell states across in vivo and ex vivo settings, identifying drivers of transcriptional plasticity and manipulating cell state to target associated vulnerabilities.

Highlights d Obesity accelerates oncogenic Kras-driven pancreatic ductal tumorigenesis in mice d Genetic or dietary weight loss intercepts pancreatic cancer progression d Obesity is associated with aberrant pancreatic islet cholecystokinin expression d Islet cholecystokinin overexpression drives pancreatic ductal cancer development

Immunotherapies targeting aspects of T cell functionality are efficacious in many solid tumors, but pancreatic ductal adenocarcinoma (PDAC) remains refractory to these treatments. Deeper understanding of the PDAC immune ecosystem is needed to identify additional therapeutic targets and predictive biomarkers for therapeutic response and resistance monitoring. To address these needs, we quantitatively evaluated leukocyte contexture in 135 human PDACs at single-cell resolution by profiling density and spatial distribution of myeloid and lymphoid cells within histopathologically defined regions of surgical resections from treatment-naive and presurgically (neoadjuvant)–treated patients and biopsy specimens from metastatic PDAC. Resultant data establish an immune atlas of PDAC heterogeneity, identify leukocyte features correlating with clinical outcomes, and, through an in silico study, provide guidance for use of PDAC tissue microarrays to optimally measure intratumoral immune heterogeneity. Atlas data have direct applicability as a reference for evaluating immune responses to investigational neoadjuvant PDAC therapeutics where pretherapy baseline specimens are not available. Significance: We provide a phenotypic and spatial immune atlas of human PDAC identifying leukocyte composition at steady state and following standard neoadjuvant therapies. These data have broad utility as a resource that can inform on leukocyte responses to emerging therapies where baseline tissues were not acquired. This article is highlighted in the In This Issue feature, p. 1861

Purpose. While abundant myeloid cell populations in the pancreatic ductal adenocarcinoma (PDAC) microenvironment have been postulated to suppress anti-tumor immunity, the composition of these populations, their spatial locations, and how they relate to patient outcomes are poorly understood. Experimental Design. To generate spatially-resolved tumor and immune cell data at single cell resolution, we developed two quantitative multiplex immunofluorescence assays to interrogate myeloid cells (CD15, CD14, ARG1, CD33, HLA-DR) and macrophages [CD68, CD163, CD86, interferon regulatory factor 5 (IRF5), MRC1 (CD206)] in the PDAC tumor microenvironment. Spatial point pattern analyses were conducted to assess the degree of co-localization between tumor cells and immune cells. Multivariable-adjusted Cox proportional hazards regression was used to assess associations with patient outcomes. Results. In a multi-institutional cohort of 305 primary PDAC resection specimens, myeloid cells were abundant, enriched within stromal regions, highly heterogeneous across tumors, and differed by somatic genotype. High densities of CD15 + ARG1 + immunosuppressive granulocytic cells and M2-polarized macrophages were associated with worse patient survival. Moreover, beyond cell density, closer proximity of M2-polarized macrophages to tumor cells was strongly associated with disease-free survival, revealing the clinical significance and biologic importance of immune cell localization within tumor areas. Conclusions. A diverse set of myeloid cells are present within the PDAC tumor microenvironment and are distributed heterogeneously across patient tumors. Not only the densities but also the spatial locations of myeloid immune cells are associated with patient Research.