Genetic Analysis of EGLN1 C127S Variant in Taiwanese Parkinson’s Disease

Abstract

Parkinson’s disease (PD) is a neurodegenerative disorder related to nigrostriatal dopaminergic neuron degeneration and iron accumulation. As a cellular oxygen sensor, prolyl hydroxylase domain containing protein 2 (PHD2, encoded by egl-9 family hypoxia inducible factor 1, EGLN1) modifies hypoxia-inducible factor alpha (HIF-α) protein for proteasomal destruction under normoxic condition. In addition, 2-oxoglutarate- (OG-) dependent dioxygenase activity of PHD2 is involved in the oxygen and iron regulation of iron-responsive element binding protein 2 (IRP2) stability. Previously increased expression of EGLN1 was found in the substantia nigra of the parkinsonian brain. We investigated the possible role of c.380 G > C (p.C127S) of EGLN1 gene in Taiwanese patients with PD. 479 patients and 435 healthy controls were recruited. Polymerase chain reaction and BsmAI restriction enzyme analysis were applied for analysis. An association between CC genotype and reduced PD risk in the recessive model (CC vs. GG + GC) was found. Our study provides a link between EGLN1 c.380 G > C SNP and the development of PD.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Parkinson’s disease (PD) is a neurodegenerative disorder characterized by α-synuclein aggregation in the brain and loss of dopaminergic neurons in the substantia nigra (SN). Although the disease etiology remains to be clarified, proteins encoded by PD-causing genes such as SNCA, Parkin, DJ1, PINK1, and LRRK2 are involved in the oxidative stress pathway that contributes to neurodegeneration of PD [1]. Among the sources and mechanisms for the generation of oxidative stress, dysfunction of iron homeostasis and iron accumulation in the SN may be involved in the development of PD [2]. Previous study by Grunblatt et al. found increased gene expression of EGLN1 (egl-9 family hypoxia inducible factor 1) in the SN of parkinsonian brains using oligonucleotide microarray technique [3]. The EGLN1 gene encodes proline hydroxylase domain 2 (PHD2), which is an oxygen/iron sensor that belongs to the 2-oxoglutarate- (2-OG-) dependent dioxygenase superfamily. The upregulation of EGLN1 promotes the degradation of hypoxic-inducible factor (HIF), an iron homeostasis-related protein [4]. Furthermore, evidence has shown that the iron regulatory protein 2 (IRP2), another key iron regulatory protein in mammals, was stabilized by a 2-OG-dependent dioxygenase specific inhibitor, dimethyloxalyglycine [5, 6], suggesting that EGLN1 may downregulate IRP2. Intriguingly, in IRP2 knock-out mice, significant iron accumulation was found in the brain with neurodegeneration and adult onset ataxia, tremor, and bradykinesia [7]. Moreover, a high frequency of linked c.[12C > G; 380 G >C] or p.[Asp4Glu; Cys127Ser] variation in Tibetans [8] and Sherpas [9] was found in recent studies, which results in an adaptive function of the EGLN1 for high altitude [8], although the finding was not replicated in the study of highland Andreans [10]. Thus, adaptation mechanisms involving EGLN1 may or may not be similar among different populations. In this study, we investigated the possible role of c.380G > C (p.C127S, rs12097901) of the EGLN gene in Taiwanese patients with PD.

2. Materials and Methods

2.1. Ethics Statement

This study was performed under a protocol approved by the institutional review boards of Chang Gung Memorial Hospital (ethical license no: 95-0531B), and all examinations were performed after obtaining written informed consent.

2.2. Study Population

Patients diagnosed with PD were recruited from the neurology clinic of Chang-Gung Memorial Hospital. The diagnosis of PD was based on the UK PD Society Brain Bank clinical diagnostic criteria [11]. Unrelated healthy adult volunteers matched for age, gender, ethnic origin, and area of residences were recruited as controls.

2.3. Genetic Analysis

DNA was extracted from leukocytes by using the standard protocols. EGLN1 gene polymorphism rs12097901 (c.380G > C, a Cys-to-Ser variant at position 127) was determined using the polymerase chain reaction and BsmAI restriction enzyme analysis. Sequences of primers and polymorphic change in the enzyme recognition site were as follows: 5′-CCACACCAGCATTCCGGC-3′ (forward primer), 5′-TTGTTCATGCACGGCACGAT-3′ (reverse primer), and GTcTC (BsmAI). The lowercase letter in the BsmAI recognition site indicates the polymorphic site. The PCR was carried out in a reaction containing 100 ng genomic DNA, 0.2 μM of each primer, 0.2 mM dNTPs, 1.5 mM MgCl₂, 0.5 U Taq polymerase, and 10% DMSO. Thermal cycling conditions consisted of 94°C for 6 minutes for initial denaturation, followed by 35 cycles of 94°C for 30 seconds, 52°C for 30 seconds, and 72°C for 30 seconds. A final extension step of 72°C for 10 minutes was followed by a 4°C hold cycle. The amplified 401 bp PCR fragments, consisting of 75.3% CG content, were digested with the BsmAI (New England Biolabs) and separated on a 2.0% agarose gel (allele G, 401-bp fragment; allele C, 233- and 168-bp fragments).

2.4. Statistical Analysis

The chi-square test was used to compare the frequency of the allele and genotypes in both cases and controls. Logistic regression analysis was performed to verify the interaction between age, sex, and the risk of PD. The genotypes of the subjects followed the Hardy–Weinberg equilibrium.

3. Results

We recruited 914 subjects, including 479 patients with PD (45.4% females) and 435 healthy controls (48.3% female). Only one proband with familial PD in the same family was recruited. The mean age at onset (AAO) of PD symptoms was 62.9 ± 10.7 years (range: 19–86), and the mean age of recruitment of 435 controls was 59.2 ± 12.7 years (range: 24–89). The frequencies of EGLN1 c.380 G > C (p.C127S) alleles were similar in both PD patients and controls. However, the CC genotype was notably lower in the PD patient group ($p=0.056$), and the results of the recessive model applied in genotype analysis (CC versus GG + GC) showed an association between CC and reduced PD risk (OR = 0.59, $p=0.002$) (Table 1). Stratification by age at onset (<50 and ≧ 50 years) and sex also did not show differences in the minor C allele frequency in each case control cohort (data not shown).

Table 1

Distributions and association of EGLN1 p.C127S (TGT > TCT) polymorphism in patients and controls.

	PD (%)	NC (%)	OR (95% CI)	$p$
	(n = 479)	(n = 435)
G	56.6	53.2	1.00
C	43.4	46.8	0.87 (0.72–1.05)	0.150
GG	28.8	30.3	1.00
GC	55.5	45.7	1.27 (0.95–1.72)	0.110
CC	15.7	23.9	0.69 (0.47–1.01)	0.056
Recessive model
GC + GG	84.3	76.1	1.00
CC	15.7	23.9	0.59 (0.42–0.82)	0.002

4. Discussion

The development of idiopathic PD is considered to be related with a complicated interplay of genetic and environmental factors affecting numerous fundamental cellular processes [12]. Among the environmental factors involved, hypoxic challenge can trigger neuronal cell death to exacerbate disease progression [13]. In addition, dysfunction in brain iron metabolism, specifically, iron accumulation in the SN, has been implicated in the pathogenesis of PD through increasing oxidative stress, impairing the ubiquitin proteasome system, and α-synuclein aggregation [2]. In this study, we studied the polymorphism of rs12097901 (c.380 G > C; p.C127S) of EGLN1 gene in patients with PD and healthy controls. We found an association between CC and reduced PD risk in the recessive model of gene analysis (CC versus GG + GC). The role of this genetic variation in PD is still unknown but may be associated with HIF-related pathways. The product of EGLN1 gene, PHD2, contains 426 residues, with the conserved catalytic domain located in the C-terminal region (∼181–426) [14]. The upregulation of EGLN1 promotes the degradation of HIF, and 2-OG-dependent-dioxygenases inhibitor stabilized IRP2 [5, 6]. We used protein structure prediction server [15] to confirm that this variant (p.C127S) did not change the protein structure at 3D level (Figure 1). Additionally, the pathogenicity of p.C127S was presumed to be benign by using SIFT and PolyPhen predictors. Although we did not check if this variant affects the gene expression level, kinetic analysis and functional experiments have shown that p.C127S variant has a lower Michaelis constant (Km) value for cosubstrate oxygen, suggesting that it promotes HIF degradation under hypoxic environment, hence blunting the hypoxic response in the literature [8]. The functional relevance of p.C127S (rs2790859) (an intergenic SNP downstream of EGLN1) alleles in an acute hypobaric hypoxia experiment was indicated as this variant was associated with percutaneous arterial oxygen saturation variations (SpO₂) and SpO₂ latency in a Japanese cohort [16]. Furthermore, loss of DJ-1, an autosomal recessive gene for familial PD, was shown to provoke hypoxic condition and increase reactive oxygen species production in human neuroblastoma cell, which in turn stabilized the HIF-1α. The posttranscriptional stabilization of HIF-1α was required to downregulate one of the receptors for the glial cell line-derived neurotrophic factor, ret proto-oncogene (RET). RET is strongly expressed in the dopaminergic (DA) neurons, and loss of RET expression in DA neurons generates PD-like dysfunctions in mouse models, thus indicating its crucial role in mesencephalic DA cell survival. Since RET is prosurvival, this study implicates that increased HIF-1α-associated response may contribute to neurodegeneration in PD [17].

[figure omitted; refer to PDF]

IRP2 is a posttranscriptional regulatory protein that binds to the iron-responsive element (IREs) in the mRNA of various iron regulatory proteins to facilitate iron uptake and prohibit iron storage during an iron deficient state [18]. Studies have suggested that 2-OG-dependent dioxygenase inhibitor and hypoxia or hypoxia mimetics stabilize IRP2 [5, 6], and N-terminally truncated form of PHD2 (residues 181-426) exhibits high affinity for, and copurifies with the endogenous level of, its Fe²⁺ cofactor and 2-oxoglutarate cosubstrate [19], linking 2-OG-dependent dioxygenase activity to IRP2 degradation. Therefore, based on our study result, it is speculated that 2-OG-dependent dioxygenase activity of PHD2 may be influenced by p.C127S to affect brain IRP2 and alters iron metabolism, contributing to the observed reduced PD risk. However, whether p.C127S affects the PHD2 activity is not known. Future studies in the functional consequence of p.C127S on the 2-OG-dependent dioxygenase activity and IRP2 degradation are necessary to further clarify how p.C127S affects PD risk. In conclusion, our study provides evidence that the protective C127S variant of EGLN1 is associated with reduced risk of PD. The limitation of our study is that only one SNP in EGLN1 was analyzed and that future studies analyzing other possible genetic variations should be warranted for further investigation.

Authors’ Contributions

Y. -R. W., G. -J. L. -C., and C. -M. C. were involved in conceptualization. Y. -C. C and C. -Y. C. carried out the methodology. Y. -R. W., C. -M. C.; Y. -C. C provided resources. H. -L. C. wrote the original draft. Y. -R. W., G. -J. L. -C., and C. -M. C. were responsible for reviewing and editing the manuscript. Y. -R. W carried out funding acquisition.

Acknowledgments

This research was funded by the National Science Council, Executive Yuan, Taiwan (MOST 95-2314-B-182A-061- and 107-2314-B-182A-045-MY2). The authors also thank Dr. Wen-Lang Fan, the assistant research fellow at Genomic Medicine Core Laboratory of Chang Gung Memorial Hospital, Linkou, Taiwan, for the protein structure prediction.

References

[1] V. Dias, E. Junn, M. M. Mouradian, "The role of oxidative stress in Parkinson’s disease," Journal of Parkinson’s Disease, vol. 3 no. 4, pp. 461-491, DOI: 10.3233/JPD-130230, 2013.

[2] H. Jiang, J. Wang, J. Rogers, J. Xie, "Brain iron metabolism dysfunction in Parkinson’s disease," Molecular Neurobiology, vol. 54 no. 4, pp. 3078-3101, DOI: 10.1007/s12035-016-9879-1, 2017.

[3] E. Grunblatt, S. Mandel, J. Jacob-Hirsch, "Gene expression profiling of parkinsonian substantia nigra pars compacta; alterations in ubiquitin-proteasome, heat shock protein, iron and oxidative stress regulated proteins, cell adhesion/cellular matrix and vesicle trafficking genes," Journal of Neural Transmission, vol. 111, pp. 1543-1573, DOI: 10.1007/s00702-004-0212-1, 2004.

[4] E. Berra, E. Benizri, A. Ginouves, V. Volmat, D. Roux, J. Pouyssegur, "HIF prolyl-hydroxylase 2 is the key oxygen sensor setting low steady-state levels of HIF-1 in normoxia," The EMBO Journal, vol. 22 no. 16, pp. 4082-4090, DOI: 10.1093/emboj/cdg392, 2003.

[5] E. S. Hanson, M. L. Rawlins, E. A. Leibold, "Oxygen and iron regulation of iron regulatory protein 2," Journal of Biological Chemistry, vol. 278 no. 41, pp. 40337-40342, DOI: 10.1074/jbc.M302798200, 2003.

[6] J. Wang, G. Chen, M. Muckenthaler, B. Galy, M. W. Hentze, K. Pantopoulos, "Iron-mediated degradation of IRP2, an unexpected pathway involving a 2-oxoglutarate-dependent oxygenase activity," Molecular and Cellular Biology, vol. 24 no. 3, pp. 954-965, DOI: 10.1128/mcb.24.3.954-965.2004, 2004.

[7] T. LaVaute, S. Smith, S. Cooperman, "Targeted deletion of the gene encoding iron regulatory protein-2 causes misregulation of iron metabolism and neurodegenerative disease in mice," Nature Genetics, vol. 27 no. 2, pp. 209-214, DOI: 10.1038/84859, 2001.

[8] F. R. Lorenzo, C. Huff, M. Myllymäki, "A genetic mechanism for Tibetan high-altitude adaptation," Nature Genetics, vol. 46 no. 9, pp. 951-956, DOI: 10.1038/ng.3067, 2014.

[9] S. Bhandari, X. Zhang, C. Cui, "Sherpas share genetic variations with Tibetans for high-altitude adaptation," Molecular Genetics & Genomic Medicine, vol. 5 no. 1, pp. 76-84, DOI: 10.1002/mgg3.264, 2017.

[10] E. C. Heinrich, L. Wu, E. S. Lawrence, "Genetic variants at the EGLN1 locus associated with high‐altitude adaptation in Tibetans are absent or found at low frequency in highland Andeans," Annals of Human Genetics, vol. 83 no. 3, pp. 171-176, DOI: 10.1111/ahg.12299, 2019.

[11] A. J. Hughes, S. E. Daniel, L. Kilford, A. J. Lees, "Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: a clinico-pathological study of 100 cases," Journal of Neurology, Neurosurgery & Psychiatry, vol. 55 no. 3, pp. 181-184, DOI: 10.1136/jnnp.55.3.181, 1992.

[12] L. V. Kalia, A. E. Lang, "Parkinson’s disease," The Lancet, vol. 386 no. 9996, pp. 896-912, DOI: 10.1016/S0140-6736(14)61393-3, 2015.

[13] O. O. Ogunshola, X. Antoniou, "Contribution of hypoxia to Alzheimer’s disease: is HIF-1 α a mediator of neurodegeneration?," Cellular and Molecular Life Sciences, vol. 66 no. 22, pp. 3555-3563, DOI: 10.1007/s00018-009-0141-0, 2009.

[14] R. K. Bruick, S. L. McKnight, "A conserved family of prolyl-4-hydroxylases that modify HIF," Science, vol. 294 no. 5545, pp. 1337-1340, DOI: 10.1126/science.1066373, 2001.

[15] C.-C. Chen, J.-K. Hwang, J.-M. Yang, "(PS)2-v2: template-based protein structure prediction server," BMC Bioinformatics, vol. 10 no. 1,DOI: 10.1186/1471-2105-10-366, 2009.

[16] Y. Yasukochi, T. Nishimura, M. Motoi, S. Watanuki, "Association of EGLN1 genetic polymorphisms with SpO2 responses to acute hypobaric hypoxia in a Japanese cohort," Journal of Physiological Anthropology, vol. 37 no. 1,DOI: 10.1186/s40101-018-0169-7, 2018.

[17] R. Foti, S. Zucchelli, M. Biagioli, "Parkinson disease-associated DJ-1 is required for the expression of the glial cell line-derived neurotrophic factor receptor RET in human neuroblastoma cells," Journal of Biological Chemistry, vol. 285 no. 24, pp. 18565-18574, DOI: 10.1074/jbc.M109.088294, 2010.

[18] O. Weinreb, S. Mandel, M. B. H. Youdim, T. Amit, "Targeting dysregulation of brain iron homeostasis in Parkinson’s disease by iron chelators," Free Radical Biology and Medicine, vol. 62, pp. 52-64, DOI: 10.1016/j.freeradbiomed.2013.01.017, 2013.

[19] L. A. McNeill, E. Flashman, M. R. G. Buck, "Hypoxia-inducible factor prolyl hydroxylase 2 has a high affinity for ferrous iron and 2-oxoglutarate," Molecular BioSystems, vol. 1 no. 4, pp. 321-324, DOI: 10.1039/b511249b, 2005.