Abstract A049: Spatially resolved microenvironment profiling reveals distinct follicular lymphoma archetypes

Abstract

Recent single-cell transcriptomic studies of follicular lymphoma (FL) have revealed heterogeneity among malignant B (MB) cells and selected components of the tumor microenvironment (TME), but have been constrained by relatively small sample sizes and a lack of spatial context. We performed comprehensive cellular profiling and spatial dissection of over five million cells by integrating single-nucleus RNA sequencing (snRNA-seq) from 182 samples (FL, n = 167; control, n = 15) with single-cell spatial transcriptomics (ST) from 386 samples (FL, n = 365; control, n = 21). After stringent quality control, 1,289,073 and 3,968,427 cells from snRNA-seq and ST data, respectively, were analyzed. Unsupervised clustering of snRNA-seq data identified 58 TME cell subtypes. To extend this classification to the ST dataset, we integrated snRNA-seq and ST data by co-embedding them into a shared principal component space, enabling robust batch correction and cell-type label transfer across modalities. Using this supervised annotation framework, ST data recapitulated 54 TME cell subtypes. MB, proliferating MB (MBprolif), and healthy B (HB) cells were classified using sample-level unsupervised clustering in ST. Macrophages exhibited previously uncharacterized heterogeneity in transcriptomic profiles and spatial distribution. We identified CXCL13-expressing macrophage subtypes with high intra- and peri-follicular abundance, termed follicular macrophages (FMs). Non-negative matrix factorization of intrafollicular cell fractions identified five TME archetypes: FM, follicular stroma, two follicular T, and HB archetypes. The FM archetype was enriched for intrafollicular macrophages, including FMs, and showed a significantly higher proportion of grade 3A/3B cases (53.4%) compared to other archetypes (≤16.5%). It was also associated with a higher frequency of high tumor burden requiring treatment initiation rather than a watch-and-wait (WW) approach. Among patients managed by WW, those with the FM archetype showed significantly shorter event-free survival (EFS) compared to those with non-FM archetypes, whereas this association was not observed among treated patients. The percentage of intrafollicular macrophages (IFMP) showed excellent predictive performance for the FM archetype (area under the curve: 0.90). IFMP calculated in a published multiplex spatial protein imaging dataset from an independent cohort of 242 FL patients validated the associations between high IFMP and both pathology grade 3A/3B and shorter EFS in WW patients. Consistently, spatial neighborhood and cell–cell communication analyses revealed co-localization of FMs and MBprolif cells, forming a distinct niche with activated signaling mediated by FM-derived CXCL13, APRIL, and BAFF. Our integrative multi-modal analysis defines previously unrecognized cellular heterogeneity and spatial architecture in FL and highlights macrophage-enriched follicular ecosystems associated with aggressive clinicopathological features, providing a framework for improved biological understanding and patient management.