(A) Stacked bar plot showing the percentage composition of cell types in benign (12,400 cells) and premetastatic niche bone marrow (18,259 cells). (B) Scatterplot of myeloid cell fraction of benign bone marrow (N = 7), premetastatic niche bone marrow (N = 8), tumor-adjacent bone marrow (N = 8), and metastatic tumor (N = 9). (C) Hypergeometric test of significantly (P < 0.05) upregulated genes in myeloid cells of the premetastatic niche compared with benign bone marrow. Heatmaps showing (D) IFNG and receptor expression (IFNGR1, IFNGR2) across all cells and (E) IFNG expression of lymphoid cells grouped by tissue type. Scatterplots of average expression of (F) IFNGR1, (G) IL10RA, and (H) IL10 in myeloid cells across bone marrow and metastatic tumor samples. Bars represent mean ± SEM. Cell proportions and differential gene expression analysis were conducted using RAISIN (version 1) to generate log fold-change and P values to determine significance.