| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
1 School of Public Health, University of California, Berkeley, Berkeley, California and 2 Laboratory for Molecular Epidemiology, Department of Epidemiology and Biostatistics, University of California, San Francisco, San Francisco, California
Requests for reprints: Joseph L. Wiemels, Department of Epidemiology and Biostatistics, University of California, San Francisco, 500 Parnassus Avenue, 420 MU-W, San Francisco, CA 94143-0560. Phone: 415-476-3059; Fax: 415-476-6014. E-mail: joe.wiemels{at}ucsf.edu
| Abstract |
|---|
|
|
|---|
Methods: Buccal cell DNA from ALL cases (n = 294) and controls (n = 369) individually matched on gender, date of birth, Hispanic status, and maternal race were whole genome amplified and genotyped for four MDR1 SNPs, T129C (rs3213619), C1236T (rs1128503), G2677T/A (rs2032582), and C3435T (rs1045642). Detailed and time-specific information on household pesticide use was obtained using in-home interviews with the mother.
Results: Allele frequencies in non-Hispanic White and Hispanic controls were similar, and with the exception of T129C, seemed to be in strong linkage disequilibrium. Overall, the SNPs considered individually or within haplotypes (C1236T-G2677T/A-C3435T) were not significantly associated with childhood ALL. However, we observed strong evidence of a differential effect of indoor insecticide exposure (interaction odds ratio, 0.31; 95% confidence interval, 0.15-0.64; P = 0.025) on risk of ALL between carriers and noncarriers of haplotype CGC.
Conclusion: These preliminary data suggest that children carrying the haplotype CGC may be less susceptible to the leukemogenic effects of indoor insecticide exposures. Future studies are needed to confirm these findings. (Cancer Epidemiol Biomarkers Prev 2007;16(6):11727)
| Introduction |
|---|
|
|
|---|
The human multidrug resistance 1 (MDR1 or ABCB1) gene encodes P-glycoprotein, a 170-kDa membrane transport protein that extrudes a wide variety of lipophilic compounds, including chemotherapeutic agents, naturally occurring xenobiotics, pesticides, and cellular metabolites (5, 6). P-glycoprotein uses an ATP-dependent efflux transport mechanism to minimize the exposure of potentially toxic compounds to the intracellular environment and is expressed primarily in regions that act as epithelial barriers or perform excretory functions, including the gastrointestinal tract, blood-tissue barrier, liver, kidney, testis, and placenta (7). Interestingly, MDR1 is also highly expressed in several subclasses of bone marrow and peripheral leukocytes (8-11), and experimental studies have shown a decrease in natural killer and CD8+ T-cell activity after suppression of MDR1 (12, 13). Furthermore, there is evidence demonstrating the involvement of P-glycoprotein in the release of interleukin-2, interleukin-4, and IFN-
from lymphocytes (14-16). Accordingly, P-glycoprotein could serve a role in leukemia etiology based on its transport of xenobiotics or modulation of immune function.
Among the several putative environmental risk factors for childhood leukemia transported by P-glycoprotein are household pesticides (17, 18). In a study comparing normal mdr1a (+/+) mice (mdr1a is the mouse equivalent to the human MDR1 gene) to constructed mdr1a disrupted mice (/), Schinkel et al. (19) found that the absence of the gene resulted in increased toxicity by the pesticide, ivermectin, and decreased elimination of this compound. P-glycoprotein is also capable of interacting with a wide range of other structurally diverse pesticides (17). Several epidemiologic studies, including the Northern California Childhood Leukemia Study, have consistently reported an increased risk of childhood leukemia associated with indoor use of pesticides (4, 18) and have shown the importance of timing of the exposure during the pregnancy period (18, 20). Studies examining the effect of outdoor use of pesticides have led to inconsistent results (4).
To date, genetic analyses have identified over 50 single nucleotide polymorphisms (SNP) along the MDR1 gene (21-25). In one of the first systematic screens of this gene, Hoffmeyer et al. (22) described the distribution of several polymorphisms in a group of healthy Caucasian volunteers and reported a significant correlation of P-glycoprotein expression and function with a synonymous SNP located in exon 26 at position 3435 (C3435T). Subsequent studies have reported mixed results leading to debate over the true functional nature of this SNP (23, 25-28). One previous case-control study conducted in Poland reported a statistically significant increased risk of childhood ALL associated with the 3435TT genotype compared with the CC and CT genotypes combined (29). Similar associations have been found with other disease outcomes, including Parkinson's disease (30, 31), adult glioma (32), renal cell carcinoma (33), colon cancer (34), and other gastrointestinal diseases (35).
Recent studies have found C3435T to be in linkage disequilibrium with two other common SNPs, the synonymous C1236T (exon 12) and nonsynonymous triallelic G2677T/A (exon 21; refs. 23-25, 36, 37). It has been suggested that the observed functionality of C3435T may be a result of its strong linkage disequilibrium with the nonsynonymous G2677T/A (25). This evidence of linkage disequilibrium, together with observed interethnic differences in allele frequencies (37), emphasizes the need for a haplotype-based approach by race/ethnicity when assessing the functional nature of MDR1 polymorphisms. The T129C (exon 1b) SNP is located in the promoter region 7 bp downstream from the transcription initiation site and has previously been associated with P-glycoprotein expression levels (25, 38).
Although MDR1 polymorphisms, particularly T129C, C1236T, G2677T/A, and C3435T, have been widely investigated in pharmacologic studies, little is known about their potential role in the etiology of cancer, including childhood leukemia. We evaluated the role of these four MDR1 SNPs individually and within predicted haplotypes in a study of childhood ALL cases and matched controls and examined whether they modify the association between indoor insecticide exposure and the risk of childhood ALL.
| Materials and Methods |
|---|
|
|
|---|
Data Collection
A detailed in-home personal interview with the primary caretaker of each subject was conducted as soon as consent was obtained. The interview collected information on household pesticide use, including the name of the product, intended purpose, frequency of use, and time period of use (3 months before pregnancy; pregnancy; and years 1, 2, and 3). A calendar with years corresponding to the time windows of interest was used during the interview to assist the recall of time-specific information. Pesticide exposure data were censored at 1 year before the diagnosis date for cases and corresponding reference date for controls, because exact dates of pesticide use were not available to determine with certainty whether pesticide exposure occurred before or after the reference date. For this analysis, we focused on indoor insecticide use, which included professional pest control services, insect repellents, indoor flea foggers, and a variety of insect control products commonly used indoors. Buccal cell specimens collected from 95% of cases and 95% of controls at the hospital or during the in-home interview using a soft cytobrush were shipped to the study laboratory where it was used to extract a DNA sample.
MDR1 Genotyping
ALL case and control DNA samples extracted from buccal cells were whole genome amplified (40). Genotype assays done on a subset of nonwhole genome amplified DNA specimens were 100% concordant with whole genome amplified products from the same specimens. Whole genome amplified products were genotyped using multiplex PCR followed by a multiplex SNP analysis technique involving a single-base extension reaction (41). All genotyping was done blinded to case or control status.
PCR Amplification. Flanking primers were used in multiplexed PCR reactions to amplify the regions containing the four MDR1 SNPs, T129C (rs3213619), C1236T (rs1128503), G2677T/A (rs2032582), and C3435T (rs1045642). Standard PCR reactions (25 µL) contained 1.5 mmol/L MgCl2, 25 pmol of each deoxynucleotide triphosphate, 20 pmol of each primer, 1x Buffer B (Promega), 1.25 units Taq (Promega) and were cycled 38 times (30 s at 94°C, 230 s at 55°C, 30 s at 72°C). Subsequent to PCR, a subset of PCR reactions was assessed by agarose gel electrophoresis.
Single-Base Extension Procedure. In the multiplexed single-base extension reaction, primers of varied lengths (4 bp differences) were designed to hybridize to a region directly adjacent to the SNP and were allowed to single-base extend with a fluorescently labeled dideoxynucleotide to diagnose the SNP allele. PCR products were pooled into one solution so that each product had a concentration of
0.15 pmol/µL. Primers and unincorporated deoxynucleotide triphosphates were removed with shrimp alkaline phosphatase and exonuclease 1 (U.S. Biochemical). The single-base extension reaction (10.0 µL) contained 2.5 µL of SNaPshot Ready Reaction Mix (ABI Prism), 0.15 pmol/µL single-base extension primers, and 0.15 pmol/µL pooled PCR products. A total of 25 cycles of 96°C for 10 s, 50°C for 5 s, and 60°C for 30 s were done. Another shrimp alkaline phosphatase treatment removed unincorporated fluorescently labeled deoxynucleotide triphosphates, and the products were electrophoresed on a DNA sequencer (ABI 377). Automated genotyping software (GeneScan 5.0) was used to call genotypes.
Statistical Analysis
To account for the matched design of the study, conditional logistic regression analysis was used to examine the association between the four MDR1 SNPs and risk of childhood ALL. Odds ratios (OR) and 95% confidence intervals (95% CI) were calculated for the heterozygous variant and homozygous variant genotype (codominant model), and these two genotypes combined (dominant model) with the homozygous wild-type genotype serving as the reference group. Haplotype frequency and haplotype pair probability estimation of phase-unknown genotypes of C1236T, G2677T/A, and C3435T was done using the expectation-maximization algorithm (42) as implemented in SAS/Genetics PROC HAPLOTYPE (43). The haplotype determination was limited to these three SNPs due to the relatively low frequency of the T129C variant in our study population. A likelihood ratio test was done over all haplotypes to examine the association between MDR1 haplotypes and risk of ALL. Estimates of exact P values were computed using Monte Carlo methods with 10,000 permutations to address limitations (due to sparse data) of the
2 approximation for the distribution of the likelihood ratio test (44). Haplotype-specific ORs were estimated using conditional logistic regression to model the log odds of disease as a function of the individuals' haplotype probabilities (45, 46). Diplotype analyses were conducted using the most likely haplotype pair inferred for each subject by the expectation-maximization algorithm. To evaluate the joint effects of MDR1 haplotypes and indoor insecticide exposure, an OR for interaction (ORGE / ORGE' x ORG'E) and 95% CIs were calculated using a conditional logistic regression model that included a multiplicative interaction term and adjusted for household income. Analyses were done among all race/ethnicities combined and in non-Hispanic White and Hispanic children separately, the two largest ethnic groups.
| Results |
|---|
|
|
|---|
The allele frequencies of T129C, C1236T, and C3435T in non-Hispanic White control children were 5.1%, 44.6%, and 53.8%, respectively (Table 1 ). The frequencies for G2677T/A were 43.3% for the variant T-allele and 3.1% for the variant A-allele. Hispanic control children had similar allele frequencies for the four polymorphisms except for C3435T where non-Hispanic Whites showed a higher frequency of the variant allele (P = 0.032). The allele frequencies among other race/ethnicities are not presented due to their small numbers in our sample population. Our data showed evidence of pairwise linkage disequilibrium between the three common SNPs, C1236T, G2677T/A, and C3435T, among both non-Hispanic White and Hispanic controls. All observed genotype frequencies do not deviate significantly from the estimated frequencies expected assuming Hardy-Weinberg equilibrium.
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
In contrast to the only other study of MDR1 genetic variation and childhood ALL conducted among 113 cases and 175 controls (29), our study did not provide consistent evidence of an association between MDR1 SNPs individually and childhood ALL risk. In their study conducted in Poland, Jamroziak et al. (29) observed a nearly 2-fold increased risk of childhood ALL associated with the 3435TT genotype compared with the CC and CT genotypes combined. We were not able to confirm these results in our analysis of total ALL conducted in a slightly larger sample of non-Hispanic White cases (n = 135) and controls. However, in an analysis composed of the hyperdiploid ALL (>50 chromosomes) cases (n = 27) and individually matched controls (n = 35), we did observe a statistically significant association with the 1236TT, 2677TA/TT/AA, and 3435TT genotypes compared with the wild-type genotypes. We did these analyses in homogeneous ethnic populations due to emerging differences in risk factor associations and tumor genetic subtype frequencies in our population (47, 48). This association among non-Hispanic White hyperdiploid ALL cases and controls was also evident in the haplotype analysis where the two most common haplotypes, CGC and TTT, seemed to be associated with a reduced and increased risk, respectively. Because it is possible that these preliminary results may reflect a chance finding due to the limited numbers of subjects available for this subgroup analysis, the hyperdiploid ALL-specific effect should be interpreted with caution and will be reevaluated in future analyses with a larger sample size.
The inconsistency in results between the ALL analysis of this study and the Polish study may be related to possible genetic differences between the two Caucasian populations as reflected by the greater allele frequency of 3435T observed in our California population compared with their Polish population, as well as likely differences in linkage disequilibrium patterns with functional variants. The allele frequencies of the four SNPs observed in our non-Hispanic White study population are consistent with those reported in previous studies conducted in large European White populations (21, 22). The synonymous SNP, C3435T, has recently received tremendous attention, which is exemplified by an already large, yet controversial body of literature that suggest it be associated with an altered expression and/or functionality of the protein product (26). However, as C3435T does not lead to an amino acid change, it is believed that a yet unidentified causal SNP in strong linkage disequilibrium located in the MDR1 gene or elsewhere in another gene may be the true risk-conferring SNP. Alternatively, a new study reports potential differences in folding as a result of translational rate differences between the various haplotypes of MDR1 (49). Several studies, including ours, have found evidence of two common MDR1 SNPs, C1236T and G2677T/A, to be in linkage disequilibrium with C3435T (23-25, 36, 37). Tang et al. (37) did a haplotype analysis of the three common MDR1 SNPs, C1236T, G2677T/A, and C3435T, and observed clear race/ethnicity-specific differences in frequency between Chinese, Malayan, and Indians. The Chinese and Malayans exhibited three major haplotypes (CGC, TGC, and TTT) whereas the Indians exhibited only two of these three (CGC and TTT). According to the haplotype profile of our control population, Hispanics and non-Hispanic Whites had similar high frequencies of CGC and TTT haplotypes, but significantly differed with respect to CGT and TGC haplotypes. Thus, it is likely that the inconsistent reports on the functional nature of the various MDR1 SNPs may be explained by differences in linkage disequilibrium patterns observed between specific populations (26, 50). These differences will also affect substrate specificity differences observed among synonymous SNP haplotypes of MDR1, which were stronger than single SNP variants (49).
Our data suggest that variation in the MDR1 gene may modify the effect of household pesticide exposure on risk of childhood ALL. Specifically, compared with noncarriers of haplotype CGC, being a carrier of the haplotype seemed to significantly dampen the effect of indoor insecticide exposure on risk childhood ALL. This analysis was motivated by previous findings from a Northern California Childhood Leukemia Study analysis in which Ma et al. (18) reported an increased risk of childhood leukemia associated with household pesticides, particularly for insecticides used indoors during the prenatal period. Previous studies generally support an association between household pesticides and childhood leukemia (4), and one other study has also reported evidence of effect modification of this association, but by genetic variants of CYP1A1, a phase I xenobiotic metabolism gene (20). Although the specific mechanisms differ in the way these two genes, CYP1A1 and MDR1, interact with pesticide exposure, both studies provide additional support for a role of household pesticides in the etiology of childhood ALL and support the hypothesis that the risk associated with pesticides may be influenced by individual genetic susceptibility governed, in part, by genes involved in chemical transport and metabolism.
Previous studies evaluating the effects of these specific MDR1 SNPs and haplotypes on P-glycoprotein function and mRNA expression support this association. Briefly, Salama et al. (51) studied the effects of C1236T, G2677T/A, and C3435T on P-glycoprotein functionality in vitro using validated stable recombinant epithelial cells and observed reduced P-glycoprotein activity associated with all three variants both individually and as haplotypes when compared with the wild-type. Song et al. (52) evaluated whether genotypes and haplotypes of G2677T/A and C3435T were associated with hepatic MDR1 mRNA expression using human liver biopsy samples. Both homozygous variant genotypes individually (2677TT and 3435TT) were associated with significantly lower MDR1 mRNA expression levels compared with the wild-type genotype, and carriers of the GC haplotype of these two loci had a significantly higher mRNA expression level compared with noncarriers. Results of both these studies on functionality and mRNA expression provide evidence in support of a protective effect of the wild-type alleles or a potentially adverse effect of variant alleles of the SNP loci.
The ethnic diversity of our California study population offered us the unique opportunity to perform subgroup analyses in non-Hispanic White and Hispanic children separately. Because allele frequencies and the results of the association analyses did not seem to differ significantly between these two race/ethnicity groups, we reported the results of the analyses that were conducted among the entire case and control series. Our current sample size did not allow us to pursue analyses separately among the less common race/ethnicities in the study population. The individual matching of cases and controls on maternal race and child's Hispanic status should minimize the potential bias from population substructure and ethnic admixture differences between cases and controls. A limitation of this study is the reduced power to detect interactions between indoor insecticide exposure and the rare haplotypes. For this reason, we focused on the evaluation of the two most common haplotypes, CGC and TTT, which each had a frequency of
40% in our study population. Another limitation is the inability to evaluate the influence of specific time windows of exposure to indoor insecticides in the interaction analysis because a large proportion of subjects have used indoor insecticides across multiple time periods. For example, >50% of cases used indoor insecticides during both the prenatal and postnatal periods, whereas only 9.2% used indoor insecticides during the prenatal period only and 20.8% during the postnatal period only. When using a crude assessment of timing by categorizing subjects based on whether indoor insecticides were used during a certain time period (yes or no), there was strong evidence of interaction between haplotype CGC and indoor insecticide use during the prenatal period, but not for the postnatal period (data not shown). Furthermore, an evaluation of the importance of the prenatal period should include data on maternal MDR1 genotypes because studies have shown P-glycoprotein to be highly expressed in human placenta (25). The relative efficiency of maternal cells of the placenta in extruding hydrophobic xenobiotics, including pesticides, could affect the level of exposure to the fetus.
As the first to examine the role of MDR1 haplotypes in the risk of childhood leukemia, this study provides preliminary evidence of a differential effect of indoor insecticide exposure on childhood ALL risk between carriers and noncarriers of haplotype CGC. These results need to be replicated in an independent population, and future studies should additionally consider other membrane-bound transport proteins in combination with an environmental assessment of relevant substrates in a study of gene-environment interaction in childhood ALL etiology.
| Acknowledgments |
|---|
| Footnotes |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Note: J.L. Wiemels is a Scholar of the Leukemia and Lymphoma Society of America.
Received 1/11/07; revised 3/ 6/07; accepted 3/13/07.
| References |
|---|
|
|
|---|
in normal human T lymphocytes. Blood 1996;88:174754.This article has been cited by other articles:
![]() |
A. P. Chokkalingam and P. A. Buffler Genetic susceptibility to childhood leukaemia Radiat Prot Dosimetry, December 1, 2008; 132(2): 119 - 129. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Metayer and P. A. Buffler Residential exposures to pesticides and childhood leukaemia Radiat Prot Dosimetry, December 1, 2008; 132(2): 212 - 219. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |