Novak, Anne J.

Abstract Background: Waldenström Macroglobulinemia (WM) is an indolent B-cell lymphoplasmacytic lymphoma characterized by heterogeneous clinical behavior and variable progression patterns. While MYD88L265Pand CXCR4 mutations are frequently encountered, the broader genomic landscape underpinning disease evolution from smoldering WM (SWM) to active treatment naïve WM (TNWM), post-treatment remission WM (PTWM) and relapsed/refractory WM (RRWM) remains incompletely understood. To address this gap, we studied the mutational profile across the WM disease spectrum. Methods We performed whole exome sequencing (WES) on sorted (CD19 and/or CD138) bone marrow aspirates of patients across the spectrum of WM. Paired germline sequencing was performed using peripheral blood mononuclear cells (PBMCs). Unique molecular identifier-tagged reads were deduplicated using fgbio. WES data were processed through the Sarek pipeline. Any valid mutation from either caller, present on the index or in COSMIC was considered. Copy number variations (CNVs) were inferred via CNVkit and summarized at the gene and arm level (≥50% directional change). Mean tumor and germline coverages were 174X and 74X, respectively. Fisher's exact test was used for enrichment and hypergeometric test for overlap analyses. Results A total of 78 high-quality WM samples [15 SWM, 44 TNWM, 13 RRWM, 6 PTWM in remission (all with paired TNWM samples)] were sequenced. The mutational burden varied across disease states (median SWM=1.48 m/Mb, TNWM= 1.12 m/Mb, RRWM = 1.59 m/Mb; p=0.04). Density plots revealed bimodal VAF distributions, consistent with clonal and subclonal architectures in WM. Overall, 77% harbored a MYD88L265P mutation by WES. Among MYD88^WT^ patients, AS-PCR identified MYD88L265P in 4 (25%) not detected by WES, raising the overall mutation rate to 83%, plus one with a non-L265P variant. Low VAF (<10%) MYD88L265P mutations were seen in 6% of matched PBMCs. MYD88 copy number gain was seen in 2 (17%) MYD88WT TNWM patients. CXCR4 mutations were detected in 38% by WES (18% each frameshift and nonsense). Two (4%) CXCR4WT samples were positive by PCR, giving an overall rate of 40%. CXCR4 mutations were more frequent in RRWM (54%) vs TNWM (41%) and SWM (27%), though not statistically significant. Del(6q) and Del(17p) were noted in 16% (n=12) and 7% (n=5) of patient samples. All del(17p) events were noted in TNWM samples and in one paired post-treatment sample. Overall rate of del(17p) and TP53 mutation was 14% (n=11). RRWM was enriched in ARID1A alterations, a gene implicated in chromatin remodeling and plasma cell identity formation [38% (n=5) in RRWM vs. 11% (n=6) in TNWM and SWM], compared to TNWM and SWM p=0.03]. Alterations in NF-kB signaling gene MAP3K14 were more frequent in RRWM (31%, n=4) compared to the SWM/TNWM cohort [(7%), n=4; p=0.036). We noted an enrichment in the TP53 mutation rate in RRMW [n=3, 23%, OR 7.7] compared to TNWM and SWM cohorts [n=2, 4%), p=0.04]. Mutations in FUBP1 (p=0.03), an upstream regulator of MYC, and CARD11 (p=.043), a B-cell signaling scaffold protein were also more frequent in RRWM. Notably, none of the SWM patients harbored an ARID1A alteration (p=0.057), suggesting acquisition during progression. CNV analysis revealed that MAP3K14 amplifications were frequent in RRWM (n=3, 23%) compared to SWM/TNWM [n=2 (4%), p=0.043]. RRWM samples showed increased CD79B mutations (15%) and significant overlap with the C5/ABC diffuse large B-cell lymphoma (DLBCL) genomic cluster (Chapuy et al., Nat. Med. 2018; p=0.007). This overlap remained significant after excluding MYD88L265P (p=0.019), driven by CARD11, CD79B, and IGLL5, suggesting that a subset of RRWM may acquire an ABC DLBCL-like mutational profile. These findings warrant further investigation using transcriptional profiling. Conclusions This study provides a genomic map of WM evolution from precursor to relapse. Increasing somatic mutations and CNVs support stepwise progression. Enrichment of genes involved in chromatin remodeling and plasma cell identity formation (ARID1A), NF-kB or BCR signaling (CARD11, MAP3K14), and DNA damage repair (TP53) were noted at relapse, all of which have been shown to drive aggressive disease biology. Emergence of an ABC-like genotype in a proportion of RRWM indicates convergent evolution with aggressive lymphomas. Transcriptional studies are warranted to explore this further.

Abstract Despite a characteristic indolent course, a substantial subset of follicular lymphoma (FL) patients has an early relapse with a poor outcome. Thus far, efforts to identify factors that predict survival have been unsuccessful. However, we and others have demonstrated the prognostic relevance of CD4+ T cells in the tumor microenvironment (TME) of FL and developed a bio-clinical risk model (called BioFLIPI) that offers improved risk stratification (Mondello et al, BCJ 2021). However, the mechanisms defining different BioFLIPI scores remain to be fully understood. To identify the programs associated with the BioFLIPI scores, we first interrogated the RNA-seq of purified B cells from 67 newly diagnosed FL1-3A (BioFLIPI 1 n=21; BioFLIPI 2-3 n=36; BioFLIPI 4 n=10). We found 1,156 upregulated genes and downregulation of 930 genes in BioFLIPI 4 vs 1, including well-known immune molecules such as CD70, IL10RA, CD69, and CCR6 (FDR<0.05). Gene set enrichment analysis (GSEA) using the MSigDB confirmed a negative enrichment in multiple immune response and interferon gamma signatures and revealed striking positive enrichment for cell proliferation and antigen presentation pathways in BioFLIPI 4 compared to BioFLIPI 1 tumors. The same pathways were dysregulated in BioFLIPI 2-3 but a lower magnitude. To investigate the somatic mutations associated with different BioFLIPI scores, we performed WES of the same 67 FLs and further expanded this cohort for a total of 90 patients (BioFLIPI 1 n=26; BioFLIPI 2-3 n=51; BioFLIPI 4 n=13). We found enrichment for mutations in several regulators of the germinal center in BioFLIPI 4 and 2-3, such as MEF2B (p=0.029), FOXO1 (p=0.038) and IRF8 (p=0.035), suggesting the cooperation of multiple oncogenic events promoting lymphomagenesis. In contrast, BioFLIPI 1 tumors harbored a higher rate of mutations in TNFRSF14 gene, known to enhance anti-lymphoma immune response. To investigate whether BioFLIPI modulates immune signatures of the microenvironment, we performed TME deconvolution by applying CIBERSORTx to bulk RNA-seq of these 67 FLs. We found that FLs with BioFLIPI 4 and 2-3 were enriched for CD4+ Tcm exh (p=0.02), Tem exh (p=0.04), Tfhexh (p=0.07) and NK cells (p=0.01) compared to BioFLIPI 1. In line with this data, CyTOF performed on a subset of these patients (BioFLIPI 1 n=12; BioFLIPI 2-3 n=25; BioFLIPI 4 n=6) confirmed progressive increase of CD4+ Tcm exh (p=0.006), Tem exh (p=0.009) and Tfh exh (p=0.025) cells from BioFLIPI 1 to BioFLIPI 4. Single cell RNA/BCR-seq of 21 FLs (BioFLIPI 1 n=6; BioFLIPI 2-3 n=12; BioFLIPI 4 n=3) corroborated a higher abundance of Tem exh, cytotoxic (CTL) Tem exh, Tfhexh and NK cells in BioFLIPI 4 and 2-3. Conversely, CD8+ naïve T cells were prevalent in BioFLIPI 1 tumors. Differential gene expression in T cell clusters from BioFLIPI 4 vs BioFLIPI 1 FL revealed a marked upregulation of multiple immune genes involved in antigen presentation (e.g., HLA-DRB, HLA-DQA, HLA-DPA/B) and immune response (e.g., ETS1, TNFSF8, IL7R). GSEA showed a positive enrichment for antigen presentation signatures and negative enrichment for TNFα signaling and Oxphos metabolism. Similar findings were observed in BioFLIPI 2-3 FLs. Using CellChat analysis, we found a stronger interaction between malignant B cells and Tem exh and CTL Tem exh cells in BioFLIPI 4 and 2-3 vs BioFLIPI 1 tumors. This was linked to an increased strenght of MHCI-CD8 interaction and concordant decrease of MHCII-CD4. Notably, the suppressive BTLA-TNFRSF14 signaling axis was present in BioFLIPI 4 and 2-3 but not in BioFLIPI 1 FLs. This was accompanied by a decreased signaling through TNF/LTA-TNFRSF1B, CD80/CD86-CD28/CTLA4. Linking the BioFLIPI-specific immunogenicity to lymphomagenesis, we found that FL patients with BioFLIPI 4 had two-three folds higher risk of transformation to an aggressive lymphoma (12% vs 7%) or death (21% vs 3%) than those with BioFLIPI 1, suggesting that this immune dysfunction favors the expansion of a dominant tumor clone and promotes a more aggressive form of FL. Collectively, our findings demonstrate that BioFLIPI scores in FL reflect distinct patterns of tumor immunogenicity that shape the TME. Notably, high BioFLIPI scores are paradoxically associated with enhanced antigen presentation alongside a suppressed immune response. This suggests a sustained antigen-specific T cell response that may drive immune exhaustion, promote tumor progression, and ultimately lead to histologic transformation

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 Description The efficacy of immunotherapy for splenic marginal zone lymphoma (SMZL) remains limited due to the strong immune suppression present in the tumor microenvironment (TME). Tregs play a key role in maintaining TME homeostasis by suppressing effector T cells and malignant B cells. Previously, we identified two Treg subsets CD26+ and CD161+ in SMZL. CD161+ Tregs were more effective at suppressing both CD8+T cells and malignant B cells, whereas CD26+Tregs primarily inhibited CD8+T cells. Patients with higher levels of CD26+ Tregs and fewer CD161+ Tregs tend to have worse outcomes. In this study, we found that cytokine-inducible SH2-containing protein (CISH), a known negative regulator of immune activation, is highly expressed in CD26+Tregs compared to CD161+ Tregs. Co-culturing B cells enhanced the expression of CISH and CD26 in CD4+CD25- Treg cells, knockdown of CISH in Tregs reduced CD26 expression and their immune suppression activity of CD8+T cells, while CRISPR knockout of CISH in Karpas 299 cells led to higher expression CD70, a molecule associated with less immunesuppression. Our findings suggest that increased CISH expression in Tregs could result in increased immune suppression and unfavorable outcome in SMZL, suggesting CISH is a potential therapeutic tagrget in SMZL by dempening immunosuppressive property of the Treg cells hence boosting the antitumor immunity. Topic Categories Tumor Immunology: Cellular Responses and Tumor Microevironment (TIME)

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 Although follicular lymphoma (FL) typically follows an indolent course, patients with FL who experience early events, such as transformation or progression, have increased risk of death related to lymphoma. The FL24Cx is an algorithm based on a 45-target gene expression profiling (GEP) assay, which was developed and trained using 265 formalin-fixed, paraffin-embedded tissue samples on a reliable platform to predict, at the time of diagnosis, whether a patient will experience an event within 24 months. The modeling also confirmed and relied upon previously reported synergy between immune response (IR) gene expression signatures IR1 and IR2. Once locked, the 5-factor logistic regression FL24Cx model was independently validated in a retrospectively assessed cohort of 232 patients from 2 immunochemotherapy-treated arms of SWOG Cancer Research Network S0016 phase 3 clinical trial, in which it assigned 169 patients to the low-risk group with 29 events before 24 months (17.2%) and 63 patients to the high-risk group with 24 events before 24 months (38.1%). The relative risk of an event within 24 months after registration among patients who were classified into the high-risk group relative to patients who were classified into the low-risk group was 2.2 (95% confidence interval, 1.41 to 3.51). An up-front GEP biomarker, such as the FL24Cx, rigorously validated in a clinical laboratory and with a clinically relevant turnaround time, could identify and steer enrollment of patients at high risk for early events in clinical trials, thus enabling timely interpretation of such trials and increasing the pace of innovation.

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.