Fig. 2

The interplay between ITH and TME in the primary PDACs and matched liver metastases. A Box plots indicating the percentage differences in T cell subclusters among origins (two-sided Wilcoxon rank sum test: *P < 0.05, **P < 0.01, ***P < 0.001). B Area plots displaying the changes in T cell subcluster composition by origin for each patient. C Dot plots illustrating ligand-receptor interactions between malignant ductal cells and Tregs. The size of a circle indicates an interaction score, and the color of a circle represents the origin. D Multiplex immunohistochemistry (IHC) showing the interaction between ICAM1 (magenta)- or IGF2R (orange)-expressing FOXP3⁺ (green) Tregs (arrows) and AREG (cyan)- or IGF2 (yellow)-expressing cytokeratin (CK)⁺ (red) tumor cells. Nuclei are counterstained with DAPI (blue). E and F Scatter plot displaying the correlation between the fraction of basal-like in ductal cells and the fraction of cytotoxic T cells in T cells (E), and between the expression level of S100A9 in ductal cells and the fraction of cytotoxic T cells in T cells (F). G and H Pearson correlation between the proportion of basal-like among ductal cells and the proportion of Tregs among the T cell population (G), and between the expression level of S100A9 in ductal cells and the fraction of Tregs among T cells (H). I and J Mapping of major cell types (I) and NMF subtypes (J) to spatial transcription spots from treatment naïve PDAC patient published by Zhou et al. using a robust cell type decomposition (RCTD) method. K The spots on the spatial transcriptome slide were colored by NMF subtypes and overlaid with contour plots of Treg enriched spots. L and M Multiplex IHC showing the expression of S100A9 (yellow) and the distribution of T cells in basal-like dominant (L) and classical dominant (M) PDAC tissues. CD8 (green) for cytotoxic T cells, FOXP3 (red) for Tregs, CK (white) for ductal cancer cells, S100A2 (magenta) for basal-like ductal cells and DAPI (blue) for nuclei were co-stained. Scale bar, 50 μm. N Box plots indicating the percentage differences in myeloid subclusters among origins (two-sided Wilcoxon rank sum test: *P < 0.05, **P < 0.01, ***P < 0.001). Samples from the same patients were connected by solid lines. O Area plots showing the change in the composition of the myeloid subclusters by origin for each patient. P and Q Scatter plot displaying the Pearson correlation between the fraction of basal-like in ductal cells and the fraction of Mono-FCN1 (P) and Mp-TGFBI (Q)