Virtually all patients with multiple myeloma become unresponsive to treatment over time. Relapsed/refractory multiple myeloma (RRMM) is accompanied by the clonal evolution of myeloma cells with heterogeneous genomic aberrations and profound changes of the bone marrow microenvironment (BME). However, the molecular mechanisms that drive drug resistance remain elusive. Here, we analyze the heterogeneous tumor cell population and its complex interaction network with the BME of 20 RRMM patients by single cell RNA-sequencing before/after treatment. Subclones with chromosome 1q-gain express a specific transcriptomic signature and frequently expand during treatment. Furthermore, RRMM cells shape an immune suppressive BME by upregulation of inflammatory cytokines and close interaction with the myeloid compartment. It is characterized by the accumulation of PD1⁺γδ T-cells and tumor-associated macrophages as well as the depletion of hematopoietic progenitors. Thus, our study resolves transcriptional features of subclones in RRMM and mechanisms of microenvironmental reprogramming with implications for clinical decision-making.

Intratumor heterogeneity as a clinical challenge becomes most evident after several treatment lines, when multidrug-resistant subclones accumulate. To address this challenge, the characterization of resistance mechanisms at the subclonal level is key to identify common vulnerabilities. In this study, we integrate whole-genome sequencing, single-cell (sc) transcriptomics (scRNA sequencing), and chromatin accessibility (scATAC sequencing) together with mitochondrial DNA mutations to define subclonal architecture and evolution for longitudinal samples from 15 patients with relapsed or refractory multiple myeloma. We assess transcriptomic and epigenomic changes to resolve the multifactorial nature of therapy resistance and relate it to the parallel occurrence of different mechanisms: (1) preexisting epigenetic profiles of subclones associated with survival advantages, (2) converging phenotypic adaptation of genetically distinct subclones, and (3) subclone-specific interactions of myeloma and bone marrow microenvironment cells. Our study showcases how an integrative multiomics analysis can be applied to track and characterize distinct multidrug-resistant subclones over time for the identification of molecular targets against them.

Early alterations within the bone marrow microenvironment that contribute to the progression of multiple myeloma (MM) from its precursor stages could be the key to identifying novel therapeutic approaches. However, the intrinsic variability in cellular populations between patients and the differences in sample processing and analysis methods have made it difficult to identify consistent changes between data sets. Here, we used single-cell RNA sequencing of bone marrow cells from precursor stages, monoclonal gammopathy of unknown significance, smoldering MM, and newly diagnosed MM and analyzed our data in combination with a previously published data set that used a similar patient population and sample processing. Despite the vast interpatient heterogeneity, some alterations were consistently observed in both data sets. We identified changes in immune cell populations as the disease progressed, which were characterized by a substantial decrease in memory and naïve CD4 T cells, and an increase in CD8+ effector T cells and T-regulatory cells. These alterations were further accompanied by an enrichment of nonclonal memory B cells and an increase in CD14 and CD16 monocytes in MM compared with its precursor stages. These results provide crucial information on the immune changes associated with the progression to clinical MM and can help to develop immune-based strategies for patient stratification and early therapeutic intervention.

In multiple myeloma spatial differences in the subclonal architecture, molecular signatures and composition of the microenvironment remain poorly characterized. To address this shortcoming, we perform multi-region sequencing on paired random bone marrow and focal lesion samples from 17 newly diagnosed patients. Using single-cell RNA- and ATAC-seq we find a median of 6 tumor subclones per patient and unique subclones in focal lesions. Genetically identical subclones display different levels of spatial transcriptional plasticity, including nearly identical profiles and pronounced heterogeneity at different sites, which can include differential expression of immunotherapy targets, such as CD20 and CD38. Macrophages are significantly depleted in the microenvironment of focal lesions. We observe proportional changes in the T-cell repertoire but no site-specific expansion of T-cell clones in intramedullary lesions. In conclusion, our results demonstrate the relevance of considering spatial heterogeneity in multiple myeloma with potential implications for models of cell-cell interactions and disease progression.

Multiple Myeloma (MM) is a malignant plasma cell disorder with an unmet medical need, in particular for relapsed and refractory patients. Molecules within deregulated signaling pathways, including the RAS/RAF/MEK/ERK, but also the PI3K/AKT-pathway belong to the most promising evolving therapeutic targets. Rationally derived compounds hold great therapeutic promise to target tumor-specific abnormalities rather than general MM-associated vulnerabilities. This paradigm is probably best depicted by targeting mutated BRAF: while well-tolerated, remarkable responses have been achieved in selected patients by inhibition of BRAFV600E alone or in combination with MEK. Targeting of AKT has also shown promising results in a subset of patients as monotherapy or to resensitize MM-cells to conventional treatment. Approaches to target transcription factors, convergence points of signaling cascades such as p53 or c-MYC, are emerging as yet another exciting strategy for pathway-directed therapy. Informed by our increasing knowledge on the impact of signaling pathways in MM pathophysiology, rationally derived Precision-Medicine trials are ongoing. Their results are likely to once more fundamentally change treatment strategies in MM.

The bone marrow microenvironment plays a crucial role in the development of multiple myeloma. As the disease progresses, malignant myeloma cells can evolve to survive outside the bone marrow. However, the processes underlying bone marrow independence and their consequences for immune control remain poorly understood. Here, we conducted single-cell and spatial multiomics analyses of bone marrow–confined intramedullary disease and paired breakout lesions that disrupt the cortical bone. These analyses revealed a distinct cellular microenvironment and architectural features of breakout lesions, characterized by extensive areas of malignant plasma cells interspersed with lesion-specific solitary natural killer and macrophage populations, as well as focal accumulations of immune cell agglomerates. Within these agglomerates, spatially confined T cell clones expanded alongside various immune cells, coinciding with the local genomic evolution of tumor cells. These analyses identify breakout lesions as a hotspot for tumor-immune cell interactions and diversification, representing a key event in myeloma pathogenesis.