The bone marrow immune microenvironment (BMME) is a critical factor in multiple myeloma (MM) pathogenesis. Despite higher incidence of MM among self-identified African Americans, the relationship between African genetic similarity and the BMME remains unclear. We present the largest analysis to date assessing the impact of Ancestry on the BMME in myeloma (N=320, including 64 with ≥50% African genetic similarity, AFR-High) leveraging our recently generated Immune Atlas. AFR-High was associated with enrichment of late T-effector cells, driven by rare cytotoxic CD4⁺ T cells, suggesting a distinct adaptive BMME. AFR-High exhibited expansion of late-differentiated adaptive NK cells with high cytotoxic GNLY and GZMH expression. AFR-High depicted reduced CD14⁺ monocytes, increased CD16⁺ non-classical monocytes, and expansion of macrophage-like monocyte populations with enrichment of metabolic adaptation, efferocytosis and phagocytic clearance programs, and immunoregulatory pathways. These findings reveal a shift from an innate inflammatory immune landscape toward adaptive cytotoxic and immunoregulatory remodeling associated with African genetic similarity in MM.
Haplotype-resolved genome assemblies of BJ and IMR-90 human fibroblast cell lines reveal extensive structural variation and enable reanalysis of historical sequencing data
Nucleic Acids Research | April 29, 2026 | Volume: 54 | Issue: 8 | Pages: No page available
Abstract We present chromosome-level, phased diploid genome assemblies of two widely used human fibroblast cell lines: BJ (46,XY) and IMR-90 (46,XX). Using Oxford Nanopore, PacBio HiFi, and Hi-C sequencing data, we generated assemblies spanning 5.9 and 6.0 Gbp with diploid quality values exceeding QV 60. To validate structural integrity, we developed KaryoScope, an alignment-free tool for generating computational karyograms from k-mer feature databases. We identify >50 000 structural variants relative to T2T-CHM13v2.0, the majority of which are heterozygous and cell-line-specific. Combining reference-based and de novo gene annotation, we uncover a previously unreported 1 Mbp homozygous duplication at the 16p11.2 locus in BJ, demonstrating that even karyotypically normal cell lines can harbor clinically relevant submicroscopic rearrangements. We show that mapping publicly available short-read, RNA-seq, and ChIP-seq data to sample-matched diploid assemblies substantially improves read alignment and enables haplotype phasing of 23%–28% of short reads. The BJ and IMR-90 assemblies and associated variant calls are publicly available as a resource for the research community.
Abstract 3887: Multiple myeloma risk classification of CoMMpass using Consensus Genomic Staging identifies differences between patients with high and low African genetic ancestry.
Cancer Research | April 3, 2026 | Volume: 86 | Issue: 7_Supplement | Pages: 3887-3887
Abstract Background: Multiple myeloma is the second most common adult blood cancer and remains incurable, despite recent advances. The IMWG/IMS recently published updated risk criteria, Consensus Genomic Staging (CGS), combining genomic events with clinical measures to identify high-risk patients (Avet-Loiseau, et al. JCO 2025). We applied CGS risk criteria to baseline samples from the Multiple Myeloma Research Foundation CoMMpass Study (NCT01454297), a prospective, longitudinal observational study that enrolled more than1100 participants between 2011-2017. Methods: We assessed patient risk following CGS criteria. Criteria 1: del(17p), with a cutoff of >20% clonal fraction, and/or TP53 mutation. Criteria 2: an IgH translocation including t(4;14), t(14;16), or t(14;20) along with 1q+ and/or del(1p32). Criteria 3: monoallelic del(1p32) along with 1q+ or biallelic del(1p32). Criteria 4: beta-2 microglobulin >= 5.5 mg/L with normal creatinine. Patients meeting one or more criteria were classified as High Risk. To test if risk showed association with genetic ancestry, we grouped patients into those with high African genetic ancestry and those with low African genetic ancestry. Genetic ancestry was measured using ADMIXTURE in comparison to 1000 Genomes. Results: We evaluated the CGS risk status of 877 CoMMpass patients, finding 69.2% (573/828) to be Standard Risk (SR) and 30.8% (255/828) to be High Risk (HR). HR patients showed worse progression-free survival (median PFS 26.3 months vs. 42.5 months; log rank test, p < 0.0001) and overall survival (median OS 63.9 months vs. median not reached; p < 0.0001). 10.6% (91/858) of patients met Criteria 1, 10.7% (93/873) met Criteria 2, 5.6% (49/877) met Criteria 3, and 8.7% (73/835) met Criteria 4. We found 14.2% (123/867) of patients to have high African-like genetic ancestry (AFR-high) and 85.8% (744/867) to have low African-like genetic ancestry (AFR-low). We observed significantly fewer HR patients within the AFR-high group (20%, 23/115) than in the AFR-low group (32.4%, 231/712) (Fisher’s exact test, p = 0.0065). The AFR-high group met HR criteria at a lower rate across all criteria, but only with Criteria 3 was the difference significant (p-value = 0.03). Higher levels of creatinine in AFR-high patients may have impacted meeting the normal creatinine requirement of Criteria 4. Conclusions: HR patients comprise 30.8% of CoMMpass subjects and show significantly worse PFS and OS than SR patients. AFR-high patients met HR criteria at a lower rate than AFR-low patients according to the CGS criteria. Clinical trials that utilize CGS criteria to determine patient eligibility may be challenged to enroll a representative patient population, including AFR-high patients. Screening and eligibility criteria should be carefully considered to mitigate potential enrollment disparities. Citation Format: Steven M. Foltz, Chaitanya R. Acharya, Alexander Gout, Yuxin Jin, David E. Avigan, Ravi Vij, Shaji Kunnathu Kumar, Sagar Lonial, Hearn Cho, Katie Wozniak, Jonathan J. Keats, Craig Cole, John D. Carpten, George J. Mulligan. Multiple myeloma risk classification of CoMMpass using Consensus Genomic Staging identifies differences between patients with high and low African genetic ancestry [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 3887.
Longitudinal Profiling of Tumor and Immune Compartments Uncovers Patterns of Dysregulation and Associations with Response in Multiple Myeloma
Blood Cancer Discovery | December 9, 2025 | Volume: 7 | Issue: 2 | Pages: 266-286
Abstract Multiple myeloma is a malignancy of clonally expanded plasma cells shaped by complex interactions with the immune microenvironment (IME). To investigate immune correlates of treatment response and disease progression, we conducted multi-omics profiling including CD138neg single-cell RNA sequencing of 243 bone marrow samples from 102 patients (631,226 cells) and CD138pos bulk RNA and whole-genome sequencing from 209 samples. In longitudinal analyses, interferon-γ signaling associated with markers of impaired T-cell memory after autologous stem cell transplant, whereas naïve B-cell abundance and immunoglobulin diversity correlated with improved progression-free survival (HR = 0.48; P = 2.3e−4). At disease progression, multiple myeloma cells upregulated cancer–testis antigens (CTAg) and immune effector genes, with concurrent B-cell depletion, enrichment of myeloid-derived suppressor cell expression, and phenotypic T-cell exhaustion. These findings highlight dynamic immune–tumor interactions, identifying naïve B-cell reconstitution as a biomarker of durable response and CTAgs as potential targets for high-risk disease at progression. Significance: Longitudinal profiling of multiple myeloma and the IME revealed dynamic immune–tumor interactions across the disease course. Altered expression in CD8+ T cells limited memory phenotype after transplant, whereas naïve B-cell recovery associated with sustained treatment response. At progression, CTAg expression associated with immunosuppression, revealing novel mechanisms of immune dysregulation.
International Myeloma Society/International Myeloma Working Group Consensus Recommendations on the Definition of High-Risk Multiple Myeloma
Journal of Clinical Oncology | June 9, 2025 | Volume: 43 | Issue: 24 | Pages: 2739-2751
Despite significant improvements in survival of patients with multiple myeloma (MM), outcomes remain heterogeneous, and a significant proportion of patients experience suboptimal outcomes. Importantly, traditional prognostic factors based on data from patients treated with older therapies no longer capture prognosis accurately in the contemporary era of novel triplet or quadruplet therapies. Therefore, risk stratification requires refinement in the context of available and investigational treatment options in routine practice and clinical trials, respectively. The current identification of high-risk MM (HRMM) in routine practice is based on the Revised International Staging System, which stratifies patients using a combination of widely available serum biomarkers and chromosomal abnormalities assessed via fluorescence in situ hybridization. In recent years, a substantial body of evidence concerning additional clinical, biological, and molecular/genomic prognostic factors has accumulated, along with new MM risk stratification tools and consensus reports. The International Myeloma Society, along with the International Myeloma Working Group, convened an Expert Panel with the primary aim of revisiting the definition of HRMM and formulating a practical and data-driven consensus definition, based on new evidence from molecular/genomic assays, updated clinical data, and contemporary risk stratification concepts. The Panel proposes the following Consensus Genomic Staging (CGS) of HRMM which relies upon the presence of at least one of these abnormalities: (1) del(17p), with a cutoff of >20% clonal fraction, and/or TP53 mutation; (2) an IgH translocation including t(4;14), t(14;16), or t(14;20) along with 1q+ and/or del(1p32); (3) monoallelic del(1p32) along with 1q+ or biallelic del(1p32); or (4) β2 microglobulin ≥5.5 mg/L with normal creatinine (<1.2 mg/dL).
Longitudinal multi-omic profiling uncovers immune escape and predictors of response in multiple myeloma
No container title available | June 1, 2025 | Volume: No volume available | Issue: No issue available | Pages: No page available
Abstract Multiple myeloma (MM) is an incurable malignancy of clonally expanded plasma cells shaped by complex interactions with the immune microenvironment. To investigate immune factors driving treatment response and resistance, we conducted multi-omics profiling including CD138 neg single-cell RNA sequencing of 243 bone marrow samples from 102 patients (631,226 cells) and CD138 pos bulk RNA and whole-genome sequencing from 209 samples. Longitudinal analyses revealed that interferon gamma signaling impairs T cell memory after autologous stem cell transplant, while naïve B cell abundance and immunoglobulin diversity correlated with improved progression-free survival (HR = 0.48, p = 2.3e -4 ). At disease progression, MM cells upregulated cancer-testis antigens and immune effector genes, with concurrent B cell depletion, enrichment of myeloid-derived suppressor cell genes in monocytes, and T cell exhaustion. These findings highlight dynamic immune-tumor interactions, identifying naïve B cell reconstitution as a biomarker of durable response, and cancer-testis antigens as potential targets for high-risk disease at progression. Statement of Significance Longitudinal profiling of multiple myeloma and the immune microenvironment revealed dynamic immune-tumor interactions across the disease course. Dysfunctional CD8⁺ T cells limited memory formation post-transplant, while naïve B recovery associated with sustained treatment response. At progression, cancer-testis antigen expression associated with immunosuppression, revealing novel mechanisms of immune escape.
ETV4-Dependent Transcriptional Plasticity Maintains
MYC
Expression and Results in IMiD Resistance in Multiple Myeloma
Blood Cancer Discovery | November 7, 2023 | Volume: 5 | Issue: 1 | Pages: 56-73
Abstract Immunomodulatory drugs (IMiD) are a backbone therapy for multiple myeloma (MM). Despite their efficacy, most patients develop resistance, and the mechanisms are not fully defined. Here, we show that IMiD responses are directed by IMiD-dependent degradation of IKZF1 and IKZF3 that bind to enhancers necessary to sustain the expression of MYC and other myeloma oncogenes. IMiD treatment universally depleted chromatin-bound IKZF1, but eviction of P300 and BRD4 coactivators only occurred in IMiD-sensitive cells. IKZF1-bound enhancers overlapped other transcription factor binding motifs, including ETV4. Chromatin immunoprecipitation sequencing showed that ETV4 bound to the same enhancers as IKZF1, and ETV4 CRISPR/Cas9-mediated ablation resulted in sensitization of IMiD-resistant MM. ETV4 expression is associated with IMiD resistance in cell lines, poor prognosis in patients, and is upregulated at relapse. These data indicate that ETV4 alleviates IKZF1 and IKZF3 dependency in MM by maintaining oncogenic enhancer activity and identify transcriptional plasticity as a previously unrecognized mechanism of IMiD resistance. Significance: We show that IKZF1-bound enhancers are critical for IMiD efficacy and that the factor ETV4 can bind the same enhancers and substitute for IKZF1 and mediate IMiD resistance by maintaining MYC and other oncogenes. These data implicate transcription factor redundancy as a previously unrecognized mode of IMiD resistance in MM. See related article by Welsh, Barwick, et al., p. 34. See related commentary by Yun and Cleveland, p. 5. This article is featured in Selected Articles from This Issue, p. 4
Perspectives on the Risk-Stratified Treatment of Multiple Myeloma
Blood Cancer Discovery | June 2, 2022 | Volume: 3 | Issue: 4 | Pages: 273-284
Summary: The multiple myeloma treatment landscape has changed dramatically. This change, paralleled by an increase in scientific knowledge, has resulted in significant improvement in survival. However, heterogeneity remains in clinical outcomes, with a proportion of patients not benefiting from current approaches and continuing to have a poor prognosis. A significant proportion of the variability in outcome can be predicted on the basis of clinical and biochemical parameters and tumor-acquired genetic variants, allowing for risk stratification and a more personalized approach to therapy. This article discusses the principles that can enable the rational and effective development of therapeutic approaches for high-risk multiple myeloma.
Transcriptional Plasticity Compensates for Ikaros and Aiolos Proteasomal Degradation and Mediates Resistance to IMiDs in Multiple Myeloma (MM)
Abstract Background: Immunoglobulin (IGH, IGL, IGK) and non-immunoglobulin (PVT1, TXNDC5, FAM46C, DUSP22, etc.) enhancers hijacking by variable genes (MYC, MAF, MAFB, CCND1/2/3, MMSET, IRF4) is a recognized oncogenic driver event in MM. However, the identity of the transcription factors (TFs) or transcriptional regulatory complexes binding and regulating the activity of these enhancers remains to be fully elucidated and may yield valuable therapeutic targets. As such the discovery of the BET family member BRD4 as the master histone acetyl mark reader at enhancers loci regulating MYC lead to promising therapeutic developments in MM and numerous other cancers. Immunomodulatory drugs (IMiDs) promote the proteasomal degradation of IKAROS (IKZF1) and AIOLOs (IKZF3) leading to the transcriptional repression of MYC and the suppression of MM cells survival and proliferation. However, acquired resistance to IMIDs and the loss of the transcriptional repression of MYC are nearly universal and occur in spite of sustained IKZF1/3 degradation suggesting that transcriptional rewiring may be sustaining hijacked enhancers activity and transcription of driver oncogenes. Methods and Results: In order to define how IMiDs repress MYC transcription, we first defined IKZF1, BRD4, the lysine acetyl transferase P300 and the mediator complex subunit MED1 mapping within the MM genome using ChIPseq. In MM cell lines (MM1S, RPMI8226, ARP1 and AMO1), IKZF1 predominantly mapped to intronic and intergenic loci which are typically enriched with enhancer and superenhancer elements. Indeed, IKZF1 mapping to the genome nearly completely (96.5%) overlapped that of P300, MED1 and BRD4 co-occupied enhancer and superenhancer loci. We also confirmed that in the MM1S sensitive cell lines IMiDs (lenalidomide 10 μM, 24h) exposure efficiently depleted IKZF1, BRD4, P300 and MED1 at enhancer loci with ensuing MYC (and MAF) downregulation. In contrast, in resistant cell lines (RPMI8226) and in spite efficient IKZF1 displacement, BRD4, P300 and MED1 were retained at the oncogenic enhancer (IGLL5) driving MYC (and MAF). These findings lead us to postulate that in IMiDs resistant cells retention of BRD4 and MED1 at oncogenic enhancers in the absence of IKZF1 likely results from rewiring of the TFs regulating MYC. To identify TFs that may co-localize with BRD4 and IKZF1, we analyzed the enrichment of DNA motifs at IKZF1and BRD4 co-occupied enhancers using the MEME suite motif-finding algorithms. This computational analysis revealed a strong enrichment at these MM enhancers of the GGAA motif recognized by the ETS family of TFs (P= 3.2 e-743) and other motifs boxes for the RUNX (P= 9.6 e-725), MYC/MYB ( P= 8.8 e-52) and interferon regulatory (IRF) (P= 3.1 e-293) TFs. We next confirmed that the ETS family TF ETV4 was indeed expressed in IMiDs resistant, but not sensitive, MM cell lines. ChiPseq occupancy profiles in IMiDs resistant RPMI8226 cell line revealed co-localization of ETV4 with IKZF1, P300 and BRD4. As predicted, lenalidomide treatment induced global depletion of IKZF1 but not ETV4 at BRD4 occupied enhancers in resistant cell lines (RPMI8226 and ARP1). Importantly, Cas9-mediated knock out of ETV4 in RPMI8226 cells sensitized them to lenalidomide with ensuing MYC downregulation and cell death. Confirming its role in MM, ETV4 transcript was indeed detectable in primary patients' samples in the CoMMpass data repository (ETV4 FPKM >1.0 in 112/724) and its expression was associated with significantly reduced survival outcomes (HR 0.64; P=0.0008). Similarly, high expression (top quartiles) of RUNX2 or MYB, TFs with enriched motifs at IKZF1 co-occupied enhancer loci, was also associated with decreased survival. Of note RNAseq analysis of paired patient samples pre- and post-IMiDs treatment (n=14 pairs) revealed significant upregulation of ETV4 at the time of acquired IMiDs resistance (7/14). Lastly transcriptome analysis of 101 patients enrolled in the RD arm (lenalidomide and dexamethasone) of the POLLUX trial (NCT02076009) confirmed the reduced survival of patients with top quartiles expression of ETV4 as well as MYB and RUNX2 (Fig.1) Conclusion: Transcriptional plasticity with expression of extra-lineage TFs such as the ETS family member ETV4 sustains oncogenic enhancers in MM overcoming IKAROS and AIOLOS dependency and promoting IMiDs resistance. Figure 1 Figure 1. Disclosures Neri: Celgene: Consultancy, Honoraria, Research Funding; Janssen: Consultancy, Honoraria, Research Funding. Soong: Jannsen: Employment. Chiu: Janssen: Employment. Bahlis: Takeda: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees; Amgen: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding, Speakers Bureau; Janssen: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding, Speakers Bureau; Celgene: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding, Speakers Bureau.
Abstract 5622: scRNA-seq reveals restructuring of the immune microenvironment of multiple myeloma following treatment with an IAP antagonist
Cancer Research | September 4, 2020 | Volume: 80 | Issue: 16_Supplement | Pages: 5622-5622
Abstract Previously, the IAP antagonist LCL161 was shown to have robust anti-myeloma activity in a transgenic mouse model where an acute inflammatory response accompanied a reduction in tumor burden. In order to fully examine alterations to the immune microenvironment, murine splenocytes from twenty mice transplanted with Vk12598 myeloma and administered single-agent treatment regimens of LCL161, cyclophosphamide, or bortezomib were harvested after 24 hours and subjected to single cell RNA-seq sequencing. After initial removal of erythrocytic and mitochondrial expressing cells, count matrices of over 50,000 total cells span twenty murine whole spleen samples. In an examination of the myeloid component, we identified 13 clusters of immune cells with unique compositions across samples due to differing treatment regimens and initial tumor burden. While minimally present in all seven non-transplanted mice, macrophages with elevated expression of complement activation system genes (C1q) were the most abundant myeloid cell population in each of our untreated, Vk12598 transplanted mice. Alternatively, populations of cDC2, pDC, moDC, and monocytes lacking expression of MHC class II markers were fixtures for healthy controls while markedly less abundant in tumor-bearing mice. Tumor burdens in three samples subjected to LCL161 therapy were significantly reduced 24-hours post treatment and accompanied by an expansion of late-stage neutrophils and migratory DCs with expression features mirroring those of DCs from mesenteric lymph nodes. This composition of immune cells was distinct from all other sample including those treated with cyclophosphamide or bortezomib. In addition, relative proportions of B, T, and NK cells were nearly fully reconstituted, matching levels observed in healthy controls, in two of these samples. In a direct functional assessment, expression levels of Fscn1 and Ccr7 (key over-expressed genes in these migratory DCs) markedly increased following treatment of bone marrow derived DCs with LCL161. Overall, the transcriptional profile of these migratory DCs indicates their pivotal role in maturation of DCs (Fscn1, Tmem123) as well as the recruitment of T and NKT cells (Cxcl16, Ccr7) and increase in antigen presentation (Cd63)–all of which may serve to ‘warm' tumors, increasing T cell infiltration necessary for immunotherapy approaches. With the increased resolution of scRNA-seq, we observe unique alterations in the composition of the immune microenvironment of MM induced by treatment with LCL161 that further detail its unique mechanism of anti-MM activity. Citation Format: Caleb K. Stein, Erin Meermeier, Seth Welsh, Victoria Garbitt, Meaghen Sharik, P. Leif Bergsagel, Nicholas Banovich, Jonathan Keats, Marta Chesi. scRNA-seq reveals restructuring of the immune microenvironment of multiple myeloma following treatment with an IAP antagonist [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5622.