Skip to main content
  • Research article
  • Open access
  • Published:

Genetic associations of Nrf2-encoding NFE2L2 variants with Parkinson’s disease – a multicenter study

Abstract

Background

The transcription factor Nrf2, encoded by the NFE2L2 gene, is an important regulator of the cellular protection against oxidative stress. Parkinson’s disease is a neurodegenerative disease highly associated with oxidative stress. In a previously published study, we reported associations of NFE2L2 haplotypes with risk and age at onset of idiopathic Parkinson’s disease in a Swedish discovery material and a Polish replication material. Here, we have extended the replication study and performed meta-analyses including the Polish material and four new independent European patient-control materials. Furthermore, all SNPs included in the haplotype windows were investigated individually for associations with Parkinson’s disease in meta-analyses including all six materials.

Methods

Totally 1038 patients and 1600 control subjects were studied. Based on previous NFE2L2 haplotype associations with Parkinson’s disease, five NFE2L2 tag SNPs were genotyped by allelic discrimination and three functional NFE2L2 promoter SNPs were genotyped by sequencing. The impact of individual SNPs and haplotypes on risk and age at onset of Parkinson’s disease were investigated in each material individually and in meta-analyses of the obtained results.

Results

Meta-analyses of NFE2L2 haplotypes showed association of haplotype GAGCAAAA, including the fully functional promoter haplotype AGC, with decreased risk (OR = 0.8 per allele, p = 0.012) and delayed onset (+1.1 years per allele, p = 0.048) of Parkinson’s disease. These results support the previously observed protective effect of this haplotype in the first study. Further, meta-analyses of the SNPs included in the haplotypes revealed four NFE2L2 SNPs associated with age at onset of Parkinson’s disease (rs7557529 G > A, −1.0 years per allele, p = 0.042; rs35652124 A > G, −1.1 years per allele, p = 0.045; rs2886161 A > G, −1.2 years per allele, p = 0.021; rs1806649 G > A, +1.2 years per allele, p = 0.029). One of these (rs35652124) is a functional SNP located in the NFE2L2 promoter. No individual SNP was associated with risk of Parkinson’s disease.

Conclusion

Our results support the hypothesis that variation in the NFE2L2 gene, encoding a central protein in the cellular protection against oxidative stress, may contribute to the pathogenesis of Parkinson’s disease. Functional studies are now needed to explore these results further.

Peer Review reports

Background

Parkinson’s disease (PD) is a neurodegenerative disease affecting the central nervous system, resulting in motor symptoms such as rigidity, slowness of movement, postural instability and a characteristic resting tremor. The motor symptoms derive from a decrease in the neurotransmitter dopamine (DA), due to death of dopaminergic (DAergic) neurons in the substantia nigra (SN). Even though several PD-causing genes have been identified, the majority of PD patients have an idiopathic form without known cause and the disease process may involve a combination of genetic and environmental factors (reviewed in [1]). PD pathogenesis is known to involve oxidative stress and the SN seems to be especially vulnerable. Incomplete intracellular oxidation of DA in the SN, favored by the presence of ferrous iron, may result in the formation of reactive DA(semi)-quinones, which together with decreased levels of reduced glutathione and chronic inflammation adds to the oxidative stress observed in the PD brain (reviewed in [2],[3]).

Nuclear factor erythroid 2 (NF-E2) related factor 2 (Nrf2), a transcription factor encoded by the NF-E2-like 2 (NFE2L2) gene, has a key role in the cellular protection against oxidative and electrophilic insults [4]. In the normal state, Nrf2 is kept largely inactive by its repressor protein kelch-like ECH-associated protein 1 (Keap1) [5], which targets Nrf2 for ubiquitin ligation and subsequent degradation by the proteasome [6],[7]. When Nrf2 is activated by oxidative or electrophilic stress, it induces transcription of a battery of cytoprotective genes by binding to a specific region in their promoters – the antioxidant response element (ARE) [8],[9]. Nrf2 can also be activated by dietary factors, such as sulforaphane or curcumin, and Nrf2-activating substances have been used in a number of pre-clinical PD models to study the effect of Nrf2 upregulation on PD progression with a pharmacological perspective (reviewed in [10],[11]. Interestingly, one study has shown that an increased oxidative stress observed in olfactory neurosphere-derived cells from PD patients could be restored by activation of Nrf2 with sulforaphane [12] and curcumin has been shown to protect DAergic SH-SY5Y neurons from 6-Hydroxydopamine toxicity [13]. Another study has shown that upregulation of Nrf2 using potent synthetic Nrf2 activators protects DAergic neurons from degeneration in a 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) mouse model of PD [14].

Nrf2 is essential in regulation of the cellular redox homeostasis, as it controls rate-limiting steps in the neo-synthesis of glutathione [15], as well as the induction of the antioxidant and neuroprotective enzyme heme oxygenase 1 (HO-1) [16],[17]. Another protein regulated by Nrf2 is the highly inducible NAD(P)H dehydrogenase, quinone 1 (NQO1) [18]. Besides its broad general antioxidant activity, NQO1 prevents toxic redox cycling of DA-quinones through a two-electron reduction into stable hydroquinones [19]. NQO1 has been found to be expressed in the SN of PD patients but not in age-matched control subjects [20], suggesting an increased Nrf2 activity in the PD brain. Nrf2 nuclear translocation is enhanced in the PD affected SN [21] and a recent study showing an increased Nrf2 activity in neurons derived from induced pluripotent stem cells (iPSCs) from PD patients [22] indicates a role for Nrf2 activation in the PD pathogenesis. Considering all this, it is feasible that a decreased radical protection due to genetic variation in the NFE2L2 gene could affect the pathogenesis of PD.

We previously reported that a haplotype in the NFE2L2 gene, GAAAA consisting of five tag single nucleotide polymorphisms (SNPs), was associated with delayed age at onset (AAO) in a Swedish discovery material and with decreased risk of PD in an independent Polish replication material [23]. In the same study, investigation of three functional SNPs in the human NFE2L2 promoter, previously shown to influence Nrf2 protein expression [24],[25], revealed that the protective haplotype GAAAA was in linkage disequilibrium (LD) with the promoter haplotype AGC, which is part of a fully functional promoter [23].

Here, based on our previous NFE2L2 haplotype associations with PD in the Swedish Gothenburg (PD-Goth) discovery material, seven NFE2L2 haplotypes and one individual SNP were investigated for associations with risk and/or AAO of idiopathic PD in meta-analyses including five independent patient-control replication materials. Furthermore, with an exploratory approach despite lack of association with risk in our previous study, all genotyped SNPs were also independently investigated in meta-analyses including the PD-Goth discovery material.

Methods

Patient-control materials

In total, 2638 individuals (1038 PD patients and 1600 controls) were included in this study. Subjects originated from six independent research centers in Europe: Italy, Malta, Poland and Germany, as well as two independent patient-control sets from Sweden: Gothenburg (PD-Goth) and Linköping (PD-Link).

PD-Goth

PD-Goth was the discovery material in our previous study [23] and included 165 PD patients and 190 control individuals. All participants were of Caucasian origin and patients and control subjects were of similar age. The patients were diagnosed according to the Parkinson’s Disease Society Brain Bank criteria for idiopathic PD, except that presence of more than one relative with PD was not considered a criterion for exclusion. PD patients with an AAO of < 50 years were screened to exclude mutations in the recognized PD-causing genes DJ-1, Parkin, PINK1 and LRRK2 [26],[27].

Italy

The Italian material included 329 PD patients and 450 control subjects. All participants were of Caucasian origin and the patients were diagnosed according to the Parkinson’s Disease Society Brain Bank criteria for idiopathic PD. The controls were unrelated subjects originating from the same geographical area as the patients and all underwent neurological examination for exclusion of neurological diseases.

PD-Link

The PD-Link material included 195 PD patients and 379 control subjects collected as part of the Geoparkinson Study [28]. Blood samples were collected from L-dopa positive PD patients in southeast Sweden visiting the Clinic from Geriatrics and Neurology, University Hospital in Linköping, Sweden. The patients were classified as having Parkinson’s disease or parkinsonism using the United Kingdom Parkinson’s Disease Society Brain Bank clinical diagnostic criteria and a neurologist confirmed the diagnosis at recruitment. Individuals with vascular or drug-induced parkinsonism were excluded from the study, as were those with dementia. In this study, only samples from patients with confirmed PD were used. The controls were unrelated Swedish subjects randomly collected from the normal population in southeast Sweden, the same study base as the patients. The patient and control groups were frequency-matched by age and sex.

Malta

The Malta material included 101 PD patients and 313 control subjects collected as part of the Geoparkinson Study [28]. All participants were of Caucasian (Maltese) origin. The patients were classified as having Parkinson’s disease or parkinsonism using the United Kingdom Parkinson’s Disease Society Brain Bank clinical diagnostic criteria. In this study, only samples from patients with confirmed PD were used. The controls were from the community or from out-patients at St Luke’s Hospital, G’Mangia, Malta. They were group-matched to the patients by age and sex.

Poland

The Polish material was the replication material in our previous study. It included 192 PD patients and 192 control subjects. All participants were of Caucasian origin from the same geographic area, and were matched by sex. Age at sampling (AAS) was significantly higher in the control subjects than in the PD patients (to minimize the risk of PD development in the controls later in life). The patients fulfilled the Parkinson’s Disease Society Brain Bank criteria for idiopathic PD except for the presence of more than one relative with PD, i.e. all patients with a family history of PD were excluded from the study.

Germany

The German material included 56 PD patients and 76 control subjects. All participants were of Caucasian origin and were matched by age and sex. PD was diagnosed by movement disorders specialists according to the Parkinson’s Disease Society Brain Bank criteria for idiopathic PD. The controls were also personally examined and free of any symptoms or signs suggestive of PD.

All materials contained information regarding sex and AAS of patients and control subjects, and AAO of patients (i.e. age at diagnosis (AAD) for Germany, Malta and PD-Link; and age at first symptom (AAFS) for PD-Goth, Poland and Italy). Some materials also included information on family history (FH; one or more 1st degree family member(s) with PD) and smoking status. Demographics are summarized in Table 1.

Table 1 Demographic characteristics of PD patients and control subjects

SNP selection

Selection of common SNP genotyping data covering NFE2L2 for the European material CEU (Utah residents with ancestry from Northern and Western Europe) were performed when designing our previous study [23]. In brief, SNP genotyping data was downloaded from the International Haplotype Mapping Project web site (http://www.hapmap.org) [29] and processed using the Haploview software [30]. LD blocks were constructed according to Gabriel et al. [31] and tag SNPs assigned using the tagger function [30]. A minor allele frequency of ≥ 5% and pair-wise tagging with a minimum r2 of 0.80 were applied to capture the common SNPs within the block covering the NFE2L2 gene. In the same study promoter SNPs were chosen due to their reported functional effect on Nrf2 transcription [24],[25]. The analyzed SNPs, included in the haplotype window previously associated with PD [23], were: five tag SNPs (rs7557529, rs2886161, rs1806649, rs2001350 and rs10183914, here SNPs 2 – 6) and three functional promoter SNPs [24],[25] (rs35652124, rs6706649 and rs6721961, here SNPs P1–P3 [23]), (Figure 1, Table 2).

Figure 1
figure 1

Schematic NFE2L2 gene positions of the studied SNPs and their LD-patterns, showing r 2 -values as 100r 2 .

Table 2 Overview of the studied SNPs

Tag SNP genotyping and promoter sequencing

A detailed description of tag SNP genotyping and sequencing of promoter SNPs is provided in our previous study [23]. In brief, all SNPs were genotyped using genomic DNA extracted from blood. All tag SNPs were genotyped using TaqMan Allelic Discrimination [32] with TaqMan® Pre-Designed SNP genotyping assays or TaqMan® Custom Made SNP genotyping assays (Applied Biosystems, Foster City, CA, USA). The promoter SNPs were genotyped by sequencing, amplifying a 423 bp region of the NFE2L2 promoter (forward primer 5’-GACCACTCTCCGACCTAAAGG-3’, reverse primer 5’-CGAGATAAAGAGTTGTTTGCGAA-3’, annealing temperature 59° C and 34 cycles on a PTC-200 ThermalCycler (Biorad, Hercules, CA, USA)). Purification of the PCR product was performed with an Illustra™ GFX™ PCR Purification Kit (GE Healthcare, Little Chalfont, Buckinghamshire, UK). Sequence reactions were performed using BigDye v3.1 (Applied Biosystems, Forster City, CA, USA) and analyzed on an ABI PRISM 3100 Automated Sequencer (Applied Biosystems, Forster City, CA, USA). Sequence data were analyzed with the DNASTAR SeqMan® software (DNASTAR Inc., Madison, WI, USA).

Statistical analyses

Demographic statistics were performed with SYSTAT11 (SYSTAT Software GmbH, Erkrath, Germany). Pearson χ2-statistics were used for sex, FH and smoking history; and Mann–Whitney U test for AAS.

Genetic association analyses were performed using HelixTree 6.3 (Golden Helix, Bozeman, MT, USA). All tag SNPs were analyzed for deviation from Hardy–Weinberg equilibrium. Based on the significant results in the PD-Goth material [23], haplotypes were analyzed for associations with PD (GAAAA, GAGGG and GAAAG in the haplotype window consisting of tag SNPs 2 – 6, AGC in the window consisting of promoter SNPs 1 – 3 and GAGCAAAA, GAGAAGGG and GAGCAAAG in the window consisting of all eight SNPs combined, see Figure 1 and Table 2 for SNP overview). A detailed description of how the haplotype windows were identified is given in our previous study [23], where the identified windows were corrected using permutations tests with 10 000 permutations and pc-values of ≤ 0.05 were considered statistically significant.

With meta-analysis sample sizes > 2000 individuals (after correction for unequal sample sizes) and a desired power of 80%, we can detect standardized differences down to approximately 0.13 in this study [33], which means that all the detected significances in risk and AAO reported in our previous paper [23] can be detected here at the applied significance level of 0.05.

SNP and haplotype associations were analyzed using logistic or linear haplotype regression, where each haplotype was analyzed individually relative to all the other haplotypes together. Sex was the only covariate available in all materials. In order to maximize sample size and power of the analyses and to make the results comparable between the independent materials, only sex was included as covariate in our model (except for risk of PD in the Polish material where the number of males and females were identical in the patient and control groups). However, when data was available and statistically relevant (see our previous study [23] for a detailed description of how statistically relevant covariates were selected), the effects of sex, FH and smoking habits were evaluated in the meta-analyses to assure that exclusion of covariates did not affect the outcomes. Haplotype frequencies were estimated in each material individually using the EM algorithm [34] allowing imputation to compensate for missing genotypes.

Meta-analyses were performed using inverse-variance weighting. In attempts to replicate our previously observed associations, the five replication materials were analyzed for associations of NFE2L2 with PD risk and AAO. One-sided p-values ≤ 0.05 were considered statistically significant. In the exploratory meta-analyses for individual SNP associations with PD risk and AAO, all six materials were included in the analyses. Two-sided p-values ≤ 0.05 were considered statistically significant.

Ethics

The study was approved by ethics committees in each respective country (Sweden Gothenburg: the regional ethics committee at University of Gothenburg, Sweden; Italy: the ethical committee at the University Hospital Mater Domini, Catanzaro, Italy; Sweden Linköping: the regional ethics committee at University of Linköping, Sweden; Malta: the research ethics committee at the University of Malta, Malta; Poland: the ethics committee of the Pomeranian Medical University, Szczecin, Poland; Germany: the ethics committee of the University of Lübeck). This study was in compliance with the Helsinki Declaration of 1975 and written informed consent was obtained from all participants.

Results

Demographics

PD-Goth: Patients and control subjects were similar in age, FH and current smoking status, but differed significantly in the distributions of sex and ever-smoker status. Italy: Patients and control subjects differed significantly in AAS and sex. No information was available regarding FH or smoking habits for the control subjects. PD-Link: Patients and control subjects did not differ with respect to FH, but differed significantly in AAS, sex and smoking status. Malta: Patients and control subjects were similar in AAS, sex and current-smoker status, but differed significantly with respect to FH and ever-smoker status. Poland: Patients and control subjects were matched in sex, but differed significantly in AAS due to the study design. No information was available regarding PD FH or smoking habits for either PD patients or control subjects. Germany: Patients and control subjects were similar regarding sex frequencies, but differed significantly in AAS. No information was available regarding PD FH or smoking habits for either PD patients or control subjects.

In general, all centers had a higher percentage of males relative to females amongst the patients (57.0 – 62.4%). Amongst the controls there was a wider range (36.8 – 62.9%) due to different study designs. Similarly, the average AAS was more similar between the patient groups (63.7 – 71.4 years) than between the control groups (57.0 – 72.9 years). The Polish material stood out, since the controls were older than the patients as a result of the study design. The average AAO varied between 55.2 and 64.1 years and the average time from AAO until AAS ranged from 6.4 to 9.2 years between the different patient groups. Where data was provided, FH of PD was always more common in the patient groups than in the controls, though this difference did not always reach statistical significance. As expected, smoking (ever smokers) were more common in all patient groups than in the control groups, whenever data was available. All demographics are given in Table 1. Frequencies for the studied NFE2L2 SNPs and haplotypes in patients and controls are summarized in Additional file 1: Table S1 and Additional file 2: Table S2, respectively.

Tag SNP genotyping

None of the studied SNPs had a Bonferroni-corrected Hardy-Weinberg equilibrium p-value of < 0.001. The overall call rate was > 95%.

Haplotype associations

In accordance with the findings in the first study, meta-analyses of NFE2L2 haplotypes including five independent replication patient-control materials showed a protective effect of the two haplotypes GAGCAAAA and GAAAA on PD. Haplotype GAGCAAAA, consisting of the five tag SNPs and the three functional promoter SNPs (SNPs 2, P1, P2, P3, 3, 4 ,5 and 6 according to gene direction, for rs numbers see Table 2), was associated with both decreased risk (OR = 0.8 per allele, p = 0.012) (Table 3, Figure 2A) and later PD onset (+1.1 years per allele, p = 0.048) (Table 4, Figure 2B). Haplotype GAAAA, consisting of tag SNPs 2 – 6 only, was associated with later PD onset (+1.3 years per allele, p = 0.024) (Table 4) and showed a tendency towards association with decreased PD risk (OR = 0.9 per allele, p = 0.052) (Table 3).

Table 3 Replication of haplotype associations with risk of PD
Figure 2
figure 2

Forest plot showing the NFE2L2 haplotype GAGCAAAA associations with decreased risk (A) and later onset (B) of PD in meta-analysis including the five replication materials: Italy, Sweden PD-Link, Malta, Poland and Germany (filled squares). The previous finding from the Swedish PD-Goth discovery study is also shown (empty squares), though the material was not included in the meta-analysis. The areas of the squares are proportional to the meta-analysis weights (inverse variance).

Table 4 Replication of haplotype associations with AAO of PD

None of the haplotypes GAGAAGGG, GAGGG, GAGCAAAG or GAAAG in the same two windows, all associated with increased risk in the PD-Goth material in the first study, was significantly associated with PD in the meta-analyses. The promoter haplotype AGC alone was not associated with PD in meta-analysis (Tables 3 and 4).

Individual SNP associations

Our previous protective association of the A allele of tag SNP 6 (rs10183914) in the PD-Goth discovery study [23] could not be replicated in the meta-analysis including the five replication materials (+0.3 year per A allele, p = 0.23).

Exploratory meta-analyses of NFE2L2 individual SNPs, including all six materials, showed no significant associations with risk of PD (Table 5).

Table 5 Individual SNP associations with risk of PD

Exploratory meta-analyses of NFE2L2 individual SNPs with AAO, including all six materials, showed associations with four SNPs; tag SNP 2 (rs7557529, −1.0 year per A allele, p = 0.042), promoter SNP 1 (rs35652124, −1.1 year per G allele, p = 0.045), tag SNP 3 (rs2886161, −1.2 year per G allele, p = 0.021) and tag SNP 4 (rs1806649, +1.2 year per A allele, p = 0.029) (Table 6).

Table 6 Individual SNP associations with AAO of PD

Discussion

Since oxidative stress is implicated in the pathogenesis of PD [2],[3] and Nrf2 is a central protein in the cellular defense against oxidative stress, genetic variation affecting the efficiency of Nrf2 could contribute to the disease. Here, based on the previous results from our group [23], we performed meta-analyses on five independent patient-control materials to investigate the effects of seven haplotypes in the Nrf2-encoding gene NFE2L2 on risk and AAO of PD. In addition, the included SNPs were analyzed individually in exploratory meta-analyses including the PD-Goth discovery material. To our knowledge, with 1038 PD patients and 1600 control subjects included, this is the most extensive attempt yet to investigate the impact of genetic variation in NFE2L2 on PD risk and AAO.

Meta-analyses of three NFE2L2 haplotypes containing five consecutive tag SNPs and three functional promoter SNPs showed association of haplotype GAGCAAAA with both decreased risk and later onset of PD (Tables 3 and 4, Figure 2). This supports our previously observed protective effect of this haplotype, which was associated with later onset of PD in the PD-Goth discovery material and decreased risk of PD in the Polish replication material in the first study [23]. In addition, meta-analyses of the three previously associated NFE2L2 haplotypes in the window including the five tag SNPs only (tag SNPs 2 – 6) showed association of haplotype GAAAA with later onset of PD. This haplotype also had a tendency towards association with decreased risk of PD. Obviously, GAGCAAAA and GAAAA are highly influenced by one another and should not be considered separate haplotypes. The promoter haplotype AGC is in LD with the PD-associated haplotype GAAAA (r2 = 0.4) [23] and the PD-associated combined haplotype GAGCAAAA contains the fully functional promoter haplotype AGC [24],[25]. This suggests that the associations could be driven by genetic variation in the NFE2L2 promoter. However, the fact that the promoter haplotype AGC alone showed no association with PD in meta-analysis indicates that other parts of the gene are involved as well. This is also in accordance with a previous study on a Taiwanese material, where none of the NFE2L2 promoter haplotypes were associated with PD [35].

Haplotypes GAGAAGGG, GAGGG, GAGCAAAG and GAAAG in the same two haplotype windows were all associated with increased risk of PD in the PD-Goth discovery material, but none of the associations could be replicated here. Likewise, the association of tag SNP 6 (rs10183914) with later PD onset in the PD-Goth discovery study could not be replicated in meta-analysis.

In line with our previous results, individual exploratory meta-analyses of the genotyped NFE2L2 SNPs on patient and control subjects from six independent European centers showed no associations with risk of PD. However, in this study four NFE2L2 SNPs were associated with AAO of PD in the exploratory meta-analyses. The minor alleles of tag SNP 2 (rs7557529), promoter SNP P1 (rs35652124) and tag SNP 3 (rs2886161) were all associated with one year earlier PD onset per allele. Tag SNP 3 of our study was previously included in a multiple candidate gene study, but was not significantly associated with PD [36]. Tag SNP 2 is in LD (r2 = 0.9) with SNP rs6726395 [23], which was previously found to increase the risk of PD in a first-tier, sib-pair whole-genome study of PD, though the association was not replicated in the second-tier case–control study [37]. The minor G allele of promoter SNP P1 has been shown to significantly decrease Nrf2 protein levels [24],[25]. Tag SNP 2 and especially tag SNP 3 of our study are both in LD with the minor G allele of P1 (r2 = 0.4 and r2 = 0.9, respectively) (Figure 1), which could possibly explain the observed associations of these SNPs with earlier PD onset in our study. The minor A allele of tag SNP 4 (rs1806649) was instead protective and delayed onset of PD with one year per allele in meta-analysis. Consistent with the discussion above, this SNP was not in LD with SNP P1 (r2 = 0.1). Activation of Nrf2 with natural or synthetic substances have been shown to reduce oxidative stress and neurodegeneration in different PD models [10],[11], which suggests that Nrf2 could be a suitable target for pharmacological intervention of PD. In this context, the associations of NFE2L2 variants with AAO of PD are interesting.

Though haplotype GAGCAAAA showed a trend towards increased risk of PD in the Italian material when the materials were analyzed individually, this haplotype decreased the risk of PD in meta-analysis (Table 3, Figure 2A). Idiopathic PD is a multifactorial disease caused by several risk genes, possibly in combination with environmental factors. This may explain why the effect of a risk gene could be more obvious in some materials than in others. Also, study design as well as different genetic background in both PD patients and controls could affect the outcome of a multicenter study. Considering this, despite the relatively small effects in the meta-analyses, the replicated associations between genetic variation in the NFE2L2 gene and idiopathic PD are interesting and should be investigated further in functional studies.

Conclusions

Meta-analyses including five independent patient-control materials from different European centers confirmed protective effects of NFE2L2 haplotypes GAGCAAAA and GAAAA against PD. In addition, exploratory meta-analyses revealed four individual SNPs affecting the AAO of PD, including a functional SNP located in the NFE2L2 promoter. In summary, our data support the hypothesis that variation in the Nrf2-encoding gene NFE2L2 may contribute to the pathogenesis of idiopathic PD.

Additional files

Abbreviations

AAD:

Age at diagnosis

AAFS:

Age at first symptom

AAO:

Age at onset

AAS:

Age at sampling

ARE:

Antioxidant response element

DA:

Dopamine

DAergic:

Dopaminergic

FH:

Family history

HO-1:

Heme oxygenase 1

iPSCs:

Induced pluripotent stem cells

Keap1:

Kelch-like ECH-associated protein 1

LD:

Linkage disequilibrium

MPTP:

1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine

NFE2L2 :

NF-E2-like 2

Nrf2:

Nuclear factor erythroid 2 (NF-E2) related factor 2

NQO1:

NAD(P)H dehydrogenase, quinone 1

PD:

Parkinson’s disease

SNP:

Single nucleotide polymorphism

SN:

Substantia nigra

References

  1. Lees AJ, Hardy J, Revesz T: Parkinson’s disease. Lancet. 2009, 373 (9680): 2055-2066. 10.1016/S0140-6736(09)60492-X.

    Article  CAS  PubMed  Google Scholar 

  2. Zhou C, Huang Y, Przedborski S: Oxidative stress in Parkinson’s disease: a mechanism of pathogenic and therapeutic significance. Ann N Y Acad Sci. 2008, 1147: 93-104. 10.1196/annals.1427.023.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Chinta SJ, Andersen JK: Redox imbalance in Parkinson’s disease. Biochim Biophys Acta. 2008, 1780 (11): 1362-1367. 10.1016/j.bbagen.2008.02.005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Copple IM, Goldring CE, Kitteringham NR, Park BK: The Nrf2-Keap1 defence pathway: role in protection against drug-induced toxicity. Toxicology. 2008, 246 (1): 24-33. 10.1016/j.tox.2007.10.029.

    Article  CAS  PubMed  Google Scholar 

  5. Itoh K, Wakabayashi N, Katoh Y, Ishii T, Igarashi K, Engel JD, Yamamoto M: Keap1 represses nuclear activation of antioxidant responsive elements by Nrf2 through binding to the amino-terminal Neh2 domain. Genes Dev. 1999, 13 (1): 76-86. 10.1101/gad.13.1.76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Cullinan SB, Gordan JD, Jin J, Harper JW, Diehl JA: The Keap1-BTB protein is an adaptor that bridges Nrf2 to a Cul3-based E3 ligase: oxidative stress sensing by a Cul3-Keap1 ligase. Mol Cell Biol. 2004, 24 (19): 8477-8486. 10.1128/MCB.24.19.8477-8486.2004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Furukawa M, Xiong Y: BTB protein Keap1 targets antioxidant transcription factor Nrf2 for ubiquitination by the Cullin 3-Roc1 ligase. Mol Cell Biol. 2005, 25 (1): 162-171. 10.1128/MCB.25.1.162-171.2005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Itoh K, Chiba T, Takahashi S, Ishii T, Igarashi K, Katoh Y, Oyake T, Hayashi N, Satoh K, Hatayama I, Yamamoto M, Nabeshima Y: An Nrf2/small Maf heterodimer mediates the induction of phase II detoxifying enzyme genes through antioxidant response elements. Biochem Biophys Res Commun. 1997, 236 (2): 313-322. 10.1006/bbrc.1997.6943.

    Article  CAS  PubMed  Google Scholar 

  9. Magesh S, Chen Y, Hu L: Small molecule modulators of Keap1-Nrf2-ARE pathway as potential preventive and therapeutic agents. Med Res Rev. 2012, 32 (4): 687-726. 10.1002/med.21257.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Kumar H, Koppula S, Kim IS, More SV, Kim BW, Choi DK: Nuclear factor erythroid 2-related factor 2 signaling in Parkinson disease: a promising multi therapeutic target against oxidative stress, neuroinflammation and cell death. CNS & neurological disorders drug targets. 2012, 11 (8): 1015-1029. 10.2174/1871527311211080012.

    Article  CAS  Google Scholar 

  11. Yang Y, Jiang S, Yan J, Li Y, Xin Z, Lin Y, Qu Y: An overview of the molecular mechanisms and novel roles of Nrf2 in neurodegenerative disorders. Cytokine Growth Factor Rev 2014, Epub ahead of print.,

  12. Cook AL, Vitale AM, Ravishankar S, Matigian N, Sutherland GT, Shan J, Sutharsan R, Perry C, Silburn PA, Mellick GD, Whitelaw ML, Wells CA, Mackay-Sim A, Wood SA: NRF2 activation restores disease related metabolic deficiencies in olfactory neurosphere-derived cells from patients with sporadic Parkinson’s disease. PLoS ONE. 2011, 6 (7): e21907-10.1371/journal.pone.0021907.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Meesarapee B, Thampithak A, Jaisin Y, Sanvarinda P, Suksamrarn A, Tuchinda P, Morales NP, Sanvarinda Y: Curcumin I mediates neuroprotective effect through attenuation of quinoprotein formation, p-p38 MAPK expression, and caspase-3 activation in 6-hydroxydopamine treated SH-SY5Y cells. Phytother Res. 2014, 28 (4): 611-616. 10.1002/ptr.5036.

    Article  CAS  PubMed  Google Scholar 

  14. Kaidery NA, Banerjee R, Yang L, Smirnova NA, Hushpulian DM, Liby KT, Williams CR, Yamamoto M, Kensler TW, Ratan RR, Sporn MB, Beal MF, Gazaryan IG, Thomas B: Targeting Nrf2-mediated gene transcription by extremely potent synthetic triterpenoids attenuate dopaminergic neurotoxicity in the MPTP mouse model of Parkinson’s disease. Antioxidants & redox signaling. 2013, 18 (2): 139-157. 10.1089/ars.2011.4491.

    Article  CAS  Google Scholar 

  15. Wild AC, Moinova HR, Mulcahy RT: Regulation of gamma-glutamylcysteine synthetase subunit gene expression by the transcription factor Nrf2. The Journal of biological chemistry. 1999, 274 (47): 33627-33636. 10.1074/jbc.274.47.33627.

    Article  CAS  PubMed  Google Scholar 

  16. Alam J, Stewart D, Touchard C, Boinapally S, Choi AM, Cook JL: Nrf2, a Cap‘n’Collar transcription factor, regulates induction of the heme oxygenase-1 gene. The Journal of biological chemistry. 1999, 274 (37): 26071-26078. 10.1074/jbc.274.37.26071.

    Article  CAS  PubMed  Google Scholar 

  17. Chen K, Gunter K, Maines MD: Neurons overexpressing heme oxygenase-1 resist oxidative stress-mediated cell death. J Neurochem. 2000, 75 (1): 304-313. 10.1046/j.1471-4159.2000.0750304.x.

    Article  CAS  PubMed  Google Scholar 

  18. Venugopal R, Jaiswal AK: Nrf1 and Nrf2 positively and c-Fos and Fra1 negatively regulate the human antioxidant response element-mediated expression of NAD(P)H:quinone oxidoreductase1 gene. Proc Natl Acad Sci U S A. 1996, 93 (25): 14960-14965. 10.1073/pnas.93.25.14960.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Segura-Aguilar J, Lind C: On the mechanism of the Mn3(+)-induced neurotoxicity of dopamine:prevention of quinone-derived oxygen toxicity by DT diaphorase and superoxide dismutase. Chem Biol Interact. 1989, 72 (3): 309-324. 10.1016/0009-2797(89)90006-9.

    Article  CAS  PubMed  Google Scholar 

  20. van Muiswinkel FL, de Vos RA, Bol JG, Andringa G, Jansen Steur EN, Ross D, Siegel D, Drukarch B: Expression of NAD(P)H:quinone oxidoreductase in the normal and Parkinsonian substantia nigra. Neurobiol Aging. 2004, 25 (9): 1253-1262. 10.1016/j.neurobiolaging.2003.12.010.

    Article  CAS  PubMed  Google Scholar 

  21. Ramsey CP, Glass CA, Montgomery MB, Lindl KA, Ritson GP, Chia LA, Hamilton RL, Chu CT, Jordan-Sciutto KL: Expression of Nrf2 in neurodegenerative diseases. J Neuropathol Exp Neurol. 2007, 66 (1): 75-85. 10.1097/nen.0b013e31802d6da9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Imaizumi Y, Okada Y, Akamatsu W, Koike M, Kuzumaki N, Hayakawa H, Nihira T, Kobayashi T, Ohyama M, Sato S, Takanashi M, Funayama M, Hirayama A, Soga T, Hishiki T, Suematsu M, Yagi T, Ito D, Kosakai A, Hayashi K, Shouji M, Nakanishi A, Suzuki N, Mizuno Y, Mizushima N, Amagai M, Uchiyama Y, Mochizuki H, Hattori N, Okano H: Mitochondrial dysfunction associated with increased oxidative stress and alpha-synuclein accumulation in PARK2 iPSC-derived neurons and postmortem brain tissue. Mol Brain. 2012, 5: 35-10.1186/1756-6606-5-35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. von Otter M, Landgren S, Nilsson S, Celojevic D, Bergstrom P, Hakansson A, Nissbrandt H, Drozdzik M, Bialecka M, Kurzawski M, Blennow K, Nilsson M, Hammarsten O, Zetterberg H: Association of Nrf2-encoding NFE2L2 haplotypes with Parkinson’s disease. BMC Med Genet. 2010, 11: 36-10.1186/1471-2350-11-36.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Marczak ED, Marzec J, Zeldin DC, Kleeberger SR, Brown NJ, Pretorius M, Lee CR: Polymorphisms in the transcription factor NRF2 and forearm vasodilator responses in humans. Pharmacogenet Genomics. 2012, 22 (8): 620-628. 10.1097/FPC.0b013e32835516e5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Marzec JM, Christie JD, Reddy SP, Jedlicka AE, Vuong H, Lanken PN, Aplenc R, Yamamoto T, Yamamoto M, Cho HY, Kleeberger SR: Functional polymorphisms in the transcription factor NRF2 in humans increase the risk of acute lung injury. Faseb J. 2007, 21 (9): 2237-2246. 10.1096/fj.06-7759com.

    Article  CAS  PubMed  Google Scholar 

  26. Bergman O, Hakansson A, Westberg L, Nordenstrom K, Carmine Belin A, Sydow O, Olson L, Holmberg B, Eriksson E, Nissbrandt H: PITX3 polymorphism is associated with early onset Parkinson’s disease. Neurobiol Aging. 2010, 31 (1): 114-117. 10.1016/j.neurobiolaging.2008.03.008.

    Article  CAS  PubMed  Google Scholar 

  27. Westerlund M, Belin AC, Anvret A, Hakansson A, Nissbrandt H, Lind C, Sydow O, Olson L, Galter D: Cerebellar alpha-synuclein levels are decreased in Parkinson’s disease and do not correlate with SNCA polymorphisms associated with disease in a Swedish material. Faseb J. 2008, 22 (10): 3509-3514. 10.1096/fj.08-110148.

    Article  CAS  PubMed  Google Scholar 

  28. Dick FD, De Palma G, Ahmadi A, Scott NW, Prescott GJ, Bennett J, Semple S, Dick S, Counsell C, Mozzoni P, Haites N, Wettinger SB, Mutti A, Otelea M, Seaton A, Söderkvist P, Felice A, Geoparkinson study group: Environmental risk factors for Parkinson’s disease and parkinsonism: the Geoparkinson study. Occup Environ Med. 2007, 64 (10): 666-672. 10.1136/oem.2006.027003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. International HapMap Consortium: The international HapMap project. Nature. 2003, 426 (6968): 789-796. 10.1038/nature02168.

    Article  Google Scholar 

  30. Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics (Oxford, England). 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.

    Article  CAS  Google Scholar 

  31. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.

    Article  CAS  PubMed  Google Scholar 

  32. Livak KJ: Allelic discrimination using fluorogenic probes and the 5′ nuclease assay. Genet Anal. 1999, 14 (5–6): 143-149. 10.1016/S1050-3862(98)00019-9.

    Article  CAS  PubMed  Google Scholar 

  33. Altman DG: Practical statistics for medical research. 1991, Chapman and Hall, New York

    Google Scholar 

  34. Excoffier L, Slatkin M: Maximum-likelihood estimation of molecular haplotype frequencies in a diploid population. Mol Biol Evol. 1995, 12 (5): 921-927.

    CAS  PubMed  Google Scholar 

  35. Chen YC, Wu YR, Wu YC, Lee-Chen GJ, Chen CM: Genetic analysis of NFE2L2 promoter variation in Taiwanese Parkinson’s disease. Parkinsonism Relat Disord. 2013, 19 (2): 247-250. 10.1016/j.parkreldis.2012.10.018.

    Article  PubMed  Google Scholar 

  36. Mizuta I, Satake W, Nakabayashi Y, Ito C, Suzuki S, Momose Y, Nagai Y, Oka A, Inoko H, Fukae J, Saito Y, Sawabe M, Murayama S, Yamamoto M, Hattori N, Murata M, Toda T: Multiple candidate gene analysis identifies alpha-synuclein as a susceptibility gene for sporadic Parkinson’s disease. Hum Mol Genet. 2006, 15 (7): 1151-1158. 10.1093/hmg/ddl030.

    Article  CAS  PubMed  Google Scholar 

  37. Maraganore DM, de Andrade M, Lesnick TG, Strain KJ, Farrer MJ, Rocca WA, Pant PV, Frazer KA, Cox DR, Ballinger DG: High-resolution whole-genome association study of Parkinson disease. Am J Hum Genet. 2005, 77 (5): 685-693. 10.1086/496902.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are indebted to the patients who participated in this study. The authors are grateful for the excellent technical assistance of Mrs. Mona Seibt Palmér. This work was supported by grants from the Swedish Research Council, the Knut and Alice Wallenberg Foundation, the Sahlgrenska University Hospital, the West Sweden RUN fundings, Edith Jacobsson Foundation, Axel Linder’s Foundation, the Göteborg Medical Society, the Swedish Medical Society, Swedish Brain Power, Stiftelsen för Gamla Tjänarinnor, the Swedish Parkinson Foundation, the Foundation for Parkinson Research at Linköping University (Stiftelsen för Parkinsonforskning), Sweden, Gun and Bertil Stohne’s Foundation, the Åhlén Foundation, the Alzheimer’s Foundation, Sweden, the Assar Gabrielsson Foundation, The Swedish Cancer Foundation, The Swedish Pain Foundation and the Herman and Lilly Schilling Foundation. The samples from Malta were from the Malta BioBank, a founding partner of the Eurobiobank. They were collected by the Laboratory of Molecular Genetics, University of Malta, as part of the Geoparkinson Study funded by the Fifth Framework programme of the European Commission, project number QLK4-CT-1999-01133. Funding sources played no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; or preparation, review, or approval of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Petra Bergström.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MVO: Research project – Organization, design; Laboratory work – Execution; Statistical analysis – Design and execution; Manuscript – Review and critique. PB: Statistical analysis – Execution; Manuscript – Writing. AQ: Laboratory work – Supplier of patient material; Manuscript – Review and critique. VDM: Laboratory work – Supplier of patient material; Manuscript – Review and critique. GA: Laboratory work – Supplier of patient material; Manuscript – Review and critique. PS: Laboratory work – Supplier of patient material; Manuscript – Review and critique. SBW: Laboratory work – Supplier of patient material; Manuscript – Review and critique. MD: Laboratory work – Supplier of patient material; Manuscript – Review and critique. MB: Laboratory work – Supplier of patient material; Manuscript – Review and critique. HN: Laboratory work –Supplier of patient material; Manuscript – Review and critique. CK: Laboratory work – Supplier of patient material; Manuscript – Review and critique. MN: Research project – Conception; Manuscript – Review and critique. OH: Research project – Conception; Manuscript – Review and critique. SN: Statistical analysis – Design, review and critique; Manuscript – Review and critique. HZ: Research project – Conception; Manuscript – Review and critique. All authors have read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

von Otter, M., Bergström, P., Quattrone, A. et al. Genetic associations of Nrf2-encoding NFE2L2 variants with Parkinson’s disease – a multicenter study. BMC Med Genet 15, 131 (2014). https://0-doi-org.brum.beds.ac.uk/10.1186/s12881-014-0131-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12881-014-0131-4

Keywords