Braggio, Esteban

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.

Abstract Extramedullary disease (EMD) is associated with an inferior prognosis and lower response rates to conventional multiple myeloma (MM)–directed therapies compared to MM without EMD. A deeper understanding of the molecular landscape and underlying drivers of EMD is essential to identify potential targets for novel therapeutic strategies. To address this, we performed whole-exome sequencing on EMD tumor tissue from 18 unique patients and bone marrow aspirates (BMAs) from 20 patients at the time of EMD development. Notably, paired EMD and BMA samples were collected from 6 patients at the point of EMD diagnosis, allowing for direct comparison of molecular profiles. Our analysis revealed a near-universal presence of mutations within the MAPK pathway in EMD samples (94%), compared to BMAs (60%; odds ratio, 10.7; P = .02; q < 0.1). Additionally, mutations in established driver genes (NRAS, KRAS, and BRAF) were common and frequently clonal, suggesting their central role in EMD pathogenesis. We also identified alterations in genes associated with cell adhesion and migration (ROBO1, ROBO2, and FAT1) and the SWI/SNF complex and epigenetic regulators (ARID1A, KMT2C, KMT2D, and EP300), although these were predominantly subclonal. Furthermore, we frequently detected biallelic alterations in the tumor suppressor genes MAX (22%), a binding partner for MYC, and CDKN2C (17%). Genomic complexity was significantly higher in EMD samples than BMAs, as evidenced by increased tumor mutational burden and the enrichment of 1q gain/amplifications. These findings underscore the distinct molecular profile of EMD compared to BMA and highlight the genomically complex and heterogeneous nature of extramedullary disease in MM.

Abstract Most disease-associated single-nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWASs) reside in noncoding regions. There is increasing evidence that they play regulatory roles through epigenetic mechanisms by altering local chromatin state, 3D chromatin structure, and looping interactions. We and others have identified 42 CLL risk loci from GWASs. At 12 risk loci located outside of gene promoters, we previously demonstrated that functional SNPs are highly enriched in enhancers and super-enhancers that mostly target BCL2 family genes and transcriptional regulators. However, the study was limited by using only public epigenomic data generated in lymphoblastoid cell lines (LCL) and normal immune related cell types. To gain a better understanding of how the functional SNPs alter local chromatin state and modulate target gene expression in CLL, we generated RNA-seq data (n=588) and transposase-accessible chromatin with sequencing (ATAC-seq) data (n=69) in CLL tumor B cells. We performed molecular quantitative trait loci (QTL) analysis with index and correlated SNPs (R2 0.5) at 41 risk loci (excluding the HLA locus). The analyses identified cis expression QTLs (eQTLs), within 1Mb of genes’ transcription start site, at 25 risk loci and chromatin accessibility QTLs (caQTLs, within peak) at 11 risk loci, with both eQTLs and caQTLs at 8 risk loci. We identified a total of 54 eGenes, for which the expression is affected by nearby SNPs within 1Mb, including 40 novel eGenes at 16 CLL risk loci not previously reported. Notably, based on publicly available chromosome interaction data from CLL, LCL, and normal B cells, for over 70% of the novel eGenes, their promoters showed looping interactions with the regions carrying the eQTLs. Further, 13 SNPs were identified as both caQTL and eQTL, with the associated novel eGenes including CD81 at 11p15.5, IPCEF1 at 6q25.2, as well as CXXC1 and MBD1 at 18q21.1. The result suggests that these SNPs play a role in altering chromatin accessibility contributed to gene expression variation. At 11p15.5 risk locus, SNP rs2651823 (A/G) represents a novel caQTL and eQTL. The alternative G allele is associated with increased chromatin accessibility and CD81 expression. Hi-C data from GM12878 and CD20+ B cells supported interaction between the open chromatin region carrying rs2651823 and CD81 promoter. CD81 facilitates the trafficking of CD19 to the surface of activated B cells, and CD81-null mice showed impaired immune response. In a mouse model, CD81 contributes to tumor progression and metastasis in Lewis lung carcinoma. In summary, our molecular QTL analysis revealed that functional SNPs, located in CLL risk loci, modulate chromatin accessibility and expression of downstream target genes that could impact CLL risk. Citation Format: Huihuang Yan, Zhiquan Wang, Mingma Hoel, Cristine Allmer, Chantal McCabe, Daniel O’Brien, Dennis Robinson, James Cerhan, Sameer Parikh, Esteban Braggio, Susan Slager. Molecular QTL analyses identify functional variants at inherited risk loci in chronic lymphocytic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 2763.

Abstract Despite the advent of novel targeted therapies with high overall response rates, chronic lymphocytic leukemia (CLL) is still largely incurable due to the development of resistance. The bone marrow microenvironment (BMME) plays an important role in the survival and drug resistance in CLL. BM stromal cells (BMSCs) are considered a key component of BMME. Our long-term goal is to gain further insights of BMSCs-CLL crosstalk influencing both survival and drug resistance in CLL. To understand the global biologic changes that occur in leukemic CLL cells when in contact with BMSCs, we assessed RNA-seq in paired CLL cells isolated from untreated blood and BM and in untreated patients’ CLL cells cultured alone or with BMSCs. CLL cells from blood/BM and cultured alone/with BMSCs were also examined by Immunoprecipitation (IP), Western blot (WB), untargeted metabolomics (LCMS+GCMS), and preclinical drug sensitivity assays. We detected only 13 upregulated genes that overlap between CLL cells from BM and cocultured with BMSCs (p&lt0.05, fold change &gt1.5). When we analyzed the expression of these 13 genes in CLL cells from a cohort of 162 untreated patients by RNA-seq, we found associations between 4 genes (PNP, C16orf54, MOB3A, CDK2AP2) and clinical outcomes including overall survival (OS), progression free survival (PFS), and time to first treatment (TTFT) [multivariable Cox proportional hazards models, p&lt0.05]. Of the 4 genes, purine nucleoside phosphorylase (PNP), an enzyme in the purine pathway, had the most significant association with OS, PFS, and TTFT. Metabolomic profiling indicated an increased level of purine pathway metabolites adenosine, adenine, inosine, and hypoxanthine in cocultured CLL cells. WB analysis showed variable baseline PNP protein expression with 2 cohorts: low /no (L) and high (H) PNP using CLL blood cells. An induction in PNP protein levels were noticed in CLL cells expressing only LPNP when cocultured. Given the association of HPNP with adverse clinical outcomes, we tested for L vs. HPNP cohorts drug sensitivity to known BTK inhibitors(is) or Bcl2is. Treatment of CLL cells from untreated patients with Bcl2is venetoclax, S55746, and LP-118 (Bcl2-Bclxli, Newave) showed increased sensitivity to the Bcl2/Bclxlis in HPNP but not in LPNP cells. In addition, BTKis showed no difference in drug sensitivity associated with CLL PNP levels. To explain this discrepancy, WB and IP analyses were conducted in L vs. HPNP leukemic cells for apoptotic profiles, where we found increased phosphorylated-Bcl2(S70) and an increased association of Bcl2 with Bax in LPNP-expressing CLL cells compared to HPNP. These findings support increased resistance to apoptosis afforded by inhibiting Bcl2 for LPNP vs. HPNP-expressing CLL cells. Our future goals are to continue to understand the prognostic significance of the 4 genes with adverse clinical outcomes and to investigate how best to exploit the peculiar, enhanced sensitivity to Bcl2/Bclxlis found in the HPNP CLL cells. Citation Format: Sutapa Sinha, Weiguo Han, Zhiquan Wang, Brianna M. Finke, Heather C. Darby, Kari G. Rabe, Susan L. Slager, Chantal E. McCabe, Daniel R. O'Brien, Sameer A. Parikh, Esteban Braggio, Yi Chen, Fenlai Tan, Stephen P. Anthony, Yu Chen, Bing Dai, Yue Shen, Neil E. Kay. CLL drug resistance- A multifaceted mechanism [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 4423.

Abstract Monoclonal gammopathy of undetermined significance (MGUS) is a precursor to multiple myeloma (MM). Clonal hematopoiesis of indeterminate potential (CHIP) is a precursor to either myeloid or lymphoid malignancy, depending on whether the mutation is in a gene associated with lymphoid (L-CHIP) or myeloid (M-CHIP) outcome. Little is known about the association between CHIP and MGUS and their combined effect on risk of MM. We examined this question within Mayo Clinic Biobank cohort of participants aged 50 years or older with no prior history of hematological malignancies. Using whole blood, individuals were screened for MGUS using a matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry assay. M-CHIP and L-CHIP were defined based on mutations in 56 genes associated with myeloid malignancies or mutations in 235 genes associated with lymphoid malignancies, respectively, using whole-exome-sequencing. Incident MM was identified using International Classification of Diseases (ICD) codes and confirmed via medical record review. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals for the association of CHIP and MGUS. Cox proportional hazards regression was used to estimate hazard ratios (HR), with time defined as date between blood sample and the first incident MM, date of death, or 12/31/2023. All models were adjusted for age and sex. Of the 16, 479 individuals (median age 65 years, 45% male) screened, 12.3% were positive for MGUS, 9.8% were positive for M-CHIP, and 1.5% were positive for L-CHIP. There was no evidence of an association between M-CHIP and MGUS (OR=0.97, 95% CI:0.83-1.12) nor L-CHIP and MGUS (OR=1.12, 95% CI:0.78-1.57), including with any MGUS isotype (IgA, IgG, or IgM; P-values>0.22). Next, we investigated the association of these precursors on risk of MM. Median follow-up was 11.7 years (range 0.1-15.1), and 46 individuals developed incident MM. Both MGUS (HR=29.7, 95% CI: 14.1-62.4) and L-CHIP (HR=4.1, 95% CI: 1.3-13.2) were significantly associated with incident MM, but not M-CHIP (HR-1.0, 95% CI: 0.4-2.5). Investigating the combined effect of MGUS and L-CHIP, we found those with MGUS alone had a 27.5-fold (95% CI: 13.0-58.2) increased risk of incident MM compared to those with no precursors. No individuals with L-CHIP only developed MM, but those with both L-CHIP and MGUS had a 129.4-fold (95% CI: 34.4-486.4) increased risk of incident MM. In the largest cohort with screened MGUS and CHIP, we found no evidence of an association between either definition of CHIP with MGUS. However, individuals with both L-CHIP and MGUS were at the highest risk of developing MM. Citation Format: Nicholas J. Boddicker, Cristine Allmer, Danelle H. Moonen, Aaron D. Norman, Angelica Macauda, Alyssa Clay-Gilmour, Janet E. Olson, Mrinal S. Patnaik, Vincent Rajkumar, Esteban Braggio, David L. Murray, Shaji Kumar, Susan L. Slager, Celine M. Vachon. The association between MGUS and CHIP and their impact on multiple myeloma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 4917.

Genome-wide association studies (GWAS) of multiple myeloma (MM) in individuals of genetic European ancestry (EA) have identified 35 susceptibility loci. Co-heritability analyses have revealed a strong genetic correlation between MM and its precursor condition, monoclonal gammopathy of undetermined significance (MGUS). Previous research validated the association between a polygenic risk score (PRS) based on 23 MM risk loci and the risk of MGUS. Went and colleagues recently reported 12 new MM risk loci. We constructed a MM-PRS using all 35 MM risk variants and assessed its association with MGUS. Additionally, among individuals with MGUS, we assessed associations of MM-PRS with established prognostic factors [M-protein (<0.2 g/dL; 0.2-1.5 g/dL; >1.5 g/dL); isotype (IgG, IgA, IgM); free light chain (FLC) ratio (Normal (0.26-1.65), Abnormal (<0.26 or >1.65))]. Our study population included 1,723 individuals with MGUS seen at Mayo Clinic, Rochester, MN and 14,175 controls from the Mayo Clinic Biobank who lived in the 27-county region around Rochester, MN, and screened negative for MGUS using MALDI-TOF MS. The majority of MGUS patients and controls were of EA. We used logistic or multinomial regression to calculate odds ratios (ORs) and 95% confidence intervals (CIs). All models were adjusted for age, sex, and the first three principal components. The MM-PRS was significantly associated with MGUS when assessed continuously (OR, 1.28 per standard deviation (SD); 95% CI, 1.22-1.35). When comparing to the middle quintile of the MM-PRS distribution, individuals in the highest quintile had 1.41-fold increased risk of MGUS (95% CI, 1.21-1.64), and individuals in the lowest quintile had 0.68 decreased risk of MGUS (95% CI, 0.57-0.81). Of the 35 SNPs evaluated, 27 demonstrated consistent directions of effect for risk of MGUS with published risk estimates for MM. The effect of MM-PRS differed by isotype (P = 2.2 x 105), with higher MM-PRS associated with IgA MGUS (OR, 1.40 per SD; 95% CI, 1.08-1.81) and a lower MM-PRS associated with IgM (OR, 0.66 per SD; 95% CI, 0.52-0.84) compared to IgG MGUS. The MM-PRS was higher in MGUS patients with M-protein > 1.5 g/dL compared to M-protein of size 0.2-1.5 g/dL, although not statistically significant (OR, 1.41 per SD; 95% CI, 0.95-2.08); there were no differences in the MM-PRS between MGUS patients with M-protein size <0.2 g/dL compared to 0.2-1.5 g/dL (OR, 1.01 per SD;95% CI, 0.81-1.26). There were also no differences in MM-PRS by FLC ratio, (OR, 0.95 per SD; 95% CI, 0.852-1.10 for abnormal vs. normal). Our findings of an association of the MM-PRS with risk of MGUS provides further evidence for shared heritability between MGUS and MM. Future studies are needed to examine whether the expanded MM-PRS adds to existing models for MGUS progression. Citation Format: Angelica Macauda, Alyssa Clay-Gilmour, Cristine Allmer, Danelle Moonen, Aaron D. Normann, Nicholas J. Boddicker, Joselle M. Cook, Linda Baughn, Elizabeth E. Brown, Vincent S. Rajkumar, Esteban Braggio, David L. Murray, Susan L. Slager, Shaji Kumar, Celine M. Vachon. Implications of multiple myeloma polygenic risk scores (PRS) for MGUS [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 2274.

Abstract Introduction: Monitoring minimal residual disease (MRD) is critical in multiple myeloma (MM) to predict outcomes and guide therapy. Traditional bone marrow (BM) aspirates for MRD detection are invasive and limited by sample quality. We therefore explored cell-free DNA (cfDNA) fragmentation as a less invasive alternative for MRD detection. Methods: We performed 30-40X whole-genome sequencing (WGS) on peripheral blood cfDNA from 45 MM patients, collecting baseline (n = 45) and follow-up samples (n = 98) from eight Canadian sites (M4 and IMMAGINE studies) and one U.S. site (SPORE study), plus 25 healthy controls. Samples were collected post-induction therapy (n = 13), 100 days post-autologous stem cell transplantation (ACST, n = 36), and during maintenance therapy (n = 49). Multiparameter flow cytometry (MFC) MRD testing was conducted on 77 samples. We evaluated insert size metrics and MM-specific chromatin accessibility (from Ordoñez et al., 2020) using Griffin (described in Doebley et al., 2022). Results: Baseline samples had a higher proportion of short fragments (20-150 bp) than follow-up samples (p = 0.0037) and healthy controls (p = 0.019). A fragment score (FS) based on weighted fragment distribution (per Vessies et al., 2022) was highest at baseline (mean = -0.384), followed by MRD-positive (mean = -0.495) and MRD-negative (mean = -0.550) samples. MRD-negative FS was lower than baseline samples (p < 0.001) but higher than healthy controls (mean = -0.639; p = 0.016). Baseline samples had lower coverage at MM-specific chromatin regions than MRD timepoints (p < 0.05; baseline mean = 0.982; MRD-positive = 0.989; MRD-negative = 0.994; healthy = 0.993), suggesting higher transcriptional activity pre-treatment. MRD-negative samples had higher coverage than MRD-positive cases (p < 0.01). Using logistic regression, we assessed the predictive performance of FS and coverage at MM-specific sites. FS alone achieved an AUC of 0.632 (sensitivity [SN] 39.4%, specificity [SP] 92.1%, accuracy [AC] 67.6%). Coverage yielded an AUC of 0.681 (SN 75.8%, SP 60.5%, AC 67.6%). Combining FS and coverage at MM sites improved performance (AUC = 0.734; SN 48.5%, SP 92.1%, AC 71.8%). We next evaluated whether adjusting the BM tumor cell percentage cutoff by MFC could improve performance. At 0.017% (1 in 5, 789 cells) for MRD positivity, the combined model’s accuracy reached 83% with an AUC of 0.757 (SN 17%, SP 97%). Conclusions: cfDNA fragmentomic analysis achieved high specificity (92-97%) in detecting MRD, providing a less invasive alternative to facilitate serial monitoring. While sensitivity remains an area for improvement, the high specificity of cfDNA makes it valuable for confirming MRD negativity and reducing invasive procedures. Future efforts aim to enhance sensitivity by integrating additional fragmentomic features. Citation Format: Dor D. Abelman, Jenna Eagles, Aimee Wong, Saumil Shah, Stephanie Pedersen, Stephenie Prokopec, David S. Scott, Sarah Bridges, Darrell White, Irwindeep Sandhu, Kevin Song, Esteban Braggio, Alli Murugesan, Anthony Reiman, Suzanne Trudel, Trevor J. Pugh. Tumor-independent monitoring of minimal residual disease in multiple myeloma using cfDNA fragmentomics [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 3250.

Introduction The major primary genetic events in multiple myeloma (MM) include translocations of immunoglobulin heavy chain (IgH) and a hyperdiploidy (HRD) karyotype (Kuehl, 2002). The main drivers of pathogenesis in HRD MM are unknown (Fonseca, 2003), therefore further research is needed in this subgroup. Translational research in HRD MM has relied on analysis of genomic and transcriptomic alterations, however changes at the genome and transcriptome level do not always correlate directionally with protein abundance. Because proteins are considered the main and final effectors of most cellular processes, proteogenomics can improve our understanding of the drivers of HRD MM. We employed it here to identify potential biomarkers in HRD MM and determine whether any given trisomy may lead to a selective advantage. Methods: Protein and mRNA were extracted from CD138+ cells from 47 HRD MM and 14 t(11;14) MM samples. RNA-seq data was processed using MAPR-Seq pipeline using default parameters and gene-wise raw counts were extracted for each sample. Proteomics data was processed using Specranaut and protein-wise intensities were extracted for each sample. A generalized linear model configured with negative binomial distribution was used to analyze gene-wise count data in edgeR software. Genes with an adjusted p-value of <=0.05 were considered significantly differentially expressed. Protein intensity data were normalized using trimmed mean of M-values and log2 transformed. A generalized linear model configured with Gaussian distribution was used to compare protein intensities between the groups. Proteins with an adjusted p-value <= 0.05 were considered significantly differentially expressed. Two separate rank lists for protein and mRNA data were generated for pathway analysis. Results When comparing our small cohort of HRD MM patient samples to those with t(11;14) MM, out of approximately 16,600 total transcripts and 9,700 proteins evaluated we identified approximately 2,100 significantly variant transcripts and 500 significantly variant proteins in HRD MM. GSEA pathway analysis of the transcriptome and proteome revealed significantly upregulated translational and ribosomal processing pathways in HRD MM. Transcriptional MYC target pathways were also significantly upregulated. RNAseq revealed significantly higher expression of MYC (Rank 24.5, p=5.3E-6) in HRD MM. It is known that patients with t(11;14) MM have increased expression of CCND1 (Kuehl and Bergsagel, 2012). As expected, we found that the HRD MM subgroup had significantly lower expression of CCND1 compared to the t(11;14) MM subgroup (Rank -79.9, p=2.8E-21) and this was also true at the protein level (Rank -21.7, p=2.2E-5). The most significantly abundant protein in the HRD MM subgroup compared to t(11;14) MM was BAHD1, which is a protein that targets H3K27me3 and recruits transcriptional co-repressors (Fan, 2020). BAHD1 has been shown to repress erythroid genes and may recruit HDAC5, which plays a role in osteoblast differentiation (Rank 72.3, p=2.9E-20; Xu, 2021). The second most abundant protein in HRD MM was CKS2, which binds to the catalytic subunit of cyclin-dependent kinases (CDKs) and therefore positively regulates cell proliferation, invasion, and migration and has been shown to be implicated in tumor progression in several solid tumors (Rank 72.3, p=2.9E-20; You, 2015). TRIM proteins promote tumorigenesis and cancer development. TRIM14 mediates cancer development through mediating the JAK/STAT, PI3K/AKT, NFkB, and p53 pathways and was found to be significantly higher in HRD MM (Rank 70.6, p=9.0E-20; Huang, 2022). We also found increased protein abundance of LAMP5 in HRD MM (Rank 20.1, p=6.0E-5). LAMP5 may play an important role in MM progression and has been shown to be highly expressed in patients with mixed-lineage leukemia and is associated with a poor prognosis (Wang, 2023; Gracia-Maldonado, 2022). We identified many proteins (not all discussed here) that have not been well described in HRD MM and thus may be biomarkers/targets of disease pathogenesis of this subgroup. Conclusion Proteins are the final effectors of most cellular processes. Proteogenomics can improve our understanding of the drivers of HRD MM. Using a proteogenomic approach we identified proteins that may lead to selective advantage and/or pathogenesis of trisomies in HRD MM.