Unsupervised hierarchical clustering of CNAs of driver genes and its survival analysis in DLBCL and PMBCL
Based on array CGH, Lenz G et al. study previously identified specific CNAs in PMBCL which were different from ABC and GCB of DLBCL[6]. ABC DLBCLs often have CNAs in FOXP1, NFKBIZ, CDKN2A, CDKN2B, INF4a, BCL2, NFATC1 and SPIB, while GCB DLBCLs frequently harbor CNAs in REL, PTEN, MDM2, MIHG1 and ING1. PMBCL often demonstrate CNAs of JAK2 and PD-L2[6]. Using unsupervised hierarchical clustering, we explored the CNA-based pattern of these genes in DLBCL and PMBCL. The result showed that a small group of DLBCL (11.7%, 9/77) was clustered together with PMBCL as Cluster_2, with amplification of JAK2 (100%,9/9) and PD-L2 (77.8%,7/9)(Fig. 1a). This subgroup of DLBCL occurred at the site of cervical lymph node (4 cases), gastrointestinal tract (3 cases), nasal cavity (1 case) and spleen(1 cases) (Fig. 1a, Table 1). The frequency of JAK2 and PD-L2 amplification in the whole cohort of DLBCL were 11.7% (9/77) and 9.1% (7/77), while both of them were 100% in PMBCL(Fig. 1a) (see Additional file 3).
With respect to OS and PFS, there were no significant differences among three CNA-based clusters (Fig. 1b). Of note, we found that Cluster_2 enriched for JAK2/PD-L2 amplification had less frequency of MYD88_L265P mutation (11.1%, 1/9) (Fig. 1a, Table 4), which was supported by the Cancer Genome Atlas data (TCGA, PanCancer Atlas) from cBioPortal[8, 9] (Fig. 1c).
JAK2/PD-L2 amplification identify a distinctive CNA-based pattern of DLBCL similar to that of PMBCL
Since DLBCL with JAK2/PD-L2 amplification had less frequency of MYD88 L265P mutation, our study separated DLBCL patients into three subgroups: DLBCL with JAK2/PD-L2 amplification (DLBCL_JAK2/PD-L2_amp), DLBCL with MYD88 L265P mutation (DLBCL_MYD88_L265P), and DLBCL without JAK2/PD-L2 amplification nor MYD88_L265P mutation (DLBCL_others) (Fig. 2a). One patient who had both JAK2/PD-L2 amplification and MYD88 L265P mutation was put into DLBCL_JAK2/PD-L2_amp subgroup.
Unlike DLBCL_MYD88_L265P and DLBCL_others, DLBCL_JAK2/PD-L2_amp showed a distinctive pattern similar to that of PMBCL, with high frequency of REL and NFKBIZ amplifications, but no amplification of BCL2 and NFATC1 and no deletion of PRDM1 was found (Fig. 2a).
With respect to clinicopathologic characteristics, DLBCL_JAK2/PD-L2_amp tend to be younger than DLBCL_MYD88_L265P (p=0.040) (Table 2). While, Hans model, international prognostic index (IPI) risk category and MYC break-apart didn’t show any significant differences (Table 2).
PD-L1 expression in DLBCL with JAK2/PD-L2 amplification was significantly higher than that in DLBCL with MYD88 L265P mutation
Totally, 32 cases were performed PD-L1 (22C3) IHC detection, including DLBCL_MYD88_L265P (14 cases), DLBCL_JAK2/PD-L2_amp (9 cases), DLBCL_others (5 cases) and PMBCL (4 cases). The result showed that PD-L1 expression in DLBCL_JAK2/PD-L2_amp was significantly higher than that in DLBCL_MYD88_L265P (p=0.011) and DLBCL_others (p=0.026) (Fig. 2b and 2d). While no significant difference was found between DLBCL_JAK2/PD-L2_amp and PMBCL (p=0.876) (Fig. 2b).
JAK2/PD-L2 amplification identify a subgroup of DLBCL with dismal survival similar to that of MYD88 L265P mutation
Trying to explore the survival indication of JAK2/PD-L2 amplification and MYD88 L265P mutation, 50 cases of DLBCLs who received RCHOP-like regiment with or without surgical resection were enrolled to performed cox proportional regression analysis of OS and PFS. The median follow-up time for patients alive was 108.07 months (range, 3.07-191.93 months).
In the univariate analysis, as compared with DLBCL_others, DLBCLs with MYD88 L265P mutation had significantly dismal OS and PFS (p=0.026 and 0.008, respectively), and the same to DLBCLs with JAK2 amplification (p=0.018 and 0.012, respectively). Meanwhile, IPI risk category were significantly associated with OS and PFS (Fig. 2c, Tables 3 and 4).
In the multivariate analysis, IPI risk category and three subgroups of DLBCL were put into analysis. As compared with DLBCL_others, DLBCL with MYD88 L265P mutation still showed poor OS and PFS (p=0.006 and 0.001, respectively), and the same to DLBCL with JAK2 amplification for PFS and OS (p=0.012 and 0.006, respectively). Meanwhile, IPI risk category was still an independent risk predictors for OS and PFS (Fig. 2c, Tables 3 and 4).
Either JAK2/PD-L2 amplification or MYD88 L265P mutation are frequently seen in relapse/refractory DLBCL with PFS less than 2 years.
DLBCL with PFS less than 2 years was defined as primary relapse/refractory cases. Among these cases who treated by RCHOP-like regime, the frequency of JAK2 and PD-L2 amplification were 20%(3/15) and 13.3% (2/15). Meanwhile, the frequency of MYD88 L265P mutation were 33.3% (5/15). DLBCL with either JAK2/PD-L2 amplification or MYD88 L265P accounted for 46.7% (7/15).