Citations

Splenic marginal zone lymphoma (SMZL) is a B-cell malignancy with a typically indolent course. However, some patients have more aggressive disease and initial treatment can range from observation to chemoimmunotherapy. The reason for this spectrum of disease behavior is incompletely understood but the tumor microenvironment (TME) likely plays a role as recent studies have shown that an optimized antitumor immune response in SMZL patients would allow the immune system to better target malignant B cells and irradicate the disease. Regulatory CD4+T cells (Tregs) are crucial in maintaining TME homeostasis by suppressing effector T cells and malignant B cells and in this study, we therefore explored the key factors influencing the function and development of intratumoral Tregs in SMZL patients. We identified two distinct Treg subsets in the TME of SMZL patients: CD26+Tregs and CD161+Tregs, each with unique phenotypic profiles, functions, and prognostic significance. CD161+Tregs suppress CD8+T cells and malignant B cells more effectively than CD26+Tregs, which only inhibit CD8+T cells. Clinically, patients with inferior outcomes have significantly higher numbers of CD26+Tregs and fewer CD161+Tregs. To delineate differences between these subsets, we performed Cellular Indexing of Transcriptomes and Epitopes Sequencing (CITE-seq) on biopsies from 7 SMZL patients. We found that cytokine-inducible SH2-containing protein (CISH) is the most highly expressed gene in CD26+Tregs compared to CD161+Tregs, aside from FOXP3. CISH has been reported to be a negative regulator of immune effector cell activation in tumor models. Although CISH deficiency leads to a loss of Foxp3 expression and reduced Treg suppressive function, the mechanism by which CISH modulates Treg function, especially immune suppression in hematologic malignancies, remains largely unknown. To further investigate the mechanism of CISH in Treg function in MZL patients, we conducted CITE-seq analysis on a larger cohort of MZL patients (n=12), divided into groups based on outcomes (EFS24 failed, n=4; EFS24 achieved, n=8). The CITE-seq analysis results showed significantly higher CISH expression in Tregs from patients with inferior outcomes (P=1.07E-06). We then explored factors influencing CISH expression and found that TCR engagement and treatment with cytokines (IL-15, IL-2, and IL-7) elevated CISH expression in Tregs. While we previously observed that lymphoma B cells could enhance Foxp3 expression in CD4+CD25-T cells, this study we found that B cells could also enhance CISH and CD26 expression in CD4+CD25-T cells in vitro, indicating that a high frequency of B cells in the environment may lead to increased CISH expression in Treg cells leading to an unfavorable outcome in patients. We further investigated the role of CISH in the function of Treg cells. Following CISH knockdown in Treg cells, we observed decreased CD26 expression and reduced suppressive function on CD8+ T cell proliferation. Using CRISPR-Cas9, we knocked out CISH in the Karpas 299 cell line, a human T cell lymphoma line that constitutively expresses CD25, Foxp3, and CISH. CISH-KO Karpas 299 cells showed significantly less apoptosis and Foxp3 expression compared to CISH-WT cells. Bulk RNA-seq analysis of CISH-KO and CISH-WT Karpas 299 cells revealed that CISH-KO cells showed enhanced DNA replication and cell cycle pathway activity. Taken together, our results suggest that high expression of CISH in Tregs may predict an unfavorable outcome in SMZL patients. Malignant B-cells in the TME of SMZL patients may promote expression of CD26 and CISH in Tregs, thereby modulating Treg function. Targeting CISH could decrease the immunosuppressive activity of Tregs in the TME, potentially optimizing antitumor immunity in SMZL patients.

Background: FL3B is the least common among FL subtypes, accounting for 5-10% of cases. Almost half of FL3B tumors have concurrent diffuse large B cell lymphoma (DLBCL) in the diagnostic tissue biopsy, while FL3B rarely coexists with FL1-2 or FL3A. Several lines of evidence suggest that FL3B displays distinct biological and clinical characteristics resembling aggressive DLBCL rather than the indolent FL1-3A. Accordingly, DLBCL anthracycline-containing regimens, such as R-CHOP, are the mainstay of therapy for FL3B with or without coexisting DLBCL. However, while some studies show that R-CHOP therapy in FL3B was associated with 5-year overall survival (OS) rate of 54%, others report outcomes comparable to FL3A with 79-90% survival. Additional controversy involves the clinical implication of pure FL3B (FL3Bp) vs FL3B with concurrent DLBCL (FL3Bc). Thus, it remains unclear if FL3B clinically resembles DLBCL or lower-grade FL and whether FL3Bp has a different prognosis from FL3Bc. To address these questions, we performed a pooled study of two large prospective multicenter study cohorts to compare the clinical characteristics and outcomes of newly diagnosed FL3B with FL1-2, FL3A and DLBCL patients, all of whom were uniformly treated with frontline R-CHOP. Patients and Methods: Weevaluated 469 newly diagnosed patients with FL1-2 (n=216), FL3A (n=170) and FL3B (n=83) who received R-CHOP within 6 months of diagnosis and 739 de novo DLBCL (GCB n=455 and non-GCB n=284) who were prospectively enrolled into the Molecular Epidemiology Resource (MER) and the Lymphoma Epidemiology of Outcomes (LEO) cohorts between 2002 and 2021. FL3B cases were designated as FL3Bc and FL3Bp based on the presence or absence of coexisting DLBCL, respectively. Histologic diagnosis was determined by an expert pathologist at each participating center and centrally reviewed. Event free survival (EFS) and OS were defined as time from diagnosis to progression, relapse, retreatment or death, and early progression defined as failing to achieve EFS at 24 months (EFS24). Associations between clinical features and EFS or OS were evaluated using Kaplan Meier curves and Cox proportional hazards models unadjusted and adjusted for age and FLIPI (comparisons with FL1-2 and FL3A) or IPI (comparisons with DLBCL). Results: The median age was 61 years (range, 18-94 years). Among FL3B patients, 21 (25%) had FL3Bc and 62 (75%) FL3Bp. Baseline characteristics and outcomes were broadly similar between the two FL3B subtypes. We then compared all FL3B patients to those with lower-grade FL. Patients with FL3B showed similar clinical features to those with FL1-2 and FL3A, with exception for a lower tumor burden. After R-CHOP, FL1-2 patients had an inferior EFS (HR 2.13 95%CI 1.31-3.46) than those with FL3B, whereas there was no difference compared to FL3A. OS was similar across all FL grades. Although FL1-2 patients failed to achieve EFS24 more frequently than FL3B and FL3A (34% vs 18% and 15%, p<0.001), FL3B patients who did fail EFS24 had almost three folds higher risk of subsequent mortality (HR 2.84 95%CI 1.35-6.01) than other FLs. At 5-year follow-up, FL3B patients had twice the risk of relapse with an aggressive subtype (i.e., DLBCL, high-grade B-cell lymphoma) than those with FL1-2 and FL3A, who instead often recurred with indolent disease (i.e., FL1-2, FL3A; p<0.001). We further compared FL3B patients to those with DLBCL classified based on cell-of-origin into GCB and non-GCB subtypes. At diagnosis, patients with FL3B had more favorable clinical features than those with GCB and non-GCB DLBCL. Compared to FL3B, patients with non-GCB DLBCL had an inferior EFS (HR 1.73 95%CI 1.07-2.79) and OS (HR 2.04 95%CI 1.11-3.75), while those with GCB DLBCL did not show a significant difference in EFS (HR 1.41 95%CI 0.88-2.25) or OS (HR 1.69, 95%CI 0.93-3.06) likely due to insufficient statistical power. Patients with FL3B failed EFS24 less frequently than patients with GCB or non-GCB DLBCL (17% vs 25% and 32%; p=0.01). Among patients who failed EFS24, those with FL3B had a similar subsequent OS as both GCB (HR 1.41 95%CI 0.73-2.72) and non-GCB DLBCL (HR 1.27 95%CI 0.65-2.47). Conclusions: We report that FL3Bp and FL3Bc appear clinically similar. While most FL3B patients have a prognosis comparable to FL1-2 and FL3A, 18% exhibit aggressive behavior with poor outcome reminiscent of DLBCL. Future studies with a larger population will be needed to confirm our findings.

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.

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.

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 Lymphoma is one of the leading causes of cancer and cancer deaths and yet has not been amenable to population screening. The role of methylated DNA markers (MDMs) in the detection of lymphoma has not been extensively studied. We aimed to discover, validate, and test tissue–derived MDMs of lymphoma in archival plasma specimens. Reduced representation bisulfite sequencing (RRBS) was performed on a discovery set of frozen tissues. MDMs identified were converted to methylation–specific PCR assays and validated on independent formalin–fixed, paraffin‐embedded (FFPE) tissues. Target enrichment long–probe quantitative‐amplified signal (TELQAS) assays were developed and assayed in plasma–extracted, bisulfite‐converted DNA from independent treatment‐naïve lymphoma patients and healthy controls. Prediction of cancer status was modeled using random forest model with in silico cross‐validation. After discovery and validation in tissue, 16 TELQAS assays ( ZNF503 , VWA5B1 , HOXA9 , GABRG3 , ITGA5 , MAX.chr17.7190 , BNC1 , CDK20 , MAX.chr4.4069 , TPBG , DNAH14 , SYT6 , CACNG8 , FAM110B , ADRA1D , and NRN1 ) were selected for testing in plasma. These detected 78% (95% CI, 74%–82%) of lymphoma cases at 90% specificity. Excluding marginal zone and T‐cell lymphomas, sensitivity increased to 84% (80%–88%). MDMs in plasma show promise to detect lymphoma and are candidates for inclusion in multi‐cancer detection studies.

SummaryDespite the clinical and molecular heterogeneity of follicular lymphoma (FL), there remains a lack of biomarker‐directed therapeutic approaches in routine clinical practice, with the notable exception of the EZH2 inhibitor tazemetostat in EZH2‐mutant FL. Here we examined whether gene mutation status predicts response to clinical mTOR inhibitors (mTORi) in FL, by performing targeted mutational profiling of biopsies from 21 relapsed/refractory FL patients treated with mTORi everolimus or temsirolimus within clinical trials. We observed an enrichment of mutations within the catalytic histone acetyltransferase (HAT) domain of CREBBP in mTORi‐responders, and describe distinct transcriptional characteristics and co‐occurring mutations of FL harbouring these mutations; reinforcing the growing appreciation of CREBBPHAT mutation as a key biological determinant and its promise as a therapeutic biomarker in FL.

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.