Diffuse large B-cell lymphoma (DLBCL), the most common type of lymphoma, arises from various pathogenic mechanisms including gene translocations and fusions. Detection of gene rearrangements using fluorescence in situ hybridization (FISH) is a standard practice for DLBCL diagnosis and guides treatment decisions. High-throughput chromosome conformation capture (Hi-C) DNA sequencing is a next-generation sequencing-based technology to identify genome-wide rearrangements using a single assay. In this study, Hi-C sequencing was performed using FFPE tissues from 159 patients with DLBCL, and identified 746 cancer genes at or proximal to the rearrangement breakpoints with a total of 1903 occurrences. Focusing on clinically important rearrangements, Hi-C detected 102 rearrangements of MYC, BCL2, and/or BCL6 in 92 patients and revealed the fusion partners, including 25 rearrangements missed by FISH. Causes of FISH negative results included FISH-cryptic breakpoints, complex or faint FISH signals, low percentage of FISH positivity, and issues related to tissue fixation. Moreover, in 20 patients (12.6%), Hi-C sequencing detected 22 rearrangements characteristic of other lymphoma types, including CCND1 rearrangements that could lead to reclassification as mantle cell lymphoma. Survival analysis for genome-wide rearrangements using machine learning models identified MYC, PD-L1, CCND1, BCL2, NTRK1/PRCC, RRAS, and FANCE rearrangements with significant prognostic effects in the DLBCL cohort. In conclusion, Hi-C sequencing detects gene rearrangements crucial for diagnosis in an unbiased and molecular manner and showed high sensitivity and specificity in our study. These advantages of Hi-C sequencing offer help to improve the workflow of clinical pathology laboratories, diagnostic precision, and treatment of large B-cell lymphoma.
Proinflammatory GSDMD activation in live macrophages and DLBCL cells marks cognate interactions and better prognosis
Abstract Pyroptosis is a form of programmed cell death characterized by the cleavage of the gasdermin (GSDM) family proteins that form pores in the plasma membrane, cell rupture, and the release of proinflammatory cytokines. In this study, we performed immunohistochemistry for cleaved gasdermin D (GSDMD), gasdermin E (GSDME) N-terminus, and gasdermin C (GSDMC) in 2 different cohorts of diffuse large B-cell lymphoma (DLBCL) and analyzed the impact of GSDM expression on prognosis and immunity. The results showed frequent cleaved GSDMD (N-terminal) expression. Only cytoplasmic GSDMD N-terminal expression correlated with significantly better patient survival in the 2 cohorts. In contrast, GSDME was mainly expressed in the vascular endothelium, and correlated with significantly adverse prognostic effect. Correlating with the multiplex fluorescent immunohistochemistry results, we found that cytoplasmic GSDMD N-terminal expression was associated with increased CD38+ (activated) M1 macrophages in both cohorts, cognate interactions between live DLBCL cells and activated M1 macrophages (and T cells), and lower PD-1/PD-L1 expression in the analyzed cases. In contrast, T-cell pyroptosis, lymphoma cell-resistance to cell death, and phagocytosis by M2 macrophages were observed in tissues with predominantly nuclear GSDMD N-terminal expression. Bulk gene expression profiling and deconvolution analysis revealed associations of cytoplasmic GSDMD N-terminal expression with downregulation of “Don’t eat me”–signaling genes, upregulation of many RNA genes, decreased frequency of “inflammatory” lymphoma microenvironment subtype, increased frequencies of prognostically favorable cell states and ecotypes, and decreased frequency of T-cell exhaustion state. In summary, this study showed distinct cellular and subcellular patterns of 3 GSDM proteins and their associated immune response phenotypes and prognostic effects, with implications for novel therapeutic strategies for B-cell lymphoma.
Prognostic gene expression and microRNA profiling signatures and genetic alterations in primary testicular diffuse large B-cell lymphoma
Abstract Primary testicular diffuse large B-cell lymphoma (PT-DLBCL) is a rare and aggressive lymphoma with molecular heterogeneity not well characterize. In this study, we performed next-generation sequencing analysis for a large number of DNA and RNA samples from patients with PT-DLBCL. DNA sequencing analysis identified ≥ 3 chromosomes with copy number variations (CNVs) and microsatellite instability as prognostic biomarkers, rather than MYD88 mutations and genetic subtypes. Remarkably, targeted RNA-seq analysis in 195 patients revealed that TP53 mutations with a ≥ 40% variant allele frequency had significantly adverse prognostic impact, and that a 150-gene expression signature subdivided PT-DLBCL into two distinct clusters, termed as testicular lymphoma tumor (TLT) and microenvironment (ME) subtypes. The TLT subtype featured upregulation of genes involved in B-cell receptor signaling, cell cycle, DNA damage and repair, higher frequencies of CNVs and MYD88 mutations, elder ages, larger tumor sizes, and significantly poorer survival. Genomic microRNA profiling analysis identified significantly differentially expressed microRNAs between 113 PT-DLBCL and 180 systemic DLBCL patients, and further subdivided the PT-DLBCL cohort by microRNA signatures. The subcohort with upregulation of 16 microRNAs associated with PT-DLBCL and testicular tissue expression had significantly better survival. This study revealed characteristic genetic, gene expression, and microRNA profiles and heterogeneity in PT-DLBCL.
DNA mismatch repair defect and intratumor heterogeneous deficiency differently impact immune responses in diffuse large B-cell lymphoma
Diffuse Large B-Cell Lymphoma Has a Low Frequency of dMMR and High Frequencies of DNA Mismatch Repair Protein High Expression Associated with Lower T-Cell Infiltration
Background: DNA mismatch repair (MMR) deficiency (dMMR), which leads to genomic and microsatellite instability, is a biomarker that predicts response to immunotherapies and has variable prognostic effects for chemotherapies in solid tumors. The dMMR can be detected using gene sequencing or immunohistochemistry for four essential MMR proteins: MSH6, MSH2, MLH1, and PMS2. Diffuse large B-cell lymphoma (DLBCL) is the most common aggressive lymphoma. Previous studies have suggested the role of MMR gene variants in DLBCL lymphomagenesis, yet the prognostic role of dMMR in DLBCL has not been well studied. Methods: We performed targeted next-generation sequencing and immunohistochemistry for MSH6, MSH2, MLH1, and PMS2 in a large cohort of DLBCL patients treated with standard chemoimmunotherapy. Gene expression profiling (GEP) was performed using the Affymetrix GeneChip Human Genome HG-U133 Plus 2.0 microarray (data in GSE31312). Fluorescent multiplex immunohistochemistry (mIHC) was performed using MultiOmyx multiplexing immunofluorescence staining protocols and antibodies against 13 immune markers. We investigated the frequencies of dMRR and MMR protein expression and correlated the mutation status and expression levels with patient survival and DLBCL biology, including the number of mutated genes, GEP data, previously analyzed DLBCL biomarkers, and immune cell immunophenotypes. Results: MMR gene mutations and loss of MMR protein expression were infrequent in DLBCL (Figure A) and did not show a significant prognostic impact. MMR proteins were commonly expressed in DLBCL samples and the germinal centers of reactive tonsil controls (Figure B), with higher mean and median percentages of tumor cells expressing MSH6 and MLH1 proteins than in the solid tumor samples used for comparison. High expression of MMR proteins was associated with Ki-67, MYC, and p53 overexpression as determined by immunohistochemistry. GEP analysis identified significantly upregulated genes involved in the mitotic cell cycle and DNA metabolism and prominent downregulated immune gene signatures in DLBCL with high MMR protein expression. Fluorescent mIHC confirmed the associations with decreased T cell abundance in MSH6/MSH2/MLH1/PMS2 highly expressing DLBCL (Figure C) independent of p53 and MYC expression status, whereas MSH6/ MLH1 mutations were associated with increased T cell frequencies. High versus low expression of MSH6, MLH1, and PMS2 showed significant unfavorable prognostic effects in DLBCL when the optimal cutoffs were used for high expression. However, when the median percentages were used as cutoffs, the prognostic effects were not significant; nonetheless, MSH6 and PMS2 high expression showed significant adverse prognostic effects in the MYC¯ DLBCL subset. Interestingly, when we used Ki-67, MYC + or BCL2 + percentages in DLBCL assessed by immunohistochemistry to subtract MMR protein percentages and then correlated the differences to survival, we found significant favorable prognostic impact by higher percentage expression differences between MMR proteins and BCL2, or by intermediately higher differences between MMR proteins and MYC. For MSH2/MLH1/PMS2, the higher percentages compared with MYC/BCL2 showed no significant correlations with T cell infiltration levels, whereas the higher differences between MSH6 and MYC/BCL2 were associated with significantly lower T cell frequencies (Figure D). Conclusions: This study revealed that high expression of MMR proteins is a common feature of DLBCL associated with lower tumor-infiltrating T cells, MYC and p53 overexpression, and context-dependent prognostic effects. In contrast, dMMR and loss of MMR proteins are infrequent and have no significant prognostic effects in DLBCL treated with standard chemoimmunotherapy. These results may have implications for understanding DLBCL biology and the low efficacy of PD-1 blockade immunotherapy in DLBCL.
Tumor-Infiltrating Normal B Cells Revealed by Immunoglobulin Repertoire Clonotype Analysis Are Highly Prognostic and Crucial for Antitumor Immune Responses in DLBCL
Abstract Purpose: Tumor-infiltrating B lymphocytes (TIL-B) have demonstrated prognostic and predictive significance in solid cancers. In this study, we aimed to distinguish TIL-Bs from malignant B-cells in diffuse large B-cell lymphoma (DLBCL) and determine the clinical and biological significance. Experimental Design: A total of 269 patients with de novo DLBCL from the International DLBCL R-CHOP Consortium Program were studied. Ultra-deep sequencing of the immunoglobulin genes was performed to determine B-cell clonotypes. The frequencies and numbers of TIL-B clonotypes in individual repertoires were correlated with patient survival, gene expression profiling (GEP) data, and frequencies of DLBCL-infiltrating immune cells quantified by fluorescent multiplex IHC at single-cell resolution. Results: TIL-B abundance, evaluated by frequencies of normal B-cell clonotypes in the immunoglobulin repertoires, remarkably showed positive associations with significantly better survival of patients in our sequenced cohorts. DLBCLs with high versus low TIL-B abundance displayed distinct GEP signatures, increased pre-memory B-cell state and naïve CD4 T-cell state fractions, and higher CD4+ T-cell infiltration. TIL-B frequency, as a new biomarker in DLBCL, outperformed the germinal center (GC) B-cell–like/activated B-cell–like classification and TIL-T frequency. The identified TIL-B–high GEP signature, including genes upregulated during T-dependent B-cell activation and those highly expressed in normal GC B cells and T cells, showed significant favorable prognostic effects in several external validation cohorts. Conclusions: TIL-B frequency is a significant prognostic factor in DLBCL and plays a crucial role in antitumor immune responses. This study provides novel insights into the prognostic determinants in DLBCL and TIL-B functions with important therapeutic implications.
Genetic Subtyping and Phenotypic Characterization of the Immune Microenvironment and MYC/BCL2 Double Expression Reveal Heterogeneity in Diffuse Large B-cell Lymphoma
Abstract Purpose: Diffuse large B-cell lymphoma (DLBCL) is molecularly and clinically heterogeneous, and can be subtyped according to genetic alterations, cell-of-origin, or microenvironmental signatures using high-throughput genomic data at the DNA or RNA level. Although high-throughput proteomic profiling has not been available for DLBCL subtyping, MYC/BCL2 protein double expression (DE) is an established prognostic biomarker in DLBCL. The purpose of this study is to reveal the relative prognostic roles of DLBCL genetic, phenotypic, and microenvironmental biomarkers. Experimental Design: We performed targeted next-generation sequencing; IHC for MYC, BCL2, and FN1; and fluorescent multiplex IHC for microenvironmental markers in a large cohort of DLBCL. We performed correlative and prognostic analyses within and across DLBCL genetic subtypes and MYC/BCL2 double expressors. Results: We found that MYC/BCL2 double-high-expression (DhE) had significant adverse prognostic impact within the EZB genetic subtype and LymphGen-unclassified DLBCL cases but not within MCD and ST2 genetic subtypes. Conversely, KMT2D mutations significantly stratified DhE but not non-DhE DLBCL. T-cell infiltration showed favorable prognostic effects within BN2, MCD, and DhE but unfavorable effects within ST2 and LymphGen-unclassified cases. FN1 and PD-1–high expression had significant adverse prognostic effects within multiple DLBCL genetic/phenotypic subgroups. The prognostic effects of DhE and immune biomarkers within DLBCL genetic subtypes were independent although DhE and high Ki-67 were significantly associated with lower T-cell infiltration in LymphGen-unclassified cases. Conclusions: Together, these results demonstrated independent and additive prognostic effects of phenotypic MYC/BCL2 and microenvironment biomarkers and genetic subtyping in DLBCL prognostication, important for improving DLBCL classification and identifying prognostic determinants and therapeutic targets.
Genomic complexity is associated with epigenetic regulator mutations and poor prognosis in diffuse large B-cell lymphoma
Abstract Diffuse large B-cell lymphoma (DLBCL) is the major type of aggressive B-cell lymphoma. High-grade B-cell lymphoma (HGBCL) with MYC/BCL2 double-hit (DH) represents a distinct entity with dismal prognosis after standard immunochemotherapy in the current WHO lymphoma classification. However, whether TP53 mutation synergizes with MYC abnormalities (MYC rearrangement and/or Myc protein overexpression) contributing to HGBCL-like biology and prognosis is not well investigated. In this study, patients with DLBCL with MYC/TP53 abnormalities demonstrated poor clinical outcome, high-grade morphology, and distinct gene expression signatures. To identify more effective therapies for this distinctive DLBCL subset, novel MYC/TP53/BCL-2–targeted agents were investigated in DLBCL cells with MYC/TP53 dual alterations or HGBCL-MYC/BCL2-DH. A BET inhibitor INCB057643 effectively inhibited cell viability and induced apoptosis in DLBCL/HGBCL cells regardless of MYC/BCL2/TP53 status. Combining INCB057643 with a MDM2-p53 inhibitor DS3032b significantly enhanced the cytotoxic effects in HGBCL-DH without TP53 mutation, while combining with the BCL-2 inhibitor venetoclax displayed potent therapeutic synergy in DLBCL/HGBCL cells with and without concurrent TP53 mutation. Reverse-phase protein arrays revealed the synergistic molecular actions by INCB057643, DS3032b and venetoclax to induce cell-cycle arrest and apoptosis and to inhibit AKT/MEK/ERK/mTOR pathways, as well as potential drug resistance mechanisms mediated by upregulation of Mcl-1 and RAS/RAF/MEK/ERK pathways. In summary, these findings support subclassification of DLBCL/HGBCL with dual MYC/TP53 alterations, which demonstrates distinct pathobiologic features and dismal survival with standard therapy, therefore requiring additional targeted therapies. Implications: The clinical and pharmacologic studies suggest recognizing DLBCL with concomitant TP53 mutation and MYC abnormalities as a distinctive entity necessary for precision oncology practice. Visual Overview: http://mcr.aacrjournals.org/content/molcanres/19/2/249/F1.large.jpg.
A refined cell-of-origin classifier with targeted NGS and artificial intelligence shows robust predictive value in DLBCL
Abstract Diffuse large B-cell lymphoma (DLBCL) is a heterogeneous entity of B-cell lymphoma. Cell-of-origin (COO) classification of DLBCL is required in routine practice by the World Health Organization classification for biological and therapeutic insights. Genetic subtypes uncovered recently are based on distinct genetic alterations in DLBCL, which are different from the COO subtypes defined by gene expression signatures of normal B cells retained in DLBCL. We hypothesize that classifiers incorporating both genome-wide gene-expression and pathogenetic variables can improve the therapeutic significance of DLBCL classification. To develop such refined classifiers, we performed targeted RNA sequencing (RNA-Seq) with a commercially available next-generation sequencing (NGS) platform in a large cohort of 418 DLBCLs. Genetic and transcriptional data obtained by RNA-Seq in a single run were explored by state-of-the-art artificial intelligence (AI) to develop a NGS-COO classifier for COO assignment and NGS survival models for clinical outcome prediction. The NGS-COO model built through applying AI in the training set was robust, showing high concordance with COO classification by either Affymetrix GeneChip microarray or the NanoString Lymph2Cx assay in 2 validation sets. Although the NGS-COO model was not trained for clinical outcome, the activated B-cell–like compared with the germinal-center B-cell–like subtype had significantly poorer survival. The NGS survival models stratified 30% high-risk patients in the validation set with poor survival as in the training set. These results demonstrate that targeted RNA-Seq coupled with AI deep learning techniques provides reproducible, efficient, and affordable assays for clinical application. The clinical grade assays and NGS models integrating both genetic and transcriptional factors developed in this study may eventually support precision medicine in DLBCL.