Targeted sequencing of Parkinson’s disease loci genes highlights SYT11, FGF20 and other associations
======================================================================================================

* Uladzislau Rudakou
* Eric Yu
* Lynne Krohn
* Jennifer A. Ruskey
* Yves Dauvilliers
* Dan Spiegelman
* Lior Greenbaum
* Stanley Fahn
* Cheryl H. Waters
* Nicolas Dupré
* Guy A. Rouleau
* Sharon Hassin-Baer
* Edward A. Fon
* Roy N. Alcalay
* Ziv Gan-Or

## Abstract

Genome-wide association studies (GWAS) have identified numerous loci associated with Parkinson’s disease. The specific genes and variants that drive the associations within the vast majority of these loci are unknown. We aimed to perform a comprehensive analysis of selected genes to determine the potential role of rare and common genetic variants within these loci. We fully sequenced 32 genes from 25 loci previously associated with Parkinson’s disease in 2,657 patients and 3,647 controls from three cohorts. Capture was done using molecular inversion probes targeting the exons, exon-intron boundaries and untranslated regions (UTRs) of the genes of interest, followed by sequencing. Quality control was performed to include only high-quality variants. We examined the role of rare variants (minor allele frequency < 0.01) using optimized sequence Kernel association tests (SKAT-O). The association of common variants was estimated using regression models adjusted for age, sex and ethnicity as required in each cohort, followed by a meta-analysis. After Bonferroni correction, we identified a burden of rare variants in *SYT11, FGF20* and *GCH1* associated with Parkinson’s disease. Nominal associations were identified in 21 additional genes. Previous reports suggested that the *SYT11* GWAS association is driven by variants in the nearby *GBA* gene. However, the association of *SYT11* was mainly driven by a rare 3’ UTR variant (rs945006601) and was independent of *GBA* variants (*p* = 5.23E-05 after exclusion of all *GBA* variant carriers). The association of *FGF20* was driven by a rare 5’ UTR variant (rs1034608171) located in the promoter region. The previously reported association of *GCH1* with Parkinson’s Disease is driven by rare nonsynonymous variants, some of which are known to cause dopamine-responsive dystonia. We also identified two *LRRK2* variants, p.Arg793Met and p.Gln1353Lys, in ten and eight controls, respectively, but not in patients. We identified common variants associated with Parkinson’s disease in *MAPT, TMEM175, BST1*, *SNCA* and *GPNMB* which are all in strong linkage disequilibrium (LD) with known GWAS hits in their respective loci. A common coding *PM20D1* variant, p.Ile149Val, was nominally associated with reduced risk of Parkinson’s disease (OR 0.73, 95% CI 0.60–0.89, *p* = 1.161E-03). This variant is not in LD with the top GWAS hits within this locus and may represent a novel association. These results further demonstrate the importance of fine mapping of GWAS loci, and suggest that *SYT11, FGF20*, and potentially *PM20D1, BST1* and *GPNMB* should be considered for future studies as possible Parkinson’s disease-related genes.

Keywords
*   Parkinson’s disease
*   genetic fine mapping
*   risk loci
*   association study
*   sequencing

## Introduction

Genome-wide association studies (GWASs) are an important tool for identifying genetic associations with complex disorders such as Parkinson’s Disease (Nalls et al., 2014; Chang et al., 2017; Nalls et al., 2019). GWASs examine the association of multiple common variants across the genome with specific traits. Typically, each associated locus includes multiple genes, and very often GWAS alone cannot identify the specific gene or variant within each locus that drives the association. Furthermore, GWASs generally overlook rare variants, while accumulation of rare variants on a specific common allele may drive some of these associations. Even if not located on a specific common allele, rare variants in GWAS genes may contribute to disease. For example, rare variants in *SNCA* (Polymeropoulos et al., 1997), *GBA* (Sidransky et al., 2009; Ran et al., 2016; Emelyanov et al., 2018; Amaral et al., 2019)*, LRRK2* (Paisan-Ruiz et al., 2004; Khan et al., 2005), *VPS13C* (Lesage et al., 2016; Jansen et al., 2017; Darvish et al., 2018; Schormair et al., 2018; Rudakou et al., 2020), *GCH1* (Mencacci et al., 2014; Guella et al., 2015; Lewthwaite et al., 2015; Xu et al., 2017; Yoshino et al., 2018; Rudakou et al., 2019), all GWAS loci genes, also cause or increase the risk of Parkinson’s disease.

Recently, a large-scale GWAS has identified 78 loci associated with risk of Parkinson’s disease (Nalls et al., 2019), yet only for a few of them we have identified the specific genes and variants that drive the association. In some of these loci, there are multiple hits, suggesting that more than one variant, or perhaps more than one gene within the locus, can be associated with Parkinson’s disease. Fine-mapping of GWAS loci may help identifying these genes and variants, as was done for *SNCA* (Soldner et al., 2016; Krohn et al., 2019b), *VPS13C* (Lesage et al., 2016; Jansen et al., 2017; Darvish et al., 2018; Schormair et al., 2018; Rudakou et al., 2020) and *TMEM175* (Krohn et al., 2020). In the *SNCA* region, fine mapping has identified independent 5’ and 3’ variants that drive the associations with different synucleinopathies (Krohn et al., 2019b). Similarly, a potentially protective haplotype in *VPS13C* (Rudakou et al., 2020) and common nonsynonymous variants in *TMEM175* have been associated with Parkinson’s Disease through fine mapping, and in the case of *TMEM175* also by functional studies (Jinn et al., 2019; Krohn et al., 2020).

In the current study, we fully sequenced 32 genes from 25 different Parkinson’s disease GWAS loci in a total of 2,657 cases and 3,647 controls. We examined the association of common and rare variants with Parkinson’s Disease, including the role of rare bi-allelic variants.

## Materials and Methods

### Study Population

The study population included, consecutively recruited, unrelated patients with Parkinson’s Disease (n = 2,657) and controls (n = 3,647) from three cohorts, collected at McGill University (Quebec, Canada and Montpellier, France), Columbia University (New York, NY, USA), and the Sheba Medical Center (Israel). The McGill cohort was recruited in Quebec, Canada (in part with the assistance of the Quebec Parkinson’s Network(Gan-Or et al., 2020)) and in Montpellier, France. It includes 1,026 Parkinson’s Disease patients (mean age 61.5±11.3, 62.6% male) and 2,588 controls (mean age 54.5±14.2, 47.4% male), all unrelated of French-Canadian/French ancestry. The Columbia cohort includes 1,026 Parkinson’s Disease patients (mean age 59.4±11.7, 64.2% male) and 525 controls (mean age 64.2±10.1, 37.1% male). This cohort has been previously described in detail, and is mainly composed of participants with European descent, but 21.2% of Parkinson’s Disease patients and 18.7% controls are of Ashkenazi Jewish (AJ) origin (Alcalay et al., 2016). The Sheba cohort included 605 Parkinson’s Disease patients (mean age 60.7±11.9, 62.3% male) and 534 controls (mean age 34.0±7.0, 55.8% male), all of full AJ origin; previously described in more details (Ruskey et al., 2019). Movement disorder specialists diagnosed Parkinson’s Disease with UK brain bank criteria (Hughes et al., 1992) or the MDS clinical diagnostic criteria (Postuma et al., 2015).

As detailed in the Statistical Analyses subsection, due to the differences in age and sex, statistical analyses were adjusted and included age and sex as co-variates. To account for ethnical heterogeneity of the Columbia cohort, an ethnicity covariate was also introduced for the analysis of this cohort (GWAS data was not available, therefore the reported ethnicity was used and not principal components). All three cohorts were sequenced in the same lab (McGill University), using the same protocol.

### Standard Protocol Approvals, Registrations, and Patient Consents

The institutional review boards of McGill University, Columbia University and Sheba Medical Center had approved the study protocols and informed consent was obtained from all individual participants before entering the study.

### Target genes

A recent Parkinson’s Disease GWAS identified 78 loci, however the current study was designed and performed before its publication and the target genes were taken from earlier literature (Nalls et al., 2014; Chang et al., 2017). Our target genes included: *ACMSD, BST1, CCDC62, DDRGK1, FGF20, GAK, GCH1, GPNMB, HIP1R, INPP5F, ITGA8, LAMP3, LRRK2, MAPT, MCCC1, PM20D1, RAB25, RAB29, RIT2, SCARB2, SETD1A, SIPA1L2, SLC41A1, SNCA, SREBF1, STK39, STX1B, SYT11, TMEM163, TMEM175, USP25* and *VPS13C*. These genes included several known Parkinson’s disease related genes (*LRRK2, GCH1, SNCA*) found in GWAS loci in order to perform full sequencing assessment. Other genes were selected for sequencing based on assessment that included effects of quantitative trait loci, brain expression, potential association with known Parkinson’s disease-causing genes and involvement in pathways potentially involved in Parkinson’s disease, such the autophagy-lysosomal pathway, mitochondria quality control and endolysosomal recycling.

### Sequencing

The coding regions, exon-intron boundaries and the 5’ and 3’ untranslated regions (UTRs) of the genes listed above were targeted using molecular inversion probes (MIPs), as previously described (O’Roak et al., 2012). All MIPs used to sequence the genes of the study are detailed in Supplementary Table 1. Targeted DNA capture and amplification was done as previously described (Ross et al., 2016), and the full protocol is available upon request. The library was sequenced using Illumina HiSeq 2500/4000 platform at the McGill University and Genome Quebec Innovation Centre.

View this table:
[Table 1:](http://medrxiv.org/content/early/2020/05/29/2020.05.29.20116111/T1)

Table 1: Results of SKAT-O analysis in the merged cohorts at 30x and 50x minimum depth of coverage

### Data processing and quality control

Reads were mapped to the human reference genome (hg19) with Burrows-Wheeler Aligner (Li and Durbin, 2009). Post-alignment quality control and variant calling was performed using Genome Analysis Toolkit (GATK, v3.8) (McKenna et al., 2010). Variant annotation was done with ANNOVAR (Wang et al., 2010). Data on the frequency of each variant in different populations were extracted from the public database Genome Aggregation Database (GnomAD) (Lek et al., 2016).

The quality control pipeline that was used in this study is available on GitHub ([https://github.com/gan-orlab/MIPVar.git](https://github.com/gan-orlab/MIPVar.git)). All calls with less than 25% reads with the variant have been removed. The remaining quality control (QC) was performed using GATK and PLINK software v1.9 (Purcell et al., 2007). Genotyping quality (GQ) < 30 was used as a threshold for variant exclusion. Different depths of coverage thresholds were used for common and rare variants. Common variants threshold for inclusion was ≥ 15x, while rare variants (minor allele frequency [MAF] < 0.01) thresholds were more stringent and repeated twice, at ≥ 30x and ≥ 50x. Variants with genotype rate lower than 90% were excluded. Variants that deviated from Hardy-Weinberg equilibrium set at *p* ≥ 0.001 were filtered out. Threshold for missingness difference between cases and controls was set at *p =* 0.05 and the filtration script adjusted it with Bonferroni correction. Samples with average genotyping rate of less than 90% were excluded.

### Statistical analysis

The burden of rare variants in each of the target genes was tested using the optimized sequence Kernel association test (SKAT-O, R package) (Lee et al., 2012). In addition, we examined the burden of specific variant subgroups including: i) Funct.: Potentially functional included all nonsynonymous mutations, stop gain/loss, frameshift mutations, splicing variants located within two base pairs of exon-intron junctions, and variants that were flagged by ENCODE (explained below). ii) NS: nonsynonymous variants. iii) LOF: loss-of-function variants included stop gain/loss, frameshift, and splicing variants located within two base pairs of exon-intron junctions. iv) CADD: variants with the Combined Annotation Dependent Depletion (CADD) score ≥ 12.37; this threshold indicates the top 2% variants predicted to be the most damaging in the genome (Amendola et al., 2015). v) ENCODE: variants that are predicted to have regulatory activity, such as promoter, enhancer, activator and repressor, based on the chromatin profiling studies (Ernst et al., 2011). To increase power, rare variants were assessed in the three merged cohorts (since rare variants are less affected by ethnicity). Analysis of the individual cohorts was also performed post-hoc. Similarly to exome-wide burden analyses of rare variants (Kenna et al., 2016), Bonferroni correction for statistical significance was performed based on the number of genes that were tested, and the thresholds are *p* < 1.85E-03 for genes with variant coverage of ≥ 50x and *p* < 1.56E-03 at ≥ 30x

The association between common variants and Parkinson’s Disease was examined by logistic regression models using PLINK v1.9, with the status (patient or control) as a dependent variable, age and sex as covariates in all cohorts. Ethnicity as an additional covariate in the Columbia University cohort due to the ethnical heterogeneity. Meta-analysis of our three cohorts was conducted using METAL (R package) (Willer et al., 2010) with a fixed effect model. Cochran’s Q-test was applied to test for residual heterogeneity. LD of common variants and their respective GWAS top marker was examined using the LDlink online tool Dmatrix ([https://ldlink.nci.nih.gov/?tab=ldmatrix](https://ldlink.nci.nih.gov/?tab=ldmatrix)) (Machiela and Chanock, 2015). Effects of common and rare variants on expression levels (eQTL) or splicing (sQTL) was evaluated using GTEx Portal online tool ([https://www.gtexportal.org/home/](https://www.gtexportal.org/home/)).

### Data availability

Anonymized data is available upon request by any qualified investigator.

## Results

### Coverage and variants identified

The average coverage for all 32 genes included in the current analysis was 642X, with a range of 86-1,064 (median 688) across all genes. An average of 97% of the targeted regions was covered at > 15X, 96% at > 30X, and 92% at > 50X. Supplementary Table 2 details the average coverage and the percent coverage for each gene. At coverage of > 30X, we found 3,327 rare variants and at > 50X, 1,713 rare were identified (Supplementary Tables 3 and 4). A total of 311 common variants were identified across all genes in all three cohorts (Supplementary Table 5) with a coverage of > 15X.

View this table:
[Table 2.](http://medrxiv.org/content/early/2020/05/29/2020.05.29.20116111/T2)

Table 2. The common variants with statistically significant association with Parkinson’s Disease in the meta-analysis of three cohorts

### Rare variants in *SYT11, FGF20* and *GCH1* are associated with Parkinson’s disease

Table 1 details the SKAT-O results of the analysis in the three combined cohorts, and Supplementary Table 6 details the results per cohort. After Bonferroni correction, *SYT11, FGF20* and *GCH1* were associated with Parkinson’s Disease. Nominal associations were identified in 21 genes (Table 1, Supplementary Table 6). When considering all variants with MAF < 0.01 in *SYT11*, there was a strong association with Parkinson’s Disease (*p* = 5.3E-06, Table 1), mainly driven by a *SYT11* 3’ UTR variant (rs945006601) with MAF of 0.004 in Parkinson’s Disease patients and 0.00015 in controls (OR = 27.0, 95%CI 3.6–202.6, *p* = 0.0013). As was previously reported (Blauwendraat et al., 2018), LD analysis showed that this variant is in some LD with the known Parkinson’s disease-associated variants in *GBA* p.Glu326Lys (p.E326K, r2 = 1.5E-05, D’ = 1) and p.Asn370Ser (p.N370S, r2 = 4.23E-05, D’ = 1). The low r2 for both variants are due to the different allele frequencies, yet the D’ of 1 indicates an LD. To test if the identified association in *SYT11* is driven by these two pathogenic *GBA* variants, we removed the carriers of these variants and reran the *SYT11* SKAT-O analysis (*p =* 5.58E-05). Since we have full sequencing data on *GBA* for all participants, we further excluded carriers of any other *GBA* variant that was ever reported to be potentially pathogenic and the association remained (*p =* 5.23E-05).

When analyzing *FGF20*, a burden of all rare variants was demonstrated (*p* = 0.0002), mainly driven by a rare 5’ UTR variant, rs1034608171, which is located in the promoter region of the gene, with MAF of 0.002 in Parkinson’s Disease patients and 0.00015 in controls (OR = 13.4, 95% CI 1.7–106.2, *p* = 0.014). The association of *GCH1* with Parkinson’s Disease is driven by rare nonsynonymous variants (Supplementary Tables 3 and 4), some of which are known to cause dopamine responsive dystonia as we have previously reported(Rudakou et al., 2019).

After removing the *LRRK2* p.Gly2019Ser variant from the analysis, there was no statistically significant burden of other rare *LRRK2* variants in our cohorts. However, we did identify two variants, p.Arg793Met (rs35173587) and p.Gln1353Lys (rs200526782), in ten and eight controls, respectively, but not in Parkinson’s Disease. Neither of these two variants had a statistically significant association after Bonferroni correction.

### Association of common variants with potential functional effects in known Parkinson’s disease loci

To examine the association of common variants in the 32 genes, we performed logistic regression adjusted for age and sex in all three cohorts, and also adjusted for ethnicity in the Columbia University cohort. Additionally, our three cohorts were meta-analyzed. After Bonferroni correction we found 24 variants from five genes associated with risk of Parkinson’s Disease (Table 2). Forest plots in Fig. 1 show the meta-analysis results for the tagging variants of these six genes. We further examined the LD of these variants with the known Parkinson’s Disease-associated variants from the most recent Parkinson’s Disease GWAS, and their potential functional effects (detailed in Methods) in ENCODE, GTEx and UCSC genome browser. In the *MAPT* locus, 17 common variants at full LD, all represent the *MAPT* H2 haplotype, were associated with reduced risk of Parkinson’s Disease (Table 2). The variants identified in *TMEM175, BST1, SNCA* and *GPNMB* are all in strong LD with known GWAS hits in their respective loci (Table 2). Two *SNCA* 3’ UTR variants were in full LD between themselves, and in strong LD with rs356219, a known 3’ GWAS hit in the *SNCA* locus. These two variants are located within a strong enhancer region and are associated with altered splicing of *SNCA* in the Cortex (Fig. 2). Three *GPNMB* variants were identified, one of which (rs5850) is located at the 3’ UTR of *GPNMB* and associated with reduced expression of *GPNMB* in numerous brain regions including the substantia nigra, putamen and cortex (Fig. 3). A coding variant in *PM20D1*, p.Ile149Val (rs1891460), was nominally associated with reduced risk of Parkinson’s Disease (OR = 0.73, 95% CI 0.60–0.89, *p* = 1.161E-03, specific frequencies in each cohort are in Supplementary Table 5). This variant is located within the *PARK16* region which also includes *RAB29* (previously named *RAB7L1*), *SLC41A1* and other genes. However, this variant is not in LD with the top GWAS hits within this locus (r2 = 0.04, D’ = 0.46).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/29/2020.05.29.20116111/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/05/29/2020.05.29.20116111/F1)

Figure 1. Forest plots, meta-analyses of common variants in the three cohorts.
Each panel presents the individual odds ratio and 95% confidence interval for each cohort, and the results of the meta-analysis for the tagging variant of six loci identified in our study. Meta-analysis *p-*values: A. *MAPT* (*p =* 8.28E-10); B. *TMEM175* (*p =* 2.11E-04); C. *BST1* (*p =* 3.14E-05); D. *SNCA* (*p =* 4.08E-05); E. *GPNMB* (*p =* 2.19E-05)

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/29/2020.05.29.20116111/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/05/29/2020.05.29.20116111/F2)

Figure 2. Effect of 3’UTR *SNCA* variants on splicing.
Figure is processed from GTeX ([https://www.gtexportal.org](https://www.gtexportal.org)), and represent the splicing quantitative trait locus (sQTL) effects of: **A**. rs1045722 in the 3’ UTR of *SNCA* and **B**. rs3857053 in the 3’ UTR of *SNCA*

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/29/2020.05.29.20116111/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/05/29/2020.05.29.20116111/F3)

Figure 3. The 3’UTR *GPNMB* variant (rs5850) is associated with reduced expression in numerous brain regions.
Figure is processed from GTeX ([https://www.gtexportal.org](https://www.gtexportal.org)), and shows the expression quantitative trait locus (eQTL) effect of the 3’ UTR *GPNMB* variant rs5850 on GPNMB expression in different brain regions.

## Discussion

In the current study, we examined the role of common and rare variants in 32 genes located within known Parkinson’s Disease risk loci, by performing targeted next generation sequencing of the coding and 5’ and 3’ untranslated regions, followed by association tests. We found a statistically significant burden, after Bonferroni correction, of rare variants in *SYT11, FGF20* and *GCH1*, and nominal associations (*p* < 0.05) with other genes and variants (Table 1, Supplementary Tables 3,4,6). In addition, we have identified common coding and putative regulatory variants associated with Parkinson’s Disease in *MAPT, TMEM175, BST1, SNCA* and *GPNMB*, including a potential novel association with a nonsynonymous variant in *PM20D1*, p.Ile149Val, which is nominally significant and requires further confirmation.

Three genes, *SYT11, FGF20* and *GCH1* were identified with a statistically significant burden of rare variants comparing Parkinson’s Disease to controls. *SYT11* is located close to the *GBA* locus, and it was reported that the GWAS signal from *SYT11-GBA* locus in Parkinson’s Disease can be attributed to coding *GBA* variants (Blauwendraat et al., 2018). The *SYT11* variant mainly responsible for the association in our data is in weak LD with two coding *GBA* variants, p.Glu326Lys (p.E326K) and p.Asn370Ser (p.N370S). After exclusion of carriers of these variants, and further exclusion of all other *GBA* variants, the association of *SYT11* remained strong. These results suggest that *SYT11* may independently be associated with risk of PD. The Protein encoded by *SYT11*, Synaptotagmin 11, belongs to the synaptotagmin family, known calcium sensors mediating calcium-dependent regulation of membrane trafficking in synaptic transmission (Xie et al., 2017). Synaptotagmin 11 is involved in regulation of endocytosis and vesicle recycling process, which is essential for neurotransmission (Wang et al., 2016; Xie et al., 2017). Overexpression of Synaptotagmin 11 in mice led to impairment of DA transmission (Wang et al., 2018). Synaptotagmin 11 was also reported to be a physiological substrate for E3 ubiquitin ligase – parkin, playing a critical role in parkin-related neurotoxicity (Wang et al., 2018). Involvement of Synaptotagmin 11 in other mechanisms potentially related to Parkinson’s Disease such as autophagy-lysosomal pathway (Bento et al., 2016), phagocytosis of α-synuclein fibrils (Du et al., 2017), and repair of injured astrocytes via lysosome exocytosis regulation (Sreetama et al., 2016), suggests that *SYT11* needs to be considered as a target for further research in PD. Fibroblast growth factor 20 (*FGF20*) encodes a neurotrophic factor preferentially expressed in the substantia nigra pars compacta. It regulates central nervous system development and function (International Parkinson Disease Genomics Consortium, 2011) and plays a role in dopaminergic neurons differentiation and survival (Itoh and Ohta, 2013). There is a number of conflicting reports about FGF20 overexpression increasing *a*-synuclein levels in dopaminergic neurons (Wang et al., 2008; Wider et al., 2009; Sekiyama et al., 2014; Tarazi et al., 2014). The rare variant (rs1034608171) that drives the association in our data is located in a regulatory region and its potential functional effects would require additional studies. *GCH1* encodes GTP-cyclohydrolase 1, an essential enzyme for the dopamine synthesis pathway in the nigrostriatal cells (Kurian et al., 2011). Rare *GCH1* variants may cause DOPA-responsive dystonia, and the same variants have also been associated with Parkinson’s Disease (Mencacci et al., 2014; Rudakou et al., 2019). Since these mutations are rare, they only reached statistical significance in the McGill cohort, as we have previously described (Rudakou et al., 2019). The two *LRRK2* variants, p.Arg793Met and p.Gln1353Lys, that were found only in controls require additional studies to examine their potential protective roles in Parkinson’s Disease. *LRRK2* is a well-validated Parkinson’s Disease related gene(Bandres-Ciga et al., 2020).

All 17 *MAPT* variants that were found to be significantly associated with Parkinson’s Disease are in perfect LD with each other, and represent the H2 haplotype, known to be associated with reduced risk of Parkinson’s Disease (Li et al., 2018). The identified coding variants in *MAPT* are located within exons that are likely not expressed in the CNS (Caillet-Boudin et al., 2015), suggesting that the association in this locus may be due to regulatory variants. *TMEM175* is located in the fourth most significant GWAS locus in Parkinson’s Disease (Nalls et al., 2019), and encodes a transmembrane potassium channel responsible for K+ conductance in endosomes and lysosomes (Cang et al., 2015). TMEM175 regulates lysosomal function through maintenance of lysosomal membrane potential and luminal pH stability, as well as its role in autophagosome-lysosome fusion (Lee et al., 2017). In a recent genetic and functional study, two coding *TMEM175* variants, including the variant identified here, were shown to be associated with risk of Parkinson’s Disease, risk of REM-sleep behavior disorder, and glucocerebrosidase activity (Krohn et al., 2019a).

The *BST1* and *SNCA* variants associated with risk of Parkinson’s Disease in the current study are in strong LD with the known GWAS markers in their respective loci. The *BST1* variant is intronic without known function, and the *SNCA* 3’UTR variants are located within strong *SNCA* enhancers. Whether these *SNCA* variants drive the association with Parkinson’s Disease and their potential functional effects require additional studies. *GPNMB* encodes the glycoprotein nonmetastatic melanoma protein B a type 1, which is a transmembrane glycoprotein shown to have a role in melanoma (Taya and Hammes, 2018). Patients with melanoma have an increased risk of Parkinson’s Disease (Liu et al., 2011), and Parkinson’s Disease patients have an increased risk for melanoma (Bertoni et al., 2010). GPNMB may have a role in the lysosome (Tomihari et al., 2009) and in regulation of microglial inflammation (Dzamko et al., 2015; Herrero et al., 2015). While the expression of GPNMB in brain is overall low, the two variants identified in the current study are associated with reduced GPNMB expression in various brain tissues (Fig. 3). One of the variants, rs5850, is located within the 3’ UTR of *GPNMB*, but there is no clear regulatory function for this variant.

Peptidase M20 Domain Containing 1 (*PM20D1*) encodes a protein that regulates the production of N-fatty-acyl amino acids which have a function of chemical uncouplers of mitochondrial respiration (Long et al., 2016). Mitochondrial quality control is likely important in Parkinson’s Disease, and genes involved in mitophagy are involved in familial forms of Parkinson’s Disease (Park et al., 2018). *PM20D1* is located within the *PARK16* locus which also contains *SLC45A3, NUCKS1, RAB29* (previously named *RAB7L1*) and *SLC41A1* (Satake et al., 2009; Simon-Sanchez et al., 2009; Nalls et al., 2019).Several studies have highlighted *RAB29* as the potential gene associated with Parkinson’s Disease in this locus (Gan-Or et al., 2012),as it was shown to interact with *LRRK2* (MacLeod et al., 2013) and affect lysosomal function (Eguchi et al., 2018). In the current study we identified a common nonsynonymous variant in *PM20D1*, p.Ile149Val (rs1891460), which is not in LD with the Parkinson’s Disease GWAS marker in this locus. Since this variant is not reported in ClinVar, has a low CADD score (1.01), and the statistical association is of borderline significance when correcting by the number of genes, it is possible that this association was found by chance. Replications in other cohorts will be required to determine whether this variant is associated with Parkinson’s Disease.

Our study has several limitations. The Parkinson’s disease patients and controls are not matched by sex and age, hence we adjusted for these variables in the statistical analysis. We found no significant difference between allele frequencies of young (<40 years old) and old controls, which allowed us to use all of the controls in our analyses. In addition, the coverage of some of the genes we analyzed was suboptimal (Supplementary Table 2), and we did not target intronic regions and other potentially regulatory regions outside the 5’ and 3’ UTRs, which may also affect gene function. Whole genome sequencing studies will be required to study these regions and overcome the coverage limitations.

Our results, combined with previous studies, highlight the importance of performing fine mapping of GWAS loci to identify the variants and genes within each locus that drive the reported associations. Thus far, the causative genes within most Parkinson’s Disease GWAS loci are unknown, with the exception of *SNCA, GBA, LRRK2*, and likely *VPS13C, MAPT, TMEM175* and *GCH1*. The current study suggests that *SYT11, FGF20*, and potentially *PM20D1, BST1* and *GPNMB* should be considered as possible Parkinson’s Disease-related genes and targets for further research.

## Data Availability

Anonymized data is available upon request by any qualified investigator.

## Funding

This work was financially supported by grants from the Michael J. Fox Foundation, the Canadian Consortium on Neurodegeneration in Aging (CCNA), the Canada First Research Excellence Fund (CFREF), awarded to McGill University for the Healthy Brains for Healthy Lives initiative (HBHL), and Parkinson Canada. The Columbia University cohort is supported by the Parkinson’s Foundation, the National Institutes of Health (K02NS080915, and UL1 TR000040) and the Brookdale Foundation.

## Competing interests

Ziv Gan-Or has received consulting fees from Lysosomal Therapeutics Inc., Idorsia, Prevail Therapeutics, Denali, Ono Therapeutics, Deerfield and Inception Sciences (now Ventus). None of these companies were involved in any parts of preparing, drafting and publishing this review. Cheryl H. Waters disclosed the following: Research support (Biogen, Roche, Sanofi); consulting fees (Acorda, Amneal, Impel, Kyowa, Mitsubishi, and Neurocrine); speaker’s honoraria (Acadia, Acorda, Adamas, Amneal, and US WorldMeds)

## Supplementary material

Supplementary Table 1 Detailed information on molecular inversion probes (MIPs)

Supplementary Table 2 Coverage statistics

Supplementary Table 3 All identified rare variants with 30x minimum depth of coverage

Supplementary Table 4 All identified rare variants with 50x minimum depth of coverage

Supplementary Table 5 All identified common variants

Supplementary Table 6 SKAT-O results for individual cohorts

## Acknowledgements

We thank the participants for contributing to the study. GAR holds a Canada Research Chair in Genetics of the Nervous System and the Wilder Penfield Chair in Neurosciences. EAF is supported by a Foundation Grant from the Canadian Institutes of Health Research (FDN grant – 154301). ZGO is supported by the Fonds de recherche du Québec – Santé (FRQS) Chercheursboursiers award, in collaboration with Parkinson Quebec, and by the Young Investigator Award by Parkinson Canada. The access to part of the participants for this research has been made possible thanks to the Quebec Parkinson’s Network ([http://rpq-qpn.ca/en/](http://rpq-qpn.ca/en/)). We thank Daniel Rochefort, Helene Catoire, Sandra B. Laurent, Clotilde Degroot and Vessela Zaharieva for their assistance.

## Abbreviations

AJ
:   Ashkenazi Jewish
CADD
:   Combined Annotation Dependent Depletion
CI
:   confidence interval
CNS
:   central nervous system
DOPA
:   dihydroxyphenylalanine
DP
:   depth of coverage
eQT
:   expression quantitative trait loci
GATK
:   Genome Analysis Toolkit
GnomAD
:   Genome Aggregation Database
GTEx
:   Genotype-tissue expression
GWAS
:   Genome-wide association study
hg19
:   human genome version 19
LD
:   linkage disequilibrium
LOF
:   loss-of-function
MAF
:   minor allele frequency
MIPs
:   molecular inversion probes
NS, Nonsyn
:   non-synonymous
OR
:   odds ratio
QC
:   quality control
REM
:   Rapid Eye Movement
SKAT-O
:   optimized sequence kernel association test
sQTL
:   splicing quantitative trait loci
UTR
:   untranslated region
yo
:   years of age

*   Received May 29, 2020.
*   Revision received May 29, 2020.
*   Accepted May 29, 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.  Alcalay RN, Levy OA, Wolf P, Oliva P, Zhang XK, Waters CH, et al. SCARB2 variants and glucocerebrosidase activity in Parkinson’s disease. NPJ Parkinsons Dis 2016; 2.
    
    

2.  Amaral CEM, Lopes PF, Ferreira JCC, Alves EAC, Montenegro MVB, Costa ETD, et al. GBA mutations p.N370S and p.L444P are associated with Parkinson’s disease in patients from Northern Brazil. Arq Neuropsiquiatr 2019; 77(2): 73–9.
    
    

3.  Amendola LM, Dorschner MO, Robertson PD, Salama JS, Hart R, Shirts BH, et al. Actionable exomic incidental findings in 6503 participants: challenges of variant classification. Genome Res 2015; 25(3): 305–15.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjI1LzMvMzA1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMjkvMjAyMC4wNS4yOS4yMDExNjExMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

4.  Bandres-Ciga S, Diez-Fairen M, Kim JJ, Singleton AB. Genetics of Parkinson’s disease: An introspection of its journey towards precision medicine. Neurobiol Dis 2020; 137: 104782.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.nbd.2020.104782&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

5.  Bento CF, Ashkenazi A, Jimenez-Sanchez M, Rubinsztein DC. The Parkinson’s disease-associated genes ATP13A2 and SYT11 regulate autophagy via a common pathway. Nat Commun 2016; 7: 11803.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms11803&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27278822&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

6.  Bertoni JM, Arlette JP, Fernandez HH, Fitzer-Attas C, Frei K, Hassan MN, et al. Increased melanoma risk in Parkinson disease: a prospective clinicopathological study. Arch Neurol 2010; 67(3): 347–52.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archneurol.2010.1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20212233&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275308600013&link_type=ISI) 

7.  Blauwendraat C, Bras JM, Nalls MA, Lewis PA, Hernandez DG, Singleton AB, et al. Coding variation in GBA explains the majority of the SYT11-GBA Parkinson’s disease GWAS locus. Mov Disord 2018; 33(11): 1821–3.
    
    

8.  Caillet-Boudin ML, Buee L, Sergeant N, Lefebvre B. Regulation of human MAPT gene expression. Mol Neurodegener 2015; 10: 28.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13024-015-0025-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26170022&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

9.  Cang C, Aranda K, Seo YJ, Gasnier B, Ren D. TMEM175 Is an Organelle K(+) Channel Regulating Lysosomal Function. Cell 2015; 162(5): 1101–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2015.08.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26317472&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

10. Chang D, Nalls MA, Hallgrimsdottir IB, Hunkapiller J, van der Brug M, Cai F, et al. A meta-analysis of genome-wide association studies identifies 17 new Parkinson’s disease risk loci. Nat Genet 2017; 49(10): 1511–6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3955&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28892059&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

11. Darvish H, Bravo P, Tafakhori A, Azcona LJ, Ranji-Burachaloo S, Johari AH, et al. Identification of a large homozygous VPS13C deletion in a patient with early-onset Parkinsonism. Mov Disord 2018; 33(12): 1968–70.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mds.27516&link_type=DOI) 

12. Du C, Wang Y, Zhang F, Yan S, Guan Y, Gong X, et al. Synaptotagmin-11 inhibits cytokine secretion and phagocytosis in microglia. Glia 2017; 65(10): 1656–67.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/glia.23186&link_type=DOI) 

13. Dzamko N, Geczy CL, Halliday GM. Inflammation is genetically implicated in Parkinson’s disease. Neuroscience 2015; 302: 89–102.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroscience.2014.10.028&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25450953&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

14. Eguchi T, Kuwahara T, Sakurai M, Komori T, Fujimoto T, Ito G, et al. LRRK2 and its substrate Rab GTPases are sequentially targeted onto stressed lysosomes and maintain their homeostasis. Proc Natl Acad Sci U S A 2018; 115(39): E9115-E24.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE1LzM5L0U5MTE1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMjkvMjAyMC4wNS4yOS4yMDExNjExMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

15. Emelyanov AK, Usenko TS, Tesson C, Senkevich KA, Nikolaev MA, Miliukhina IV, et al. Mutation analysis of Parkinson’s disease genes in a Russian data set. Neurobiol Aging 2018; 71: 267 e7-e10.
    
    

16. Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, et al. Mapping and analysis of chromatin state dynamics in nine human cell types. Nature 2011; 473(7345): 43–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature09906&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21441907&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000290218300027&link_type=ISI) 

17. Gan-Or Z, Bar-Shira A, Dahary D, Mirelman A, Kedmi M, Gurevich T, et al. Association of sequence alterations in the putative promoter of RAB7L1 with a reduced parkinson disease risk. Arch Neurol 2012; 69(1): 105–10.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archneurol.2011.924&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22232350&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000298944500015&link_type=ISI) 

18. Gan-Or Z, Rao T, Leveille E, Degroot C, Chouinard S, Cicchetti F, et al. The Quebec Parkinson Network: A Researcher-Patient Matching Platform and Multimodal Biorepository. J Parkinsons Dis 2020; 10(1): 301–13.
    
    

19. Guella I, Sherman HE, Appel-Cresswell S, Rajput A, Rajput AH, Farrer MJ. Parkinsonism in GTP cyclohydrolase 1 mutation carriers. Brain 2015; 138(Pt 5): e349.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awu341&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25497597&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

20. Herrero MT, Estrada C, Maatouk L, Vyas S. Inflammation in Parkinson’s disease: role of glucocorticoids. Front Neuroanat 2015; 9: 32.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnana.2015.00032&link_type=DOI) 

21. Hughes AJ, Daniel SE, Kilford L, Lees AJ. Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: a clinico-pathological study of 100 cases. J Neurol Neurosurg Psychiatry 1992; 55(3): 181–4.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiam5ucCI7czo1OiJyZXNpZCI7czo4OiI1NS8zLzE4MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzI5LzIwMjAuMDUuMjkuMjAxMTYxMTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

22. International Parkinson Disease Genomics Consortium (IPDGC). Imputation of sequence variants for identification of genetic risks for Parkinson’s disease: a meta-analysis of genome-wide association studies. Lancet 2011; 377(9766): 641–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(10)62345-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21292315&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000287776600033&link_type=ISI) 

23. Itoh N, Ohta H. Roles of FGF20 in dopaminergic neurons and Parkinson’s disease. Front Mol Neurosci 2013; 6: 15.
    
    

24. Jansen IE, Ye H, Heetveld S, Lechler MC, Michels H, Seinstra RI, et al. Discovery and functional prioritization of Parkinson’s disease candidate genes from large-scale whole exome sequencing. Genome Biol 2017; 18(1): 22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-017-1147-9&link_type=DOI) 

25. Jinn S, Blauwendraat C, Toolan D, Gretzula CA, Drolet RE, Smith S, et al. Functionalization of the TMEM175 p.M393T variant as a risk factor for Parkinson disease. Hum Mol Genet 2019; 28(19): 3244–54.
    
    

26. Kenna KP, van Doormaal PT, Dekker AM, Ticozzi N, Kenna BJ, Diekstra FP, et al. NEK1 variants confer susceptibility to amyotrophic lateral sclerosis. Nat Genet 2016; 48(9): 1037–42.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3626&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27455347&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

27. Khan NL, Jain S, Lynch JM, Pavese N, Abou-Sleiman P, Holton JL, et al. Mutations in the gene LRRK2 encoding dardarin (PARK8) cause familial Parkinson’s disease: clinical, pathological, olfactory and functional imaging and genetic data. Brain 2005; 128(Pt 12): 2786–96.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awh667&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16272164&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233667500006&link_type=ISI) 

28. Krohn L, Ozturk TN, Vanderperre B, Ouled Amar Bencheikh B, Ruskey JA, Laurent SB, et al. Genetic, Structural, and Functional Evidence Link TMEM175 to Synucleinopathies. Ann Neurol 2019a.
    
    

29. Krohn L, Ozturk TN, Vanderperre B, Ouled Amar Bencheikh B, Ruskey JA, Laurent SB, et al. Genetic, Structural, and Functional Evidence Link TMEM175 to Synucleinopathies. Ann Neurol 2020; 87(1): 139–53.
    
    

30. Krohn L, Wu RYJ, Heilbron K, Ruskey JA, Laurent SB, Blauwendraat C, et al. Fine-mapping of SNCA in REM sleep behavior disorder and overt synucleinopathies. BioRxiv 2019b: 756528.
    
    

31. Kurian MA, Gissen P, Smith M, Heales S, Jr.., Clayton PT. The monoamine neurotransmitter disorders: an expanding range of neurological syndromes. Lancet Neurol 2011; 10(8): 721–33.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1474-4422(11)70141-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21777827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000293612600014&link_type=ISI) 

32. Lee C, Guo J, Zeng W, Kim S, She J, Cang C, et al. The lysosomal potassium channel TMEM175 adopts a novel tetrameric architecture. Nature 2017; 547(7664): 472–5.
    
    

33. Lee S, Emond MJ, Bamshad MJ, Barnes KC, Rieder MJ, Nickerson DA, et al. Optimal unified approach for rare-variant association testing with application to small-sample case-control whole-exome sequencing studies. Am J Hum Genet 2012; 91(2): 224–37.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2012.06.007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22863193&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

34. Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 2016; 536(7616): 285–91.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature19057&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27535533&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000381804900026&link_type=ISI) 

35. Lesage S, Drouet V, Majounie E, Deramecourt V, Jacoupy M, Nicolas A, et al. Loss of VPS13C Function in Autosomal-Recessive Parkinsonism Causes Mitochondrial Dysfunction and Increases PINK1/Parkin-Dependent Mitophagy. Am J Hum Genet 2016; 98(3): 500–13.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2016.01.014&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26942284&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

36. Lewthwaite AJ, Lambert TD, Rolfe EB, Olgiati S, Quadri M, Simons EJ, et al. Novel GCH1 variant in Dopa-responsive dystonia and Parkinson’s disease. Parkinsonism Relat Disord 2015; 21(4): 394–7.
    
    

37. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009; 25(14): 1754–60.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btp324&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19451168&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267665900006&link_type=ISI) 

38. Li J, Ruskey JA, Arnulf I, Dauvilliers Y, Hu MTM, Hogl B, et al. Full sequencing and haplotype analysis of MAPT in Parkinson’s disease and rapid eye movement sleep behavior disorder. Mov Disord 2018; 33(6): 1016–20.
    
    

39. Liu R, Gao X, Lu Y, Chen H. Meta-analysis of the relationship between Parkinson disease and melanoma. Neurology 2011; 76(23): 2002–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e31821e554e&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21646627&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

40. Long JZ, Svensson KJ, Bateman LA, Lin H, Kamenecka T, Lokurkar IA, et al. The Secreted Enzyme PM20D1 Regulates Lipidated Amino Acid Uncouplers of Mitochondria. Cell 2016; 166(2): 424–35.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2016.05.071&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27374330&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

41. Machiela MJ, Chanock SJ. LDlink: a web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics 2015; 31(21): 3555–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv402&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26139635&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

42. MacLeod DA, Rhinn H, Kuwahara T, Zolin A, Di Paolo G, McCabe BD, et al. RAB7L1 interacts with LRRK2 to modify intraneuronal protein sorting and Parkinson’s disease risk. Neuron 2013; 77(3): 425–39.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2012.11.033&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23395371&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000317030800007&link_type=ISI) 

43. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 2010; 20(9): 1297–303.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjk6IjIwLzkvMTI5NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzI5LzIwMjAuMDUuMjkuMjAxMTYxMTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

44. Mencacci NE, Isaias IU, Reich MM, Ganos C, Plagnol V, Polke JM, et al. Parkinson’s disease in GTP cyclohydrolase 1 mutation carriers. Brain 2014; 137(Pt 9): 2480–92.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awu179&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24993959&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000342996000017&link_type=ISI) 

45. Nalls MA, Blauwendraat C, Vallerga CL, Heilbron K, Bandres-Ciga S, Chang D, et al. Identification of novel risk loci, causal insights, and heritable risk for Parkinson’s disease: a meta-analysis of genome-wide association studies. Lancet Neurol 2019; 18(12): 1091–102.
    
    [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%2F05%2F29%2F2020.05.29.20116111.atom) 

46. Nalls MA, Pankratz N, Lill CM, Do CB, Hernandez DG, Saad M, et al. Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson’s disease. Nat Genet 2014; 46(9): 989–93.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3043&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25064009&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

47. O’Roak BJ, Vives L, Fu W, Egertson JD, Stanaway IB, Phelps IG, et al. Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders. Science 2012; 338(6114): 1619–22.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMzgvNjExNC8xNjE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMjkvMjAyMC4wNS4yOS4yMDExNjExMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

48. Paisan-Ruiz C, Jain S, Evans EW, Gilks WP, Simon J, van der Brug M, et al. Cloning of the gene containing mutations that cause PARK8-linked Parkinson’s disease. Neuron 2004; 44(4): 595–600.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2004.10.023&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15541308&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225338100006&link_type=ISI) 

49. Park JS, Davis RL, Sue CM. Mitochondrial Dysfunction in Parkinson’s Disease: New Mechanistic Insights and Therapeutic Perspectives. Curr Neurol Neurosci Rep 2018; 18(5): 21.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11910-018-0829-3&link_type=DOI) 

50. Polymeropoulos MH, Lavedan C, Leroy E, Ide SE, Dehejia A, Dutra A, et al. Mutation in the alphasynuclein gene identified in families with Parkinson’s disease. Science 1997; 276(5321): 2045–7.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIyNzYvNTMyMS8yMDQ1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMjkvMjAyMC4wNS4yOS4yMDExNjExMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

51. Postuma RB, Berg D, Stern M, Poewe W, Olanow CW, Oertel W, et al. MDS clinical diagnostic criteria for Parkinson’s disease. Mov Disord 2015; 30(12): 1591–601.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mds.26424&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26474316&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

52. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007; 81(3): 559–75.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/519795&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17701901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

53. Ran C, Brodin L, Forsgren L, Westerlund M, Ramezani M, Gellhaar S, et al. Strong association between glucocerebrosidase mutations and Parkinson’s disease in Sweden. Neurobiol Aging 2016; 45: 212 e5-e11.
    
    

54. Ross JP, Dupre N, Dauvilliers Y, Strong S, Ambalavanan A, Spiegelman D, et al. Analysis of DNAJC13 mutations in French-Canadian/French cohort of Parkinson’s disease. Neurobiol Aging 2016; 45: 212 e13-e17.
    
    

55. Rudakou U, Ouled Amar Bencheikh B, Ruskey JA, Krohn L, Laurent SB, Spiegelman D, et al. Common and rare GCH1 variants are associated with Parkinson’s disease. Neurobiol Aging 2019; 73: 231 e1-e6.
    
    

56. Rudakou U, Ruskey JA, Krohn L, Laurent SB, Spiegelman D, Greenbaum L, et al. Analysis of common and rare VPS13C variants in late-onset Parkinson disease. Neurol Genet 2020; 6(1): 385.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoibm5nIjtzOjU6InJlc2lkIjtzOjc6IjYvMS8zODUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNS8yOS8yMDIwLjA1LjI5LjIwMTE2MTExLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

57. Ruskey JA, Greenbaum L, Ronciere L, Alam A, Spiegelman D, Liong C, et al. Increased yield of full GBA sequencing in Ashkenazi Jews with Parkinson’s disease. Eur J Med Genet 2019; 62(1): 65–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejmg.2018.05.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

58. Satake W, Nakabayashi Y, Mizuta I, Hirota Y, Ito C, Kubo M, et al. Genome-wide association study identifies common variants at four loci as genetic risk factors for Parkinson’s disease. Nat Genet 2009; 41(12): 1303–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.485&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19915576&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000272144900011&link_type=ISI) 

59. Schormair B, Kemlink D, Mollenhauer B, Fiala O, Machetanz G, Roth J, et al. Diagnostic exome sequencing in early-onset Parkinson’s disease confirms VPS13C as a rare cause of autosomalrecessive Parkinson’s disease. Clin Genet 2018; 93(3): 603–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/cge.13124&link_type=DOI) 

60. Sekiyama K, Takamatsu Y, Waragai M, Hashimoto M. Role of genomics in translational research for Parkinson’s disease. Biochem Biophys Res Commun 2014; 452(2): 226–35.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbrc.2014.06.028&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24950403&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

61. Sidransky E, Nalls MA, Aasly JO, Aharon-Peretz J, Annesi G, Barbosa ER, et al. Multicenter analysis of glucocerebrosidase mutations in Parkinson’s disease. N Engl J Med 2009; 361(17): 1651–61.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa0901281&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19846850&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000270977400007&link_type=ISI) 

62. Simon-Sanchez J, Schulte C, Bras JM, Sharma M, Gibbs JR, Berg D, et al. Genome-wide association study reveals genetic risk underlying Parkinson’s disease. Nat Genet 2009; 41(12): 1308–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.487&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19915575&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000272144900012&link_type=ISI) 

63. Soldner F, Stelzer Y, Shivalila CS, Abraham BJ, Latourelle JC, Barrasa MI, et al. Parkinson-associated risk variant in distal enhancer of alpha-synuclein modulates target gene expression. Nature 2016; 533(7601): 95–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature17939&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27096366&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 

64. Sreetama SC, Takano T, Nedergaard M, Simon SM, Jaiswal JK. Injured astrocytes are repaired by Synaptotagmin XI-regulated lysosome exocytosis. Cell Death Differ 2016; 23(4): 596–607.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/cdd.2015.124&link_type=DOI) 

65. Tarazi FI, Sahli ZT, Wolny M, Mousa SA. Emerging therapies for Parkinson’s disease: from bench to bedside. Pharmacol Ther 2014; 144(2): 123–33.
    
    

66. Taya M, Hammes SR. Glycoprotein Non-Metastatic Melanoma Protein B (GPNMB) and Cancer: A Novel Potential Therapeutic Target. Steroids 2018; 133: 102–7.
    
    

67. Tomihari M, Hwang SH, Chung JS, Cruz PD, Jr., Ariizumi K. Gpnmb is a melanosome-associated glycoprotein that contributes to melanocyte/keratinocyte adhesion in a RGD-dependent fashion. Exp Dermatol 2009; 18(7): 586–95.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1600-0625.2008.00830.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19320736&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266699400003&link_type=ISI) 

68. Wang C, Kang X, Zhou L, Chai Z, Wu Q, Huang R, et al. Synaptotagmin-11 is a critical mediator of parkin-linked neurotoxicity and Parkinson’s disease-like pathology. Nat Commun 2018; 9(1): 81.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-02593-y&link_type=DOI) 

69. Wang C, Wang Y, Hu M, Chai Z, Wu Q, Huang R, et al. Synaptotagmin-11 inhibits clathrin-mediated and bulk endocytosis. EMBO Rep 2016; 17(1): 47–63.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiZW1ib3IiO3M6NToicmVzaWQiO3M6NzoiMTcvMS80NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzI5LzIwMjAuMDUuMjkuMjAxMTYxMTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

70. Wang G, van der Walt JM, Mayhew G, Li YJ, Zuchner S, Scott WK, et al. Variation in the miRNA-433 binding site of FGF20 confers risk for Parkinson disease by overexpression of alpha-synuclein. Am J Hum Genet 2008; 82(2): 283–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2007.09.021&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18252210&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253223900002&link_type=ISI) 

71. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 2010; 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%2F05%2F29%2F2020.05.29.20116111.atom) 

72. Wider C, Dachsel JC, Soto AI, Heckman MG, Diehl NN, Yue M, et al. FGF20 and Parkinson’s disease: no evidence of association or pathogenicity via alpha-synuclein expression. Mov Disord 2009; 24(3): 455–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mds.22442&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19133659&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000263747200024&link_type=ISI) 

73. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 2010; 26(17): 2190–1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq340&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20616382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F29%2F2020.05.29.20116111.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281738900017&link_type=ISI) 

74. Xie Z, Long J, Liu J, Chai Z, Kang X, Wang C. Molecular Mechanisms for the Coupling of Endocytosis to Exocytosis in Neurons. Front Mol Neurosci 2017; 10: 47.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnmol.2017.00047&link_type=DOI) 

75. Xu Q, Li K, Sun Q, Ding D, Zhao Y, Yang N, et al. Rare GCH1 heterozygous variants contributing to Parkinson’s disease. Brain 2017; 140(7): e41.
    
    

76. Yoshino H, Nishioka K, Li Y, Oji Y, Oyama G, Hatano T, et al. GCH1 mutations in dopa-responsive dystonia and Parkinson’s disease. J Neurol 2018; 265(8): 1860–70.