Abstract
Multiple myeloma (MM) progression is driven by immune dysregulation within the tumor microenvironment (TME). However, myeloma-intrinsic mechanisms underlying immune dysfunction remain poorly defined, and current immunotherapies show limited efficacy. Using RNA-seq data from 859 MM patients (MMRF-CoMMpass), we integrated xCELL, CIBERSORT, and ESTIMATE algorithms to deconvolute immune-stromal dynamics. Consensus clustering identified immune subtypes, followed by differential gene analysis and LASSO-Cox regression to construct a prognostic model validated in an independent cohort (GSE19784, N = 328). Immune Subtype Classification: Two subgroups emerged: Multiple myeloma-associated immune-related cluster 1 (N = 482): Immune-dysfunctional TME with Th2 cell enrichment, preadipocyte accumulation, and CXCL family suppression, linked to poor survival (P < 0.001). Multiple myeloma-associated immune-related cluster 2 (N = 377): Immune-active TME with cytotoxic CD8 + T/NK cell infiltration and favorable outcomes. Prognostic Gene Signature: Ten immune-related genes (UBE2T, E2F2, EXO1, SH2D2A, DRP2, WNT9A, SHROOM3, TMC8, CDCA7, and GPR132) predicted survival (The One-year AUC = 0.682 and The Over 5-years AUC = 0.714). We define a myeloma-intrinsic immune classification system and a 10-gene prognostic index, offering a framework for risk-stratified immunotherapy. Integration with flow cytometry could optimize precision treatment in MM.