- Rapid communication
- Open Access
Mutations of CREBBP and SOCS1 are independent prognostic factors in diffuse large B cell lymphoma: mutational analysis of the SAKK 38/07 prospective clinical trial cohort
Journal of Hematology & Oncologyvolume 10, Article number: 70 (2017)
Recently, the mutational background of diffuse large B cell lymphoma (DLBCL) has been revealed, identifying specific genetic events that drive lymphomagenesis. However, the prognostic value of these mutations remains to be determined. Prognostic biomarkers in DLBCL are urgently needed, since the current clinical parameter-based factors (e.g., International Prognostic Index (IPI)) are insufficient, particularly in identifying patients with poor prognosis who might benefit from alternative treatments.
We investigated the prognostic value of somatic mutations in DLBCL in a clinical trial (NCT00544219) patient cohort homogenously treated with six cycles of rituximab, cyclophosphamide, hydroxydaunorubicin, vincristine, and prednisone (R-CHOP), followed by two cycles of R (R-CHOP-14). The primary endpoint was event-free survival (EFS) at 2 years. Secondary endpoints included progression-free survival (PFS) and overall survival (OS). Targeted high-throughput sequencing (HTS) of tumor genomic DNA was performed on all exons or hotspots of 68 genes frequently mutated in B cell lymphomas. Mutational data was correlated with the endpoints to identify prognostic associations.
Targeted HTS detected somatic mutations in 71/76 (93%) of investigated cases. The most frequently mutated genes were KMT2D, SOCS1, GNA13, and B2M. Survival analysis revealed that CREBBP- and EP300-mutated cases had significantly worse OS, PFS, and EFS. In addition, ATM mutations predicted worse outcomes for all three clinical endpoints in germinal center B cell-like DLBCL. In contrast, SOCS1 mutations were associated with better PFS. On multivariable analysis taken into account IPI and failure to achieve complete remission, CREBBP and EP300 mutations remained significant to predict worse OS, PFS, and EFS.
Targeted mutation analysis of a uniformly treated prospective clinical trial DLBCL cohort identifies tumor-based genetic prognostic markers that could be useful in the clinical management of such patients.
Diffuse large B cell lymphoma (DLBCL) is a heterogeneous and aggressive lymphoid neoplasm, the treatment of which has significantly improved in the last decade with addition of the anti-CD20 monoclonal antibody rituximab (R) to the chemotherapy regimen consisting of cyclophosphamide, hydroxydaunorubicin, vincristine, and prednisone (R-CHOP) . Since then, R-CHOP has become a standard treatment for the vast majority of primary DLBCL cases with cure rates at about 60%. The remaining 40% of DLBCL either relapses after a period of remission or are refractory to the applied first-line therapy. Patients in this group are treated with aggressive salvage regimens supported by autologous stem cell transplantation (ASCT) , but success rates are modest, particularly in primary refractory DLBCL [3, 4]. Therefore, it is important to identify high-risk patients before administration of first-line therapy so that the potentially more aggressive tumors can be treated with alternative regimens .
Risk stratification of DLBCL has relied for more than 20 years on the International Prognostic Index (IPI), which is based on evaluation of multiple clinical parameters . Following the addition of rituximab, IPI was revised (R-IPI)  and, most recently, further enhanced (NCCN-IPI), to better identify high-risk patients . According to NCCN-IPI, the high-risk group has 5-year overall survival (OS) probability of 33% compared to 54% predicted by IPI, although this remains to be confirmed by the datasets of prospective trials.
In addition to patient-based factors, tumor-based prognostic markers have been proposed. It is known that ~10% of DLBCL cases have MYC rearrangements that are strongly associated with worse outcomes, especially if linked to MYC protein overexpression . It has also been shown that the activated B cell (ABC) cell-of-origin (COO) DLBCL subtype has a worse outcome compared to the germinal center B cell (GCB) subtype . Because original classification based on gene expression profiling on a transcriptome level proved to be technically too challenging for implementation in clinical practice, surrogate immunohistochemistry-based algorithms were developed . However, their utility was limited by suboptimal concordance to the gene expression-based gold standard. Recently, a Lymph2Cx assay, which can measure expression of 20 genes and can be also applied on formalin-fixed paraffin-embedded (FFPE) material, was proposed for reproducible classification of DLBCL into COO subgroups . It remains to be seen if this new technology will be widely accepted in lymphoma centers worldwide. In addition to COO classification, immunohistochemical (IHC) studies have identified multiple protein markers, such as CD5, Ki-67, FOXP1, HLA-I, p21, and CD40, that prospectively showed prognostic value for R-CHOP-treated DLBCL [13–17].
During the last decade, substantial progress has been made toward understanding the genetic basis of DLBCL [18–20]. Shared and COO subtype-specific DNA lesions have been identified, converging into several most frequently dysregulated cellular pathways . Based on these findings, a handful of molecular prognostic markers have been identified, among which, TP53, FOXP1, and MYD88 mutations and CDKN2A deletions were associated with inferior outcomes in R-CHOP-treated DLBCL [22–24]. The prognostic role of such molecular markers is likely to increase in the near future as high-throughput sequencing (HTS) enters routine practice in many institutions. However, more studies, particularly prospective analyses, are required to validate the existing molecular prognostic markers and to discover new ones that would add power to the existing prognostication algorithms.
In this study, we employed targeted HTS to identify somatic mutations in tumors of a well-documented prospective clinical cohort consisting of uniformly treated primary DLBCL patients. By correlating gene mutation status to the robust survival data, we aimed to discover new, and validate known, prognostic markers in DLBCL.
The clinical trial SAKK 38/07 (NCT00544219), active between 2007 and 2010, included 138 eligible patients with primary untreated DLBCL to prospectively determine the prognostic value of interim PET/CT scans by standardized treatment and evaluation criteria. The main clinical results have been published previously .
Tissue specimens of patients who consented for additional translational research were used for a subsequent study investigating the prognostic value of phenotypic and genotypic profiles by IHC and fluorescence in situ hybridization (FISH), the results of which have been recently published in this journal .
For the current study, FFPE tissues of 84 primary untreated de novo DLBCL patients with adequate amounts of remaining material were selected. Tumor content was determined by morphological evaluation and was at least 50% in all samples.
Updated clinical data (last follow-up on 31 January 2017) were used for evaluation of the prognostic role of genetic mutations. All patients in the study cohort were uniformly treated with six cycles of R-CHOP, followed by two cycles of R (R-CHOP-14). The primary endpoint was event-free survival (EFS) at 2 years (for definition, see the “Statistical analysis” section), and the secondary endpoints were progression-free survival (PFS) and OS at 2 and 5 years as well as objective responses according to international criteria .
DNA extraction and quantification
Genomic DNA was extracted with the GeneRead DNA FFPE kit (Qiagen, Nussloch, Germany) following manufacturers’ recommendations with minor modifications. Briefly, one to three 10–25 μm thick tissue sections were deparaffinized by several xylene washes, rehydrated and digested with proteinase K overnight at 56 °C in a shaking heat block. Following digestion, the samples were incubated for 1 h at 90 °C to reverse fixation-induced DNA crosslinks and inactivate proteinase K. Thereafter, uracil-N-glycosylase (UNG) enzyme was added to remove artificially (formalin) induced uracils and reduce the number of false-positive C > T transitions. After incubation at 37 °C for 1 h, the samples were loaded into a DNA purification column and were washed and eluted in 40 μl of nuclease-free water. DNA yields were quantified with the Qubit High sensitivity DNA assay (Life Technologies, Eugene, OR, USA).
Targeted HTS sequencing variant calling and filtering
A target enrichment panel was designed to cover mutational hotspots or all exons of genes most frequently mutated in B cell lymphoid neoplasms according to the COSMIC database (release v70) and manual review of the literature . Sequencing libraries were constructed exactly as described previously . One microliter of the prepared library was used for the Bioanalyzer High Sensitivity DNA assay (Agilent, USA) to confirm expected DNA fragment length distribution. Quantification was performed with the Ion Library Quantitation kit (Thermo Fisher Scientific, Carlsbad, CA, USA) following the original protocol. Libraries were diluted to 40 pM, loaded to the Ion540 sequencing chips by automated IonChef instrument, and sequenced with the Ion Torrent S5 XL machine (Thermo Fisher Scientific, Carlsbad, CA, USA). The depth of coverage, coverage uniformity, and number of variants called per sample are summarized in Additional file 1: Table S1. Mutation identification was performed by the Variant caller plug-in v5.0 of the Torrent Suite (Thermo Fisher Scientific) using default low stringency parameters for somatic mutation calling. Mutations were annotated using the Ion Reporter variant annotation workflow v5.0 and dbNSFP v3.0 database . MetaLR rank score was used to predict the functional impact of non-synonymous point mutations to the encoded protein . After annotation, the variants were subjected to additional, more stringent, and quality- and relevance-based filtering by criteria as detailed in Table 1. All variants with variant allelic frequency >5% were used in downstream analysis.
Finally, aligned BAM files were manually inspected at sites of all remaining variants to exclude false-positive mutations or other artifacts introduced during library preparation .
All statistical analyses were performed using the Statistical Package of Social Sciences (IBM SPSS version 22.0, Chicago, IL, USA) for Windows. EFS was calculated from registration to progressive disease or relapse, death of any cause, and initiation of any non-protocol anti-cancer treatment because of lymphoma symptoms or need of concomitant radiotherapy. PFS was calculated from registration to progressive disease or relapse, and death of any cause. OS was calculated from registration to death. Patients not experiencing an event were censored at last follow-up. The survival probabilities were determined using the Kaplan–Meier method, and groups were compared using the log-rank test. Factors of prognostic significance in univariable models underwent multivariable analysis using the Cox proportional hazards model. For other endpoints, differences between groups were tested either with t test, Wilcoxon rank-sum test, or Fisher’s exact test, as appropriate. In all tests, p values are two-sided, considered significant if <0.05, and not corrected for multiple testing. For survival analysis within the COO subgroups, p values were corrected for multiple testing and were considered significant if <0.017.
Patients and clinico–pathologic characteristics
In total, 84 patient samples were included. Four samples yielded insufficient DNA for HTS library preparation, and another four samples failed quality control after sequencing (low coverage uniformity, low percentage of reads mapped to target), leaving 76 cases for further study. Patient characteristics are presented in Table 2. Survival data was available for all 76 patients. The 2-year OS, PFS, and EFS of the entire cohort were 91% (95% CI 84 to 98%), 78% (95% CI 68 to 88%), and 63% (95% CI 53 to 75%), respectively. There was no significant difference in any of the survival endpoints between the two different immunohistochemically determined COO subtypes. IHC and FISH data on all samples were collected as part of the previous translational study focusing on the prognostic value of genotypic and phenotypic characteristics . “Double-hit score,” representing IHC-based evaluation of tumor cell expression of BCL2 and MYC (Table 2), showed that 10/76 (13%) cases were positive for both MYC and BCL2. None of the investigated cases had chromosomal translocations affecting both MYC and BCL2 genes. Thirty of 76 (39%) and 1/76(1%) of cases were positive for FOXP1 and CD5, respectively, both associated with worse DLBCL outcomes in our previous study . All investigated cases were Epstein–Barr virus (EBV)-negative as determined by in situ hybridization for EBV early RNA (EBER).
In total, 317 somatic alterations were detected. Missense mutations were the most frequent, 232/317 (73%), followed by 40 (13%) nonsense mutations and 37 (12%) frameshift insertions/deletions (Fig. 1a–b and Additional file 1: Table S2). Sixty-eight percent (46/68) of the investigated genes were mutated at least once, and at least one somatic mutation was identified in 71/76 (93%) of investigated tumor samples. The most frequently mutated gene was KMT2D (34% of cases), followed by cytokine signaling 1 (SOCS1) (28%), GNA13 (17%), and B2M (15%) (Fig. 1c). EZH2 (p = 0.0176) along with GNA13 (p = 0.0181) and SGK1 (p = 0.028) were more frequently mutated in GCB–DLBCL. On average, GCB–DLBCL cases had more mutated genes than ABC-DLBCL cases (mean 3.9 vs 2.8, p = 0.018). However, there was no significant difference between COO subtypes considering the cumulative number of mutations per case (some genes had multiple mutations in the same tumor).
We used the ensemble MetaLR score that integrates nine different functional prediction algorithms and population allele frequencies to predict the functional impact of the detected missense mutations. Ninety-eight percent of all missense mutations were successfully annotated, and 86/227 (38%) were reliably predicted to be deleterious. Functional effects did not distribute evenly among investigated genes. In some (GNA13, TP53, CREBBP, TNFAIP), all or nearly all substitutions were either nonsense or predicted damaging missense mutations, whereas in others (PIM1, BCL2), the majority of mutations were “neutral,” and thus likely played a passenger role in lymphomagenesis (Fig. 2).
We performed Kaplan–Meier survival analysis of all genes that were mutated in at least 4 cases within our cohort and, also, on combinations of mutated genes that act on the same pathway. Only significant findings are reported.
SOCS1 mutated cases had better PFS compared to wild-type cases (p = 0.022), but the difference was not significant for OS and EFS (p = 0.11 and p = 0.182, respectively) (Fig. 3a). In total, 47 mutations of SOCS1 were detected, affecting 21 cases (12 non-GCB and 9 GCB). Mutations were distributed evenly along the coding sequence with two hotspots at Ala17 and Phe79 (four and three mutations, respectively) (Additional file 2: Figure S1). In all but 2 cases, at least one truncating (frameshift insertion/deletion or nonsense) or predicted damaging missense mutation was present, suggesting a deleterious effect on protein synthesis and/or function. Nineteen of 47 (40%) SOCS1 mutations were potentially induced by aberrant somatic hypermutation (aSHM), as they occurred within the RGYW/WRCY DNA sequence motif known to be targeted by activation-induced deaminase.
Since CREBBP and EP300 have close functional interactions in acetylating histone and nonhistone proteins, we decided to evaluate their combined prognostic value. Cases with mutations in CREBBP or EP300 had significantly worse OS, EFS, and PFS (p = 0.0015, p = 0.005, and p = 0.017, respectively) (Fig. 3b–d). In total, 14 cases were affected, with one heterozygous mutation in each case. All detected variants were missense substitutions localized in the acetyltransferase HAT domain with predicted deleterious effects on protein function (Additional file 2: Figure S1). However, the existence of more mutations cannot be excluded, since the applied gene panel only partly investigated the coding sequence of both CREBBP and EP300. CREBBP mutations alone had a negative impact on OS, but this did not reach statistical significance (p = 0.061). Analysis of potential linkage of these mutations to aSHM showed no such evidence.
In our previous study, we reported that overexpression of FOXP1 evaluated by IHC was prognostic of worse OS . We combined mutation- and IHC-based prognostic markers to determine whether this improved identification of DLBCL cases with poor outcomes. Fifty-one percent (39/76) of cases overexpressed FOXP1 and/or had mutations of CREBBP/EP300. Five tumors shared both of these features. Survival analysis showed significantly worse OS, PFS, and EFS in cases with either FOXP1 overexpression or CREBBP or EP300 mutations (p = 0.0003, p = 0.006, and p = 0.0002, respectively) (Fig. 3e–f).
ATM mutations emerged as prognostic factor for worse EFS, PFS, and OS within the GCB–DLBCL subgroup (p = 0.000002, p = 0.00041, and p = 0.0065, respectively; Fig. 3g–i). In total, 11 mutations affecting 10 DLBCL cases were found (5 in GCB and 5 in non-GCB instances). Only in one case ATM mutations were reliably predicted to be deleterious, while the functional significance of other mutations remained unclear. Also, there was no evidence of targeting of this gene by aberrant somatic hypermutation as only 2/11 mutations occurred within the RGYW motif.
In multivariable analysis that included IPI and failure to achieve complete remission (FACR), CREBBP and EP300 mutational status (but not SOCS1 status) emerged as the most significant factor predicting worse OS (p = 0.013, HR 4.12, CI 1.38–14.15) (Table 3). While these mutations were also significant for predicting worse EFS and PFS (p = 0.021 and p = 0.041, respectively), FACR was a superior predictor of the latter two endpoints (p = 0.0006 and p = 0.001).
There is a critical need for additional tumor-based prognostic markers in DLBCL to identify high-risk DLBCL patients prior to first treatment that might benefit from alternative risk-adjusted therapies. Tumor mutations represent promising candidates for outcome prognostication due to the unequivocal nature of results obtained by respective mutational detection techniques and the increased availability and robustness of HTS. We have demonstrated that targeted mutational analysis of a relatively small but uniformly treated and prospectively followed up patient cohort can reveal significant prognostic associations in DLBCL. This was only possible (a lesson learned) because of the proper trial design contemplating central collection of tissue for translational analysis as an integral study part. Additional lessons learned from this prospective trial such as the importance of a central diagnostic pathology review as well as handling of biological entities and subentities in the spectrum of so-called high grade B cell lymphomas have been discussed in a paper of ours recently published in this journal .
We found that deleterious mutations in two acetyltransferase genes, CREBBP and EP300, which belong to the KAT3 family of histone/protein lysine acetyltransferases, predict worse OS, PFS, and EFS in DLBCL independent of IPI and FACR. Point mutations or deletions of CREBBP/EP300 reportedly affect 39% of all DLBCL cases . In line with other studies, we confirm that CREBBP mutations are more frequent than EP300 and tend to occur more often in GCB–DLBCL . In the context of cancerogenesis, CREBBP and EP300 act as tumor suppressors. It has been shown that loss of one CREBBP allele leads to reduced acetylation and inactivation of p53, impaired expression of glucocorticoid-receptor-responsive genes, and upregulation of BCL6 [31, 32]. Finally, recent data suggest that heterozygous deleterious CREBBP mutations lead to decreased global histone H3 lysine 14 (H3K14), K18, and K27 acetylation and reduced MHC class II expression (Hashwah et al., currently under peer review). Association with poorer outcomes were suggested in studies that found CREBBP mutations in 20% of relapsed/refractory GCB–DLBCL , and in a large proportion of relapsed acute lymphoblastic leukemia patients [31, 34]. Despite this important clue, however, the prognostic value of CREBBP/EP300 mutations in DLBCL has not been previously reported.
SOCS1 mutations predicted excellent PFS in our cohort, but differences in OS and EFS were not significant. Suppressor of cytokine signaling 1 (SOCS1 ) is a known inhibitor of JAK/STAT-dependent signal transduction, which binds to phosphorylated JAK and marks it for proteosomal degradation [35, 36]. Mutations of SOCS1 have been previously shown to be associated with favorable survival in DLBCL  and are also frequently detected in Hodgkin lymphoma and primary mediastinal B cell lymphoma (PMBCL), both with a relatively favorable prognosis . We also previously showed that SOCS1 mutations occurred exceptionally in non-relapsing primary DLBCL whereas being completely absent in relapsing DLBCL cases, supporting the association with favorable prognosis . Schif et al. reported that DLBCL bearing truncating SOCS1 mutations have excellent OS, whereas those with only missense mutations have markedly worse prognosis . Our data, however, does not confirm such distinction: 6 of 19 SOCS1-mutated cases had only missense mutations, but their PFS was equally favorable as those with truncating mutations. The SOCS1 mutational frequency in our cohort (28%) was higher than reported average (~13%) . All detected mutations had relatively high variant allelic frequencies (median 28%, range 6–78%); therefore, our higher mutation rate cannot be explained by comparably high sequencing depth and sensitive detection of subclonal mutations that could have been missed by exome-scale sequencing studies. A potential explanation might be bias of our cohort toward better than average survival, e.g., due to study protocol exclusion of patients with performance status >2 on the ECOG scale, with symptomatic central nervous system disease, or with HIV, and/or hepatitis infection, thus being potentially enriched for cases with better prognosis that consequently more often bear SOCS1 mutations.
The prognostic value of multiple other gene mutations such as TP53, MYD88, FOXP1, and FOXP2 in DLBCL has been suggested previously [22, 23, 40, 41]. In our cohort, TP53-mutated cases had worse OS, but this difference did not reach statistical significance. Also, the low number (n = 5) of MYD88 L265P-mutated cases did not allow detection of any reliable prognostic associations. Analogously, despite ATM mutations were consistently associated with worse outcomes in GCB–DLBCL in our study, this observation was based on a small number of events and cases bearing mutations and therefore remains to be validated on other larger collectives. While we were unable to identify mutations of FOXP1, we previously reported that overexpression of the FOXP1 protein is associated with worse OS in the investigated DLBCL cohort . Here, we show that combination of HTS- and IHC-based prognostic markers (CREBBP/EP300 mutations and FOXP1 overexpression) enables even better stratification between patients with good and worse prognosis.
It is unlikely that any single biomarker will significantly improve risk stratification in DLBCL due to the profound heterogeneity of this disease. This point is illustrated by 3 cases in our cohort with mutations of both SOCS1 and CREBBP (Fig. 2, cases UPN57, UPN112, and UPN115) that clinically behaved like CREBBP-mutated cases and had worse prognosis. It is more likely that a combination of several different types of markers established on uniformly treated prospective cohorts and selected by robust statistical methods would provide models that could be prospectively validated on larger DLBCL collections. A good example of such an updated composite prognostic model is the m7-FLIPI, a recently developed score for follicular lymphoma, which incorporates mutations and clinical factors and provides superior prognostication compared to traditionally used clinical factor-based FLIPI .
Deleterious mutations in the HAT domain of the acetyltransferases CREBBP and EP300 are associated with worse OS, PFS, and EFS in DLBCL. ATM mutations are prognostic for worse survival at all clinical endpoints in GCB–DLBCL, but due to low case number remains to be verified on larger collectives. The previously reported beneficial prognostic role of SOCS1 mutations in DLBCL is valid for predicting better PFS.
Autologous stem cell transplantation
Aberrant somatic hypermutation
Binary sequence alignment/map format
Diffuse large B cell lymphoma
Epstein-Barr virus early RNA
Eastern Cooperative Oncology Group
Failure to achieve complete remission
Fluorescence in situ hybridization
Follicular Lymphoma International Prognostic Index
Germinal center B cell
International Prognostic Index
Positron emission tomography/computed tomography
Primary mediastinal B cell lymphoma
Rituximab, cyclophosphamide, hydroxydaunorubicin, vincristine, and prednisone
Coiffier B, Thieblemont C, Van Den Neste E, Lepeu G, Plantier I, Castaigne S, et al. Long-term outcome of patients in the LNH-98.5 trial, the first randomized study comparing rituximab-CHOP to standard CHOP chemotherapy in DLBCL patients: a study by the Groupe d’Etudes des Lymphomes de l’Adulte. Blood. 2010;116:2040–5.
Friedberg JW. Relapsed/refractory diffuse large B-cell lymphoma. Hematology. 2011;2011:498–505.
Costa LJ, Maddocks K, Epperla N, Reddy NM, Karmali R, Umyarova E, et al. Diffuse large B-cell lymphoma with primary treatment failure: ultra-high risk features and benchmarking for experimental therapies. Am J Hematol. 2017;92:161–70.
Hitz F, Connors JM, Gascoyne RD, Hoskins P, Moccia A, Savage KJ, et al. Outcome of patients with primary refractory diffuse large B cell lymphoma after R-CHOP treatment. Ann Hematol. 2015;94:1839–43.
Gisselbrecht C, Glass B, Mounier N, Singh Gill D, Linch DC, Trneny M, et al. Salvage regimens with autologous transplantation for relapsed large B-cell lymphoma in the rituximab era. J Clin Oncol. 2010;28:4184–90.
A predictive model for aggressive non-Hodgkin’s lymphoma. The International Non-Hodgkin’s Lymphoma Prognostic Factors Project. N Engl J Med. 1993;329:987–94.
Sehn LH, Berry B, Chhanabhai M, Fitzgerald C, Gill K, Hoskins P, et al. The revised International Prognostic Index (R-IPI) is a better predictor of outcome than the standard IPI for patients with diffuse large B-cell lymphoma treated with R-CHOP. Blood. 2007;109:1857–61.
Zhou Z, Sehn LH, Rademaker AW, Gordon LI, LaCasce AS, Crosby-Thompson A, et al. An enhanced International Prognostic Index (NCCN-IPI) for patients with diffuse large B-cell lymphoma treated in the rituximab era. Blood. 2014;123:837–42.
Tzankov A, Xu-Monette ZY, Gerhard M, Visco C, Dirnhofer S, Gisin N, et al. Rearrangements of MYC gene facilitate risk stratification in diffuse large B-cell lymphoma patients treated with rituximab-CHOP. Mod Pathol. 2014;27:958–71.
Lenz G, Wright G, Dave SS, Xiao W, Powell J, Zhao H, et al. Stromal gene signatures in large-B-cell lymphomas. N Engl J Med. 2008;359:2313–23.
Meyer PN, Fu K, Greiner TC, Smith LM, Delabie J, Gascoyne RD, et al. Immunohistochemical methods for predicting cell of origin and survival in patients with diffuse large B-cell lymphoma treated with rituximab. J Clin Oncol. 2011;29:200–7.
Scott DW, Wright GW, Williams PM, Lih C-J, Walsh W, Jaffe ES, et al. Determining cell-of-origin subtypes of diffuse large B-cell lymphoma using gene expression in formalin-fixed paraffin embedded tissue. Blood. 2014;123:1214–7.
Tzankov A, Leu N, Muenst S, Juskevicius D, Klingbiel D, Mamot C, et al. Multiparameter analysis of homogeneously R-CHOP-treated diffuse large B cell lymphomas identifies CD5 and FOXP1 as relevant prognostic biomarkers: report of the prospective SAKK 38/07 study. J Hematol Oncol. 2015;8:70.
Miyazaki K, Yamaguchi M, Suzuki R, Kobayashi Y, Maeshima AM, Niitsu N, et al. CD5-positive diffuse large B-cell lymphoma: a retrospective study in 337 patients treated by chemotherapy with or without rituximab. Ann Oncol Off J Eur Soc Med Oncol. 2011;22:1601–7.
Salles G, de Jong D, Xie W, Rosenwald A, Chhanabhai M, Gaulard P, et al. Prognostic significance of immunohistochemical biomarkers in diffuse large B-cell lymphoma: a study from the Lunenburg Lymphoma Biomarker Consortium. Blood. 2011;117:7070–8.
Winter JN, Li S, Aurora V, Variakojis D, Nelson B, Krajewska M, et al. Expression of p21 protein predicts clinical outcome in DLBCL patients older than 60 years treated with R-CHOP but not CHOP: a prospective ECOG and Southwest Oncology Group correlative study on E4494. Clin Cancer Res. 2010;16:2435–42.
Song G, Ni H, Zou L, Wang S, Tian F, Liu H, et al. Expression of CD40 is a positive prognostic factor of diffuse large B-cell lymphoma treated with R-CHOP (rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone). Onco Targets Ther. 2016;9:3799–805.
Pasqualucci L, Trifonov V, Fabbri G, Ma J, Rossi D, Chiarenza A, et al. Analysis of the coding genome of diffuse large B-cell lymphoma. Nat Genet. 2011;43:830–7.
Morin RD, Mungall K, Pleasance E, Mungall AJ, Goya R, Huff RD, et al. Mutational and structural analysis of diffuse large B-cell lymphoma using whole-genome sequencing. Blood. 2013;122:1256–65.
Zhang J, Grubor V, Love CL, Banerjee A, Richards KL, Mieczkowski PA, et al. Genetic heterogeneity of diffuse large B-cell lymphoma. Proc Natl Acad Sci U S A. 2013;110:1398–403.
Pasqualucci L, Dalla-Favera R. The genetic landscape of diffuse large B-cell lymphoma. Semin Hematol. 2015;52:67–76.
Xu-Monette ZY, Wu L, Visco C, Tai YC, Tzankov A, Liu W, et al. Mutational profile and prognostic significance of TP53 in diffuse large B-cell lymphoma patients treated with R-CHOP: report from an International DLBCL Rituximab-CHOP Consortium Program Study. Blood. 2012;120:3986–96.
Fernández-Rodríguez C, Bellosillo B, García-García M, Sánchez-González B, Gimeno E, Vela MC, et al. MYD88 (L265P) mutation is an independent prognostic factor for outcome in patients with diffuse large B-cell lymphoma. Leukemia. 2014;28:2104–6.
Trinh DL, Scott DW, Morin RD, Mendez-Lago M, An J, Jones SJM, et al. Analysis of FOXO1 mutations in diffuse large B-cell lymphoma. Blood. 2013;121:3666–74.
Mamot C, Klingbiel D, Hitz F, Renner C, Pabst T, Driessen C, et al. Final results of a prospective evaluation of the predictive value of interim positron emission tomography in patients with diffuse large B-cell lymphoma treated with R-CHOP-14 (SAKK 38/07). J Clin Oncol. 2015;33:2523–9.
Cheson BD, Pfistner B, Juweid ME, Gascoyne RD, Specht L, Horning SJ, et al. Revised response criteria for malignant lymphoma. J Clin Oncol. 2007;25:579–86.
Juskevicius D, Lorber T, Gsponer J, Perrina V, Ruiz C, Stenner-Liewen F, et al. Distinct genetic evolution patterns of relapsing diffuse large B-cell lymphoma revealed by genome-wide copy number aberration and targeted sequencing analysis. Leukemia. 2016;30:2385–95.
Liu X, Wu C, Li C, Boerwinkle E. dbNSFP v3.0: a one-stop database of functional predictions and annotations for human nonsynonymous and splice-site SNVs. Mutat. 2016;37:235–41.
Dong C, Wei P, Jian X, Gibbs R, Boerwinkle E, Wang K, et al. Comparison and integration of deleteriousness prediction methods for nonsynonymous SNVs in whole exome sequencing studies. Hum Mol Genet. 2015;24:2125–37.
McCall CM, Mosier S, Thiess M, Debeljak M, Pallavajjala A, Beierl K, et al. False positives in multiplex PCR-based next-generation sequencing have unique signatures. J Mol Diagn. 2014;16:541–9.
Pasqualucci L, Dominguez-Sola D, Chiarenza A, Fabbri G, Grunn A, Trifonov V, et al. Inactivating mutations of acetyltransferase genes in B-cell lymphoma. Nature. 2011;471:189–95.
Mullighan CG, Zhang J, Kasper LH, Lerach S, Payne-Turner D, Phillips LA, et al. CREBBP mutations in relapsed acute lymphoblastic leukaemia. Nature. 2011;471:235–9.
Morin RD, Assouline S, Alcaide M, Mohajeri A, Johnston RL, Chong L, et al. Genetic landscapes of relapsed and refractory diffuse large B-cell lymphomas. Clin Cancer Res. 2016;22:2290–300.
Inthal A, Zeitlhofer P, Zeginigg M, Morak M, Grausenburger R, Fronkova E, et al. CREBBP HAT domain mutations prevail in relapse cases of high hyperdiploid childhood acute lymphoblastic leukemia. Leukemia. 2012;26:1797–803.
Kamizono S, Hanada T, Yasukawa H, Minoguchi S, Kato R, Minoguchi M, et al. The SOCS box of SOCS-1 accelerates ubiquitin-dependent proteolysis of TEL-JAK2. J Biol Chem. 2001;276:12530–8.
Endo TA, Masuhara M, Yokouchi M, Suzuki R, Sakamoto H, Mitsui K, et al. A new protein containing an SH2 domain that inhibits JAK kinases. Nature. 1997;387:921–4.
Schif B, Lennerz JK, Kohler CW, Bentink S, Kreuz M, Melzner I, et al. SOCS1 mutation subtypes predict divergent outcomes in diffuse large B-cell lymphoma (DLBCL) patients. Oncotarget. 2013;4:35–47.
Steidl C, Gascoyne RD. The molecular pathogenesis of primary mediastinal large B-cell lymphoma. Blood. 2011;118:2659–69.
Forbes SA, Beare D, Gunasekaran P, Leung K, Bindal N, Boutselakis H, et al. COSMIC: exploring the world’s knowledge of somatic mutations in human cancer. Nucleic Acids Res. 2015;43:D805–11.
Gascoyne DM, Banham AH. The significance of FOXP1 in diffuse large B-cell lymphoma. Leuk Lymphoma. 2016;58:1037–51.
Wong KK, Gascoyne DM, Soilleux EJ, Lyne L, Spearman H, Roncador G, et al. FOXP2-positive diffuse large B-cell lymphomas exhibit a poor response to R-CHOP therapy and distinct biological signatures. Oncotarget. 2016;7:52940–56.
Pastore A, Jurinovic V, Kridel R, Hoster E, Staiger AM, Szczepanowski M, et al. Integration of gene mutations in risk prognostication for patients receiving first-line immunochemotherapy for follicular lymphoma: a retrospective analysis of a prospective clinical trial and validation in a population-based registry. Lancet Oncol. 2015;16:1111–22.
We are grateful to Prof. Anne Müller for her insights on the functional relevance of CREBBP/EP300 mutations and fruitful scientific discussions. Also, special thanks to Valeria Perrina and Sibylle Tschumi for their help in performing high-throughput sequencing.
The study was supported by the Oncosuisse grant Nr. OCS 02072-04-2007 and Krebsliga beider Basel.
Availability of data and materials
Targeted high-throughput sequencing data was deposited in the Sequence Read Archive (SRA) Ref Nr.: SRP098666.
DJus performed the sequencing experiments and data analysis and wrote the manuscript. DJuc performed the sequencing experiments and analysis. DK and KM performed the statistical analysis and revised the manuscript. SD performed the immunohistochemistry analysis and revised the manuscript. AT designed the study; performed the immunohistochemistry, FISH, and statistical analysis; patrially wrote and revised the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Approval for the clinical trial was obtained from the ethics committee of each participating institution.
Consent for translational study was obtained from all patients, and it was approved by the Ethics Committee of Northwestern and Central Switzerland (EKNZ 2014-252).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Statistics of high-throughput sequencing and variant calling. Table S2. Somatic mutations detected by high-throughput sequencing in 76 cases of primary DLBCL. (XLSX 384 kb)
Mapping of the detected mutations to the protein sequences within respective domains in the ten most frequently mutated and relevant genes is this study. Mutations in CREBBP and EP300 always occurred within the HAT domain responsible for acetylation. In SOCS1, mutations were distributed evenly along the protein sequence. Y-axis represents the number of mutations detected at a particular amino acid position. Green–missense mutation; black–frameshift insertion/deletion; purple–nonsense mutation. (PDF 442 kb)