Dynamic single-cell RNA-seq analysis reveals distinct tumor program associated with microenvironmental remodeling and drug sensitivity in multiple myeloma

Background Multiple myeloma (MM) is a hematological malignancy characterized by clonal proliferation of malignant plasma cells. Despite extensive research, molecular mechanisms in MM that drive drug sensitivity and clinic outcome remain elusive. Results Single-cell RNA sequencing was applied to study tumor heterogeneity and molecular dynamics in 10 MM individuals before and after 2 cycles of bortezomib–cyclophosphamide–dexamethasone (VCD) treatment, with 3 healthy volunteers as controls. We identified that unfolded protein response and metabolic-related program were decreased, whereas stress-associated and immune reactive programs were increased after 2 cycles of VCD treatment. Interestingly, low expression of the immune reactive program by tumor cells was associated with unfavorable drug response and poor survival in MM, which probably due to downregulation of MHC class I mediated antigen presentation and immune surveillance, and upregulation of markers related to immune escape. Furthermore, combined with immune cells profiling, we uncovered a link between tumor intrinsic immune reactive program and immunosuppressive phenotype in microenvironment, evidenced by exhausted states and expression of checkpoint molecules and suppressive genes in T cells, NK cells and monocytes. Notably, expression of YBX1 was associated with downregulation of immune activation signaling in myeloma and reduced immune cells infiltration, thereby contributed to poor prognosis. Conclusions We dissected the tumor and immune reprogramming in MM during targeted therapy at the single-cell resolution, and identified a tumor program that integrated tumoral signaling and changes in immune microenvironment, which provided insights into understanding drug sensitivity in MM.