Original Article |
1 Servicios de Hematología, Hospital Infanta Leonor, Madrid;
2 IBMCC, Centro de Investigación del Cáncer, Universidad de Salamanca-CSIC;
3 Hospital Clínico Universitario de Salamanca;
4 Hospital Virgen Blanca, León;
5 Hospital del Río Hortega, Valladolid;
6 Hospital Virgen del Puerto, Plasencia, Cáceres;
7 Hospital Clínico Universitario, Valladolid;
8 Hospital General de Segovia;
9 Hospital del Bierzo, León;
10 Hospital Nuestra Señora de Sonsoles, Ávila and
11 Hospital Morales Meseguer, Murcia, Instituto de Estudios de Ciencias de la Salud de Castilla y León (IECSCYL), Unidad de Investigación, Hospital Universitario de Salamanca, Spain
Correspondence: Jesús María Hernández, Servicio de Hematología y, Departamento de Medicina, Hospital Universitario de, Salamanca Paseo San Vicente, 58, 37007 Salamanca, Spain., E-mail:jmhr{at}usal.es
|
|
|---|
Design and Methods: Three hundred and fifty patients with B-cell chronic lymphoid leukemia were studied. Clinical data were collected and fluorescence in situ hybridization and molecular studies were carried out. In addition, a gene expression profile was obtained by microarray-based analysis.
Results: In 109 out of the 350 cases (31.1%) loss of 13q was the sole cytogenetic aberration at diagnosis. In the subgroup of patients with 80% or more of cells with loss of 13q (18 cases), the overall survival was 56 months compared with not reached in the 91 cases in whom less than 80% of cells had loss of 13q (p< 0.0001). The variables included in the multivariate analysis for overall survival were the percentage of losses of 13q14 (p=0.001) and B symptoms (p=0.007). The time to first therapy in the group with 80% or more vs. less than 80% of losses was 38 months vs. 87 months, respectively (p=0.05). In the multivariate analysis the variables selected were unmutated status of IgVH (p=0.001) and a high level of β2microglobulin (p=0.003). Interestingly, these differences regarding overall survival and time to first therapy were also present when other cut-offs were considered. The gene expression profile of patients with a high number of losses in 13q14 showed a high proliferation rate, downregulation of apoptosis-related genes, and dysregulation of genes related to mitochondrial functions.
Conclusions: Patients with B-cell chronic lymphoid leukemia with a high number of losses in 13q14 as the sole cytogenetic aberration at diagnosis display different clinical and biological features: short overall survival and time to first therapy as well as more proliferation and less apoptosis. A quantification of the number of cells showing a genetic abnormality should, therefore, be included in the study of the prognostic factors of B-cell chronic lymphoid leukemia.
Key words: B chronic lymphoid leukemia, 13q14 deletion, outcome, proliferation, apoptosis.
|
|
|---|
B-CLL patients show several cytogenetic aberrations, mainly in the regions of chromosomes 13q, 12, 11q, 17p, 14q and 6q. Some of these abnormalities can be better assessed by means of fluorescence in situ hybridization (FISH), which has shown that 62–80% of patients with B-CLL have cytogenetic abnormalities.10,12 These cytogenetic changes are strongly correlated with the prognosis in terms of overall survival and time to progression (defined as the time to first therapy).12,16–19 Patients with a deletion in 13q14 have a better outcome while patients with deletions in 11q23 or 17p13 have a shorter survival and shorter time to progression.12 Classically, patients with B-CLL and a normal karyotype or trisomy 12 have been considered to have an intermediate prognosis.12 It should, however, be noted that, in some series with a long-term follow-up, patients with B-CLL and a normal karyotype showed a better survival from 12 years on, as compared to patients with 13q-.9 In addition, several studies have demonstrated that the percentage of cells displaying a particular cytogenetic abnormality (e.g. loss of P53)20 or antigenic markers (e.g. CD38 or ZAP-70)7 can be related to prognosis.
For these reasons, we decided to perform an analysis of patients diagnosed with B-CLL and deletion in 13q14, as the sole cytogenetic abnormality. The clinical features, including outcome, and the biological features of the patients displaying different degrees of infiltration by 13q- cells were assessed. Moreover, to gain further insights into the molecular mechanisms involved in 13q14 deletion B-CLL, a gene expression profile study using a microarray-based analysis was also carried out.
|
|
|---|
and
23. In addition, FISH studies, including specific probes for the regions 11q21, 12q13, 13q14, 14q32, and 17p13, were carried out. The study protocol was approved by local ethical committees and written informed consent was obtained from the patients.
Clinical data
Clinical data were recorded by reviewing the clinical histories of patients included in the study. In most cases (283 patients; 81%) the FISH study was performed at the time of diagnosis. In more than 95% (61 patients) of the remaining cases, (patients with a long follow-up), the FISH study was normal (25 patients) or showed alterations in 13q14 as the sole cytogenetic aberration (36 patients). Only six patients showed other cytogenetic alterations: two patients had 11q deletions, one had t(14q32) and three patients had two cytogenetic alterations (13q14 deletion plus t(14q32) in two cases and 13q14 deletion plus 11q deletion in the other one). Progression was defined according to previously reported criteria:24 the presence of disease-related symptoms, massive or progressive organomegaly, bone marrow failure or recurrent infections.
Fluorescence in situ hybridization
Interphase FISH was performed on bone marrow samples using commercially available probes for the following regions: 13q14, 12q13, 11q22/ATM, 17p13/P53, and 14q32/IGH (Vysis/Abbott Co, Downers Grove IL, USA). The methods used for the FISH analysis have been described elsewhere.25 14q32 translocations, trisomy 12 and deletions were considered to be present when
5%,
3% and
8% interphase cells showed a split signal, three signals and one signal, respectively. Dual-color FISH using differently-labeled control probes and test probes was performed and signal screening was carried out on at least 200 cells with well-delineated signals. Hybridization was repeated on those slides with less than 80% cells showing two control-probe signals.
Mutation status of IGVH genes
Amplification and sequencing of IGVH genes was performed according to the ERIC recommendations on IGHV gene mutational status analysis in B-CLL.26 Cases were classified as IGVH unmutated if there was at least 98% concordance between the tumor DNA and the respective family sequence, and IGVH mutated if there was less than 98% concordance.
Statistical analysis
Statistical tests were performed with SPSS 13.0 (SPSS, Chicago, IL, USA). The
2 test was used to assess associations between categorical variables, while continuous variables were analyzed with the Kruskal-Wallis test. The variables with statistical significance related to overall survival and time to first therapy were calculated by the Kaplan-Meier method (log-rank). Results were considered statistically significant for p values
0.05. Multivariate analysis of survival and time to first therapy was performed using the Cox regression method.
Gene expression profile analysis
Patients and samples
Bone marrow samples were obtained from 37 patients with B-CLL and deletion of 13q14 as the sole cytogenetic aberration at diagnosis. Fifteen had more than 80% of 13q- cells, while the remaining 22 cases had less than 80% of 13q- cells in the bone marrow. Mononuclear cells from all samples were isolated using Ficoll gradient, snap frozen and stored at –80ºC. Both groups of patients showed more than 80% of clonal B-cell lymphocytes. RNA isolation, labeling and microarray hybridization were performed as previously reported.27
Normalization, signal calculation, significant differential expression, and sample/gene profile clustering
A robust microarray analysis algorithm was used for background correction, intra- and inter-microarray normalization, and expression signal calculation.28–30 Once the absolute expression signal for each gene (i.e., the signal value for each probe set) had been calculated for each microarray, a method, called significance analysis of microarray,31 was applied to calculate significant differential expression and find the gene probe sets that characterized the samples of each compared state. This method uses permutations to provide robust statistical inference of the most significant genes and provides p values adjusted to multiple testing using the false discovery rate (FDR).32 A FDR of less than 0.05 was used for all the differential expression calculations. Finally, the resulting lists of candidate genes associated to a high degree with 13q14 band deletion were tested using another algorithm, the so-called global test,33 which reveals the group of genes that has a global expression pattern most significantly related to the clinical feature studied. We applied all these methods using R and Bioconductor.
The function of the genes included in the expression signature of CLL with a high degree of 13q14 was assigned by applying the GeneCodis program,34 which finds concurrent annotations in GO and KEGG, and thereby derives several groups of genes with functional significance. The functional analysis to identify the most relevant biological mechanisms, pathways and functional categories in the data sets of genes selected by statistical analysis was generated through the use of Ingenuity Pathways Analysis (Ingenuity Systems, Mountain View, CA, USA).
|
|
|---|
|
View this table: [in a new window] [Download PPT slide] |
Table 1. Incidence of genomic aberrations assessed by FISH in the global series of 350 patients with B-CLL.
|
![]() View larger version (17K): [in a new window] [Download PPT slide] |
Figure 1. (A) Overall survival and (B) time to first therapy in the 350 patients with B-CLL. (C) Overall survival of patients with 80% or <80% cells showing the 13q14 deletion against that of patients in the other cytogenetic groups. The median survival and median time to first therapy is shown in months for each cytogenetic group.
|
Regarding the percentage of cells with 13q-, the majority of cases (83.5%) had this abnormality in less than 80% of the analyzed cells. No relation was found between monoallelic and biallelic losses (22% in patients with <80% of 13q losses vs. 18% in the
80% group) in the 13q- patients. No major differences regarding clinical, biological, immunophenotypic or mutational status features were found in the cases with low (<80%) vs. high infiltration (
80%) of 13q- cells, except for a high lymphocyte count (median of 14 vs. 19.7x109/L, respectively) (p=0.007) and a trend for an association with a diffuse pattern of bone marrow infiltration (17% vs. 40%; p=0.07) and splenomegaly (11% vs. 28%; p=0.07) in the group with high 13q- (Table 2). Of 109 cases considered, all 33 treated received, fludarabine-based therapies; in three cases allogeneic transplantation with reduced intensity conditioning was also performed.
|
View this table: [in a new window] [Download PPT slide] |
Table 2. Characteristics of 109 patients with 13q14 deletion as the sole cytogenetic aberration at diagnosis, divided according to the percentage of losses detected by FISH: < 80% (n=91) or 80% (n=18).
|
![]() View larger version (13K): [in a new window] [Download PPT slide] |
Figure 2. (A) Overall survival and (B) time to first therapy of patients with B-CLL and 13q14 deletion as the sole aberration and <80% or 80% FISH losses.
|
40%, 41–69%,
70%, data not shown, and
50%, 51–79% and
80%) (Online Supplementary Figure S1).
Gene expression profiles of the subsets of B-CLL patients displaying different degrees of 13q loss
For the gene expression profile analysis two groups of patients with 13q- were compared: those in whom 80% or more of cells showed 13q- (group A) and those in whom less than 80% of cells showed 13q losses (group B). The comparative analysis of the gene expression profile of both groups identified a set of 1755 differentially expressed genes: 1073 genes were upregulated in group A, and 682 were downregulated (Figure 3). The gene function analysis revealed that most of these genes are involved in apoptosis, cellular growth and proliferation, mitochondrial, endoplasmic reticulum and calcium mediated activity (Table 3). Thus, the patients with more than 80% of 13q- cells had overexpression of genes related to cell proliferation and MAP kinase activity, such as ATF2, FGF3, FGF10, FGFR2, RRAS2, GRB2, JUN, MAPT and MAP3K1. By contrast, this group of patients had downregulation of genes related to apoptosis (PDCD10, EGLN3, CASP6, CL, and DAPK1) and cell cycle arrest, such as CDKN2C, GAS2L1, UHMK1, GAS1, GAS2L3, ZAK, and GAS7. Most of the genes involved in proteasome function were downregulated in the group with 13q losses in more than 80% of cells. The proteasome (PSM) genes dysregulated in our study codify proteasomal subunits, both catalytic (PSMA 1–7) and regulatory (PSMD1, PSMD7) ones. The ubiquitin process was also altered. Thus, genes for several enzymes that participate in different steps of the ubiquitin cycle were downregulated in the group with more losses in 13q: peptidases (USP27X), ubiquitin-conjugating enzymes (UBE2), thiolesterases (UCHL3) and ligas-es (HECTD1) (Table 3).
![]() View larger version (66K): [in a new window] [Download PPT slide] |
Figure 3. Unsupervised analysis of patients displaying 13q- in more than 80% of the cells (group A) and the cases with less than 80% of 13q-. Overall, 1755 genes were deregulated, most of them (n=1073) were upregulated in group A.
|
|
View this table: [in a new window] [Download PPT slide] |
Table 3. Most relevant differentially expressed genes in patients with 80% or more cells showing 13q- (Up: overexpressed. Down: underexpressed).
|
|
|
|---|
Focusing on patients with loss of 13q as the sole cytogenetic aberration at diagnosis, we observed that the number of malignant cells carrying this genetic lesion influences the disease outcome. Thus, according to the percentage of cells with 13q- two prognostic groups could be established: the patients with a high proportion (
80%) of 13q- cells had both a shorter overall survival than that of patients with <80% 13q- cells (56 months vs. not reached, p<0.0001) and a shorter time to first therapy (38 months vs. 87 months, p=0.05). These differences persisted when other cut-offs were considered (Online Supplementary Figure S1). The clinical relevance of the percentage of cells displaying a specific genetic abnormality has been recently demonstrated in CLL patients. Thus, the presence of more than 20% of cells with loss of TP53 has been associated with an adverse prognosis, while patients with less than 20% of cells with loss of TP53 had a prognosis similar to that of the global series.20,36 In the present study, a high number of 13q- cells together with the presence of B-symptoms were the only independent adverse factors for a short survival in the multivariate analysis in the group of B-CLL patients with 13q-.
Our data support the concept that patients with 13q-do not constitute a homogeneous group. Thus, patients with a high proportion of 13q losses had a high lymphocyte count and a trend to have more frequently a diffuse pattern of bone marrow infiltration and splenomegaly. In addition, overall survival and time to first therapy in this group were shorter than in patients with a normal karyotype or trisomy 12, being similar to that in patients with 11q deletion, although better than in the rest of cytogenetic subtypes. By contrast, CLL patients with a low number of malignant cells carrying the 13q deletion had a better prognosis than patients with a normal karyotype. It could, therefore, be useful to include quantification of the number of cells showing a genetic abnormality as part of the study of the prognostic factors in this disease.
To test our hypothesis that differences in several cellular functions such as more proliferation and less apoptosis could justify the clinical and prognostic differences found in the subsets of patients with 13q- divided according to the percentage of cells with this loss, we carried out a gene expression profiling analysis. This methodology has been applied in B-CLL showing that patients with specific genomic aberrations, gene expression phenotype or IGVH mutation status have distinct gene expression profiles.37–39
In the present study, we identified important functions differentially deregulated in the two groups of patients displaying 13q-. The main differential pathways involved are related to apoptosis, proliferation, mitochondrial and endopolasmic reticulum function, as well as ubiquitin metabolism. Patients with more 13q- cells had overexpression of genes involved in proliferation. Thus, the MAPK signaling pathway was affected since GRBS, RRAS, JUN or SOS1 were upregulated in the group of B-CLL patients with a high number of losses in 13q. In addition, some genes related to cell cycle arrest, such as CDKN2C, GAS2L1, GAS1, GAS7, ZAK, and GAS2L3, were downregulated in the group of patients with more 13q. Both the upregulation of MAPK and the underexpression of genes related to cell cycle arrest would lead to greater cell proliferation.40 By contrast, genes related to apoptosis (CASP6, CLU, DAPK1, and E2F1) were found to be downregulated in the group of cases with a high proportion of 13q- cells, leading to a decrease in apoptotic activity. The accumulation of mature B cells that have escaped programmed cell death and undergone cell-cycle arrest in the G0/G1 phase is the hallmark of CLL.41 Moreover, in our study we found alterations in other apoptosis-related pathways such as those mediated by mitochondria, endoplasmic reticulum and calcium metabolism. Thus, endoplasmic reticulum and protein-vesicular transport are inhibited because of an upregulation of genes involved in calcium-mediated activities, such as calmodulin binding and calcium ion binding.42,43 Furthermore, mitochondrial oxi-doreductase activity was inhibited in the group of B-CLL patients with a high number of cells with 13q-. The ATP dependency of apoptosis is well-known and both the respiratory chain and ATP synthase itself have been attributed central roles in the apoptotic process.44,45 Therefore, patients with a high number of cells with 13q- not only had more proliferation, but also less apoptosis, which could be related to the more aggressive disease that we observed in the clinical study.
In recent years several hematologic malignancies have been treated with proteasome inhibitors.46 Our study supports the idea that proteasome function is inhibited in the group of B-CLL patients with a high proportion of 13q-.47 It is, therefore, logical to expect that such patients would not benefit from the use of proteasome inhibitors.
The PI3K/NF
B pathway seems to play a pivotal role in B-CLL cell survival and growth.48 In our study, we found deregulation in PI3K β and
catalytic subunits (PI3KCB and PI3KCG) as well as in other genes involved in this pathway (RAS, c-JUN, PKC and several integrins). This could be essential for mitogenic and antiapoptotic functions. In addition Ras-MAPK and Ras-PI3K signaling pathways are related and both of them lead cells to greater proliferation.
In summary, our findings suggest that the number of malignant cells with 13q- can influence the outcome of B-CLL patients in whom 13q- is the single genetic abnormality. A high number of cells with 13q deletion, as quantified by FISH, is associated with short overall survival and time to progression. To the best of our knowledge, this is the first time that the percentage of losses in 13q14 has been related to survival and progression. In addition, several cellular functions, such as more proliferation and less apoptosis, were found to be altered in the subsets of 13q- patients.
Our results need to be confirmed by additional studies, preferably in the context of large, randomized clinical trials.
J-AH designed and performed the research and statistical analyses, recorded and analyzed data and wrote the paper; A-ER: designed and performed the microarray experiments and wrote the paper; CF, JLG, N-CG and JdR performed and interpreted the gene expression profiling studies; MG and J-MH designed and performed research and wrote the paper; VR, MR, GM-N, AGC, RF, JG, IR and FO provided patients data and wrote the paper; J-FSM corrected and approved the final version of the manuscript.
The authors reported no potential conflicts of interest.
Received for publication August 27, 2008. Revision received October 17, 2008. Accepted for publication November 11, 2008.
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||