Over the last few years, several proposals for a genetic classification of DLBCL have been published. Hereunder, the main contributions will be summarized and discussed based on the technical approach used.
4.1. Whole-Exome Sequencing (WES)-Based Studies
In 2017, Reddy and coworkers
[49] reported on the whole-exome sequencing (WES) of 1001 FF DLBCLs and 400 paired germline DNAs. They found 150 driver genes to be recurrently mutated. The 60 top genes frequently exhibited a pattern of either predominant missense and/or copy number gains consistent with an oncogene or truncating mutations and/or copy number losses consistent with a tumor suppressor gene. When the mutational pattern was matched with the COO, 20 genes were differentially mutated between the two groups, including
EZH2, SGK1, GNA13, SOCS1, STAT6, and
TNFRSF14, which were mutated in GCB tumors, and
ETV6, MYD88, PIM1 and
TBL1XR1, which were mutated in ABC tumors. Interestingly,
MLL2 mutations were associated with those of
MYC, while
TP53 mutations occurred in a mutually exclusive fashion with
KLHL6. CRISPR screening revealed that knockout of
EBF1,
IRF4,
CARD11,
MYD88, and
IKBKB was selectively lethal in ABC DLBCL cell lines, as was knockout of
ZBTB7A, XPO1, TGFBR2, and
PTPN6 in the GCB lines. On prognostic grounds,
MYC mutations were strongly associated with poorer survival, as were mutations in
CD79B and
ZFAT. Mutations in
NF1 and
SGK1 were associated with more favorable survival. Furthermore, in ABC DLBCLs, genetic alterations in
KLHL14, BTG1, PAX5, and
CDKN2A were associated with significantly poorer survival, while those in
CREBBP were associated with favorable outcomes. In the GCB-DLBCL group, genetic alterations in
NFKBIA and
NCOR1 were associated with poorer prognosis, while alterations in
EZH2, MYD88, and
ARID5B were all associated with a significantly better prognosis. The authors developed a multivariate supervised learning approach for defining the association of survival with combinations of genetic markers (150 genetic driver genes) and gene expression markers (cell of origin, MYC, and BCL2). This led to the proposal of a three-subgroup molecular risk model that was found to outperform all existing predictors (i.e., COO, MYC/BCL2 DE, and IPI). However, the recent application of this model to 499 DLBCLs by Bolen et al.
[50] did not provide independent validation. This might reflect the technical differences between the two studies (WES of FF samples by Reddy et al. vs. targeted NGS of DNA extracted from FFPE biopsies by Bolen et al.).
Two studies published in 2018 proposed a molecular subclassification of DLBCLs that had potential prognostic and therapeutic implications
[51][52]. They both were based on WES and copy-number analysis of a large series of FF DLBCLs (304 and 574, respectively)
[51][52].
Chapuy et al.
[51] described five clusters characterized by different genetic lesions that were capable of identifying subgroups within the COO categories showing different behaviors.
Most cases included in clusters (Cs) 1 and 5 were classified as ABC. However, they showed important differences on molecular and prognostic grounds. C1 cases were thought to derive from marginal-zone B-cells, as they showed a stable mutational pattern, structural variants (SVs) of BCL6, and mutations of genes involved in the NOTCH2 and NF-kB pathways (NOTCH2, SPEN, BCL10, TNFAIP3, and FAS). Besides the multiple genetic lesions of genes involved in immune escape (BM2, CD70, FAS, PD-L1, PD-L2), these C1 cases carried MYD88 mutations which were non-L265P, unlike what was observed in the cases included in C5. Notably, C1 cases had a rather favorable course and revealed potential therapeutic targets related to NOTCH2 and BCL6 signaling and immune evasion mechanisms. C5 tumors, which behaved more aggressively than the C1 ones, showed mutations of MYD88L265P, CD79B, PIM1, TBL1XR1, GRHPR, and BTG1, SV of 18q, and activation of the NF-kB pathway. In addition, they carried ongoing mutations, being at least in part under the effect of AID. Potential targets for C5 cases corresponded to BCR/TLR signaling and BCL2.
Cs 3 and 4 were significantly enriched in GCB cases but were characterized by different genetic lesions and responses to chemoimmunotherapy. The majority of DLBCLs in C3 harbored BCL2 mutations with concordant SVs. They also exhibited frequent mutations in chromatin modifiers, KMT2D, CREBBP, and EZH2, and increased transcriptional abundance of EZH2 targets by gene set enrichment analisys (GSEA). These tumors also had alterations in the B-cell transcription factors MEF2B and IRF8, and indirect modifiers of BCR and PI3K signaling (TNFSF14(HVEM), HCNV1, and GNA13). In addition, C3 tumors had two alternative mechanisms of inactivating PTEN: focal 10q23.31/PTEN loss and predominantly truncating PTEN mutations, events that play a role in the process of lymphomagenesis. C4 DLBCLs were characterized by mutations in four linker and four core histone genes, multiple immune evasion molecules (CD83, CD58, and CD70), BCR/PI3K signaling intermediates (RHOA, GNA13, and SGK1), NF-kB modifiers (CARD11, NFKBIE, and NFKBIA), and RAS/JAK/STAT pathway members (BRAF and STAT3). Comparison of the C3 and C4 genetic signatures further revealed that these GCB-DLBCLs utilized distinct mechanisms to perturb common pathways such as PI3K signaling. In contrast to C3 DLBCLs, C4 tumors rarely exhibited PTEN alterations but harbored more frequent RHOA mutations. In addition, C4 DLBCLs rarely exhibited BCL2 alterations and had higher mutational density. The distinct genetic features of C3 and C4 GCB-DLBCLs led Chapuy et al. to suggest specific targeted therapies including inhibition of BCL2, PI3K, and the epigenetic modifiers EZH2 and CREBBP in C3 GCB tumors, and JAK/STAT and BRAF/MEK1 blockade in C4 GCB-DLBCLs. Last but not least, C3 cases had a far worse prognosis.
C2 DLBCLs harbored frequent biallelic inactivation of
TP53 by mutations and
17p copy loss. In addition, they often exhibited copy loss of
9p21.13/CDKN2A and
13q14.2/RB1, perturbing chromosomal stability and cell cycle. C2 tumors also had significantly more driver somatic copy number alterations (SCNAs) and a higher proportion of genome doubling events. This cluster included both GCB- and ABC-DLBCLs, as did prior DLBCL cohorts with
TP53 mutations in targeted analyses
[53]. Prognostically significant SCNAs, including
13q31.31/miR-17-92 copy gain and
1q42.12 copy loss, were also more common in these DLBCLs, which were characterized by a rather unfavorable prognosis.
A further cluster, termed 0, was also detected, which apparently lacked significant genetic alterations. However, as the C0 group consisted almost exclusively of T-cell rich/histiocyte-rich B-cell lymphomas, the obtained results might have been largely influenced by the small number of neoplastic cells.
The authors further evaluated BCL2 and MYC alterations. Tumors with cooccurring BCL2 and MYC SVs were significantly more frequent in C3 DLBCLs.
Importantly, the coordinate genetic signatures reported by Chapuy et al. predicted outcomes independent of IPI which could suggest new combination treatment strategies and, more broadly, provide a roadmap for actionable DLBCL classification
[51].
By their integrated approach, Schmitz et al.
[52] identified four prominent genetic subtypes among 574 DLBCLs which they termed MCD (based on the co-occurrence of
MYD88L265P and
CD79B mutations), BN2 (based on
BCL6 fusions and
NOTCH2 mutations), N1 (based on
NOTCH1 mutations), and EZB (based on
EZH2 mutations and
BCL2 translocations). Interestingly, Schmitz and co-workers enriched their series with unclassified DLBCLs. The latter turned out to frequently carry mutations affecting
SPEN and
NOTH2 as well as
BCL6 fusions. ABC cases were enriched in
MYD88L265P and
CD79B or
NOTCH1 mutations, with the two conditions being mutually exclusive. GCB tumors showed the co-occurrence of
EZH2 mutations and
BCL2 translocations. The MCD and N1 subtypes were dominated by ABC cases, while EZB included mostly GCB tumors, and BN2 had contributions from all GEP subgroups. Overall, about 45% of the samples were classified into the genetically pure subtypes of DLBCL.
The MCD subtype displayed 82% of cases carrying MYD88L265P or CD79B aberrations (mutation or amplification), with 42% bearing both abnormalities. The MCD subtype showed a frequent gain or amplification of SPIB, encoding a transcription factor that, with IRF4, defines the ABC phenotype and promotes plasmacytic differentiation. Known tumor suppressors in MCD include CDKN2A, ETV6, BTG1, and BTG2, and putative tumor suppressors include TOX, SETD1B, FOXC1, TBL1XR1, and KLHL14. The tumor suppressor TP53 was mutated significantly less often in MCD as compared to other subtypes. Immune editing appeared prominent in MCD genomes, with 76% acquiring a mutation or deletion of HLA-A, HLA-B, or HLA-C and 30% acquiring truncating mutations targeting CD58.
BN2 was dominated by NOTCH pathway aberrations, with 73% acquiring a NOTCH2 mutation or amplification, SPEN mutation, or mutation in DTX1, a NOTCH target gene. BCL6 fusion, the other BN2 hallmark, occurred in 73% of cases. BCL6 fusions were enriched in cases with NOTCH2, SPEN, or DTX1 lesions to a significantly greater extent in BN2 than in non-BN2 cases. Genetic aberrations (mutations or amplifications) targeting regulators of the NF-kB pathway were a prominent feature of BN2. These more often affected TNFAIP3, PRKCB, and BCL10. Other likely gain-of-function events included mutations targeting cyclin D3 and CXCR5, whereas inactivating lesions targeting the immune regulator CD70 suggested immune escape.
N1 was characterized by NOTCH1 mutations and aberrations targeting transcriptional regulators of B-cell differentiation (IRF4, ID3, and BCOR), which may contribute to its plasmacytic phenotype. TNFAIP3 mutations in N1 could reinforce this phenotype by fostering NF-KB-induced IRF4 expression.
EZB was enriched for most of the genetic events previously ascribed to GCB-DLBCL, including BCL2 translocation, EZH2 mutation, and REL amplification, as well as inactivation of the tumor suppressors TNFRSF14, CREBBP, EP300, and KMT2D. The germinal-center homing pathway involving S1PR2 and GNA1314 was disrupted in 38% of EZB cases. JAK-STAT signaling was promoted in about half cases by a STAT6 mutation or amplification or by a mutation or deletion targeting SOCS1. PI3K target of rapamycin signaling turned out to be activated in 23% of cases by MTOR mutations or the amplification of MIR17HG. Immune editing was of interest in EZB genomes since 39% acquired lesions in the major histocompatibility complex class II pathway genes CIITA and HLA-DMA.
The four subtypes differed significantly in PFS and OS, with the BN2 and EZB subtypes having much more favorable outcomes than the MCD and N1 subtypes. The predicted 5-year OS rates for the MCD, N1, BN2, and EZB subtypes were 26%, 36%, 65%, and 68%, respectively. Within ABC DLBCL, patients with MCD had significantly inferior survival as compared with those with BN2, and patients with either MCD or N1 had significantly inferior survival as compared with patients with ABC tumors that were not genetically classified. Within GCB-DLBCL, there was a trend toward inferior OS among patients with EZB as compared with those with other GCB tumors. The COO subgroups and genetic subtypes independently contributed to survival in a multivariate analysis. Conversely, the IPI score did not vary significantly among the genetic subtypes, but the latter significantly added to IPI. A trend toward increased extranodal involvement (e.g., CNS) was a feature of MCD, which reflected the frequent CD79B and MYD88L265P mutations.
On therapeutic grounds, constitutive BCR signaling activation was most frequent in MCD and least frequent in EZB, but genetic alterations involving the BCR cascade occurred in all genetics subtypes, suggesting that constitutive BCR signaling is a pervasive aspect of DLBCL pathogenesis. BN2 was notably enriched for BCR–NF-kB and IKK regulator aberrations. In addition to NF-kB, survival of DLBCL cells turned out to be promoted by antiapoptotic BCL2 family members, which were targeted by genomic amplification or translocation in 17.4% of cases. As expected, BCL2 mRNA levels were significantly higher in EZB tumors with BCL2 translocations than in other EZB tumors. MCD tumors also had high BCL2 mRNA expression as compared with other cases, a finding due to mechanisms other than translocation or amplification.