Maura, Francesco

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.

Abstract Background Myeloproliferative neoplasms (MPNs) can progress from chronic phase to accelerated or blast phase (MPN-AP/BP), a transition with poor prognosis and limited treatment options. While some abnormalities have been described, the genomic and clonal drivers of transformation remain poorly defined. We characterized mutational and structural alterations in MPN-AP/BP, leveraging paired samples to trace clonal evolution. Methods We profiled 186 samples from 160 individuals, including 159 MPN-AP/BP and 27 chronic-phase MPN samples (26 MF, 1 PV), using targeted sequencing of a 588 gene panel (mean depth: 709x). Paired pre- and post-transformation samples were available for 26 patients (n= 52). A consensus pipeline was used to detect somatic mutations, copy number variations (CNVs), and copy-neutral loss of heterozygosity (cnLOH). Clinical cytogenetic data were integrated when available. Single-cell validation is ongoing. Results Genomic complexity was widespread: 99% of patients harbored ≥1 oncogenic mutation (median 4; range 1–10). In total, 859 mutations were identified, most commonly in JAK2 (68%), ASXL1 (28%), TET2 (26%), SRSF2 (24%), TP53 (22%), and RUNX1 (22%). 684 structural alterations were identified in 86% of patients (median 3; range 1-20), including arm-level (22.5%), segmental (>5 Mb; 60%), and focal (<2 Mb; 9%) alterations. Recurrent events included 9p cnLOH (16%), del(5q) (13%), del(17p) (10%), and focal deletion of 21q22.12 involving RUNX1 (10%). Complex karyotypes (≥3 abnormalities) were observed in 55% of patients. Biallelic mutations accounted for 21% (183/859) of mutations, affecting 58% of patients across 21 genes, most frequently JAK2 V617F (28%), TP53 (17%) and TET2 (14%). Patients with biallelic mutations in these genes had significantly higher blast percentages compared to those with monoallelic mutations (p<0.05), highlighting an association with aggressive disease. In 26 patients with longitudinal samples (chronic phase MPN and MPN-AP/BP), 87% (125/144) of mutations were already detectable during the chronic phase, on average more than 7 years prior to transformation. Only 11% (16/144) of mutations were newly acquired, including additional RUNX1 and TP53 hits, or structural changes that converted pre-existing JAK2 V617F and TP53 R248Q mutations from monoallelic to biallelic via 9p cnLOH or 17p loss, respectively. These findings suggest transformation is often driven by outgrowth of pre-existing high-risk clones, indicating that these clones likely have a competitive advantage and that the mutations within these clones have a cooperative biological interaction which underlies this competitive advantage. To investigate this hypothesis, we examined whether patterns of recurrent co-mutations were identifiable in this cohort. We identified recurrent co-occurring mutations in SRSF2-IDH2, ASXL1-SRSF2, and ASXL1-EZH2, which were enriched at transformation (p < 0.05; Fisher exact test). Importantly, presence of these co-mutations was consistently associated with a higher transformation risk in two independent cohorts of MF patients from Mayo Clinic (n= 405) and University of Florence (n= 518). Clonal architecture analysis revealed that these co-mutations frequently arose within the same malignant clone: 100% for SRSF2-IDH2, 80% for ASXL1-SRSF2, and 85% for ASXL1-EZH2, and are often within the dominant clone. Prior studies have demonstrated mechanistic interplay between SRSF2 and IDH2 (Yoshimi et al.), as well as ASXL1 and SRSF2 (Sui et al.) mutations in myeloid malignancies. We assessed functional synergy of ASXL1-EZH2 using a murine model and identified that dual deletion caused a rapidly fatal myeloid neoplasm with shorter survival than single-gene loss. These results reinforce that functional synergy emerges within a shared cellular context, supporting clonal co-dependence as a driver of progression. Single-cell RNA sequencing of ASXL1-EZH2 mutant cases is ongoing and will be presented. Conclusions Transformation to MPN-AP/BP is driven by gradual clonal remodeling and expansion of pre-existing high-risk abnormalities. Biallelic and multi-hit alterations in TP53, JAK2, TET2, and RUNX1 promote clonal dominance and blast progression. Recurrent co-mutations such as IDH2-SRSF2 and ASXL1-EZH2 arise within the same clone, exhibit functional synergy, and are linked to poor outcomes. These findings support the clinical utility of early clonal and co-mutational profiling to guide risk-adapted intervention.

Abstract Cyclin D1, encoded by CCND1, is a cell cycle regulator, overexpressed in ~15% of multiple myeloma (MM) cases via CCND1::IGH super-enhancer juxtaposition through t(11;14). Cyclin D1 also plays a non-canonical role in transcriptional regulation, but in t(11;14) MM its function in gene regulation remains undefined. We investigated the global transcriptional and epigenomic consequences of CCND1overexpression in t(11;14) MM. To define a CCND1 signature, RNA-seq from two independent primary MM patient cohorts, CoMMpass and Mayo, was utilized to compare high CCND1 t(11;14) cases (≥300 transcripts-per-million [TPM], n=271) vs. low CCND1 hyperdiploid cases (<30 TPM, n=233). Two CCND1 knockout (KO) clones in t(11;14) MM U266B1 were also generated by CRISPR/Cas9, introducing a frameshift insertion (A) after chr11:69641478 (hg38) in exon 1. ATAC-seq and ChIP-seq (H3K4me3, H3K4me1, and H3K27ac) were performed in duplicates, and RNA-seq in triplicates in U266B1 wild-type (WT) and KO. Libraries were sequenced on Illumina NextSeq 2000 (2x51 bp for ChIP/ATAC, 2x150 bp for RNA-seq). Reads were mapped to hg38 by BWA (ChIP/ATAC) or STAR (RNA-seq), peaks called by MACS2 (FDR< 0.01), peak and gene raw counts estimated by featureCounts. Differential peaks (fold-change ≥1.5) and differentially expressed genes (fold-change ≥2) were identified by edgeR (FDR<0.05), including peaks detected in ≥2 samples and genes with TPM ≥0.1 in ≥1 sample. Interactions between promoters and distal regulatory regions were identified using publicly available chromatin interaction data in blood cell types. Pathway analysis was performed using Enrichr. Co-immunoprecipitation was performed in t(11;14) KMS12BM cells. Differential analysis comparing CCND1-high t(11;14) and CCND1-low MM cases in the two patient cohorts identified 21 upregulated and 21 downregulated genes associated with CCND1 expression which were also captured in the wild-typevs. CCND1 KO U266 cells. Enriched pathways (p<0.05) included Notch signaling, early/late estrogen response, epithelial-mesenchymal transition and reduced mitotic spindle and myogenesis. Dysregulated genes included those related to apoptosis, ubiquitin-proteosome, extracellular matrix, adhesion, and immune response, demonstrating a broader role for cyclin D1 in MM beyond cell cycle control. To explore the epigenomic basis of this transcriptional regulation, for each of the 42 differentially expressed genes we examined the changes of histone marks and chromatin accessibility in promoter and the interacting distal regulatory regions in the U266 KO model. Ten of the 21 genes upregulated in CCND1-high (downloaded in KO; e.g., argininosuccinate synthase 1 [ASS1], PMAIP1, MYEOV, CCR7, and DTX1) and 12 of the 21 genes downregulated in CCND1-high (upregulated in KO; e.g., UACA, RASAL2, PTPRZ1, PPT2, and IL23R) showed concordant changes in histone marks and/or chromatin accessibility, predominantly in distal regulatory regions, suggesting a chromatin-mediated mechanism of gene regulation by cyclin D1. The most prominent changes involved H3K27ac, a mark defining active enhancers. As H3K27ac is deposited by histone acetyltransferases CBP, we tested for the possible cyclin D1-CBP interaction. Co-IP confirmed this interaction, suggesting a role for this complex in chromatin remodeling and transcriptional control in t(11;14) MM. One gene notably upregulated in association with CCND1 is ASS1, which encodes a urea cycle enzyme responsible for endogenous arginine synthesis. Metabolomic studies demonstrated KO cells had elevated aspartate and citrulline and lowered arginine and ornithine compared to WT cells (p<0.01) demonstrating impaired urea cycle function. Further, KO cells displayed significantly reduced ASS1 protein expression, cell growth and S-phase entry compared to WT cells in arginine-replete conditions, with both WT and KO cell types further inhibited by arginine depletion (p<0.05). Together, a cyclin D1-associated gene signature was identified supporting its role in t(11;14) MM beyond the CDK4/6 mediated cell cycle control, including a role in metabolic plasticity and survival under arginine depleted conditions. These findings highlight cyclin D1 as a therapeutic vulnerability and open new avenues for personalized treatment strategies.

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 Motivation Efforts to address health disparities are often limited by the lack of robust computational tools for inferring genetic ancestry by calculating an individual’s genetic similarity to continental groups. We have already shown that a preferred alternative to self-described race is using ancestry-informative markers (AIMs) that can be classified into ancestral components and used to estimate their similarity to those of known populations to identify continental groups. However, real-world genomic data can present challenges, including limited availability of germline DNA, a small number of AIMs for each sample, and the use of different variant calling software, limiting the application of existing solutions. Results Here, we describe a novel supervised machine-learning tool AncestryGeni, which infers genetic ancestry for samples with even a hundred markers and is applicable to any genomic data, including whole exome sequencing (WES) and RNA sequencing (RNA-Seq) data. Applying AncestryGeni to a real-world genomic dataset obtained from the Multiple Myeloma Research Foundation (MMRF) CoMMpass study, we show that it is more accurate than the commonly used FastNGSadmix when using nonstandard genomic material. We also demonstrate that when using AncestryGeni, the tumor-derived sequence obtained from WES and RNA-Seq can be a robust data source to accurately estimate an individual’s genetic similarity to a continental group. Availability and implementation AncestryGeni pipeline is available at https://github.com/eelhaik/AncestryGeni/tree/main.

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.

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).