Parkinson’s disease causality and heterogeneity: a proteogenomic view ======================================================================= * Sergio Kaiser * Luqing Zhang * Brit Mollenhauer * Jaison Jacob * Simonne Longerich * Jorge Del-Aguila * Jacob Marcus * Neha Raghavan * David Stone * Olumide Fagboyegun * Douglas Galasko * Mohammed Dakna * Bilada Bilican * Mary Dovlatyan * Anna Kostikova * Jingyao Li * Brant Peterson * Michael Rotte * Vinicius Sanz * Tatiana Foroud * Samantha J. Hutten * Mark Frasier * Hirotaka Iwaki * Andrew Singleton * Ken Marek * Karen Crawford * Fiona Elwood * Mirko Messa * Pablo Serrano-Fernandez ## Abstract The pathogenesis and clinical heterogeneity of Parkinson’s disease have been evaluated from molecular, pathophysiological, and clinical perspectives. High-throughput proteomic analysis of CSF has opened new opportunities for scrutinizing this heterogeneity. To date, this is the most comprehensive CSF-based proteomics profiling study in Parkinson’s disease (1103 patients, 4135 proteins). Combining CSF aptamer-based proteomics with genetics we determined protein quantitative trait loci (pQTLs). Analyses of pQTLs together with summary statistics from the largest Parkinson’s disease genome wide association study (GWAS) identified 68 potential causal proteins by Mendelian randomization. The top causal protein, GPNMB was previously reported to be upregulated in the substantia nigra of Parkinson’s disease patients. We also compared the CSF proteomes of patients and controls. The Parkinson’s disease cohort comprised not only *LRRK2*+ and *GBA*+ mutation carriers but also idiopathic patients. Proteome differences between *GBA+* patients and unaffected *GBA*+ controls suggest degeneration of dopaminergic neurons, altered dopamine metabolism and increased brain inflammation. The proteins discriminating *LRRK2+* patients from unaffected *LRRK2+* controls, revealed dysregulated lysosomal degradation, as well as altered alpha-synuclein processing, and neurotransmission. Proteome differences between idiopathic patients and controls suggest increased neuroinflammation, mitochondrial dysfunction / oxidative stress, altered iron metabolism and potential neuroprotection mediated by vasoactive substances. Finally, we used proteomic data to stratify idiopathic patients into “endotypes”. The identified endotypes show differences in cognitive and motor disease progression based on the use of previously reported protein-based risk scores. In summary, we: i) identified causal proteins for Parkinson’s disease, ii) assessed CSF proteome differences in Parkinson’s disease patients of genetic and idiopathic etiology, and. iii) stratified idiopathic patients into robust clinically relevant subtypes. Our findings not only contribute to the identification of new therapeutic targets but also to shaping personalized medicine in CNS neurodegeneration. Keywords * proteomics * genetics * proteogenomics * CSF * pQTL ## Introduction Parkinson’s disease is the second most prevalent neurodegenerative disorder.1 Parkinson’s disease patients experience selective degeneration and loss of dopaminergic neurons in the substantia nigra *pars compacta*. In most cases, Parkinson’s disease is classified as idiopathic, but a growing set of genetic variants increase Parkinson’s disease risk or accelerate its onset. Many of the identified genes are involved either in mitochondrial or endo-lysosomal biology.2 The two most common Parkinson’s disease risk genes are leucine rich kinase 2 (*LRRK2*) and glucosidase beta acid (*GBA*). Mutations in these genes are linked to ∼10% of sporadic cases and up to 30% in specific ethnic subgroups and familial disease.3 Some of the genetic variants in *LRRK2* and *GBA* have been associated with specific clinical phenotypes.4,5 For both genes, the pathological mutations are thought to exacerbate the toxicity of alpha-synuclein, which – in an aggregated form – contributes to neuronal death and amplifies the neuroinflammatory response. Clinical heterogeneity of Parkinson’s disease has motivated many disease stratification efforts. Some of those have focused on clinical variables, mostly hypothesis-driven, while others have focused on molecular data, mostly hypothesis-free.6 An association between such strata and the Parkinson’s disease risk mutations remains elusive.7 Proteins hold great potential as predictors, causal biomarkers and surrogates of disease progression and/or stratification. However, the biological and pathophysiological complexity of Parkinson’s disease, the difficulties of collecting standardized biological samples (especially CSF) from large cohorts throughout the course of disease, and the technical limitations of high-throughput proteomic analyses hamper the identification of biomarkers at a proteomic level. The multicenter Parkinson Progression Marker Initiative (PPMI) was initiated to overcome some of these limitations, particularly in terms of number of samples and clinical data 8. In this collaboration, 1103 baseline (not longitudinal) CSF samples from patients and control participants with known status of *LRRK2* and *GBA* pathogenic variants were analyzed using the SomaScan® aptamer-based proteomics platform.9 There is also whole genome sequencing data from 804 patients out of the 1103, after quality control. To our knowledge, this is to date the largest proteomic and genetic data set for interrogating causal proteins for Parkinson’s disease in a neurologically relevant biofluid. The main goals of this study are summarized in Figure 1: i. Parkinson’s disease causal protein identification using mendelian randomization based on proteomics and genetics, ii. identification of differences between Parkinson’s disease patients and controls within and between subcohorts (*LRRK2+, GBA+* and idiopathic), and iii. hypothesis-free stratification of idiopathic Parkinson’s disease patients into clinically relevant endotypes and then compared by means of protein-based risk scores reflecting cognitive and motor progression 10. ![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/14/2022.03.09.22272131/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/06/14/2022.03.09.22272131/F1) Figure 1. Full analytical design. 1103 subjects were analyzed with SomaScan for 4135 unique proteins in CSF. The comparison between *GBA+* patients and *GBA+* unaffected controls (UC) retrieved six differentially expressed (DE) proteins. The comparison between *LRRK2+* patients and *LRRK2+* UC retrieved seven DE proteins. The comparison between idiopathic patients and healthy controls (HC) non-mutation carriers retrieved 23 DE proteins. Patients and controls were also combined and compared, which retrieved 122 DE proteins. Idiopathic patients were further analyzed, and two endotypes were identified based on CSF proteomics. 1264 subjects were sequenced genome wide to detect a total of 9743041 SNPs. For the 804 patients that had both genomic and proteomic data, a pQTL analysis was performed that identified 744 unique proteins with a significant cis-pQTL. The pQTLs combined with a meta GWAS for Parkinson’s disease performed by Nalls *et al*. (2019), led to the proposal of 68 unique CSF proteins presumed to be causal for Parkinson’s disease. ## Materials and methods The clinical data and samples used in this study were obtained from the PPMI ([http://www.ppmi-info.org/data](http://www.ppmi-info.org/data)) on October 1, 2020. PPMI samples were collected under a standardized protocol over 33 centers and includes clinical and imaging data as well as plasma and CSF samples. Study protocol and manuals are available online ([http://www.ppmi-info.org/study-design](http://www.ppmi-info.org/study-design)). Separate subcohorts of patients with Parkinson’s disease and their respective controls were enrolled following inclusion and exclusion criteria.11 One subcohort is comprised of recently diagnosed, drug-naïve, idiopathic Parkinson’s disease patients and healthy controls, while the second and third subcohorts are comprised of Parkinson’s disease patients, carriers of a severe *GBA* or *LRRK2* mutation, either Parkinson’s disease patients or unaffected controls. Parkinson’s disease patients from the genetic subcohorts had a higher disease duration, were partially under Parkinson’s disease medication (*n* = 203), were over-represented for individuals of Ashkenazi Jewish descent and differed by sex distribution from the idiopathic Parkinson’s disease patients (higher proportion of men among idiopathic). The study was approved by the Institutional Review Board at each site, and participants provided written informed consent. Genetic testing was done by the centralized PPMI genetic testing core. Non-manifesting carriers received pre-testing and post-testing genetic counselling by phone from certified genetic counsellors at the University of Indiana or site-qualified personnel. The *LRRK2* genetic testing battery includes G2019S and R1441G mutations. *GBA* genetic testing includes N370S (for all participants), and L483P, L444P, IVS2+1, and 84GG (for a subset of participants) mutations. Dual mutation carriers *(LRRK2* and *GBA)* were considered as *LRRK2* carriers for simplicity (*n* = 1). Six patients were diagnosed as idiopathic Parkinson’s disease at enrolment but were re-classified during follow-up (two patients were diagnosed as multiple system atrophy and four patients did not have a final diagnose but Parkinson’s disease had been excluded). These patients were removed from the analysis. Four patients were initially diagnosed as genetic Parkinson’s disease, but the diagnose changed to prodromal during follow up. These patients were considered as unaffected controls in their corresponding genetic subcohort (Five *LRRK2*+ unaffected controls and one *GBA*+ unaffected control). One subject originally classified as healthy, but later shown to have an unclear health status, was removed from the analysis. Subjects recruited into the subcohort of idiopathic Parkinson’s disease patients and healthy controls but identified as carriers of a severe *GBA* or *LRRK2* mutation, were moved to the corresponding genetic subcohort (*GBA n* = 15; *LRRK2 n* = 7). The genetic screening also detected *GBA* mutations of unknown or moderate risk: A459P, E365K, T408M. Carriers of these mutations were removed from analysis (*n* = 38). Finally, 10 carriers of a mutation in *SCNA* (eight Parkinson’s disease patients and two unaffected controls), were also removed due to lack of statistical power for analysis. The original data set was comprised of 1190 samples out of which 32 samples were pools, which were discarded for this study, and as described above, additional six Parkinson’s disease patients and one healthy control were removed due to change in diagnose, 38 subjects were removed due to non-severe *GBA* mutations and 10 patients were removed for being carriers of a mutation in *SCNA*. Hence, the final data set used for analysis was comprised of 1103 proteomic samples divided into three subcohorts: no mutation (idiopathic Parkinson’s disease patients and healthy controls with no severe mutation in *GBA* or *LRRK2*), *GBA*+ (Parkinson’s disease patient carriers of severe *GBA* mutations and unaffected controls carrying the same mutations) and *LRRK2*+ (Parkinson’s disease patients carriers of severe *LRRK2* mutations and unaffected controls carrying the same mutations). The exact composition is summarized in Table I. View this table: [Table I](http://medrxiv.org/content/early/2022/06/14/2022.03.09.22272131/T1) Table I Parkinson’s disease subcohorts We make an explicit distinction between “healthy” and “unaffected” controls, because there is evidence of prodromal pathophysiology in *LRRK2*+ and *GBA*+ controls when compared to healthy controls that are non-mutation carriers.12 ### Proteomics Proteomic profiling was performed using SomaScan® in a platform version that is proprietary to Novartis and includes 4785 SOMAmers® (Slow off-rate modified aptamers) targeting 4135 human proteins. SOMAmer levels were determined and standardized at SomaLogic Inc. (Boulder, US) including hybridization normalization (controls for variability in the readout of individual microarrays), plate scaling (accounts for plate-by-plate variation), median signal normalization (controls for total signal differences between individual samples) and calibration (removes the variation between assay runs within and across experiments). Relative fluorescence units are transformed to log2 scale, normalized to the median separately by dilution level across all plates. Finally, the data set is adjusted for batch effects between plates using an empirical Bayes method as implemented in the R package *sva*.13 ### Genetics PPMI whole genome sequencing results were lifted over to hg19 coordinates. Biallelic SNPs on autosomes were extracted. Standard GWAS quality control (QC) was applied at both individual and SNP level. 22 patients with outlying heterozygosity and 93 patients with high identity-by-descent were excluded after QC. 306031 SNPs were removed due to missing genotype and 35907596 SNPs removed due to minor allele count less than 20. Finally, 9743041 variants and 1264 subjects passed QC. ### pQTL calculation Among the 1264 subjects who passed QC for genetics and the 1103 who passed QC for proteomics, 804 subjects overlapped. Protein expression values were ranked and inverse normal transformed. For pQTL calculation, each protein level was regressed with each independent genetic variant (SNP; MAF > 0.05), adjusted for age, sex, subcohort, protein principal components 1-4 and genetic principal components 1-10 using R package *MatrixEQTL*.14 Cis-pQTLs were defined as SNPs located inside the +/- 1Mb region flanking the gene that encodes the given protein. A genome-wide threshold of *P* < 5×10−8 defined a significant cis-pQTL. ### Causal Analysis The two-sample mendelian randomization method implemented in the R package *TwoSampleMR**15,16* was applied to find causal proteins for Parkinson’s disease in CSF. Although PPMI is a well-controlled study with genetic data, to avoid weak instrumental variable bias and take the advantage of a larger Parkinson’s disease GWAS we relied on the meta-analysis from Nalls *et al*.17, which includes 17 datasets with Parkinson’s disease cases ranging from 363 to 33674 individuals, and healthy controls ranging from 165 to 449056 individuals. Instrumental variables were selected for each SNP with MAF > 0.05, F-statistics larger than 10 and significant cis-pQTL *P* < 5×10−8. Shared SNPs in both cis-pQTL and Parkinson’s disease GWAS were harmonized and then clumped using a linkage disequilibrium (LD) threshold of either *r**2* < 0.01 or *r**2* < 0.3. The more stringent LD threshold of r2<0.01 resulted in only one instrumental variable for most proteins, therefore the results we present are from the less stringent threshold of r2<0.3. The Wald ratio was used when only one instrument survived clumping, while the inverse variance weighted meta-analysis method was used when more than one instrumented SNP was available. Horizontal pleiotropy was tested using the R package *MRPRESSO*.18 Colocalization probability was calculated using the R package *coloc*.19 Default priors of *p**1* = 10−4, *p**2* = 10−4, and *p**12* = 10−5 were used, where p1 is the prior probability of a SNP being associated with Parkinson’s disease, p2 is the prior probability of a SNP being associated with CSF pQTL, and *p**12* is the prior probability of a SNP being associated with both Parkinson’s disease and CSF pQTL. We considered *PPH**4* > 0.75 as strong evidence for colocalization. *PPH**4* is the posterior probability of one shared SNP being associated with both Parkinson’s disease and CSF pQTL. ### Clinical variables The clinical assessment battery is described on the PPMI website ([http://www.ppmi-info.org](http://www.ppmi-info.org)). Parkinson’s disease status was assessed with the Unified Parkinson’s Disease Rating Scale in the revised version published by the Movement Disorder Society (MDS-UPDRS) scores 1, 2, and 3.20 Cognitive testing comprised screening with the Montreal Cognitive Assessment (MoCA).21 High resolution xy-weighted 3 tesla MRI was available for 545 Parkinson’s disease patients and 177 controls. Caudate, putamen and striatum thicknesses were calculated as the arithmetical mean between the right and left brain hemispheres. CSF was collected using standardized lumbar puncture procedures. Sample handling, shipment and storage were carried out as described in the PPMI biologics manual ([http://ppmi-info.org](http://ppmi-info.org)). Besides the SomaScan analysis described earlier, data from immunoassay kits were also used for measuring CSF total alpha-synuclein, amyloid-beta 1-42, total tau and phospho-tau (p-tau 181) protein as described previously.22,23 Phospho-tau was measured with the Elecsys® assay run on the fully automated Roche Cobas® system. Use of medications for Parkinson’s disease was recorded at each visit after baseline assessment. For simplicity, we used this as a binary variable (medication present / absent). ### Protein risk scores The Parkinson’s disease protein risk scores used in this study were generated as previously described 10. They were defined as “non-genetic Parkinson’s Disease-associated Proteomic Score (ngPD-ProS). ### Statistical analysis We first compared the protein profiles of Parkinson’s disease patients with controls within each subcohort (idiopathic, *GBA*+ and *LRRK2*+). A linear regression model was applied using the Bioconductor R package *limma*.24 The model included the following covariates: age, sex, study center and proteomic principal components 1-4. The genetic subcohorts also included levodopa treatment (yes/no) as a covariate in the model to exclude treatment effects. This was skipped for idiopathic patients as they were drug-naïve. Additionally, a linear model with an interaction term was tested. The interaction term was between the disease status (case /control) and the mutation status (no mutation / *GBA+* / *LRRK2+*). The covariates were the same as in the models above. Comparisons of clinical variables between endotypes of idiopathic patients were performed using a chi-squared test for categorical variables or a generalized linear model adjusted for age and sex for quantitative variables. To test differences in age between endotypes a Mann Whitney U-test was used. Predictive modeling for the idiopathic classes was performed using a partition tree with pruning as implemented in the R package *rpart*.25 The model was defined on a training set (70% of the idiopathic patients) and tested on an independent test set (30% of the idiopathic patients). The pruning was based upon a 10-fold cross validation, the default for *rpart*. All p-values were adjusted for multiple testing using false discovery rate (FDR). ### Cluster analysis Network analysis of the CSF proteome of idiopathic patients was carried out using the R package *WGCNA*.26 Co-expressed proteins (SOMAmers) were grouped into modules. Consensus clustering as implemented in the R package *ConsensusClusterPlus*27 used the SOMAmer modules to identify idiopathic patient subclasses. To avoid confounders being responsible for the differences between patient subclasses, network analysis was performed on the SOMAmer residuals of a linear regression on age, sex and study center. Heatmaps were generated using the R package *Heatplus*.28 ### Data availability The data used for this study is publicly available in the PPMI web page [https://www.ppmi-info.org/access-data-specimens/download-data](https://www.ppmi-info.org/access-data-specimens/download-data). The free access requires registration. The clinical data snapshot used here is kept under the tab “Archived PPMI data” >> “Publication Associated Archives” >> “2022-0001 Serrano-Fernandez: Parkinson’s Disease Proteogenomics (Version: 2022-05-18)”. The proteomic data is available under “Biospecimen” >> “Proteomic Analysis” >> “Project 151 Identification of proteins & protein networks & pQTL analysis in CSF x of 7 (Batch Corrected)” (7 files in total). The original adat files are also available under “Biospecimen” >> “Proteomic Analysis” >> “Project 151 Identification of proteins & protein networks & pQTL analysis in CSF - ADAT files”. ## Results ### Causal analysis Our analysis reported significant cis-pQTLs for 856 SOMAmers – corresponding to 744 unique proteins (Supplementary Table 5). From these 856 cis-pQTLs, we identified statistically significant evidence for causation for 68 proteins in CSF (Table 2). Out of those proteins, GPNMB, FCGR2A and FCGR2B also had a strong colocalization signal (see Methods), indicating the same SNP is both associated with protein level and Parkinson’s disease risk (Figure 2). The full tables with nominal *P* ≤ 0.05 are included in Supplementary Tables 1 and 2 for less or more stringent clumping, respectively. View this table: [Table 2](http://medrxiv.org/content/early/2022/06/14/2022.03.09.22272131/T2) Table 2 Parkinson’s disease causal proteins from CSF ![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/14/2022.03.09.22272131/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2022/06/14/2022.03.09.22272131/F2) Figure 2. Causal Analysis (best protein candidate). (**A**) Locus visualization of GPNMB pQTL hits suggest a strong association between GPNMB SOMAmer levels and its cis-SNPs. Colors indicate the linkage disequilibrium (LD) correlation of other SNPs with chr7:23294144 (rs858275). (**B**) Locus visualization of GWAS hits in the GPNMB locus for the risk of developing Parkinson’s disease. The y axis is the -log10 nominal *P*-value of the GWAS results. ### Differential protein expression in subcohorts of Parkinson’s disease patients To identify proteins differentially expressed in Parkinson’s disease, we compared SOMAmers in each of the subcohorts (*GBA+, LRRK2+* and idiopathic) to their corresponding controls. Our statistical analyses revealed six differentially expressed SOMAmers for *GBA+* patients, seven SOMAmers for *LRRK2+* patients and 23 SOMAmers for idiopathic patients. Directionality of the change and adjusted *P*-values for each of these markers are reported in Table 3. View this table: [Table 3](http://medrxiv.org/content/early/2022/06/14/2022.03.09.22272131/T3) Table 3 **SOMAmers differentially expressed between Parkinson’s disease patient vs. controls, divided by subcohort** For each subcohort, several identified markers confirmed previously reported proteins dysregulated in Parkinson’s disease. Interestingly, there was little overlap between proteins dysregulated in *GBA+, LRRK2+* and idiopathic subcohorts (SEMG2 and DLK1 were shared by *GBA+* and the idiopathic subcohort) though in each list there is a high percentage (4/6 in *GBA+*, 4/7 in *LRRK2+* and 10/23 in the idiopathic subcohort) of markers previously reported in relation to Parkinson’s disease (Table 3). Only one protein, CTSB, passed the FDR significance threshold for the interaction between disease status and mutation status. As shown in Table 3, CTSB was also differentially expressed in the *LRRK2+* subcohort. Additional analysis comparing all patient subcohorts with all controls, using subcohort membership (no mutation/*GBA*+/*LRRK2*+) and treatment status (yes/no) as additional covariates resulted in 129 SOMAmers, tagging 122 distinct proteins, passing FDR correction (Supplementary Table 3). ### Identification of subtypes of idiopathic Parkinson’s disease patients #### Identification of endotypes To determine if distinct endotypes were present in the idiopathic subcohort we performed a network analysis on the CSF proteome of idiopathic patients. Two modules of co-expressed proteins were identified. They comprised 889 and 600 SOMAmers, respectively. Applying consensus clustering on these two protein modules split the idiopathic subcohort (350 patients) in proteome-based patient endotypes 1 (85 patients) and 2 (165 patients) (Figure 3A). As seen in the tracking plot (Figure 3B) these endotypes suffer only negligible changes as the number of modeled subclasses increases. Moreover, a predictive model for the endotypes was built based on clinical parameters, avoiding the re-use of the same proteomic data involved in the definition of the endotypes. Patients with CSF phospho-tau ≥ 11pg/mL (as measured with the Elecsys® assay) were enriched for endotype 2, and patients with CSF phospho-tau < 11pg/mL were enriched for endotype 1 (Figure 3C). The model accuracy in the training (244 patients) and in the independent test (106 patients) sets, was 0.82 and 0.73, respectively. The estimated area under the curve (AUC) for the test set was 0.77. ![Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/14/2022.03.09.22272131/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2022/06/14/2022.03.09.22272131/F3) Figure 3. Idiopathic Parkinson’s disease endotypes. (**A**) Heatmap of z-scores of the protein values as measured with SomaScan, corresponding to the two modules identified using Weighted Gene Co-expression Network Analysis (WGCNA). The proteins in these modules are used for cluster analysis using Consensus Clustering, which retrieves two clusters (endotypes) of idiopathic patients. Patients are shown in the x-axis, separated by endotype, while proteins are shown in the y-axis, separated by module. (**B**) Tracking plot depicting how the idiopathic patients are assigned to specific endotypes by Consensus Clustering as the number of potential endotypes increases. (**C**) Partition Tree predicting endotype membership of the idiopathic patients based on clinical variables only. One node suffices to separate patients into endotypes based on phospho-tau levels (p-tau) in CSF as measured with a clinical assay, the cutoff being 11 pg/mL. Endotypes did not significantly differ in age or sex. There were neither significant differences in caudate, striatum or putamen thickness, nor in UPDRS score parts II and III or MoCA scores. Significant differences were found for UPDRS score part I (higher for endotype 2; *P* = 0.044), as well as CSF levels of amyloid beta (lower for endotype 2; *P* = 1.95×10−15.), phospho-tau (lower for endotype 2; *P* =1.25×10−15.), total tau (lower for endotype 2; *P* < 2×10−16.) and alpha-synuclein (lower for endotype 2; *P* < 2×10−16). Endotypes also showed significant differences in ngPD-ProS scores (higher for endotype 2; *P* = 0.034). ### Proteins differentially expressed in endotypes (CSF SomaScan) To identify the unique proteins significantly dysregulated in each endotype, a linear model was used to identify the differences between each of these two endotypes and the healthy controls to the idiopathic subcohort. For endotype 1, five markers were significantly different compared to the control group (CNTFR, LPO, MMP10, RIPK2 and VEGFA). LPO, RIPK2 and VEGFA were also part of the differences between healthy controls and the whole idiopathic group (see above). Endotype 2, however, showed 200 differentially expressed SOMAmers, 197 unique proteins (see Supplementary Table 4). Among those proteins, AK1, CCL14, FRZB, GPI, HAMP, LPO, NETO1, PTPRR, RAB31, RELT, RIPK2, ROBO3, RSPO4, SHANK1, SPINK9, VEGFA and VIP were dysregulated for the whole idiopathic group as compared to healthy controls. Differentially expressed SOMAmers between endotypes added up to 155, 153 unique proteins (see Supplementary Table 4). ## Discussion ### Causal analysis Independently of its etiology (genetic or idiopathic), multiple cellular and metabolic alterations (e.g. iron metabolism), inflammation, and oxidative stress underly neurodegeneration in Parkinson’s disease. In this regard, the identification of causal proteins not only enhances its understanding, but also assists the search for druggable targets. This study identified 68 CSF proteins as causal and GPNMB, FCGR2A, FCGR2B and CTSB were among the top ones. GPNMB - expressed by myeloid cells - is found at high levels in the substantia nigra of Parkinson’s disease patients.29 It has been suggested that GPNMB modulates immune response 30,31 having primarily a protective role.32-35 CTSB cleaves alpha-synuclein fibrils with the potential for decreasing alpha-synuclein aggregation.36 FCGR2A and FCGR2B are involved in phagocytosis and modulate inflammatory responses.37 Supporting our findings, independent reports suggest that GPNMB, CTSB and FCGR2A are causal to Parkinson’s disease.38-40 FCGR2B - aggregated alpha-synuclein interaction inhibits its microglial phagocytosis, promoting neurodegeneration.41 However, blockade of FCGR2B in neurons may suppress Lewy body-like inclusion body formation.42 ### Genetic patients GBA and LRRK2 are enzymes involved in ceramide metabolism and lysosomal degradation of aggregated alpha-synuclein. Genetic variants of the *GBA* and *LRRK2* genes are known risk factors for the development of Parkinson’s disease and dementia associated with accumulation of Lewy bodies.43-45 To better understand the differences between diseased and unaffected mutation carriers we compared their CSF proteomes (see Table 3). Protein differences between *GBA*+ patients and unaffected *GBA*+ controls are related to brain dopaminergic neurons (DLK1), dopamine metabolism (CALCA, GCH1) and inflammation effector cells (IL17A). All these proteins have been previously linked to Parkinson’s disease, though not specifically to *GBA*+. DLK1 is coupled to both tyrosine hydroxylase expression and neurotrophic signaling.46 Significantly lower levels of DLK1 in *GBA*+ patients suggest increased degeneration of dopaminergic neurons. GCH1 and CALCA are involved in nigrostriatal dopamine synthesis,47 and dopamine release and metabolism,48,49 respectively. We found both elevated in *GBA*+ patients. *GCH1* variants contribute to the risk and earlier age-at-onset of Parkinson’s disease.50 Moreover, CALCA was reported as significantly elevated in the CSF of depressed Parkinson’s disease patients relative to major depressive disorder patients.51 One of the inflammatory mechanisms proposed for neurodegeneration in Parkinson’s disease involves Th17 cells.52 Our study found increased IL17A levels in *GBA+* patients. It has been recently proposed that elevated IL17A plays a key role in neurodegenerative diseases.53 Among the CSF proteins differentially expressed between *LRRK2*+ patients and unaffected *LRRK2*+ controls, several stood out for their relevance in Parkinson’s disease: ARSA, SMPD1, CTSB and TENM4.ARSA is a lysosomal chaperone that prevents alpha-synuclein aggregation, secretion and cell-to-cell propagation.54 We found ARSA elevated in *LRRK2*+ patients, which could be interpreted as a protection mechanism to prevent the formation of alpha-synuclein aggregates. CTSB (causal protein, see causal analysis) and SMPD1 play important roles in Parkinson’s disease autophagy and lysosomal degradation processes.55 *CTSB* and *SMPD1* genetic variants are known to be associated with Parkinson’s disease risk.55 In this study, higher CTSB and SMPD1 levels in *LRRK2*+ patients indicate dysregulation of autophagy-endolysosomal pathway and potentially increased macroautophagy.56,57 Brain TENM4 is involved in axon guidance and myelination.58 In our study, it was significantly reduced in *LRRK2*+ patients. Loss-of-function and missense variants in *TENM4* are associated with early onset Parkinson’s disease and essential tremor, a potential risk factor for developing Parkinson’s disease.60-59 ### Idiopathic patients Worldwide, ∼90% of Parkinson’s disease cases are idiopathic. Proteome differences between idiopathic patients and controls comprise not only markers of inflammation, but also of mitochondrial dysfunction / oxidative stress, iron metabolism and other pathological processes (see Table 3). The AK1 kinase is expressed by neurons and astrocytes.60 At advanced Parkinson’s disease stages, AK1 is downregulated in the substantia nigra probably due to mitochondrial dysfunction and dopaminergic neuronal death.60 The observed elevation of CSF AK1 levels may be associated with Parkinson’s disease progression stages, frontal cortex primary alteration or compensation of altered purine metabolism.60 SHANK may be regulated by the mitochondrial kinase PINK1, for which variants are known to be causal for Parkinson’s disease.61 It has been reported that knockdown in neurons of *PINK* decreases PSD95 and SHANK1.61 SHANK1 was decreased in idiopathic patients suggesting impaired synaptic plasticity. TXN promotes cell proliferation, protection against oxidative stress and anti-apoptotic functions in the brain, which makes it a good candidate for a neurodegeneration marker. Here we find TXN decreased in Parkinson’s disease patients. Iron dysregulation is associated with oxidative stress and lipid peroxidation.62 It has been proposed that LPO - heme peroxidase - in the substantia nigra is involved in neurodegeneration.63 Lower LPO levels in idiopathic patients found here contrasts with previously reported elevated CSF LPO levels.63 The high CSF levels of HAMP could help explaining this discrepancy. HAMP reduces iron accumulation and neuroinflammation by decreasing mitochondrial dysfunction, oxidative stress, and ultimately dopamine neuronal loss.64 Moreover, HAMP overexpression – as seen here – promotes alpha-synuclein clearance through autophagy.65 It has been also reported that dopamine and levodopa reduce LPO levels.63 Given that the idiopathic patients recruited were drug-naïve early Parkinson’s disease patients, dopamine levels in this subpopulation may have helped maintaining low levels of CSF LPO. The glucose metabolism enzyme GPI was elevated in idiopathic patients and may be protective. Its overexpression in dopaminergic neurons protects against alpha-synuclein-induced neurotoxicity.66 As seen in *GBA*+ patients, reduced DLK1 levels in idiopathic patients may suggest neurodegeneration. Lower levels of RAB31, a small GTPase involved in exosome biogenesis67 and potentially in alpha-synuclein spreading in Parkinson’s disease68 were observed, too. RIPK2, a LRRK2 substrate, is lower in idiopathic patients, which matches the fact that LRRK2 deficiency leads to reduced activation of RIPK2.69 The neurotrophic factor VEGFA is neuroprotective and has genetic variant associated with Parkinson’s disease risk.70 VIP enhances striatal plasticity and prevents dopaminergic cell loss in parkinsonian rats.71 VEGFA and VIP at lower levels in idiopathic patients may reflect ongoing neurodegeneration. ### Heterogeneity of idiopathic patients Disease heterogeneity challenges the development of disease modifying therapies. In this study, we used proteomic data to stratify idiopathic patients into clinically relevant endotypes. The lack of differences in the CSF levels of causal proteins, suggests that endotype molecular differences may be downstream from causal effects and affect the specific characteristics of Parkinson’s disease phenotypes rather than Parkinson’s disease risk. And yet, those molecular differences can affect progression as reflected by differences in ngPD-ProS and might be the basis of personalized therapy approaches. The endotype robustness was reflected by the identified patient cluster stability (Figure 3B) and the high performance of the endotype predictive model (see Figure 3C). The fact that CSF phospho-tau levels sufficed to predict endotypes, suggests that targeted assessment of CSF proteins may be appropriate for idiopathic patient stratification in a clinical setting. It is worth noting that, while in Tsukita *et al*. 10 the ngPD-ProS perform similarly for both the idiopathic and the combined genetic subcohorts, here we see little overlap in the differentially expressed proteins between the idiopathic and each of the genetic subcohorts separately. This apparent discrepancy could just be reflecting a common CSF proteomic signature rising from the combination of genetic subcohorts even though each genetic subcohort may have different etiology. Supporting this hypothesis, when we compared the proteomes of the combined genetic subcohorts (*GBA+* and *LRRK2+*) with their respective controls (data not shown) we found seven proteins shared with the 24 in the ngPD-ProS (i.e. DLK1, LPO, NEFH, RIPK2, SEMG2, VIP and LRFN2). These seven proteins intersected almost completely with the ten proteins that overlapped between the 24 proteins in the ngPD-ProS and the proteins differentially expressed between idiopathic patients and controls (i.e. DLK1, LPO, NEFH, RIPK2, SEMG2, VIP, CCL14, HAMP, RSPO4, and TXN). This similarity between the idiopathic subcohort and the combination of the genetic subcohorts could help explaining the apparent discrepancy mentioned above. In summary, we: i) identified causal proteins for Parkinson’s disease, ii) assessed CSF proteome differences in Parkinson’s disease patients of genetic and idiopathic etiology, and. iii) stratified idiopathic patients into robust subtypes. Our findings not only contribute to the identification of new therapeutic targets but also to shaping personalized medicine in CNS neurodegeneration. ## Supporting information Supplementary Material [[supplements/272131_file03.pdf]](pending:yes) ## Data Availability The results of the current study are publicly available on the LONI database and all data produced in the present work are contained in the manuscript ## Funding The study was funded by the Novartis Institutes for Biomedical Research and Merck. Protein measurements were performed at SomaLogic. ## Competing interests The authors report no competing interests. ## Supplementary material Supplementary material is available at *Brain* online. ## Acknowledgements The authors would like to thank Rose Case (Indiana University Genetics Biobank, Indiana University School of Medicine, Indianapolis, USA) for her key role in sample management and Myung Shin (Genome and Biomarker Sciences. Merck Exploratory Science Center. Cambridge, USA), Karla Gonzalez (Department of Medical and Molecular Genetics, Indiana University School of Medicine. Indianapolis, USA) and Faraz Faghri (Laboratory of Neurogenetics, National Institute on Aging, National Institutes of Health. Bethesda, USA) for constructive feedback. ## Abbreviations AUC : Area Under the Curve FDR : False Discovery Rate GWAS : Genome-wide Association Study LD : Linkage Disequilibrium MAF : Major Allele Frequency MDS : Movement Disorder Society MoCA : Montreal Cognitive Assessment ngPD-ProS : non-genetic Parkinson’s Disease-associated Proteomic Score PPMI : Parkinson Progression Marker Initiative pQTL : Proteomic Quantitative Trait Locus QC : Quality Control SOMAmer : Slow Off-rate Modified Aptamer SNP : Single Nucleotide Polymorphism UPDRS : Unified Parkinson’s Disease Rating Scale WGCNA : Weighted Gene Co-expression Network Analysis * Received March 9, 2022. * Revision received June 14, 2022. * Accepted June 14, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Group GBDNDC. Global, regional, and national burden of neurological disorders during 1990-2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet Neurol. 2017;2017(11):877–897. 2. 2.Blauwendraat C, Nalls MA, Singleton AB. The genetic architecture of Parkinson’s disease. Lancet Neurol. 2020;2020(2):170–178. 3. 3.Bonifati V. Genetics of Parkinson’s disease--state of the art, 2013. Parkinsonism Relat Disord. 2014;20 Suppl 1:S23–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1353-8020(13)70009-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24262182&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F14%2F2022.03.09.22272131.atom) 4. 4.Kestenbaum M, Alcalay RN. Clinical Features of LRRK2 Carriers with Parkinson’s Disease. Adv Neurobiol. 2017;14:31–48. 5. 5.Thaler A, Bregman N, Gurevich T, et al. Parkinson’s disease phenotype is influenced by the severity of the mutations in the GBA gene. Parkinsonism Relat Disord. 2018;55:45–49. 6. 6.Qian E, Huang Y. Subtyping of Parkinson’s Disease - Where Are We Up To? Aging Dis. 2019;2019(5):1130–1139. 7. 7.Ma LY, Chan P, Gu ZQ, Li FF, Feng T. Heterogeneity among patients with Parkinson’s disease: cluster analysis and genetic association. J Neurol Sci. 2015;351(1-2):41–45. 8. 8.Parkinson Progression Marker I. The Parkinson Progression Marker Initiative (PPMI). Prog Neurobiol. 2011;2011(4):629–35. 9. 9.Gold L, Ayers D, Bertino J, et al. Aptamer-based multiplexed proteomic technology for biomarker discovery. PLoS One. 2010;2010(12):e15004. 10. 10.Tsukita K, Sakamaki-Tsukita H, Kaiser S, et al. Changes in the cerebrospinal fluid proteome precede and stratify the course of Parkinson’s Disease. medRxiv. 2022:2022.06.08.22276035. 11. 11.Marek K, Chowdhury S, Siderowf A, et al. The Parkinson’s progression markers initiative (PPMI) - establishing a PD biomarker cohort. Ann Clin Transl Neurol. 2018;2018(12):1460–1477. 12. 12.Simuni T, Uribe L, Cho HR, et al. Clinical and dopamine transporter imaging characteristics of non-manifest LRRK2 and GBA mutation carriers in the Parkinson’s Progression Markers Initiative (PPMI): a cross-sectional study. Lancet Neurol. 2020;2020(1):71–80. 13. 13.Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;2012(6):882–3. 14. 14.Shabalin AA. Matrix eQTL: ultra fast eQTL analysis via large matrix operations. Bioinformatics. 2012;2012(10):1353–8. 15. 15.Hemani G, Zheng J, Elsworth B, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7 16. 16.Hemani G, Tilling K, Davey Smith G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;2017(11):e1007081. 17. 17.Nalls MA, Blauwendraat C, Vallerga CL, 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;2019(12):1091–1102. 18. 18.Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;2018(5):693–698. 19. 19.Giambartolomei C, Vukcevic D, Schadt EE, et al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. 2014;2014(5):e1004383. 20. 20.Goetz CG, Tilley BC, Shaftman SR, et al. Movement Disorder Society-sponsored revision of the Unified Parkinson’s Disease Rating Scale (MDS-UPDRS): scale presentation and clinimetric testing results. Mov Disord. 2008;2008(15):2129–70. 21. 21.Nasreddine ZS, Phillips NA, Bedirian V, et al. The Montreal Cognitive Assessment, MoCA: a brief screening tool for mild cognitive impairment. J Am Geriatr Soc. 2005;2005(4):695–9. 22. 22.Kang JH, Irwin DJ, Chen-Plotkin AS, et al. Association of cerebrospinal fluid beta-amyloid 1-42, T-tau, P-tau181, and alpha-synuclein levels with clinical features of drug-naive patients with early Parkinson disease. JAMA Neurol. 2013;2013(10):1277–87. 23. 23.Mollenhauer B, Caspell-Garcia CJ, Coffey CS, et al. Longitudinal analyses of cerebrospinal fluid alpha-Synuclein in prodromal and early Parkinson’s disease. Mov Disord. 2019;2019(9):1354–1364. 24. 24.Ritchie ME, Phipson B, Wu D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;2015(7):e47. 25. 25.rpart: Recursive Partitioning and Regression Trees. R package version 4.1-15. 2019. [https://github.com/bethatkinson/rpart](https://github.com/bethatkinson/rpart) 26. 26.Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-9-559&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19114008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F14%2F2022.03.09.22272131.atom) 27. 27.Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010;2010(12):1572–3. 28. 28.Heatplus: Heatmaps with row and/or column covariates and colored clusters. R package. Bioconductor; 2015. 29. 29.Moloney EB, Moskites A, Ferrari EJ, Isacson O, Hallett PJ. The glycoprotein GPNMB is selectively elevated in the substantia nigra of Parkinson’s disease patients and increases after lysosomal stress. Neurobiol Dis. 2018;120:1–11. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.nbd.2018.08.013&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30149180&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F14%2F2022.03.09.22272131.atom) 30. 30.Ripoll VM, Irvine KM, Ravasi T, Sweet MJ, Hume DA. Gpnmb is induced in macrophages by IFN-gamma and lipopolysaccharide and acts as a feedback regulator of proinflammatory responses. J Immunol. 2007;2007(10):6557–66. 31. 31.Chung JS, Sato K, Dougherty, II, Cruz PD, Jr.., Ariizumi K. DC-HIL is a negative regulator of T lymphocyte activation. Blood. 2007;2007(10):4320–7. 32. 32.Nagahara Y, Shimazawa M, Tanaka H, et al. Glycoprotein nonmetastatic melanoma protein B ameliorates skeletal muscle lesions in a SOD1G93A mouse model of amyotrophic lateral sclerosis. J Neurosci Res. 2015;2015(10):1552–66. 33. 33.Nagahara Y, Shimazawa M, Ohuchi K, et al. GPNMB ameliorates mutant TDP-43-induced motor neuron cell death. J Neurosci Res. 2017;2017(8):1647–1665. 34. 34.Saade M, Araujo de Souza G, Scavone C, Kinoshita PF. The Role of GPNMB in Inflammation. Front Immunol. 2021;12:674739. 35. 35.Neal ML, Boyle AM, Budge KM, Safadi FF, Richardson JR. The glycoprotein GPNMB attenuates astrocyte inflammatory responses through the CD44 receptor. J Neuroinflammation. 2018;2018(1):73. 36. 36.McGlinchey RP, Lee JC. Cysteine cathepsins are essential in lysosomal degradation of alpha-synuclein. Proc Natl Acad Sci U S A. 2015;2015(30):9322–7. 37. 37.Bournazos S, Gupta A, Ravetch JV. The role of IgG Fc receptors in antibody-dependent enhancement. Nat Rev Immunol. 2020;2020(10):633–643. 38. 38.Kia DA, Zhang D, Guelfi S, et al. Identification of Candidate Parkinson Disease Genes by Integrating Genome-Wide Association Study, Expression, and Epigenetic Data Sets. JAMA Neurol. 2021;2021(4):464–472. 39. 39.Storm CS, Kia DA, Almramhi MM, et al. Finding genetically-supported drug targets for Parkinson’s disease using Mendelian randomization of the druggable genome. Nat Commun. 2021;2021(1):7342. 40. 40.Schilder BM, Raj T. Fine-mapping of Parkinson’s disease susceptibility loci identifies putative causal variants. Hum Mol Genet. 2022;2022(6):888–900. 41. 41.Choi YR, Kang SJ, Kim JM, et al. FcgammaRIIB mediates the inhibitory effect of aggregated alpha-synuclein on microglial phagocytosis. Neurobiol Dis. 2015;83:90–9. 42. 42.Choi YR, Cha SH, Kang SJ, Kim JB, Jou I, Park SM. Prion-like Propagation of alpha-Synuclein Is Regulated by the FcgammaRIIB-SHP-1/2 Signaling Pathway in Neurons. Cell Rep. 2018;2018(1):136–148. 43. 43.Bastien J, Menon S, Messa M, Nyfeler B. Molecular targets and approaches to restore autophagy and lysosomal capacity in neurodegenerative disorders. Mol Aspects Med. 2021;82:101018. 44. 44.Pang SY, Lo RCN, Ho PW, et al. LRRK2, GBA and their interaction in the regulation of autophagy: implications on therapeutics in Parkinson’s disease. Transl Neurodegener. 2022;2022(1):5. 45. 45.Rivero-Rios P, Romo-Lozano M, Fasiczka R, Naaldijk Y, Hilfiker S. LRRK2-Related Parkinson’s Disease Due to Altered Endolysosomal Biology With Variable Lewy Body Pathology: A Hypothesis. Front Neurosci. 2020;14:556. 46. 46.Bossers K, Meerhoff G, Balesar R, et al. Analysis of gene expression in Parkinson’s disease: possible involvement of neurotrophic support and axon guidance in dopaminergic cell death. Brain Pathol. 2009;2009(1):91–107. 47. 47.Kapatos G. The neurobiology of tetrahydrobiopterin biosynthesis: a model for regulation of GTP cyclohydrolase I gene transcription within nigrostriatal dopamine neurons. IUBMB Life. 2013;2013(4):323–33. 48. 48.Deutch AY, Roth RH. Calcitonin gene-related peptide in the ventral tegmental area: selective modulation of prefrontal cortical dopamine metabolism. Neurosci Lett. 1987;1987(2):169–74. 49. 49.Drumheller A, Menard D, Fournier A, Jolicoeur FB. Neurochemical effects of CGRP. Ann N Y Acad Sci. 1992;657:546–8. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1637121&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F14%2F2022.03.09.22272131.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992JL49900078&link_type=ISI) 50. 50.Pan HX, Zhao YW, Mei JP, et al. GCH1 variants contribute to the risk and earlier age-at-onset of Parkinson’s disease: a two-cohort case-control study. Transl Neurodegener. 2020;2020(1):31. 51. 51.Svenningsson P, Palhagen S, Mathe AA. Neuropeptide Y and Calcitonin Gene-Related Peptide in Cerebrospinal Fluid in Parkinson’s Disease with Comorbid Depression versus Patients with Major Depressive Disorder. Front Psychiatry. 2017;8:102. 52. 52.Bolte AC, Lukens JR. Th17 Cells in Parkinson’s Disease: The Bane of the Midbrain. Cell Stem Cell. 2018;2018(1):5–6. 53. 53.Chen J, Liu X, Zhong Y. Interleukin-17A: The Key Cytokine in Neurodegenerative Diseases. Front Aging Neurosci. 2020;12:566922. 54. 54.Angelopoulou E, Paudel YN, Villa C, Piperi C. Arylsulfatase A (ASA) in Parkinson’s Disease: From Pathogenesis to Biomarker Potential. Brain Sci. 2020;10(10) 55. 55.Navarro-Romero A, Montpeyo M, Martinez-Vicente M. The Emerging Role of the Lysosome in Parkinson’s Disease. Cells. 2020;9(11) 56. 56.Manzoni C. The LRRK2-macroautophagy axis and its relevance to Parkinson’s disease. Biochem Soc Trans. 2017;2017(1):155–162. 57. 57.Plowey ED, Cherra SJ, 3rd., Liu YJ, Chu CT. Role of autophagy in G2019S-LRRK2-associated neurite shortening in differentiated SH-SY5Y cells. J Neurochem. 2008;2008(3):1048–56. 58. 58.Hor H, Francescatto L, Bartesaghi L, et al. Missense mutations in TENM4, a regulator of axon guidance and central myelination, cause essential tremor. Hum Mol Genet. 2015;2015(20):5677–86. 59. 59.Algarni M, Fasano A. The overlap between Essential tremor and Parkinson disease. Parkinsonism Relat Disord. 2018;46 Suppl 1:S101–S104. 60. 60.Garcia-Esparcia P, Hernandez-Ortega K, Ansoleaga B, Carmona M, Ferrer I. Purine metabolism gene deregulation in Parkinson’s disease. Neuropathol Appl Neurobiol. 2015;2015(7):926–40. 61. 61.Hernandez CJ, Baez-Becerra C, Contreras-Zarate MJ, Arboleda H, Arboleda G. PINK1 Silencing Modifies Dendritic Spine Dynamics of Mouse Hippocampal Neurons. J Mol Neurosci. 2019;2019(4):570–579. 62. 62.Ndayisaba A, Kaindlstorfer C, Wenning GK. Iron in Neurodegeneration - Cause or Consequence? Front Neurosci. 2019;13:180. 63. 63.Fernandez-Espejo E, Rodriguez de Fonseca F, Suarez J, Martin de Pablos A. Cerebrospinal fluid lactoperoxidase level is enhanced in idiopathic Parkinson’s disease, and correlates with levodopa equivalent daily dose. Brain Res. 2021;1761:147411. 64. 64.Liang T, Qian ZM, Mu MD, Yung WH, Ke Y. Brain Hepcidin Suppresses Major Pathologies in Experimental Parkinsonism. iScience. 2020;2020(7):101284. 65. 65.Urrutia PJ, Borquez DA, Nunez MT. Inflaming the Brain with Iron. Antioxidants (Basel). 2021;10(1) 66. 66.Knight AL, Yan X, Hamamichi S, et al. The glycolytic enzyme, GPI, is a functionally conserved modifier of dopaminergic neurodegeneration in Parkinson’s models. Cell Metab. 2014;2014(1):145–57. 67. 67.Wei D, Zhan W, Gao Y, et al. RAB31 marks and controls an ESCRT-independent exosome pathway. Cell Res. 2021;2021(2):157–177. 68. 68.Kumar R, Donakonda S, Muller SA, Botzel K, Hoglinger GU, Koeglsperger T. FGF2 Affects Parkinson’s Disease-Associated Molecular Networks Through Exosomal Rab8b/Rab31. Front Genet. 2020;11:572058. 69. 69.Yan R, Liu Z. LRRK2 enhances Nod1/2-mediated inflammatory cytokine production by promoting Rip2 phosphorylation. Protein Cell. 2017;2017(1):55–66. 70. 70.Wu Y, Zhang Y, Han X, Li X, Xue L, Xie A. Association of VEGF gene polymorphisms with sporadic Parkinson’s disease in Chinese Han population. Neurol Sci. 2016;2016(12):1923–1929. 71. 71.Korkmaz O, Ay H, Ulupinar E, Tuncel N. Vasoactive intestinal peptide enhances striatal plasticity and prevents dopaminergic cell loss in Parkinsonian rats. J Mol Neurosci. 2012;2012(3):565–73.