Bhardwaj, Vaishali

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.

Abstract Waldenstrom macroglobulinemia (WM) is a lymphoplasmacytic lymphoma associated with an increase in serum immunoglobulin M (IgM) monoclonal protein. The underlying causes of the disease are still not fully understood, but cells in the bone marrow (BM) microenvironment (TME) are anticipated to play a significant role in disease progression. We previously reported that myeloid-derived suppressor cells (MDSCs), a phenotypically heterogeneous group of cells known for their ability to suppress immune responses, are expanded in the bone marrow of patients with WM and that the expansion of MDSCs is particularly notable for increased numbers of polymorphonuclear MDSCs (PMN-MDSCs). In the present study, we evaluated the ability of WM cells to attract MDSCs to the tumor microenvironment (TME). We also evaluated whether PMN-MDSCs suppress the cytotoxic function of T-cells and thereby promote the survival of WM cells. Our study confirms the expansion of MDSCs, notably PMN-MDSCs in BM of WM patients as compared to normals. We also found that PMN-MDSCs expressing CD66b was a predominant population in WM patients. To determine the phenotype and transcriptome of both PMN-MDSCs and monocytic MDSCs (M-MDSCS) in WM, we compared BM aspirates from WM patients to normal bone marrow (NBM) and performed cellular indexing of transcriptomes and epitopes by sequencing (CITEseq) analysis (WM=3, NBM=2) followed by qPCR validation. We found that MDSCs in WM express significantly higher levels of HIF1-α overall (p<0.02) compared to NBM. Furthermore, PMN-MDSCs are more abundant in WM bone marrow than M-MDSCs and express significantly higher levels of CXCR4 (p<0.001), myeloperoxidase (MPO; p<0.0003), and azuricidin (AZU1; p<0.0004), suggesting an activated phenotype. To determine whether WM cells attract MDSCs to the BM, we performed a transwell assay and showed that WM cells (BCWM.1 cell line) attracted MDSCs, particularly PMN-MDSCs, to a greater extent than M-MDSCs (n=3; p<0.001). We next evaluated whether MDSCs prevent cytotoxic T-cells from lysing WM cells. In-vitro co-culture studies with BCWM.1 cells and T-cytotoxic cells (B:CD8-T) showed that T-cytotoxic cells were able to lyse the WM cell line by 46.4% , (median range 46.4 - 47.8%, p=0.001). However, when BCWM.1 cells and T-cytotoxic cells were co-cultured in the presence of MDSCs (i.e B:CD8-T:M), the specific lysis of BCWM.1 cells was reduced by 12.6% respectively (range 2.2- 12.6%, p<0.001). Furthermore, we observed that B:CD8-T:M co-culture significantly reduced CD8-T cell activation as compared to activated T cells alone (17.9%) (median range 17.8- 18%, p<0.001). In contrast, no affect was observed with B: CD8-T coculture by 91% (range 90.6- 91.4% p<.001) confirming that MDSCs suppress CD8-T cell activation. Our study therefore provides evidence that MDSCs in WM have an activated phenotype and that PMN-MDSCs in WM are attracted by malignant cells in the bone marrow. Furthermore, we find that the increased presence of PMN-MDSCs in the WM bone marrow environment suppresses T-cell lysis of malignant cells and may thereby contribute to disease progression. These findings provide a strong rationale for targeting MDSCs as a potential therapeutic strategy in WM.

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.

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.

Reprogramming of cellular metabolism is a hallmark of cancer leading to improved cellular fitness and providing a selective advantage to malignant cells. Recently, we found that follicular lymphoma (FL) patients with poor prognosis display increased transcriptional and protein expression of IRF4 (IRF4 high), altered immune signaling and an immunosuppressive tumor microenvironment (TME). Furthermore, high expression of IRF4 rewired B cell metabolism towards glycolysis ( Mondello P et al, Blood 2022;140 Suppl 1:168-169). However, it's unknown whether this metabolic reprogramming may be responsible in part for the impairment of immune cells in the TME, which also rely on glucose for their proliferation and function. Using RNA-seq on 88 FL patients, we found a negative correlation between IRF4-regulated genes associated with glycolysis and tumor infiltration by major immune cell subsets in IRF4 high compared to IRF4 low tumors ( Fig 1A), suggesting that IRF4-driven tumor glycolysis may reinforce an immune suppressed phenotype. Especially, we observed that ALDOA and GPI genes were inversely correlated with expression of CD4 + T cells (r=0.34, p=0.003; and r=0.38, p=0.015, respectively) and T reg cells (r=0.3, p=0.05; and r=0.35, p=0.027, respectively) in IRF4 high but not in IRF4 low tumors. As radioactive glucose is the cornerstone of PET assessment, we investigated the correlation between baseline SUV values and IRF4 expression. Remarkably, we found significantly higher levels of FDG uptake in FL IRF4 high compared to IRF4 low cases (p= 0.003), supporting the glycolytic features of Iymphoma expressing IRF4. To define the metabolic activity and function of lymphoma cells with different IRF4 status, we performed seahorse assay in TMD8 and SUDHL2 cells (lymphoma cells with IRF4 high expression) transfected with siIRF4 or siCtr for 72 hours. In line with the IRF4-dependence on glycolysis, TMD8 transfected with siCtr had a low baseline glycolytic ECAR, while the same value doubled after silencing IRF4. Notably, after adding oligomycin, which stimulates glycolysis, TMD8 control cells showed a lower glycolytic capability compared with those deprived of IRF4, probably due to the already maximal addiction to glycolysis ( Fig 1B left). We then investigated the glycolytic demand of germinal center (GC) B cells differently expressing IRF4. Using irf4 fl/fl mice ( Klein et al, Nat Immunol 2006) crossed with Cγ 1cre animals ( Cγ 1cre; irf4 fl/fl) to induce conditional deletion of irf4 in GC B cells, we tested ex-vivo sorted GC B cells from Cγ 1cre; irf4 fl/fl, Cγ 1cre; irf4 fl/+ and WT mice immunized for 10 days. This time point is critical for studying IRF4 as it corresponds to the GC reaction peak when GC B cells express elevated IRF4 levels promoting terminal differentiation to plasma cells. As expected, GC B cells minimally relied on glycolysis and showed a low ability to respond to oligomycin in WT mice ( Weiser et al, Nat Immunol 2020). However, GC B cells with partial or complete loss of irf4 demonstrated a superior ability to response to increased glycolytic demand ( Fig 1B right), further supporting the conclusion that IRF4 influences B cell dependency on glycolysis. As the malignant cell metabolism can alter nutrient availability in the TME, we explored if lymphoma cells differently expressing IRF4 may change glucose levels in the media and in turn impact T cells activity. We found low glucose and high lactate levels in the media of cell lines with IRF4 high expression, while the opposite was observed for those that were IRF4 low. We then cocultured primary human CD4 + T cells with TMD8 cells transfected with either siIRF4 or siCtr for 3 days at three different concentrations of D-glucose: high (22.2 mM, Gluc hi), normal (5.5 mM, Gluc nor) or absent (0 mM, Gluc nil). In line with our prior findings, coculture of B cells in which IRF4 was knocked out showed a significant decrease in abundance of T reg cells and an increase in T FH cells compared to control in presence of Gluc hi or Gluc nor. In contrast, there was no difference in the overall amount of either T cell subtype irrespective of IRF4 status with Gluc nil, supporting the prerequisite of microenvironmental glucose as primary fuel for T cell survival and function. In conclusion, IRF4 rewires the metabolism in B cells with profound functional implications on the TME. Our data suggest that glucose metabolism represents a functional vulnerability that might be exploited therapeutically.

A substantial subset of follicular lymphoma (FL) patients has an early relapse with a poor outcome. We showed that this unfavorable group displays increased expression of IRF4, dysregulated immune signaling and an immunosuppressive microenvironment. At the molecular level, IRF4 controls expression of antigen presentation and immune molecules (e.g. HLA-DR, CD40, PDL1). Furthermore, silencing of IRF4 triggers anti-lymphoma immunity ( Mondello P et al, Blood 2022;140 Suppl 1:168-169). However, little is known on the molecular mechanisms whereby IRF4 controls the crosstalk between B and T cells. To address this question, first we performed RNA-seq, ATAC-seq and CUT&TAG of active (H3K27Ac, H3K4me3) and repressive (H3K27me3, H3K4me1) histone marks in TMD8 and HBL1 (lymphoma cells with high IRF4 expression) transfected with two pairs of siIRF4 or siCtr for 72 hours. Gene expression revealed a broad change in transcription, with 432 upregulated genes and 136 downregulated genes in siIRF4 compared to siCtr. Notably, the genes with increased expression following IRF4 silencing coded for antigen presentation, immune molecules and cytokines (e.g. HLA-DR, CD86, IL6, IL10). Accordingly, we observed significant enrichment of antigen presentation, interferon and immune response signatures ( Fig 1A). Integrating RNA-seq with chromatin accessibility showed a significant shift towards transcriptional upregulation at genes nearest to upregulated ATAC peaks, with significant enrichment at enhancers (n=771). Conversely, transcriptional suppression was enriched at enhancers (n=2,458) of downregulated peaks. k-means analyses of all union peaks, generated by merging ATAC with CUT&TAG peaks, revealed a significant gain of H3K4me3, with concordant loss of H3K4me1 at enhancers with upregulated gene expression. In contrast, H3K4me3 was severely decreased at suppressed peaks. TF motif analysis revealed an enrichment for PU.1 and E2A at IRF4-bound regions, suggesting their cooperation in regulating immune signaling. Since antigen presentation and other immune molecules are critical components of the immune synapse, we explored the role of IRF4 on the interaction between germinal center (GC) B and T cells by crossing irf4 fl/fl mice with conditional deletion of irf4 ( Klein et al, Nat Immunol 2006) and Cγ 1cre mice ( Cγ 1cre; irf4 fl/fl) to induce recombination in GC B cells. After 10 days of immunization (GC reaction peak), Cγ 1cre; irf4 fl/fl and Cγ 1cre; irf4 fl/+ GC B cells showed significant increase of HLA-DR and CD40 expression compared to WT. As expected, Cγ 1cre; irf4 fl/fl mice displayed a significantly skewed light-zone:dark-zone ratio in favor of increased dark-zone, which are the cells that most reflect the action of BCL6, with concordant suppression of plasma cells. Additionally, we observed a significant increase of CD4 + T cells after partial or complete deletion of irf4, with upregulation of T FH cells and decrease of T reg cells ( Fig 1B). Notably, mice with loss of irf4 showed a significant decrease in exhausted T FH cells. To test if this effect was dependent on MHC:TCR interaction, we cocultured CD4 + T cells with TMD8 cells transfected with two siIRF4 or siCtr in presence of blocking antibodies for MHC class I, MHC class II or both. Blocking one or the other MHC class I or II partially rescued the IRF4 effect on T reg and T FH cells respectively as well as reduced anti-lymphoma cytotoxicity. Further experiments blocking IRF4-dependent IL6 and IL10 identified these two cytokines as responsible for the additional effects on T cells. To confirm this as an antigen-dependent immune response we coculture sorted GC B cells from NP-OVA-immunized Cγ 1cre; irf4 fl/fl, Cγ 1cre; irf4 fl/+ and WT mice with CSFE-labeled OVA-specific OT-II CD4 + T cells at a ratio of 2:1 in an organoid system for 5 days. We observed induced antigen specific proliferation of CD4 + T cells in the presence of Cγ 1cre; irf4 fl/fl and Cγ 1cre; irf4 fl/+ GC B cellscompared to WT,as indicated by decrease in CSFE expression . We confirmed an increase of T FH cells and a decrease of T reg cells, supporting the conclusion that IRF4 directly controls the antigen-dependent immune response . In summary, IRF4 expression in GC and lymphoma B cells regulates immune signaling with PU.1 and E2A, thereby modulating antigen-dependent immune responses with potential implications for lymphomagenesis. Future studies investigating targeting of IRF4 in lymphoma are warranted.

Abstract T follicular helper (TFH) cells are a subset of CD4+ T cells that reside within the follicles of secondary lymphoid organs and play an essential role in the germinal center (GC) reaction. To examine the role TFH cells in follicular lymphoma (FL), we first characterized TFH cells using single cell analysis with the high-throughput technologies of mass cytometry (CyTOF) and the Indexing of Transcriptomes and Epitopes by Sequencing (CITE-seq). Clustering analysis revealed that TFH cells form subsets with unique phenotypes, and the expression of many common markers including CXCR5, ICOS, IL-21, TOX2, CD57, GZMK, XCL1, XCL2, CXCL13 and CD40L varies on TFH subsets. We found that TFH cells are phenotypically more heterogenous in FL than benign tonsil tissues by both proteomic (n=82) and transcriptomic (n=8) profiling. Each FL patient showed different dominant TFH subsets, in contrast to tonsil donors which were highly homogenous. Some TFH subsets were only present in FL and absent in tonsil tissues, suggesting the existence of FL-specific TFH cells. Cytokines (IL-2, IL12, IL-17 and IL-21) and activated B cells upregulated TFH markers including PD-1, CXCR5 and IL-21 to induce the development of TFH subsets such as IL-21+ or IL-21+IFN-g+ TFH cells, which was inhibited by rapamycin, an immunosuppressor. We observed that the TFH cell frequency significantly correlated with the number of lymphoma B cells. Futhermore, while TFH cells as a whole showed no correlation with prognosis, TFH subsets such as CD57+ TFH cells were significantly associated with an inferior patient outcome. Taken together, our results indicate that TFH cells are highly heterogenous with different tumor-specific TFH subsets present in FL and associated with patient outcomes in FL.