Tumor immune microenvironmental alterations occur early in multiple myeloma (MM) development. In this study, we aim to systematically characterize the tumor immune microenvironment (TME) and the tumor-immune interactions from precursor stages, i.e., monoclonal gammopathy of undetermined significance (MGUS) and smoldering MM (SMM), to newly diagnosed MM, comparing these to healthy donors. CIBERSORT analysis revealed proportions of 27 cell types, including 10 innate immune cells (monocytes, M0/M1/M2 macrophages, resting/activated dendritic cells, resting/activated mast cells, eosinophils, and neutrophils), 7 T cell subsets (CD8+ T cells, naive CD4+ T cells, resting/activated memory CD4+ T cells, follicular helper T cells, regulatory T cells, and gamma delta T cells), resting/activated NK cells, naive B cells, memory B cells, plasma cells, myeloma plasma cells, memory plasma cells, osteoblasts, osteoclasts and adipocytes. We removed plasma cells (PC) from the 100% total cell proportions to avoid bias resulting from different tumor burdens. We found that the proportions of neutrophils, mast cells, and monocytes account for the majority of innate immune cells and were decreased in NDMM and its precursor stages when compared with NBM. We noticed a significant abundance of M2 macrophages in NDMM, but not in the precursor stages. The proportion of CD8+ T cells was increased at the SMM and MM stages, and the proportion of activated memory CD4+ T cells continued to decrease from NBM to MGUS, SMM, and MM. We performed correlation analysis of immune cell proportions with the time of progression of MGUS and SMM patients. We observed neutrophil proportions were negatively correlated with the time to progression in SMM patients, indicating more neutrophils predict inferior outcomes for SMM patients. However, in NDMM patients, we observed that the neutrophil percentage was decreased in patients with high-risk status based on the 70-gene Prognostic Risk Score (GEP-70) or at stage III of International Staging System (ISS), and patients with high proportions of neutrophils had a significantly superior overall survival (OS) and event-free survival (EFS). For T cell populations, we observed that the proportions of CD8+ T cells were negatively correlated with the time of progression in SMM patients and the γδ T cells were positively correlated. We observed the proportions of γδ T cells were also positively correlated with the time of progression in MGUS patients. Consistent with the time-to-progression data, high-risk SMM patients showed higher proportions of CD8+ T cells and lower proportions of naïve CD4+ T cells and γδ T cells.
Overall design
Bone marrow biopsies were processed for RNA extraction and hybridization on Affymetrix microarrays. CIBERSORT-based analysis was used to deconvolution the bone marrow cell types. Specifically, the 27-myeloma genome signature matrix (MGSM27), a gene-expression–based computational technique, was implemented.