Pharmacogenetics of outcome in children with acute lymphoblastic leukemia
http://www.100md.com
《血液学杂志》
the Hematologic Malignancies Program, Department of Pharmaceutical Sciences, the Department of Hematology-Oncology
the Department Biostatistics, Saint Jude Children's Research Hospital of the University of Tennessee, Memphis
the Department of Psychiatry, Pediatrics
Human Genetics, University of Chicago, Chicago, IL.
Abstract
Acquired genetic characteristics of acute lymphoblastic leukemia (ALL) cells are used to individualize therapy, whereas germ line genetic characteristics generally are not. We determined whether ALL outcome was related to 16 genetic polymorphisms affecting the pharmacodynamics of antileukemic agents. Of 246 children, 116 were treated on the lower-risk (LR) and 130 on the higher-risk (HR) arms of a St Jude protocol. Patients in the HR group with the glutathione S-transferase (GSTM1) nonnull genotype had greater risk of hematologic relapse (P = .03), which was further increased by the thymidylate synthetase (TYMS) 3/3 genotype (P = .03). These genotypes remained predictive in multivariate analyses (P < .001 and .003, respectively). No genotypes were predictive in the LR arm. expression of these 2 genes in ALL blasts was lower in those with low-activity genotypes. For central nervous system relapse, among the HR group, the vitamin D receptor start site (P = .02) and intron 8 genotypes (P = .04) predisposed, whereas for LR patients the TYMS 3/3 genotype predisposed (P = .04). The GSTM1 non-null and TYMS 3/3 genotypes are plausibly linked to drug resistance. Polymorphisms interact to influence antileukemic outcome and represent determinants of response that can be used to optimize therapy. (Blood. 2005;105:4752-4758)
Introduction
Childhood acute lymphoblastic leukemia (ALL) is cured in approximately 80% of patients.1-4 One treatment strategy for improving cure has been to modify the intensity of therapy based on the acquired genetic characteristics of the leukemia, using more intensive therapy for ALL with molecular markers (eg, the t(9;22) or MLL translocations), indicating a resistant leukemia.5 Other than age, host characteristics are usually not used for assigning patients with ALL to risk groups or to determine the intensity of therapy.
ALL regimens consist of combination chemotherapy (eg, vincristine, glucocorticoid, methotrexate, a thiopurine, and asparaginase with or without anthracyclines, topoisomerase II inhibitors, cytarabine, and cyclophosphamide) administered continuously over 2 to 3 years.1-4,6 Outcome of ALL may be influenced by modest changes in drug dose or exposure.7-10 Thus, if the determinants of interpatient variability in drug pharmacodynamics were better defined, tailoring drug therapy based on these factors might further improve outcome.
Germ line polymorphisms in genes that code for the proteins involved in the pharmacodynamics of antileukemic agents are common, with the frequency of the "variant" allele ranging from 5% to 50%.11,12 Although we and others have explored whether germ line polymorphisms may relate to ALL outcomes,13-19 such studies have examined only a few polymorphisms at a time, have failed to adjust for other risk factors, have not distinguished among relapse risk versus other types of adverse outcome, or have potential for selection bias by including only a small fraction of patients enrolled on the treatment protocols. Most prior studies have not assessed whether gene-gene interactions could influence the relationships between germ line polymorphisms and outcomes.
Herein, we have used a candidate gene approach to assess whether ALL outcome is related to 16 common polymorphisms in genes plausibly linked to the pharmacodynamics of the drug therapy, accounting for known prognostic factors, types of ALL outcome, competing risks, and gene-gene interactions.
Patients, materials, and methods
Study population
Of the 247 children with newly diagnosed ALL treated on the St Jude Children's Research Hospital Total XIIIB study for ALL,20 246 were included; one patient died within a few weeks after diagnosis and a germ line DNA sample could not be obtained. The research was approved by our Institutional Review Board, and informed consent was obtained from parents, guardians, or patients (as appropriate). Patients were prospectively assigned to risk groups that determined the intensity of therapy; those with at least one of the following were assigned to the higher-risk group: age younger than 1 year or older than 10 years (except cases with DNA index 1.16 and 1.60), initial leukocyte count greater than 50 x 109/L (except cases with DNA index 1.16 and 1.60), central nervous system (CNS) involvement (defined as 5 white cells/mm3 in cerebrospinal fluid and the presence of leukemic blasts on cytocentrifuge preparation from a nontraumatic spinal tap, or cranial palsies), testicular involvement, greater than 5% leukemic blasts in bone marrow on day 19, T-cell ALL, presence of BCR-ABL [t(9;22)], MLL-AF4 [t(4;11)], e2A-PBX1 [t(1;19)] in pre-B ALL, or MLL rearrangement by karyotype and/or molecular probes. All others were enrolled on the lower-risk arm.21 Minimal residual disease on day 43 was assessed by flow cytometric assays as previously described.21,22
Race/ethnicity was classified by a hospital staff member at the time of diagnosis (using 1994-1998 policies). Parents (or patients if older than 18 years of age) were asked which of the race categories (American Indian, black, Chinese, Filipino, Hispanic, Japanese, white, or unknown) that they considered applied to their child. Race and/or ethnicity was assessed in this study, because of the prior association of race with ALL outcome in some studies,23,24 and because of the considerable racial differences in frequencies of genotypes.11,12,25
In summary, treatment consisted of high-dose methotrexate, mercaptopurine, or both, followed by therapy with prednisone, vincristine, daunorubicin, asparaginase, cytarabine, and etoposide.20 At complete remission, patients received 2 weeks of high-dose methotrexate and mercaptopurine, followed by 120 weeks of risk-directed postremission rotating pairs of chemotherapy. For lower-risk cases, postremission therapy consisted of daily mercaptopurine and weekly methotrexate, high-dose methotrexate and mercaptopurine every 8 weeks for the first year, and dexamethasone plus vincristine every 4 weeks. For higher-risk cases, postremission therapy consisted of drug pairs (etoposide, cyclophosphamide, mercaptopurine, methotrexate, cytarabine, dexamethasone, vincristine, and mercaptopurine) administered in weekly rotation, interrupted by reinduction therapy from weeks 16 to 21.21 events were defined as any relapse, death, or secondary malignancies, using criteria previously published.26
Genotyping
DNA was extracted from normal blood cells. Candidate genes exhibited polymorphisms and encoded proteins that were involved in the pharmacokinetics or pharmacodynamics of antileukemic agents (Figure 1). Priority was given to polymorphisms that were demonstrated to be associated with phenotypes in both clinical and preclinical studies.9,11-14,27-33 Vincristine, prednisone, dexamethasone, etoposide, and daunorubicin are substrates for cytochrome P4503A4 (CYP3A4), cytochrome P4503A5 (CYP3A5), and/or p-glycoprotein (MDR1), which are partially regulated by the vitamin D receptor (VDR).27,34,35 The parent drug or metabolites of etoposide, daunorubicin, and cyclophosphamide are substrates for glutathione S-transferases (GSTM1, GSTP1, and GSTT1).28 etoposide is glucuronidated by UDP-glucuronosyltransferase 1A1 (UGT1A1).36 Mercaptopurine is a substrate for thiopurine S-methyltransferase (TPMT),9 and methotrexate interacts with methylenetetrahydrofolate reductase (MTHFR), thymidylate synthetase (TYMS), and the reduced folate carrier (RFC, or SLC19A1).11 Prednisone and dexamethasone bind to the glucocorticoid receptor (NR3C1).29
Genotyping was performed for the 16 polymorphic loci using methods as previously described: CYP3A41B (A>G at position -392) and CYP3A53 (G>A at position 22893)37; GSTP1 313A>G,37 GSTM1 deletion and GSTT1 deletion16; MDR1 exon 21 (2677G>T/A) and MDR1 exon 26 (3435C>T)37; MTHFR 677C>T38 and MTHFR 1298A>C30; NR3C1 1088A>G31; RFC 80G>A38; TPMT 238G>C, TPMT 460G>A, and TPMT 719A>G32; TYMS enhancer repeat33; UGT1A1 promoter repeat polymorphism37; VDR intron 8 G>A, and VDR FokI (start-site) T>C.37
Gene expression
In all children with available RNA in their diagnostic ALL blasts, gene expression levels were analyzed using the Affymetrix GeneChip array HG_U95Av2, as described.39 expression signals were scaled to a target average intensity and analyzed using Affymetrix Microarray Suite (MAS) version 5.0 (Affymetrix, Santa Clara, CA). The expression levels of GSTM1 and TYMS were based on probe sets 39054_at and 37899_at, respectively. Patients whose blasts had numeric or structural abnormalities of chromosomes 1 and 18 were excluded from the analysis of blast expression of the GSTM1 and TYMS genes, located on 1p13.3 and 18p11.32, respectively.
Statistical analysis
The distribution of genotypes between racial groups or between risk groups was compared using chi-square test with Yates correction. The cumulative incidences of hematologic and CNS relapse by risk group were compared using the Gray estimator with incorporation of competing events.40 Death in remission and second malignancy were treated as competing risks. Survival time was the time between diagnosis and the date of last follow-up.
To assess possible pharmacogenetic determinants of relapse, a classification and regression tree (CART) approach was used.41 Statistical analysis was performed on each polymorphism individually to make a decision tree for prediction of relapses. For each outcome, the most significant polymorphism determined the second node of the CART (after treatment arm), and polymorphisms for the remaining 15 loci were then analyzed, and the process was repeated. At each node, all the loci that had not yet been assessed for association between genotypes with outcome were analyzed, and the polymorphism with the most significant association (by the smallest P value using Gray test) was used to define that node and the downward branches. The CARTs were stopped either when no further significant genotypes were observed or when the number of patients in the terminal nodes was equal to 10. In these univariate analyses, Gray test was used to compare the cumulative incidences of relapse by genotypic groups.40
To assess whether genetic variants had independent prognostic significance for outcome, we performed multiple regression analyses adjusting for the established risk factors in this protocol21: initial leukocyte count (< 100 x 109/L versus 100 x 109/L), DNA index (1.16-1.6 versus others), presence versus absence of BCR-ABL [t(9;22)] or of MLL-AF4 [t(4;11)], percentage of blasts in the bone marrow on day 19 (< 5% blasts versus 5% blasts), and minimal residual disease on remission date (< 0.01% versus 0.01%). For the analysis of CNS relapse, CNS status at diagnosis was also included as a possible prognostic factor (no blasts in cerebrospinal fluid in an atraumatic or traumatic lumbar puncture versus other CNS status). Fine and Gray estimator with incorporation of competing events was used to compare the cumulative incidence of relapse by genotypic groups.42
To assess the possible correlation between genotypes at the 16 polymorphic loci, pair-wise associations were performed using Fisher exact or chi-square tests. We used Wilcoxon rank-sum test to compare the expression levels of GSTM1 and TYMS among genotypic groups.
Results
Among those on the higher-risk arm, the nonnull GSTM1 genotype (Figure 3A) was associated with an increased risk of hematologic relapse (isolated or combined) (5-year cumulative incidence, 17.1% ± 4.5% versus 5.1% ± 2.9% for GSTM1 null genotype, P = .03). Among those with the GSTM1 nonnull genotype, the presence of the TYMS 3/3 genotype was associated with a higher-risk of hematologic relapse (5-year cumulative incidence, 29.2% ± 9.5% versus 10.9% ± 4.7% for TYMS 2/3 or 2/2 genotypes, P = .02) (Figure 3B), whereas among those with GSTM1 null genotype, the NR3C1 AG genotype was associated with an increased risk of hematologic relapse (5-year cumulative incidence, 33.3% ± 33.3% versus 3.6% ± 2.5% for NR3C1 AA genotype, P = .02). In univariate analyses, the only other significant predictors of hematologic relapse were the presence of BCR-ABL [t(9;22)] at diagnosis (P < .001), the percentage of blasts in the day 19 bone marrow (P < .001), and minimal residual disease status on remission date (P < .001). In multivariate analyses, including other features predictive of hematologic relapse, either the GSTM1 nonnull genotype alone (P < .001) or the combined GSTM1 nonnull and TYMS 3/3 genotypes (P = .003), were independent prognostic factors for hematologic relapse (Tables 2 and 3). No genotype was associated with hematologic relapse among those patients on the lower-risk arm.
To explore whether genotypes were related to expression of those genes in the target tissue, we compared gene expression levels in the same set of patients in whom GSTM1 and TYMS genotypes were prognostic. In diagnostic ALL blasts from children treated with the higher-risk therapy, those with the GSTM1 nonnull genotype had higher GSTM1 gene expression than those with the null genotype (P < .001), and among those with the GSTM1 nonnull genotype, there was a trend for higher TYMS expression in those with the TYMS 3/3 compared with 2/3 or 2/2 genotypes (P = .09) (Figure 4).
Because different polymorphisms might be important for drug penetration into the CNS, we performed a separate analysis of genetic determinants of CNS relapse (isolated or combined), although only 6 and 2 CNS relapses occurred on the higher- and lower-risk arms, respectively. On the higher-risk arm, the VDR FokI T allele was associated with an increased risk of CNS relapse (5-year cumulative incidence, 8.6% ± 3.4% versus 0% for VDR FokI CC genotype, P = .02). All CNS relapses occurred among those with at least 1 VDR FokI T allele, but the frequency of VDR FokI genotypes did not differ (P > .14 for all) between patients with high- versus low-risk features (by initial leukocyte count, CNS status at diagnosis, DNA index, BCR-ABL [t(9;22)], MLL-AF4 [t(4;11)], day 19 bone marrow blasts, and minimal residual disease). Among those with at least 1 VDR FokI T allele, the VDR intron 8 GG genotype was associated with a greater risk of CNS relapse (5-year cumulative incidence, 15.6% ± 6.6% versus 2.7% ± 2.7% for VDR intron 8 AA genotype, P = .04), and the combined VDR FokI and VDR intron 8 genotypes were significantly associated with CNS relapse (P = .01) in a multivariate analysis including other prognostic features. On the lower-risk arm, the TYMS 3/3 genotype was associated with an increased risk of CNS relapse (5-year cumulative incidence, 5.7% ± 4% versus 0% for TYMS 2/3 or 2/2 genotypes, P = .04).
We compared the frequencies of significantly predictive genotypes in the higher- and lower-risk arms: no predictive polymorphisms differed in frequency between risk groups. Within risk groups, the only association of predictive genotypes with other patient features was among the higher-risk group: GSTM1 null genotype was more common (P = .03) among those with initial leukocyte count at least 100 x 109/L (57.9%) than those with less than 100 x 109/L (36.8%), but GSTM1 remained a significant predictor of hematologic relapse even after adjustment for presenting leukocyte counts and other prognostic features (Table 2). To assess the possibility that the association of genotypes with outcome was confounded by associations among genotypes, we assessed the correlations among genotypes at the 16 polymorphic loci. Although there was the expected linkage disequilibrium within some loci (eg, CYP3A4 with CYP3A5, MDR1 exon 21 with MDR1 exon 26 polymorphisms), there was no significant association between any of the predictive genotypes (eg, GSTM1 was not associated with TYMS or NR3C1 genotypes) (data not shown).
Discussion
The most common reason for failure in childhood ALL is hematologic relapse. We previously found that the only prognostic features in this protocol were risk group (treatment arm), initial leukocyte count, DNA index, BCR-ABL [t(9;22)], MLL-AF4 [t(4;11)], and early treatment response.21 Herein, after accounting for the difference in relapse hazard because of treatment arms, risk group, and other prognostic features (Tables 2 and 3), the incidence of hematologic relapse was significantly associated with common germ line genetic polymorphisms. In the higher-risk arm, the GSTM1 deletion was a determinant for the risk of hematologic relapse: 13 of the 16 children who experienced a hematologic relapse had the nonnull genotype. Our analysis allowed us to identify gene-gene interactions that modified the relapse risk. For patients with the GSTM1 nonnull genotype, the presence of the higher activity TYMS 3/3 genotype was associated with greater risk of hematologic relapse (P = .03). That GSTM1 and TYMS genotypes remained prognostic even after adjustment for factors reflective of early treatment response suggests that these genotypes continue to have an impact on prognosis for treatment administered beyond the remission induction period. This finding is plausible, in that drugs that are substrates for glutathione S-transferases (GSTs; eg, cyclophosphamide, etoposide), or that target thymidylate synthetase (TYMS; eg, methotrexate) were important components of the therapy given after induction. Although confidence intervals are large, the hazard ratios associated with unfavorable GSTM1 or combined GSTM1/TYMS genotypes are on par with those associated with features such as MLL translocations or the presence of minimal residual disease (Tables 2 and 3). Among those with the GSTM1 null genotype, the NR3C1 G allele was associated with an increased relapse risk, although a very small number of patients carried the G allele.
No polymorphisms were associated with hematologic relapse among patients treated on the lower-risk arm. That the GSTM1 germ line polymorphism was important for relapse among patients in the higher-risk arm, but not in the lower-risk arm, is consistent with more intensive therapy given to the former patients.20 GSTs catalyze the inactivation of many antileukemic agents and their metabolites (eg, cyclophosphamide, anthracyclines, etoposide, and steroids)28 and protect against oxidative stress. Null GST genotypes have been linked with improved outcome in some ALL trials14,16 but not in others.17 Our study is distinguished by the fact that virtually all patients on the trial were included in this analysis. Because cyclophosphamide and etoposide constituted a major proportion of therapy in the higher-risk patients, it is plausible that GST genotype was of greater importance in the higher-risk than in the lower-risk arm. Using pretreatment ALL blasts, we confirmed that the germ line nonnull GSTM1 genotype is associated with higher GSTM1 expression, in support of the hypothesis that the inferior outcome associated with the nonnull genotype is related to the increased detoxification conferred by GSTM1 expression.
The TYMS enhancer repeat 3/3 genotype causes higher expression and activity of one of the major targets of methotrexate, thymidylate synthetase.11,33 Higher thymidylate synthetase requires higher concentrations of methotrexate and its metabolites for inhibition33 and cytotoxicity.44,45 Although one single-gene study reported shorter event-free survival in children with ALL with the 3/3 genotype,18 another case-control study of 80 patients failed to find such an association.15 Such disparate findings may be due to differences in intensity of methotrexate therapy between studies. In our study, among patients on the higher-risk arm with an unfavorable GSTM1 status, the TYMS expression tended to be higher (Figure 4), and the TYMS 3/3 genotype was associated with a higher risk of hematologic relapse. Our findings illustrate the potential for gene-gene interactions: in the patients with low glutathione transferase activity, it is possible that methotrexate was a less critical element of therapy than among those with high glutathione transferase; thus, a polymorphism in the target of methotrexate (thymidylate synthetase) was of less importance in those with a favorable (low activity) than an unfavorable (high activity) GSTM1 genotype.
Only 2 CNS relapses occurred in the lower-risk group, and both patients were homozygous for the high-activity TYMS 3/3 genotype. For those on the lower-risk arm, methotrexate (given as systemic high doses and as intrathecal injections) was the primary form of CNS prophylaxis; thus, it is plausible that this unfavorable genotype might constitute a therapeutic disadvantage for this subset of patients.
In the higher-risk arm, 2 polymorphisms in the VDR locus were associated with risk of CNS relapse, which remained significant predictors after adjusting for other prognostic features. The VDR regulates the expression of CYP3A4 and possibly p-glycoprotein,34,35 which are important in the disposition of vincristine, etoposide, daunorubicin, prednisone, and dexamethasone. Both VDR polymorphisms have been associated with many clinical phenotypes.46-48 If the VDR FokI T allele and intron 8 GG genotype were associated with higher expression of p-glycoprotein, their association with an increased risk of relapse might be consistent with worse penetration of active drugs into the CNS (due to higher p-glycoprotein in the blood-brain barrier) and with higher risk of relapse because of greater p-glycoprotein-mediated biliary excretion of antileukemic drugs. Whether polymorphisms in VDR directly affect drug sensitivity or expression of downstream targets (CYP3A4 and/or p-glycoprotein) is not clear.
Relationships between genotypes and outcome could be influenced by a causative relationship between the polymorphism and response to therapy, or by an a priori relationship between the polymorphism and the development of specific subtypes of ALL. Our data favor the former rather than the latter. First, no polymorphisms differed in allele frequency between risk groups. Second, among the higher-risk group, although the GSTM1 genotype was associated with a high presenting leukocyte count, GSTM1 remained a significant predictor of relapse after adjusting for this feature.
This single-institution protocol was an ideal platform in which to assess how genotypes predicted outcomes in a racially diverse patient group, because, although genotypes differed by race, this is one of the first studies in which ALL cure rates did not differ by race.21;43 Thus, the genotype/outcome relationships we identified maintained prognostic significance in a therapeutic setting in which race itself was unimportant. Moreover, the same genotypes (GSTM1 nonnull and TYMS 3/3) were predictive for hematologic relapse within the major subgroups of whites and blacks (data not shown). Dosing based on pharmacogenetics holds the promise for delivery of "color-blind" therapy: genotyping allows for individualized dosing according to genetic rather than racial characteristics.49
Improved outcomes in patients with ALL depend on interactions among drugs, ALL blast sensitivity, and host factors. Currently, patients whose ALL blasts have unfavorable acquired genetic characteristics receive more intensive therapy. In our study of 16 target polymorphic loci, several common polymorphisms were found to be prognostic, most of which have pharmacologic plausibility. However, we highlight that prognostic genotypes are highly treatment dependent; hence, caution should be used in attempting to generalize these results to other protocols. Both of the primary polymorphisms associated with relapse, the high activity TYMS 3/3 and the GSTM1 nonnull genotypes, would be expected to cause drug resistance and were linked with higher gene expression in the target tissue of ALL blasts. We acknowledge that it is likely that additional pharmacogenetic prognostic factors will be identified in other protocols and that genotypes that were not prognostic in this study might be important in other protocols: as with all prognostic features, they depend on treatment. In the future, patients may be classified not only on the genetic characteristics of their blasts but also on host genetic characteristics, so that therapy could be intensified according to a pharmacogenetic index of drug resistance. We conclude that pharmacogenetics influences the outcome of antileukemic therapy and may provide a tool to improve individualization of therapy.
Acknowledgements
We thank our protocol coinvestigators, clinical staff, research nurses (Sheri Ring, Lisa Walters, Margaret edwards, Terri Kuehner, and Paula Condy), the patients and their parents for their participation. We also thank Jean Cai, Pam McGill, and Nancy Duran for laboratory assistance and Wenjian Yang, Nancy Kornegay, Carl Panetta, and Mark Wilkinson for computing assistance.
All the authors contributed substantially to the manuscript. M.V.R. and W.e.e. designed the research; J.C.C.R., S.K., S.D., e.H.C., J.T.S., J.R., R.R., D.C., C.-H.P., W.e.e., and M.V.R. performed the research; J.C.R., C.C., and W.L. analyzed the data; and J.C.R. and M.V.R. were primary authors of the paper.
Footnotes
Prepublished online as Blood First edition Paper, February 15, 2005; DOI 10.1182/blood-2004-11-4544.
Supported by the National Cancer Institute (CA 51001, CA 78224, CA 60419, CA21 765) and the National Institute of General Medical Sciences (NIGMS) Pharmacogenetics Research Network and Database (U01 1GM61374) from the National Institutes of Health, by a Center of excellence grant from the State of Tennessee, and by American Lebanese Syrian Associated Charities (ALSAC).
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.
References
Pui CH, Relling MV, Downing JR. Acute lymphoblastic leukemia. N engl J Med. 2004;350: 1535-1548.
Nachman JB, Sather HN, Sensel MG, et al. Augmented post-induction therapy for children with high-risk acute lymphoblastic leukemia and a slow response to initial therapy. N engl J Med. 1998;338: 1663-1671.
Schrappe M, Reiter A, Ludwig WD, et al. Improved outcome in childhood acute lymphoblastic leukemia despite reduced use of anthracyclines and cranial radiotherapy: results of trial ALL-BFM 90 German-Austrian-Swiss ALL-BFM Study Group. Blood. 2000;95: 3310-3322.
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.
Pui CH, Gaynon PS, Boyett JM, et al. Outcome of treatment in childhood acute lymphoblastic leukaemia with rearrangements of the 11q23 chromosomal region. Lancet. 2002;359: 1909-1915.
Pui CH, evans We. Acute lymphoblastic leukemia. N engl J Med. 1998;339: 605-615.
evans We, Relling MV, Rodman JH, et al. Conventional compared with individualized chemotherapy for childhood acute lymphoblastic leukemia. N engl J Med. 1998;338: 499-505.
evans We, Crom WR, Abromowitch M, et al. Clinical pharmacodynamics of high-dose methotrexate in acute lymphocytic leukemia. N engl J Med. 1986;314: 471-477.
Relling MV, Hancock ML, Boyett JM, Pui C-H, evans We. Prognostic importance of 6-mercaptopurine dose intensity in acute lymphoblastic leukemia. Blood. 1999;93: 2817-2823.
Camitta B, Mahoney D, Leventhal B, et al. Intensive intravenous methotrexate and mercaptopurine treatment of higher-risk non-T, non-B acute lymphocytic leukemia: a Pediatric Oncology Group study. J Clin Oncol. 1994;12: 1383-1389.
evans We, McLeod HL. Pharmacogenomics—drug disposition, drug targets, and side effects. N engl J Med. 2003;348: 538-549.
Weinshilboum R. Inheritance and drug response. N engl J Med. 2003;348: 529-537.
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.
Stanulla M, Schrappe M, Brechlin AM, Zimmermann M, Welte K. Polymorphisms within glutathione S-transferase genes (GSTM1, GSTT1, GSTP1) and risk of relapse in childhood B-cell precursor acute lymphoblastic leukemia: a case-control study. Blood. 2000;95: 1222-1228.
Lauten M, Asgedom G, Welte K, Schrappe M, Stanulla M. Thymidylate synthase gene polymorphism and its association with relapse in childhood B-cell precursor acute lymphoblastic leukemia. Haematologica. 2003;88: 353-354.
Chen C-L, Liu Q, Pui C-H, et al. Higher frequency of glutathione S-transferase deletions in black children with acute lymphoblastic leukemia. Blood. 1997;89: 1701-1707.
Davies SM, Bhatia S, Ross JA. et al. Glutathione S-transferase genotypes, genetic susceptibility, and outcome of therapy in childhood acute lymphoblastic leukemia. Blood. 2002;100: 67-71.
Krajinovic M, Costea I, Chiasson S. Polymorphism of the thymidylate synthase gene and outcome of acute lymphoblastic leukaemia. Lancet. 2002;359: 1033-1034.
Krajinovic M, Labuda D, Mathonnet G, et al. Polymorphisms in genes encoding drugs and xenobiotic metabolizing enzymes, DNA repair enzymes, and response to treatment of childhood acute lymphoblastic leukemia. Clin Cancer Res. 2002;8: 802-810.
Relling MV, Boyett JM, Blanco JG, et al. Granulocyte-colony stimulating factor and the risk of secondary myeloid malignancy after etoposide treatment. Blood. 2003;101: 3862-3867.
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.
Coustan-Smith e, Sancho J, Hancock ML, et al. Clinical importance of minimal residual disease in childhood acute lymphoblastic leukemia. Blood. 2000;96: 2691-2696.
Kadan-Lottick NS, Ness KK, Bhatia S, Gurney JG. Survival variability by race and ethnicity in childhood acute lymphoblastic leukemia. JAMA. 2003;290: 2008-2014.
Pollock BH, DeBaun MR, Camitta BM, et al. Racial differences in the survival of childhood B-precursor acute lymphoblastic leukemia: a Pediatric Oncology Group Study. J Clin Oncol. 2000;18: 813-823.
Blanco JG, edick MJ, Hancock ML, et al. Genetic polymorphisms in CYP3A5, CYP3A4 and NQO1 in children who developed therapy-related myeloid malignancies. Pharmacogenetics. 2002;12: 605-611.
Pui CH, Cheng C, Leung W, et al. extended follow-up of long-term survivors of childhood acute lymphoblastic leukemia. N engl J Med. 2003;349: 640-649.
Pastan I, Gottesman MM. Multidrug resistance. Annu Rev Med. 1991;42: 277-286.
Townsend DM, Tew KD. The role of glutathione-S-transferase in anti-cancer drug resistance. Oncogene. 2003;22: 7369-7375.
Tissing WJ, Meijerink JP, den Boer ML, Pieters R. Molecular determinants of glucocorticoid sensitivity and resistance in acute lymphoblastic leukemia. Leukemia. 2003;17: 17-25.
Weisberg IS, Jacques PF, Selhub J, et al. The 1298AC polymorphism in methylenetetrahydrofolate reductase (MTHFR): in vitro expression and association with homocysteine. Atherosclerosis. 2001;156: 409-415.
Huizenga NA, Koper JW, De Lange P, et al. A polymorphism in the glucocorticoid receptor gene may be associated with and increased sensitivity to glucocorticoids in vivo. J Clin endocrinol Metab. 1998;83: 144-151.
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.
Wang W, Marsh S, Cassidy J, McLeod HL. Pharmacogenomic dissection of resistance to thymidylate synthase inhibitors. Cancer Res. 2001;61: 5505-5510.
Schmiedlin-Ren P, Thummel Ke, Fisher JM, et al. expression of enzymatically active CYP3A4 by Caco-2 cells grown on extracellular matrix-coated permeable supports in the presence of 1alpha,25-dihydroxyvitamin D3. Mol Pharmacol. 1997;51: 741-754.
Thummel Ke, Brimer C, Yasuda K, et al. Transcriptional control of intestinal cytochrome P-4503A by 1alpha,25-dihydroxy vitamin D3. Mol Pharmacol. 2001;60: 1399-1406.
Watanabe Y, Nakajima M, Ohashi N, Kume T, Yokoi T. Glucuronidation of etoposide in human liver microsomes is specifically catalyzed by UDP-glucuronosyltransferase 1A1. Drug Metab Dispos. 2003;31: 589-595.
Kishi S, Yang W, Boureau B, et al. effects of prednisone and genetic polymorphisms on etoposide disposition in children with acute lymphoblastic leukemia. Blood. 2004;103: 67-72.
Kishi S, Griener JC, Cheng C, et al. Homocysteine, pharmacogenetics, and neurotoxicity in children with leukemia. J Clin Oncol. 2003;21: 3084-3091.
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.
Gray RJ. A class of K-sample tests for comparing the cumulative incidence of a competing risk. Ann Statistics. 1988;16: 1141-1154.
Breiman L, Friedman JH, Olshen RA, Stone CJ. Classification and Regression Trees. Wadsworth and Brooks, Pacific Grove, CA 1984.
Fine JP, Gray RJ. A proportional hazards model for the subdistribution of a competing risk. J A Stat Assoc. 1999;94: 496-509.
Pui CH, Sandlund JT, Pei D, et al. Results of therapy for acute lymphoblastic leukemia in black and white children. JAMA. 2003;290: 2001-2007.
Johnston PG, Mick R, Recant W, et al. Thymidylate synthase expression and response to neoadjuvant chemotherapy in patients with advanced head and neck cancer. J Natl Cancer Inst. 1997;89: 308-313.
Villafranca e, Okruzhnov Y, Dominguez MA, et al. Polymorphisms of the repeated sequences in the enhancer region of the thymidylate synthase gene promoter may predict downstaging after preoperative chemoradiation in rectal cancer. J Clin Oncol. 2001;19: 1779-1786.
Morrison NA, Qi JC, Tokita A, et al. Prediction of bone density from vitamin D receptor alleles. Nature. 1994;367: 284-287.
eisman JA. Pharmacogenetics of the vitamin D receptor and osteoporosis. Drug Metab Dispos. 2001;29: 505-512.
Sainz J, Van Tornout JM, Loro ML, et al. Vitamin D-receptor gene polymorphisms and bone density in prepubertal American girls of Mexican descent. N engl J Med. 1997;337: 77-82.
Phillips KA, Veenstra DL, Oren e, Lee JK, Sadee W. Potential role of pharmacogenomics in reducing adverse drug reactions: a systematic review. JAMA. 2001;286: 2270-2279.(Jose Claudio C. Rocha, Ch)
the Department Biostatistics, Saint Jude Children's Research Hospital of the University of Tennessee, Memphis
the Department of Psychiatry, Pediatrics
Human Genetics, University of Chicago, Chicago, IL.
Abstract
Acquired genetic characteristics of acute lymphoblastic leukemia (ALL) cells are used to individualize therapy, whereas germ line genetic characteristics generally are not. We determined whether ALL outcome was related to 16 genetic polymorphisms affecting the pharmacodynamics of antileukemic agents. Of 246 children, 116 were treated on the lower-risk (LR) and 130 on the higher-risk (HR) arms of a St Jude protocol. Patients in the HR group with the glutathione S-transferase (GSTM1) nonnull genotype had greater risk of hematologic relapse (P = .03), which was further increased by the thymidylate synthetase (TYMS) 3/3 genotype (P = .03). These genotypes remained predictive in multivariate analyses (P < .001 and .003, respectively). No genotypes were predictive in the LR arm. expression of these 2 genes in ALL blasts was lower in those with low-activity genotypes. For central nervous system relapse, among the HR group, the vitamin D receptor start site (P = .02) and intron 8 genotypes (P = .04) predisposed, whereas for LR patients the TYMS 3/3 genotype predisposed (P = .04). The GSTM1 non-null and TYMS 3/3 genotypes are plausibly linked to drug resistance. Polymorphisms interact to influence antileukemic outcome and represent determinants of response that can be used to optimize therapy. (Blood. 2005;105:4752-4758)
Introduction
Childhood acute lymphoblastic leukemia (ALL) is cured in approximately 80% of patients.1-4 One treatment strategy for improving cure has been to modify the intensity of therapy based on the acquired genetic characteristics of the leukemia, using more intensive therapy for ALL with molecular markers (eg, the t(9;22) or MLL translocations), indicating a resistant leukemia.5 Other than age, host characteristics are usually not used for assigning patients with ALL to risk groups or to determine the intensity of therapy.
ALL regimens consist of combination chemotherapy (eg, vincristine, glucocorticoid, methotrexate, a thiopurine, and asparaginase with or without anthracyclines, topoisomerase II inhibitors, cytarabine, and cyclophosphamide) administered continuously over 2 to 3 years.1-4,6 Outcome of ALL may be influenced by modest changes in drug dose or exposure.7-10 Thus, if the determinants of interpatient variability in drug pharmacodynamics were better defined, tailoring drug therapy based on these factors might further improve outcome.
Germ line polymorphisms in genes that code for the proteins involved in the pharmacodynamics of antileukemic agents are common, with the frequency of the "variant" allele ranging from 5% to 50%.11,12 Although we and others have explored whether germ line polymorphisms may relate to ALL outcomes,13-19 such studies have examined only a few polymorphisms at a time, have failed to adjust for other risk factors, have not distinguished among relapse risk versus other types of adverse outcome, or have potential for selection bias by including only a small fraction of patients enrolled on the treatment protocols. Most prior studies have not assessed whether gene-gene interactions could influence the relationships between germ line polymorphisms and outcomes.
Herein, we have used a candidate gene approach to assess whether ALL outcome is related to 16 common polymorphisms in genes plausibly linked to the pharmacodynamics of the drug therapy, accounting for known prognostic factors, types of ALL outcome, competing risks, and gene-gene interactions.
Patients, materials, and methods
Study population
Of the 247 children with newly diagnosed ALL treated on the St Jude Children's Research Hospital Total XIIIB study for ALL,20 246 were included; one patient died within a few weeks after diagnosis and a germ line DNA sample could not be obtained. The research was approved by our Institutional Review Board, and informed consent was obtained from parents, guardians, or patients (as appropriate). Patients were prospectively assigned to risk groups that determined the intensity of therapy; those with at least one of the following were assigned to the higher-risk group: age younger than 1 year or older than 10 years (except cases with DNA index 1.16 and 1.60), initial leukocyte count greater than 50 x 109/L (except cases with DNA index 1.16 and 1.60), central nervous system (CNS) involvement (defined as 5 white cells/mm3 in cerebrospinal fluid and the presence of leukemic blasts on cytocentrifuge preparation from a nontraumatic spinal tap, or cranial palsies), testicular involvement, greater than 5% leukemic blasts in bone marrow on day 19, T-cell ALL, presence of BCR-ABL [t(9;22)], MLL-AF4 [t(4;11)], e2A-PBX1 [t(1;19)] in pre-B ALL, or MLL rearrangement by karyotype and/or molecular probes. All others were enrolled on the lower-risk arm.21 Minimal residual disease on day 43 was assessed by flow cytometric assays as previously described.21,22
Race/ethnicity was classified by a hospital staff member at the time of diagnosis (using 1994-1998 policies). Parents (or patients if older than 18 years of age) were asked which of the race categories (American Indian, black, Chinese, Filipino, Hispanic, Japanese, white, or unknown) that they considered applied to their child. Race and/or ethnicity was assessed in this study, because of the prior association of race with ALL outcome in some studies,23,24 and because of the considerable racial differences in frequencies of genotypes.11,12,25
In summary, treatment consisted of high-dose methotrexate, mercaptopurine, or both, followed by therapy with prednisone, vincristine, daunorubicin, asparaginase, cytarabine, and etoposide.20 At complete remission, patients received 2 weeks of high-dose methotrexate and mercaptopurine, followed by 120 weeks of risk-directed postremission rotating pairs of chemotherapy. For lower-risk cases, postremission therapy consisted of daily mercaptopurine and weekly methotrexate, high-dose methotrexate and mercaptopurine every 8 weeks for the first year, and dexamethasone plus vincristine every 4 weeks. For higher-risk cases, postremission therapy consisted of drug pairs (etoposide, cyclophosphamide, mercaptopurine, methotrexate, cytarabine, dexamethasone, vincristine, and mercaptopurine) administered in weekly rotation, interrupted by reinduction therapy from weeks 16 to 21.21 events were defined as any relapse, death, or secondary malignancies, using criteria previously published.26
Genotyping
DNA was extracted from normal blood cells. Candidate genes exhibited polymorphisms and encoded proteins that were involved in the pharmacokinetics or pharmacodynamics of antileukemic agents (Figure 1). Priority was given to polymorphisms that were demonstrated to be associated with phenotypes in both clinical and preclinical studies.9,11-14,27-33 Vincristine, prednisone, dexamethasone, etoposide, and daunorubicin are substrates for cytochrome P4503A4 (CYP3A4), cytochrome P4503A5 (CYP3A5), and/or p-glycoprotein (MDR1), which are partially regulated by the vitamin D receptor (VDR).27,34,35 The parent drug or metabolites of etoposide, daunorubicin, and cyclophosphamide are substrates for glutathione S-transferases (GSTM1, GSTP1, and GSTT1).28 etoposide is glucuronidated by UDP-glucuronosyltransferase 1A1 (UGT1A1).36 Mercaptopurine is a substrate for thiopurine S-methyltransferase (TPMT),9 and methotrexate interacts with methylenetetrahydrofolate reductase (MTHFR), thymidylate synthetase (TYMS), and the reduced folate carrier (RFC, or SLC19A1).11 Prednisone and dexamethasone bind to the glucocorticoid receptor (NR3C1).29
Genotyping was performed for the 16 polymorphic loci using methods as previously described: CYP3A41B (A>G at position -392) and CYP3A53 (G>A at position 22893)37; GSTP1 313A>G,37 GSTM1 deletion and GSTT1 deletion16; MDR1 exon 21 (2677G>T/A) and MDR1 exon 26 (3435C>T)37; MTHFR 677C>T38 and MTHFR 1298A>C30; NR3C1 1088A>G31; RFC 80G>A38; TPMT 238G>C, TPMT 460G>A, and TPMT 719A>G32; TYMS enhancer repeat33; UGT1A1 promoter repeat polymorphism37; VDR intron 8 G>A, and VDR FokI (start-site) T>C.37
Gene expression
In all children with available RNA in their diagnostic ALL blasts, gene expression levels were analyzed using the Affymetrix GeneChip array HG_U95Av2, as described.39 expression signals were scaled to a target average intensity and analyzed using Affymetrix Microarray Suite (MAS) version 5.0 (Affymetrix, Santa Clara, CA). The expression levels of GSTM1 and TYMS were based on probe sets 39054_at and 37899_at, respectively. Patients whose blasts had numeric or structural abnormalities of chromosomes 1 and 18 were excluded from the analysis of blast expression of the GSTM1 and TYMS genes, located on 1p13.3 and 18p11.32, respectively.
Statistical analysis
The distribution of genotypes between racial groups or between risk groups was compared using chi-square test with Yates correction. The cumulative incidences of hematologic and CNS relapse by risk group were compared using the Gray estimator with incorporation of competing events.40 Death in remission and second malignancy were treated as competing risks. Survival time was the time between diagnosis and the date of last follow-up.
To assess possible pharmacogenetic determinants of relapse, a classification and regression tree (CART) approach was used.41 Statistical analysis was performed on each polymorphism individually to make a decision tree for prediction of relapses. For each outcome, the most significant polymorphism determined the second node of the CART (after treatment arm), and polymorphisms for the remaining 15 loci were then analyzed, and the process was repeated. At each node, all the loci that had not yet been assessed for association between genotypes with outcome were analyzed, and the polymorphism with the most significant association (by the smallest P value using Gray test) was used to define that node and the downward branches. The CARTs were stopped either when no further significant genotypes were observed or when the number of patients in the terminal nodes was equal to 10. In these univariate analyses, Gray test was used to compare the cumulative incidences of relapse by genotypic groups.40
To assess whether genetic variants had independent prognostic significance for outcome, we performed multiple regression analyses adjusting for the established risk factors in this protocol21: initial leukocyte count (< 100 x 109/L versus 100 x 109/L), DNA index (1.16-1.6 versus others), presence versus absence of BCR-ABL [t(9;22)] or of MLL-AF4 [t(4;11)], percentage of blasts in the bone marrow on day 19 (< 5% blasts versus 5% blasts), and minimal residual disease on remission date (< 0.01% versus 0.01%). For the analysis of CNS relapse, CNS status at diagnosis was also included as a possible prognostic factor (no blasts in cerebrospinal fluid in an atraumatic or traumatic lumbar puncture versus other CNS status). Fine and Gray estimator with incorporation of competing events was used to compare the cumulative incidence of relapse by genotypic groups.42
To assess the possible correlation between genotypes at the 16 polymorphic loci, pair-wise associations were performed using Fisher exact or chi-square tests. We used Wilcoxon rank-sum test to compare the expression levels of GSTM1 and TYMS among genotypic groups.
Results
Among those on the higher-risk arm, the nonnull GSTM1 genotype (Figure 3A) was associated with an increased risk of hematologic relapse (isolated or combined) (5-year cumulative incidence, 17.1% ± 4.5% versus 5.1% ± 2.9% for GSTM1 null genotype, P = .03). Among those with the GSTM1 nonnull genotype, the presence of the TYMS 3/3 genotype was associated with a higher-risk of hematologic relapse (5-year cumulative incidence, 29.2% ± 9.5% versus 10.9% ± 4.7% for TYMS 2/3 or 2/2 genotypes, P = .02) (Figure 3B), whereas among those with GSTM1 null genotype, the NR3C1 AG genotype was associated with an increased risk of hematologic relapse (5-year cumulative incidence, 33.3% ± 33.3% versus 3.6% ± 2.5% for NR3C1 AA genotype, P = .02). In univariate analyses, the only other significant predictors of hematologic relapse were the presence of BCR-ABL [t(9;22)] at diagnosis (P < .001), the percentage of blasts in the day 19 bone marrow (P < .001), and minimal residual disease status on remission date (P < .001). In multivariate analyses, including other features predictive of hematologic relapse, either the GSTM1 nonnull genotype alone (P < .001) or the combined GSTM1 nonnull and TYMS 3/3 genotypes (P = .003), were independent prognostic factors for hematologic relapse (Tables 2 and 3). No genotype was associated with hematologic relapse among those patients on the lower-risk arm.
To explore whether genotypes were related to expression of those genes in the target tissue, we compared gene expression levels in the same set of patients in whom GSTM1 and TYMS genotypes were prognostic. In diagnostic ALL blasts from children treated with the higher-risk therapy, those with the GSTM1 nonnull genotype had higher GSTM1 gene expression than those with the null genotype (P < .001), and among those with the GSTM1 nonnull genotype, there was a trend for higher TYMS expression in those with the TYMS 3/3 compared with 2/3 or 2/2 genotypes (P = .09) (Figure 4).
Because different polymorphisms might be important for drug penetration into the CNS, we performed a separate analysis of genetic determinants of CNS relapse (isolated or combined), although only 6 and 2 CNS relapses occurred on the higher- and lower-risk arms, respectively. On the higher-risk arm, the VDR FokI T allele was associated with an increased risk of CNS relapse (5-year cumulative incidence, 8.6% ± 3.4% versus 0% for VDR FokI CC genotype, P = .02). All CNS relapses occurred among those with at least 1 VDR FokI T allele, but the frequency of VDR FokI genotypes did not differ (P > .14 for all) between patients with high- versus low-risk features (by initial leukocyte count, CNS status at diagnosis, DNA index, BCR-ABL [t(9;22)], MLL-AF4 [t(4;11)], day 19 bone marrow blasts, and minimal residual disease). Among those with at least 1 VDR FokI T allele, the VDR intron 8 GG genotype was associated with a greater risk of CNS relapse (5-year cumulative incidence, 15.6% ± 6.6% versus 2.7% ± 2.7% for VDR intron 8 AA genotype, P = .04), and the combined VDR FokI and VDR intron 8 genotypes were significantly associated with CNS relapse (P = .01) in a multivariate analysis including other prognostic features. On the lower-risk arm, the TYMS 3/3 genotype was associated with an increased risk of CNS relapse (5-year cumulative incidence, 5.7% ± 4% versus 0% for TYMS 2/3 or 2/2 genotypes, P = .04).
We compared the frequencies of significantly predictive genotypes in the higher- and lower-risk arms: no predictive polymorphisms differed in frequency between risk groups. Within risk groups, the only association of predictive genotypes with other patient features was among the higher-risk group: GSTM1 null genotype was more common (P = .03) among those with initial leukocyte count at least 100 x 109/L (57.9%) than those with less than 100 x 109/L (36.8%), but GSTM1 remained a significant predictor of hematologic relapse even after adjustment for presenting leukocyte counts and other prognostic features (Table 2). To assess the possibility that the association of genotypes with outcome was confounded by associations among genotypes, we assessed the correlations among genotypes at the 16 polymorphic loci. Although there was the expected linkage disequilibrium within some loci (eg, CYP3A4 with CYP3A5, MDR1 exon 21 with MDR1 exon 26 polymorphisms), there was no significant association between any of the predictive genotypes (eg, GSTM1 was not associated with TYMS or NR3C1 genotypes) (data not shown).
Discussion
The most common reason for failure in childhood ALL is hematologic relapse. We previously found that the only prognostic features in this protocol were risk group (treatment arm), initial leukocyte count, DNA index, BCR-ABL [t(9;22)], MLL-AF4 [t(4;11)], and early treatment response.21 Herein, after accounting for the difference in relapse hazard because of treatment arms, risk group, and other prognostic features (Tables 2 and 3), the incidence of hematologic relapse was significantly associated with common germ line genetic polymorphisms. In the higher-risk arm, the GSTM1 deletion was a determinant for the risk of hematologic relapse: 13 of the 16 children who experienced a hematologic relapse had the nonnull genotype. Our analysis allowed us to identify gene-gene interactions that modified the relapse risk. For patients with the GSTM1 nonnull genotype, the presence of the higher activity TYMS 3/3 genotype was associated with greater risk of hematologic relapse (P = .03). That GSTM1 and TYMS genotypes remained prognostic even after adjustment for factors reflective of early treatment response suggests that these genotypes continue to have an impact on prognosis for treatment administered beyond the remission induction period. This finding is plausible, in that drugs that are substrates for glutathione S-transferases (GSTs; eg, cyclophosphamide, etoposide), or that target thymidylate synthetase (TYMS; eg, methotrexate) were important components of the therapy given after induction. Although confidence intervals are large, the hazard ratios associated with unfavorable GSTM1 or combined GSTM1/TYMS genotypes are on par with those associated with features such as MLL translocations or the presence of minimal residual disease (Tables 2 and 3). Among those with the GSTM1 null genotype, the NR3C1 G allele was associated with an increased relapse risk, although a very small number of patients carried the G allele.
No polymorphisms were associated with hematologic relapse among patients treated on the lower-risk arm. That the GSTM1 germ line polymorphism was important for relapse among patients in the higher-risk arm, but not in the lower-risk arm, is consistent with more intensive therapy given to the former patients.20 GSTs catalyze the inactivation of many antileukemic agents and their metabolites (eg, cyclophosphamide, anthracyclines, etoposide, and steroids)28 and protect against oxidative stress. Null GST genotypes have been linked with improved outcome in some ALL trials14,16 but not in others.17 Our study is distinguished by the fact that virtually all patients on the trial were included in this analysis. Because cyclophosphamide and etoposide constituted a major proportion of therapy in the higher-risk patients, it is plausible that GST genotype was of greater importance in the higher-risk than in the lower-risk arm. Using pretreatment ALL blasts, we confirmed that the germ line nonnull GSTM1 genotype is associated with higher GSTM1 expression, in support of the hypothesis that the inferior outcome associated with the nonnull genotype is related to the increased detoxification conferred by GSTM1 expression.
The TYMS enhancer repeat 3/3 genotype causes higher expression and activity of one of the major targets of methotrexate, thymidylate synthetase.11,33 Higher thymidylate synthetase requires higher concentrations of methotrexate and its metabolites for inhibition33 and cytotoxicity.44,45 Although one single-gene study reported shorter event-free survival in children with ALL with the 3/3 genotype,18 another case-control study of 80 patients failed to find such an association.15 Such disparate findings may be due to differences in intensity of methotrexate therapy between studies. In our study, among patients on the higher-risk arm with an unfavorable GSTM1 status, the TYMS expression tended to be higher (Figure 4), and the TYMS 3/3 genotype was associated with a higher risk of hematologic relapse. Our findings illustrate the potential for gene-gene interactions: in the patients with low glutathione transferase activity, it is possible that methotrexate was a less critical element of therapy than among those with high glutathione transferase; thus, a polymorphism in the target of methotrexate (thymidylate synthetase) was of less importance in those with a favorable (low activity) than an unfavorable (high activity) GSTM1 genotype.
Only 2 CNS relapses occurred in the lower-risk group, and both patients were homozygous for the high-activity TYMS 3/3 genotype. For those on the lower-risk arm, methotrexate (given as systemic high doses and as intrathecal injections) was the primary form of CNS prophylaxis; thus, it is plausible that this unfavorable genotype might constitute a therapeutic disadvantage for this subset of patients.
In the higher-risk arm, 2 polymorphisms in the VDR locus were associated with risk of CNS relapse, which remained significant predictors after adjusting for other prognostic features. The VDR regulates the expression of CYP3A4 and possibly p-glycoprotein,34,35 which are important in the disposition of vincristine, etoposide, daunorubicin, prednisone, and dexamethasone. Both VDR polymorphisms have been associated with many clinical phenotypes.46-48 If the VDR FokI T allele and intron 8 GG genotype were associated with higher expression of p-glycoprotein, their association with an increased risk of relapse might be consistent with worse penetration of active drugs into the CNS (due to higher p-glycoprotein in the blood-brain barrier) and with higher risk of relapse because of greater p-glycoprotein-mediated biliary excretion of antileukemic drugs. Whether polymorphisms in VDR directly affect drug sensitivity or expression of downstream targets (CYP3A4 and/or p-glycoprotein) is not clear.
Relationships between genotypes and outcome could be influenced by a causative relationship between the polymorphism and response to therapy, or by an a priori relationship between the polymorphism and the development of specific subtypes of ALL. Our data favor the former rather than the latter. First, no polymorphisms differed in allele frequency between risk groups. Second, among the higher-risk group, although the GSTM1 genotype was associated with a high presenting leukocyte count, GSTM1 remained a significant predictor of relapse after adjusting for this feature.
This single-institution protocol was an ideal platform in which to assess how genotypes predicted outcomes in a racially diverse patient group, because, although genotypes differed by race, this is one of the first studies in which ALL cure rates did not differ by race.21;43 Thus, the genotype/outcome relationships we identified maintained prognostic significance in a therapeutic setting in which race itself was unimportant. Moreover, the same genotypes (GSTM1 nonnull and TYMS 3/3) were predictive for hematologic relapse within the major subgroups of whites and blacks (data not shown). Dosing based on pharmacogenetics holds the promise for delivery of "color-blind" therapy: genotyping allows for individualized dosing according to genetic rather than racial characteristics.49
Improved outcomes in patients with ALL depend on interactions among drugs, ALL blast sensitivity, and host factors. Currently, patients whose ALL blasts have unfavorable acquired genetic characteristics receive more intensive therapy. In our study of 16 target polymorphic loci, several common polymorphisms were found to be prognostic, most of which have pharmacologic plausibility. However, we highlight that prognostic genotypes are highly treatment dependent; hence, caution should be used in attempting to generalize these results to other protocols. Both of the primary polymorphisms associated with relapse, the high activity TYMS 3/3 and the GSTM1 nonnull genotypes, would be expected to cause drug resistance and were linked with higher gene expression in the target tissue of ALL blasts. We acknowledge that it is likely that additional pharmacogenetic prognostic factors will be identified in other protocols and that genotypes that were not prognostic in this study might be important in other protocols: as with all prognostic features, they depend on treatment. In the future, patients may be classified not only on the genetic characteristics of their blasts but also on host genetic characteristics, so that therapy could be intensified according to a pharmacogenetic index of drug resistance. We conclude that pharmacogenetics influences the outcome of antileukemic therapy and may provide a tool to improve individualization of therapy.
Acknowledgements
We thank our protocol coinvestigators, clinical staff, research nurses (Sheri Ring, Lisa Walters, Margaret edwards, Terri Kuehner, and Paula Condy), the patients and their parents for their participation. We also thank Jean Cai, Pam McGill, and Nancy Duran for laboratory assistance and Wenjian Yang, Nancy Kornegay, Carl Panetta, and Mark Wilkinson for computing assistance.
All the authors contributed substantially to the manuscript. M.V.R. and W.e.e. designed the research; J.C.C.R., S.K., S.D., e.H.C., J.T.S., J.R., R.R., D.C., C.-H.P., W.e.e., and M.V.R. performed the research; J.C.R., C.C., and W.L. analyzed the data; and J.C.R. and M.V.R. were primary authors of the paper.
Footnotes
Prepublished online as Blood First edition Paper, February 15, 2005; DOI 10.1182/blood-2004-11-4544.
Supported by the National Cancer Institute (CA 51001, CA 78224, CA 60419, CA21 765) and the National Institute of General Medical Sciences (NIGMS) Pharmacogenetics Research Network and Database (U01 1GM61374) from the National Institutes of Health, by a Center of excellence grant from the State of Tennessee, and by American Lebanese Syrian Associated Charities (ALSAC).
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.
References
Pui CH, Relling MV, Downing JR. Acute lymphoblastic leukemia. N engl J Med. 2004;350: 1535-1548.
Nachman JB, Sather HN, Sensel MG, et al. Augmented post-induction therapy for children with high-risk acute lymphoblastic leukemia and a slow response to initial therapy. N engl J Med. 1998;338: 1663-1671.
Schrappe M, Reiter A, Ludwig WD, et al. Improved outcome in childhood acute lymphoblastic leukemia despite reduced use of anthracyclines and cranial radiotherapy: results of trial ALL-BFM 90 German-Austrian-Swiss ALL-BFM Study Group. Blood. 2000;95: 3310-3322.
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.
Pui CH, Gaynon PS, Boyett JM, et al. Outcome of treatment in childhood acute lymphoblastic leukaemia with rearrangements of the 11q23 chromosomal region. Lancet. 2002;359: 1909-1915.
Pui CH, evans We. Acute lymphoblastic leukemia. N engl J Med. 1998;339: 605-615.
evans We, Relling MV, Rodman JH, et al. Conventional compared with individualized chemotherapy for childhood acute lymphoblastic leukemia. N engl J Med. 1998;338: 499-505.
evans We, Crom WR, Abromowitch M, et al. Clinical pharmacodynamics of high-dose methotrexate in acute lymphocytic leukemia. N engl J Med. 1986;314: 471-477.
Relling MV, Hancock ML, Boyett JM, Pui C-H, evans We. Prognostic importance of 6-mercaptopurine dose intensity in acute lymphoblastic leukemia. Blood. 1999;93: 2817-2823.
Camitta B, Mahoney D, Leventhal B, et al. Intensive intravenous methotrexate and mercaptopurine treatment of higher-risk non-T, non-B acute lymphocytic leukemia: a Pediatric Oncology Group study. J Clin Oncol. 1994;12: 1383-1389.
evans We, McLeod HL. Pharmacogenomics—drug disposition, drug targets, and side effects. N engl J Med. 2003;348: 538-549.
Weinshilboum R. Inheritance and drug response. N engl J Med. 2003;348: 529-537.
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.
Stanulla M, Schrappe M, Brechlin AM, Zimmermann M, Welte K. Polymorphisms within glutathione S-transferase genes (GSTM1, GSTT1, GSTP1) and risk of relapse in childhood B-cell precursor acute lymphoblastic leukemia: a case-control study. Blood. 2000;95: 1222-1228.
Lauten M, Asgedom G, Welte K, Schrappe M, Stanulla M. Thymidylate synthase gene polymorphism and its association with relapse in childhood B-cell precursor acute lymphoblastic leukemia. Haematologica. 2003;88: 353-354.
Chen C-L, Liu Q, Pui C-H, et al. Higher frequency of glutathione S-transferase deletions in black children with acute lymphoblastic leukemia. Blood. 1997;89: 1701-1707.
Davies SM, Bhatia S, Ross JA. et al. Glutathione S-transferase genotypes, genetic susceptibility, and outcome of therapy in childhood acute lymphoblastic leukemia. Blood. 2002;100: 67-71.
Krajinovic M, Costea I, Chiasson S. Polymorphism of the thymidylate synthase gene and outcome of acute lymphoblastic leukaemia. Lancet. 2002;359: 1033-1034.
Krajinovic M, Labuda D, Mathonnet G, et al. Polymorphisms in genes encoding drugs and xenobiotic metabolizing enzymes, DNA repair enzymes, and response to treatment of childhood acute lymphoblastic leukemia. Clin Cancer Res. 2002;8: 802-810.
Relling MV, Boyett JM, Blanco JG, et al. Granulocyte-colony stimulating factor and the risk of secondary myeloid malignancy after etoposide treatment. Blood. 2003;101: 3862-3867.
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.
Coustan-Smith e, Sancho J, Hancock ML, et al. Clinical importance of minimal residual disease in childhood acute lymphoblastic leukemia. Blood. 2000;96: 2691-2696.
Kadan-Lottick NS, Ness KK, Bhatia S, Gurney JG. Survival variability by race and ethnicity in childhood acute lymphoblastic leukemia. JAMA. 2003;290: 2008-2014.
Pollock BH, DeBaun MR, Camitta BM, et al. Racial differences in the survival of childhood B-precursor acute lymphoblastic leukemia: a Pediatric Oncology Group Study. J Clin Oncol. 2000;18: 813-823.
Blanco JG, edick MJ, Hancock ML, et al. Genetic polymorphisms in CYP3A5, CYP3A4 and NQO1 in children who developed therapy-related myeloid malignancies. Pharmacogenetics. 2002;12: 605-611.
Pui CH, Cheng C, Leung W, et al. extended follow-up of long-term survivors of childhood acute lymphoblastic leukemia. N engl J Med. 2003;349: 640-649.
Pastan I, Gottesman MM. Multidrug resistance. Annu Rev Med. 1991;42: 277-286.
Townsend DM, Tew KD. The role of glutathione-S-transferase in anti-cancer drug resistance. Oncogene. 2003;22: 7369-7375.
Tissing WJ, Meijerink JP, den Boer ML, Pieters R. Molecular determinants of glucocorticoid sensitivity and resistance in acute lymphoblastic leukemia. Leukemia. 2003;17: 17-25.
Weisberg IS, Jacques PF, Selhub J, et al. The 1298AC polymorphism in methylenetetrahydrofolate reductase (MTHFR): in vitro expression and association with homocysteine. Atherosclerosis. 2001;156: 409-415.
Huizenga NA, Koper JW, De Lange P, et al. A polymorphism in the glucocorticoid receptor gene may be associated with and increased sensitivity to glucocorticoids in vivo. J Clin endocrinol Metab. 1998;83: 144-151.
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.
Wang W, Marsh S, Cassidy J, McLeod HL. Pharmacogenomic dissection of resistance to thymidylate synthase inhibitors. Cancer Res. 2001;61: 5505-5510.
Schmiedlin-Ren P, Thummel Ke, Fisher JM, et al. expression of enzymatically active CYP3A4 by Caco-2 cells grown on extracellular matrix-coated permeable supports in the presence of 1alpha,25-dihydroxyvitamin D3. Mol Pharmacol. 1997;51: 741-754.
Thummel Ke, Brimer C, Yasuda K, et al. Transcriptional control of intestinal cytochrome P-4503A by 1alpha,25-dihydroxy vitamin D3. Mol Pharmacol. 2001;60: 1399-1406.
Watanabe Y, Nakajima M, Ohashi N, Kume T, Yokoi T. Glucuronidation of etoposide in human liver microsomes is specifically catalyzed by UDP-glucuronosyltransferase 1A1. Drug Metab Dispos. 2003;31: 589-595.
Kishi S, Yang W, Boureau B, et al. effects of prednisone and genetic polymorphisms on etoposide disposition in children with acute lymphoblastic leukemia. Blood. 2004;103: 67-72.
Kishi S, Griener JC, Cheng C, et al. Homocysteine, pharmacogenetics, and neurotoxicity in children with leukemia. J Clin Oncol. 2003;21: 3084-3091.
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.
Gray RJ. A class of K-sample tests for comparing the cumulative incidence of a competing risk. Ann Statistics. 1988;16: 1141-1154.
Breiman L, Friedman JH, Olshen RA, Stone CJ. Classification and Regression Trees. Wadsworth and Brooks, Pacific Grove, CA 1984.
Fine JP, Gray RJ. A proportional hazards model for the subdistribution of a competing risk. J A Stat Assoc. 1999;94: 496-509.
Pui CH, Sandlund JT, Pei D, et al. Results of therapy for acute lymphoblastic leukemia in black and white children. JAMA. 2003;290: 2001-2007.
Johnston PG, Mick R, Recant W, et al. Thymidylate synthase expression and response to neoadjuvant chemotherapy in patients with advanced head and neck cancer. J Natl Cancer Inst. 1997;89: 308-313.
Villafranca e, Okruzhnov Y, Dominguez MA, et al. Polymorphisms of the repeated sequences in the enhancer region of the thymidylate synthase gene promoter may predict downstaging after preoperative chemoradiation in rectal cancer. J Clin Oncol. 2001;19: 1779-1786.
Morrison NA, Qi JC, Tokita A, et al. Prediction of bone density from vitamin D receptor alleles. Nature. 1994;367: 284-287.
eisman JA. Pharmacogenetics of the vitamin D receptor and osteoporosis. Drug Metab Dispos. 2001;29: 505-512.
Sainz J, Van Tornout JM, Loro ML, et al. Vitamin D-receptor gene polymorphisms and bone density in prepubertal American girls of Mexican descent. N engl J Med. 1997;337: 77-82.
Phillips KA, Veenstra DL, Oren e, Lee JK, Sadee W. Potential role of pharmacogenomics in reducing adverse drug reactions: a systematic review. JAMA. 2001;286: 2270-2279.(Jose Claudio C. Rocha, Ch)