Statistics from Altmetric.com
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.
Parkinson’s disease is the second most common neurodegenerative disease after Alzheimer’s disease. It is characterised by bradykinesia, rigidity, resting tremor, and postural instability.1 It is a genetically heterogeneous disorder. Pathogenic mutations in several genes—including α-synuclein, Parkin, UCH-L1 (ubiquitin-C terminal hydrolase-L1) and DJ-1—have previously been identified in rare monogenic forms of this disease showing autosomal dominant, autosomal recessive, or maternal transmission, with or without genetic anticipation.2,3 The more common, sporadic form of Parkinson’s disease appears to result from an interaction between genetic and environmental factors.4 Polymorphisms in several genes, including those implicated in familial forms of the disease such as α-synuclein5 and Parkin,6,7 are also reported to be associated with the sporadic form.8
Genetic susceptibility to sporadic Parkinson’s disease was also found to be modulated by genes involved in xenobiotic management. A meta-analysis of 84 association studies of 14 genes showed that polymorphisms in four genes are significantly associated with the disease.9 These genes are either responsible for xenobiotic metabolism, such as NAT210,11 and GSTT1,12 or may interact with environmental agents, such as monoamine oxidase (MAOB).13 Poor metaboliser alleles of the cytochrome P450 xenobiotic metabolism enzyme, CYP2D6, may also be associated with increased risk of Parkinson’s disease.14–20 Furthermore, there may be sex effects in the association of CYP2D6 mutant alleles with Parkinson’s disease.21
These genetic association studies corroborate epidemiological studies, which have long suggested that Parkinson’s disease is associated with exposure to certain environmental xenobiotics. Although most of the specific agents remain to be identified, rural living, well water consumption, industrialisation, and herbicide/pesticide exposure have been implicated as potential risk factors.1,22,23
Another category of genes that may influence susceptibility to Parkinson’s disease is the ATP binding cassette (ABC) superfamily of transporter genes which regulate the bioavailability of xenobiotics within critical tissues and cells in the body, of which the MDR1 multidrug transporter or P-glycoprotein is the best characterised member. Unlike drug metabolising enzymes, whose major drug metabolising functions occur in the liver, the MDR1 transporter is expressed at the interface of major organs. This pattern of distribution suggests that the MDR1 transporter regulates the traffic of drugs and xenobiotics in the body at two levels: its expression in the epithelial cells of the gut serves as a first initial barrier regulating the absorption of xenobiotics into the body, while its expression at the blood–brain and blood–germ cell/fetal interface serves as a second barrier controlling the uptake of xenobiotics into these sensitive tissues.24
The importance of the MDR1 transporter as a component of the blood–brain barrier is evident in knockout mouse studies. Mdr1a(−/−) mice were found to accumulate toxic levels of the anticancer drug, vinblastine, in the brain.25 Also, loperamide—an antidiarrhoeal narcotic analogue that normally does not enter the central nervous system (CNS)—was found to enter the brain of mdr1a(−/−) mice, causing them to develop abnormal behaviour characteristic of toxicity to CNS permeable opiates (for example morphine).26 Hence, we hypothesised that functional polymorphisms in the MDR1 gene may compromise its blood–brain barrier transporter function, increase accessibility of neurotoxic xenobiotics to the brain, and result in increased susceptibility to Parkinson’s disease.
Seven single nucleotide polymorphisms (SNPs) spanning ~100 kb of the MDR1 gene were examined in 206 Chinese patients with Parkinson’s disease and 224 matched normal controls.
Three SNPs—e12/1236(C/T), e21/2677(G/T/A), and e26/3435(C/T)—showed a significant association with Parkinson’s disease. In particular, e12/1236T, e21/2677T, and e26/3435T, or haplotypes containing these alleles, were found to be overrepresented in the matched normal controls compared with the Parkinson patients.
The significant effects of these SNPs were primarily observed in men and in patients with age of onset ⩾60 years; they were not associated with significant risk for Parkinson’s disease in women or in patients with a younger age of onset (⩽55 years).
It appears that the MDR1 transporter is a significant modulator of susceptibility to Parkinson’s disease among male ethnic Chinese ⩾60 years of age.
Several single nucleotide polymorphisms (SNPs) have been identified in the MDR1 gene, of which two (e21/2677(G/T/A) and e26/3435(C/T)) have been reported to be associated with differences in MDR1 expression and function, although the functional significance remains unclear. The non-synonymous SNP e21/2677(G/T/A) was reported to change the efflux of digoxin in cells in vitro in one study,27 but did not alter the efflux of several substrates in another study that used a different experimental system.28 The synonymous SNP e26/3435(C/T) has variously been associated with differences in MDR1 protein expression and plasma drug concentration,27,29–31 with drug induced side effects,32 and with drug response.33 Recently, these two SNPs and a third one, e1/-129(T/C), were examined in two case–control studies of approximately 100 patients with Parkinson’s disease and matched normal controls.34,35 No statistical significance was found between any of these SNPs and Parkinson’s disease.
In this study, we examined seven SNPs as well as haplotypes of these SNPs spanning ~100 kb in potentially functional regions of the MDR1 gene (that is, promoter region, coding regions, and 3′UTR) for an association with Parkinson’s disease. We found a significant association between Parkinson’s disease and the SNPs e12/1236(C/T), e21/2677(G/T/A), and e26/3435(C/T) (p values between 0.0367 and 0.00067), or haplotypes of these SNPs (p<0.05), in the Chinese population.
All patients with Parkinson’s disease and controls in this study were ethnic Chinese from Singapore. The Chinese in Singapore are predominantly descendents of migrants from south China. Individuals identified from the health screening programme in Singapore with no evidence of neurodegenerative disease on clinical examination were selected to serve as controls for the study. The diagnosis of Parkinson’s disease was made by neurologists specialising in movement disorders according to the United Kingdom Parkinson’s disease brain bank criteria.36 DNA was isolated from blood samples collected from 206 patients with Parkinson’s disease and 224 controls matched for age, sex, and ethnic group (table 1).
Ethical approval was obtained from the Singapore General Hospital research ethics committee.
The seven SNPs spanning ~100 kb of the MDR1 gene are located in five potentially functional genomic regions (promoter, exons 12, 21, 26, and 28) (fig 1). The five genomic segments were amplified in a single polymerase chain reaction, and all seven SNPs were genotyped by multiplex minisequencing as previously described.37
Genotype frequencies for the various SNPs in Parkinson’s disease patients and controls were assessed for deviation from Hardy–Weinberg equilibrium using Pearson’s χ2 test.38 A log-linear model embedded within the EM algorithm was used to estimate haplotype frequencies and haplotype–disease association.39,40 The analyses assumed Hardy-Weinberg equilibrium but allowed for linkage disequilibrium. A likelihood ratio test was used to assess whether haplotype–disease association models fitted better than models assuming no haplotype–disease association. As the likelihood ratio test assessed models rather than particular haplotypes, we also estimated odds ratios (OR) for each haplotype to quantify the strength and direction of the association of individual haplotypes, using the more prevalent haplotypes as reference. We obtained 95% confidence intervals (CI) of the odds ratios by the profile likelihood approach; a 95% CI that excluded the value of 1 indicated a significant relation between a particular haplotype and Parkinson’s disease risk.41,42 The EM algorithm estimation was carried out using the Stata program.41 All probability (p) values were two sided, and a p value smaller than 0.05 was considered significant.
SNPs with frequencies below 5% were excluded from the haplotype–disease association studies. In supplementary analyses, we examined the conditional independency of the excluded SNPs from Parkinson’s disease given the flanking SNPs by a likelihood ratio test,39 to determine whether the inclusion of these SNPs could improve the haplotype–Parkinson’s disease (haplotype–PD) association models given the flanking SNPs.
In subset analyses we further explored whether the association of the various alleles/haplotypes in the MDR1 gene with Parkinson’s disease differed between categories of sex and age of onset. As the average age of onset of Parkinson’s disease is around 60 years (table 1), early onset was defined as developing the disease at or before the age of 55, while late onset was defined as developing the disease at or after the age of 60. A gap of four years between 56 and 59 was not analysed, to allow for uncertainty in the ascertainment of the exact age of onset of some of the patients. Odds ratios and their confidence intervals were estimated separately in the different sex and age of onset groups. A sensitivity analysis was also carried out whereby we restricted the analysis of haplotype–disease association to subjects with phase-known haplotypes only. A logistic regression was used to estimate the odds ratio of disease.
As the genetic basis for complex disorders including Parkinson’s disease is still unclear, there could be extensive allelic variation at any disease locus, resulting in multiple susceptibility alleles of independent origin present in the population.43–45 It has been suggested that analysis of haplotypes rather than individual SNPs may be more advantageous in the presence of multiple susceptibility alleles at a single disease locus.43 In this study, we examined the association of individual SNPs as well as SNP haplotypes with Parkinson’s disease in ethnic Chinese.
Pearson’s χ2 test showed that all seven SNPs in our study population were consistent with the Hardy–Weinberg equilibrium assumption (each p>0.05).
Association of MDR1 SNPs and their haplotypes with Parkinson’s disease
As shown in table 2, the C allele of SNP e12/1236(C/T) (OR 1.341 (95% CI, 1.022 to 1.766); p = 0.0367), the G allele of SNP e21/2677(G/T/A) (OR 1.765 (1.317 to 2.365); p = 0.00067), and the C allele of SNP e26/3435(T/C) (OR 1.622 (1.223 to 2.161); p = 0.00074) were individually significantly associated with a higher risk of developing Parkinson’s disease. These three SNPs have previously been shown to be in tight linkage disequilibrium in the Chinese population.46 Calculated p values for all the possible haplotypes containing the above SNPs showed significant associations between these SNP combinations and Parkinson’s disease (p = 0.04321 to 0.00147), except for three combinations containing SNP i-1/-41(A/G) (i-1/-41(A/G)-e12/1236(C/T) (p = 0.1511), i-1/-41(A/G)-e12/1236(C/T)-e21/2677(G/T/A)-e26/3435(C/T) (p = 0.1041), and i-1/-41(A/G)-e12/1236(C/T)-e21/2677(G/T/A)-e26/3435(C/T)-e28/4036(A/G) (p = 0.4751)) (table 2). Even so, some specific haplotypes within these three SNP combinations were individually found to be associated with an increased risk of Parkinson’s disease (table 2).
SNPs e1/-145(C/G) and e1/-129(T/C) were excluded from the haplotype–association analyses as the minor alleles of these SNPs occur at less than 5% frequency. To evaluate whether the inclusion of these two SNPs would improve the haplotype–PD association models, we undertook conditional independence tests of the two SNPs from Parkinson’s disease, given the flanking SNPs by the likelihood ratio test. It was found that these two SNPs did not improve the haplotype–PD association model significantly (each p>0.05).
A sensitivity test using only phase-known haplotypes yielded similar results as EM estimated haplotype frequencies (data not shown), suggesting that the EM estimated haplotype frequencies were reliable.
Sex differences in risk determination
The characteristics of male and female Parkinson’s disease patients in our study population were found to be different. The women tended to be older and to have a later age of disease onset than the men (p<0.05) (table 1). We proceeded to examine whether there are sex specific associations between SNPs/haplotypes of the MDR1 gene and Parkinson’s disease. Our results showed that only haplotypes e12/1236C-e21/2677G-e26/3435C (OR 1.835 (95% CI, 1.082 to 3.175) and e21/2677G-e26/3435C-e28/4036A (OR 1.739 (1.012 to 2.996) were significantly associated with Parkinson’s disease in women (table 3). In contrast, most of the MDR1 SNPs and haplotypes that were significant in table 2 were also significant in men (table 3). Only SNP e12/1236C, and haplotypes i-1/-41A-e12/1236C, e12/1236C-e21/2677G-e26/3435C-e28/4036A, and i-1/-41A-e12/1236C-e21/2677G-e26/3435C-e28/4036A were not significantly associated with Parkinson’s disease in men, although their association with the disease in the overall population was significant. In addition, haplotypes e21/2677G-e26/3435C-e28/4036G (OR 3.644 (1.652 to 7.269), e12/1236T-e21/2677G-e26/3435C-e28/4036G (OR 2.715 (1.182 to 6.888), e12/1236C-e21/2677G-e26/3435C-e28/4036G (OR 5.778 (1.286 to 93.813), and i-1/-41A-e12/1236T-e21/2677G-e26/3435C-e28/4036G (OR 2.804 (1.090 to 7.126) were significantly associated with Parkinson’s disease in men but not overall (table 3).
Role of SNPs/haplotypes in the MDR1 gene in later onset of Parkinson’s disease
Interesting observations were made when we examined the age of onset specific association of SNPs/haplotypes in the MDR1 gene with Parkinson’s disease. While the promoter SNP i-1/-41(A/G) was found not to be associated with Parkinson’s disease in our overall or sex specific analyses, the low frequency G allele of this SNP was found to be significantly associated (p = 0.01), with a decreased risk of developing Parkinson’s disease at or before the age of 55 years (OR 0.307 (95% CI, 0.125 to 0.758) (table 4). Conversely, SNPs e21/2677(G/T/A) (p = 0.0102), e26/3435(C/T) (p = 0.0061), and SNP combinations e26/3435(C/T)-e28/4036(A/G) (p = 0.0423) and e21/2677(G/T/A)-e26/3435(C/T)-e28/4036(A/G) (p = 0.0225) were associated with increased risk of developing Parkinson’s disease at or after age 60, with SNPs e21/2677G (OR 1.748 (1.209 to 2.534)) and e26/3435C (OR 1.642 (1.148 to 2.354)), and haplotypes e26/3435C-e28/4036A (OR 1.657 (1.081 to 2.583)) and e21/2677G-e26/3435C-e28/4036A (OR 1.963 (1.250 to 3.106)) being associated with the increased risk (table 4). Some haplotypes that include either or both of the SNPs e21/2677(G/T/A) and e26/3435(C/T) were also associated with an increased risk of developing Parkinson’s disease (table 4). Curiously, although SNPs i1/-41(A/G) and e12/1236(C/T) were not individual risk factors, the haplotype i-1/-41A-e12/1236C (OR 1.470 (1.005 to 2.143)) was significantly associated with increased risk of late onset Parkinson’s disease (table 4).
Overall, the results from table 4 suggest that SNP i-1/-41(A/G) may be associated with decreased risk for developing Parkinson’s disease at or before the age of 55, while SNPs e21/2677(G/T/A) and e26/3435(C/T) and haplotypes containing these SNPs are associated with later onset disease (⩾60 years).
Environmental xenobiotics have been implicated in the development of Parkinson’s disease, a complex genetically heterogeneous disorder.1,22,23 The blood–brain barrier plays an important role in regulating the traffic of environmental xenobiotics in the brain, and individual differences in the “quality” of this barrier may influence the susceptibility to Parkinson’s disease. The MDR1 multidrug transporter represents an important component of the blood–brain barrier and has been shown to regulate the uptake of drugs and xenobiotics into this sensitive organ.25,26,47 It is conceivable that polymorphisms which alter the expression levels or transport ability of this transporter could result in altered susceptibility to neurotoxic substances and thus alter the genetic threshold for the development of Parkinson’s disease.
Two recent case–control studies have examined the role of MDR1 gene polymorphisms (SNPs e1/-129(T/C), e21/2677(G/T/A), and e26/3435(C/T)) in Parkinson’s disease development. The studies involved approximately 100 white Italian and Polish patients and 100 controls from the same geographical regions.34,35 No significant associations between these SNPs and Parkinson’s disease were detected. However, our present study of 206 Chinese patients and 224 controls showed that three SNPs—e12/1236(C/T) (p = 0.0367), e21/2677(G/T/A) (p = 0.00067), and e26/3435(C/T) (p = 0.00074), all in tight linkage disequilibrium with each other46—are significantly associated with an altered risk of developing Parkinson’s disease (table 2). The odds ratios of the haplotypes that were associated with Parkinson’s disease were not very high. These observations are, however, consistent with the widely held view that Parkinson’s disease is a complex disorder involving the interaction of multiple genes with different environmental factors, whereby the individual contribution of each causative gene may not be large.
We recently found strong evidence of positive selection for the e21/2677T and e26/3435T alleles in the Chinese, but only marginal evidence for this in white Americans (Tang K, Wong L, Lee E, et al, Human Molecular Genetics (in press)). The Chinese samples in that study were from anonymised umbilical cord blood from Chinese neonates, and allele frequencies of the seven SNPs were found to be very similar to those in the present study. When we used cord blood DNA samples as controls to compare against the Parkinson’s disease samples, we obtained a similar, statistically significant association between Parkinson’s disease and these two SNPs (data not shown). The strong evidence of a recent positive selection for the T alleles of these two SNPs supports our current observation that these alleles are significantly underrepresented in patients with Parkinson’s disease compared with unaffected controls, suggesting that the T alleles of these SNPs may confer better protection for the brain against xenobiotic insults in the Chinese population.
It is possible that the earlier Italian and Polish association studies did not detect a significant statistical association because of their limited sample size. There may be another reason why neither study was able to detect a significant association between any MDR1 SNPs and Parkinson’s disease. If we assume that the Italian and Polish subjects34,35 were genetically similar to white Americans, their MDR1 haplotype and LD profiles may not favour the detection of associations. Our observation of only marginal evidence of recent positive selection in white Americans compared with the Chinese supports this hypothesis. Nonetheless, it remains to be determined whether the white Italians and Poles are in fact similar to white Americans in their underlying genetic architecture at this locus.
It is possible that either SNP e21/2677(G/T/A) or e26/3435(C/T) could be potential causal SNPs as they had much lower p values than SNP e12/1236(C/T). Consistent with our observation that individuals carrying the G allele at the non-synonymous SNP e21/2677(G/T/A) have a higher risk of developing Parkinson’s disease, the MDR1 transporter carrying the e21/2677G allele—coding for Ala at amino acid position 893—has been shown to be a less effective transporter than one carrying the T allele (Ser 893).27 The synonymous SNP e26/3435(C/T) appears to be associated with altered MDR1 transporter expression and function. While several reports found that the T allele is associated with lower MDR1 expression,29,30,33,48 resulting in lower efflux or higher plasma levels of drugs and xenobiotics,29,30 others have reported lower drug plasma concentration in individuals carrying the T allele.27,31,33 Most of these studies examined only SNP e26/3435(C/T) without taking into account the underlying haplotype and linkage disequilibrium architecture of the study population. Detailed characterisation of the genetic and evolutionary history of the entire MDR1 gene in each study population, and the influence of recent events in the history of each population on linkage disequilibrium and the likelihood of detecting an association, could resolve these conflicting reports. Our data showing an association between e26/3435T and a lower risk of developing Parkinson’s disease support observations that the T allele alters MDR1 function, resulting in a greater efflux of drugs or xenobiotics. Although SNP e26/3435(C/T) is a synonymous SNP and does not result in an amino acid change, there are several possible explanations for this observation. The observed correlation with e26/3435T could reflect either differential codon usage of the C or T allele at the wobble position of the isoleucine codon, or allele specific differences in RNA folding,49 sometimes influencing RNA processing50 or splicing,51,52 or differences in translation control53 and regulation.54 It is also possible that neither SNP e21/2677(G/T/A) nor e26/3435(C/T) represents the causal SNP, but that they are merely in strong linkage disequilibrium with an unobserved causal SNP. A strong association of these two SNPs with Parkinson’s disease could suggest that the linked causal variant resides within a region defined by strong LD.
An interesting observation was made when male and female patients with Parkinson’s disease were investigated independently—the MDR1 gene appears to play a more important role in determining risk of developing the disease in men than in women (table 3). This is consistent with the view that the MDR1 transporter regulates the accumulation of neurotoxic xenobiotics in the brain to modulate the risk of developing Parkinson’s disease. As older women in urban Singapore are primarily home makers while men often work out of doors, it is conceivable that the observed greater risk for Parkinson’s disease in men compared with women is related to increased exposure to environmental susceptibility factors among men, given the same genetic risk factors in the two sexes.
When patients with Parkinson’s disease were compared on the basis of their age at disease onset, we found that several polymorphisms in the MDR1 gene seemed to play a greater role in later onset disease (⩾60 years) (table 4). One hypothesis is that, in individuals with particular MDR1 genotypes (for example, e12/1236C, e21/2677G, e26/3435C) and haplotypes, the blood–brain barrier allows neurotoxic xenobiotics easier access and gradual accumulation in the brain, eventually leading to Parkinson’s disease. Conversely, individuals with the alternative alleles (that is, e12/1236T, e21/2677T and e26/3435T) are better protected from xenobiotic insults and hence from Parkinson’s disease. In contrast, early onset Parkinson’s disease is probably a result of other genetic factors and hence is less dependent on genetic variation at the MDR1 locus.
The promoter SNP i-1/-41(A/G), which resides in a putative CCAAT box, was found to influence the risk of Parkinson’s disease in patients with a younger age of onset (p = 0.01) (table 4). The G allele of this SNP appeared to protect individuals from Parkinson’s disease (OR 0.307 (95% CI, 0.125 to 0.758)). This observation, however, should be interpreted cautiously, given the low frequency (<10%) of i-1/-41G in the general population and the resultant sample sizes in this comparison.
We have produced strong statistical evidence that particular alleles and haplotypes of MDR1 SNPs—e12/1236(C/T), e21/2677(G/T/A), and e26/3435(C/T)—are important risk factors for the development of Parkinson’s disease in ethnic Chinese, especially in men, through sex associated lifestyle differences, and in individuals with a later age of onset (⩾60 years). The wide variations in allele frequencies of the MDR1 SNPs (especially SNP e12/1236(C/T), e21/2677(G/T/A), and e26/3435(C/T)) among different ethnic populations46 may account for the differences in the ability to detect an association between MDR1 and Parkinson’s disease in other ethnic groups, especially if the increase in relative risk is small.
This study was supported by a grant from the National Medical Research Council, Singapore (NMRC/0657/2002) to CGL, SSC, and EJDL.
Conflicts of interest: none declared