Fig. 3

Immune cell landscape of TNBC with respect to age at diagnosis. a CIBERSORTx estimated B cell proportions per sample versus stratified age groups in all cases (left), HRDetect-high, HRDetect-low, and PAM50 basal-like cases (right). Top axes indicate group sizes. Two-sided p values calculated using Kruskal-Wallis test. Linear regression modeling showing p value and slope coefficient (k) when using B cell proportion and continuous age in the model. b Heatmap of 102 immune cell marker genes in 232 SCAN-B cases using Pearson correlation and ward.D linkage. Hierarchical clustering of cases (columns) and genes (rows) was performed using Pearson correlation as distance metric and ward.D linkage on mean-centered log2 transformed data with an offset of 0.1. c TIL percentage estimated from whole-slide hemotoxylin and eosin-stained sections versus stratified age groups in all cases (left), HRDetect-high, HRDetect-low, and PAM50 basal-like cases (right). Top axes indicate group sizes. Two-sided p values calculated using Kruskal-Wallis test. Linear regression modeling showing p value and slope coefficient (k) when using TIL percentage and continuous age in the model. For age group definitions, “[” equals ≥, “)” equals <, and “]” equals ≤ for the value specified next to it. In panels a and c, separate results of a sensitivity analysis for trend are reported in red due to very small sample numbers across the full stratified age range