Pancreatic ductal adenocarcinoma (PDAC) has been classified into classical and basal-like transcriptional subtypes by bulk RNA measurements. However, recent work has uncovered greater complexity to transcriptional subtypes than was initially appreciated using bulk RNA expression profiling. To provide a deeper understanding of PDAC subtypes, we developed a multiplex immunofluorescence (mIF) pipeline that quantifies protein expression of six PDAC subtype markers (CLDN18.2, TFF1, GATA6, KRT17, KRT5, and S100A2) and permits spatially resolved, single-cell interrogation of pancreatic tumors from resection specimens and core needle biopsies. Both primary and metastatic tumors displayed striking intratumoral subtype heterogeneity that was associated with patient outcomes, existed at the scale of individual glands, and was significantly reduced in patient-derived organoid cultures. Tumor cells co-expressing classical and basal markers were present in > 90% of tumors, existed on a basal-classical polarization continuum, and were enriched in tumors containing a greater admixture of basal and classical cell populations. Cell–cell neighbor analyses within tumor glands further suggested that co-expressor cells may represent an intermediate state between expression subtype poles. The extensive intratumoral heterogeneity identified through this clinically applicable mIF pipeline may inform prognosis and treatment selection for patients with PDAC. Significance: A high-throughput pipeline using multiplex immunofluorescence in pancreatic cancer reveals striking expression subtype intratumoral heterogeneity with implications for therapy selection and identifies co-expressor cells that may serve as intermediates during subtype switching.

Purpose: Neoadjuvant chemotherapy is increasingly administered to patients with resectable or borderline resectable pancreatic ductal adenocarcinoma (PDAC), yet its impact on the tumor immune microenvironment is incompletely understood. Experimental: Design: We employed quantitative, spatially resolved multiplex immunofluorescence and digital image analysis to identify T-cell subpopulations, macrophage polarization states, and myeloid cell subpopulations in a multi-institution cohort of up-front resected primary tumors (n = 299) and in a comparative set of resected tumors after FOLFIRINOX-based neoadjuvant therapy (n = 36) or up-front surgery (n = 30). Multivariable-adjusted Cox proportional hazards models were used to evaluate associations between the immune microenvironment and patient outcomes. Results: In the multi-institutional resection cohort, immune cells exhibited substantial heterogeneity across patient tumors and were located predominantly in stromal regions. Unsupervised clustering using immune cell densities identified four main patterns of immune cell infiltration. One pattern, seen in 20% of tumors and characterized by abundant T cells (T cell–rich) and a paucity of immunosuppressive granulocytes and macrophages, was associated with improved patient survival. Neoadjuvant chemotherapy was associated with a higher CD8:CD4 ratio, greater M1:M2–polarized macrophage ratio, and reduced CD15+ARG1+ immunosuppressive granulocyte density. Within neoadjuvant-treated tumors, 72% showed a T cell–rich pattern with low immunosuppressive granulocytes and macrophages. M1-polarized macrophages were located closer to tumor cells after neoadjuvant chemotherapy, and colocalization of M1-polarized macrophages and tumor cells was associated with greater tumor pathologic response and improved patient survival. Conclusions: Neoadjuvant chemotherapy with FOLFIRINOX shifts the PDAC immune microenvironment toward an anti-tumorigenic state associated with improved patient survival.