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.
4.2. Targeted NGS and Bioinformatic-Based Studies
Lacy et al.
[54] applied a 293-gene chip to DNA extracted from FFPE tissue samples by using a Covaris LE220. The authors sequenced a large, unselected cohort consisting of 928 DLBCL patients all treated with R-CHOP and provided with full clinical follow-up. Bernoulli mixture-model clustering was applied, and the resulting subtypes analyzed in relation to their clinical characteristics and outcomes. Five molecular subtypes were resolved, termed MYD88, BCL2, SOCS1/SGK1, TET2/SGK1, and NOTCH2, along with an unclassified group. The subtypes characterized by genetic alterations of
BCL2,
NOTCH2, and
MYD88 recapitulated the above-mentioned studies showing good, intermediate, and poor prognosis, respectively. The SOCS1/SGK1 subtype showed biological overlap with primary mediastinal B-cell lymphoma and conferred excellent prognosis. Although not identified as a distinct cluster,
NOTCH1 mutation was associated with poor prognosis. The impact of
TP53 mutation varied with genomic subtypes, conferring no effect in the NOTCH2 subtype and poor prognosis in the MYD88 subtype. The results obtained by Lacy et al. are summarized in , where they are also compared with the subtypes reported by Chapuy et al.
[51], and Schmitz et al.
[52].
Table 2.
Molecular subtypes of DLBCL according to Lacy et al. in comparison with those of Chapuy et al. and Schmitz et al.
BCL2/BCL6 rearrangement data. Schmitz’s cohort was used as training set, while those from Chapuy et al. and Ennishi et al. were used for validation. Wright et al. developed a model that is summarized in , which also includes information on potential therapeutic targets related to the genetic subtype of DLBCL.
Table 3.
Implications of genetic subtypes of DLBCL for therapy (from Wright et al., modified).
Ennishi et al.
[55] performed an integrative genomic and transcriptomic analysis of DLBCL using a British Columbia population-based registry. They uncovered recurrent biallelic
TMEM30A loss-of-function mutations which were associated with a favorable outcome and were uniquely observed in DLBCL. Using
TMEM30A-knockout systems, increased accumulation of chemotherapy drugs was observed in
TMEM30A-knockout cell lines and
TMEM30A-mutated primary cells, accounting for the improved treatment outcome. Furthermore, they found increased tumor-associated macrophages and an enhanced effect of anti-CD47 blockade limiting tumor growth in
TMEM30A-knockout models. By contrast, Ennishi et al. showed that TMEM30A loss-of-function increased B-cell signaling following antigen stimulation—a mechanism conferring selective advantage during B-cell lymphoma development. These findings suggested intrinsic and extrinsic vulnerabilities of cancer cells that can be therapeutically exploited.
Finally, Wright et al.
[56] developed an algorithm to determine the probability of a patient’s lymphoma belonging to one of seven genetic subtypes based on its genetic features. This represented a probabilistic classification tool (LymphGen) using any combination of mutational, copy number, and
5. Open Issues and Perspectives
The paper of Lacy et al.
[54] was accompanied by a commentary from Morin and Scott
[57], who concluded that comprehensive sequencing of a larger number of tumors with the combination of whole-genome and transcriptome sequencing is warranted to develop a new molecular taxonomy which may be concretely translated into clinical benefits. In fact, between 7.5% and 55% of the cases reported by Chapuy et al., Schmitz et al., and Lacy et al. did not fit into any of the major genetic categories they identified
[51][53][58]. The fact that the genomic studies hitherto reported show a certain variability in terms of results may depend on different factors, such as the size of the analyzed cohort or heterogeneity of the techniques used (e.g., FF vs. FFPE tissue, whole exome vs. targeted sequencing, and the statistical approach applied), but also on the actual heterogeneity of the lesions occurring in these tumors. For instance, divergent evolution within the same biopsy, which corresponded to different morphologic, phenotypic, and COO features
[59], has been reported. Although the distinct components had a common clonal origin and shared the bulk of genetic aberrations, each revealed private mutations, in keeping with the above-mentioned morpho-phenotypic and molecular differences.
The heterogeneity of genetic lesions is much greater than was thought until a couple of years ago. This has been highlighted by liquid biopsy (LB)
[59][60][61]. By ultradeep sequencing of the cell-free circulating tumoral DNA (cfDNA) released by neoplastic cells undergoing apoptosis, it has been shown that the global mutational landscape of DLBCL is indeed wider than that observed in diagnostic biopsies, which means that different mutations can occur at different anatomic sites. Once a standardized methodology is developed and the cost per test is reduced, LB can represent a real-time noninvasive tool for disease monitoring. In fact, patients achieving early molecular response (a 2-log decrease of ctDNA after one cycle of standard chemoimmunotherapy) and major molecular response (a 2.5-log decrease after two cycles) show a significantly superior outcome at 24 months independently of IPI and interim positron emission tomography. Conversely, among treatment-resistant subjects, new mutations are acquired in cfDNA, marking resistant clones selected during the clonal evolution.
The continuous development of sequencing and bioinformatic techniques will allow us to achieve the long searched-for goal of using customized therapies based on the molecular characteristics of each individual tumor. Some approaches do appear to be more easily and cheaply applicable to daily life. Nevertheless, the more comprehensive the bioinformatic approach is, the higher the likelihood of overcoming today’s standard chemoimmunotherapy and designing chemotherapy-free protocols capable of curing most (if not all) patients, with minimal or no toxic effects.