Common and rare variants in HFE are not associated with Parkinson’s disease in Europeans ========================================================================================== * Prabhjyot Saini * Sara Bandres-Ciga * Jose Luis Alcantud * Clara Ruz * Ziv Gan-Or * for the International Parkinson’s disease Genomics Consortium ## Abstract A recent study suggested that the p.H63D variant in *HFE*, a gene involved in iron homeostasis, may modify α-synuclein pathology, the pathological hallmark of Parkinson’s disease (PD). If indeed this gene and specific variant are involved in PD, we expect to find differential distribution of *HFE* variants when comparing PD patients and controls. We analyzed genome-wide association study (GWAS) data from 14,671 PD patients and 17,667 controls and full sequencing data from additional 1,647 PD patients and 1,050 controls, using logistic regression models, and burden and Kernel tests. The *HFE* p.H63D variant was not associated with PD, nor did all the other common variants in the *HFE* locus. We did not find association of rare *HFE* variants with PD as well in all types of burden and Kernel tests. Our results do not support a role for *HFE* in PD. ## Introduction Parkinson’s disease (PD) is a complex neurodegenerative disorder, caused in most cases by multiple genetic and environmental factors, and the normal process of aging. Previous studies have suggested that iron homeostasis may have a role in PD; pathological and imaging studies showed that iron accumulation in the substantia nigra may contribute to PD pathogenesis (Ward et al., 2014). Multiple genes are involved in iron homeostasis, one of the most important of them being *HFE*, encoding the homeostatic iron regulator protein (also called hereditary hemochromatosis protein). This protein is a membrane protein thought to regulate iron absorption, and bi-allelic *HFE* mutations may lead to hereditary haemochromatosis, an iron storage disorder. Despite several negative genetic studies (Biasiotto et al., 2008; Duan et al., 2016; Guerreiro et al., 2006; Rhodes et al., 2014; Xia et al., 2015), a recent study suggested that the p.H63D variant in *HFE* may contribute to Parkinson’s disease pathogenesis by modifying α□synuclein pathology in cell models (Kim et al., 2020). The same group has previously suggested, using a mouse model with the homolog of the same variant (H67D), that *HFE* is a disease modifier in PD (Nixon et al., 2018). In order to examine whether this and other *HFE* variants are associated with PD, we performed a comprehensive analysis of common and rare *HFE* variants in large case-control cohorts of PD patients and controls of European origin. ## Methods Two cohorts were included in the current study: 1) a total of 14,671 PD patients and 17,667 controls collected through collaborators from the International PD Genomics Consortium (IPDGC), and 2) 1,647 PD patients and 1,050 controls from AMP-PD ([https://amp-pd.org/](https://amp-pd.org/)). Demographic details on these cohorts can be found in Supplementary Table S1A and S1B. We performed standard genome-wide association study (GWAS) quality control (QC) on the individual and variant level data as previously described (Nalls et al., 2019). We performed a similar QC process for the AMP-PD whole genome sequencing data, as detailed in the AMP-PD portal ([https://amp-pd.org/whole-genome-data](https://amp-pd.org/whole-genome-data)). Genotype data was extracted from both datasets for all variants within 100kb upstream and downstream of *HFE*, and was annotated using ANNOVAR (Wang et al., 2010). To examine whether common variants (allele frequency >0.01) in *HFE* are associated with PD, we performed logistic regression adjusted for age, sex, 10 principal components and sites of cohorts. To determine whether rare variants (allele frequency <0.01) we performed a set of burden and Kernel association tests included in the RVtests R package (Zhan et al., 2016). Bonferroni correction for multiple comparisons was applied, and all code used in the current study is available at the IPDGC GitHub at [https://github.com/ipdgc/IPDGC-Trainees/blob/master/](https://github.com/ipdgc/IPDGC-Trainees/blob/master/). The institutional review board (McGill University Health Center Research Ethics Board -MUHC REB) approved the study protocols. ## Results None of the common variants in the *HFE* region was associated with risk of PD (Figure 1). The specific *HFE* variant recently reported to be associated with decreased aggregation of α□synuclein, p.H63D, had similar allele frequencies in patients and controls (0.15 and 0.14, respectively), and was not associated with risk of PD (*p*=1 after correction for multiple comparisons). Another *HFE* variant previously reported to be associated with PD, p.C282Y, had similar frequencies in patients and controls as well (0.049 in both) and was not associated with PD (corrected *p*=1). Similarly, none of the burden and Kernel association tests we used to examine the association of rare *HFE* variants with PD yielded statistically significant results (Table 1). The list of all common and rare variants found in *HFE* is available in Supplementary Table S2 and S3, respectively. ![Figure1](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/26/2020.10.22.20217190/F1.medium.gif) [Figure1](http://medrxiv.org/content/early/2020/10/26/2020.10.22.20217190/F1) View this table: [Table 1.](http://medrxiv.org/content/early/2020/10/26/2020.10.22.20217190/T1) Table 1. Burden tests for *HFE* in the AMP-PD cohorts ## Discussion In the current study we performed a thorough genetic analysis of *HFE* in large case-control cohorts of PD patients and controls. Our results do not support a role for common or rare *HFE* variants in PD. Previous smaller genetic association studies focusing on *HFE* reported contradicting results, as some studies reported associations of different *HFE* variants with PD, while other studies reported lack of association (Biasiotto et al., 2008; Borie et al., 2002; Buchanan et al., 2002; Duan et al., 2016; Guerreiro et al., 2006; Halling et al., 2008; Rhodes et al., 2014; Xia et al., 2015). Specifically for p.H63D, all the previous studies, including two meta-analyses (Duan et al., 2016; Xia et al., 2015), reported lack of association. Together with our study, there is no support from human genetic studies in PD patients for the involvement of *HFE* p.H63D in PD. Furthermore, there is no convincing evidence to suggest an association between other *HFE* variants and risk of PD. We cannot rule out that a very small effect does exist, but since it was not detected by the recent GWAS with over 50,000 PD patients and proxy-cases (Nalls et al., 2019), it is likely to be too small to be clinically meaningful, if at all. The main limitation of the current study is that it includes only participants of European descent, and since PD is a complex disorder, it is possible that in other populations *HFE* variants may still contribute to the disease when found on another genetic background. In a recent GWAS of Asian population, however, no GWAS-significant association has been reported for *HFE* (Foo et al., 2020). Despite the lack of association in the current study, the role of iron and iron homeostasis in PD should be further studied, whether by genetic studies or by biochemical, functional, clinical and pathological studies. We recommend that genetic data from the IPDGC and other available studies should be consulted before studying specific genetic variants in different models. ## Supporting information Supplementary Tables S1-S3 [[supplements/217190_file04.xlsx]](pending:yes) ## Data Availability Data used for the analysis can be found in the supplementary tables. The code used for the analysis is available at https://github.com/ipdgc/IPDGC-Trainees ## Conflict of interests ZGO has received consulting fees from Lysosomal Therapeutics Inc., Idorsia, Prevail Therapeutics, Denali, Ono Therapeutics, Neuron23, Handl Therapeutics, Deerfield and Inception Sciences (now Ventus). None of these companies were involved in any parts of preparing, drafting and publishing this study. Other authors have no additional disclosures to report. ## Acknowledgements We thank the participants for contributing to the study. We would like to also thank all members of the International Parkinson Disease Genomics Consortium (IPDGC). For a complete overview of members, acknowledgements and funding, please see [http://pdgenetics.org/partners](http://pdgenetics.org/partners). This research was supported in part by the Intramural Research Program of the NIH, National institute on Aging. This work was also supported in part by grants from the Michael J. Fox Foundation, the Canadian Consortium on Neurodegeneration in Aging (CCNA), and the Canada First Research Excellence Fund (CFREF), awarded to McGill University for the Healthy Brains for Healthy Lives initiative (HBHL). ZGO is supported by the Fonds de recherche du Québec - Santé (FRQS) Chercheurs-boursiers award, in collaboration with Parkinson Quebec, and by the Young Investigator Award by Parkinson Canada. * Received October 22, 2020. * Revision received October 22, 2020. * Accepted October 26, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. Biasiotto, G., Goldwurm, S., Finazzi, D., Tunesi, S., Zecchinelli, A., Sironi, F., Pezzoli, G., Arosio, P., 2008. HFE gene mutations in a population of Italian Parkinson’s disease patients. Parkinsonism Relat Disord 14(5), 426–430. 2. Borie, C., Gasparini, F., Verpillat, P., Bonnet, A.M., Agid, Y., Hetet, G., Brice, A., Durr, A., Grandchamp, B., French Parkinson’s disease genetic study, g., 2002. Association study between iron-related genes polymorphisms and Parkinson’s disease. J Neurol 249(7), 801–804. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00415-002-0704-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12140659&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177159700001&link_type=ISI) 3. Buchanan, D.D., Silburn, P.A., Chalk, J.B., Le Couteur, D.G., Mellick, G.D., 2002. The Cys282Tyr polymorphism in the HFE gene in Australian Parkinson’s disease patients. Neurosci Lett 327(2), 91–94. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0304-3940(02)00398-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12098643&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom) 4. Duan, C., Wang, M., Zhang, Y., Wei, X., Huang, Y., Zhang, H., Cheng, L., Gai, Z., 2016. C282Y and H63D Polymorphisms in Hemochromatosis Gene and Risk of Parkinson’s Disease: A Meta-Analysis. Am J Alzheimers Dis Other Demen 31(3), 201–207. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1533317515602220&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26340960&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom) 5. Foo, J.N., Chew, E.G.Y., Chung, S.J., Peng, R., Blauwendraat, C., Nalls, M.A., Mok, K.Y., Satake, W., Toda, T., Chao, Y., Tan, L.C.S., Tandiono, M., Lian, M.M., Ng, E.Y., Prakash, K.M., Au, W.L., Meah, W.Y., Mok, S.Q., Annuar, A.A., Chan, A.Y.Y., Chen, L., Chen, Y., Jeon, B.S., Jiang, L., Lim, J.L., Lin, J.J., Liu, C., Mao, C., Mok, V., Pei, Z., Shang, H.F., Shi, C.H., Song, K., Tan, A.H., Wu, Y.R., Xu, Y.M., Xu, R., Yan, Y., Yang, J., Zhang, B., Koh, W.P., Lim, S.Y., Khor, C.C., Liu, J., Tan, E.K., 2020. Identification of Risk Loci for Parkinson Disease in Asians and Comparison of Risk Between Asians and Europeans: A Genome-Wide Association Study. JAMA Neurol 77(6), 746–754. 6. Guerreiro, R.J., Bras, J.M., Santana, I., Januario, C., Santiago, B., Morgadinho, A.S., Ribeiro, M.H., Hardy, J., Singleton, A., Oliveira, C., 2006. Association of HFE common mutations with Parkinson’s disease, Alzheimer’s disease and mild cognitive impairment in a Portuguese cohort. BMC Neurol 6, 24. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2377-6-24&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16824219&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom) 7. Halling, J., Petersen, M.S., Grandjean, P., Weihe, P., Brosen, K., 2008. Genetic predisposition to Parkinson’s disease: CYP2D6 and HFE in the Faroe Islands. Pharmacogenet Genomics 18(3), 209–212. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/FPC.0b013e3282f5106e&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18300942&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom) 8. Kim, Y., Stahl, M.C., Huang, X., Connor, J.R., 2020. H63D variant of the homeostatic iron regulator (HFE) gene alters alpha-synuclein expression, aggregation, and toxicity. J Neurochem. 9. Nalls, M.A., Blauwendraat, C., Vallerga, C.L., Heilbron, K., Bandres-Ciga, S., Chang, D., Tan, M., Kia, D.A., Noyce, A.J., Xue, A., Bras, J., Young, E., von Coelln, R., Simon-Sanchez, J., Schulte, C., Sharma, M., Krohn, L., Pihlstrom, L., Siitonen, A., Iwaki, H., Leonard, H., Faghri, F., Gibbs, J.R., Hernandez, D.G., Scholz, S.W., Botia, J.A., Martinez, M., Corvol, J.C., Lesage, S., Jankovic, J., Shulman, L.M., Sutherland, M., Tienari, P., Majamaa, K., Toft, M., Andreassen, O.A., Bangale, T., Brice, A., Yang, J., Gan-Or, Z., Gasser, T., Heutink, P., Shulman, J.M., Wood, N.W., Hinds, D.A., Hardy, J.A., Morris, H.R., Gratten, J., Visscher, P.M., Graham, R.R., Singleton, A.B., andMe Research, T., System Genomics of Parkinson’s Disease, C., International Parkinson’s Disease Genomics, C., 2019. Identification of novel risk loci, causal insights, and heritable risk for Parkinson’s disease: a meta-analysis of genome-wide association studies. Lancet Neurol 18(12), 1091–1102. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1474-4422(19)30320-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom) 10. Nixon, A.M., Meadowcroft, M.D., Neely, E.B., Snyder, A.M., Purnell, C.J., Wright, J., Lamendella, R., Nandar, W., Huang, X., Connor, J.R., 2018. HFE Genotype Restricts the Response to Paraquat in a Mouse Model of Neurotoxicity. J Neurochem 145(4), 299–311. 11. Rhodes, S.L., Buchanan, D.D., Ahmed, I., Taylor, K.D., Loriot, M.A., Sinsheimer, J.S., Bronstein, J.M., Elbaz, A., Mellick, G.D., Rotter, J.I., Ritz, B., 2014. Pooled analysis of iron-related genes in Parkinson’s disease: association with transferrin. Neurobiol Dis 62, 172–178. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.nbd.2013.09.019&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24121126&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom) 12. Wang, K., Li, M., Hakonarson, H., 2010. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 38(16), e164. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkq603&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20601685&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom) 13. Ward, R.J., Zucca, F.A., Duyn, J.H., Crichton, R.R., Zecca, L., 2014. The role of iron in brain ageing and neurodegenerative disorders. Lancet Neurol 13(10), 1045–1060. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1474-4422(14)70117-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25231526&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom) 14. Xia, J., Xu, H., Jiang, H., Xie, J., 2015. The association between the C282Y and H63D polymorphisms of HFE gene and the risk of Parkinson’s disease: A meta-analysis. Neurosci Lett 595, 99–103. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neulet.2015.04.010&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25863172&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom) 15. Zhan, X., Hu, Y., Li, B., Abecasis, G.R., Liu, D.J., 2016. RVTESTS: an efficient and comprehensive tool for rare variant association analysis using sequence data. Bioinformatics 32(9), 1423–1426. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw079&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27153000&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F26%2F2020.10.22.20217190.atom)