Figure 6 from Spatial Architecture of Myeloid and T Cells Orchestrates Immune Evasion and Clinical Outcome in Lung Cancer
Authors: Enfield Katey S.S. , Colliver Emma , Lee Claudia , Magness Alastair , Moore David A. , Sivakumar Monica , Grigoriadis Kristiana , Pich Oriol , Karasaki Takahiro , Hobson Philip S. , Levi Dina , Veeriah Selvaraju , Puttick Clare , Nye Emma L. , Green Mary , Dijkstra Krijn K. , Shimato Masako , Akarca Ayse U. , Marafioti Teresa , Salgado Roberto , Hackshaw Allan , Jamal-Hanjani Mariam , van Maldegem Febe , McGranahan Nicholas , Glass Benjamin , Pulaski Hanna , Walk Eric , Reading James L. , Quezada Sergio A. , Hiley Crispin T. , Downward Julian , Sahai Erik , Swanton Charles , Angelova Mihaela 
Published: June 2024

<p>Neutrophil infiltration is associated with recent subclonal expansion and poorer disease-free survival. <b>A,</b> Comparison of neutrophil cell density in primary tumors with metastasis-seeding clones detectable at time of surgery or during follow-up (metastasizing tumors) to tumors from patients who were metastasis-free and recurrence-free for more than 3 years of follow-up time (control) in LUAD (<i>n</i> = 28) and LUSC (<i>n</i> = 15) within the TRACERx discovery cohort. Maximum neutrophil density taken for tumors with multiple tumor cores. <i>P</i> values derived from one-tailed Wilcoxon test. <b>B,</b> Kaplan–Meier curves for DFS according to tumor-level TAN score in the validation cohort (<i>n</i> = 332 patients). <i>P</i> value derived from univariate Cox model adjusted for histology. <b>C,</b> Multivariable Cox proportional hazard regression analysis of DFS using tumor-level TAN score and tumor-level necrosis evaluation from H&E images of diagnostic tumor blocks. <b>D,</b> Recent subclonal expansion score, measured as the maximum cancer cell fraction of subclones at the terminus of the phylogenetic tree, compared between TAN-high and TAN-low tumor regions in LUAD and LUSC patients from the TRACERx 421 cohort. <i>P</i> values were derived from a linear mixed effects model with patient as random effect. <b>E,</b> Spatial ITH score of cell types and communities and ITH probabilities of TME classes in multiregion analysis (pan-immune <i>n</i> = 41 tumors, 112 cores). ITH score was calculated as the average standard deviation of the cell/community density in multiple regions per tumor and z-score transformed. ITH score of αSMA<sup>+</sup> cells was derived from T cells and stroma panel (<i>n</i> = 39 tumors, 105 cores). ITH probability was calculated as 1 − probability of all regions having the same indicated TME class. <b>F,</b> Example phylogenetic tree depicting a stage IIIA LUSC case with a clonal <i>PIK3CA</i> mutation, high TAN scores, and recent subclonal expansion, including in a region (R5) that seeded a lymph node (LN) metastasis (FLN, FFPE LN). The metastasis-seeding lineage is highlighted in orange. Tumor-level TAN scores and regional subclonal expansion (SubExp) scores are reported for primary tumor regions. The reported TAN score represents the maximum of the tumor nest and stroma scores. Each cluster in the phylogenetic tree is assigned a color that is also represented in the region clone maps. The clone maps illustrate the prevalence of each clone within a region. <b>G,</b> IMC images shown for R5 and the LN metastasis. Scale bar, 200 μm. <b>H,</b> Summary schematic of the link between tumors with a neutrophil enriched microenvironment with tumor progression. <b>·</b>, <i>P</i> < 0.1; *, <i>P</i> < 0.05; **, <i>P</i> < 0.01; ***, <i>P</i> < 0.001; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; ITH, intratumor heterogeneity.</p>