Mondello, Patrizia

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 Chimeric antigen receptor (CAR) T-cell therapy achieves high response rates in relapsed/refractory aggressive B-cell non-Hodgkin lymphoma, yet only 40% of patients achieve durable remission. Identifying biological programs linked to long-term response remains critical. As metabolic fitness underlies T-cell persistence and effector function, we profiled metabolic pathways associated with clinical outcomes following commercial CD19 CAR T-cell therapy. Single-cell RNA sequencing of peripheral blood mononuclear cells collected at baseline (BL), peak CAR-T expansion (PK), and one-month post-infusion (M1) was analyzed using CellRanger v7.0.1, immunopipe. Seurat v4.3.0 was applied for unsupervised clustering to delineate cell subsets based on top differentially expressed genes. Patients were categorized as durable complete remission (CR ≥6 months; n = 16), primary refractory (PD1; n = 4), or relapse after initial response (PD2; n = 12). Gene-set enrichment analysis defined metabolic pathway activity across T-cell, monocyte, dendritic cell (DC), and natural killer (NK) subsets. Oxidative phosphorylation (OXPHOS) emerged as the dominant metabolic program distinguishing clinical outcomes. At BL and PK, OXPHOS was consistently enriched in PD1 relative to CR across T-cell, monocyte, and DC subsets, suggesting early oxidative activation in non-responders. By M1, this pattern inverted, with higher OXPHOS activity in CR. PD2 largely paralleled PD1, but several subsets including CD8 T central memory, classical monocytes TGFβ, intermediate monocytes CD38, monocytic myeloid-derived suppressor cell (mMDSC) HIF1A, and mMDSC SIRPA, showed OXPHOS enrichment in CR at BL and/or PK that persisted through M1. In contrast, DC and NK subsets exhibited the opposite pattern: OXPHOS was enriched in CR at PK (conventional DC 2 (cDC2), plasmacytoid DC (pDC), NK, proliferating NK), but shifted toward enrichment in PD2 at M1 (cDC2, NK, NK CD56bright). Glycolysis (GLY) followed a similar trajectory in monocytes and NK cells, with enrichment in PD1/PD2 at BL and PK, followed by enrichment in CR at M1. In T-cells, PD1 maintained GLY enrichment from BL through PK, with no significant differences observed at M1. The inositol-phosphate metabolism pathway showed a more static pattern and was consistently enriched in multiple PD2 effector and memory T-cell subsets without reversal at M1. Taken together, these data reveal distinct, lineage-specific metabolic signatures that differentiate durable remission from early and late progression. Dynamic OXPHOS and GLY programming, characterized by lower activity early and enhanced activity at M1 in CR, may reflect adaptive metabolic programming that supports sustained antitumor immunity. Immune-metabolic profiling may therefore serve as a biomarker of response and highlight actionable metabolic pathways to enhance CAR T-cell durability. Citation Format: Melinda S.Y. Tan, Panwen Wang, Patrizia Mondello, Jacqueline Turner, Andre de Menezes Silva Corraes, Chen Wu, Zuoyi Shao, Kevin Regan, Ma Audrey, Arushi Khurana, Nora N. Benanni, Yucai Wang, Paul Hampel, Jonas Paludo, Saad J. Kenderian, Urshila Durani, Patrick B. Johnston, Jose Caetano Villasboas, Stephen M. Ansell, Ying Li, Haidong Dong, Hu Zeng, Yi Lin. Metabolic pathway signatures defining response to CD19 CAR T-Cell therapy in aggressive B-cell non-hodgkin lymphoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 3278.

Abstract Introduction: Therapeutic outcomes in B-cell non-Hodgkin’s lymphoma (B-NHL) have substantially improved with the introduction of chimeric antigen receptor T-cell (CAR-T) therapy. Yet, persistent immune dysfunction and underlying metabolic dysregulation remain major obstacles to achieving durable clinical responses in B-NHL. We hypothesize that mitochondrial pathways play a critical role in mediating anti-tumor immunity during CAR-T therapy. Methods: Peripheral blood mononuclear cells were collected from healthy controls (CNTRL, n=5) and patients with advanced stage lymphoma (LYM) who received FDA approved CAR-T (n=32). Samples were collected before lymphodepletion (BL), at peak CAR-T expansion (PK) and one-month post-infusion (M1). Single-cell RNA sequencing was used to interrogate differential mitochondrial gene expression. Samples were obtained from LYM patients with complete remission for >6 months (CR), primary refractory (PD1), or relapsed (PD2) disease. Results: Compared to CNTRLs, LYM patients exhibited reduced expression of cytochrome c oxidase. At BL, CR patients demonstrated broader mitochondrial differential gene expression across all T-cell subsets relative to PD1 and PD2 and showed higher expression of MT-ATP-6, MT-ATP-8, and MT-CO3. Following CAR-T infusion, global mitochondrial gene expression increased across patients corresponding with expansion of CD8 effector T cells (Tem tumor circulating) and activated peripheral memory T cells (Tpm). Certain mitochondrial genes, such as MT-CO1, were durably expressed across LYM patients at BL, PK, and M1 in peripheral and central memory T cells (Tpm, Tcm). At PK, cytochrome c oxidase and ATP synthase genes were upregulated in CR Tpm and Tcm. At M1, MT-ATP-6, MT-ATP-8, MT-CO3 were significantly upregulated in CD8 memory precursor effector T-cells (MPECs) of CR patients compared to both PD1 and PD2. While CR and PD2 patients appeared clinically responsive at M1, CR samples maintained significantly higher MT-ATP-8 expression across T cell subsets (CD8 MPECS, Tpm, Tcm) compared to PD2 (p values 0.008; <0.0001; <0.0001). Conclusions: We report a temporal shift in immunometabolism with CAR-T therapy in LYM. Clinical responders exhibited sustained mitochondrial activity characterized by persistent expression of ATP synthase in CD8 MPECs and cytochrome c oxidase in Tpm and Tcm. In sum, mitochondrial programming may distinguish effective versus dysfunctional T cells and identify targets for therapeutic development. Citation Format: Jacqueline Turner, Panwen Wang, Patrizia Mondello, Melinda Tan, Christoph Schaefers, Chen Wu, Andre de Menezes Silva Corraes, Kevin Regan, Zuoyi Shao, Ma Audrey, Arushi Khurana, Nora N. Benanni, Yucai Wang, Paul Hampel, Saad J. Kenderian, Jonas Paludo, Urshila Durani, Patrick B. Johnston, Jose Caetano Villasboas, Stephen M. Ansell, Haidong Dong, Ying Li, Zeng Hu, Yi Lin. Mitochondrial transcriptome profiles associated with clinical responses to CAR-T therapy in aggressive lymphoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 1213.

Abstract Waldenström macroglobulinemia (WM) is characterized by recurrent MYD88 and CXCR4 mutations, whose prognostic value in patients treated with chemoimmunotherapy remains unclear. Moreover, the typically prolonged progression-free survival (PFS) correlates inconsistently with overall survival (OS), underscoring the importance of examining other surrogates. Progression of disease within 24 months (POD24), an established early end point, delineates functionally high-risk patients in other indolent lymphomas. This international study evaluated 253 patients receiving frontline fixed-duration bendamustine-rituximab (BR), a common chemoimmunotherapy for WM. At a median follow-up of 5.9 years, 5-year PFS and OS were 65% and 87%, respectively; 5-year PFS was similar between MYD88L265P (90%) and MYD88wild-type (WT) subcohorts (64% each; P = .4). Among 89 patients with known CXCR4 status, the subcohort with CXCR4mutation (28%) had shorter PFS (median 3.3 vs 8.8 years; hazard ratio [HR], 2.8; P = .0036) and OS (HR, 2.6; P = .036) compared with CXCR4WT. POD24 occurred in 11.5% of patients who demonstrated inferior subsequent OS (5-year OS, 71% vs 86%; HR, 3.1; P = .005) and higher mortality (standardized mortality ratio [SMR], 3.7), unlike the non-POD24 group, whose mortality was comparable to the matched general population (SMR, 1.1). In conclusion, BR is effective irrespective of the MYD88 status, but CXCR4 mutations and POD24 portend worse outcomes. Patients without POD24 represent a cohort with distinctly favorable outcome.

Abstract Chimeric antigen receptor (CAR) T-cell therapy has demonstrated remarkable response rates in relapsed/refractory aggressive B-cell non-Hodgkin's lymphoma. However, durable remissions occur in only 40% of treated patients, highlighting mechanisms of therapeutic resistance and immune escape. While cytotoxicity is principally executed by CAR T-cells, monocytes shape CAR-T cell efficacy and orchestrate cytokine-driven immunity. In this study, we defined distinct metabolic profiles of T-cell and monocyte subsets based on clinical response. Single-cell RNA sequencing (scRNA-seq) and T-cell receptor sequencing (scTCR-seq) were performed on peripheral blood mononuclear cells at three timepoints: before lymphodepletion (BL), at peak CAR-T expansion (PK), and one-month post-infusion (M1). Data were processed with CellRanger v7.0.1 and analyzed with the immunopipe pipeline (Wang et al., 2025). Seurat v4.3.0 was utilized for unsupervised hierarchical clustering analysis; subsets identified included classical monocytes (CM), intermediate monocytes (IM), non-classical monocytes (NM), and monocytic myeloid derived suppressor cells (mMDSC). Metabolic profiles were compared across patients achieving ≥6-month complete remission (CR, n=16), primary refractory disease (PD1, n=4), or relapsed after a prior response (PD2, n=12). A total of 107035 T cells and 129981 monocytes from 77 serial samples in 32 patients were analyzed. CR patients had a higher T-cell:monocyte ratio than PD1/PD2 at all time points, especially at BL and PK (BL 1.2 vs 0.4, p=0.011; PK 1.3 vs 0.3, p=0.027). Clinical responses were compared with metabolic pathway analysis and identified oxidative phosphorylation (OXPHOS) as the dominant pathway. OXPHOS was consistently enriched in PD cohorts across T-cell subsets, particularly in PD1 vs CR at BL (mucosal-associated invariant T-cell (MAIT) and T peripheral memory (Tpm)), and at PK (CD4/CD8 memory precursor effector cells (MPEC) and activated MPECs (MPECa), CD4 T-helper 1 proliferating (Th1p), CD4 Th2 activated (Th2a), regulatory T-cell (Treg), CD8 naive senescent (CD8NS), CD8 T effector memory (Tem) tumor recirculating (TR), CD8 Tpm, CD8 T resident memory (Trm) TR and MAIT Tpm). By M1, OXPHOS enrichment in CD4/CD8MPECa, CD4Th2a, CD8NS, CD8TemTR reversed, now predominant in CR although CD8Tpm and CD8TrmTR remained enriched in PD1. In the monocyte subsets at BL, PD2 showed OXPHOS enrichment in CD163+ monocytes (mono), IM IL1, mMDSC SIRPA, mMDSC TGFβ, and CM IL6 subsets compared to PD1 and CR, while PD1 showed OXPHOS enrichment in IM CD38, CM TGFβ, and mMDSC HIF1A vs CR and PD2. This persisted at PK. By M1, the pattern reversed and mirrored that of T-cells, with CR patients showing OXPHOS enrichment in CD163+mono, CM IL6, CM TGFβ, IM CD38, mMDSC HIF1A, mMDSC SIRPA and mMDSC TGFβ, and additionally in CM EMT, CM p53, NM, NM CD36 vs PD1/PD2. Glycolysis/gluconeogenesis (GLY) was also enriched across subsets. In T-cells, PD1 vs CR showed GLY enrichment in CD4Th2a, CD8Tpm, CD8Tpm proliferating and MAIT Tpm at PK. In monocytes, GLY was enriched in PD1/2 vs CR at BL (CD163+mono, CM IL6) and PK (CD163+mono, IM CD18, mMDSC SIRPA). By M1, GLY metabolism was enriched in CR vs PD1/PD2 in CD163+mono, CM IL6, mMDSC SIRPA. Additional GLY enriched subsets at CR vs PD1 included CM, CM p53, EMT, IM CD38, IM IL1, mMDSC HIF1A, NM, NM CD36. Inositol phosphate metabolism (IP) was predominantly enriched in PD2 vs CR at BL (CD8NS, CD8Tcm, CD8Tpm activated), PK (CD4Th2, CD8NS, CD8Tem TR, CD8Tpm proliferating, γδT Tem) and M1 (CD4/CD8 MPECa, CD4 Th1/Th2, CD8NS, CD8Tpm, CD8 Tpm proliferating, Treg and MAIT Tpm). Notably, Treg was the only T-cell subset enriched for IP at PK in CR vs PD1/PD2. Additionally, the TCA cycle was enriched in PD1 Tregs and Th2 cells at BL and PK vs CR/PD2. A higher T-cell:monocyte ratio in CR vs PD1/PD2 supports prior evidence that elevated pre-lymphodepletion lymphocyte:monocyte ratios predict better outcomes. We identified distinct metabolic signatures in T-cells and monocytes between CR and PD1/PD2. OXPHOS enrichment in T-cells and monocytes at M1 in CR patients may reflect oxidative programming as a hallmark of sustained anti-tumor immunity. Conversely, GLY enrichment in PD1 T-cells at PK may reflect short-lived effector states lacking persistence. Overall, these findings highlight metabolic profiling as a promising predictive biomarker and therapeutic programming target to enhance CAR T-cell outcomes.

Abstract Introduction: Outcomes for FL have markedly improved, and the 10-year (yr) overall survival (OS) for patients (pts) diagnosed in the rituximab era is about 80%. Historically, FL is the leading cause of death in the first decade after diagnosis (dx); however, modern outcomes beyond 10 yrs after dx remain unknown. Here we report long-term outcomes of FL from a large prospective cohort of pts with FL from the rituximab era. Methods: Pts with biopsy-proven FL 1-3A at dx were prospectively enrolled within 9 months of dx from 2002-2015 in the MER. Outcomes of interest were OS, event-free survival (EFS), histologic transformation (HT), and cause of death. HT was defined as subsequent biopsy-proven FL 3B, diffuse large B cell lymphoma, or high-grade B cell lymphoma. Analyses were performed using Cox models and KM curves; cause of death was evaluated using a competing risk approach. A subset analysis was performed in younger (18-59 yrs) vs older (60+ yrs) pts. A landmark analysis was conducted in pts alive at 10 yrs. Results: 1,155 pts with FL 1-3A were enrolled in the MER from 2002-2015. Median age at dx was 60 yrs (IQR 51-69). Clinical characteristics were: 67% stage III/IV, 12% grade 3A, 4% ECOG PS 2-4, and 19% elevated LDH. Initial management was 35% observation (obs), 12% R-monotherapy (R-mono), 39% immunochemotherapy (IC; 43% R-CHOP, 28% R-CVP, 27% BR), and 14% other. FLIPI scores were 40% low, 33% intermediate (int), and 27% high; FLIPI24 risk categories were 49% low, 25% int, and 26% high. At a median follow-up of 14 yrs (IQR 11-18 yrs), 751 events, 148 HTs, and 402 deaths were observed. EFS estimates (95% CI) were: 10-yr: 40% (37-43), 15-yr: 31% (28-35), and 18-yr: 27% (24-31). 15 yr & 18 yr EFS by initial txs were 22% (18-27) & 18% (14-25) for obs, 30% (23-40) & 18% (11-30) for R-mono, and 40% (35-46) & 38% (33-45) for IC. The incidence (95% CI) of HT was 12% (10-14) at 10 yrs, 13% (11-16) at 15 yrs, and 14% (12-17) at 18 yrs. OS (95% CI) for the entire cohort was 77% (75-80%) at 10 yrs, 63% (60-66) at 15 yrs, and 55% (51-59) at 18 yrs. The incidence of non-lymphoma-related death (non-LRD) surpassed lymphoma-related death (LRD) after 12 yrs in pts with known cause of death. Estimates (95% CI) for LRD were 11% (10-13) at 10 yrs, 14% (12-17) at 15 yrs, and 16% (14-19) at 18 yrs, while estimates for non-LRD were at 8% (7-10) at 10 yrs, 17% (15-20) at 15 yrs, and 22% (19-26) at 18 yrs. Among pts diagnosed before age 60, LRD remained the predominant cause of mortality into the second decade following dx, with incidences (95% CI) of LRD vs. non-LRD of 7% (5-9) vs 3% (2-4) at 10 yrs, 10% (8-13) vs 5% (3-8) at 15 yrs, and 13% (9-17) vs 6% (4-10) at 18 yrs, respectively. In patients age >60 non-LRD was the predominant cause of mortality in the second decade. Incidence of LRD vs. non-LRD in those aged ≥60 at dx were 15% (13-19) vs 13% (10-16) at 10 yrs, 18% (15-22) vs 28% (24-32) at 15 yrs, and 20% (16-24) vs 35% (31-41) at 18 yrs, respectively. In the landmark analysis of pts alive at 10 yrs from dx (n=720), median age was 69 yrs (IQR 60-77). 63 pts (9%) had HT prior to landmark. The median number of prior treatments was 1 (IQR 1-2, range 1-8). Risk of death (95% CI) was 19% (15-22) in the next 5 yrs (i.e. 15 yrs from dx) and 29% (24-33) in the next 8 yrs (18 yrs from dx). The risk of LRD vs non-LRD was 4% (3-6) vs 12% (9-15) at 5 yrs and 7% (3-9) vs 18% (10-23) at 8 yrs.Conclusions: Long-term outcomes for pts with FL remain excellent into the second decade from dx. HT risk decreases significantly after 10 yrs, though late HTs still occur. At 15 yrs from dx, 40% of pts remained alive and event-free after 1L IC. Nearly a quarter of pts observed at dx remained treatment-free. The most common cause of death after 10 yrs from dx was from non-lymphoma causes, notably in older pts aging into their 70s and 80s. While younger pts have a lower overall risk of LRD compared to older pts, FL represents the leading cause of death in this population during the second decade following dx. These results provide critical insights into the natural history of FL diagnosed in the modern era.

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 A subset of follicular lymphoma 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. However, little is known on the molecular mechanisms whereby IRF4 controls the crosstalk between B and T cells. To this aim, we leveraged Cγ1cre;irf4fl/fl mice that yield conditional deletion of irf4 (cKO) in germinal center B cells (GCB). After 10 days of immunization, cKO GCB showed significant increase of HLA-DR and CD40 expression, a skewed LZ:DZ ratio towards increased DZ, and suppression of plasma cells compared to WT. Upon irf4 deletion, we also found a significant increase of total CD4 and TFH cells and a decrease in Treg cells. To test if this effect was antigen-dependent, we cocultured sorted GCB from NP-OVA-immunized cKO and WT mice with CSFE-labeled OVA-specific OT-II CD4 T cells at a ratio of 2:1 for 5 days. We observed increased antigen specific proliferation of CD4 T cells with cKO GCB compared to WT. We then performed seahorse assay and mass spectrometry of the same GCB to define if IRF4-dependent B cell metabolic reprogramming influences immune responses. GCB cKO had lower glycolytic ECAR and increased intracellular glucose and non-essential amino acids than WT cells, suggesting a metabolic shift from glycolysis to the TCA cycle. Thus, IRF4 controls immune and metabolic signaling with functional effects on the B-T cell crosstalk. Funding Sources NCI K08CA279652, LRF Clinical Investigator Award (CDA 1020588), the Gerstner Family Career Development Award, the ASH Junior Faculty Scholar Award. Topic Categories Tumor Immunology: Cellular Responses and Tumor Microevironment (TIME)