Abstract Granulocytic myeloid-derived suppressor cells (G-MDSCs) promote tumor growth and immunosuppression in multiple myeloma (MM). However, their phenotype is not well established for accurate monitoring or clinical translation. We aimed to provide the phenotypic profile of G-MDSCs based on their prognostic significance in MM, immunosuppressive potential, and molecular program. The preestablished phenotype of G-MDSCs was evaluated in bone marrow samples from controls and MM patients using multidimensional flow cytometry; surprisingly, we found that CD11b+CD14−CD15+CD33+HLADR− cells overlapped with common eosinophils and neutrophils, which were not expanded in MM patients. Therefore, we relied on automated clustering to unbiasedly identify all granulocytic subsets in the tumor microenvironment: basophils, eosinophils, and immature, intermediate, and mature neutrophils. In a series of 267 newly diagnosed MM patients (GEM2012MENOS65 trial), only the frequency of mature neutrophils at diagnosis was significantly associated with patient outcome, and a high mature neutrophil/T-cell ratio resulted in inferior progression-free survival (P < .001). Upon fluorescence-activated cell sorting of each neutrophil subset, T-cell proliferation decreased in the presence of mature neutrophils (0.5-fold; P = .016), and the cytotoxic potential of T cells engaged by a BCMA×CD3-bispecific antibody increased notably with the depletion of mature neutrophils (fourfold; P = .0007). Most interestingly, RNA sequencing of the 3 subsets revealed that G-MDSC–related genes were specifically upregulated in mature neutrophils from MM patients vs controls because of differential chromatin accessibility. Taken together, our results establish a correlation between the clinical significance, immunosuppressive potential, and transcriptional network of well-defined neutrophil subsets, providing for the first time a set of optimal markers (CD11b/CD13/CD16) for accurate monitoring of G-MDSCs in MM.

The historical lack of preclinical models reflecting the genetic heterogeneity of multiple myeloma (MM) hampers the advance of therapeutic discoveries. To circumvent this limitation, we screened mice engineered to carry eight MM lesions (NF-κB, KRAS, MYC, TP53, BCL2, cyclin D1, MMSET/NSD2 and c-MAF) combinatorially activated in B lymphocytes following T cell-driven immunization. Fifteen genetically diverse models developed bone marrow (BM) tumors fulfilling MM pathogenesis. Integrative analyses of ∼500 mice and ∼1,000 patients revealed a common MAPK–MYC genetic pathway that accelerated time to progression from precursor states across genetically heterogeneous MM. MYC-dependent time to progression conditioned immune evasion mechanisms that remodeled the BM microenvironment differently. Rapid MYC-driven progressors exhibited a high number of activated/exhausted CD8+ T cells with reduced immunosuppressive regulatory T (Treg) cells, while late MYC acquisition in slow progressors was associated with lower CD8+ T cell infiltration and more abundant Treg cells. Single-cell transcriptomics and functional assays defined a high ratio of CD8+ T cells versus Treg cells as a predictor of response to immune checkpoint blockade (ICB). In clinical series, high CD8+ T/Treg cell ratios underlie early progression in untreated smoldering MM, and correlated with early relapse in newly diagnosed patients with MM under Len/Dex therapy. In ICB-refractory MM models, increasing CD8+ T cell cytotoxicity or depleting Treg cells reversed immunotherapy resistance and yielded prolonged MM control. Our experimental models enable the correlation of MM genetic and immunological traits with preclinical therapy responses, which may inform the next-generation immunotherapy trials.