To elucidate interpatient variability in thioguanine nucleotide (TGN) concentrations in acute lymphoblastic leukemia (ALL) cells, we determined the TGN concentrations in leukemic blasts from 82 children with newly diagnosed ALL after intravenous administration of mercaptopurine (MP). Patients treated with MP alone achieved higher TGN concentrations than those treated with the combination of methotrexate plus mercaptopurine (MTX + MP). Analysis of the expression of approximately 9600 genes in ALL cells obtained at diagnosis identified 60 gene probes significantly associated with TGN accumulation in patients treated with MP alone and 75 gene probes in patients treated with MTX + MP, with no overlap between the 2 sets of genes. Genes significantly associated with intracellular TGN accumulation after MP alone included those encoding MP metabolic enzymes and transporters (eg, SLC29A1). Inhibition of SLC29A1 by nitrobenzylmercaptopurine ribonucleoside (NBMPR) caused a 33% to 45% reduction of TGN in ALL cells in vitro (P < .006), consistent with the gene expression findings. Genes associated with TGN concentration after combination therapy included those involved in protein and adenosine triphosphate (ATP)-biosynthesis. Together, these in vivo and in vitro data provide new insight into the genomic basis of interpatient differences in intracellular TGN accumulation and reveal significant differences between treatment with MP alone and treatment with MP and MTX. (Blood. 2005;106:1778-1785)

Mercaptopurine (MP) is a purine antimetabolite widely used in the treatment of acute lymphoblastic leukemia (ALL).1  MP is an inactive prodrug that requires intracellular activation catalyzed by multiple enzymes to exert cytotoxicity (Figure A in Supplemental Document S1, available on the Blood website; see the Supplemental Document link at the top of the online article).2-5  MP is converted by hypoxanthine guanine phosphoribosyltransferase (HGPRT) to thioinosine monophosphate (TIMP) and subsequently to thioguanosine monophosphate (TGMP) by a 2-step process involving inosine monophosphate dehydrogenase (IMPDH) and guanosine monophosphate synthetase (GMPS).2  Subsequently, TGMP can be converted to the active thioguanine nucleotides (TGNs) by adenosine kinase,2,6,7  Cytotoxic effects of MP are achieved primarily through the incorporation of TGN into RNA and DNA.8,9  The incorporation of TGN, as deoxy-TGN triphosphate, inhibits the function of several enzymes involved in DNA replication and repair10,11  and induces DNA damage, such as single-strand breaks, DNA protein cross-links, interstrand cross-links, and sister chromatid exchanges.12-15  The pathway leading to TGN synthesis is in competition with inactivation pathways catalyzed by xanthine oxidase (XDH) and thiopurine methyltransferase (TPMT).16 

Numerous clinical studies have reported substantial interpatient variability in intracellular TGN concentrations during continuation therapy for ALL.5,17-19  Determinants of this variability include differences in MP absorption, genetic polymorphism of TPMT, and patient compliance with treatment.20-23  Interpatient differences in intracellular TGN concentrations have been associated with treatment outcome in children with ALL, whereas variability in MP plasma concentrations has not. Patients with TGN accumulation greater than the population average had better 5-year relapse-free survival compared with patients with TGN concentrations below the population average when they were treated with antimetabolite-based chemotherapy regimens.5,24 

In most treatment protocols for childhood ALL, MP is administered in combination with methotrexate (MTX). The pharmacologic rationale for this combination is based in part on a direct effect of MTX on enzymes involved in TGN metabolism25  or on an indirect effect on plasma purines.17,26,27  Recently, it has been shown that patients treated with this combination in the presence of circulating ALL cells accumulate lower intracellular TGN concentrations in bone marrow leukemic lymphoblasts, peripheral blasts, and erythrocytes than patients treated with MP alone.28  For reasons that are largely unknown, however, there were substantial interindividual differences in TGN accumulation in ALL cells after each treatment (MP alone, 0.21-16.6; MTX + MP, 0.08-1.74; range, 10th-90th percentiles [pmol/5 × 106 cells]).

The current investigation was undertaken to identify genes significantly associated with the accumulation of TGN in ALL cells in vivo and to identify gene expression patterns that differed in leukemia cells that accumulated high compared with low intracellular concentrations of TGN after the administration of MP alone or in combination with MTX.

Patients and treatment

Eighty-two patients aged 18 years or younger with newly diagnosed ALL, enrolled on the St Jude Children's Research Hospital Total Therapy XIIIB protocol,29  were included in the study. The institutional review board approved the study, and informed consent was obtained from patients, parents, or legal guardians, as appropriate. The diagnosis of ALL was based on morphology, cytochemical staining, immunophenotyping, cytogenetic analysis, and genetic characterization using methods previously described.1,30  Patients were randomly assigned to receive 1 of 3 initial therapies (Figure B in Document S1). Patients were administered intravenous MP alone (200 mg/m2 over a period of 20 minutes plus 800 mg/m2 over a period of 5 hours and 40 minutes); low-dose oral MTX (30 mg/m2 every 6 hours for 6 doses) followed by intravenous MP as described (LDMTX + MP); or high-dose intravenous MTX (200 mg/m2 push plus 800 mg/m2 over a 24-hour period) followed by intravenous MP as described (HDMTX + MP). Patients with contraindications to urate oxidase or rasburicase31,32  were administered allopurinol orally. As previously described,28  the use of allopurinol at least 72 hours before bone marrow aspiration did not alter thiopurine metabolite levels among treatment groups. For TGN analysis, we divided patients into 3 groups according to the randomized treatment arms: MP (n = 33), HDMTX + MP (n = 29), and LDMTX + MP (n = 20).

Isolation of ALL blasts from bone marrow aspirate

For all 82 patients, ALL cells at the time of diagnosis and 20 hours after initiation of the MP infusion (after treatment) were obtained and processed as previously described.33  The percentages of leukemia blasts in bone marrow aspirates at diagnosis and after treatment were similar (median, 97% and 95%, respectively). The final cell yield was determined by hemocytometer, and viability was determined by trypan blue exclusion.

Determination of intracellular thiopurine nucleotide concentrations in bone marrow ALL cells

We used a reverse-phase high-performance liquid chromatography (HPLC) method reported previously34  to determine the amount of TGNs in bone marrow cells 20 hours after in vivo MP treatment of patients. Briefly, 5 × 106 ALL cells were sonicated for 10 seconds, and proteins were removed by ultrafiltration. Thiopurine nucleotides were calculated by subtracting the thiopurine ribonucleotides from thiopurine ribonucleosides after conversion with acid phosphatase. The concentration of TGN is expressed as picomoles per 5 × 106 cells.

Mercaptopurine pharmacokinetic analysis

To assess whether plasma MP and MP metabolite concentrations were associated with intracellular TGN concentrations, we measured the plasma concentrations of mercaptopurine, mercaptopurine riboside, thioguanine, thioguanine riboside, and thiouric acid at 3 and 6 hours after the start of MP in 63 evaluable patients (MP, n = 23; HDMTX + MP, n = 17; LDMTX + MP, n = 23) using the HPLC method we have previously reported.34 

TPMT phenotype and genotype

TPMT activity was measured in erythrocytes of 60 evaluable patients at diagnosis, as previously described.35 TPMT genotype was determined using polymerase chain reaction (PCR) methods for the TPMT*2, TPMT*3A, TPMT*3B, and TPMT*3C alleles on somatic cell DNA, as previously described.36 

RNA extraction, gene expression profiling, and quantitative real-time-polymerase chain reaction

Patients for whom gene expression was determined (all 82) had sufficient blasts in their bone marrow aspirates to permit RNA isolation from 5 to 10 × 106 leukemia cells. We extracted high-quality total RNA with TriReagent (MRC, Cincinnati, OH) from mononuclear cells of bone marrow aspirates at diagnosis. As previously described,37,38  high-quality total RNA was processed and hybridized to the HG-U95Av2 oligonucleotide microarray (Affymetrix; Santa Clara, CA; see manufacturer's manual for details).

To further establish the validity of gene expression determined by microarray analysis, we performed quantitative real-time-polymerase chain reaction (RT-PCR) in 4 patient samples for 3 genes (GAB1, GRB2-associated binding protein 1; HNRPF, heterogeneous nuclear ribonucleoprotein F; DPM2, dolichol-phosphate mannosyltransferase polypeptide 2, regulatory subunit).

Determination of TGN accumulation in leukemia cells in the presence or absence of a specific inhibitor of the SLC29A1 transporter

We used the reverse-phase HPLC method reported previously34  to determine intracellular TGN concentration after treatment with 10 μM MP (Sigma, St Louis, MO) in the presence or absence of a specific SLC29A1 inhibitor, nitrobenzylmercaptopurine ribonucleoside (NBMPR) (Sigma). CEM, Nalm-6, and MOLT-4 human leukemia cells (American Type Culture Collection [ATCC], Manassas, VA) (50 × 106 cells in 50 mL media) were incubated for 24 hours at 37°C, and the effect of SLC29A1 transport inhibition was assessed by adding 100 nM NBMPR 20 minutes before MP exposure. As a negative control, we treated CEM cells with 1 nM adenine (Sigma), which is not a substrate for SLC29A1,39,40  and measured the intracellular adenine concentration after 24-hour exposure in the presence or absence of 100 nM NBMPR. These experiments were performed 3 times with 3 replicates.

Statistical analysis

Pairwise comparisons using the Wilcoxon rank sum test were used to assess differences in posttreatment TGN concentrations among patients in the 3 treatment arms. The Holm method was used to adjust for multiple testing.41  The Wilcoxon rank sum test was used to compare plasma MP and MP-metabolite concentrations between patients randomly assigned to HDMTX + MP and LDMTX + MP. Multiple regression analysis was used to assess the association between intracellular TGN concentrations and plasma MP and MP-metabolite concentrations separately in the group of patients treated with MP alone and in that treated with the combination of MTX + MP. Spearman rank correlation was used to assess the association between TPMT activity and TGN concentration. Default settings of Affymetrix Microarray Suite software version 5 (MAS 5.0; Affymetrix) were used to calculate gene expression values. Expression values were scaled to the target intensity of 2500 and were log-transformed. From the total of 12 599 probe sets, those expressed in less than 5% of the 82 patients were omitted, leaving 5013 probe sets for subsequent analyses. Analyses were performed separately in the 2 groups of patients treated either with MP alone or with the combination of MTX plus MP. Spearman rank correlation and false discovery rate (FDR) using the Storey q-value42  were computed to identify genes associated with TGN concentration. Leave-one-out cross-validation using the support vector machine (SVM) based on Spearman rank correlation was used to assess the validity of association between gene expression profiles and TGN concentrations.

To assess whether the expression of gene probe sets associated with intracellular TGN concentration in patients treated with MP alone was also associated with MP sensitivity in the 60 human tumor cell line panel (NCI60), we correlated the expression levels of those genes with the growth inhibition values (GI50) obtained after treatment with MP in these 60 cell lines. The GI50 values and the HG-U95Av2 gene expression data are publicly available through the National Cancer Institute's Developmental Therapeutics Program (http://dtp.nci.nih.gov/mtargets/download.html). Spearman rank correlation was used to assess the association between gene expression and GI50 values.

R 1.8.0 statistical software (http://www.r-project.org) was used to perform these analyses. Principal component analysis (PCA) and 2-dimensional hierarchical clustering were performed using GeneMaths 2.1 software (AppliedMaths, St Martens-Latem, Belgium).

Differences in TGN accumulation in ALL cells of patients randomly assigned to 3 initial treatments

As shown in Figure 1, 33 patients who received MP alone (2.46 ± 5.36 pmol/5 × 106 cells; median ± SD) accumulated significantly higher concentrations of TGN in ALL cells than 29 patients treated with either LDMTX + MP (0.51 ± 0.64 pmol/5 × 106 cells; P < .001) or 20 patients treated with HDMTX + MP (0.57 ± 0.54 pmol/5 × 106 cells; P < .001). There were no differences in TGN concentrations between patients randomly assigned to receive LDMTX + MP and those randomly assigned to receive HDMTX + MP (P = .30). Patient characteristics (age, sex, race, white blood cell [WBC] count, ALL subtype) were similar among the 3 groups of patients (Figure B and Table A in Document S1).

Therefore, because of the results of our current analysis and our previous findings that showed no difference in TGN accumulation when MP was given after low-dose and high-dose MTX,28  in subsequent analyses we combined patients treated with either combination of HDMTX + MP and LDMTX + MP (MTX + MP; n = 49) into one group. Multiple regression analyses did not show any association between TGN concentration and patient age, sex, initial WBC count, or ALL subtype (Table B in Document S1) in either group of patients.

Figure 1.

Intracellular concentrations of TGNs in patients after treatment with MP alone or in combination with MTX. Box plot of TGN concentrations in ALL cells for patients randomly assigned to treatment with MP alone (n = 33), HDMTX + MP (n = 20), or LDMTX + MP (n = 29). Horizontal line indicates the median for each group; box depicts the 25th and 75th percentile ranges. ⋄ indicates outlier values. P was determined by pairwise comparison using the Wilcoxon rank sum test adjusted for multiple testing.

Figure 1.

Intracellular concentrations of TGNs in patients after treatment with MP alone or in combination with MTX. Box plot of TGN concentrations in ALL cells for patients randomly assigned to treatment with MP alone (n = 33), HDMTX + MP (n = 20), or LDMTX + MP (n = 29). Horizontal line indicates the median for each group; box depicts the 25th and 75th percentile ranges. ⋄ indicates outlier values. P was determined by pairwise comparison using the Wilcoxon rank sum test adjusted for multiple testing.

Close modal

Correlation between MP pharmacokinetic parameters and intracellular TGN accumulation

Plasma concentrations of MP and MP metabolites at 3 and 6 hours after MP treatment were similar between patients randomly assigned to treatment with HDMTX + MP or LDMTX + MP.28  Plasma concentrations of thioguanine and thioguanine riboside were negligible (data not shown). Multiple regression analyses did not show any significant association between plasma concentrations of mercaptopurine, mercaptopurine riboside, or thiouric acid and intracellular TGN accumulation in ALL cells after treatment with MP alone (P > .29; n = 23) or with the combination of MTX + MP (P > .20; n = 40) (Table C in Document S1). The power to detect an association of R2 = 0.20 at the α= 0.10 level for 23 or 40 persons was 68.7% or 90%, respectively.

Association between TPMT phenotype/genotype and TGN accumulation in bone marrow ALL cells

In patients treated with MP alone (n = 24) or in patients treated with the combination of MTX + MP (n = 36), intracellular TGN concentrations were not associated with TPMT activity measured at the time of diagnosis (MP, R2 = 0.01; MTX + MP, R2 = 0.02). However, among the 60 patients who underwent genotyping, only 2 had TPMT heterozygosity and none were deficient in TPMT.

Relation between gene expression and intracellular TGN accumulation

Unsupervised hierarchical clustering, which groups patients based on predominant similarities in overall gene expression, did not cluster patients according to ALL TGN concentrations (Figure C in Document S1). However, Spearman rank correlation and leave-one-out cross-validation using SVM identified 60 gene probe sets (ρ= 0.60; P < .001; FDR = 45%) (Table D in Document S1) that were associated with TGN concentrations in ALL cells after treatment with MP alone (n = 33). Hierarchical clustering using these selected gene probe sets distinctly separated the 33ALL patients treated with MP alone into 2 major groups with either a low (less than 2.46 pmol/5 × 106 cells) or a high (more than 2.46 pmol/5 × 106 cells) TGN concentration (Figure 2A). PCA illustrates the degree of discrimination between the 2 groups of patients (Figure 2B).

In a separate analysis, 75 gene probe sets (ρ= 0.65; P < .001; FDR = 25%) (Table E in Document S1) were significantly associated with TGN concentrations in ALL cells after treatment MTX + MP (n = 49). Two major branches were also found for patients treated with MTX + MP (less than 0.55 pmol/5 × 106 cells and more than 0.55 pmol/5 × 106 cells; Figure 3A), and the degree of discrimination between these 2 groups of patients is evident by PCA (Figure 3B).

Interestingly, there was no overlap between the 2 sets of genes. Among genes associated with TGN accumulation in ALL cells after treatment with MP alone, 31 were positively and 29 were negatively correlated with TGN concentration. In the group of patients treated with MTX + MP, 50 genes were positively correlated and 25 were negatively correlated with TGN accumulation. Genes significantly associated with TGN concentration in patients treated with MP alone included xanthine oxidase (XDH)(ρ= -0.51; P = .002), solute carrier family 29 member 1 (SLC29A1) (ρ= 0.54; P = .001), aldehyde dehydrogenase 1 family member A2 (ALDH1A2)(ρ= 0.56; P < .001), adenosine deaminase (ADA)(ρ= 0.48; P = .004), and genes related to cell proliferation and apoptosis (caspase 7 [CASP7], DNA topoisomerase II-binding protein [TOPBP1], anaphase-promoting complex subunit 5 [ANAPC5], and chaperonin-containing TCP1 subunit 4 [CCT4]). In the group of patients treated with MTX + MP, we found genes involved in adenosine triphosphate (ATP) synthesis (solute carrier family 25, mitochondrial carrier member 3 [SLC25A3], ATP synthetase, H+ transporting, mitochondrial complex, O subunit [ATP50], cytochrome c oxidase subunit Vb [COX5B], and VIIa polypeptide 2-like [COX7A2L]) and several genes implicated in protein synthesis, such as ribosomal proteins (RPS19, RPL18, RPS25, RPL23) and translation factors (EEF1G, EIF3S5, eIF3k).

Figure 2.

Clustering of genes associated with intracellular concentrations of TGN in leukemia cells of patients treated with MP alone. (A) Hierarchical clustering of 33 ALL patients using 60 gene probe sets significantly associated with TGN accumulation in ALL cells after treatment with MP alone. Each row represents a probe set and each column a patient. Standardized expression values are shown according to the scale. Red and green symbols indicate patients with TGN concentrations higher and lower than 2.46 pmol/5 × 106 cells, respectively. Listed are the Affymetrix probe set IDs, gene names, GenBank accession numbers, and coefficients of correlation (rho). A positive rho indicates a positive correlation, and a negative rho indicates negative correlation between gene expression and intracellular TGN concentration. (B) PCA of ALL cells with high (red) and low (green) TGN concentrations based on the selected probe sets. Similarities are visualized in 3-dimensional space.

Figure 2.

Clustering of genes associated with intracellular concentrations of TGN in leukemia cells of patients treated with MP alone. (A) Hierarchical clustering of 33 ALL patients using 60 gene probe sets significantly associated with TGN accumulation in ALL cells after treatment with MP alone. Each row represents a probe set and each column a patient. Standardized expression values are shown according to the scale. Red and green symbols indicate patients with TGN concentrations higher and lower than 2.46 pmol/5 × 106 cells, respectively. Listed are the Affymetrix probe set IDs, gene names, GenBank accession numbers, and coefficients of correlation (rho). A positive rho indicates a positive correlation, and a negative rho indicates negative correlation between gene expression and intracellular TGN concentration. (B) PCA of ALL cells with high (red) and low (green) TGN concentrations based on the selected probe sets. Similarities are visualized in 3-dimensional space.

Close modal

Correlation between gene expression and GI50 using NCI60 cancer cell lines

Using Spearman rank correlation, we found that among the 60 gene probe sets associated with TGN accumulation in ALL cells after in vivo treatment with MP alone, 10 genes were significantly associated with the 50% growth inhibitory concentration (GI50) of MP in vitro sensitivity (P < .05) measured in a panel of 60 human tumor cell lines. Nine genes were negatively and one was positively correlated with GI50. Six of the negatively correlated genes were positively correlated with TGN accumulation in ALL cells in our patients. These genes encode for transporters (SLC29A1, ATP2A3), proteins involved in oxidative phosphorylation (NADH dehydrogenase flavoprotein 1 [NDUFV1]), and carbohydrate metabolism (isocitrate dehydrogenase 3 β [IDH3B] and DNA replication ([TOPBP1]) (Table 1; Figure D in Document S1).

Table 1.

Genes significantly associated with GI50for MP in NCI60 human tumor cell lines


Probe set ID

Gene name

P GI50

ρ GI50

P TGN

ρ TGN
35835_at   Anaphase promoting complex subunit 5*  < .001   −.45   .007   .47  
41544_at   Serum-inducible kinase   .001   .41   .004   .49  
37329_at   NADH dehydrogenase flavoprotein 1*  .003   −.39   .001   .55  
33901_at   Solute carrier family 29, member 1*  .003   −.38   .003   .54  
32007_at  Homo sapiens-transcribed sequences   .004   −.37   .002   −.52  
38834_at   Topoisomerase (DNA) II binding protein*  .012   −.32   .001   .58  
37825_at   Galactokinase 2   .016   −.31   .007   −.46  
36482_s_at   ATPase, Ca++ transporting, ubiquitous*  .022   −.30   .004   .50  
40111_g_at   Isocitrate dehydrogenase 3 (NAD+) β*  .031   −.28   .003   .51  
1937_at
 
Retinoblastoma 1
 
.046
 
−.26
 
.004
 
−.49
 

Probe set ID

Gene name

P GI50

ρ GI50

P TGN

ρ TGN
35835_at   Anaphase promoting complex subunit 5*  < .001   −.45   .007   .47  
41544_at   Serum-inducible kinase   .001   .41   .004   .49  
37329_at   NADH dehydrogenase flavoprotein 1*  .003   −.39   .001   .55  
33901_at   Solute carrier family 29, member 1*  .003   −.38   .003   .54  
32007_at  Homo sapiens-transcribed sequences   .004   −.37   .002   −.52  
38834_at   Topoisomerase (DNA) II binding protein*  .012   −.32   .001   .58  
37825_at   Galactokinase 2   .016   −.31   .007   −.46  
36482_s_at   ATPase, Ca++ transporting, ubiquitous*  .022   −.30   .004   .50  
40111_g_at   Isocitrate dehydrogenase 3 (NAD+) β*  .031   −.28   .003   .51  
1937_at
 
Retinoblastoma 1
 
.046
 
−.26
 
.004
 
−.49
 
*

Genes with negative correlation with NCI60 GI50 and positive correlation with TGN concentration.

TGN accumulation in leukemia cells differed after in vitro treatment with MP alone or in the presence of NBMPR

NBMPR is an inhibitor of the SLC29A1 nucleoside transporter.43  As depicted in Figure 4, in CEM human leukemia cells treated for 24 hours with MP + NBMPR, concentrations of TGN were 45% lower in cells treated with MP alone (672 ± 160 vs 1204 ± 262 pmol/5 × 106 cells; P < .002), 33% lower in Nalm-6 cells (764 ± 153 vs 1133 ± 146 pmol/5 × 106 cells; P = .001), and 42% lower in MOLT-4 cells (156 ± 62 vs 267 ± 44 pmol/5 × 106 cells; P = .006), comparing use of inhibitor with no inhibitor, respectively, for each cell line. The intracellular accumulation of adenine, which is not a substrate of SLC29A1,39  was not different in CEM cells treated with adenine in the presence or the absence of NBMPR (47.44 ± 8.6 vs 44.6 ± 12.6 pmol/5 × 106 cells; P = .65).

Validation of gene expression

Gene expression levels determined by microarray analysis were validated using quantitative real-time-polymerase chain reaction (RT-PCR) in 4 patient samples for 3 genes (GAB1, HNRPF, DPM2). Correlation (R2) between mRNA levels determined by RT-PCR and microarray was 0.97 for GAB1, 0.77 for HNRPF, and 0.83 for DPM2, confirming expression levels determined by the gene expression array (Figure E in Document S1).

Intracellular accumulation of TGNs is a critical event for the cytotoxic effects of mercaptopurine.8,9  Previous studies have shown that several factors, including the activity of enzymes involved in MP metabolism, MP dose intensity, drug-drug interactions, and patient compliance, significantly influence interpatient differences in TGN concentrations.18,44,45  Furthermore, these interpatient differences have been associated with treatment efficacy3,5,24  and MP toxicity.23,46 

Well-defined sources of variability include inherited differences in TPMT36,46-48  and inhibition of xanthine oxidase by allopurinol.49  Consistent with previous reports from our laboratory and others,26,28  the current study revealed a significant in vivo effect of MTX on intracellular accumulation of TGN compared with treatment with MP given alone (P < .001). The current study also documented substantial interpatient variability in TGN accumulation not explained by these drug interactions or by inherited differences in TPMT activity.

Figure 3.

Hierarchical clustering using genes associated with intracellular concentrations of TGN in ALL cells of patients treated with MP + MTX. (A) Hierarchical clustering of 49 ALL patients using 75 probe sets significantly associated with TGN accumulation in ALL cells after treatment with MTX + MP. Each row represents a probe set and each column a patient. Standardized expression values are shown according to the scale. Red symbols depict patients with TGN concentrations greater than 0.55 pmol/5 × 106 cells, and green symbols depict those with concentrations less than 0.55 pmol/5 × 106 cells. Listed are the Affymetrix probe set IDs, gene names, GenBank accession numbers, and coefficients of correlation (rho). A positive rho indicates a positive correlation, and a negative rho indicates a negative correlation between gene expression and intracellular TGN concentration. (B) PCA of ALL cells with high (red) and low (green) TGN concentration based on the selected probe sets. Similarities are visualized in 3-dimensional space.

Figure 3.

Hierarchical clustering using genes associated with intracellular concentrations of TGN in ALL cells of patients treated with MP + MTX. (A) Hierarchical clustering of 49 ALL patients using 75 probe sets significantly associated with TGN accumulation in ALL cells after treatment with MTX + MP. Each row represents a probe set and each column a patient. Standardized expression values are shown according to the scale. Red symbols depict patients with TGN concentrations greater than 0.55 pmol/5 × 106 cells, and green symbols depict those with concentrations less than 0.55 pmol/5 × 106 cells. Listed are the Affymetrix probe set IDs, gene names, GenBank accession numbers, and coefficients of correlation (rho). A positive rho indicates a positive correlation, and a negative rho indicates a negative correlation between gene expression and intracellular TGN concentration. (B) PCA of ALL cells with high (red) and low (green) TGN concentration based on the selected probe sets. Similarities are visualized in 3-dimensional space.

Close modal

We evaluated MP with or without MTX as initial therapy because that is the only time it is possible to assess MP dispositions and effects in primary cells after in vivo treatment (by week 4, most patients have no ALL cells in their peripheral blood or bone marrow). In our study, we were focused on treatment-naive ALL and its ability to accumulate active MP metabolites in vivo. This necessitated that we assess gene expression at diagnosis and TGN accumulation after in vivo MP treatment during the first 4 days of treatment.

To better understand the biologic basis underlying interpatient variability in TGN accumulation in ALL cells after MP treatment, we used oligonucleotide microarrays to perform a genomewide determination of pretreatment gene expression in ALL cells to identify genes and gene expression patterns significantly associated with high or low TGN accumulation in leukemia cells. This revealed 60 gene probe sets (49 genes, 11 cDNA clones) highly related to the extent of TGN accumulation in leukemia cells of patients treated with MP alone and 75 gene probe sets (68 genes, 7 cDNA clones) of patients treated with MTX + MP. Notably, there was no overlap between the 2 sets of genes, indicating that different genes influence the accumulation of TGN when MP is given alone and when it is given in combination with MTX. Previous studies have shown that in vivo gene expression changes are treatment specific in ALL cells and that the nature of these cellular drug responses is different when medications are given alone and when they are given in combination with a second agent.37,50  Of note, gene expression was not biased by low-dose or high-dose MTX + MP shown by “unsupervised” clustering analysis using all genes and not by using genes related to TGN accumulation (Figure C in Document S1; P = .5).

Figure 4.

Effect of concomitant treatment with nucleoside transporter inhibitor NBMPR on TGN and adenine concentrations in ALL cells. (Left) Shown are intracellular TGN concentrations in CEM, Nalm-6, and MOLT-4 cells after treatment with MP in the presence or absence of nitrobenzylmercaptopurine ribonucleoside (NBMPR), a nucleoside transport (SLC29A1) inhibitor. (Right) Shown is the intracellular adenine concentration in CEM cells after incubation with adenine, in the presence or absence of NBMPR. Columns represent the averages of 3 independent experiments, each with 3 replicates; error bars represent range.

Figure 4.

Effect of concomitant treatment with nucleoside transporter inhibitor NBMPR on TGN and adenine concentrations in ALL cells. (Left) Shown are intracellular TGN concentrations in CEM, Nalm-6, and MOLT-4 cells after treatment with MP in the presence or absence of nitrobenzylmercaptopurine ribonucleoside (NBMPR), a nucleoside transport (SLC29A1) inhibitor. (Right) Shown is the intracellular adenine concentration in CEM cells after incubation with adenine, in the presence or absence of NBMPR. Columns represent the averages of 3 independent experiments, each with 3 replicates; error bars represent range.

Close modal

After treatment with MP alone, genes significantly related to TGN accumulation in ALL cells included those encoding transporters (SLC29A1, ATP2A3, NDUFS1), enzymes involved in purine metabolism (ADA, XDH), regulation of cell cycle (ANAPC5, CCT4), and apoptosis (CASP7, PRKR). Interestingly, only 2 (XDH, ADA) of the genes identified in the current study have been previously associated with MP metabolism and disposition. Expression of SLC29A1 (or ENT1), a member of the family of the equilibrative nucleoside transporters that mediates the uptake of nucleosides and chemotherapeutic agents,39,51  has been shown to be positively associated with cytotoxicity of nucleoside analogs (MP, azacytidine, cytarabine, inosine-glycodialdehyde) in human cancer cell lines52  and clinical response to cytarabine in patients with acute myeloid leukemia.53  Consistent with this, in our patients, SLC29A1 expression was high in ALL cells that accumulated high posttreatment TGN concentrations (ρ= 0.54; P = .0015).

We measured TGN accumulations in CEM, Nalm-6, and MOLT-4 human leukemia cells after exposure to MP in the presence or absence of NBMPR, a specific inhibitor of the SLC29A1 transporter. Inhibition of SLC29A1 by NBMPR reduced the concentrations of TGN at 24 hours by 45%, 33%, and 42%, respectively (P < .006). Together, these in vivo and in vitro data provide the first evidence that the inhibition of SLC29A1 reduced steady state TGN accumulation by 45% to 33%, and it also indicates the existence of other MP and TGN transport and metabolic mechanisms in ALL cells.

XDH encodes xanthine oxidase, an important enzyme involved in the first-pass metabolism of MP.16  This enzyme, predominantly expressed in the intestinal mucosa and liver, converts MP to inactive thiouric acid,54,55  a principal pathway for MP catabolism. In ALL, XDH was expressed at a relatively low level (median, 1575, Affymetrix scaled signal). The significant negative correlation between XDH expression and TGN concentration in ALL cells (ρ= -0.52; P = .003) is consistent with the role of this enzyme in the catabolism of MP, but it may be related to concordant variability of XDH in other tissues (eg, liver). The nonsignificant correlation of XDH expression and TGN (ρ= 0.08; P = .54) in patients treated with MTX + MP may be attributed to the inhibitory effects of MTX on xanthine oxidase activity.25 

Another gene that showed a positive relation between gene expression and TGN accumulation (ρ= 0.48; P = .005) was human adenosine deaminase (ADA). ADA encodes a purine-metabolizing enzyme that indirectly influences the metabolism of MP through perturbation of the intracellular pool of enzymes involved in nucleotide synthesis. ADA has been extensively studied in ALL patients, and low ADA activity has been associated with poor outcome after in vivo treatment with the structurally similar 6-thioguanine.56  It is postulated that high ADA expression results in the accumulation of inosine and subsequently inosine monophosphate (IMP),57  which inhibits phosphoribosyl pyrophosphate aminotransferase (PPAT), thereby increasing intracellular phosphoribosyl pyrophosphate (PRPP), a rate-limiting precursor for the formation of TGN.58 

Interestingly, among the 60 gene probe sets significantly associated with TGN accumulation in patients treated with MP alone, 10 were significantly correlated with MP sensitivity (GI50) in NCI60 human tissue cell lines. Notably, 6 of these genes were negatively correlated with GI50 in cell lines and positively correlated with TGN accumulation in our patients. These findings indicate that increased expression of genes associated with TGN accumulation is associated with increased sensitivity to MP.

Among genes related to TGN accumulation after MTX + MP treatment were several involved in mitochondrial activities (eg, ATP biosynthesis and fatty acid metabolism) and protein synthesis (eg, ribosomal proteins and translation factors). To our knowledge, none of those 68 genes has been previously implicated in MP metabolism, TGN accumulation, or thiopurine resistance. Four of these genes, SLC25A3, ATP50, COX5B, and COX7A2L, are involved in the synthesis of the structurally similar nucleotide, ATP. SLC25A3 encodes the phosphate carrier PCH, an integral protein of the mitochondrial inner membrane involved in the terminal process of oxidative phosphorylation.59 ATP50 encodes for the oligomycin sensitivity conferring protein (OSCP) subunit of the stalk of the mitochondrial respiratory chain F1F0-ATP synthase.60  Finally, COX5B and COX7A2L encode for subunits of the cytochrome c oxidase, the terminal enzyme of the mitochondrial electron transport chain. ATP is a cofactor for GMPS,61  an important enzyme involved in the conversion of thioxanthosine monophosphate (TXMP) to TGMP, an essential step in TGN formation. Low ATP levels have been associated with significantly lower intracellular TGN accumulation in leukemia cells treated in vitro with MP.62  MTX, through the inhibition of de novo purine synthesis, decreases intracellular levels of ATP.44,63  Therefore, it is plausible that patients with increased expression of SLC25A3, ATP50, COX5B, and COX7A2L have higher ATP levels and are, therefore, less susceptible to the inhibition of ATP synthesis by MTX and, as a consequence, achieve relatively higher intracellular TGN concentrations than other patients treated with MTX + MP. Expression of these genes was not significantly associated with TGN concentrations in the group of patients treated with MP alone. Together, these in vivo and in vitro data provide new insight into the genomic basis of interpatient differences in TGN accumulation in ALL cells and reveal significant differences when MP is given alone or in combination with MTX.

Prepublished online as Blood First Edition Paper, May 19, 2005; DOI 10.1182/blood-2005-01-0143.

Supported in part by National Institutes of Health grants R37 CA36401 (W.E.E., M.V.R., C.H.P.), R01 CA78224 (W.E.E., M.V.R., C.H.P.), RO1 CA51001 (M.V.R., C.H.P.), and U01 GM61393 (M.V.R., W.E.E.) and Cancer Center Support Grant CA21765; by an F.M. Kirby Clinical Research Professorship from the American Cancer Society (C.H.P.); and by the American Lebanese Syrian Associated Charities.

The online version of the article contains a data supplement.

An Inside Blood analysis of this article appears in the front of this issue.

The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 U.S.C. section 1734.

We thank the following people for their outstanding support: Dr Natalia Krynetskaia, Dr Evgeny Krynetsky, Yaqin Chu, Margaret Needham, Nancy Kornegay, and Dr Clayton Naeve and his staff in the Hartwell Center for Bioinformatics and Biotechnology at St. Jude Children's Research Hospital. We also thank Paxton Baker, Cong Ding, Nancy Duran, Emily Melton, and Siamac Salehy for excellent technical assistance.

1
Pui CH, Evans WE. Acute lymphoblastic leukemia.
N Engl J Med.
1998
;
339
:
605
-615.
2
Elion GB. The purine path to chemotherapy.
Science.
1989
;
244
:
41
-47.
3
Lennard L, Lilleyman JS. Variable mercaptopurine metabolism and treatment outcome in childhood lymphoblastic leukemia.
J Clin Oncol.
1989
;
7
:
1816
-1823.
4
Evans WE, Horner M, Chu YQ, Kalwinsky D, Roberts WM. Altered mercaptopurine metabolism, toxic effects, and dosage requirement in a thiopurine methyltransferase-deficient child with acute lymphocytic leukemia.
J Pediatr.
1991
;
119
:
985
-989.
5
Lennard L, Lilleyman JS, Van Loon J, Weinshilboum RM. Genetic variation in response to 6-mercaptopurine for childhood acute lymphoblastic leukaemia.
Lancet.
1990
;
336
:
225
-229.
6
Keuzenkamp-Jansen CW, De Abreu RA, Bokkerink JP, Trijbels JM. Determination of extracellular and intracellular thiopurines and methylthiopurines by high-performance liquid chromatography.
J Chromatogr.B Biomed.Appl.
1995
;
672
:
53
-61.
7
Zimmerman TP, Chu LC, Bugge CJ, et al. Identification of 6-methylmercaptopurine ribonucleoside 5'-diphosphate and 5'-triphosphate as metabolites of 6-mercaptopurine in man.
Cancer Res.
1974
;
34
:
221
-224.
8
Warren DJ, Andersen A, Slordal L. Quantitation of 6-thioguanine residues in peripheral blood leukocyte DNA obtained from patients receiving 6-mercaptopurine-based maintenance therapy.
Cancer Res.
1995
;
55
:
1670
-1674.
9
LePage GA. Basic biochemical effects and mechanism of action of 6-thioguanine.
Cancer Res.
1963
;
23
:
1202
-1206.
10
Krynetskaia NF, Krynetski EY, Evans WE. Human RNase H-mediated RNA cleavage from DNA-RNA duplexes is inhibited by 6-deoxythioguanosine incorporation into DNA.
Mol Pharmacol.
1999
;
56
:
841
-848.
11
Krynetskaia NF, Cai X, Nitiss JL, Krynetski EY, Relling MV. Thioguanine substitution alters DNA cleavage mediated by topoisomerase II.
FASEB J.
2000
;
14
:
2339
-2344.
12
Ling YH, Chan JY, Beattie KL, Nelson JA. Consequences of 6-thioguanine incorporation into DNA on polymerase, ligase, and endonuclease reactions.
Mol Pharmacol.
1992
;
42
:
802
-807.
13
Bodell WJ. Molecular dosimetry of sister chromatid exchange induction in 9L cells treated with 6-thioguanine.
Mutagenesis.
1991
;
6
:
175
-177.
14
Christie NT, Drake S, Meyn RE, Nelson JA. 6-Thioguanine-induced DNA damage as a determinant of cytotoxicity in cultured Chinese hamster ovary cells.
Cancer Res.
1984
;
44
:
3665
-3671.
15
Pan BF, Nelson JA. Characterization of the DNA damage in 6-thioguanine-treated cells.
Biochem Pharmacol.
1990
;
40
:
1063
-1069.
16
Bostrom B, Erdmann G. Cellular pharmacology of 6-mercaptopurine in acute lymphoblastic leukemia.
Am J Pediatr Hematol Oncol.
1993
;
15
:
80
-86.
17
Balis FM, Holcenberg JS, Poplack DG, et al. Pharmacokinetics and pharmacodynamics of oral methotrexate and mercaptopurine in children with lower risk acute lymphoblastic leukemia: a joint children's cancer group and pediatric oncology branch study.
Blood.
1998
;
92
:
3569
-3577.
18
Relling MV, Hancock ML, Boyett JM, Pui CH, Evans WE. Prognostic importance of 6-mercaptopurine dose intensity in acute lymphoblastic leukemia.
Blood.
1999
;
93
:
2817
-2823.
19
Schmiegelow K, Bruunshuus I. 6-Thioguanine nucleotide accumulation in red blood cells during maintenance chemotherapy for childhood acute lymphoblastic leukemia, and its relation to leukopenia.
Cancer Chemother.Pharmacol.
1990
;
26
:
288
-292.
20
Lennard L, Van Loon JA, Lilleyman JS, Weinshilboum RM. Thiopurine pharmacogenetics in leukemia: correlation of erythrocyte thiopurine methyltransferase activity and 6-thioguanine nucleotide concentrations.
Clin Pharmacol Ther.
1987
;
41
:
18
-25.
21
Lennard L, Lilleyman JS. Individualizing therapy with 6-mercaptopurine and 6-thioguanine related to the thiopurine methyltransferase genetic polymorphism.
Ther Drug Monit.
1996
;
18
:
328
-334.
22
Evans WE, Relling MV. Moving towards individualized medicine with pharmacogenomics.
Nature.
2004
;
429
:
464
-468.
23
Evans WE, Hon YY, Bomgaars L, et al. Preponderance of thiopurine S-methyltransferase deficiency and heterozygosity among patients intolerant to mercaptopurine or azathioprine.
J Clin Oncol.
2001
;
19
:
2293
-2301.
24
Lilleyman JS, Lennard L. Mercaptopurine metabolism and risk of relapse in childhood lymphoblastic leukaemia.
Lancet.
1994
;
343
:
1188
-1190.
25
Lewis AS, Murphy L, McCalla C, Fleary M, Purcell S. Inhibition of mammalian xanthine oxidase by folate compounds and amethopterin.
J Biol Chem.
1984
;
259
:
12
-15.
26
Balis FM, Holcenberg JS, Zimm S, et al. The effect of methotrexate on the bioavailability of oral 6-mercaptopurine.
Clin Pharmacol.Ther.
1987
;
41
:
384
-387.
27
Innocenti F, Danesi R, Di Paolo A, et al. Clinical and experimental pharmacokinetic interaction between 6-mercaptopurine and methotrexate.
Cancer Chemother Pharmacol.
1996
;
37
:
409
-414.
28
Dervieux T, Hancock ML, Pui CH, et al. Antagonism by methotrexate on mercaptopurine disposition in lymphoblasts during up-front treatment of acute lymphoblastic leukemia.
Clin Pharmacol.Ther.
2003
;
73
:
506
-516.
29
Pui CH, Sandlund JT, Pei D, et al. Improved outcome for children with acute lymphoblastic leukemia: results of Total Therapy Study XIIIB at St Jude Children's Research Hospital.
Blood.
2004
;
104
:
2690
-2696.
30
Silverman LB, Gelber RD, Dalton VK, et al. Improved outcome for children with acute lymphoblastic leukemia: results of Dana-Farber Consortium Protocol 91-01.
Blood.
2001
;
97
:
1211
-1218.
31
Pui CH, Relling MV, Lascombes F, et al. Urate oxidase in prevention and treatment of hyperuricemia associated with lymphoid malignancies.
Leukemia.
1997
;
11
:
1813
-1816.
32
Pui CH, Mahmoud HH, Wiley JM, et al. Recombinant urate oxidase for the prophylaxis or treatment of hyperuricemia in patients with leukemia or lymphoma.
J Clin Oncol.
2001
;
19
:
697
-704.
33
Dervieux T, Brenner TL, Hon YY, et al. De novo purine synthesis inhibition and antileukemic effects of mercaptopurine alone or in combination with methotrexate in vivo.
Blood.
2002
;
100
:
1240
-1247.
34
Dervieux T, Chu Y, Su Y, et al. HPLC determination of thiopurine nucleosides and nucleotides in vivo in lymphoblasts following mercaptopurine therapy.
Clin Chem.
2002
;
48
:
61
-68.
35
Dervieux T, Blanco JG, Krynetski EY, et al. Differing contribution of thiopurine methyltransferase to mercaptopurine versus thioguanine effects in human leukemic cells.
Cancer Res.
2001
;
61
:
5810
-5816.
36
Yates CR, Krynetski EY, Loennechen T, et al. Molecular diagnosis of thiopurine S-methyltransferase deficiency: genetic basis for azathioprine and mercaptopurine intolerance.
Ann Intern Med.
1997
;
126
:
608
-614.
37
Cheok MH, Yang W, Pui CH, et al. Treatment-specific changes in gene expression discriminate in vivo drug response in human leukemia cells.
Nat Genet.
2003
;
34
:
85
-90.
38
Yeoh EJ, Ross ME, Shurtleff SA, et al. Classification, subtype discovery, and prediction of outcome in pediatric acute lymphoblastic leukemia by gene expression profiling.
Cancer Cell.
2002
;
1
:
133
-143.
39
Kong W, Engel K, Wang J. Mammalian nucleoside transporters.
Curr Drug Metab.
2004
;
5
:
63
-84.
40
Masson E, Synold TW, Relling MV, et al. Allopurinol inhibits de novo purine synthesis in lymphoblasts of children with acute lymphoblastic leukemia.
Leukemia.
1996
;
10
:
56
-60.
41
Westfall PH, Zaykin DV, Young SS. Multiple tests for genetic effects in association studies.
Methods Mol Biol.
2002
;
184
:
143
-168.
42
Storey JD, Tibshirani R. Statistical significance for genomewide studies.
Proc Natl Acad Sci U S A.
2003
;
100
:
9440
-9445.
43
Mangravite LM, Badagnani I, Giacomini KM. Nucleoside transporters in the disposition and targeting of nucleoside analogs in the kidney.
Eur J Pharmacol.
2003
;
479
:
269
-281.
44
Bokkerink JP, Stet EH, De Abreu RA, et al. 6-Mercaptopurine: cytotoxicity and biochemical pharmacology in human malignant T-lymphoblasts.
Biochem Pharmacol.
1993
;
45
:
1455
-1463.
45
Krynetski EY, Krynetskaia NF, Yanishevski Y, Evans WE. Methylation of mercaptopurine, thioguanine, and their nucleotide metabolites by heterologously expressed human thiopurine S-methyltransferase.
Mol Pharmacol.
1995
;
47
:
1141
-1147.
46
Krynetski EY, Evans WE. Pharmacogenetics of cancer therapy: getting personal.
Am J Hum Genet.
1998
;
63
:
11
-16.
47
Relling MV, Hancock ML, Rivera GK, et al. Mercaptopurine therapy intolerance and heterozygosity at the thiopurine S-methyltransferase gene locus.
J Natl Cancer Inst.
1999
;
91
:
2001
-2008.
48
Krynetski EY, Schuetz JD, Galpin AJ, et al. A single point mutation leading to loss of catalytic activity in human thiopurine S-methyltransferase.
Proc Natl Acad Sci U S A.
1995
;
92
:
949
-953.
49
Zimm S, Collins JM, O'Neill D, Chabner BA, Poplack DG. Inhibition of first-pass metabolism in cancer chemotherapy: interaction of 6-mercaptopurine and allopurinol.
Clin Pharmacol Ther.
1983
;
34
:
810
-817.
50
Winick NJ, Carroll WL, Hunger SP. Childhood leukemia—new advances and challenges.
N Engl J Med.
2004
;
351
:
601
-603.
51
Griffiths M, Beaumont N, Yao SY, et al. Cloning of a human nucleoside transporter implicated in the cellular uptake of adenosine and chemotherapeutic drugs.
Nat Med.
1997
;
3
:
89
-93.
52
Huang Y, Anderle P, Bussey KJ, et al. Membrane transporters and channels: role of the transportome in cancer chemosensitivity and chemoresistance.
Cancer Res.
2004
;
64
:
4294
-4301.
53
Galmarini CM, Thomas X, Calvo F, et al. In vivo mechanisms of resistance to cytarabine in acute myeloid leukaemia.
Br J Haematol.
2002
;
117
:
860
-868.
54
Hashimoto H, Kubota M, Shimizu T, et al. Biochemical basis of the prevention of 6-thiopurine toxicity by the nucleobases, hypoxanthine and adenine.
Leuk Res.
1990
;
14
:
1061
-1066.
55
Zimm S, Collins JM, Riccardi R, et al. Variable bioavailability of oral mercaptopurine: is maintenance chemotherapy in acute lymphoblastic leukemia being optimally delivered?
N Engl J Med.
1983
;
308
:
1005
-1009.
56
Pieters R, Huismans DR, Loonen AH, et al. Adenosine deaminase and purine nucleoside phosphorylase in childhood lymphoblastic leukemia: relation with differentiation stage, in vitro drug resistance and clinical prognosis.
Leukemia.
1992
;
6
:
375
-380.
57
Agarwal RP, Crabtree GW, Parks RE Jr, et al. Purine nucleoside metabolism in the erythrocytes of patients with adenosine deaminase deficiency and severe combined immunodeficiency.
J Clin Invest.
1976
;
57
:
1025
-1035.
58
Lennard L, Hale JP, Lilleyman JS. Red blood cell hypoxanthine phosphoribosyltransferase activity measured using 6-mercaptopurine as a substrate: a population study in children with acute lymphoblastic leukaemia.
Br J Clin Pharmacol.
1993
;
36
:
277
-284.
59
Jabs EW, Thomas PJ, Bernstein M, et al. Chromosomal localization of genes required for the terminal steps of oxidative metabolism: alpha and gamma subunits of ATP synthase and the phosphate carrier.
Hum Genet.
1994
;
93
:
600
-602.
60
Chen H, Morris MA, Rossier C, Blouin JL, Antonarakis SE. Cloning of the cDNA for the human ATP synthase OSCP subunit (ATP5O) by exon trapping and mapping to chromosome 21q22.1-q22.2.
Genomics.
1995
;
28
:
470
-476.
61
Nakamura J, Lou L. Biochemical characterization of human GMP synthetase.
J Biol Chem.
1995
;
270
:
7347
-7353.
62
Liliemark J, Pettersson B, Engberg B, et al. On the paradoxically concentration-dependent metabolism of 6-mercaptopurine in WEHI-3b murine leukemia cells.
Cancer Res.
1990
;
50
:
108
-112.
63
Bokkerink JP, Bakker MA, Hulscher TW, De Abreu RA, Schretlen ED. Purine de novo synthesis as the basis of synergism of methotrexate and 6-mercaptopurine in human malignant lymphoblasts of different lineages.
Biochem Pharmacol.
1988
;
37
:
2321
-2327.
Sign in via your Institution