Citations

Abstract Diffuse Large B Cell Lymphoma (DLBCL), the most common subtype of aggressive B-cell Non-Hodgkin Lymphoma (NHL), accounts for 30-40% of newly diagnosed NHL cases. Despite recent progress improving our understanding of molecular features of DLBCL, there is still a high relapse rate (~45%) among patients who undergo standard treatment (R-CHOP). Several groups have identified reproducible molecular clusters of DLBCL although, the biologic significance of many of the specific gene mutations within these clusters is poorly understood. DTX1 mutations are a feature of the BN2/C1 DLBCL classifiers and are among the mutations with unknown significance. BN2/C1 DLBCL are a mixture of ABC, GCB, and Unclassified cell of origin (COO) with a common signature being the presence of genetic alterations in genes known to regulate the NOTCH pathway; yet, these tumors lack transcriptional signatures of traditional NOTCH activation, suggesting a unique and undefined biology driven through non-canonical NOTCH pathway activation. DTX1 is an E3 ubiquitin ligase which directly binds the intracellular domain of the NOTCH receptor (NICD), targeting it for degradation, resulting in suppression of NOTCH pathway activity. We therefore hypothesize that DTX1 mutations lead to NOTCH pathway dysregulation and are a critical driver for a subset of BN2/C1 DLBCLs. Using whole exome sequencing (WES) data from 404 newly diagnosed DLBCL patients, we found 8% of patients with mutations in DTX1. Alterations in DTX1 were localized to the NICD binding domain, WWE1. In silico Alphafold modeling of DTX1 WWE1-interaction with NICD revealed that patient mutations clustered within predicted interacting surfaces with NICD. These findings lead to the hypothesis that mutations in the DTX1 WWE1 domain suppress its ability to bind to and degrade NICD, enhancing NOTCH pathway activation in DLBCL. To address this question, we first compared DTX1 protein expression by immunoblot in DTX1 WWE1 mutant (OCI-Ly3 and SUDHL2) and wildtype (WT) (OCI-Ly7, OCI-Ly19, SUDHL6, and U2932) DLBCL cell lines. We found higher expression of NICD1 (p=0.03) in OCI-Ly3 and NICD2 (p<0.0001) in SUDHL2 cell lines compared to WT. This increase was accompanied by an increase in NOTCH surface expression with higher NOTCH1 (p<0.0001) and NOTCH2 (p<0.0001) expression on OCI-Ly3 and SUDHL2 cells, respectively, compared to WT. To remove variability across cell lines we next used CRISPR to generate DTX1 knock out (KO) models using WT DTX1 DLBCL cell lines U2932 and OCI-Ly19. Analysis of NICD and NOTCH receptor protein expression confirmed our initial finding, showing an 8.5 and 4.8-fold increase of NICD1 and a 6 and 2.5-fold increase of NICD2 expression in the U2932 and OCI-Ly19 DTX1 KO cell lines, respectively, compared to empty vector (EV) controls. An increase in surface NOTCH1, and to a lesser extent NOTCH2, was seen in U2932 (p=0.01 and p=0.07, respectively) and OCI-Ly19 (p=0.001 and p=0.2, respectively) KO cells compared to EV controls. To identify pathways affected by DTX1 mutations, we next performed RNAseq on U2932 EV and KO cell lines. Of interest, we did not detect any significant changes in classical NOTCH driven genes (Hes1, Hey1, etc.), supporting prior studies on BN2/C1 tumors. However, we did identify potential non-canonical NOTCH-target genes that may contribute to lymphomagenesis. Gene expression in KO cells showed upregulation of genes involved in germinal center (GC) processes (AICDA, IRF4; p<0.001) and cellular growth and survival (IL6R, CCND2, STAT5B, BCL2; p<0.0001). Measurement of IL6R by flow cytometry and CyclinD2 by immunoblot further verified DTX1-mediated upregulation of IL6R (p<0.01) and CyclinD2 (p<0.01) for U2932 KO cells compared to EV control. We also observed downregulation of genes contributing to cell death (ex. BAX, CASP8, FAS; p<0.05) in KO cells compared to EV. Together, these novel findings support our hypothesis that DTX1 WWE1 mutations lead to dysregulation of the NOTCH pathway by driving overexpression of NICD and NOTCH, resulting in dysregulated expression of proteins involved in cell cycle, GC reaction, and apoptosis. We further provide new biologic insight into BN2/C1 DLBCL, ultimately highlighting a path toward novel therapeutic targets and the advancement of personalized medicine, as patients harboring these mutations may benefit from new therapies targeting inflammatory pathways, such as CELMoDs.

Abstract Description Follicular Lymphoma (FL) is a germinal center-derived malignancy representing approximately 20% of all non-Hodgkin’s lymphomas. Despite a relatively indolent disease course, FL is a biologically and clinically heterogeneous disease which poses significant clinical challenges due to its propensity for relapse and transformation into aggressive and untreatable forms of lymphoma. Here, we have leveraged computational techniques with WES, RNA-sequencing, and immune profiling data from FL tumors to gain further insight into the biologic basis for FL heterogeneity. We performed computational reconstruction of the B cell receptor repertoire on B cell sorted RNA-sequencing data from normal lymph nodes (n = 5) and FL tumors (n = 89). To gain a better understanding of FL B cell biology at the single cell level, we further investigated the clonal architecture of dominant malignant clonotypes in FL tumors using computational analysis of paired single cell V(D)J and 5’ GEX sequencing data from FL tumors (n = 21) and non-malignant lymph nodes (n = 3), allowing us to characterize the germinal center programs in malignant versus non-malignant B cells. Our data suggest that a subset of FL tumors consist of complex polyclonal B cell populations, some of which represent a pre-malignant FL precursor. Our data provide a deepened understanding of FL biology, including precursor cells, and contributes a refined transcriptomic characterization of malignant and non-malignant B cells in Follicular Lymphoma. Funding Sources This work was supported by the National Institutes of Health (R01 CA212162-01A1 to A.J.N and J.R.C., and P50 CA097274 to J.R.C., and A.J.N) Topic Categories Computational and Systems Immunology (COMP)

Abstract Twenty percent of patients with follicular lymphoma relapse early with poor outcomes; however, the molecular mechanisms underlying this aggressive behavior are unknown. Using a multiomic approach, we show that patients with follicular lymphoma with elevated IRF4 expression (IRF4hi) have increased transformation risk, dysregulated immune signaling, and a suppressive tumor microenvironment. Loss- and gain-of-function experiments in IRF4hi lymphoma cells, along with chromatin profiling, demonstrate that IRF4 impairs their interaction with T cells by repressing antigen presentation and co-receptor gene modules while promoting the expression of cytokines that antagonize T-follicular helper cell and regulatory T-cell functions. Additionally, IRF4 rewires tumor metabolism and restricts glucose availability to immune cells. Silencing of IRF4 inhibits tumor cell growth and restores immune surveillance mechanisms, thus representing a promising target for therapy. Our data suggest that IRF4hi lymphoma cells co-opt a developmental mechanism used to exit the germinal center response in promoting a more aggressive cancer via engagement of multiple immune-evasive mechanisms. Significance: We show that IRF4 controls immune and metabolic signaling with functional effects on the cross-talk between B and T cells, resulting in a suppressive immune microenvironment when IRF4 is overexpressed in lymphoma. This underscores the importance of IRF4 as a pivotal therapeutic target to restore antitumor immunosurveillance for IRF4-associated lymphomas.

Introduction: Diffuse large B cell lymphoma (DLBCL) patients experience favorable outcomes following frontline immunochemotherapy, however, 30-40% of patients will exhibit primary refractory disease or relapse (rrDLBCL). Recent studies have expanded upon well-defined characteristics of newly diagnosed DLBCL (ndDLBCL) to elucidate genomic landscapes underlying rrDLBCL; yet, high-throughput examination of resistance and relapse mechanisms within rrDLBCL subtypes, including the newly defined dark zone signature (DZsig) addition to cell of origin (COO), have yet to be described. To expand upon the current tumor etiology and define subtype-specific driving mechanisms, we describe a multi-omic analysis leveraging clinical and biologic variables in rrDLBCL. Methods: Tumor biopsies (N=228) collected at the time of relapse were obtained from n=85 participants in the Mayo Clinic/University of Iowa Lymphoma Molecular Epidemiology Resource (MER), n=73 in the Mayo Clinic Biobank, and n=38 from CC-122-DLBCL-001 or n=32 CC-122-ST-001 trials. Tumor samples were analyzed by WES (n=158) or WGS (n=70) and by RNAseq (n=145). Comparative analysis against ndDLBCL MER samples (N=444) utilized WES (n=404) and RNAseq (n=321). Chi Square or Fisher's exact test was used for enrichment analysis of categorical variables, while Wilcoxon or Kruskal-Wallis test was used for continuous variables. Results: rrDLBCL samples were enriched for ABC (46% vs 33% in ndDLBCL, p=0.01), and fewer GCB (27% vs 40%, p=0.006). Amongst rrDLBCL patients, variance in overall survival (from time of original diagnosis; log-rank p=0.001) was observed between COO subtypes, with poor outcomes in patients with DZsig (median 15.0 mo) and GCB (36.5 mo) compared to ABC tumors (74.4 mo). Patients with progression or relapse within 6 months featured a more even COO distribution (ABC - 22%, GCB - 34%, DZsig - 28%) than patients with early (6-24 mo; 51%, 22%, 22%) or late relapses (>24 mo; 62%, 19%, 4%), which were predominantly ABC. Comparing mutation frequency within subtypes, rrDLBCL GCB tumors were enriched for B2M (41% vs 15% in ndDLBCL, p=0.002), EP300 (19% vs 7%, p=0.059), CIITA (19% vs 5%, p=0.02), PRKCB (11% vs 2%, p=0.035), CD79B (11% vs 0%, p=0.004), and CD83 (8% vs 1%, p=0.049) mutations, and CNV events including gains at 8q24.21 (MYC locus; 48% vs 14%, p<0.001) and 18q21.33 (BCL2, TCF4 locus; 36% vs 15%, p=0.021), and losses at 17p13.1 (TP53 locus; 40% vs 19%, p=0.035). rrDLBCL DZsig tumors featured enrichment for BCL6 (22% vs 3%, p=0.025), and SPEN (17% vs 3%, p=0.062) mutations and 2q22.2 losses (LRP1B locus; 20% vs 0%, p=0.019). Despite differential mutation and CNV profiles, gene expression analysis revealed minimal variance within COO subtypes in rrDLBCL and ndDLBCL, however upregulation of lymphoma driving genes including MYBPC2 (LogFC=1.68, FDR<0.001), MYOM2 (LogFC=1.33, FDR=0.007), MYC (LogFC=0.61, FDR=0.04), and CELSR2 (LogFC=0.53, FDR=0.02) were observed in ABC rrDLBCL. Increased frequency of LymphGen classifier A53 was observed in ABC (43% vs 17%, p=0.005) and GCB (40% vs 14%, p=0.009) rrDLBCL. Using CIBERSORTx to infer lymphoma microenvironment (LME) cell type abundance, an LME-depleted class was evident in rrDLBCL tumors, including reduction of CD8+ T cells (effect size (r)=-0.78, p<0.001) and T follicular helper cells (r=-0.48, p=0.012) in GCB rrDLBCL and reduction of T follicular helper cells (r=-0.43, p=0.023) and Tregs (r=-0.32, p=0.065) in DZsig rrDLBCL. Significance: We present an update to the genomic landscape of rrDLBCL, highlighting GCB and DZsig tumors with differential biology and advanced clinical course. Our data supports increased implementation of sequencing efforts in rrDLBCL to more broadly evaluate relapse mechanisms, with optimism that a more detailed understanding of r/r disease may guide trials for targeted treatment strategies. For example, rrDLBCL patients may be strong candidates for MYC pathway targeting agents (BETi), or a Lenalidomide-based regimen (R2 or Tafa-Len), as Lenalidomide has been shown to reprogram the LME via PD-L1 expression in addition to downregulation of MYC and MYC target genes. Or, due to enrichment of variants affecting anti-apoptotic proteins including BCL6 or BCL2 in DZsig and GCB rrDLBCL, a venetoclax-based therapeutic may be considered. Further targeted strategies should be considered as development of a more personalized approach remains an unmet need.

Introduction MYD88 L265P mutations (MUT) are found in 90% of patients with Waldenstrom macroglobulinemia (WM) and lead to the activation of pro-survival signaling. Limited literature has demonstrated that patients with MYD88 wild-type (WT) WM have poorer responses to conventional therapies and a higher likelihood of histological transformation. Increasingly the role of non-coding genomics in lymphomagenesis has demonstrated importance and currently, the epigenetic landscape of MYD88 WT WM remains underexplored. This study utilized a multi-omics approach to assess DNA methylation in MYD88 WT WM compared to MYD88 MUT patients and controls. Methods Bone marrow samples were prospectively collected from 20 subjects [3 MYD88 WT/CXCR4 WT WM, 12 MYD88 MUT WM and 5 normal controls] and positive selection for CD19+ and/or CD138+ B cells was performed. The analysis focused on the comparison between MYD88 WT WM to normal controls and MYD88-WT to MUT WM. Genome-wide DNA methylation analysis was performed, and CpG methylation ratios were segmented into 200-bp regions. Differentially methylated regions (DMRs) were determined using a Q-value <0.05 and absolute methylation difference ≥10%. Next, RNA-seq data was interrogated and differentially expressed genes (log2fold change>0.5 or <-0.5, false discovery rate<0.05) with corresponding promoter hypo- or hypermethylation were assessed. Ingenuity Pathway Analysis was used to determine potentially epigenetically regulated pathways using a corrected p<0.05. Results Comparing MYD88 WT WM samples to controls, a statistically greater number of hypermethylated DMRs were observed across genomic regions (hyper- vs hypomethylated): promoter (29129 vs 6300), 3'untranslated region (UTR) (2068 vs 1551), 5'UTR (10169 vs 735), intron (30344 vs 24511), exon (2224 vs 6855), CpG island (30372 vs 1240) and CpG shore (14477 vs 5872) (all p<0.001). Next, mRNA analysis demonstrated that 214 genes were differentially up-/downregulated with concordant promoter hypo-/hypermethylation. On the assessment of promoter-methylation-based pathways, multiple pathways involved in cell signaling, proliferation, gene regulation, cytokine and immune signaling and metabolism were found to be deregulated. Most pathways were found to be negatively enriched in MYD88 WT patients compared to controls including several members of interleukin (IL) signaling (IL-2, IL-3, IL,7) and regulators of cell proliferation and signaling including mTOR, PI3K/AKT, and p53. Positive enrichment of the PTEN, FGFR3 and ERK/MAPK pathways were observed. Underlying multiple pathways was promoter hypermethylation with corresponding mRNA downregulation of AKT3 and PIK3CA. Furthermore, promoter hypomethylation with transcript upregulation of PTEN was observed. Next, in comparison of MYD88 WT to MUT samples, a statistically greater number of hypermethylated DMRs were observed across genomic regions (hyper- vs hypomethylated): promoter (24829 vs 2544), 3'UTR (3534 vs 378), 5'UTR (6111 vs 605), intron (52292 vs 5360), exon (21066 vs 2187), CpG island (18334 vs 1551) and CpG shore (16351 vs 1957) (all p<0.001). Assessing mRNA expression, 354 genes were up-/downregulated with concordant promoter hypo-/hypermethylation. Pathway analysis revealed positive enrichment of the mitotic G1 phase to G1/S transition pathway with upregulation of the CDK4, LYN, and PSMA6 transcripts in MYD88 WT compared to MUT samples. Further, negative enrichment of cell signaling regulators including the RHO GTPase family, IL pathways (IL-4, IL-8, IL-13, IL-17) and integrin signaling were found in MYD88 WT compared to MUT samples. Underlying multiple pathways were downregulation with promoter hypermethylation of BCL2, an important regulator of apoptosis. Conclusion Genome-wide methylation analysis of MYD88 WT WM revealed a unique DNA methylation profile, with potential epigenetic alteration of key oncogenic transcripts and pathways. The negative enrichment of cell signaling and cytokine/inflammatory pathways in MYD88 WT samples compared to controls and MYD88 MUT WM suggests an epigenetic mechanism for immune evasion and potentially underlies the reduced response to conventional WM-directed therapies. Overall, this study reveals that MYD88 WT WM is epigenetically distinct from MYD88 MUT WM and regulators such as methylation may underlie the unique clinical and phenotypic properties observed.

ABSTRACTDiffuse large B‐cell lymphoma (DLBCL) patients that fail to achieve a complete metabolic response with frontline immunochemotherapy have a poor prognosis. Genomic profiling has led to a broader understanding of the molecular drivers in DLBCL, but it is unknown how well current classifiers identify patients that will experience primary treatment resistance (PTR). Using whole exome and RNA sequencing data from newly diagnosed DLBCL patients, we evaluated the genomic landscape of PTR and compared it to that of non‐PTR DLBCL. We found a significant increase in the frequency of TP53 (34% vs. 15%, p = 0.005) and ARID1A mutations (21% vs. 7%, p = 0.007) in PTR cases, with pathway analysis further demonstrating a downregulation of TP53 and an increase in chromatin modifying pathways. These results suggest that TP53 and ARID1A may be key mediators of PTR and important pathways contributing to the poor outcomes. We found that the current molecular classifiers were unable to identify PTR cases at diagnosis. However, our newly identified high‐risk signature identified 46% of PTR cases at diagnosis. Overall, these results contribute to our understanding of the genomic landscape of patients with primary treatment resistance.

Abstract Currently, the role of DNA methylation in the immunoglobulin M (IgM) monoclonal gammopathy disease spectrum remains poorly understood. In the present study, a multiomics prospective analysis was conducted integrating DNA methylation, RNA sequencing (RNA-seq), and whole-exome sequencing data in 34 subjects (23 with Waldenström macroglobulinemia [WM], 6 with IgM monoclonal gammopathy of undetermined significance [MGUS], and 5 normal controls). Analysis was focused on defining differences between IgM gammopathies (WM/IgM-MGUS) compared with controls, and specifically between WM and IgM-MGUS. Between groups, genome-wide DNA methylation analysis demonstrated a significant number of differentially methylated regions that were annotated according to genomic region. Next, integration of RNA-seq data was performed to identify potentially epigenetically deregulated pathways. We found that pathways involved in cell cycle, metabolism, cytokine/immune signaling, cytoskeleton, tumor microenvironment, and intracellular signaling were differentially activated and potentially epigenetically regulated. Importantly, there was a positive enrichment of the CXCR4 signaling pathway along with several interleukin (interleukin 6 [IL-6], IL-8, and IL-15) signaling pathways in WM compared with IgM-MGUS. Further assessment of known tumor suppressor genes and oncogenes uncovered differential promoter methylation of several targets with concordant change in gene expression, including CCND1 and CD79B. Overall, this report defines how aberrant DNA methylation in IgM gammopathies may play a critical role in the epigenetic control of oncogenesis and key cellular functions.

Abstract The role of the bone marrow (BM) microenvironment in regulating the antitumor immune response in Waldenstrom macroglobulinemia (WM) remains limited. Here we transcriptionally and phenotypically profiled non-malignant (CD19− CD138−) BM cells from WM patients with a focus on myeloid derived suppressive cells (MDSCs) to provide a deeper understanding of their role in WM. We found that HLA-DRlowCD11b+CD33+ MDSCs were significantly increased in WM patients as compared to normal controls, with an expansion of predominantly polymorphonuclear (PMN)-MDSCs. Single-cell immunogenomic profiling of WM MDSCs identified an immune-suppressive gene signature with upregulated inflammatory pathways associated with interferon and TNF signaling. Gene signatures associated with an inflammatory and immune suppressive environment were predominately expressed in PMN-MDSCs. In vitro WM PMN-MDSCs demonstrated robust T-cell suppression, and their viability and expansion was notably enhanced by G-CSF and TNFα. Furthermore, BM malignant B-cells attracted PMN-MDSCs to a greater degree than monocytic MDSCs. Collectively, these data suggest that malignant WM B cells actively recruit PMN-MDSCs which promote an immunosuppressive BM microenvironment through a direct T cell inhibition, while release of G-CSF/TNFα in the microenvironment further promotes PMN-MDSC expansion and in turn immune suppression. Targeting PMN-MDSCs may therefore represent a potential therapeutic strategy in patients with WM.

Introduction: Molecular classifications have been developed to characterize the heterogeneity of diffuse large B cell lymphoma (DLBCL), with activation of BCL6 and Notch signaling pathways acting as driving pathogenic features in C1/BN2 subtype tumors. Among C1/BN2 defining variants are truncating/inactivating SPEN mutations. SPEN regulates the Notch signaling pathway by initiating formation of a repressive complex that facilitates transcriptional repression of Notch target genes. Our group has identified SPEN mutations enriched in DLBCL cases failing to achieve event-free status at 24 months (EFS24). This contrasted with generally favorable outcomes in C1/BN2 tumors, suggesting a potential high-risk molecular feature with unexplored significance. Presenting an opportunity to expand underlying mechanisms and risk stratification in C1/BN2 tumors, we aimed to use a functional genomic approach, combined with clinical data, to characterize the alternative biology stemming from SPEN truncating mutations. Methods: FFPE tumor biopsies from newly diagnosed DLBCL cases (N=444) enrolled in the Mayo Clinic/Iowa Lymphoma Molecular Epidemiology Resource (MER) program were analyzed by WES (n=404) and RNA-seq (n=321). Biologic and clinical variables were annotated for all cases. CRISPR/Cas9 introduced truncating mutations to the endogenous SPEN locus in OCI-LY3 & SU-DHL-2 ABC-type human cell lines to establish models of SPEN depletion. Functional in vitro analysis utilized standard methods for RNA-seq, cell proliferation, western blot, and flow cytometry. Results: WES identified 29/404 ndDLBCL (7.2%) harboring SPEN mutations; n=17 missense and n=12 nonsense or frameshift insertions/deletions (truncating variants - [trunc]), the later considered as the focal point of our study. Against all cases, SPEN trunc cases displayed enrichment for high-risk factors including non-GCB COO (p=0.06, OR=3.7, 95% CI [0.9, 22.2]) and primary refractory disease (PTR; p=0.02, OR=5.0, 95% CI [1.05, 19.9]), and demonstrated inferior OS (p=0.06, HR=2.2), where a majority of clinical progression events occurred within 12 months of diagnosis. Leveraging RNA-seq, we observed 11 of 11 SPEN trunc cases exhibiting a double-hit negative gene expression signature, 9/11 cases with inflammatory or depleted microenvironment (LME), and 8/11 with a signature associating with intermediate/high risk of EFS24 failure. Applying the LymphGen algorithm, 5/12 SPEN trunc cases received the BN2 classification, observing co-occurring mutations with CD70 (p<0.01, OR=9.68), BRCA2 (p=0.03, OR=5.84), TET2 (p=0.02, OR=4.91), IRF4 (p=0.02, OR=7.27), and PTPRD (p<0.01, OR=8.31) genes, and were mutually exclusive with ARID1A, BCL2, EZH2, FAS, HLA-B, IRF8, and TMEM30A mutations ( Figure 1A). Isolating SPEN trunc against BN2 cases without SPEN variants (n=14), despite similar clustering of non-GCB COO (p=0.7), PTR was exclusively observed in SPEN trunc cases, accompanied again by inferior OS (p=0.07, HR=4.2; Figure 1B). In the OCI-LY3 SPEN depletion model, GSEA analysis from RNA-seq revealed activation of E2F (NES=2.7; p=9.4E-20) and MYC target genes (NES=2.3; p=4.1E-6), while both OCI-LY3 and SU-DHL-2 models featured depletion of antigen presentation genes (NES= -1.8, p=5.7E-4; NES= -2.2, p=3.8E-7, respectively). A decrease of MHC class II proteins (HLA-DR, HLA-DP) was confirmed by flow cytometry in both models (p<0.01). Hypothesizing loss of SPEN enhances reliance on the BCL6-HDAC3 axis featured in BN2 tumors, we treated both SPEN depletion models with a class-I HDAC inhibitor highly selective for HDAC3, BRD-3308, and indeed observed a 1.4 and 2.4-fold increase in sensitivity compared to WT cells (p=0.07, p<0.01, respectively). Significance: We find SPEN trunc mutations in poor acting DLBCL cases with clinical and biologic characteristics that diverge from C1/BN2 biology. The enhanced rate of clinical progression highlights urgency to better characterize these tumors. Our data suggests loss of SPEN regulatory mechanisms may act cooperatively with C1/BN2 biology to facilitate tumor immune evasion and drive progression through Notch-mediated pathways. Importantly, in vitro evidence shows SPEN deficient lymphoma cells are susceptible to targeted therapeutics against dependent pathways, providing rationale for advancing individualized care strategies for affected patients.

Introduction: Waldenström macroglobulinemia (WM) is a rare B-cell lymphoma that may be preceded by an asymptomatic IgM-monoclonal gammopathy of undetermined significance (MGUS). The mechanisms underlying the progression from IgM-MGUS to WM remain to be fully elucidated; however, the stepwise accumulation of genomic abnormalities has been suggested as a potential driver. Currently, the role of epigenetic regulation of DNA remains underexplored in this disease spectrum. Aberrant DNA methylation is a hallmark of oncogenesis that plays an important role in hematologic malignancies. The aim of this study was to assess differential DNA methylation in WM compared to IgM-MGUS and identify hyper- and hypomethylated pathways that may alter key cellular functions. Methods: This study included 15 patients with WM and 6 with IgM-MGUS. Bone marrow samples were prospectively collected and CD19+ and/or CD138+ B cells were positively selected. Subsequently, genome-wide DNA methylation analysis was performed. CpG methylation ratios were segmented into 200-bp regions and differentially methylated regions (DMRs) comparing WM to IgM-MGUS were identified. DMRs with a Q-value <0.05 and absolute methylation difference of ≥20% were considered significant and utilized for analysis. DMRs were annotated according to genomic regions of interest (promotors, 3'UTR, 5'UTR, introns, exons, CpG islands, CpG shores). Next, genes with significant promotor hyper- or hypomethylation were assessed using the Ingenuity Pathway Analysis. Pathways with a corrected p-value<0.05 were considered statistically significant and a Z-value was calculated based on the degree of differential methylation. Results: Genome-wide DNA methylation analysis revealed a significant number of DMRs (n= 15691, hyper: 8248; hypo: 7443) in WM compared to IgM-MGUS. Assessing genomic regions, promotors (hyper: 1379; hypo: 2403), introns (hyper: 1548; hypo: 2292), CpG shores (hyper: 587; hypo: 849) and 3'UTR (hyper: 66; hypo: 193) regions were observed to have a significantly greater number of hypomethylated DMRs (all comparisons, p<0.001) in WM compared to IgM-MGUS (Figure A). While exons (hyper: 1686; hypo: 1378), CpG islands (hyper: 2128; hypo: 82), and 5'UTR (hyper: 854; hypo: 246) regions had a significantly greater number of hypermethylated DMRs (all comparisons, p<0.001) (Figure A). Next, promotor methylation-based pathway analysis revealed differential methylation of multiple cellular pathways important in WM pathogenesis (Figure B). These included pathways involved in intracellular signaling and the regulation of cellular proliferation, such as the STAT3 and CXCR4 pathways which were hypomethylated, and the PTEN signaling pathway which was hypermethylated. Of interest, we have previously demonstrated the transcription factor, STAT3, to have an increased expression in WM and to be implicated in IgM secretion (Elsawa, 2011). Next, multiple cytokine signaling pathways were observed to be hypomethylated in WM compared to IgM-MGUS including members of the interleukin (IL) family (IL-1, IL-8, IL-15, IL-23), VEGF, and TGF-β. Of relevance, previous studies have demonstrated the elevation of IL-8, IL-15, and increased expression of angiogenic cytokines, including VEGF, in patients with WM (Lewicki, 2010; Anagnostopoulos, 2007). Furthermore, our group has also shown miRNA-based epigenetic regulation of IL-8 and IL-15 in WM compared to IgM-MGUS (Chohan, 2023). Lastly, several pathways critical in the regulation of the cell structure and the extracellular matrix (ECM) were found to be hypomethylated including the FGF, actin cytoskeleton, PAK, and FAK signaling pathways. Of note, the FGF/FGFR axis has been demonstrated to be dysregulated in WM, and its targeting has been shown to halt WM cell growth (Sacco, 2021). Conclusions: This study demonstrates differential DNA methylation of key genomic regions in WM compared to IgM-MGUS. Promotor methylation-based analysis reveals hypomethylation of several important cellular pathways involved in cytokine signaling, cell structure/ECM, and intracellular signaling. These results underscore the influence of epigenetics in the IgM-gammopathy disease spectrum and suggest the potential role of aberrant DNA methylation underlying the malignant progression from IgM-MGUS to WM.