Diamond, Benjamin

Abstract Studies have reported conflicting findings regarding the contribution of germline variants or somatic genomic drivers to racial disparities in multiple myeloma. To comprehensively investigate somatic drivers in relation to inherited genetics in multiple myeloma, we combined newly sequenced whole-genome sequencing data with publicly available datasets (total n = 1,286). Overall, we did not identify germline or somatic genomic differences that explain the different risk of developing multiple myeloma between patients with genetic similarity to African (AFR) or European (EUR) reference populations. A difference in the detectability and timing of APOBEC-associated and germinal center mutational activity was observed. Integrating epidemiologic data and mutational signature-based temporal estimates, we challenge the assumption that individuals in the AFR group develop multiple myeloma at a younger age. Finally, we demonstrate that, with equal access to efficacious therapies, patients in the AFR and EUR groups have equivalent clinical outcomes. Significance: Multiple myeloma is reported to occur at higher rates in individuals who self-identify as non-Hispanic Black. In this large dataset, genomic drivers occur at the same rate among ancestry groups, except for APOBEC mutagenesis. With equivalent therapy, clinical outcomes did not differ for patients grouped by genetic ancestry similarity.

Abstract Multiple myeloma (MM) with t(11;14)(CCND1;IGH) remains the only subset sensitive to the BCL2 inhibitor venetoclax. However, not all patients with t(11;14)(CCND1;IGH) respond to treatment, and some progress early after initial response. To investigate this, we examined 44 whole-genome and whole-exome sequencing data samples from 34 patients with t(11;14) MM treated with venetoclax. The presence of mutations in the RAS pathway was strongly associated with shortened progression-free survival (PFS) and was validated in an independent cohort of 21 patients with MM. The presence of 1q gain was also associated with shorter PFS in patients without RAS mutations. In 10 patients with paired prevenetoclax and postvenetoclax treatment samples, postvenetoclax progression was recurrently driven by the selection of genomic events in the BCL2/MCL1 and RAS pathways and of high-risk features (eg, loss of TP53 and CDKN2C). Overall, our study shows that comprehensive genomic profiling can identify most mechanisms underlying resistance to BCL2 inhibition in t(11;14)(CCND1;IGH) MM.

ABSTRACT Chimeric antigen receptor T-cell (CART) and T-cell engager (TCE) therapies targeting B-cell maturation antigen (BCMA) are transforming the treatment landscape for relapsed multiple myeloma (MM). However, despite impressive initial response rates, most patients eventually relapse. To investigate this unmet medical need, we applied whole-genome sequencing (WGS) to MM cells from cohorts of 102 relapsed patients treated with anti-BCMA CART and TCE therapies. Several genomic alterations were associated with clinical outcomes, particularly primary refractoriness, including high genomic complexity and mutations in genes regulating plasma cell identity, which predicted resistance to therapy. Single-cell RNA sequencing further revealed that MM cells from refractory patients exhibited high proliferation signatures and reduced expression of TNFRSF17 (encoding BCMA), while were less enriched for plasma cell–associated transcriptional programs, a phenomenon we term “plasma cell identity escape.” This profile was strongly associated with immune dysregulation of CD8 T cells including increased activation and exhaustion. This evolution of MM toward a more proliferative and lineage-divergent state, refractory to the anti-BCMA T-cell redirecting therapies, was functionally validated in preclinical MM mouse models. Collectively, our results comprehensively define the cellular and molecular mechanisms underlying primary resistance to anti-BCMA therapies.

Abstract Background: Multiple myeloma (MM) is always preceded by monoclonal gammopathy of undetermined significance (MGUS) and smoldering multiple myeloma (SMM). While these precursor conditions are asymptomatic, they are not entirely benign and carry a lifelong risk of progression to MM. Unlike other cancers defined by pathology, malignant transformation from MGUS or SMM to MM has so far relied on demonstration of clinical end-organ damage or high tumor cell burden (i.e. SLiM criteria, SV Rajkumar Lancet Onc 2014), as morphology and cytogenetics cannot reliably distinguish them. Methods: To identify genomic features underlying clinical behavior in MGUS and SMM, a total of 374 patients with SMM (n= 290) or MGUS (n=84) with available whole exome (whole exome sequencing, WES, n=190) or whole genome sequencing (whole genome sequencing WGS, n=184) were included in the study. Overall, none of the patients had evidence of SLiM or CRAB criteria at the time of the biopsy. Based on the center of origin and availability of WES and WGS, 227 and 97 patients were used to form our study training and validation set. The study was designed in two parts: 1) defining a new classification for MM precursor conditions; 2) integrating genomics to existing prognostic score to improve accuracy in identifying patients at imminent risk of progression. Overall, 68 patients in the training set were either lost to follow-up or enrolled in a clinical trial and they were only included in the classification part of the study. Results: By leveraging the distribution of established MM-defining genomic events (Maura et al. JCO 2024) in the training set, we developed a workflow that differentiates MGUS and SMM into two genomically distinct entities: one with evidence of neoplastic transformation (i.e., genomic MM) and one without (i.e., genomic MGUS). Overall, 39% of MGUS and 91% of SMM cases had genomic evidence of neoplastic transformation in the training set (i.e., genomic MM). When we applied the same workflow and classification to the validation set, 46% and 83% of MGUS and SMM had evidence of neoplastic transformation, respectively. Importantly, none of the patients classified as genomic-MGUS progressed to MM in either the training or validation cohorts, with a median follow-up of 40 and 61 months, respectively. Moreover, all patients with IMWG 2/20/20 high risk SMM were classified as genomic MM in both the training and validation set. Interestingly, 5 out of 62 (8%) patients enrolled in early intervention trials for high risk SMM were genomically classified as genomic-MGUS. Overall, these data demonstrated that the MM genomic background and neoplastic transformation can be acquired very early in time, even before the SMM phase, in line with what observed in solid tumors. To avoid overestimating the prognostic impact of myeloma-defining genomic events—which are, by definition, enriched in the genomic-MM group compared to genomic-MGUS—we restricted our prognostic analysis to patients in whom transformation has already occurred (genomic-MM). In fact, distinct genomic drivers may predict malignant transformation but not necessarily the risk of imminent progression to MM once transformation has occurred. After running 5-fold cross validations corrected for IMWG 2/20/20, only presence of RAS mutations, genomic events involving MYC, presence of neoplastic copy number variant (CNV) signatures and APOBEC mutagenesis were predictive of early progression in the training set. Adding the presence of high-risk myeloma-defining genomic events to the IMWG 2/20/20 model (i.e., genomic-IMWG) significantly improves the accuracy of predicting MM progression in both the training (c-index 0.69 vs. 0.74) and validation sets (c-index 0.74 vs. 0.79), compared to the IMWG 2/20/20 alone. Conclusion: By leveraging comprehensive genomic profiling integrated with the clinical IMWG 2/20/20 risk score, we propose a novel classification and prognostication system for MM precursor conditions, which will better define patients who may benefit from early therapeutic intervention versus watch-and-wait monitoring.

Multiple myeloma (MM) is consistently preceded by monoclonal gammopathy of undetermined significance (MGUS) and smoldering multiple myeloma (SMM). While these precursor conditions are asymptomatic, they are not entirely benign and carry a lifelong risk of progression to MM. Unlike other cancers defined by pathology, malignant transformation from MGUS or SMM to MM has so far relied on demonstration of clinical end-organ damage as morphology and cytogenetics cannot reliably distinguish them. In this study, using genomic data from 374 patients with MGUS or SMM (277 training, 97 validation), to our knowledge, we demonstrate for the first time the ability to identify malignant transformation in MGUS and SMM. We introduce the concept of genomic MM and genomic MGUS to differentiate the subsets of MGUS and SMM that are biologically malignant with genomic features indistinguishable from MM from the subset that is premalignant and unlikely to progress to malignancy. Importantly, we find that most SMM has biological features of malignant transformation indistinguishable from MM. As expected, this subset that we consider having genomic MM is associated with a high risk of progression to MM although some patients remained progression-free beyond 5 years. Conversely, 60% of MGUS and 10% of SMM have no evidence of malignant transformation (genomic MGUS), with no progression during follow-up. Integration of genomic features with the 2/20/20 International Myeloma Working Group model significantly improved the prediction of progression among genomic MM. These findings support the use of genomic criteria to refine the classification and the risk stratification in myeloma precursor conditions.

Abstract Monoclonal antibodies targeting CD38 are a therapeutic mainstay in multiple myeloma (MM). Although they have contributed to improved outcomes, most patients still experience disease relapse, and little is known about tumor-intrinsic mechanisms of resistance to these drugs. Antigen escape has been implicated as a mechanism of tumor-cell evasion in immunotherapy. Yet, it is unknown whether MM cells can develop permanent resistance to anti-CD38 antibodies by acquiring genomic events leading to biallelic disruption of the CD38 gene locus. Here, we analyzed whole-genome and whole-exome sequencing data from patients 701 newly diagnosed MM, 67 patients at relapse with naivety to anti-CD38 antibodies, and 50 patients collected at relapse after anti-CD38 antibodies. We report a loss of CD38 in 10 of 50 patients (20%) after CD38 therapy, 3 of whom exhibited a loss of both copies. Two of these cases showed convergent evolution in which distinct subclones independently acquired similar advantageous variants. Functional studies on missense mutations involved in biallelic CD38 events revealed that 2 variants, L153H and C275Y, decreased binding affinity and antibody-dependent cellular cytotoxicity of the commercial antibodies daratumumab and isatuximab. However, a third mutation, R140G, conferred selective resistance to daratumumab, while retaining sensitivity to isatuximab. Clinically, patients with MM are often rechallenged with CD38 antibodies after disease progression and these data suggest that next-generation sequencing may play a role in subsequent treatment selection for a subset of patients.

Introduction Monoclonal antibodies (MoAbs) directed against CD38 are a mainstay in multiple myeloma (MM). Downregulation of CD38 expression is a resistance mechanism of antigenic escape, however, cell surface expression may return. While this tumor cell plasticity is likely driven by epigenetic mechanisms, it is unknown if MM cells can develop permanent resistance to anti-CD38 MoAb by acquiring genomic events leading to biallelic disruption of CD38 thus precluding re-treatment. Here, we report biallelic loss of CD38 as a recurrent mechanism of resistance, often via convergent evolution where distinct subclones are advantaged towards the same selective pressure. Methods We interrogated two sets of patients (Mayo Clinic; n = 31, Calgary University; n =19) with relapse following anti-CD38 MoAbs. Two additional cases of interest were selected from Heidelberg University. Whole genome sequencing (WGS, 60-100x) or whole exome sequencing (WXS), flow cytometry (flow), and bulk RNAseq were used to validate the effect of genomic events on CD38. Missense mutations predicted to affect CD38 MoAb binding were validated with site-directed mutagenesis of a CD38 plasmid, transduction into K562 MM cells, and CD38 MoAb binding assays. Results The prevalence of biallelic inactivation of CD38 at anti-CD38 MoAb relapse was 6% (3/50). An additional selected biallelic event was noted in the Heidelberg set. Each case is described in detail: MM-19 (Calgary) received daratumumab (Dara), lenalidomide, and dexamethasone (Dex), and then had a second exposure to Dara, following relapse. At this second relapse, WGS was performed on two independent extramedullary disease localizations in the liver. Both samples had unique, large deletions (>3 Mb) on chromosome 4 coupled with additional unique focal deletions involving CD38. Flow confirmed the presence of CD138+/CD38- plasma cell populations in both sites. A second case of convergent evolution was seen in a patient (Mayo) after seventh-line Dara/Len/Dex. WGS revealed a common deletion encompassing CD38. On the remaining allele, one subclone had a structural variation-mediated loss of CD38(cancer cell fraction; CCF 40%) while another harbored an L153H missense mutation (CCF 60%). Flow confirmed presence of two CD138+ clones; one with diminished CD38 expression. L153H was functionally validated to abrogate Dara binding, while maintaining CD38 expression. The third sample was collected after fifth-line Dara and sixth-line Dara/Bortezomib/Dex (Heidelberg) and represents another case of convergent evolution with a common loss of chromosome 4 and two subclones detected by WGS (P98Lfs*12, CCF 39% resulting in biallelic inactivation; R140G, CCF 61%). Bulk RNAseq revealed, in comparison to NDMM from CoMMpass, CD38 expression below the first decile. Notably, the R140G variant reduced binding affinity to Dara but not to Isatuximab (Isa). A fourth case (Mayo) relapsing after 6 months of Dara-based therapy similarly had a focal deletion of CD38 with an L18Sfs*16 frameshift deletion on the remaining allele associated with a 4.35-fold reduction in CD38 expression by RNAseq compared to pre-treatment. Interestingly, in 3/4 cases the CD38-negative clone emerged after a second exposure to anti-CD38 MoAb, suggesting emergence with the first exposure and then selection during the second. One other monoallelic event (C275Y; Heidelberg) was detected post-Isa. Copy number of the alternate allele is not known given targeted sequencing methodology, but the variant was seen to reduce binding affinity to both Dara and Isa. The prevalence of monoallelic loss of CD38 in newly diagnosed MM from MMRF CoMMpass was 7% (50/701) and two patients had mutations in CD38, but there were no cases of biallelic loss suggesting these events are driven by treatment pressure. Conclusion Distinction of the mode of CD38 loss is valuable for predicting outcome of re-challenge to anti-CD38 MoAbs. We see here an estimated 5-10% of anti-CD38 MoAb relapses are characterized by genomic antigen escape. Results here are in line with antigen escape to BCMA- and GPRC5D-directed immunotherapies (Lee et al, Nat Med 2023), suggesting that acquired antigen escape resulting from prolonged exposure to the selective pressures of immunotherapy is a recurrent mechanism of resistance in MM. Moreover, in select cases, non-overlapping resistance mutations provide rationale for switching agents at second exposure.

Introduction The pathogenic role of somatic non-coding (nc) single nucleotide variants (SNVs) is underexplored in multiple myeloma (MM) and additional characterization of these regions has the potential to provide significant new insights into disease progression. To investigate the role of ncSNVs and their relationship with MM pathogenesis, we have characterized a novel series of WGS from 294 cases of MGUS, SMM, and MM using computational tools optimized to identify functional driver SNVs in these regions. We hypothesize the ncSNVs that co-localize to sites of super-enhancers (SEs) active in plasma cells (PC) are sites of high functional impact SNVs and these mutations lead to the dysregulated expression of genes that should be silent at the terminally differentiated stage of PCs, providing a selective advantage that shapes the fitness of the MM genomic landscape. Methods We analyzed 294 high-coverage matched tumor/normal pairs (80x tumor and 40x normal) derived from 21 MGUS, 86 SMM, 167 NDMM and 40 RRMM. In addition to previously sequenced samples, newly sequenced patients were obtained from NYU Langone Health, Memorial Sloan Kettering Cancer Center, University of Alabama at Birmingham and University of Miami following IRB-approved informed consent. DNA was extracted from CD138+ purified bone marrow plasma cells using standard approaches. We used the MGP1000 (https://github.com/pblaney/mgp1000) to process sequencing data which employs a union-consensus mechanism for determining the final set of variants. We collected enhancer elements (E/SEs) active in GC B-cells, ABC-DLBCL and GBC-DLBCL (Bal, et al. Nature, 2022) and PC specific E/SEs (Loven, et al. Cell 2013). FunSeq2 was used to score ncSNVs. GRAM was used to predict how variants alter TF binding and gene expression in a subset of 6 cases with paired WGS in our series and will expand this analysis on a validation set of 263 NDMM and 93 RRMM cases from University-Hospital Heidelberg enrolled in the GMMG-HD6 clinical trial. Results We identified 1,952,232 high-quality SNVs with 18,199 (0.93%) accumulating in 280 focal sites we termed “mutationally enriched regions” or MERs. To discern SNVs in MERs with the highest potential for imparting a selective advantage, we clustered mutations based on the density of FunSeq2 impact scores which uses reference data contexts on sequence conservation, modification of transcription factor (TF) binding, and gene linkage to regulatory elements and network hubs from pan-cancer studies. We show that MERs are enriched for ncSNVs with high impact scores (mean=2.57) compared to random similarly-sized regions of the genome (Fisher's exact test P<9.7e-12, OR=10.10 [3.88-37.40]). The high impact ncSNVs (hi-ncSNVs) in MERs (n=220) account for 8.09% of all hi-ncSNVs and are present in 46.3% of cases (136/294) with some patients having up to 7. The proportion of cases carrying a MER hi-ncSNVs was similar across disease stages (MGUS=42.9%, SMM=46.5%, MM=40.7%, and RR=47.5%). Analysis of the VAF of MER hi-ncSNVs shows a minor shift towards more clonal versus subclonal mutations (53.6% vs 46.4% VAF>=0.25) suggesting some high impact mutations are acquired in the GC but additional hi-ncSNV are acquired later in PC development, providing evidence that hi-ncSNVs accumulate in MERs throughout disease progression. We mapped these hi-ncSNVs to associated genes then intersected the loci with E/SEs active in GCBs and PCs to determine their functional context. This approach identified 12 genes affected by hi-ncSNVs including the previously identified nc-drivers NEAT1, MIR142, ZCCHC7, and IGLL5 as well as the PC differentiation genes IRF8, BACH2, proliferation genes RHOH, SPRED2, known coding MM drivers CCND1, ZFP36L1, and chr1q genes BTG2 and CROCC. We predicted the effect of ncSNVs on expression through modulation of TF binding at this set of putative nc drivers in MM and identified a trend towards increase in expression of SPRED2, CCND1, NEAT1, and IGLL5 compared to non-mutated cases (P=0.089). We will present data on the validation of this set of genes on a larger set of data derived from the Heidelberg group. Conclusions We identified a set of 12 putative nc driver genes in MM that modify expression and have the potential to significantly alter the pathogenesis of MM by providing a selective survival advantage. We will present the impact of these by clinical stage and on outcome.

PURPOSE Outcomes for patients with newly diagnosed multiple myeloma (NDMM) are heterogenous, with overall survival (OS) ranging from months to over 10 years. METHODS To decipher and predict the molecular and clinical heterogeneity of NDMM, we assembled a series of 1,933 patients with available clinical, genomic, and therapeutic data. RESULTS Leveraging a comprehensive catalog of genomic drivers, we identified 12 groups, expanding on previous gene expression–based molecular classifications. To build a model predicting individualized risk in NDMM (IRMMa), we integrated clinical, genomic, and treatment variables. To correct for time-dependent variables, including high-dose melphalan followed by autologous stem-cell transplantation (HDM-ASCT), and maintenance therapy, a multi-state model was designed. The IRMMa model accuracy was significantly higher than all comparator prognostic models, with a c-index for OS of 0.726, compared with International Staging System (ISS; 0.61), revised-ISS (0.572), and R2-ISS (0.625). Integral to model accuracy was 20 genomic features, including 1q21 gain/amp, del 1p, TP53 loss, NSD2 translocations, APOBEC mutational signatures, and copy-number signatures (reflecting the complex structural variant chromothripsis). IRMMa accuracy and superiority compared with other prognostic models were validated on 256 patients enrolled in the GMMG-HD6 (ClinicalTrials.gov identifier: NCT02495922 ) clinical trial. Individualized patient risks were significantly affected across the 12 genomic groups by different treatment strategies (ie, treatment variance), which was used to identify patients for whom HDM-ASCT is particularly effective versus patients for whom the impact is limited. CONCLUSION Integrating clinical, demographic, genomic, and therapeutic data, to our knowledge, we have developed the first individualized risk-prediction model enabling personally tailored therapeutic decisions for patients with NDMM.