Single-cell transcriptomics reveals proliferative, immune-evasive, and bone marrow–independent transcriptional programs in extramedullary multiple myeloma.

Abstract

Abstract Background Extramedullary disease (EMD) in multiple myeloma (MM) represents an aggressive phenotype with poor outcomes. We performed single-cell RNA sequencing (scRNA-seq) on clonal plasma cells (CPCs) from EMD tissue, solitary plasmacytomas (SPM), and bone marrow aspirates (BMAs) from patients with MM to identify transcriptional states unique to EMD. Methods We processed FFPE tumor samples from EMD and SPMs and cryopreserved CD138+ BMAs from MM, using the 10x Genomics FLEX platform. Cells/nuclei were fixed, hybridized with Human WTA probes, pooled, and sequenced (10,000–40,000 cells/library) on NovaSeq X. After removing low quality cells, data was normalized and integrated using RPCA in Seurat. Plasma cell identity was assigned via Cellarium and canonical markers. Clonality was inferred by light chain restrictions and refined via random forest modeling. Cell cycle effects were regressed out to reduce proliferation bias. Results We analyzed 128, 908 CPCs from 16 samples [8 EMD (62, 669 CPCs), 4 SPM (20, 607 CPCs), 4 BMA (45, 632 CPCs)] and identified eight transcriptionally distinct cell clusters. Differential expression and pathway comparisons were performed between EMD/SPM/BMA within each of these clusters. Integration of gene expression and pathway analyses revealed consistent upregulation of cell cycle and E2F pathways in EMD. We identified two highly proliferative cell cycle clusters (Cluster 2 with S-phase genes and Cluster 4 with a G2/M transcriptional program) that were expanded in EMD compared to BMA and SPM. Cluster 0 comprised a quiescent CPC population enriched for RNA processing pathways alongside downregulation of ER stress and autophagy. EMD cells demonstrated suppressed interleukin signaling [fold change (FC) –33.8), and class I MHC antigen processing (FC –23.5) relative to BMA, suggestive of an immune suppressive transcriptional program. Angiogenesis-related genes [VEGFA (log2FC 2.4), ADM2 (log2FC 2.2)] and EZH2 [log2FC 2.1; q<0.0001)] were upregulated in EMD. Compared to SPM, EMD plasma cells in this cluster demonstrated marked downregulation of CD74 (log2FC –2.6; q<0.0001), a key modulator of MHC class II antigen presentation. Cluster 1 comprised of metabolically stressed secretory CPCs and displayed downregulation of early growth response pathways in EMD compared to BMA, concurrent with upregulation of E2F signaling. Suppression of IL-2/STAT5, IL-6/JAK/STAT3, and TGF-β pathways in EMD suggests reduced reliance on bone marrow–derived cytokines and impaired immunoregulatory signaling. CD44 was upregulated in EMD versus BMA (log2FC 3.2), consistent with its role in EMD spread and marrow independence. Cluster 3 was enriched for extracellular matrix–interacting plasma cells (COL1A1, MMP2) with upregulated epithelial–mesenchymal transition (EMT) pathways in EMD compared to both BMA (FC 7.5) and SPM (FC 7.07). Cluster 6 was defined by a migratory, motility-primed plasma cell population (CD44, AHNAK, LMNA and S100A10). Notably, IL32 was overexpressed in EMD compared to BMA (log2FC 5.1; adj. p = 2.02×10⁻¹⁴), a finding linked to hypoxia-driven proliferation. Compared to SPM, EMD cells also had elevated MYC/E2F activity and DNA repair, hallmarks of a proliferative phenotype. Cluster 5 displayed a strong interferon-stimulated gene signature (ISG15, IFIT1, TNFSF10), with enrichment of type I IFN, JAK–STAT, and TNF pathways. Within this inflammatory cluster, EMD exhibited marked suppression of immune signaling, metabolism (TCA, OXPHOS), and TP53 programs, consistent with an immune-evasive, metabolically suppressed, and deregulated state. Finally, within Cluster 7, a smaller subset defined by NFκB–driven anti-apoptotic transcriptional activity, EMD exhibited significant downregulation of KLF2 (log2FC –3.7; adj. p = 1.77×10⁻²⁷), likely reflecting disruption of the KDM3A–KLF2–IRF4 axis, a transcriptional network critical for bone marrow homing and plasma cell survival. We consistently observed downregulation of the MAPK pathway in EMD compared to BMA, aligning with prior bulk RNA-seq data and contrasting with the near-universal presence of MAPK mutations, highlighting a disconnect between mutational status and pathway activation. Conclusion In the largest yet reported single-cell transcriptional atlas of EMD and SPM, we demonstrate EMD to be highly proliferative with metabolic reprogramming, EMT, and downregulation of key immunologic signaling pathways, indicating immune escape and microenvironmental independence.