Evidence for Distinct Pathomechanisms in Genetic Subgroups of Chronic Lymphocytic Leukemia Revealed by Quantitative Expression Analysis of C
http://www.100md.com
▲還散笫雖悝◎
the Department of Internal Medicine III, University of Ulm, Ulm
Department of Molecular Genetics and Department of Biostatistics, German Cancer Research Center, Heidelberg, Germany
Department of Internal Medicine I, Medical University of Vienna, Vienna, Austria
ABSTRACT
PURPOSE: In patients with chronic lymphocytic leukemia (CLL), the VH mutation status and genomic aberrations (13q每, +12q, 11q每, 17p每) identify distinct prognostic subgroups. The aim was to elucidate biologic mechanisms through which these genetic markers may exert their pathogenic influence.
PATIENTS AND METHODS: Twenty-four genes involved in apoptosis, cell cycle, B-cell activation, and B-cell receptor (BCR) signaling were analyzed by real-time quantitative reverse transcription polymerase chain reaction (RQ-PCR) in 82 CLL cases constituting prototypic genetic CLL subgroups as defined by the VH mutation status and the genomic aberrations 13q每, +12, 11q每, and 17p每.
RESULTS: The VH mutation subgroups were characterized by a differential expression of the BCR associated genes ZAP70 and PI3K. Among the subgroups defined by genomic aberrations, there was a deregulation of candidate genes from the affected critical genomic regions such as CDK4 (up), ATM (down), and TP53 (down) in the groups +12, 11q每, and 17p每, respectively. Additionally, the genomic subgroups were characterized by a significant deregulation of cell cycle and apoptosis regulators: AKT (up) in 13q, E2F1 (up) in +12, MYC (up) and BCL-2 (down) in 17p每, and CCND3 (down) in 11q每 as well as 17p每. The 17p每 subgroup showed an additional down-regulation of BCR-associated genes such as SYK and PI3K.
CONCLUSION: The characteristic gene expression patterns observed implicate a differential regulation of cell cycle, apoptosis, and BCR signaling in the genetic subgroups illustrating distinct pathomechanisms and are evidence for a gene dosage effect being operative in CLL. These findings link the biologic diversity and clinical heterogeneity of CLL.
INTRODUCTION
Chronic lymphocytic leukemia (CLL) is characterized by a variable clinical course, with survival times ranging from months to decades.1,2 This clinical heterogeneity has been shown to be related to inherent genetic features of the CLL cells such as the VH mutation status and genomic aberrations. While CLL with mutated VH displays a good prognosis, VH unmutated CLL is characterized by a poor outcome,3,4 suggesting that the different maturation stages according to the VH mutation status may define distinct CLL pathomechanisms.5 However, by genome-wide gene expression analysis using microarrays, the transcriptional profile of the VH subgroups could not be distinguished in unsupervised analyses, and only a small number of genes were found to be differentially expressed when comparing the VH subgroups by supervised analyses.6,7 Among these, ZAP70 had a prominent role, and further studies showed that B-cell receptor (BCR) crosslinking on ZAP70-positive CLL cells led to increased tyrosine phosphorylation (eg, of p72 [Syk]), indicating an increased activation after BCR stimulation as a pathomechanism.8
Additionally, genomic aberrations, which have an independent prognostic impact and may be transforming events in CLL, were not included in these studies. The most common genomic aberrations are losses at the chromosomal bands 13q14 (13q每: 55%), 11q22 to q23 (11q每: 18%), 17p13 (17p每: 7%), and gains at chromosome 12 (+12: 16%).9 17p每 and 11q每 identify poor prognostic subgroups in multivariate analysis.9-12 Recurrent genomic aberrations point to the loci of pathogenic genes such as ATM at 11q22-23 and TP53 at 17p13. However, the expression pattern of these and other candidate genes is so far unresolved in the subgroups of CLL defined by genomic aberrations.
The aim of the present study was to elucidate pathogenic events in pathways of oncogenic potential in the genetic subgroups of CLL. To investigate candidate gene expression of genetically distinct CLL subtypes, we built prototypic subgroups composed of cases defined by the VH mutation status and genomic aberrations, ie, VH mutated (VH MUT), VH unmutated (VH UM), normal karyotype by FISH (normal), 13q每 as sole abnormality (13q每), 11q每, +12, and 17p每. To obtain a highly sensitive quantitation of gene expression, we used real-time quantitative reverse transcription polymerase chain reaction (RQ-PCR) and focused on a selected set of functionally relevant genes playing central roles in the regulation of apoptosis, cell cycle, B-cell activation, and BCR signaling, namely AICL (CLECSF2), AID (AICDA), AKT1, ATM, BAX, BCL2, CCND1, CCND2, CCND3, CDK4, MYC, E2F1, P21 (CDKN1A), P27 (CDKN1B), PI3K (PIK3CB), MDM2, MCL1, NFKB1, Survivin (BIRC5), SYK, TCL (TCL1A), TP53, TRADD, and ZAP70.
PATIENTS AND METHODS
Patients
Peripheral blood was obtained for diagnostic procedures from all cases after informed consent and the studies were institutional review board每approved. Samples from 82 patients with CLL diagnosed according to established morphologic and immunophenotypic criteria were included. Staging information was complete for 76 patients; 30 patients were Binet A; 26, Binet B; and 20, Binet C. Median age was 57 years (range, 38 to 84 years), and 56% were male. All samples were from the same date as the corresponding VH and fluorescent in situ hybridization (FISH) analyses. Sixty-seven samples were collected at time of diagnosis (15 during the course of the disease), and all except seven patients were untreated. In these cases, mostly patients with refractory disease, one to three previous therapies had been administered, and the samples were collected in active recurrent disease after at least 2 treatment-free months. Clinical follow-up for survival analysis was available for 60 patients. Translocation t(11;14)每positive cases were not included.
Genetic Analysis
FISH analysis and VH sequencing was performed in all cases as previously described.9,12,13 A germline homology of 98% was used as the cutoff between VH mutated and VH unmutated cases.
Genetic Subgroups
To avoid the analysis of subclones, a high tumor load as defined by a high percentage of cells carrying the respective genomic aberration detected by FISH was the prerequisite for inclusion in the study. Median values of aberrant cells were 84% in 13q每 (range, 75% to 93%), 80% in +12 (range, 71% to 85%), 90% in 11q每 (range, 77% to 97%), and 89% in 17p每 (range, 73% to 97%). In cases with a normal karyotype by FISH, leukocyte and lymphocyte counts, in combination with CD5/CD19 expression in fluorescence-assisted cell sorter analysis, were used to determine the tumor load (median, 93%; range, 72% to 98%). The following subgroups were studied according to their genetic characteristics as defined by VH mutation status and genomic aberrations (Fig 1): VH MUT and VH UM (each equally divided into 13q- and normal); 13q每 (as single abnormality) and normal (no aberration by FISH; each equally divided into VH MUT and VH UM); +12 (and absence of 11q每 or 17p每), 11q每 (and absence of +12 or 17p每), and 17p每 (and absence of +12 or 11q-), that is, there was no occurrence of a +12, 11q每, or 17p每 outside the subgroup defined by these aberrations. If existing, additional aberrations in the groups 11q每, +12, and 17p每 consisted mainly of 13q每. In rare cases, other abnormalities (6q每, four cases; 14q每, three cases; +3q, one case; and +8q, one case) were present. The subgroups 11q每, +12, and 17p每 consisted mainly of VH UM cases: 11 of 15 with a +12, 16 of 16 with a 11q每, and 14 of 17 with a 17p每, allowing the additional comparison of homogeneously VH UM cases with +12 (+12[UM]: n = 11), 11q- (11q每[UM]: n = 16), or 17p每 (17p每[UM]: n = 14) versus VH UM cases without these aberrations (VH UM). Median VH homology in the VH MUT and VH UM subgroups was 92.8% (range, 88.4% to 95.7%) and 99.7% (range, 98.1% to 100%), respectively. For standardization purposes, CD19-positive fractions purified by magnet-activated cell separation CD19 MicroBeads (Miltenyi Biotech, Bergisch Gladbach, Germany) from buffy coat preparations of four healthy individuals were used.
Sample and RNA Preparation
The procedures were done as previously reported,14 including a DNAse I digestion of total RNA to avoid contamination with genomic DNA.
RQ-PCR
RQ-PCR analysis was carried out as previously described14 using the ABI Prism 7700 Sequence Detector (Applied Biosystems, Weiterstadt, Germany). Quantitative assessment of DNA amplification was achieved either through SYBR Green or TaqMan probes according to the manufacturer's manuals (SYBR Green PCR Core Reagents, TaqMan Universal Master Mix; Applied Biosystems). For SYBR Green reactions, a melting curve of the final RQ-PCR products was generated to avoid nonspecific fluorescence derived from side products such as primer-dimers (ABI Prism Dissociation Curve Software; PE Applied Biosystems, Foster City, CA). Primers for SYBR Green detection were designed using Primer Express software (Applied Biosystems) and are listed in Table 1. Primers for the following genes were used as previously described: BAX, CCND3, CDK4, E2F1, CMYC, P21, P27, TP53, lamin B1 (LMNB1), and phosphoglycerate kinase 1 (PGK1)14; CCND115; and Survivin.16 For beta actin (ACTB) PreDevelopedAssay-reagents, and for AID as well as ATM, Assay-on-Demand Gene Expression Probes were used (all supplied by Applied Biosystems).
Data Normalization
Eight serial dilutions of cDNA from tumor cell lines were used to obtain a calibration graph, which was measured in every single RQ-PCR run. In general, the B-cell CLL cell line EHEB was used for calibration. The T-cell line Jurkat and the Burkitt cell line Namalva were used for calibration of ZAP70 and AID expression levels, respectively, due to absent expression in EHEB. The average value of three independent endogenous control amplicons (PGK1, LMNB1, ACTB) were used for normalization.17
Control Experiments
Experiments were repeated for TP53, MYC, CDK4, BCL2, CCND3, E2F1, and for all three housekeeping genes to obtain an independent data set.
Statistical Analysis
For gene screening, we computed unadjusted and permutation-adjusted P values using standardized Wilcoxon rank sum test statistics. Permutation-adjusted P values were computed for the minP step-down multiple testing procedure described in Westfall and Young.18 This procedure provides strong control for the family-wise error rate and also incorporates the joint dependence structure between the test statistics. The number of permutations used was 10,000. In addition, we used summary measures based on receiver operating characteristic (ROC) curves as described by Pepe et al.19 The area under the ROC curve, which corresponds to the Wilcoxon rank sum test statistic, scores for discrimination ability (and equals 0.5 for a random classifier). Fold changes of expression values between compared subgroups were estimated by "back-transformation" of the Hodges-Lehmann estimate of the difference of the location parameters of the log-transformed expression data. Gene correlations were estimated by Spearman's statistic and plotted as a red and blue color image of the correlation matrix. Increasingly positive correlations are represented with reds of increasing intensity, and increasingly negative correlations are represented with blues of increasing intensity. To present correlation matrices we started with the pair of genes with the largest positive coefficient of pairwise correlation. The next gene selected is the one whose expression values show the strongest positive correlation with the mean expression of the two previously selected genes. This procedure is continued by selecting the gene whose expression values have the strongest positive correlation with the mean expression of all the genes previously selected.
RESULTS
Quantitative expression levels of 24 genes (listed in Table 2 with their chromosomal position) with central roles in the regulation of apoptosis, cell cycle, BCR signaling, and B-cell activation were analyzed by RQ-PCR in 82 cases constituting the prototypic genetic CLL subgroups (VH MUT, VH UM, 13q每, normal, +12, 11q每, and 17p每). A survival probability analysis was performed to test whether the selected cases were a representative collection with regard to the clinical impact of the genetic markers to be analyzed (Fig 2). The individual subgroups VH MUT, VH UM, +12, 11q每, and 17p每 showed the expected survival probabilities, demonstrating that the selected cases represent a valid collection for this study.
VH Mutation Subgroups
The VH MUT versus VH UM analysis was restricted to the set of normal and 13q每 cases to avoid interference of potentially distinct expression signatures of +12, 11q每, and 17p每 cases. Only two of the 24 genes analyzed showed significant differential distributions of expression levels between the VH mutation subgroups (summarized in Table 3): ZAP70 levels were 3.24-fold higher in VH UM than in VH MUT cases (Fig 3A), whereas PI3K showed lower expression levels in VH UM cases (0.68-fold). CD3 RNA expression levels, used as a calibrator for T-cell contamination, were comparable between the VH subgroups, and a correction for CD3 expression did not change the results for ZAP70 significantly (data not shown). The CCND1/CCND3 ratio was higher in VH MUT. When comparing VH mutated and VH unmutated cases in the overall cohort (n = 82; Table 4) including +12, 11q每, and 17p每 cases, similar results were obtained for ZAP70 and PI3K, but in addition, a down-regulation of ATM and CCND1 was evident, probably as a consequence of the inclusion of 11q每 and 17p每 cases, respectively (see Results section Genomic Aberrations).
Due to the potential pathogenic impact of ZAP70 expression in CLL, an additional comparison of ZAP70-high and ZAP70-low expressing cases was carried out. Therefore, 21 cases from the quartile with the highest ZAP70 expression were compared with 21 cases from the quartile with the lowest expression. In the ZAP70-high group, an upregulation of SYK (P = .02, 1.59-fold) and a down-regulation of PI3K (P = .04, 0.71-fold) were identified.
Genomic Aberrations
When comparing the different subgroups defined by genomic aberrations, a majority of the investigated genes (17 of 24) showed a differential expression pattern.
13q每 (single) was compared with normal cases (Table 5). The most characteristic finding was a low-level but highly differential overexpression of AKT (P = .001) in 13q每 cases (Fig 3B). To explore whether the differential expression of AKT was restricted to cases with a 13q每 single aberration (ie, no additional abnormality) or a general phenomenon of 13q deletion, we compared all cases with a 13q deletion (n = 42) against all cases without a 13q deletion (n = 40), irrespective of additional aberrations. The higher expression of AKT was also observed when comparing the 13q deletion (irrespective of additional aberrations) with the non-13q deletion subgroup (Fig 4).
The subgroups +12, 11q每, and 17p每 were individually compared with the cohort of normal and 13q- cases (n = 34, Fig 1). This comparison was chosen due to the occurrence of an additional 13q每 in approximately half of the cases with +12, 11q每, and 17p每. In addition, to investigate expression characteristics without interferences of VH mutation status, a restricted analysis was carried out including only VH UM cases. In this comparison, the homogeneously unmutated subgroups +12(UM), 11q每(UM), and 17p每(UM) were compared individually with the VH UM subgroup without these aberrations.
When comparing the +12 subgroup with the normal/13q每 cohort (Table 6), characteristic expression patterns were observed for E2F1 (P < .001), BAX (P = .002), P27 (P = .022), and CDK4 (P = .050). Notably, all of these differentially expressed genes showed an overexpression (1.33- to 1.73-fold) in +12. BCL2/BAX and CCND1/CCND3 ratios were significantly lower in +12. The restricted analysis of +12(UM) versus VH UM (Table 7) was statistically limited by the small number of +12(UM) cases (n = 11). However, comparable alterations of candidate gene expression levels were identified, confirming the upregulation of E2F1.
In the 11q每 subgroup (Table 8), the most differentially expressed genes were ATM (Fig 3C) and CCND3, both with a lower expression in 11q每 as compared with 13q每/normal (P < .001, 0.43-fold and P = .015, 0.66-fold, respectively). When the comparison was restricted to 11q每(UM) versus VH UM (Table 9), the lower ATM expression remained the most prominent finding (P < .001), but additional trends for upregulation of MCL and NFKB and down-regulation of SYK were noted.
The 17p每 subgroup displayed the highest number of differentially expressed genes affecting all investigated pathways (Tables 10 and 11). In comparison to normal/13q每, the 17p每 cases showed a significantly lower expression of CCND3 (P < .001), BCL2 (P = .001), SYK (P = .004), TP53 (P = .007; Fig 3D), ATM (P = .005), TCL (P = .006), PI3K (P = .021), CCND1 (P = .030), and AID (0.013) with a range of underexpression between 0.31- and 0.66-fold. Significant overexpression was observed for P21 (1.94-fold), MYC (2.18-fold), and AICL (3.57-fold). BCL2/BAX ratios were lower in the 17p每 subgroup. The comparison of 17p每(UM) versus VH UM produced similar results. In this analysis, TP53 was the most significantly down-regulated gene (P = .003).
Figure 5 gives a comparison of median expression levels of selected candidate genes involved in apoptosis, cell cycle regulation and BCR signaling in the genetic subgroups VH MUT, VH UM, +12, 11q每, and 17p每. ZAP70 expression levels were highest in VH UM, but lower in +12, 11q每, and 17p每 cases, though these were predominantly VH unmutated. Similarly, SYK expression was lower in +12, 11q每, and 17p每 as compared with VH UM. While ATM expression was lowest in 11q每 cases, TP53 levels were lowest in 17p每 cases. Both subgroups harbor a genomic loss of the respective gene locus. Surprisingly, BCL2 showed the highest expression in the VH MUT subgroup but the lowest expression in the poor-risk 17p每 subgroup. BAX expression was homogeneous except for an overexpression in +12. Whereas CCND3 expression was characteristically low in 11q每 and 17p每, MYC was overexpressed in 17p每.
Candidate genes that discriminated the genomic subgroups +12, 11q每, and 17p每 from normal/13q每 (P .05) were used for an outline of the major pathways in which the investigated genes are involved and the mode of their deregulation in the genetic CLL subgroups (Fig. 6A to C).
Gene Expression Inter-Relations
For a general comparison of subgroup-specific expression inter-relations, all genes were clustered hierarchically by their strength of correlation (Figs 7A to E), starting with the pair of genes with the largest positive coefficient of pairwise correlation (see Patients and Methods). Using this algorithm, highly distinct, subgroup-specific expression inter-relations can be illustrated. Whereas in VH MUT the strongest cluster included positive correlations between MCL, NFKB, and P27, the cluster consisted of CCND2, TCL, and BCL2 in VH UM (Figs 7A and B). Although the subgroups +12, 11q每, and 17p每 were predominantly composed of VH unmutated cases, their clustering was distinct for each of the subgroups and differed from the VH UM group. The strongest cluster in +12 consisted of PI3K, ATM, and TCL, in contrast to MCL, P21, and MDM2 in 11q每, and BAX and TRADD in 17p (Fig 7C每E).
Control Experiments
To ensure the reproducibility of the results, experiments were repeated in replicate with independent measurements of six candidate genes and all housekeeping genes. Highly comparable results were obtained. Distributions with initial P values less than .05 could be reproduced in all cases, and no additional significant results were obtained (data not shown).
DISCUSSION
Although genetic factors have a strong impact on the clinical course of CLL, little is known about the pathomechanisms of the genetic subgroups. Despite their highly distinct clinical behavior it was not possible to separate the CLL subgroups with mutated and unmutated VH genes by unsupervised cluster analysis of global gene expression data derived from microarrays.6,7 However, in these studies, genomic aberrations were not included, and the analyses were not performed by a highly quantitative method such as RQ-PCR, possibly obscuring subtle expression differences in a wealth of data points. In contrast to previous expression analyses, we used unmanipulated cell samples (ie, no CD19+ selection), but selected for a high content of genetically defined tumor cells. We thereby gained a "genetic" tumor purity of approximately 85% as compared with greater than 95% "B-cell purity" by CD19 selection. A disadvantage of this approach was a potentially remaining influence of non-CLL cells in the sample. This applies particularly for ZAP70, which is strongly expressed in T-cells. Since the current study focused on subgroup comparisons, this problem may have been reduced by the balanced content of genetically defined CLL cells and contaminating T-cells (monitored by CD3) in the different subgroups.
In line with pivotal studies6,7 only few differences in gene expression were observed between the VH mutation subgroups. ZAP70 was previously shown to be related to increased BCR signaling in CLL.8 In line with this, we identified elevated levels of SYK, a ZAP70-related kinase, in cases with high ZAP70 expression, arguing for a synergistic role of ZAP70 and SYK in initial BCR signaling. In contrast, PI3K expression, a gene known to be involved in BCR signaling by linking initial receptor signals from tyrosine kinases to downstream effectors,20-23 was lower in VH UM CLL. Similarly, PI3K was identified between the ZAP70 subgroups with a lower expression in ZAP70 high expressing cases. In contrast to Chen et al,8 these findings suggest an inactivity or suppression of physiologic downstream BCR cascades as a consequence of ZAP70 overexpression, potentially favoring atypical cellular responses or alternative pathways following BCR stimulation.
In CLL characterized by genomic aberrations, multiple gene deregulations affecting different pathways were observed, suggesting a combined model of cell-cycle deregulation, impairment of apoptosis, and dysfunctional cellular responsiveness for these subgroups. Importantly, most deregulations between genomic subgroups were reidentified when comparing the genomic subgroups in purely VH unmutated cases, clearly demonstrating the additional impact of genomic aberrations on the gene expression pattern independently of the VH mutation status.
Although 13q deletion is the most common genomic abnormality in CLL, its biologic role is still largely unresolved.24,25 In 13q每, the most characteristic finding was a low-level overexpression of AKT, which was also observed in the overall cohort of cases with 13q deletion, despite the inclusion of heterogeneous cases with secondary genetic events (11q每, +12, 17p每). Therefore, AKT overexpression seems to be a general phenomenon of 13q deletion potentially providing a new pathomechanism for this genetic abnormality, and not a phenomenon linked only to the favorable 13q每 (single) group. AKT can be activated by the BCR and promotes B-cell survival by inactivating components of the apoptotic machinery.26,27 In CLL cells, CD40-crosslinking induced cell proliferation, which was associated with an increase of AKT.28 Additionally, AKT activation was shown to protect CLL cells from chlorambucil- and radiation-induced apoptosis.29 This might explain a survival advantage for tumor clones with enhanced AKT expression, thereby favoring cells with a 13q deletion. Whether AKT overexpression could result from a deregulation of candidate genes from the critical genomic region on 13q14, has to be assessed.
Remarkably, in +12, all of the differentially expressed genes (E2F1, BAX, P27, CDK4) were overexpressed. Overexpression of CDK4 and P27 might be explained by a direct gene dosage-dependant upregulation due to their location on chromosome 12. The transcription factor E2F1 is a critical determinant of the G1/S-phase transition during cell cycle, leading to the transcription of S-phase proteins.30-32 Overexpression of E2F1 could be correlated with proliferation in human cancers and B-cell lymphomas.33,34 Therefore, E2F1 overexpression in +12 suggests enhanced proliferative activity as a potential pathomechanism in the evolution of this genetic abnormality. In addition, E2F1 activity is known to be regulated through the action of G1 cyclins and associated kinases such as CDK4.30-32 Therefore, overexpression of CDK4 might pathomechanistically contribute to E2F1 expression in +12.
Deletion in chromosome bands 11q22-23 is associated with extensive lymph node involvement and poor survival in CLL.9 ATM is located within the minimally deleted region and coordinates responses to DNA damage and is able to induce repair mechanisms and apoptosis by activating the TP53 pathway,35 suggesting that ATM might be acting as a tumor suppressor gene. However, only about one third of the CLL cases with 11q每 exhibit ATM mutations of the second allele,36,37 implicating the requirement of additional ATM inactivation mechanisms in CLL with 11q每. In this study, we demonstrate that the majority of 11q每 cases (13 of 16) exhibit reduced ATM transcript levels, pointing to a gene dosage-dependant down-regulation. Pettitt et al showed that p53 dysfunction in CLL may occur independently of TP53 mutations.38 Interestingly, these cases were associated with ATM mutations and exhibited reduced ATM protein expression. Therefore, down-regulation of ATM transcription in 11q- CLL may be a pathogenic factor contributing to dysfunctional apoptosis and dismal clinical outcome.
The largest number of deregulated genes was detected in the 17p- subgroup possibly reflecting the distinctly aggressive biology of this subgroup. This finding is in contrast to Stankovic et al39 reporting an indistinguishable native gene expression pattern when comparing wild-type and TP53-mutant CLL. The limited number of cases studied by Stankovic et al and the higher sensitivity of RQ-PCR in comparison to the microarray approach may account for this discrepancy. The characteristic two-fold reduction of TP53 transcript observed in our series indicates a gene-dosage effect as additional or synergistic mechanism contributing to TP53 dysfunction in CLL. Of note, we detected significantly reduced ATM levels in cases with 17p每, a finding that might contribute synergistically to apoptosis impairment in the 17p每 subgroup. The consensual down-regulation of ATM in the 11q每 and 17p每 subgroup argues for related pathogenic pathways involving a common mechanism of ATM inactivation as previously suggested.38,39 Overexpression of P21 despite low TP53 expression in 17p每 appears unexpected and may indicate TP53-independent counter-regulation of P21 in 17p每 CLL.
The finding of particularly low BCL2 levels in 17p每, the subgroup with the worst clinical course, was surprising due to the known antiapoptotic function of BCL2 and the association of high BCL2/BAX ratios with aggressive disease.40-42 This finding argues against a relevant role of BCL2 overexpression or elevated BCL2/BAX ratios in the pathogenesis of the genetic CLL subgroups or in CLL prognosis.
In addition to the down-regulation of apoptosis-associated genes in 17p每, there was a profound deregulation of cell-cycle regulators, including CCND3, CCND1, and MYC. CCND1 and, more pronounced, CCND3, showed reduced transcript levels in 17p每. However, as a consequence of this consensual down-regulation, the relationship between the transcript levels of both genes (CCND1/CCND3 ratio) remained largely unchanged. In contrast to the cyclins, MYC showed a differential overexpression in 17p每. MYC is known as a potent activator of cell proliferation,43 and MYC expression has been suggested as a marker for poor prognosis in B-cell lymphomas,44-46 implicating enhanced proliferative activity as an additional pathomechanism of the 17p每 subgroup.
A common phenomenon observed in the subgroups +12, 11q每, and 17p每 was the deregulation of genes from critical genomic regions implicating a gene-dosage effect being operative in CLL, which is in line with a recent report.47 In addition to the differential expression of individual genes described in this work, cluster analysis of gene expression correlations lead to subgroup-specific patterns reinforcing the concept of pathomechanisms distinct for the genetic CLL subgroups and thereby linking the biologic diversity and clinical heterogeneity of CLL.
Authors' Disclosures of Potential Conflicts of Interest
The authors indicated no potential conflicts of interest.
NOTES
Supported by the DFG (Sti 296/1-1), Deutsche Krebshilfe (70-3183-Li1), Wilhelm Sander Stiftung (2002.095.1), Jubilumsfonds of the Austrian National Bank (No. 9964; U.J.).
Presented at the Annual Meeting of the American Society of Hematology, San Diego, CA, 2003.
Authors' disclosures of potential conflicts of interest are found at the end of this article.
REFERENCES
Rozman C, Montserrat E: Chronic lymphocytic leukemia. N Engl J Med 333:1052-1057, 1995
Zwiebel JA, Cheson BD: Chronic lymphocytic leukemia: Staging and prognostic factors. Semin Oncol 25:42-59, 1998
Hamblin TJ, Davis Z, Gardiner A, et al: Unmutated Ig V(H) genes are associated with a more aggressive form of chronic lymphocytic leukemia. Blood 94:1848-1854, 1999
Damle RN, Wasil T, Fais F, et al: Ig V gene mutation status and CD38 expression as novel prognostic indicators in chronic lymphocytic leukemia. Blood 94:1840-1847, 1999
K邦ppers R, Klein U, Hansmann ML, et al: Cellular origin of human B-cell lymphomas. N Engl J Med 341:1520-1529, 1999
Rosenwald A, Alizadeh AA, Widhopf G, et al: Relation of gene expression phenotype to immunoglobulin mutation genotype in B cell chronic lymphocytic leukemia. J Exp Med 194:1639-1647, 2001
Klein U, Tu Y, Stolovitzky GA, et al: Gene expression profiling of B cell chronic lymphocytic leukemia reveals a homogeneous phenotype related to memory B cells. J Exp Med 194:1625-1638, 2001
Chen L, Widhopf G, Huynh L, et al: Expression of ZAP-70 is associated with increased B-cell receptor signaling in chronic lymphocytic leukemia. Blood 100:4609-4614, 2002
Dhner H, Stilgenbauer S, Benner A, et al: Genomic aberrations and survival in chronic lymphocytic leukemia. N Engl J Med 343:1910-1918, 2000
Geisler CH, Philip P, Egelund, et al: In B-cell chronic lymphocytic leukaemia chromosome 17 abnormalities and not trisomy 12 are the single most important cytogenetic abnormalities for the prognosis: A cytogenetic and immunophenotypic study of 480 unselected newly diagnosed patients. Leuk Res 21:1011-1023, 1997
Oscier DG, Gardiner AC, Mould SJ, Glide S, Davis ZA, Ibbotson RE, Corcoran MM, Chapman RM, Thomas PW, Copplestone JA, Orchard JA, Hamblin TJ: Multivariate analysis of prognostic factors in CLL: Clinical stage, IGVH gene mutational status, and loss or mutation of the p53 gene are independent prognostic factors. Blood 100:1177-1784, 2002
Krber A, Seiler T, Benner A, et al: V(H) mutation status, CD38 expression level, genomic aberrations, and survival in chronic lymphocytic leukemia. Blood 100:1410-1416, 2002
Kienle D, Krber A, Katzenberger T, et al: VH mutation status and VDJ rearrangement structure in mantle cell lymphoma: Correlation with genomic aberrations, clinical characteristics, and outcome. Blood 102:3003-3009, 2003
Korz C, Pscherer A, Benner A, et al: Evidence for distinct pathomechanisms in B-cell chronic lymphocytic leukemia and mantle cell lymphoma by quantitative expression analysis of cell cycle and apoptosis-associated genes. Blood 99:4554-4561, 2002
Bijwaard KE, Aguilera NS, Monczak Y, et al: Quantitative real-time reverse transcription-PCR assay for cyclin D1 expression: Utility in the diagnosis of mantle cell lymphoma. Clin Chem 47:195-201, 2001
Schmidt SM, Schag K, M邦ller MR, et al: Survivin is a shared tumor-associated antigen expressed in a broad variety of malignancies and recognized by specific cytotoxic T cells. Blood 102:571-576, 2003
Mertens D, Wolf S, Schroeter P, et al: Down-regulation of candidate tumor suppressor genes within chromosome band 13q14.3 is independent of the DNA methylation pattern in B-cell chronic lymphocytic leukemia. Blood 99:4116-4121, 2002
Westfall PH, Young SS: Resampling-Based Multiple Testing: Examples and Methods for P-Value Adjustment. New York, NY, John Wiley & Sons, 1993, p 66
Pepe MS, Longton G, Anderson G, et al: Selecting differentially expressed genes from microarray experiments. Biometrics 59:133-142, 2003
Beitz LO, Fruman DA, Kurosaki T, et al: SYK is upstream of phosphoinositide 3-kinase in B cell receptor signaling. J Biol Chem 274:32662-32666, 1999
Bone H, Williams NA: Antigen-receptor cross-linking and lipopolysaccharide trigger distinct phosphoinositide 3-kinase-dependent pathways to NF-kappa B activation in primary B cells. Int Immunol 13:807-816, 2001
Bernal A, Pastore RD, Asgary Z, et al: Survival of leukemic B cells promoted by engagement of the antigen receptor. Blood 98:3050-3057, 2001
Glassford J, Soeiro I, Skarell SM, et al: BCR targets cyclin D2 via Btk and the p85alpha subunit of PI3-K to induce cell cycle progression in primary mouse B cells. Oncogene 22:2248-2259, 2003
Calin GA, Dumitru CD, Shimizu M, et al: Frequent deletions and down-regulation of micro- RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A 99:15524-15529, 2002
Migliazza A, Bosch F, Komatsu H, et al: Nucleotide sequence, transcription map, and mutation analysis of the 13q14 chromosomal region deleted in B-cell chronic lymphocytic leukemia. Blood 97:2098-2104, 2001
Datta SR, Brunet A, Greenberg ME: Cellular survival: A play in three Akts. Genes Dev 13:2905-2927, 1999
Pogue SL, Kurosaki T, Bolen J, et al: B cell antigen receptor-induced activation of Akt promotes B cell survival and is dependent on Syk kinase. J Immunol 165:1300-1306, 2000
Cuni S, Perez-Aciego P, Perez-Chacon G, et al: A sustained activation of PI3K/NF-kappaB pathway is critical for the survival of chronic lymphocytic leukemia B cells. Leukemia 18:1391-1400, 2004
Jones DT, Ganeshaguru K, Anderson RJ, et al: Albumin activates the AKT signaling pathway and protects B-chronic lymphocytic leukemia cells from chlorambucil- and radiation-induced apoptosis. Blood 101:3174-3180, 2003
Bartek J, Lukas J: Pathways governing G1/S transition and their response to DNA damage. FEBS Lett 490:117-122, 2001
Sherr CJ, Roberts JM. CDK inhibitors: Positive and negative regulators of G1-phase progression. Genes Dev 13:1501-1512, 1999
Van der Sman J, Thomas NS, Lam EW: Modulation of E2F complexes during G0 to S phase transition in human primary B-lymphocytes. J Biol Chem 274:12009-12016, 1999
Zacharatos P, Kotsinas A, Evangelou K, et al: Distinct expression patterns of the transcription factor E2F-1 in relation to tumour growth parameters in common human carcinomas. J Pathol 203:744-753, 2004
Chiarle R, Budel LM, Skolnik J, et al: Increased proteasome degradation of cyclin-dependent kinase inhibitor p27 is associated with a decreased overall survival in mantle cell lymphoma. Blood 95:619-626, 2000
Waterman MJ, Stavridi ES, Waterman JL, et al: ATM-dependent activation of p53 involves dephosphorylation and association with 14-3-3 proteins. Nat Genet 19:175-178, 1998
Schaffner C, Stilgenbauer S, Rappold GA, et al: Somatic ATM mutations indicate a pathogenic role of ATM in B-cell chronic lymphocytic leukemia. Blood 94:748-753, 1999
Austen B, Fegan C, Moss P, et al ATM mutations define a subset of B-CLL tumors that only partially correlates with 11q status and is associated with impaired patient survival. Blood 104, 2004 (abstr 774)
Pettitt AR, Sherrington PD, Stewart G, et al: p53 dysfunction in B-cell chronic lymphocytic leukemia: Inactivation of ATM as an alternative to TP53 mutation. Blood 98:814-822, 2001
Stankovic T, Hubank M, Cronin D, et al: Microarray analysis reveals that TP53- and ATM-mutant B-CLLs share a defect in activating proapoptotic responses after DNA damage but are distinguished by major differences in activating prosurvival responses. Blood 103:291-300, 2004
Robertson LE, Plunkett W, McConnell K, et al: Bcl-2 expression in chronic lymphocytic leukemia and its correlation with the induction of apoptosis and clinical outcome. Leukemia 10:456-459, 1996
Aguilar-Santelises M, Rottenberg ME, Lewin N, et al: Bcl-2, Bax and p53 expression in B-CLL in relation to in vitro survival and clinical progression. Int J Cancer 69:114-119, 1996
Molica S, Dattilo A, Giulino C, Levato D, Levato L: Increased bcl-2/bax ratio in B-cell chronic lymphocytic leukemia is associated with a progressive pattern of disease. Haematologica 83:1122-1124, 1998
Bouchard C, Staller P, Eilers M: Control of cell proliferation by Myc. Trends Cell Biol 8:202-206, 1998
Nagy B, Lundan T, Larramendy ML, et al: Abnormal expression of apoptosis-related genes in haematological malignancies: Overexpression of MYC is poor prognostic sign in mantle cell lymphoma. Br J Haematol 120:434-441, 2003
Chang CC, Kampalath B, Schultz C, et al: Expression of p53, c-Myc, or Bcl-6 suggests a poor prognosis in primary central nervous system diffuse large B-cell lymphoma among immunocompetent individuals. Arch Pathol Lab Med 127:208-212, 2003
Aref S, Fouda M, El-Dosoky E, et al: C-Myc oncogene and Cdc25A cell activating phosphatase expression in non-Hodgkin's lymphoma. Hematology 8:183-190, 2003
Haslinger C, Schweifer N, Stilgenbauer S, et al: Microarray gene expression profiling of B-cell chronic lymphocytic leukemia subgroups defined by genomic aberrations and VH mutation status. J Clin Oncol 22:3937-3949, 2004(Dirk L. Kienle, Christian)
Department of Molecular Genetics and Department of Biostatistics, German Cancer Research Center, Heidelberg, Germany
Department of Internal Medicine I, Medical University of Vienna, Vienna, Austria
ABSTRACT
PURPOSE: In patients with chronic lymphocytic leukemia (CLL), the VH mutation status and genomic aberrations (13q每, +12q, 11q每, 17p每) identify distinct prognostic subgroups. The aim was to elucidate biologic mechanisms through which these genetic markers may exert their pathogenic influence.
PATIENTS AND METHODS: Twenty-four genes involved in apoptosis, cell cycle, B-cell activation, and B-cell receptor (BCR) signaling were analyzed by real-time quantitative reverse transcription polymerase chain reaction (RQ-PCR) in 82 CLL cases constituting prototypic genetic CLL subgroups as defined by the VH mutation status and the genomic aberrations 13q每, +12, 11q每, and 17p每.
RESULTS: The VH mutation subgroups were characterized by a differential expression of the BCR associated genes ZAP70 and PI3K. Among the subgroups defined by genomic aberrations, there was a deregulation of candidate genes from the affected critical genomic regions such as CDK4 (up), ATM (down), and TP53 (down) in the groups +12, 11q每, and 17p每, respectively. Additionally, the genomic subgroups were characterized by a significant deregulation of cell cycle and apoptosis regulators: AKT (up) in 13q, E2F1 (up) in +12, MYC (up) and BCL-2 (down) in 17p每, and CCND3 (down) in 11q每 as well as 17p每. The 17p每 subgroup showed an additional down-regulation of BCR-associated genes such as SYK and PI3K.
CONCLUSION: The characteristic gene expression patterns observed implicate a differential regulation of cell cycle, apoptosis, and BCR signaling in the genetic subgroups illustrating distinct pathomechanisms and are evidence for a gene dosage effect being operative in CLL. These findings link the biologic diversity and clinical heterogeneity of CLL.
INTRODUCTION
Chronic lymphocytic leukemia (CLL) is characterized by a variable clinical course, with survival times ranging from months to decades.1,2 This clinical heterogeneity has been shown to be related to inherent genetic features of the CLL cells such as the VH mutation status and genomic aberrations. While CLL with mutated VH displays a good prognosis, VH unmutated CLL is characterized by a poor outcome,3,4 suggesting that the different maturation stages according to the VH mutation status may define distinct CLL pathomechanisms.5 However, by genome-wide gene expression analysis using microarrays, the transcriptional profile of the VH subgroups could not be distinguished in unsupervised analyses, and only a small number of genes were found to be differentially expressed when comparing the VH subgroups by supervised analyses.6,7 Among these, ZAP70 had a prominent role, and further studies showed that B-cell receptor (BCR) crosslinking on ZAP70-positive CLL cells led to increased tyrosine phosphorylation (eg, of p72 [Syk]), indicating an increased activation after BCR stimulation as a pathomechanism.8
Additionally, genomic aberrations, which have an independent prognostic impact and may be transforming events in CLL, were not included in these studies. The most common genomic aberrations are losses at the chromosomal bands 13q14 (13q每: 55%), 11q22 to q23 (11q每: 18%), 17p13 (17p每: 7%), and gains at chromosome 12 (+12: 16%).9 17p每 and 11q每 identify poor prognostic subgroups in multivariate analysis.9-12 Recurrent genomic aberrations point to the loci of pathogenic genes such as ATM at 11q22-23 and TP53 at 17p13. However, the expression pattern of these and other candidate genes is so far unresolved in the subgroups of CLL defined by genomic aberrations.
The aim of the present study was to elucidate pathogenic events in pathways of oncogenic potential in the genetic subgroups of CLL. To investigate candidate gene expression of genetically distinct CLL subtypes, we built prototypic subgroups composed of cases defined by the VH mutation status and genomic aberrations, ie, VH mutated (VH MUT), VH unmutated (VH UM), normal karyotype by FISH (normal), 13q每 as sole abnormality (13q每), 11q每, +12, and 17p每. To obtain a highly sensitive quantitation of gene expression, we used real-time quantitative reverse transcription polymerase chain reaction (RQ-PCR) and focused on a selected set of functionally relevant genes playing central roles in the regulation of apoptosis, cell cycle, B-cell activation, and BCR signaling, namely AICL (CLECSF2), AID (AICDA), AKT1, ATM, BAX, BCL2, CCND1, CCND2, CCND3, CDK4, MYC, E2F1, P21 (CDKN1A), P27 (CDKN1B), PI3K (PIK3CB), MDM2, MCL1, NFKB1, Survivin (BIRC5), SYK, TCL (TCL1A), TP53, TRADD, and ZAP70.
PATIENTS AND METHODS
Patients
Peripheral blood was obtained for diagnostic procedures from all cases after informed consent and the studies were institutional review board每approved. Samples from 82 patients with CLL diagnosed according to established morphologic and immunophenotypic criteria were included. Staging information was complete for 76 patients; 30 patients were Binet A; 26, Binet B; and 20, Binet C. Median age was 57 years (range, 38 to 84 years), and 56% were male. All samples were from the same date as the corresponding VH and fluorescent in situ hybridization (FISH) analyses. Sixty-seven samples were collected at time of diagnosis (15 during the course of the disease), and all except seven patients were untreated. In these cases, mostly patients with refractory disease, one to three previous therapies had been administered, and the samples were collected in active recurrent disease after at least 2 treatment-free months. Clinical follow-up for survival analysis was available for 60 patients. Translocation t(11;14)每positive cases were not included.
Genetic Analysis
FISH analysis and VH sequencing was performed in all cases as previously described.9,12,13 A germline homology of 98% was used as the cutoff between VH mutated and VH unmutated cases.
Genetic Subgroups
To avoid the analysis of subclones, a high tumor load as defined by a high percentage of cells carrying the respective genomic aberration detected by FISH was the prerequisite for inclusion in the study. Median values of aberrant cells were 84% in 13q每 (range, 75% to 93%), 80% in +12 (range, 71% to 85%), 90% in 11q每 (range, 77% to 97%), and 89% in 17p每 (range, 73% to 97%). In cases with a normal karyotype by FISH, leukocyte and lymphocyte counts, in combination with CD5/CD19 expression in fluorescence-assisted cell sorter analysis, were used to determine the tumor load (median, 93%; range, 72% to 98%). The following subgroups were studied according to their genetic characteristics as defined by VH mutation status and genomic aberrations (Fig 1): VH MUT and VH UM (each equally divided into 13q- and normal); 13q每 (as single abnormality) and normal (no aberration by FISH; each equally divided into VH MUT and VH UM); +12 (and absence of 11q每 or 17p每), 11q每 (and absence of +12 or 17p每), and 17p每 (and absence of +12 or 11q-), that is, there was no occurrence of a +12, 11q每, or 17p每 outside the subgroup defined by these aberrations. If existing, additional aberrations in the groups 11q每, +12, and 17p每 consisted mainly of 13q每. In rare cases, other abnormalities (6q每, four cases; 14q每, three cases; +3q, one case; and +8q, one case) were present. The subgroups 11q每, +12, and 17p每 consisted mainly of VH UM cases: 11 of 15 with a +12, 16 of 16 with a 11q每, and 14 of 17 with a 17p每, allowing the additional comparison of homogeneously VH UM cases with +12 (+12[UM]: n = 11), 11q- (11q每[UM]: n = 16), or 17p每 (17p每[UM]: n = 14) versus VH UM cases without these aberrations (VH UM). Median VH homology in the VH MUT and VH UM subgroups was 92.8% (range, 88.4% to 95.7%) and 99.7% (range, 98.1% to 100%), respectively. For standardization purposes, CD19-positive fractions purified by magnet-activated cell separation CD19 MicroBeads (Miltenyi Biotech, Bergisch Gladbach, Germany) from buffy coat preparations of four healthy individuals were used.
Sample and RNA Preparation
The procedures were done as previously reported,14 including a DNAse I digestion of total RNA to avoid contamination with genomic DNA.
RQ-PCR
RQ-PCR analysis was carried out as previously described14 using the ABI Prism 7700 Sequence Detector (Applied Biosystems, Weiterstadt, Germany). Quantitative assessment of DNA amplification was achieved either through SYBR Green or TaqMan probes according to the manufacturer's manuals (SYBR Green PCR Core Reagents, TaqMan Universal Master Mix; Applied Biosystems). For SYBR Green reactions, a melting curve of the final RQ-PCR products was generated to avoid nonspecific fluorescence derived from side products such as primer-dimers (ABI Prism Dissociation Curve Software; PE Applied Biosystems, Foster City, CA). Primers for SYBR Green detection were designed using Primer Express software (Applied Biosystems) and are listed in Table 1. Primers for the following genes were used as previously described: BAX, CCND3, CDK4, E2F1, CMYC, P21, P27, TP53, lamin B1 (LMNB1), and phosphoglycerate kinase 1 (PGK1)14; CCND115; and Survivin.16 For beta actin (ACTB) PreDevelopedAssay-reagents, and for AID as well as ATM, Assay-on-Demand Gene Expression Probes were used (all supplied by Applied Biosystems).
Data Normalization
Eight serial dilutions of cDNA from tumor cell lines were used to obtain a calibration graph, which was measured in every single RQ-PCR run. In general, the B-cell CLL cell line EHEB was used for calibration. The T-cell line Jurkat and the Burkitt cell line Namalva were used for calibration of ZAP70 and AID expression levels, respectively, due to absent expression in EHEB. The average value of three independent endogenous control amplicons (PGK1, LMNB1, ACTB) were used for normalization.17
Control Experiments
Experiments were repeated for TP53, MYC, CDK4, BCL2, CCND3, E2F1, and for all three housekeeping genes to obtain an independent data set.
Statistical Analysis
For gene screening, we computed unadjusted and permutation-adjusted P values using standardized Wilcoxon rank sum test statistics. Permutation-adjusted P values were computed for the minP step-down multiple testing procedure described in Westfall and Young.18 This procedure provides strong control for the family-wise error rate and also incorporates the joint dependence structure between the test statistics. The number of permutations used was 10,000. In addition, we used summary measures based on receiver operating characteristic (ROC) curves as described by Pepe et al.19 The area under the ROC curve, which corresponds to the Wilcoxon rank sum test statistic, scores for discrimination ability (and equals 0.5 for a random classifier). Fold changes of expression values between compared subgroups were estimated by "back-transformation" of the Hodges-Lehmann estimate of the difference of the location parameters of the log-transformed expression data. Gene correlations were estimated by Spearman's statistic and plotted as a red and blue color image of the correlation matrix. Increasingly positive correlations are represented with reds of increasing intensity, and increasingly negative correlations are represented with blues of increasing intensity. To present correlation matrices we started with the pair of genes with the largest positive coefficient of pairwise correlation. The next gene selected is the one whose expression values show the strongest positive correlation with the mean expression of the two previously selected genes. This procedure is continued by selecting the gene whose expression values have the strongest positive correlation with the mean expression of all the genes previously selected.
RESULTS
Quantitative expression levels of 24 genes (listed in Table 2 with their chromosomal position) with central roles in the regulation of apoptosis, cell cycle, BCR signaling, and B-cell activation were analyzed by RQ-PCR in 82 cases constituting the prototypic genetic CLL subgroups (VH MUT, VH UM, 13q每, normal, +12, 11q每, and 17p每). A survival probability analysis was performed to test whether the selected cases were a representative collection with regard to the clinical impact of the genetic markers to be analyzed (Fig 2). The individual subgroups VH MUT, VH UM, +12, 11q每, and 17p每 showed the expected survival probabilities, demonstrating that the selected cases represent a valid collection for this study.
VH Mutation Subgroups
The VH MUT versus VH UM analysis was restricted to the set of normal and 13q每 cases to avoid interference of potentially distinct expression signatures of +12, 11q每, and 17p每 cases. Only two of the 24 genes analyzed showed significant differential distributions of expression levels between the VH mutation subgroups (summarized in Table 3): ZAP70 levels were 3.24-fold higher in VH UM than in VH MUT cases (Fig 3A), whereas PI3K showed lower expression levels in VH UM cases (0.68-fold). CD3 RNA expression levels, used as a calibrator for T-cell contamination, were comparable between the VH subgroups, and a correction for CD3 expression did not change the results for ZAP70 significantly (data not shown). The CCND1/CCND3 ratio was higher in VH MUT. When comparing VH mutated and VH unmutated cases in the overall cohort (n = 82; Table 4) including +12, 11q每, and 17p每 cases, similar results were obtained for ZAP70 and PI3K, but in addition, a down-regulation of ATM and CCND1 was evident, probably as a consequence of the inclusion of 11q每 and 17p每 cases, respectively (see Results section Genomic Aberrations).
Due to the potential pathogenic impact of ZAP70 expression in CLL, an additional comparison of ZAP70-high and ZAP70-low expressing cases was carried out. Therefore, 21 cases from the quartile with the highest ZAP70 expression were compared with 21 cases from the quartile with the lowest expression. In the ZAP70-high group, an upregulation of SYK (P = .02, 1.59-fold) and a down-regulation of PI3K (P = .04, 0.71-fold) were identified.
Genomic Aberrations
When comparing the different subgroups defined by genomic aberrations, a majority of the investigated genes (17 of 24) showed a differential expression pattern.
13q每 (single) was compared with normal cases (Table 5). The most characteristic finding was a low-level but highly differential overexpression of AKT (P = .001) in 13q每 cases (Fig 3B). To explore whether the differential expression of AKT was restricted to cases with a 13q每 single aberration (ie, no additional abnormality) or a general phenomenon of 13q deletion, we compared all cases with a 13q deletion (n = 42) against all cases without a 13q deletion (n = 40), irrespective of additional aberrations. The higher expression of AKT was also observed when comparing the 13q deletion (irrespective of additional aberrations) with the non-13q deletion subgroup (Fig 4).
The subgroups +12, 11q每, and 17p每 were individually compared with the cohort of normal and 13q- cases (n = 34, Fig 1). This comparison was chosen due to the occurrence of an additional 13q每 in approximately half of the cases with +12, 11q每, and 17p每. In addition, to investigate expression characteristics without interferences of VH mutation status, a restricted analysis was carried out including only VH UM cases. In this comparison, the homogeneously unmutated subgroups +12(UM), 11q每(UM), and 17p每(UM) were compared individually with the VH UM subgroup without these aberrations.
When comparing the +12 subgroup with the normal/13q每 cohort (Table 6), characteristic expression patterns were observed for E2F1 (P < .001), BAX (P = .002), P27 (P = .022), and CDK4 (P = .050). Notably, all of these differentially expressed genes showed an overexpression (1.33- to 1.73-fold) in +12. BCL2/BAX and CCND1/CCND3 ratios were significantly lower in +12. The restricted analysis of +12(UM) versus VH UM (Table 7) was statistically limited by the small number of +12(UM) cases (n = 11). However, comparable alterations of candidate gene expression levels were identified, confirming the upregulation of E2F1.
In the 11q每 subgroup (Table 8), the most differentially expressed genes were ATM (Fig 3C) and CCND3, both with a lower expression in 11q每 as compared with 13q每/normal (P < .001, 0.43-fold and P = .015, 0.66-fold, respectively). When the comparison was restricted to 11q每(UM) versus VH UM (Table 9), the lower ATM expression remained the most prominent finding (P < .001), but additional trends for upregulation of MCL and NFKB and down-regulation of SYK were noted.
The 17p每 subgroup displayed the highest number of differentially expressed genes affecting all investigated pathways (Tables 10 and 11). In comparison to normal/13q每, the 17p每 cases showed a significantly lower expression of CCND3 (P < .001), BCL2 (P = .001), SYK (P = .004), TP53 (P = .007; Fig 3D), ATM (P = .005), TCL (P = .006), PI3K (P = .021), CCND1 (P = .030), and AID (0.013) with a range of underexpression between 0.31- and 0.66-fold. Significant overexpression was observed for P21 (1.94-fold), MYC (2.18-fold), and AICL (3.57-fold). BCL2/BAX ratios were lower in the 17p每 subgroup. The comparison of 17p每(UM) versus VH UM produced similar results. In this analysis, TP53 was the most significantly down-regulated gene (P = .003).
Figure 5 gives a comparison of median expression levels of selected candidate genes involved in apoptosis, cell cycle regulation and BCR signaling in the genetic subgroups VH MUT, VH UM, +12, 11q每, and 17p每. ZAP70 expression levels were highest in VH UM, but lower in +12, 11q每, and 17p每 cases, though these were predominantly VH unmutated. Similarly, SYK expression was lower in +12, 11q每, and 17p每 as compared with VH UM. While ATM expression was lowest in 11q每 cases, TP53 levels were lowest in 17p每 cases. Both subgroups harbor a genomic loss of the respective gene locus. Surprisingly, BCL2 showed the highest expression in the VH MUT subgroup but the lowest expression in the poor-risk 17p每 subgroup. BAX expression was homogeneous except for an overexpression in +12. Whereas CCND3 expression was characteristically low in 11q每 and 17p每, MYC was overexpressed in 17p每.
Candidate genes that discriminated the genomic subgroups +12, 11q每, and 17p每 from normal/13q每 (P .05) were used for an outline of the major pathways in which the investigated genes are involved and the mode of their deregulation in the genetic CLL subgroups (Fig. 6A to C).
Gene Expression Inter-Relations
For a general comparison of subgroup-specific expression inter-relations, all genes were clustered hierarchically by their strength of correlation (Figs 7A to E), starting with the pair of genes with the largest positive coefficient of pairwise correlation (see Patients and Methods). Using this algorithm, highly distinct, subgroup-specific expression inter-relations can be illustrated. Whereas in VH MUT the strongest cluster included positive correlations between MCL, NFKB, and P27, the cluster consisted of CCND2, TCL, and BCL2 in VH UM (Figs 7A and B). Although the subgroups +12, 11q每, and 17p每 were predominantly composed of VH unmutated cases, their clustering was distinct for each of the subgroups and differed from the VH UM group. The strongest cluster in +12 consisted of PI3K, ATM, and TCL, in contrast to MCL, P21, and MDM2 in 11q每, and BAX and TRADD in 17p (Fig 7C每E).
Control Experiments
To ensure the reproducibility of the results, experiments were repeated in replicate with independent measurements of six candidate genes and all housekeeping genes. Highly comparable results were obtained. Distributions with initial P values less than .05 could be reproduced in all cases, and no additional significant results were obtained (data not shown).
DISCUSSION
Although genetic factors have a strong impact on the clinical course of CLL, little is known about the pathomechanisms of the genetic subgroups. Despite their highly distinct clinical behavior it was not possible to separate the CLL subgroups with mutated and unmutated VH genes by unsupervised cluster analysis of global gene expression data derived from microarrays.6,7 However, in these studies, genomic aberrations were not included, and the analyses were not performed by a highly quantitative method such as RQ-PCR, possibly obscuring subtle expression differences in a wealth of data points. In contrast to previous expression analyses, we used unmanipulated cell samples (ie, no CD19+ selection), but selected for a high content of genetically defined tumor cells. We thereby gained a "genetic" tumor purity of approximately 85% as compared with greater than 95% "B-cell purity" by CD19 selection. A disadvantage of this approach was a potentially remaining influence of non-CLL cells in the sample. This applies particularly for ZAP70, which is strongly expressed in T-cells. Since the current study focused on subgroup comparisons, this problem may have been reduced by the balanced content of genetically defined CLL cells and contaminating T-cells (monitored by CD3) in the different subgroups.
In line with pivotal studies6,7 only few differences in gene expression were observed between the VH mutation subgroups. ZAP70 was previously shown to be related to increased BCR signaling in CLL.8 In line with this, we identified elevated levels of SYK, a ZAP70-related kinase, in cases with high ZAP70 expression, arguing for a synergistic role of ZAP70 and SYK in initial BCR signaling. In contrast, PI3K expression, a gene known to be involved in BCR signaling by linking initial receptor signals from tyrosine kinases to downstream effectors,20-23 was lower in VH UM CLL. Similarly, PI3K was identified between the ZAP70 subgroups with a lower expression in ZAP70 high expressing cases. In contrast to Chen et al,8 these findings suggest an inactivity or suppression of physiologic downstream BCR cascades as a consequence of ZAP70 overexpression, potentially favoring atypical cellular responses or alternative pathways following BCR stimulation.
In CLL characterized by genomic aberrations, multiple gene deregulations affecting different pathways were observed, suggesting a combined model of cell-cycle deregulation, impairment of apoptosis, and dysfunctional cellular responsiveness for these subgroups. Importantly, most deregulations between genomic subgroups were reidentified when comparing the genomic subgroups in purely VH unmutated cases, clearly demonstrating the additional impact of genomic aberrations on the gene expression pattern independently of the VH mutation status.
Although 13q deletion is the most common genomic abnormality in CLL, its biologic role is still largely unresolved.24,25 In 13q每, the most characteristic finding was a low-level overexpression of AKT, which was also observed in the overall cohort of cases with 13q deletion, despite the inclusion of heterogeneous cases with secondary genetic events (11q每, +12, 17p每). Therefore, AKT overexpression seems to be a general phenomenon of 13q deletion potentially providing a new pathomechanism for this genetic abnormality, and not a phenomenon linked only to the favorable 13q每 (single) group. AKT can be activated by the BCR and promotes B-cell survival by inactivating components of the apoptotic machinery.26,27 In CLL cells, CD40-crosslinking induced cell proliferation, which was associated with an increase of AKT.28 Additionally, AKT activation was shown to protect CLL cells from chlorambucil- and radiation-induced apoptosis.29 This might explain a survival advantage for tumor clones with enhanced AKT expression, thereby favoring cells with a 13q deletion. Whether AKT overexpression could result from a deregulation of candidate genes from the critical genomic region on 13q14, has to be assessed.
Remarkably, in +12, all of the differentially expressed genes (E2F1, BAX, P27, CDK4) were overexpressed. Overexpression of CDK4 and P27 might be explained by a direct gene dosage-dependant upregulation due to their location on chromosome 12. The transcription factor E2F1 is a critical determinant of the G1/S-phase transition during cell cycle, leading to the transcription of S-phase proteins.30-32 Overexpression of E2F1 could be correlated with proliferation in human cancers and B-cell lymphomas.33,34 Therefore, E2F1 overexpression in +12 suggests enhanced proliferative activity as a potential pathomechanism in the evolution of this genetic abnormality. In addition, E2F1 activity is known to be regulated through the action of G1 cyclins and associated kinases such as CDK4.30-32 Therefore, overexpression of CDK4 might pathomechanistically contribute to E2F1 expression in +12.
Deletion in chromosome bands 11q22-23 is associated with extensive lymph node involvement and poor survival in CLL.9 ATM is located within the minimally deleted region and coordinates responses to DNA damage and is able to induce repair mechanisms and apoptosis by activating the TP53 pathway,35 suggesting that ATM might be acting as a tumor suppressor gene. However, only about one third of the CLL cases with 11q每 exhibit ATM mutations of the second allele,36,37 implicating the requirement of additional ATM inactivation mechanisms in CLL with 11q每. In this study, we demonstrate that the majority of 11q每 cases (13 of 16) exhibit reduced ATM transcript levels, pointing to a gene dosage-dependant down-regulation. Pettitt et al showed that p53 dysfunction in CLL may occur independently of TP53 mutations.38 Interestingly, these cases were associated with ATM mutations and exhibited reduced ATM protein expression. Therefore, down-regulation of ATM transcription in 11q- CLL may be a pathogenic factor contributing to dysfunctional apoptosis and dismal clinical outcome.
The largest number of deregulated genes was detected in the 17p- subgroup possibly reflecting the distinctly aggressive biology of this subgroup. This finding is in contrast to Stankovic et al39 reporting an indistinguishable native gene expression pattern when comparing wild-type and TP53-mutant CLL. The limited number of cases studied by Stankovic et al and the higher sensitivity of RQ-PCR in comparison to the microarray approach may account for this discrepancy. The characteristic two-fold reduction of TP53 transcript observed in our series indicates a gene-dosage effect as additional or synergistic mechanism contributing to TP53 dysfunction in CLL. Of note, we detected significantly reduced ATM levels in cases with 17p每, a finding that might contribute synergistically to apoptosis impairment in the 17p每 subgroup. The consensual down-regulation of ATM in the 11q每 and 17p每 subgroup argues for related pathogenic pathways involving a common mechanism of ATM inactivation as previously suggested.38,39 Overexpression of P21 despite low TP53 expression in 17p每 appears unexpected and may indicate TP53-independent counter-regulation of P21 in 17p每 CLL.
The finding of particularly low BCL2 levels in 17p每, the subgroup with the worst clinical course, was surprising due to the known antiapoptotic function of BCL2 and the association of high BCL2/BAX ratios with aggressive disease.40-42 This finding argues against a relevant role of BCL2 overexpression or elevated BCL2/BAX ratios in the pathogenesis of the genetic CLL subgroups or in CLL prognosis.
In addition to the down-regulation of apoptosis-associated genes in 17p每, there was a profound deregulation of cell-cycle regulators, including CCND3, CCND1, and MYC. CCND1 and, more pronounced, CCND3, showed reduced transcript levels in 17p每. However, as a consequence of this consensual down-regulation, the relationship between the transcript levels of both genes (CCND1/CCND3 ratio) remained largely unchanged. In contrast to the cyclins, MYC showed a differential overexpression in 17p每. MYC is known as a potent activator of cell proliferation,43 and MYC expression has been suggested as a marker for poor prognosis in B-cell lymphomas,44-46 implicating enhanced proliferative activity as an additional pathomechanism of the 17p每 subgroup.
A common phenomenon observed in the subgroups +12, 11q每, and 17p每 was the deregulation of genes from critical genomic regions implicating a gene-dosage effect being operative in CLL, which is in line with a recent report.47 In addition to the differential expression of individual genes described in this work, cluster analysis of gene expression correlations lead to subgroup-specific patterns reinforcing the concept of pathomechanisms distinct for the genetic CLL subgroups and thereby linking the biologic diversity and clinical heterogeneity of CLL.
Authors' Disclosures of Potential Conflicts of Interest
The authors indicated no potential conflicts of interest.
NOTES
Supported by the DFG (Sti 296/1-1), Deutsche Krebshilfe (70-3183-Li1), Wilhelm Sander Stiftung (2002.095.1), Jubilumsfonds of the Austrian National Bank (No. 9964; U.J.).
Presented at the Annual Meeting of the American Society of Hematology, San Diego, CA, 2003.
Authors' disclosures of potential conflicts of interest are found at the end of this article.
REFERENCES
Rozman C, Montserrat E: Chronic lymphocytic leukemia. N Engl J Med 333:1052-1057, 1995
Zwiebel JA, Cheson BD: Chronic lymphocytic leukemia: Staging and prognostic factors. Semin Oncol 25:42-59, 1998
Hamblin TJ, Davis Z, Gardiner A, et al: Unmutated Ig V(H) genes are associated with a more aggressive form of chronic lymphocytic leukemia. Blood 94:1848-1854, 1999
Damle RN, Wasil T, Fais F, et al: Ig V gene mutation status and CD38 expression as novel prognostic indicators in chronic lymphocytic leukemia. Blood 94:1840-1847, 1999
K邦ppers R, Klein U, Hansmann ML, et al: Cellular origin of human B-cell lymphomas. N Engl J Med 341:1520-1529, 1999
Rosenwald A, Alizadeh AA, Widhopf G, et al: Relation of gene expression phenotype to immunoglobulin mutation genotype in B cell chronic lymphocytic leukemia. J Exp Med 194:1639-1647, 2001
Klein U, Tu Y, Stolovitzky GA, et al: Gene expression profiling of B cell chronic lymphocytic leukemia reveals a homogeneous phenotype related to memory B cells. J Exp Med 194:1625-1638, 2001
Chen L, Widhopf G, Huynh L, et al: Expression of ZAP-70 is associated with increased B-cell receptor signaling in chronic lymphocytic leukemia. Blood 100:4609-4614, 2002
Dhner H, Stilgenbauer S, Benner A, et al: Genomic aberrations and survival in chronic lymphocytic leukemia. N Engl J Med 343:1910-1918, 2000
Geisler CH, Philip P, Egelund, et al: In B-cell chronic lymphocytic leukaemia chromosome 17 abnormalities and not trisomy 12 are the single most important cytogenetic abnormalities for the prognosis: A cytogenetic and immunophenotypic study of 480 unselected newly diagnosed patients. Leuk Res 21:1011-1023, 1997
Oscier DG, Gardiner AC, Mould SJ, Glide S, Davis ZA, Ibbotson RE, Corcoran MM, Chapman RM, Thomas PW, Copplestone JA, Orchard JA, Hamblin TJ: Multivariate analysis of prognostic factors in CLL: Clinical stage, IGVH gene mutational status, and loss or mutation of the p53 gene are independent prognostic factors. Blood 100:1177-1784, 2002
Krber A, Seiler T, Benner A, et al: V(H) mutation status, CD38 expression level, genomic aberrations, and survival in chronic lymphocytic leukemia. Blood 100:1410-1416, 2002
Kienle D, Krber A, Katzenberger T, et al: VH mutation status and VDJ rearrangement structure in mantle cell lymphoma: Correlation with genomic aberrations, clinical characteristics, and outcome. Blood 102:3003-3009, 2003
Korz C, Pscherer A, Benner A, et al: Evidence for distinct pathomechanisms in B-cell chronic lymphocytic leukemia and mantle cell lymphoma by quantitative expression analysis of cell cycle and apoptosis-associated genes. Blood 99:4554-4561, 2002
Bijwaard KE, Aguilera NS, Monczak Y, et al: Quantitative real-time reverse transcription-PCR assay for cyclin D1 expression: Utility in the diagnosis of mantle cell lymphoma. Clin Chem 47:195-201, 2001
Schmidt SM, Schag K, M邦ller MR, et al: Survivin is a shared tumor-associated antigen expressed in a broad variety of malignancies and recognized by specific cytotoxic T cells. Blood 102:571-576, 2003
Mertens D, Wolf S, Schroeter P, et al: Down-regulation of candidate tumor suppressor genes within chromosome band 13q14.3 is independent of the DNA methylation pattern in B-cell chronic lymphocytic leukemia. Blood 99:4116-4121, 2002
Westfall PH, Young SS: Resampling-Based Multiple Testing: Examples and Methods for P-Value Adjustment. New York, NY, John Wiley & Sons, 1993, p 66
Pepe MS, Longton G, Anderson G, et al: Selecting differentially expressed genes from microarray experiments. Biometrics 59:133-142, 2003
Beitz LO, Fruman DA, Kurosaki T, et al: SYK is upstream of phosphoinositide 3-kinase in B cell receptor signaling. J Biol Chem 274:32662-32666, 1999
Bone H, Williams NA: Antigen-receptor cross-linking and lipopolysaccharide trigger distinct phosphoinositide 3-kinase-dependent pathways to NF-kappa B activation in primary B cells. Int Immunol 13:807-816, 2001
Bernal A, Pastore RD, Asgary Z, et al: Survival of leukemic B cells promoted by engagement of the antigen receptor. Blood 98:3050-3057, 2001
Glassford J, Soeiro I, Skarell SM, et al: BCR targets cyclin D2 via Btk and the p85alpha subunit of PI3-K to induce cell cycle progression in primary mouse B cells. Oncogene 22:2248-2259, 2003
Calin GA, Dumitru CD, Shimizu M, et al: Frequent deletions and down-regulation of micro- RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A 99:15524-15529, 2002
Migliazza A, Bosch F, Komatsu H, et al: Nucleotide sequence, transcription map, and mutation analysis of the 13q14 chromosomal region deleted in B-cell chronic lymphocytic leukemia. Blood 97:2098-2104, 2001
Datta SR, Brunet A, Greenberg ME: Cellular survival: A play in three Akts. Genes Dev 13:2905-2927, 1999
Pogue SL, Kurosaki T, Bolen J, et al: B cell antigen receptor-induced activation of Akt promotes B cell survival and is dependent on Syk kinase. J Immunol 165:1300-1306, 2000
Cuni S, Perez-Aciego P, Perez-Chacon G, et al: A sustained activation of PI3K/NF-kappaB pathway is critical for the survival of chronic lymphocytic leukemia B cells. Leukemia 18:1391-1400, 2004
Jones DT, Ganeshaguru K, Anderson RJ, et al: Albumin activates the AKT signaling pathway and protects B-chronic lymphocytic leukemia cells from chlorambucil- and radiation-induced apoptosis. Blood 101:3174-3180, 2003
Bartek J, Lukas J: Pathways governing G1/S transition and their response to DNA damage. FEBS Lett 490:117-122, 2001
Sherr CJ, Roberts JM. CDK inhibitors: Positive and negative regulators of G1-phase progression. Genes Dev 13:1501-1512, 1999
Van der Sman J, Thomas NS, Lam EW: Modulation of E2F complexes during G0 to S phase transition in human primary B-lymphocytes. J Biol Chem 274:12009-12016, 1999
Zacharatos P, Kotsinas A, Evangelou K, et al: Distinct expression patterns of the transcription factor E2F-1 in relation to tumour growth parameters in common human carcinomas. J Pathol 203:744-753, 2004
Chiarle R, Budel LM, Skolnik J, et al: Increased proteasome degradation of cyclin-dependent kinase inhibitor p27 is associated with a decreased overall survival in mantle cell lymphoma. Blood 95:619-626, 2000
Waterman MJ, Stavridi ES, Waterman JL, et al: ATM-dependent activation of p53 involves dephosphorylation and association with 14-3-3 proteins. Nat Genet 19:175-178, 1998
Schaffner C, Stilgenbauer S, Rappold GA, et al: Somatic ATM mutations indicate a pathogenic role of ATM in B-cell chronic lymphocytic leukemia. Blood 94:748-753, 1999
Austen B, Fegan C, Moss P, et al ATM mutations define a subset of B-CLL tumors that only partially correlates with 11q status and is associated with impaired patient survival. Blood 104, 2004 (abstr 774)
Pettitt AR, Sherrington PD, Stewart G, et al: p53 dysfunction in B-cell chronic lymphocytic leukemia: Inactivation of ATM as an alternative to TP53 mutation. Blood 98:814-822, 2001
Stankovic T, Hubank M, Cronin D, et al: Microarray analysis reveals that TP53- and ATM-mutant B-CLLs share a defect in activating proapoptotic responses after DNA damage but are distinguished by major differences in activating prosurvival responses. Blood 103:291-300, 2004
Robertson LE, Plunkett W, McConnell K, et al: Bcl-2 expression in chronic lymphocytic leukemia and its correlation with the induction of apoptosis and clinical outcome. Leukemia 10:456-459, 1996
Aguilar-Santelises M, Rottenberg ME, Lewin N, et al: Bcl-2, Bax and p53 expression in B-CLL in relation to in vitro survival and clinical progression. Int J Cancer 69:114-119, 1996
Molica S, Dattilo A, Giulino C, Levato D, Levato L: Increased bcl-2/bax ratio in B-cell chronic lymphocytic leukemia is associated with a progressive pattern of disease. Haematologica 83:1122-1124, 1998
Bouchard C, Staller P, Eilers M: Control of cell proliferation by Myc. Trends Cell Biol 8:202-206, 1998
Nagy B, Lundan T, Larramendy ML, et al: Abnormal expression of apoptosis-related genes in haematological malignancies: Overexpression of MYC is poor prognostic sign in mantle cell lymphoma. Br J Haematol 120:434-441, 2003
Chang CC, Kampalath B, Schultz C, et al: Expression of p53, c-Myc, or Bcl-6 suggests a poor prognosis in primary central nervous system diffuse large B-cell lymphoma among immunocompetent individuals. Arch Pathol Lab Med 127:208-212, 2003
Aref S, Fouda M, El-Dosoky E, et al: C-Myc oncogene and Cdc25A cell activating phosphatase expression in non-Hodgkin's lymphoma. Hematology 8:183-190, 2003
Haslinger C, Schweifer N, Stilgenbauer S, et al: Microarray gene expression profiling of B-cell chronic lymphocytic leukemia subgroups defined by genomic aberrations and VH mutation status. J Clin Oncol 22:3937-3949, 2004(Dirk L. Kienle, Christian)