Chronic Lymphocytic Leukemia |
1 Unitat dHematopatologia i
2 Departament dHematologia, Hospital Clínic, Institut dInvestigacions Biomèdiques August Pi i Sunyer (IDIBAPS), Universitat de Barcelona, Barcelona
3 Departament de Bioquímica i Biologia Molecular, Institut de Biomedicina (IBUB), Universitat de Barcelona i CIBER EHD, Barcelona, Spain
Correspondence: Marcal Pastor-Anglada, Departament de Bioquímica i Biologia Molecular, Universitat de Barcelona, Diagonal 645, E-08028 Barcelona, Spain. E-mail:mpastor{at}ub.edu and Dolors Colomer, Unitat dHematopatologia, Hospital Clínic, Villarroel 170, 08036 Barcelona, Spain. E-mail:dcolomer{at}clinic.ub.es
|
|
|---|
Design and Methods: We performed gene expression profiling in cells from two fludarabine-sensitive and two fludarabine-resistant cases of chronic lymphocytic leukemia treated with fludarabine either in the presence or the absence of nitrobenzylthioinosine, a hENT1-specific blocker. Twenty selected fludarabine-inducible genes were validated using Taqman low-density arrays in cells from 20 chronic lymphocytic leukemia patients with the same experimental design.
Results: Sixteen of the twenty genes (DDB2, GADD45A, TYMS, BAX, TIGAR, FAS, TNFSF7, TNFSF9, CCNG1, CDKN1A, MDM2, SESN1, MAP4K4, PPM1D, OSBPL3 and WIG1) correlated with the ex vivo sensitivity of chronic lymphocytic leukemia cells to fludarabine, TIGAR (TP53-induced glycolysis and apoptosis regulator) being the gene that showed the strongest correlation (p<0.0001; r2= 0.6022).We observed that the transcriptomic response was weakly sensitive to the hENT1 blocker nitrobenzylthioinosine. Interestingly, we also found a correlation between hENT2 expression and induction of TIGAR after fludarabine treatment.
Conclusions: We demonstrate a correlation between the recently described p53-inducible apoptosis gene TIGAR and both sensitivity to fludarabine and hENT2 expression in chronic lymphocytic leukemia cells. These results, as well as the variability in fludarabine response among chronic lymphocytic leukemia patients with wild type p53, support the major role of hENT2 in the uptake of fludarabine into chronic lymphocytic leukemia cells.
Key words: TIGAR, hENT2, fludarabine, chronic lymphocytic leukemia.
|
|
|---|
Fludarabine is a prodrug that is converted to the free nucleoside 9-β-D-arabinosyl-2-fluoroadenine (F-ara-A), which enters cells and accumulates mainly as the 5-triphosphate, F-ara-ATP. F-ara-ATP has multiple mechanisms of action, which are mostly directed towards DNA.3 Several in vitro studies suggest that fludarabine triggers p53-mediated apoptosis in CLL cells4,5 although cell death by p53-independent mechanisms has also been described.6 Responses, both in vitro and in vivo, are very heterogeneous, there being reported cases with functional p53 but low response to fludarabine. The mechanisms implicated in these processes are not well known. Therefore, no suitable system to predict clinical outcome in CLL patients treated with fludarabine is available.
Anticancer therapy using nucleoside-derived analogs is dependent on drug transport across the plasma membrane and subsequent metabolic activation. Although some enzymes have been shown to be suitable biomarkers of nucleoside metabolism, thus modulating response to therapy,7,8 the role of plasma membrane transporters in determining nucleoside-derived drug bioavailability is less well known.
The uptake of nucleosides and nucleoside-derived drugs into cells is mediated by concentrative nucleoside transporter (CNT) and equilibrative nucleoside transporter (ENT) proteins, encoded by the gene families SLC28 and SLC29, respectively.9,10 CNT proteins mediate high-affinity, concentrative nucleoside transport into cells, whereas equilibrative transporters are responsible for the facilitative uptake of nucleosides and nucleoside-derived drugs, with broad selectivity but lower affinity than that shown by CNT proteins.11
CLL cells express hENT1, hENT2, hCNT2 and hCNT3 mRNA, although most of the biological activity responsible for fludarabine uptake is associated with hENT transporters.12 This is probably due to the fact that the hCNT3 protein is localized mostly in intracellular compartments of CLL cells,13 while fludarabine is not a suitable permeant for hCNT2,11 a transporter protein that is functional at the plasma membrane of CLL cells.12 In fact, although both hENT1 and hENT2 can mediate fludarabine internalization, hENT2 but not hENT1 protein levels have recently been correlated with ex vivo sensitivity of CLL cells to this nucleoside analog.14 The possibility that a single drug, despite having several transporter proteins as putative internalization pathways, channels its cytotoxic action via a particular route of entry has been recently addressed in the breast cancer cell line MCF7. In these cells, the transcriptomic response that follows treatment with 5-deoxy-5-fluorouridine (5-DFUR) is selectively blocked by the ENT1-specific inhibitor nitrobenzylthioinosine (NBTI),15 even though 5-DFUR is also a ENT2 substrate and MCF7 cells show significant ENT2-related activity.
Here we report a high-throughput analysis of gene expression in CLL cells treated ex vivo with fludarabine as a way to identify putative new biomarkers of fludarabine responsiveness in this disease, and we also provide evidence of transporter-mediated transcriptomic responses to nucleoside-derived drugs in primary CLL cells.
|
|
|---|
|
View this table: [in a new window] [Download PPT slide] |
Table 1. Characteristics of the patients with chronic lymphocytic leukemia.
|
Treatment and flow cytometry detection of apoptosis
Cells were incubated for 3, 24 and 48 hours with fludarabine 1 µg/mL (Schering, Berlin, Germany), a concentration that is comparable to that achieved in the blood during CLL treatment.18 When indicated, cells were preincubated for 1 hour with the hENT1 inhibitor nitrobenzylthioinosine (NBTI; Sigma, St Louis, MO, USA) at 100 nM before drug exposure. Cell viability was quantified by double staining with annexin-V conjugated to fluorescein isothiocyanate (FITC) and propidium iodide (PI) (BenderMedsystems, Vienna, Austria). Ten thousand cells per sample were acquired in a FACScan flow cytometer (Becton Dickinson, San Jose, CA, USA) and the labeled populations of cells were analyzed with the Paint-A-Gate software (Becton Dickinson).
High-density array study
Total RNA was isolated in each experimental condition by Trizol reagent (Invitrogen, Carlsbad, CA, USA), and then purified using the RNeasy MinElute Cleanup Kit (Qiagen, Hilden, Germany). In all cases the integrity of RNA was verified with a Bioanalyzer 2100 (Agilent Technologies, Palo Alto, CA, USA). Amplified biotinylated complementary RNA (2 µg) was produced with an in vitro transcription labeling reaction and was subsequently hybridized onto HU133A GeneChips (Affymetrix, Santa Clara, CA, USA), following the Affymetrix protocol for high-density arrays. Scans were carried out on an Agilent G2500A GeneArray scanner (Agilent Technologies) and the fluorescence intensities of scanned arrays were measured with the Affymetrix GeneChip software. For unsupervised clustering, dChip v1.3 software (dChip, Boston, MA, USA) was used applying a variation filter of p<0.001. For unsupervised analysis, scanned microarray data were introduced in the Gene Expression Pattern Analysis Suite v3.1 (GEPAS). First, genes with intensity levels below 20 were removed. After referencing each condition to its control, data were preprocessed as previously described.19 Briefly, ratios were log-transformed (base 2) and missing values and flat patterns were filtered. The data set was then sent to a clustering tool (Sotatree) that allowed us to select the group of genes that showed different expression profiles between fludarabine-treated and fludarabine-untreated cells, and between fludarabine-sensitive and -resistant cases. In both instances, we considered only the genes whose expression differed by at least two-fold from the respective control. Functional associations between differentially expressed genes were analyzed using the Ingenuity Pathway Analysis platform (Mountain View, CA, USA).
Low-density array study
From each CLL sample, total RNA was extracted from 107 cells using Trizol reagent (Invitrogen). One microgram of total RNA was then retrotranscribed to cDNA using random primers with the High Capacity cDNA Archive Kit (Applied Biosystems, Foster City, CA, USA) according to the manufacturers protocols. Forty genes were selected because of their differential expression shown in the high-density array study and five genes related to nucleoside transport and metabolism [hCNT2, hCNT3, hENT1, hENT2 and deoxycytidine kinase (DCK)] were also included. The selected genes (Online Supplementary Table S1) were studied by real-time polymerase chain reaction (PCR) using TaqMan Low Density Arrays (Applied Biosystems) in an additional series of CLL patients. cDNA samples were subjected to real-time PCR in duplicate in an Abi Prism 7900HT Sequence Detection System (Applied Biosystems). The relative expression of each gene was quantified by the comparative cycle threshold (Ct) method (
Ct), using β-glucuronidase (GUS) as an endogenous control. mRNA expression levels are given as arbitrary quantitative PCR units. For individual patient analysis, the control sample of each gene was taken as the calibrator. For comparisons among different patients, the calibrator was the average Ct value of each gene in all samples grouped together.
|
|
|---|
![]() View larger version (14K): [in a new window] [Download PPT slide] |
Figure 1. Effect of hENT1 inhibition on fludarabine-induced apoptosis in chronic lymphocytic leukemia primary cells. Cells from chronic lymphocytic leukemia patients were incubated with or without fludarabine (1 µg/mL) for 48 hours either in the absence or in the presence of NBTI (100 nM), added 1 hour before the fludarabine treatment. Cytotoxicity was determined by annexin-V/PI staining and results from ten representative chronic lymphocytic leukemia cases are shown. As a positive control of hENT1 inhibition, cells from a mantle cell lymphoma patient were treated for 48 hours with or without gemcitabine (25 µg/mL) in the presence or absence of NBTI (100 nM). Data (mean ± SEM) are given as percentages of the apoptotic rates relative to untreated cells.
|
![]() View larger version (32K): [in a new window] [Download PPT slide] |
Figure 2. Gene expression analysis of fludarabine-sensitive and -resistant chronic lymphocytic leukemia cells. ( A ) Unsupervised hierarchical clustering of the microarray study. The 2102 genes obtained after a variation filter (p<0.001) were visualized by hierarchical clustering and a partial image is shown. Two main branches clearly distinguished sensitive (blue) from resistant (purple) chronic lymphocytic leukemia cases. Four clusters were identified corresponding to the four chronic lymphocytic leukemia patients included in the study, and small subsets of each cluster were constituted by the different experimental conditions from the same patient (C: control, N: NBTI, F: fludarabine, FN: fludarabine + NBTI). Red and green represent up- or down-regulation of a given gene, respectively. (B) Validation of the 20 selected genes differentially expressed in the microarray analysis performed in the same four chronic lymphocytic leukemia cases by low-density arrays. ( C ) Validation in the whole set of samples. The 27 chronic lymphocytic leukemia cases were distributed into four groups according to their fludarabine-response. The SERPINB9 and SERPINE2 genes were differentially expressed depending on fludarabine cytotoxicity. Statistical significance between the most sensitive and the most resistant group was assessed by Students t-test (*p<0.05).
|
Genes implicated in the response of chronic lymphocytic leukemia cells to fludarabine
Changes in gene expression profiles after 3 and 24 hours of fludarabine treatment in the four selected CLL cases (CLL 2, 4, 24 and 26) were analyzed. By performing an unsupervised analysis using the preprocessor tool, we found 61 genes regulated by fludarabine after 24 hours of treatment. This analysis discriminated the two sensitive and the two resistant cases. After removing genes showing inconsistent changes and calculating the mean fold-change of those genes represented more than once on the array, we identified a set of 40 genes as being responsive to fludarabine treatment (Online Supplementary Table S3). Interestingly, none of these genes showed a significant change after 3 hours of treatment. A time-course analysis of the mRNA levels of selected genes (p21, BAX, FAS and TIGAR) did not show significant changes between 3 and 9 hours after treatment (data not shown).
Twenty fludarabine-regulated genes were chosen according to their differential expression (Table 2). These genes were validated in a customized low-density array in 20 cases of CLL. These selected cases are marked in Table 1. All but three genes (TGFB1, Aiolos, MMP9) were actually validated as being regulated by fludarabine (fold-change-ratio higher than 1.4). Moreover, 16 of the 20 genes showed a significant correlation between ex vivo sensitivity to fludarabine and gene-fold induction. These genes encoded proteins involved in the apoptotic machinery (BAX, FAS, TNFSF7, TNFSF9 and TIGAR), cell cycle regulation (cyclin G1, MDM2, p21 and sestrin 1), and also those encoding proteins implicated in DNA synthesis and repair (such as thymidylate synthetase, DDB2 and GADD45A), intracellular protein trafficking (oxysterol binding protein-like 3 and the p53 target zinc finger protein WIG1) and signal transduction, such as PPM1D and MAP4K4. In fact, when the mean fold-induction of the whole panel of up-regulated genes in this set of 20 CLL cases was plotted against the ex vivo cytotoxicity triggered by fludarabine, a highly significant correlation (p=0.0002, r2: 0.55) was found (Figure 3A). It is important to note that TIGAR, a gene recently described as a p53-inducible regulator of apoptosis,22 showed the strongest correlation with sensitivity to fludarabine (p<0.0001 and correlation of 0.60; Figure 3B). Functional analysis of all these probable targets of fludarabine action, using Ingenuity Pathways Analysis software, defined a p53-network of interconnecting genes, most of them (except MMP9) up-regulated after fludarabine treatment, which supports the biological logic of these findings (Online Supplementary Figure S1).
|
View this table: [in a new window] [Download PPT slide] |
Table 2. Selected fludarabine-regulated genes identified by microarray analysis and validated using low-density arrays.
|
![]() View larger version (11K): [in a new window] [Download PPT slide] |
Figure 3. Genes regulated by fludarabine in chronic lymphocytic leukemia (CLL) primary cells. Cytotoxicity of cells from 20 CLL patients to fludarabine 1 µg/mL after 48 hours of treatment was plotted against the mean fold-induction of the whole panel of genes analyzed using low-density arrays (A) and the fold-induction of TIGAR analyzed by low-density arrays (B). Correlation coefficients and p-values are shown.
|
![]() View larger version (24K): [in a new window] [Download PPT slide] |
Figure 4. Role of hENT2 in fludarabine response in chronic lymphocytic leukemia primary cells. (A) For each single case, fold-induction following fludarabine treatment (1 µg/mL) was plotted against the fold-induction after fludarabine 1 µg/mL + NBTI 100 nM treatment. Linear regression is represented for each case and slope values are shown. (B) Mean fold-change after treatment of the 20 genes analyzed using low-density arrays. Each sample is referenced to its control (untreated cells). This is shown by a horizontal line in the figure. (C) Correlation between the percentage of the gene induction inhibited by NBTI 100 nM, and cytotoxicity to fludarabine 1 µg/mL after 48 hours. (D) Correlation between the percentages of the gene induction inhibited by NBTI (100 nM) and hENT2 mRNA level. (E) Relative hENT2 mRNA level of cells from 20 selected chronic lymphocytic leukemia cases was plotted against the fold-induction of TIGAR mRNA after fludarabine treatment. Correlation coefficients and p values are given in the inset boxes.
|
|
|
|---|
Our study was focused on CLL cases bearing a wild type p53 genotype. Since we observed marked variability in the response of these cells to fludarabine, it is likely that other mechanisms participate in the response to this drug in CLL cells. Transcriptomic profiles of fludarabine-sensitive and -resistant CLL patients have been reported by others and most of the genes induced by fludarabine treatment turned out to be p53-dependent.29 Our data are consistent with this view and most of the identified genes are known transcriptional targets of p53. This is illustrated by the network of interacting genes induced by fludarabine (Online Supplementary Figure S1). Some genes, however, were to some extent unexpectedly identified as putative mediators of fludarabine action. This is the case, for instance, for the novel p53-inducible regulator of apoptosis, TIGAR,22 which was consistently up-regulated after fludarabine treatment of CLL cells, showing the highest correlation with sensitivity to fludarabine among all identified genes. TIGAR is the TP53-induced glycolysis and apoptosis regulator, whose up-regulation had been described after adriamycin treatment of various wild type p53 non-hematologic cell lines. Furthermore, p53-mediated upregulation of TIGAR resulted in inhibition of glycolysis and a decrease in cellular levels of reactive oxygen species (ROS).22 Thus, induction of TIGAR after fludarabine treatment of CLL cells may explain why fludarabine induces apoptosis in a caspase-dependent manner irrespective of ROS production.4 Our finding that fludarabine activated a p53 response in all fludarabine-sensitive cases of CLL indicated that other initial events are implicated in the response to this nucleoside analog.
Other proposed mechanisms of chemoresistance to nucleoside analogs are the dysregulation of intracellular enzymes (e.g. kinases, deaminases and nucleotidases) responsible for their metabolism,30 or the levels of proteins implicated in drug transport.12,31 Fludarabine is a hydrophilic compound that does not readily cross plasma membranes by diffusion; functional nucleoside transporters are, therefore, required for cellular entry. Nucleoside transporters may play a major role in facilitating nucleoside-derived drug uptake and their subsequent actions. An increasing body of evidence suggests that the amount of particular types of nucleoside transporters might determine drug-induced cytotoxicity and, accordingly, statistical correlations between the pharmacological effects of nucleoside-derived drugs and expression of nucleoside transporters have been found in a variety of types of cancer so far. Thus, the number of high-affinity binding sites for NBTI (a specific inhibitor of hENT1) has been correlated with cytarabine cytotoxicity in acute myeloid leukemia and acute lymphocytic leukemia cells,32,33 while similar positive correlations have been reported when quantifying hENT1-related mRNA levels.34–36 In fact, hENT1 protein has been reported to be a suitable biomarker for predicting survival in patients suffering from pancreatic adenocarcinoma and treated with gemcitabine (a good hENT1 substrate) as a single therapy.37 Similarly, sensitivity to gemcitabine ex vivo has been reported to be positively correlated with the levels of hENT1 protein in MCL cells.21 In fact, it has been reported that hENT1 expression, despite showing marked variability among individuals, is comparatively higher in most tumors compared to in their non-transformed tissue counterparts38 and that hENT1 protein expression is strongly retained in gynecological tumors, particularly when compared to the highly efficient nucleoside-derived drug transporter protein hCNT1.39
A general consensus in the field is that hENT1 is a ubiquitous, highly expressed transporter protein that can account for most of the nucleoside supply required for salvage processes, particularly during tissue growth. Uncontrolled proliferation processes do not appear to be implicated in the development of CLL and, in accordance with what has been discussed above, hENT1 expression contributes only slightly to fludarabine-induced cytotoxicity in CLL cells,12,14 even though fludarabine is a hENT1 substrate.11 In fact, the levels of hENT2 protein correlate with the ex vivo cytotoxicity induced by this nucleoside analog in CLL cells.12 In this study we were able to show that the transcriptomic response of CLL cells to fludarabine was poorly sensitive to nucleoside transport inhibition triggered by NBTI, a finding that is consistent with hENT2 playing a major role in the genomic response of CLL cells to fludarabine. This is relevant because it provides evidence, at the transcriptomic level, linking hENT2 function and fludarabine-triggered cytotoxicity, thus strongly suggesting that the correlation between hENT2 expression and fludarabine-induced apoptosis is not coincidental but rather the result of a drug-channeled action at the protein level, via hENT2-mediated fludarabine transport. Furthermore, it has also been described that CLL subjects with elevated hCNT3 expression, a concentrative nucleoside transporter that has been shown to mediate cellular entry of fludarabine, had a lower complete response rate to fludarabine therapy.13 However, the current study did not show any relationship between response to fludarabine and expression of hCNT3, although it did provide evidence of a positive correlation between mRNA expression of hENT2 and induction of TIGAR, thus supporting the concept that hENT2-related function might determine the transcriptomic response to fludarabine.
In summary, this is the first study showing a correlation between the recently described p53-inducible apoptosis gene TIGAR and sensitivity to fludarabine in CLL cells. Although all patients analyzed here had wild type p53, they showed variable sensitivity to fludarabine. This finding suggests that, besides the mutational status of p53, early events in drug action, such as its transport across the plasma membrane, might determine cytotoxicity. In fact, evidence is provided of a major role of hENT2-related biological function in the transcriptomic response of CLL cells to fludarabine, a finding that strongly supports the importance of this particular nucleoside transporter in CLL chemotherapy.
# DC and MP-A should be considered co-senior authors. ![]()
The online version of this article contains a supplementary appendix.
ML-G, LT: performed the research, analyzed the data and wrote the paper. MM-A: participated in the analysis of the data and wrote the paper. NV, FJ-C: contributed with analytical tools. EM, EC: revised the manuscript. DC, MP-A: designed and supervised the research, analyzed the data and wrote the paper. All authors revised the manuscript critically and approved the final version to be published. The authors reported no potential conflicts of interest.
Funding: this work was supported in part by grants SAF06/8850 (Ministerio de Educación y Ciencia, Spain), RED 2006-20-014 (Redes Temáticas de Investigación Cooperativa de Cáncer (RTICC, Instituto de Salud Carlos III) and Lymphoma Research Foundation to DC; SAF 2005-01259 and SAF 2008-00577 (Ministerio de Educación y Ciencia, Spain), 2005SGR00315 (Generalitat de Catalunya, Catalonia, Spain) and Fundació La Caixa to MP-A. LT was also funded by a Fundació La Caixa research grant. ML-G is the recipient of a FI pre-doctoral fellowship from Generalitat de Catalunya.
Received for publication April 4, 2008. Revision received July 24, 2008. Accepted for publication August 19, 2008.
|
|
|---|
This article has been cited by other articles:
![]() |
M. Lopez-Guerra, G. Roue, P. Perez-Galan, R. Alonso, N. Villamor, E. Montserrat, E. Campo, and D. Colomer p65 Activity and ZAP-70 Status Predict the Sensitivity of Chronic Lymphocytic Leukemia Cells to the Selective I{kappa}B Kinase Inhibitor BMS-345541 Clin. Cancer Res., April 15, 2009; 15(8): 2767 - 2776. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||