Single-cell landscape of bone marrow metastases in human neuroblastoma unraveled by deep multiplex imaging

Bone marrow commonly serves as a metastatic niche for disseminated tumor cells (DTCs) of solid cancers in patients with unfavorable clinical outcome. Single-cell assessment of bone marrow metastases is essential to decipher the entire spectrum of tumor heterogeneity in these cancers, however, has previously not been performed. Here we used multi-epitope-ligand cartography (MELC) to spatially profile 20 biomarkers and assess morphology in DTCs as well as hematopoietic and mesenchymal cells of eight bone marrow metastases from neuroblastoma patients. We developed DeepFLEX, a single-cell image analysis pipeline for MELC data that combines deep learning-based cell and nucleus segmentation and overcomes frequent challenges of multiplex imaging methods including autofluorescence and unspecific antibody binding. Using DeepFLEX, we built a single-cell atlas of bone marrow metastases comprising more than 35,000 single cells. Comparisons of cell type proportions between samples indicated that microenvironmental changes in the metastatic bone marrow are associated with tumor cell infiltration and therapy response. Hierarchical clustering of DTCs revealed multiple phenotypes with highly diverse expression of markers such as FAIM2, an inhibitory protein in the Fas apoptotic pathway, which we propose as a complementary marker to capture DTC heterogeneity in neuroblastoma. The presented single-cell atlas provides first insights into the heterogeneity of human bone marrow metastases and is an important step towards a deeper understanding of DTCs and their interactions with the bone marrow niche.