Abstract
Both the tumor and tumor microenvironment (TME) are crucial for pathogenesis and chemotherapy resistance in multiple myeloma (MM). Bortezomib, commonly used for MM treatment, works on both MM and TME cells, but innate and acquired resistance easily develop. By single-cell RNA sequencing (scRNA-seq), we investigated bone marrow aspirates of 18 treatment-naïve MM patients who later received bortezomib-based treatments. Twelve plasma and TME cell types and their subsets were identified. Suboptimal responders (SORs) to bortezomib exhibited higher copy number alteration burdens than optimal responders (ORs). Forty-four differentially expressed genes for SORs based on scRNA-seq data were further analyzed in an independent cohort of 90 treatment-naïve MMs, where 24 genes were validated. A combined model of three clinical variables (older age, low absolute lymphocyte count, and no autologous stem cell transplantation) and 24 genes was associated with bortezomib responsiveness and poor prognosis. In T cells, cytotoxic memory, proliferating, and dysfunctional subsets were significantly enriched in SORs. Moreover, we identified three monocyte subsets associated with bortezomib responsiveness and an MM-specific NK cell trajectory that ended with an MM-specific subset. scRNA-seq predicted the interaction of the GAS6-MERTK, ALCAM-CD6, and BAG6-NCR gene networks. Of note, tumor cells from ORs and SORs were the most prominent sources of ALCAM on effector T cells and BAG6 on NK cells, respectively. Our results indicate that the complicated compositional and molecular changes of both tumor and immune cells in the bone marrow (BM) milieu are important in the development and acquisition of resistance to bortezomib-based treatment of MM.