Yang, Zhi-Zhang

Abstract Germinal center (GC)-derived lymphomas arise within a specialized immune ecosystem, yet whether malignant cells direct its remodeling remains unclear. Integrating genomic, transcriptomic, spatial, and functional analyses across follicular lymphoma and diffuse large B-cell lymphoma, we define a developmental framework linking malignant B-cell differentiation state to tumor microenvironment (TME) organization. Rather than segregating by histology, lymphomas align along a shared differentiation continuum in which proliferative dark zone (DZ) states associate with immune-depleted TMEs, whereas post-GC/memory B-cell (MBC) states drive inflamed but immunosuppressed TMEs. Spatial profiling reveals collapse of GC architecture and emergence of suppressive myeloid niches along this trajectory. Single-cell lineage reconstruction uncovers a DZ-to-MBC developmental axis underlying lymphoma evolution. Genetic reprogramming in vivo demonstrates that malignant differentiation state is sufficient to instruct TME remodeling. Collectively, these findings establish malignant differentiation state as a causal determinant of immune architecture across GC-derived lymphomas and identify lineage-directed immune reprogramming as a potential therapeutic strategy.

Abstract Background: Follicular lymphoma (FL) is a germinal center–derived B cell malignancy characterized by the clonal expansion of a dominant malignant B cell population. While prior studies have identified transcriptional subtypes across patients, most analyses have treated B cells within individual tumors as a transcriptionally uniform population. However, this assumption has not been systematically tested using high-resolution, clonotype-resolved transcriptomic data. The broader landscape of intratumoral B cells—including those not belonging to the dominant clone—remains poorly defined, and their roles in disease biology are not well understood. These unexpanded, intratumoral B (ITB) cells may include residual non-malignant populations, early-stage malignant cells, or immune-modulating cells with distinct roles in the tumor microenvironment. It is also unknown whether their presence varies across patients or correlates with clinical features or immune microenvironmental interactions. Methods: We performed integrated single-cell 5′ RNA sequencing and B cell receptor (BCR) repertoire sequencing on 7 newly diagnosed FL 1-3A fresh-frozen lymph node biopsies (Dasari et al., Blood Cancer Discov 2025). Using a high-confidence approach that incorporates heavy and light chain pairing, V and J gene usage, and CDR3 nucleotide junction identity, we classified B cells into two groups: (1) cells belonging to the dominant expanded clone (FL B cells) and (2) all other B cells with detectable and non-expanded BCRs (ITB). These clonal identities were then integrated with single-cell transcriptomic profiles using Seurat and Immcantation pipelines. Differential expression analysis was conducted using edgeR to compare FL B cells to ITB cells across all FL biopsies. Results: After quality filtering, a total of 44,637 single cells were retained across all 7 samples (minimum per sample: 3,847 cells), including malignant and non-malignant B cells as well as diverse tumor microenvironment populations. The proportion of B cells assigned to the dominant FL clone varied substantially across tumors, ranging from 22% to 89%, indicating that some tumors contain a substantial population of ITB cells. In one tumor, we identified two expanded malignant B cell clones comprising 40.6% and 34.4% of all cells. These clones exhibited non-overlapping V and J gene usage (IGHV3-11/IGHJ5 vs. IGHV3-72/IGHJ3), distinct CDR3 junctions, and different immunoglobulin isotypes (IGHG1 vs. IGM). BCL2 was expressed in most cells within both clones (Clone 1: 87.8%, Clone 2: 67.0%), with mean log-normalized expression levels of 1.79 and 0.95, respectively. Both clones formed distinct UMAP clusters and showed transcriptional divergence from ITBs, with 453 differentially expressed genes (DEGs) identified for Clone 1 and 98 for Clone 2 (adjusted p<0.01, |log₂FC|>1, min.pct≥0.25). Additionally, 748 DEGs were identified between the two clones. This evidence of clonal divergence within a single FL tumor reinforces the need to reconsider intratumoral homogeneity in FL and suggests that multiple malignant clones may sometimes coexist. When comparing FL clones to ITB cells across tumors we found that FL clones exhibited a transcriptional program enriched for proliferation and cell cycle regulators, including MKI67 (p<0.005), E2F1 (p<0.0001), CDKN2A (p<0.0001), and CDKN2B (p<0.01), as well as elevated expression of AICDA (p<0.01), EZH2 (p<0.0001), and PD-L1 (p<0.05). As expected, the malignant cells also had higher expression of BCL2 (p<0.01) and CD10/MME (p<0.0001). In contrast, ITB cells demonstrated relatively higher expression of SELL (p<0.0001) and TNFRSF14 (p<0.05), which were significantly downregulated in the FL clone. Conclusions: Together, these findings challenge the prevailing notion that FL tumors are defined by a single homogeneous malignant B cell population. Instead, our data support a model in which multiple transcriptionally and functionally distinct B cell populations coexist within tumors, including dominant clones and a diverse population of ITB cells. This is, to our knowledge, the first study to comprehensively characterize these populations using integrated, lineage-aware single-cell multiomic profiling across FL patients. These insights may ultimately inform new strategies for understanding tumor evolution, immune evasion, and therapeutic targeting in follicular lymphoma.

Abstract Description Signal regulatory protein-α (SIRPα) is a key member of the “do-not-eat-me” signaling pathway and is expressed on intratumoral monocytes/macrophages in lymphoma. In this study, we found that SIRPα is also expressed by a subset of malignant B cells in lymphoma. SIRPα+ B cells have a memory phenotype with germinal center (GC) dark-zone features, and RNA velocity analysis showed that SIRPα+ B cells exhibit higher spliced transcripts confirming a more mature phenotype. SIRPα expression is substantially upregulated by treating lymphoma B cells with LPS or CpG and SIRPα+ B cells showed increased proliferation when compared to SIRPα- B cells, which was associated with reduced cellular autophagy. SIRPα+ B cells exhibited increased expression of BAFF signaling pathway molecules as the likely mechanism of increased proliferation and survival. Furthermore, in co-culture, SIRPα+ lymphoma B cells suppressed T-cell function and blocking SIRPa/CD47 signaling reversed T cell inhibition. Clinically, we found that a greater number of SIRPα+ B cells present in the diagnostic tumor biopsies of FL patients was associated with a significantly poorer overall survival. We conclude that SIRPα+ lymphoma B cells are a proliferative component of indolent lymphomas, are associated with T-cell suppression, and increased prevalence of these cells is associated with a poor prognosis. Targeting SIRPα+ B cells in lymphoma may therefore be a promising therapeutic strategy. Topic Categories Tumor Immunology: Cellular Responses and Tumor Microevironment (TIME)

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 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.