Abstract:Diffuse large B-cell lymphoma (DLBCL) is the most common aggressive B-cell lymphoma with high clinical and biological heterogeneity. Despite current effective immunochemotherapy, up to 40% of patients do not respond or develop refractory disease. DLBCL is characterized by two major cell-of-origin (COO) subtypes, germinal center B cell-like (GCB) and activated B cell-like (ABC), with recent research uncovering additional molecular subgroups based on genomic alterations or tumor microenvironment (TME) features. However, current classifiers use bulk or deconvolution profiling which blurs intra-tumor heterogeneity, hindering the discovery of true functional cell states underlying clinical diversity. Here, we constructed a multimodal single-cell atlas of 62 DLBCL cases, integrating single-cell transcriptomes, immune repertoire and genetic features, to decrypt recurrent B cell states and TME ecosystems while identifying clinically relevant prognostic biomarkers (Fig.1). We performed simultaneous 5'-end single-cell RNA, B-cell receptor (BCR) and T-cell receptor sequencing on 63 DLBCL lymph node biopsies, recovering 109,294 cells after stringent quality controls. These biopsies were obtained from the real-world, clinically-annotated, multicentric CeVi collection and included 43 patient samples at diagnosis (ABC = 17, GCB = 11, Unclassified = 6, Not Evaluated = 9; 21 of which did not achieve event-free survival at 24 months (EFS24) after R-chemotherapy regimens), 11 relapsed/refractory DLBCL and 8 DLBCL transformed from indolent lymphoma. We distinguished malignant B-cells from TME using canonical markers and BCR sequence, exploring heterogeneity in both components. We developed an original bioinformatic approach to define malignant B cell transcriptional archetypes shared across patients. To assess whether subclonal genetics are associated with transcriptional states, we inferred copy number variations (CNV) and studied BCR intraclonal variations of malignant cells. We also integrated TME data (39,674 cells) to quantify subpopulations and clonally expanded T-cell subsets, and performed robust statistical analyses to identify clinically relevant features linked to patient outcomes. To circumvent interpatient heterogeneity in malignant B cells, we defined transcriptional cellular archetypes conserved across patients by combining 3 unsupervised methods: rank-biased overlap meta-clustering, variational autoencoder and canonical correlation analysis integration. We identified 5 cellular archetypes (Arch.1-5): Arch.1 was characterized by plasma cell identity, Arch.2 by MYC and NFKB pathways, Arch.3 by MYC pathway only, Arch.4 by quiescence markers, and Arch.5 by GC markers. Several archetypes co-existed in each tumor sample, and archetype identity was superimposed to the transcriptional profile conferred by COO molecular class as an additional layer of transcriptional heterogeneity. Cell state-specific CNV and BCR subclones were very rare, suggesting that subclonal genetic changes were not the main drivers of malignant B cell plasticity and archetypes. Kaplan-Meier analysis showed a significant negative association between proportion of Arch.4 cells and continuous EFS in our diagnosis cohort. Unsupervised analysis of archetype composition further distinguished ABC patients who achieved EFS24 (enriched in Arch.3 cells) from those with a particularly poor prognosis (enriched in Arch.4 cells). Integrated analysis of the TME revealed three major ecosystems, including one enriched for effector and cytotoxic T cells. In diagnosis samples, pseudo-temporal trajectory analysis revealed the enrichment of clonally-expanded hyperactivated CD8+ T cells expressing canonical exhaustion markers in patients achieving EFS24, suggesting a positive link between the extent of CD8+ T cell differentiation and response to R-CHOP based regimens. Shedding new light on DLBCL biological and clinical heterogeneity, our comprehensive single-cell atlas of DLBCL uncovered recurrent tumor cell states and TME features crucial for stratifying non-responsive ABC patients and guiding treatment decisions. Ongoing single-cell spatial transcriptomics analyses will be used to further validate our findings and elucidate archetype-TME interactions that may drive malignant B cell plasticity, with the potential to uncover novel therapeutic targets.

Integrative genomic analysis identifies unique immune environments associated with immunotherapy response in diffuse large B cell lymphoma

Integrative genomic analysis identifies unique immune environments associated with immunotherapy response in diffuse large B cell lymphoma

Single-cell spatial analysis of tumor immune architecture in diffuse large B-cell lymphoma

Integrative analyses reveal outcome-associated and targetable molecular partnerships between TP53, BRD4, TNFRSF10B, and CDKN1A in diffuse large B-cell lymphoma

Immune Profiling of Diagnostic DLBCL Biopsies Dramatically Improves upon Cell-of-Origin Risk Stratification

Integrated Genomic and Transcriptomic Analyses of Diffuse Large B-Cell Lymphoma With Multiple Abnormal Immunologic Markers

Immune cell constitution in the tumor microenvironment predicts the outcome in diffuse large B-cell lymphoma

Abstract PR06: Comprehensive characterization of the non-tumor microenvironment of relapsed/refractory large B-cell lymphoma identifies patients with greatest benefit from CD19 CAR T-cell therapy

Immune-based therapies in diffuse large B-cell lymphoma

Immune Profiling and Quantitative Analysis Decipher the Clinical Role of Immune-Checkpoint Expression in the Tumor Immune Microenvironment of DLBCL

Spatial transcriptomics unveils immune cellular ecosystems associated with patient survival in diffuse large B-cell lymphoma

Genomic drivers of large B-cell lymphoma resistance to CD19 CAR-T therapy

Tumor-infiltrating normal B cells revealed by immunoglobulin repertoire clonotype analysis are highly prognostic and crucial for antitumor immune responses in DLBCL

Data from Tumor-Infiltrating Normal B Cells Revealed by Immunoglobulin Repertoire Clonotype Analysis Are Highly Prognostic and Crucial for Antitumor Immune Responses in DLBCL

Immune microenvironment-related gene mapping predicts immunochemotherapy response and prognosis in diffuse large B-cell lymphoma

Gene set-based identification of two immune subtypes of diffuse large B cell lymphoma for guiding immune checkpoint blocking therapy

A Single Cell Atlas of Diffuse Large B Cell Lymphomas Reveals Distinct Cellular States Predictive of Outcomes

Baseline Immune State and T-cell Clonal Kinetics are Associated with Durable Response to CAR-T Therapy in Large B-cell Lymphoma

Genetic landscapes and curative effect of CAR T-cell immunotherapy in patients with relapsed or refractory DLBCL.

Hitting back at lymphoma: How do modern diagnostics identify high-risk diffuse large B-cell lymphoma subsets and alter treatment?

B cells and tertiary lymphoid structures promote immunotherapy response