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.

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 An accurate assessment of p53’s functional statuses is critical for cancer genomic medicine. However, there is a significant challenge in identifying tumors with non-mutational p53 inactivation which is not detectable through DNA sequencing. These undetected cases are often misclassified as p53-normal, leading to inaccurate prognosis and downstream association analyses. To address this issue, we built the support vector machine (SVM) models to systematically reassess p53’s functional statuses in TP53 wild-type (TP53WT) tumors from multiple The Cancer Genome Atlas (TCGA) cohorts. Cross-validation demonstrated the good performance of the SVM models with a mean area under the receiver operating characteristic curve (AUROC) of 0.9822, precision of 0.9747, and recall of 0.9784. Our study revealed that a significant proportion (87%–99%) of TP53WT tumors actually had compromised p53 function. Additional analyses uncovered that these genetically intact but functionally impaired (termed as predictively reduced function of p53 or TP53WT-pRF) tumors exhibited genomic and pathophysiologic features akin to TP53-mutant tumors: heightened genomic instability and elevated levels of hypoxia. Clinically, patients with TP53WT-pRF tumors experienced significantly shortened overall survival or progression-free survival compared to those with predictively normal function of p53 (TP53WT-pN) tumors, and these patients also displayed increased sensitivity to platinum-based chemotherapy and radiation therapy.

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.

Abstract Extrachromosomal DNA (ecDNA) drives the evolution of cancer cells by conferring survival advantages over normal cells. However, the functional importance of ecDNA replication and maintenance is yet to be evaluated, and the molecular components that contributed to these processes remain largely unknown. Here, using CRISPR-C (circularization of genes and chromosomes by CRISPR) technology and synthetic circular DNA mimicking ecDNA identified in tumors, we generated ecDNA-positive (ecDNA+) cell models. Using these models, we identified the proteins involved in ecDNA replication and maintenance and found DNA damage response (DDR)-related pathways were significantly enriched. We further demonstrated that ecDNA replication led to increased DDR and genome instability and DDR was also essential for ecDNA maintenance, consistent with observations that ecDNA+ tumors exhibited higher DDR scores than ecDNA-negative (ecDNA–) tumors in patients. ATM, MDC1, and CHK2, but not ATR, were activated in ecDNA+ cells compared to ecDNA– cells. Moreover, ecDNA+ cells were more sensitive to ATM and CHK2 inhibition. We further found that type II topoisomerase TOP2B and DNA ligase III (LIG3) were critical regulators for ecDNA replication and maintenance and related DDR. In summary, our findings of the reciprocal interactions between ecDNA and DDR provide new insights into detection and therapy for ecDNA+ tumors.