normalization factor was applied to PLA2G2D data, and values were converted into conventional fold changes by setting the value of a control sample to 1 as indicated, scaling all other values proportionally. Primer sequences used: PLA2G2D forward: 5′-CAACCCAAAGATGCCACG-3′ and reverse: 5′-CCAGCTTCCCTTGTCAGAG-3′. ACTB forward: 5′-ATTGCCGACAGGATGCAGAA-3′ and reverse: 5′-GCTG ATCCACATCTGCTGGAA-3′. HPRT1 forward: 5′-GACCAGTCAACAGGGG ACAT-3′ and reverse: 5′-AACACTTCGTGGGGTCCTTTTC-3′. Statistical analysis For IMC data analyses, differential cell abundances were calculated using a two-sided t -test (Rstatix v.0.7.2; https://rpkgs.datanovia.com/ rstatix/). Wilcoxon signed-rank tests and unpaired t -tests were imple- mented for non-normally and normally distributed data, respectively, using the stat_compare_means() function included in ggpubr v.0.6.0. Correlations were calculated as Pearson’s correlation coefficient ( r ) as implemented by ggpubr v0.6.0 (https://rpkgs.datanovia.com/ggpubr/) using the stat_cor function. For DSP whole-transcriptome analysis, DEGs were calculated using a LMM as implemented in the GeomxTools v.3.10.0 R package. Statistical analyses of the flow cytometry, tumour growth curves, survival analysis and IHC and immunofluorescence results were performed in GraphPad Prism v.10. Error bars reported for bar graphs in the figures denote s.e.m. unless indicated otherwise. For box plots, boxes indicate the median and interquartile range, with whiskers representing the data distribution excluding outliers. Paired or unpaired t -tests were used for comparisons of paired data with Gaussian distributions, and paired or unpaired Wilcoxon tests were performed for non-normally distributed paired data. P values (and adjusted P values) smaller than 0.05 were considered statisti- cally significant; exact P values are provided in the Supplementary Data. Data points used in the figures were obtained from independ- ent biological samples. When representative data are shown, simi- lar results were obtained from one or more independent replicate experiments. Reporting summary Further information on research design is available in the Nature Port- folio Reporting Summary linked to this article. Data availability The PLA2G2D blocking antibodies developed for this study are made available to the scientific community and can be requested from L.-F.L. and K.V.L. Raw flow-cytometry, IHC and immunofluorescence data files generated and analysed for the current study are available upon reasonable request. IMC, GeoMx and single-cell transcriptomic data have been deposited in publicly accessible databases (the Gene Expression Omnibus or Zenodo): IMC data: https://doi.org/10.5281/ zenodo.20038486 (ref. 72); GeoMx data: GSE282437; bulk RNA-seq data: https://doi.org/10.5281/zenodo.20314670 (ref. 73); and single-cell RNA-seq data: GSE333464. Source data are provided with this paper. Code availability Computer code used for data analysis and figure plotting is available through GitHub: https://github.com/pulmed/Spatial-proteogenomics- of-tumour-draining-lymph-nodes-.git. 48. NSCLC Meta-analyses Collaborative Group. Adjuvant chemotherapy, with or without postoperative radiotherapy, in operable non-small-cell lung cancer: two meta-analyses of individual patient data. Lancet 375 , 1267–1277 https://doi.org/10.1016/S0140- 6736(10)60059-1 (2010). 49. Windhager, J. et al. An end-to-end workflow for multiplexed image processing and analysis. Nat. Protoc. 18 , 3565–3613 (2023). 50. Stirling, D. R. et al. CellProfiler 4: improvements in speed, utility and usability. BMC Bioinformatics 22 , 433 (2021). 51. Berg, S. et al. ilastik: interactive machine learning for (bio)image analysis. Nat. Methods 16 , 1226–1232 (2019). 52. Greenwald, N. F. et al. Whole-cell segmentation of tissue images with human-level performance using large-scale data annotation and deep learning. Nat. Biotechnol. 40 , 555–565 (2022). 53. Haghverdi, L., Lun, A. T. L., Morgan, M. D. & Marioni, J. C. Batch effects in single-cell RNA- sequencing data are corrected by matching mutual nearest neighbors. Nat. Biotechnol. 36 , 421–427 (2018). 54. Levine, J. H. et al. Data-driven phenotypic dissection of AML reveals progenitor-like cells that correlate with prognosis. Cell 162 , 184–197 (2015). 55. Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. Nat. Methods 9 , 676–682 (2012). 56. van Maldegem, F. et al. Characterisation of tumour microenvironment remodelling following oncogene inhibition in preclinical studies with imaging mass cytometry. Nat. Commun. 12 , 5906 (2021). 57. Eling, N., Damond, N., Hoch, T. & Bodenmiller, B. cytomapper: an R/Bioconductor package for visualization of highly multiplexed imaging data. Bioinformatics 36 , 5706–5708 (2020). 58. Zhou, Y. et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10 , 1523 (2019). 59. Xu, S. et al. Using clusterProfiler to characterize multiomics data. Nat. Protoc. 19 , 3292–3320 (2024). 60. Qian, J. et al. A pan-cancer blueprint of the heterogeneous tumor microenvironment revealed by single-cell profiling. Cell Res. 30 , 745–762 (2020). 61. Cheng, S. et al. A pan-cancer single-cell transcriptional atlas of tumor infiltrating myeloid cells. Cell 184 , 792–809 (2021). 62. Zhang, Q. et al. Landscape and dynamics of single immune cells in hepatocellular carcinoma. Cell 179 , 829–845 (2019). 63. Tang, Z., Kang, B., Li, C., Chen, T. & Zhang, Z. GEPIA2: an enhanced web server for large- scale expression profiling and interactive analysis. Nucleic Acids Res. 47 , W556–W560 (2019). 64. Bankhead, P. et al. QuPath: open source software for digital pathology image analysis. Sci. Rep. 7 , 16878 (2017). 65. Vuylsteke, R. J. et al. Sampling tumor-draining lymph nodes for phenotypic and functional analysis of dendritic cells and T cells. Am. J. Pathol. 161 , 19–26 (2002). 66. Elliott, B. et al. Successful live cell harvest from bisected sentinel lymph nodes research report. J. Immunol. Methods 291 , 71–78 (2004). 67. Belk, J. A. et al. Genome-wide CRISPR screens of T cell exhaustion identify chromatin remodeling factors that limit T cell persistence. Cancer Cell 40 , 768–786 (2022). 68. Miller, J. C. et al. Deciphering the transcriptional network of the dendritic cell lineage. Nat. Immunol. 13 , 888–899 (2012). 69. Kloosterman, D. J. et al. Macrophage-mediated myelin recycling fuels brain cancer malignancy. Cell 187 , 5336–5356 (2024). 70. Onrust-van Schoonhoven, A. et al. 3D chromatin reprogramming primes human memory T H 2 cells for rapid recall and pathogenic dysfunction. Sci. Immunol. 8 , eadg3917 (2023). 71. Heinz, S. et al. Simple combinations of lineage-determining transcription factors prime cis -regulatory elements required for macrophage and B cell identities. Mol. Cell 38 , 576–589 (2010). 72. van Krimpen, A., Stadhouders, R., Dammeijer, F. & van Nimwegen, M. PLA2G2D in tumour- draining lymph nodes regulates anti-tumour immunity. Zenodo https://doi.org/10.5281/ zenodo.20038486 (2026). 73. Stadhouders, R., Dammeijer, F. & Eterman, M. PLA2G2D in tumour-draining lymph nodes regulates anti-tumour immunity. Zenodo https://doi.org/10.5281/zenodo.20314670 (2026). Acknowledgements We thank members of the Pulmonary Medicine department at the Erasmus MC for discussions; T. Cupedo, I. Touw and the Haematology department of the Erasmus MC for setting up the IMC infrastructure in-house; the Erasmus MC Biomics core facility for assisting with library preparation and sequencing of GeoMx DSP samples; G. van Beek for initial 10X single-cell RNA-seq data processing; N. Reticker-Flynn for the LN-metastatic B16F0-derived LN6-987AL line; and T. van Hall for the H-2Kb TRP2 tetramer. Author contributions A.v.K., D.J.G., C.V., J.G.J.V.A., L.-F.L., K.V.L., R.S. and F.D. designed the study. J.C.J., V.G. and A.L. collected and analysed patient samples and clinical data. R.v.W. collected cLNs. J.v.d.T. and A.M. provided advice on TDLN histology, ROI selection and image data quantification. A.v.K. and M.v.N. performed IMC experiments. A.O.-V.S. designed GeoMx DSP experiments and performed data analysis. A.O.-V.S. and D.V. performed GeoMx DSP experiments. A.v.K. analysed IMC, GeoMx DSP and single-cell RNA-seq data. F.v.M. provided advice on IMC cell segmentation and data analysis. V.G. and N.R. performed flow-cytometry and cell-culture experiments. T.P.P.v.d.B created TMAs and performed histological staining. V.G. and F.D. performed cell quantification and analysis of histological images. R.S. and F.D. analysed TCGA data. M.E. performed in vivo experiments and single-cell transcriptome experiments. M.E. and H.L. performed in vitro studies with primary human T cells. E.B. performed single-cell transcriptomics experiments. M.U., J.H., X.J., K.L., A.Z.L., Y.-R.L., L.-F.L. and K.V.L. generated PLA2G2D blocking antibodies and performed in vivo experiments and in vitro studies. M.M. provided key reagents. Y.N. performed sPLA2 enzyme activity assays. R.W.H. provided advice on experimental design. C.M.B., S.V. and B.N.L. performed in vivo depletion experiments. A.v.K., L.-F.L., K.V.L., R.S. and F.D. wrote the manuscript with input from all authors. M.E., V.G., M.U. and J.C.J. made equal contributions as second authors. Funding This work was supported by grants from the Daniel den Hoed Foundation and the Support Casper Foundation; an unrestricted grant from the Van Herk Foundation; and the National Cancer Institute of the National Institutes of Health under award number R43CA278134, which partially supported J.H., L.-F.L., K.V.L. and X.J. F.D. is supported by a VENI grant (09150162510193) from the Dutch Research Council (NWO/ZonMw). R.S. is supported by