Abstract
Background The aim of this study was to determine the impact of a gain of the MALT1 gene on gene expression and clinical parameters in diffuse large B-cell lymphoma.Design and Methods We analyzed 116 cases of diffuse large B-cell lymphoma by fluorescence in situ hybridization, array-based comparative genomic hybridization, and transcriptional profiling.Results A gain of 18q21 including MALT1 was detected in 44 cases (38%) and was accompanied by a gain of BCL2 in 43 cases. All cases with a 18q21/MALT1 gain showed BCL2 protein whereas 79% in the group without a 18q21/MALT1 gain did so (p<0.001). Cases with 18q21/MALT1 gain more frequently showed an activated B-cell-like (ABC) gene expression signature (65%) than a germinal center B-cell-like (GCB) one (23%) (p<0.001). Ninety-eight genes including MALT1, BCL2, and some selected nuclear factor-κB target genes were differentially expressed between the two genetic groups of diffuse large B-cell lymphoma. By global testing of each chromosome, we identified 33 genes, all located on chromosome 18q, which were differentially expressed between the two genetic groups independently of the ABC/GCB status. In multivariate analysis, the 18q21/MALT1 status represented an independent negative prognostic factor for overall survival (p=0.03).Conclusions In diffuse large B-cell lymphoma, gain of 18q21 including MALT1 is significantly associated with differential expression of genes located on 18q, the ABC gene expression subtype, increased BCL2 gene and protein expression and might indicate an unfavorable prognosis.Introduction
Diffuse large B-cell lymphoma (DLBCL) is the most frequent lymphoma in adults, accounting for 30% to 40% of lymphoid neoplasms.1 The diversity in clinical presentation and outcome, as well as the pathological and biological heterogeneity suggest that DLBCL comprises several disease entities that may require different therapeutic approaches.2–4
Recent gene expression profiling studies have identified several prognostically different subgroups of DLBCL with gene expression patterns indicative of the different stages of B-cell differentiation, termed germinal center B-cell-like (GCB), activated B-cell-like (ABC), and primary mediastinal (PM) DLBCL,5–7 as well as activation of distinct biological pathways, especially oxidative phosphorylation, B-cell receptor/proliferation, and host response.8,9 In addition, we recently defined the gene expression signature of Burkitt’s lymphoma and showed that a small proportion of DLBCL also bears this signature.10 Patients with a “molecular Burkitt’s lymphoma signature” were younger, in a less advanced clinical stage, and a had better clinical outcome compared to patients without this signature.10 The 5-year survival rates of patients with GCB-DLBCL, ABC-DLBCL, and PM-DLBCL are 59%, 30%, and 64%, respectively.5–7 GCB-DLBCL is characterized by recurrent REL amplification, BCL2 translocations,6,11 and ongoing somatic hypermutation of the immunoglobulin genes.12 In contrast, ABC-DLBCL has constitutively activated transcription factor nuclear factor (NF)-κB, which plays a critical role in the pathogenesis of this lymphoma,13–15 and BCL2 protein overexpression in the absence of the t(14;18).16 The mechanisms underlying BCL2 overexpression in t(14;18)-negative DLBCL are largely unknown. Since BCL2 is a target gene for NF-κB17 and NF-κB is constitutively expressed in ABC-DLBCL, BCL2 upregulation may be mediated through NF-κB. Another possible mechanism resulting in BCL2 protein overexpression is a gain of BCL2, as recently shown by Iqbal et al.16
In this scenario, the role of the MALT1 gene, which lies close to BCL2 on chromosomal band 18q21 and has an important impact on NF-κB activation, was not known and was, therefore, explored in the present study.
Design and Methods
Patients
This study is part of the German network project Molecular Mechanisms of Malignant Lymphomas, for which ethic approval was obtained. Mature aggressive B-cell lymphomas with more than 70% of tumor cells were analyzed by a comprehensive approach including gene expression profiling, array-based comparative genomic hybridization (array-CGH), and interphase fluorescence in situ hybridization (FISH).10 The diagnoses were established by seven expert hematopathologists, who evaluated morphology and immunohistochemistry in a centralized panel review as previously described.10 The intensity of BCL2 expression by the tumor cells was determined in relation to that of BCL2-positive reactive T-cells. Only lymphomas with an intensity of BCL2 expression similar to or stronger than that of reactive T cells and with more than 25% of BCL2-positive tumor cells were scored as BCL2-positive.
We included all cases analyzed in the above mentioned network project until December 2005 which fulfilled the following criteria: (i) a histopathologic diagnosis of DLBCL according to the WHO classification, (ii) lack of a molecular Burkitt’s lymphoma signature, (iii) samples not from patients with relapse, human immunodeficiency virus-infection, or PM-DLBCL, (iv) interphase FISH evaluable for imbalances of the BCL2 and MALT1 loci.
From 142 DLBCL cases with available interphase FISH results, 10 cases were excluded because they carried a molecular Burkitt’s lymphoma signature and 16 further cases were excluded because relapse of the lymphoma, a human immunodeficiency virus-infection, or a PM-DLBCL was diagnosed. None of the ten DLBCL with a molecular Burkitt’s lymphoma signature showed a gain of MALT1 or BCL2. Thus, a total of 116 cases entered the study (Table 1). Clinical and survival data were available for 81 of these patients (Table 2). The median follow-up time of these patients was 5 years (range, 0 to 209 months).
Fluorescence in situ hybridization
Interphase FISH was performed on 5 μm thick frozen tissue sections according to standard methods. Two hundred interphase nuclei were analyzed in each case by at least two independent investigators. Gene amplification was defined by the presence of a cloud-like signal accumulation or more than five individual signals in the interphase cells. A gain of the respective gene was defined by the presence of three to five hybridization signals per interphase nucleus occurring in more than 5% of nuclei. In cases with a gain of MALT1 or BCL2, as shown by interphase FISH and a balanced array-CGH status for these genes, further FISH analyses were evaluated for imbalances of the BCL6, MYC, and IGH loci to differentiate polysomies of 18q21 from poly-ploid karyotypes without a real 18q21 gain. The MALT1 and the BCL2 loci were analyzed using dual color FISH with probe LSI MALT1 Dual Color, Break Apart Rearrangement Probe and probe LSI IGH/BCL2 Dual Color, Dual Fusion Translocation Probe, respectively (both from Abbott Vysis, Downers Grove, IL, USA). In addition, as described by Hummel et al.,10 breakpoints within the MYC, BCL6, and IGH genes were analyzed with the LSI MYC, LSI BCL6, and LSI IGH break apart assays (all from Abbott Vysis) and MYC partners with the LSI IGH/MYC (Abbott Vysis) and home-brewed IGK-MYC and IGL-MYC double-color, double-fusion assays.10
Array-CGH
Array-CGH was performed as previously described.10,18,19 A 2.8k array CGH chip containing 2799 fragments was applied. The previously published GoldenPath clone set20 served as the backbone for the set up of the chip. In total, 1500 DNA clones at intervals of approximately 2 Mb were selected. This compilation was enriched by 600 DNA fragments, which represent critical regions in B-cell neoplasms as well as 699 additional DNA fragments containing proto-oncogenes or tumor suppressor genes. DNA preparations of all BAC and PAC clones and subsequent degenerate oligonucleotide primed (DOP) polymerase chain reaction (PCR)-mediated DNA amplification were performed as described by Fiegler et al.21 The PCR products were spotted onto Corning CMT-Gaps II glass slides using an Omnigrid microarrayer (Gene Machines, San Carlos, CA, USA). DNA labeling, hybridization and data acquisition was done as recently described.10,18,19 The ratios of two hybridizations with reversed dye labeling were averaged and normalized. The diagnostic cut-off level for each individual experiment was determined by calculating the mean and subsequently using plus/minus three standard deviations of all clones from chromosomes 1 to 22.
Affymetrix GeneChip hybridization and statistical analysis of the microarray data
Affymetrix hybridization was carried out with the U133A GeneChips according to the manufacturer's recommendations, as previously described.10 Probe intensities were normalized using a variance stabilization method.22 Gene expression levels were estimated by fitting an additive model23 employing a median polish routine. The significance of differential gene expression was analyzed with a two-step procedure, a global step and a gene-by-gene step. Global tests for differential gene expression were performed for predefined sets of genes using the ANCOVA test of Mansmann and Meister.24 This test allows several factors to be included in a linear model and analyzed separately for differential gene expression. The resulting empirical permutation p values were adjusted for multiplicity using Bonferroni’s correction since multiple chromosomes were tested at the same time. We included the ABC/GCB type of the lymphomas in the model when testing for the impact of a MALT1 gain on gene expression. In order to determine lists of differentially expressed genes we performed a gene-by-gene analysis using the method of Smyth et al.25 GCB and ABC sample labels were computed according to Wright et al.26
The data presented are available from the Gene Expression Omnibus (GEO) of the National Center for Biotechnology Information (www.ncbi.nlm.nih.gov/geo/) through GEO accession number GSE4475 (the patients’ identification numbers (MPI-number) are shown in Online Supplementary Figure S1).
Statistical analysis of genetic, histopathological, and clinical data
Mann-Whitney U, χ, Fisher’s exact, and log-rank tests were applied to test for differences between groups. Survival was calculated from the day of diagnosis until death or the end of follow-up. Time-to-treatment-failure (TTTF) was calculated from the day of diagnosis until progression, relapse, start of salvage therapy, end of therapy without complete remission or death. A Cox proportional hazard model was employed to analyze prognostic factors. Since not all parameters of the International Prognostic Index were available for all cases, age at diagnosis and Ann Arbor stage were used to adjust for known prognostic factors.27
Results
Interphase FISH and array CGH: the minimal overrepresented region on chromosome 18 lies in 18q21 and includes MALT1
According to array-CGH, the minimal overrepresented region on chromosome 18 was delineated to a 3.9 Mb region in 18q21.2 to 18q21.31 (clones RP11-7L24 to RP11-350K6) and encompassed 15 annotated genes including MALT1, but not BCL2. The telomeric margin of this region was mainly defined by one case, which showed an amplification of 18q21 including MALT1 and a normal BCL2 status. The centromeric margin was defined by another case showing an amplification of 18q21 spanning exactly the minimal overrepresented region between clones RP11-7L24 and RP11-350K6 including MALT1. Interestingly, the region telomeric to RP11-350K6 including BCL2 was gained but not amplified.
From the 40 cases with a 18q gain, 30 cases (75%) had a gain of 18q or part of it and ten (25%) had trisomy or polysomy of the whole chromosome.18 According to FISH, the MALT1 gene was gained in 44 of the 116 cases (38%). In all cases with a MALT1 gain but one, a concomitant gain of BCL2 was detected by interphase FISH and array-CGH. The latter case had MALT1 amplification in 64% of the cells and a normal BCL2 signal constellation (see array-CGH results). In the remaining cases, the signal numbers for MALT1 and BCL2 were similar, indicating a gain of the 18q21 region or part of it including both genes (further referred to as gain of 18q21/MALT1).
Among the 44 cases with a gain of MALT1, 28 (63%) showed mainly three copies, six (14%) showed four copies, six (14%) showed five copies, and four (9%) had an amplification of MALT1. The percentage of cells with a MALT1 or BCL2 gain ranged from 11% to 90%. Array-CGH results corresponded well with the inter-phase FISH results in 40 cases. In the remaining four cases, interphase FISH revealed a trisomy of MALT1 and BCL2 together with two signals for the IGH, BCL6, and MYC probes, but array-CGH showed a normal signal constellation for the BCL2 and MALT1 genes. In these four cases, the percentage of cells bearing three hybridization signals was 11%, 18%, 22%, and 40%. The relatively low proportion of cells carrying the abnormality was possibly the reason why these gains were not detected by array-CGH.18 These cases were assigned as gain of 18q21. All cases without gains of MALT1 and/or BCL2 by FISH had normal array-CGH results for these regions.
The group without a 18q21/MALT1 gain showed a higher frequency of IGH-BCL2 fusions (p=0.06) and breaks in the BCL6 locus (p=0.07), but the difference in frequency was not statistically significant (Table 1). Three cases with 18q21/MALT1 gains showed a IGH-BCL2 fusion. In one case, array-CGH revealed a tetrasomy 18. In the other two cases, the gains involved the 18q12.1-18q21.33 and 18q21.1-18q21.33 regions, and the proximal border of the gain was telomeric of BCL2 in both cases. Breaks targeting MALT1 were not detected in any case.
Gain of 18q21 including MALT1 correlates with an ABC-DLBCL expression signature and differential gene expression overlaps with that observed in ABC- vs. GCB comparisons
Out of the 44 cases in the group with a 18q21/MALT1 gain, 29 (66%) showed an ABC signature and 10 (23%) had a GCB profile, whereas in the other 72 cases only 12 cases (17%) showed an ABC signature and 39 cases (54%) had a GCB pattern (p<0.001) (Table 1).
In order to determine differential gene expression between cases with and without a gain of 18q21/MALT1, we preceded in several steps. First, we applied a global test for differential expression between 18q21/MALT1-positive and 18q21/MALT1-negative patients, ignoring the ABC and GCB signatures. Global ANCOVA returned an empirical permutation p value of 0.004. Hence, there is differential gene expression between the two genetic subtypes of DLBCL. From a gene-by-gene test we obtained a list of 135 Affymetrix probe sets representing 98 differentially expressed genes with an estimated false discovery rate of 0.05. This list included MALT1 and BCL2 and some selected NF-κB target genes. In the next step, we included the GCB or ABC signature as a confounding factor in a full linear model using global ANCOVA. We restricted this part of the analysis to patients with an unambiguous GCB or ABC signature and excluded cases with an unclassified (type 3) signature. The global test for the influence of a 18q21/MALT1 gain on gene expression gave a negative result (p=0.41). It is not, therefore, possible to distinguish 18q21/MALT1- specific from ABC-specific effects on gene expression.
Genomic gain of 18q21 including MALT1 is significantly associated with RNA upregulation of genes located on chromosome 18q including MALT1
Up to this point, we had only considered the global gene expression profiles looking at all probe sets at the same time. This represents a massive multiple testing problem (>22,000 array features) compromising the statistical power of the screening. In the following step we, therefore, increased the power of the test by including prior knowledge of the chromosomal location of the genes. We applied a global test of gene expression for each chromosome (excluding the X and Y chromosomes, 22 tests) separately and adjusted the 22 p values for multiple testing. This test identifies chromosomes with global changes in gene expression between the two groups. Interestingly, only chromosome 18 displayed global changes in gene expression, which were independent of the ABC/GCB status (padjust-ed<0.03). To further elucidate this finding, we analyzed the genes on chromosome 18 independently using a gene-by-gene test and identified 47 probe sets representing 33 differentially expressed genes with an estimated false discovery rate of 0.05. (Online Supplementary Table S1 and Figure S1). These differentially expressed genes included MALT1, SERPINB8, NEDD4L, POLI, VPS4B, TXNL1, SMAD2, TCF4, SMAD4 and PIK3C3, but not BCL2.
Upregulation of NF-κB target genes in the group with a 18q21/MALT1 gain overlaps with ABC gene expression characteristics
Since MALT1 is known to activate the NF-κB pathway and BCL2 is transcriptionally regulated by NF-κB,17 we assessed the differential expression of NF-κB target genes in the groups with and without a 18q21/MALT1 gain. For this purpose, we selected all NF-κB target genes, which were previously reported by Feuerhake et al.14 ( Online Supplementary Table S2) and were present on the U133A-array used in this study (106 of the 116 probe sets reported, representing 68 genes). We subjected these 106 probe sets to a global test for differential gene expression24 between cases with and without a 18q21 gain. We found that this set of genes together displayed global changes in gene expression (p<0.0001), if the ABC/GCB status, which overlaps with the presence of a 18q21 gain, was ignored. However, the significance was lost when the ABC/GCB status was taken into account as a confounding factor within the test (p=0.12). We cannot, therefore, distinguish between 18q21/MALT1- related and ABC-related effects on the expression of the analyzed NF-κB target genes.
Correlation with histopathology
In the group with a 18q21/MALT1 gain, the histopathological panel review classified the DLBCL as follows: centroblastic variant (n=31; 71%), plasmablastic variant (n=4; 9%), immunoblastic variant (n=2; 4%), and anaplastic variant (n=1; 2%). In six cases (14%) a definitive assignment to one of the WHO DLBCL variants was not possible. In the cases without a 18q21/MALT1 gain the following DLBCL variants were observed: centroblastic variant (n=46; 64%), immunoblastic variant (n=7; 10%), plasmablastic (n=3; 4%), anaplastic variant (n=2; 3%). In 14 DLBCL cases without a 18q21/MALT1 gain (19%) a clear subgrouping according to the WHO classification was not possible.
All cases with a 18q21/MALT1 gain showed BCL2 protein expression, whereas 79% of the other cases did so (p<0.001) (Table 1). In addition, the Ki-67 index, a marker of cell proliferation, was significantly higher in the group with a 18q21/MALT1 gain (p=0.016) (Table 1). CD10 expression tended to be higher in the group without a 18q21/MALT1 gain (p=0.056), which might be related to the higher frequency of GCB-type DLBCL. No significant difference was found in BCL6 protein expression (p=0.1) (Table 1).
Interestingly, among the cases with a 18q21/MALT1 gain, four showed concomitant amplification of MALT1 and BCL2. These cases were diagnosed as DLBCL centroblastic variant and did not differ morphologically or immunophenotypically from the other cases with a gain of 18q21/MALT1. One further case showing a MALT1 amplification but normal BCL2 status was detected and diagnosed as DLBCL centroblastic variant without evidence that this case represented a transformed marginal zone B-cell lymphoma.
Correlation with clinical parameters and prognosis: gain of 18q21 including MALT1 is an independent negative prognostic factor in DLBCL
Concerning known clinical risk factors, the patients in the group with a 18q21/MALT1 gain showed a trend to more frequent bone marrow involvement (21% vs. 6%, p=0.097). There was no significant difference with respect to age, sex ratio, B-symptoms, bulky disease, extranodal involvement, lactate dehydrogenase level, or applied therapy including the use of rituxan (Table 2). In univariate analysis, we found a lower fraction of patients reaching a complete remission and a higher number of patients with progressive disease (p=0.014) in the group with a 18q21/MALT1 gain (Table 2). Time-to-treatment-failure was not significantly shorter in the group with a 18q21/MALT1 gain (p=0.12), but we found a trend to a shorter overall survival (5-year survival, 28% vs. 38%, p=0.09),
In the multivariate analysis incorporating age, Ann Arbor stage, a break within the MYC gene, ABC/GCB gene expression subtype, and 18q21/MALT1 status in all patients for whom information on these factors was available (n=71), the 18q21/MALT1 status remained an independent negative prognostic factor for overall survival [hazard ratio 2.13 (95% CI: 1.08–4.23), p=0.03].
Discussion
The 116 DLBCL analyzed in our study were diagnosed according to the histopathological and immunohistochemical characteristics described in the WHO classification and additionally according to novel gene expression features recently defined by our group.10
Forty-three of our 116 cases of DLBCL had a genomic gain of BCL2 and MALT1 and one case revealed an amplification of MALT1 together with a normal BCL2 gene status (together 38%), while 72 cases (62%) had normal copy numbers of the BCL2 and MALT1 genes not considering changes in ploidy. Interestingly, gain of the BCL2 gene was invariably associated with gain of the MALT1 gene in our series. Moreover, the presence of a 18q21/MALT1 gain correlated with upregulation of genes located on 18q including MALT1 at the RNA level, independently of the ABC or GCB gene expression signature, and was strongly associated with over-expression of BCL2 protein. We found a significant accumulation of cases with a 18q21/MALT1 gain among those with an ABC signature and 98 differentially expressed genes between the two genetic groups of DLBCL, including some selected NF-κB target genes. However, when we included the ABC/GCB signature as a confounding factor into the test, the influence of the 18q21/MALT1 gain on gene expression was no longer significant. The overlap of ABC and 18q21 gains complicated efforts to distinguish effects on gene expression specific to ABC status or 18q21 gains. Linear models generated in all steps of analysis allowed us to test for differential gene expression (18q21 gain vs. normal) while taking into account the presence of an additional confounding factor (ABC vs. GCB signature). We interpret this by the unequal distribution of the DLBCL cases with and without 18q21/MALT1 gains in the ABC- and GCB-categories. As a result, we could not distinguish 18q21/MALT1- specific from ABC-specific effects on gene expression.
BCL2 and MALT1 are both located on chromosome 18q21, about 5 Mb apart from each other, and are both key genes in the pathogenesis of different subtypes of non-Hodgkin’s lymphomas. MALT1 codes for a paracaspase, a scaffold protein recently proposed to be related to selective control of the c-Rel subunit in the canonical NF-κB pathway.28 The API2-MALT1 fusion protein resulting from the t(11;18)(q21;q21)29–31 and the overexpressed MALT1 protein in another MALT1-associated translocation in MALT lymphomas, the t(14;18) (q32;q21),32,33 have been shown to activate NF-κB and thereby promote cellular proliferation and resistance to apoptosis.34–36 In addition to the two translocations, amplifications and gains of MALT1 have been described in different subtypes of non-Hodgkin’s lymphomas,32,37 but their pathophysiological significance is not known.
BCL2 is a gene involved in the regulation of apoptosis and plays a major role in the response of malignant cells to a variety of stresses that may lead to apoptosis, including chemotherapy.38 One major mechanism leading to BCL2 protein overexpression is the t(14;18) (q32;q21) resulting in juxtaposition of the BCL2 gene and the IGH locus. The t(14;18) is found in most follicular lymphomas and about 20% of DLBCL.11,39 The t(14;18) occurs mainly in GCB-DLBCL, being present in approximately 35% of these cases. In GCB-DLBCL, the t(14;18) correlates with BCL2 mRNA and protein expression,39 indicating that this is the major mechanism in upregulation of BCL2 expression in this subgroup.
However, ABC-DLBCL also have high levels of BCL2 expression in the absence of t(14;18), indicating alternative mechanisms of BCL2 upregulation. Using inter-phase FISH, Iqbal et al. recently showed that BCL2 protein expression in ABC-DLBCL is significantly associated with genomic gain of BCL2.16 In another study, gain of 18q21-22 was detected by CGH in one third of ABC-DLBCL, but less frequently in GCB-DLBCL (10%) and PM-DLBCL (16%).40 Our study revealed compatible results with an accumulation of cases with a gain or amplification of the BCL2 and MALT1 genes in the ABC subtype as well as a highly significant association of these gains with BCL2 RNA and protein expression.
We additionally showed that a gain of BCL2 was accompanied by a gain of MALT1 in all of our cases of DLBCL. In line with this observation, gain of BCL2 was associated with gain of MALT1 in one study analyzing lymphoma cell lines by array-CGH and in a further study using array-CGH and quantitative PCR to evaluate DLBCL patients with chromosomal gains in 18q21-q22.32,40 Most other studies applied conventional CGH analysis or interphase FISH with BCL2-specific probes. It is not, therefore, possible to assess in these studies whether or not MALT1 was gained in addition to BCL2.16,41 Despite the scarcity of data, the frequent association of genomic gains of BCL2 and MALT1 might suggest a cooperative action of these two and possibly additional 18q21 genes in DLBCL.
Given that BCL2 is transcriptionally regulated by NF-κB17 and NF-κB is activated by MALT1, which was found to be upregulated on genomic and RNA levels in our study, we evaluated the differential expression of 68 selected NF-κB target genes in our DLBCL cases. There were significant changes in the expression of this set of genes between the two MALT1 genetic groups; however, the significance was lost when ABC/GCB status was included as a confounding factor in the test. Hence, also for the analyzed NF-κB target genes, it was not possible to distinguish between 18q21/MALT1- and ABC-related effects on their expression.
In contrast to the analyzed NF-κB target genes and the global differential gene expression between the two genetic groups, the differential expression of MALT1 itself and 32 other genes located on chromosome 18q was independent of the gene expression subtype. We focused our investigations on the MALT1 gene, but cannot exclude that any of the other 32 genes may play a pathogenetic role in DLBCL cases with a 18q21 gain. The pathogenetic mechanism of action of MALT1 and the other 18q genes in this context cannot be concluded from our data.
We found that gain of 18q21 including MALT1 and BCL2 was an independent negative prognostic factor in our cases of DLBCL, although the relatively small number of patients and the heterogeneous treatment protocols applied should be taken into account when interpreting this result. In the multivariate analysis including age, Ann Arbor stage, a break within the MYC gene, ABC/GCB gene expression subtype, and 18q21/MALT1 status, 18q21/MALT1 status was a statistically significant independent negative prognostic factor for overall survival (hazard ratio 2.13, p=0.03). In the univariate analysis, response to treatment was significantly reduced (p=0.014) and a trend to a shorter overall survival (p=0.09) was demonstrated.
Bea et al. found overrepresentation of 18q in DLBCL to be significantly associated with a primary nodal presentation, high serum lactate dehydrogenase, and high International Prognostic Index score.41 In addition, 18q gains were more commonly found in tumors with a higher number of chromosomal imbalances.41 Similar to the present study, patients with 18q gain had a significantly shorter cause-specific overall survival, higher risk of relapse, and a tendency to lack a complete response to treatment.41 This relationship between 18q gains and a more aggressive clinical behavior in DLBCL had not been well recognized previously. A recent study on follicular lymphoma identified a distinct clinical presentation in patients with 18q gains, although no differences were observed in outcome.42 In another study, dup(18q) was found to predict an unfavorable overall survival in follicular lymphomas with the t(14;18).43 In surgically resected t(11;18)-negative gastrointestinal DLBCL, trisomy of 18q21, as shown by interphase FISH, was also associated with a shorter overall survival.44
The various reports in the literature on the t(14;18) or BCL2 expression as a prognostic indicator in DLBCL have been inconclusive.45–49 Interestingly, a recent study showed that BCL2 protein expression represents a negative prognostic factor in the ABC-DLBCL subgroup, but not in the GCB subgroup.16 In the latter study, BCL2 protein expression in the ABC-DLBCL subgroup was associated with amplification and gain of the BCL2 gene,16 which was also demonstrated in the present study. Since 18q21/MALT1-positive cases are frequently found in the ABC-DLBCL subgroup, it cannot be excluded that this prognostically unfavorable gene expression subtype had an influence on the clinical course in our study.
In summary, our data demonstrate that a gain of 18q21 including MALT1 is associated with differential RNA expression of genes located on 18q, independently of the gene expression subtype, as well as BCL2 protein expression and clustering in the ABC subtype. The differences in expression of genes not located on 18q, including the analyzed NF-κB target genes, overlapped with the ABC/GCB, gene expression subtype. So, the individual influence of the gene expression signature and the 18q21/MALT1 status on gene expression cannot be viewed separately. Finally, 18q21/MALT1 gain may be a prognostic indicator in DLBCL, although this needs further evaluation in a prospective trial with homogeneously treated patients combining genomic analysis and gene expression profiling.
Footnotes
- Funding: this work was supported by the German Cancer Aid Network “Molecular Mechanisms in Malignant Lymphoma”, 70-3173-Tr3, and grant 106092 to JD, from the German Cancer Aid. A complete list of contributors of the Deutsche Krebshilfe Network Project “Molecular Mechanisms in Malignant Lymphomas” is provided in the Appendix.
- Further members of the Deutsche Krebshilfe Network Project “Molecular Mechanisms in Malignant Lymphomas“ who contributed to the present work are: TFE Barth, C Becher, H-W Bernd, C Bokemeyer, AC Feller, M-L Hansmann, L Harder, D Hasenclever, K Hinz, MD Kühn, P Lichter, JI Martin-Subero, H-K Müller-Hermelink, C Pott, M Rosolowski, B Stuerzenhofecker, M Szczepanowski, H Trautmann, and H-H Wacker.
- The online version of this article contains a supplementary appendix.
- Authorship and Disclosures JD and EMMP contributed equally to the study and should both be regarded as first authors; JD, EMMP, and RSi: conceived, designed and performed research, collected and analyzed data, wrote the manuscript; SB: analyzed and interpreted data, wrote the manuscript; SW and CS: performed research, collected and analyzed data, wrote the manuscript; HB: collected and analyzed data, wrote the manuscript; MH: provided material, collected and analyzed data, wrote the manuscript; WK and DL: performed research and collected data; AR: conceived and designed research, provided material, performed research, collected and analyzed data; EH: provided material, performed research, collected and analyzed data; GO and HS: provided material, collected data; SBG: provided material; PM: provided material, wrote the manuscript; ML and RSp: analyzed data, wrote the manuscript; LT: collected and analyzed data All authors approved the final version of the manuscript. The authors reported no potential conflicts of interest.
- Received July 30, 2007.
- Revision received December 7, 2007.
- Accepted December 31, 2007.
References
- A clinical evaluation of the International Lymphoma Study Group classification of non-Hodgkin’s lymphoma. Blood. 1997; 89:3909-18. Google Scholar
- Tumours of the Haematopoietic and Lymphoid Tissues. IARC Press: Lyon (France); 2001. Google Scholar
- Coiffier B. State-of-the-art therapeutics: diffuse large B-cell lymphoma. J Clin Oncol. 2005; 23:6387-93. Google Scholar
- A predictive model for aggressive non-Hodgkin’s lymphoma. N Engl J Med. 1993; 329:987-94. Google Scholar
- Alizadeh AA, Eisen MB, Davis RE, Ma C, Lossos IS, Rosenwald A. Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature. 2000; 403:503-11. Google Scholar
- Rosenwald A, Wright G, Chan WC, Connors JM, Campo E, Fisher RI. The use of molecular profiling to predict survival after chemotherapy for diffuse large-B-cell lymphoma. N Engl J Med. 2002; 346:1937-47. Google Scholar
- Rosenwald A, Wright G, Leroy K, Yu X, Gaulard P, Gascoyne RD. Molecular diagnosis of primary mediastinal B cell lymphoma identifies a clinically favorable subgroup of diffuse large B cell lymphoma related to Hodgkin lymphoma. J Exp Med. 2003; 198:851-62. Google Scholar
- Shipp MA, Ross KN, Tamayo P, Weng AP, Kutok JL, Aguiar RC. Diffuse large B-cell lymphoma outcome prediction by gene-expression profiling and supervised machine learning. Nat Med. 2002; 8:68-74. Google Scholar
- Monti S, Savage KJ, Kutok JL, Feuerhake F, Kurtin P, Mihm M. Molecular profiling of diffuse large B-cell lymphoma identifies robust subtypes including one characterized by host inflammatory response. Blood. 2005; 105:1851-61. Google Scholar
- Hummel M, Bentink S, Berger H, Klapper W, Wessendorf S, Barth TF. A biologic definition of Burkitt's lymphoma from transcriptional and genomic profiling. N Engl J Med. 2006; 354:2419-30. Google Scholar
- Huang JZ, Sanger WG, Greiner TC, Staudt LM, Weisenburger DD, Pickering DL. The t(14;18) defines a unique subset of diffuse large B-cell lymphoma with a germinal center B-cell gene expression profile. Blood. 2002; 99:2285-90. Google Scholar
- Lossos IS, Alizadeh AA, Eisen MB, Chan WC, Brown PO, Botstein D. Ongoing immunoglobulin somatic mutation in germinal center B cell-like but not in activated B cell-like diffuse large cell lymphomas. Proc Natl Acad Sci USA. 2000; 97:10209-13. Google Scholar
- Davis RE, Brown KD, Siebenlist U, Staudt LM. Constitutive nuclear factor kappaB activity is required for survival of activated B cell-like diffuse large B cell lymphoma cells. J Exp Med. 2001; 194:1861-74. Google Scholar
- Feuerhake F, Kutok J, Monti S, Chen W, LaCasce AS, Cattoretti G. NFkB activity, function, and target gene signatures in primary mediastinal large B-cell lymphoma and diffuse large B-cell lymphoma subtypes. Blood. 2005; 15:1392-9. Google Scholar
- Lam LT, Davis RE, Pierce J, Hepperle M, Xu Y, Hottelet M. Small molecule inhibitors of IkappaB kinase are selectively toxic for subgroups of diffuse large B-cell lymphoma defined by gene expression profiling. Clin Cancer Res. 2005; 11:28-40. Google Scholar
- Iqbal J, Neppalli VT, Wright G, Dave BJ, Horsman DE, Rosenwald A. BCL2 expression is a prognostic marker for the activated B-cell-like type of diffuse large B-cell lymphoma. J Clin Oncol. 2006; 24:961-8. Google Scholar
- Catz SD, Johnson JL. Transcriptional regulation of bcl-2 by nuclear factor kappa B and its significance in prostate cancer. Oncogene. 2001; 20:7342-51. Google Scholar
- Schwaenen C, Nessling M, Wessendorf S, Salvi T, Wrobel G, Radlwimmer B. Automated array-based genomic profiling in chronic lymphocytic leukemia: development of a clinical tool and discovery of recurrent genomic alterations. Proc Natl Acad Sci USA. 2004; 101:1039-44. Google Scholar
- Holzmann K, Kohlhammer H, Schwaenen C, Wessendorf S, Kestler HA, Schwoerer A. Genomic DNA-chip hybridization reveals a higher incidence of genomic amplifications in pancreatic cancer than conventional comparative genomic hybridization and leads to the identification of novel candidate genes. Cancer Res. 2004; 64:4428-33. Google Scholar
- Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J. Initial sequencing and analysis of the human genome. Nature. 2001; 409:860-921. Google Scholar
- Fiegler H, Carr P, Douglas EJ, Burford DC, Hunt S, Scott CE. DNA microarrays for comparative genomic hybridization based on DOP-PCR amplification of BAC and PAC clones. Genes Chromosomes Cancer. 2003; 36:361-74. Google Scholar
- Huver W, von Heydebreck A, Sultmann H, Poustka A, Vingron M. Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics. 2002; 18(Suppl 1):S96-104. Google Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003; 4:249-64. Google Scholar
- Mansmann U, Meister R. Testing differential gene expression in functional groups. Goeman's global test versus an ANCOVA approach. Methods Inf Med. 2005; 44:449-53. Google Scholar
- Smyth GK. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004; 3Google Scholar
- Wright G, Tan B, Rosenwald A, Hurt EH, Wiestner A, Staudt LM. A gene expression-based method to diagnose clinically distinct subgroups of diffuse large B cell lymphoma. Proc Natl Acad Sci USA. 2003; 100:9991-6. Google Scholar
- Shipp MA. Prognostic factors in aggressive non-Hodgkin's lymphoma: who has “high-risk” disease?. Blood. 1994; 83:1165-73. Google Scholar
- Ferch U, zum Büschenfelde CM, Gewies A, Wegener E, Rauser S, Peschel C. MALT1 directs B cell receptor-induced canonical nuclear factor-kappaB signaling selectively to the c-Rel subunit. Nat Immunol. 2007; 8:984-91. Google Scholar
- Akagi T, Motegi M, Tamura A, Suzuki R, Hosokawa Y, Suzuki H. A novel gene, MALT1 at 18q21, is involved in t(11;18) (q21;q21) found in low-grade B-cell lymphoma of mucosa-associated lymphoid tissue. Oncogene. 1999; 18:5785-94. Google Scholar
- Dierlamm J, Baens M, Wlodarska I, Stefanova-Ouzounova M, Hernandez JM, Hossfeld DK. The apoptosis inhibitor gene API2 and a novel 18q gene, MLT, are recurrently rearranged in the t(11;18) (q21;q21) associated with MALT lymphomas. Blood. 1999; 93:3601-9. Google Scholar
- Morgan JA, Yin Y, Borowsky AD, Kuo F, Nourmand N, Koontz JI. Breakpoints of the t(11;18) (q21;q21) in mucosa-associated lymphoid tissue (MALT) lymphoma lie within or near the previously undescribed gene MALT1 in chromosome 18. Cancer Res. 1999; 59:6205-13. Google Scholar
- Sanchez-Izquierdo D, Buchonnet G, Siebert R, Gascoyne RD, Climent J, Karran L. MALT1 is deregulated by both chromosomal translocation and amplification in B-cell non-Hodgkin lymphoma. Blood. 2003; 101:4539-46. Google Scholar
- Streubel B, Lamprecht A, Dierlamm J, Cerroni L, Stolte M, Ott G. t(14;18)(q32;q21) involving IGH and MALT1 is a frequent chromosomal aberration in MALT lymphoma. Blood. 2003; 101:2335-9. Google Scholar
- Uren AG, O'Rourke K, Aravind L, Pisabarro MT, Seshagiri S, Koonin EV. Identification of paracaspases and metacaspases: two ancient families of caspase-like proteins, one of which plays a key role in MALT lymphoma. Mol Cell. 2000; 6:961-7. Google Scholar
- Lucas PC, Yonezumi M, Inohara N, McAllister-Lucas LM, Abazeed ME, Chen FF. Bcl10 and MALT1, independent targets of chromosomal translocation in MALT lymphoma, cooperate in a novel NF-ÎB signaling pathway. J Biol Chem. 2001; 276:19012-9. Google Scholar
- Ho L, Davis RE, Conne B, Chappuis R, Berczy M, Mhawech P. MALT1 and the API2-MALT1 fusion act between CD40 and IKK and confer NF-kB-dependent proliferative advantage and resistance against FAS-induced cell death in B cells. Blood. 2005; 105:2891-9. Google Scholar
- Murga Penas EM, Hinz K, Röser K, Copie-Bergman C, Wlodarska I, Marynen P. Translocations t(11;18)(q21;q21) and t(14;18) (q32;q21) are the main chromosomal abnormalities involving MLT/MALT1 in MALT lymphomas. Leukemia. 2003; 17:2225-9. Google Scholar
- Reed JC. Bcl-2 family proteins: regulators of apoptosis and chemoresistance in hematologic malignancies. Semin Hematol. 1997; 34:9-19. Google Scholar
- Iqbal J, Sanger WG, Horsman DE, Rosenwald A, Pickering DL, Dave B. BCL2 translocation defines a unique tumor subset within the germinal center B-cell-like diffuse large B-cell lymphoma. Am J Pathol. 2004; 165:159-66. Google Scholar
- Bea S, Zettl A, Wright G, Salaverria I, Jehn P, Moreno V. Diffuse large B-cell lymphoma subgroups have distinct genetic profiles that influence tumor biology and improve gene-expression–based survival prediction. Blood. 2005; 106:3183-90. Google Scholar
- Bea S, Colomo L, Lopez-Guillermo A, Salaverria I, Puig X, Pinyol M. Clinicopathologic significance and prognostic value of chromosomal imbalances in diffuse large B-cell lymphomas. J Clin Oncol. 2004; 22:3498-506. Google Scholar
- Viardot A, Moller P, Hogel J, Werner K, Mechtersheimer G, Ho AD. Clinicopathologic correlations of genomic gains and losses i n follicular lymphoma. J Clin Oncol. 2002; 20:4523-30. Google Scholar
- Lestou VS, Gascoyne RD, Sehn L, Ludkovski O, Chhanabhai M, Klasa RJ. Multicolour fluorescence in situ hybridization analysis of t(14;18)-positive follicular lymphoma and correlation with gene expression data and clinical outcome. Br J Haematol. 2003; 122:745-59. Google Scholar
- Krugmann J, Tzankov A, Dirnhofer S, Fend F, Greil R, Siebert R. Unfavourable prognosis of patients with trisomy 18q21 detected by fluorescence in situ hybridisation in t(11;18) negative, surgically resected, gastrointestinal B cell lymphomas. J Clin Pathol. 2004; 57:360-4. Google Scholar
- Gascoyne RD, Adomat SA, Krajewski S, Krajewska M, Horsman DE, Tolcher AW. Prognostic significance of Bcl-2 protein expression and Bcl-2 gene rearrangement in diffuse aggressive non-Hodgkin’s lymphoma. Blood. 1997; 90:244-51. Google Scholar
- Kramer MH, Hermans J, Wijburg E, Philippo K, Geelen E, van Krieken JH. Clinical relevance of BCL2, BCL6, and MYC rearrangements in diffuse large B-cell lymphoma. Blood. 1998; 92:3152-62. Google Scholar
- Barrans SL, Carter I, Owen RG, Davies FE, Patmore RD, Haynes AP. Germinal center phenotype and bcl-2 expression combined with the International Prognostic Index improves patient risk stratification in diffuse large B-cell lymphoma. Blood. 2002; 99:1136-43. Google Scholar
- Barrans SL, Evans PA, O'Connor SJ, Owen RG, Morgan GJ, Jack AS. The t(14;18) is associated with germinal center-derived diffuse large B-cell lymphoma and is a strong predictor of outcome. Clin Cancer Res. 2003; 9:2133-9. Google Scholar
- Colomo L, Lopez-Guillermo A, Perales M, Rives S, Martinez A, Bosch F. Clinical impact of the differentiation profile assessed by immunophenotyping in patients with diffuse large B-cell lymphoma. Blood. 2003; 101:78-84. Google Scholar