While the bone marrow attracts tumor cells in many solid cancers leading to poor outcome in affected patients, comprehensive analyses of bone marrow metastases have not been performed on a single-cell level. We here set out to capture tumor heterogeneity and unravel microenvironmental changes in neuroblastoma, a solid cancer with bone marrow involvement. To this end, we employed a multi-omics data mining approach to define a multiplex imaging panel and developed DeepFLEX, a pipeline for subsequent multiplex image analysis, whereby we constructed a single-cell atlas of over 35,000 disseminated tumor cells (DTCs) and cells of their microenvironment in the metastatic bone marrow niche. Further, we independently profiled the transcriptome of a cohort of 38 patients with and without bone marrow metastasis. Our results revealed vast diversity among DTCs and suggest that FAIM2 can act as a complementary marker to capture DTC heterogeneity. Importantly, we demonstrate that malignant bone marrow infiltration is associated with an inflammatory response and at the same time the presence of immuno-suppressive cell types, most prominently an immature neutrophil/granulocytic myeloid-derived suppressor-like cell type. The presented findings indicate that metastatic tumor cells shape the bone marrow microenvironment, warranting deeper investigations of spatio-temporal dynamics at the single-cell level and their clinical relevance.