The Parkinson's Disease Genome-Wide Association Study Locus Browser

BACKGROUND

Parkinson's disease (PD) is a neurodegenerative disease with an often complex component identifiable by genome-wide association studies. The most recent large-scale PD genome-wide association studies have identified more than 90 independent risk variants for PD risk and progression across more than 80 genomic regions. One major challenge in current genomics is the identification of the causal gene(s) and variant(s) at each genome-wide association study locus. The objective of the current study was to create a tool that would display data for relevant PD risk loci and provide guidance with the prioritization of causal genes and potential mechanisms at each locus.

METHODS

We included all significant genome-wide signals from multiple recent PD genome-wide association studies including themost recent PD risk genome-wide association study, age-at-onset genome-wide association study, progression genome-wide association study, and Asian population PD risk genome-wide association study. We gathered data for all genes 1 Mb up and downstream of each variant to allow users to assess which gene(s) are most associated with the variant of interest based on a set of self-ranked criteria. Multiple databases were queried for each gene to collect additional causal data.

RESULTS

We created a PD genome-wide association study browser tool (https://pdgenetics.shinyapps.io/GWASBrowser/) to assist the PD research community with the prioritization of genes for follow-up functional studies to identify potential therapeutic targets.

CONCLUSIONS

Our PD genome-wide association study browser tool provides users with a useful method of identifying potential causal genes at all known PD risk loci from large-scale PD genome-wide association studies. We plan to update this tool with new relevant data as sample sizes increase and new PD risk loci are discovered. © 2020 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

LRRK2 mediates microglial neurotoxicity via NFATc2 in rodent models of synucleinopathies

Synucleinopathies are neurodegenerative disorders characterized by abnormal α-synuclein deposition that include Parkinson's disease, dementia with Lewy bodies, and multiple system atrophy. The pathology of these conditions also includes neuronal loss and neuroinflammation. Neuron-released α-synuclein has been shown to induce neurotoxic, proinflammatory microglial responses through Toll-like receptor 2, but the molecular mechanisms involved are poorly understood. Here, we show that leucine-rich repeat kinase 2 (LRRK2) plays a critical role in the activation of microglia by extracellular α-synuclein. Exposure to α-synuclein was found to enhance LRRK2 phosphorylation and activity in mouse primary microglia. Furthermore, genetic and pharmacological inhibition of LRRK2 markedly diminished α-synuclein-mediated microglial neurotoxicity via lowering of tumor necrosis factor-α and interleukin-6 expression in mouse cultures. We determined that LRRK2 promoted a neuroinflammatory cascade by selectively phosphorylating and inducing nuclear translocation of the immune transcription factor nuclear factor of activated T cells, cytoplasmic 2 (NFATc2). NFATc2 activation was seen in patients with synucleinopathies and in a mouse model of synucleinopathy, where administration of an LRRK2 pharmacological inhibitor restored motor behavioral deficits. Our results suggest that modulation of LRRK2 and its downstream signaling mediator NFATc2 might be therapeutic targets for treating synucleinopathies.

Differences in the Presentation and Progression of Parkinson's Disease by Sex

BACKGROUND

Previous studies reported various symptoms of Parkinson's disease (PD) associated with sex. Some were conflicting or confirmed in only one study.

OBJECTIVES

We examined sex associations to PD phenotypes cross-sectionally and longitudinally in large-scale data.

METHODS

We tested 40 clinical phenotypes, using longitudinal, clinic-based patient cohorts, consisting of 5946 patients, with a median follow-up of 3.1 years. For continuous outcomes, we used linear regressions at baseline to test sex-associated differences in presentation, and linear mixed-effects models to test sex-associated differences in progression. For binomial outcomes, we used logistic regression models at baseline and Cox regression models for survival analyses. We adjusted for age, disease duration, and medication use. In the secondary analyses, data from 17 719 PD patients and 7588 non-PD participants from an online-only, self-assessment PD cohort were cross-sectionally evaluated to determine whether the sex-associated differences identified in the primary analyses were consistent and unique to PD.

RESULTS

Female PD patients had a higher risk of developing dyskinesia early during the follow-up period, with a slower progression in activities of daily living difficulties, and a lower risk of developing cognitive impairments compared with male patients. The findings in the longitudinal, clinic-based cohorts were mostly consistent with the results of the online-only cohort.

CONCLUSIONS

We observed sex-associated contributions to PD heterogeneity. These results highlight the necessity of future research to determine the underlying mechanisms and importance of personalized clinical management. © 2020 International Parkinson and Movement Disorder Society.

Large-scale pathway specific polygenic risk and transcriptomic community network analysis identifies novel functional pathways in Parkinson disease

Polygenic inheritance plays a central role in Parkinson disease (PD). A priority in elucidating PD etiology lies in defining the biological basis of genetic risk. Unraveling how risk leads to disruption will yield disease-modifying therapeutic targets that may be effective. Here, we utilized a high-throughput and hypothesis-free approach to determine biological processes underlying PD using the largest currently available cohorts of genetic and gene expression data from International Parkinson's Disease Genetics Consortium (IPDGC) and the Accelerating Medicines Partnership-Parkinson's disease initiative (AMP-PD), among other sources. We applied large-scale gene-set specific polygenic risk score (PRS) analyses to assess the role of common variation on PD risk focusing on publicly annotated gene sets representative of curated pathways. We nominated specific molecular sub-processes underlying protein misfolding and aggregation, post-translational protein modification, immune response, membrane and intracellular trafficking, lipid and vitamin metabolism, synaptic transmission, endosomal-lysosomal dysfunction, chromatin remodeling and apoptosis mediated by caspases among the main contributors to PD etiology. We assessed the impact of rare variation on PD risk in an independent cohort of whole-genome sequencing data and found evidence for a burden of rare damaging alleles in a range of processes, including neuronal transmission-related pathways and immune response. We explored enrichment linked to expression cell specificity patterns using single-cell gene expression data and demonstrated a significant risk pattern for dopaminergic neurons, serotonergic neurons, hypothalamic GABAergic neurons, and neural progenitors. Subsequently, we created a novel way of building de novo pathways by constructing a network expression community map using transcriptomic data derived from the blood of PD patients, which revealed functional enrichment in inflammatory signaling pathways, cell death machinery related processes, and dysregulation of mitochondrial homeostasis. Our analyses highlight several specific promising pathways and genes for functional prioritization and provide a cellular context in which such work should be done.

Genetic modifiers of risk and age at onset in GBA associated Parkinson's disease and Lewy body dementia

Parkinson's disease is a genetically complex disorder. Multiple genes have been shown to contribute to the risk of Parkinson's disease, and currently 90 independent risk variants have been identified by genome-wide association studies. Thus far, a number of genes (including SNCA, LRRK2, and GBA) have been shown to contain variability across a spectrum of frequency and effect, from rare, highly penetrant variants to common risk alleles with small effect sizes. Variants in GBA, encoding the enzyme glucocerebrosidase, are associated with Lewy body diseases such as Parkinson's disease and Lewy body dementia. These variants, which reduce or abolish enzymatic activity, confer a spectrum of disease risk, from 1.4- to >10-fold. An outstanding question in the field is what other genetic factors that influence GBA-associated risk for disease, and whether these overlap with known Parkinson's disease risk variants. Using multiple, large case-control datasets, totalling 217 165 individuals (22 757 Parkinson's disease cases, 13 431 Parkinson's disease proxy cases, 622 Lewy body dementia cases and 180 355 controls), we identified 1691 Parkinson's disease cases, 81 Lewy body dementia cases, 711 proxy cases and 7624 controls with a GBA variant (p.E326K, p.T369M or p.N370S). We performed a genome-wide association study and analysed the most recent Parkinson's disease-associated genetic risk score to detect genetic influences on GBA risk and age at onset. We attempted to replicate our findings in two independent datasets, including the personal genetics company 23andMe, Inc. and whole-genome sequencing data. Our analysis showed that the overall Parkinson's disease genetic risk score modifies risk for disease and decreases age at onset in carriers of GBA variants. Notably, this effect was consistent across all tested GBA risk variants. Dissecting this signal demonstrated that variants in close proximity to SNCA and CTSB (encoding cathepsin B) are the most significant contributors. Risk variants in the CTSB locus were identified to decrease mRNA expression of CTSB. Additional analyses suggest a possible genetic interaction between GBA and CTSB and GBA p.N370S induced pluripotent cell-derived neurons were shown to have decreased cathepsin B expression compared to controls. These data provide a genetic basis for modification of GBA-associated Parkinson's disease risk and age at onset, although the total contribution of common genetics variants is not large. We further demonstrate that common variability at genes implicated in lysosomal function exerts the largest effect on GBA associated risk for disease. Further, these results have implications for selection of GBA carriers for therapeutic interventions.

Validation of Serum Neurofilament Light Chain as a Biomarker of Parkinson's Disease Progression

Background:

The objective of this study was to assess neurofilament light chain as a Parkinson’s disease biomarker.

Methods:

We quantified neurofilament light chain in 2 independent cohorts: (1) longitudinal cerebrospinal fluid samples from the longitudinal de novo Parkinson’s disease cohort and (2) a large longitudinal cohort with serum samples from Parkinson’s disease, other cognate/neurodegenerative disorders, healthy controls, prodromal conditions, and mutation carriers.

Results:

In the Parkinson’s Progression Marker Initiative cohort, mean baseline serum neurofilament light chain was higher in Parkinson’s disease patients (13 ± 7.2 pg/mL) than in controls (12 ± 6.7 pg/mL), P = 0.0336. Serum neurofilament light chain increased longitudinally in Parkinson’s disease patients versus controls (P < 0.01). Motor scores were positively associated with neurofilament light chain, whereas some cognitive scores showed a negative association.

Conclusions:

Neurofilament light chain in serum samples is increased in Parkinson’s disease patients versus healthy controls, increases over time and with age, and correlates with clinical measures of Parkinson’s disease severity. Although the specificity of neurofilament light chain for Parkinson’s disease is low, it is the first blood-based biomarker candidate that could support disease stratification of Parkinson’s disease versus other cognate/neurodegenerative disorders, track clinical progression, and possibly assess responsiveness to neuroprotective treatments. However, use of neurofilament light chain as a biomarker of response to neuroprotective interventions remains to be assessed.

Analysis of DNM3 and VAMP4 as genetic modifiers of LRRK2 Parkinson's disease

The LRRK2 gene has rare (p.G2019S) and common risk variants for Parkinson’s disease (PD). DNM3 has previously been reported as a genetic modifier of the age at onset in PD patients carrying the LRRK2 p.G2019S mutation. We analyzed this effect in a new cohort of LRRK2 p.G2019S heterozygotes (n = 724) and meta-analyzed our data with previously published data (n = 754). VAMP4 is in close proximity to DNM3, and was associated with PD in a recent study, so it is possible that variants in this gene may be important. We also analyzed the effect of VAMP4 rs11578699 on LRRK2 penetrance. Our analysis of DNM3 in previously unpublished data does not show an effect on age at onset in LRRK2 p.G2019S carriers; however, the inter-study heterogeneity may indicate ethnic or population-specific effects of DNM3. There was no evidence for linkage disequilibrium between DNM3 and VAMP4. Analysis of sporadic patients stratified by the risk variant LRRK2 rs10878226 indicates a possible interaction between common variation in LRRK2 and VAMP4 in disease risk.

•

Analysis of DNM3 in previously unpublished data does not show an interaction between DNM3 and LRRK2 G2019S for AAO.

•

There was no evidence for linkage disequilibrium between DNM3 rs2421947 and VAMP4 rs11578699.

•

Analysis of sporadic PD patients stratified by the PD risk variant rs10878226 indicates a possible interaction between LRRK2 and VAMP4.

Comprehensive assessment of PINK1 variants in Parkinson's disease

Multiple genes have been associated with monogenic Parkinson's disease and Parkinsonism syndromes. Mutations in PINK1 (PARK6) have been shown to result in autosomal recessive early-onset Parkinson's disease. In the past decade, several studies have suggested that carrying a single heterozygous PINK1 mutation is associated with increased risk for Parkinson's disease. Here, we comprehensively assess the role of PINK1 variants in Parkinson's disease susceptibility using several large data sets totalling 376,558 individuals including 13,708 cases with Parkinson's disease and 362,850 control subjects. After combining these data, we did not find evidence to support a role for heterozygous PINK1 mutations as a robust risk factor for Parkinson's disease.

Functionalization of the TMEM175 p.M393T variant as a risk factor for Parkinson disease

Multiple genome-wide association studies (GWAS) in Parkinson disease (PD) have identified a signal at chromosome 4p16.3; however, the causal variant has not been established for this locus. Deep investigation of the region resulted in one identified variant, the rs34311866 missense SNP (p.M393T) in TMEM175, which is 20 orders of magnitude more significant than any other SNP in the region. Because TMEM175 is a lysosomal gene that has been shown to influence α-synuclein phosphorylation and autophagy, the p.M393T variant is an attractive candidate, and we have examined its effect on TMEM175 protein and PD-related biology. After knocking down each of the genes located under the GWAS peak via multiple shRNAs, only TMEM175 was found to consistently influence accumulation of phosphorylated α-synuclein (p-α-syn). Examination of the p.M393T variant showed effects on TMEM175 function that were intermediate between the wild-type (WT) and knockout phenotypes, with reduced regulation of lysosomal pH in response to starvation and minor changes in clearance of autophagy substrates, reduced lysosomal localization, and increased accumulation of p-α-syn. Finally, overexpression of WT TMEM175 protein reduced p-α-syn, while overexpression of the p.M393T variant resulted in no change in α-synuclein phosphorylation. These results suggest that the main signal in the chromosome 4p16.3 PD risk locus is driven by the TMEM175 p.M393T variant. Modulation of TMEM175 may impact α-synuclein biology and therefore may be a rational therapeutic strategy for PD.

LRRK2-mediated microglial activation via NFATc2: a novel mechanism of neurotoxic inflammation in synucleinopathies

Neuroinflammation is a key pathological feature of a set of neurological disorders called synucleinopathies, which are characterized by abnormal a-synuclein (a-syn) dynamics and include Parkinson’s disease (PD), dementia with Lewy Bodies (DLB), and multiple system atrophy (MSA). Microglia, the resident macrophages of the central nervous system, are critical for orchestrating the inflammatory response in these diseases. In particular, microglia can directly bind neuron-released a-syn via receptors such as Toll-like receptor 2 (TLR2) to adopt a reactive phenotype. Here we explored the pathogenic role of Leucine-rich repeat kinase 2 (LRRK2), a PD-associated gene also known to be highly expressed in T cells, in microglial activation. We observed that neuron-released a-syn enhances LRRK2 kinase activity in microglia in a TLR2-dependent manner to induce the release of neurotoxic cytokines, TNFa and IL-6. We further determined that LRRK2 directly modulated nuclear translocation of activated T-cells cytoplasmic 2 (NFATc2) through multi-site phosphorylation. As such, modulating the LRRK2-NFATc2 cascade may be a promising therapeutic strategy for ameliorating neuroinflammation in synucleinopathies.

Genetics of Parkinson's disease: An introspection of its journey towards precision medicine

A substantial proportion of risk for Parkinson's disease (PD) is driven by genetics. Progress in understanding the genetic basis of PD has been significant. So far, highly-penetrant rare genetic alterations in SNCA, LRRK2, VPS35, PRKN, PINK1, DJ-1 and GBA have been linked with typical familial PD and common genetic variability at 90 loci have been linked to risk for PD. In this review, we outline the journey thus far of PD genetics, highlighting how significant advances have improved our knowledge of the genetic basis of PD risk, onset and progression. Despite remarkable progress, our field has yet to unravel how genetic risk variants disrupt biological pathways and molecular networks underlying the pathobiology of the disease. We highlight that currently identified genetic risk factors only represent a fraction of the likely genetic risk for PD. Identifying the remaining genetic risk will require us to diversify our efforts, performing genetic studies across different ancestral groups. This work will inform us on the varied genetic basis of disease across populations and also aid in fine mapping discovered loci. If we are able to take this course, we foresee that genetic discoveries in PD will directly influence our ability to predict disease and aid in defining etiological subtypes, critical steps for the implementation of precision medicine for PD.

Disease modification and biomarker development in Parkinson disease: Revision or reconstruction?

A fundamental question in advancing Parkinson disease (PD) research is whether it represents one disorder or many. Does each genetic PD inform a common pathobiology or represent a unique entity? Do the similarities between genetic and idiopathic forms of PD outweigh the differences? If aggregates of α-synuclein in Lewy bodies and Lewy neurites are present in most (α-synucleinopathies), are they also etiopathogenically significant in each (α-synuclein pathogenesis)? Does it matter that postmortem studies in PD have demonstrated that mixed protein-aggregate pathology is the rule and pure α-synucleinopathy the exception? Should we continue to pursue convergent biomarkers that are representative of the diverse whole of PD or subtype-specific, divergent biomarkers, present in some but absent in most? Have clinical trials that failed to demonstrate efficacy of putative disease-modifying interventions been true failures (shortcomings of the hypotheses, which should be rejected) or false failures (shortcomings of the trials; hypotheses should be preserved)? Each of these questions reflects a nosologic struggle between the lumper's clinicopathologic model that embraces heterogeneity of one disease and the splitter's focus on a pathobiology-specific set of diseases. Most important, even if PD is not a single disorder, can advances in biomarkers and disease modification be revised to concentrate on pathologic commonalities in large, clinically defined populations? Or should our efforts be reconstructed to focus on smaller subgroups of patients, distinguished by well-defined molecular characteristics, regardless of their phenotypic classification? Will our clinical trial constructs be revised to target larger and earlier, possibly even prodromal, cohorts? Or should our trials efforts be reconstructed to target smaller but molecularly defined presymptomatic or postsymptomatic cohorts? At the Krembil Knowledge Gaps in Parkinson's Disease Symposium, the tentative answers to these questions were discussed, informed by the failures and successes of the fields of breast cancer and cystic fibrosis.

Fine-Mapping of SNCA in Rapid Eye Movement Sleep Behavior Disorder and Overt Synucleinopathies

OBJECTIVE

Rapid eye movement sleep behavior disorder (RBD) is a prodromal synucleinopathy, as >80% will eventually convert to overt synucleinopathy. We performed an in-depth analysis of the SNCA locus to identify RBD-specific risk variants.

METHODS

Full sequencing and genotyping of SNCA was performed in isolated/idiopathic RBD (iRBD, n = 1,076), Parkinson disease (PD, n = 1,013), dementia with Lewy bodies (DLB, n = 415), and control subjects (n = 6,155). The iRBD cases were diagnosed with RBD prior to neurodegeneration, although some have since converted. A replication cohort from 23andMe of PD patients with probable RBD (pRBD) was also analyzed (n = 1,782 cases; n = 131,250 controls). Adjusted logistic regression models and meta-analyses were performed. Effects on conversion rate were analyzed in 432 RBD patients with available data using Kaplan-Meier survival analysis.

RESULTS

A 5'-region SNCA variant (rs10005233) was associated with iRBD (odds ratio [OR] = 1.43, p = 1.1E-08), which was replicated in pRBD. This variant is in linkage disequilibrium (LD) with other 5' risk variants across the different synucleinopathies. An independent iRBD-specific suggestive association (rs11732740) was detected at the 3' of SNCA (OR = 1.32, p = 4.7E-04, not statistically significant after Bonferroni correction). Homozygous carriers of both iRBD-specific SNPs were at highly increased risk for iRBD (OR = 5.74, p = 2E-06). The known top PD-associated variant (3' variant rs356182) had an opposite direction of effect in iRBD compared to PD.

INTERPRETATION

There is a distinct pattern of association at the SNCA locus in RBD as compared to PD, with an opposite direction of effect at the 3' of SNCA. Several 5' SNCA variants are associated with iRBD and with pRBD in overt synucleinopathies. ANN NEUROL 2020;87:584-598.

Penetrance of Parkinson's Disease in LRRK2 p.G2019S Carriers Is Modified by a Polygenic Risk Score

BACKGROUND

Although the leucine-rich repeat kinase 2 p.G2019S mutation has been demonstrated to be a strong risk factor for PD, factors that contribute to penetrance among carriers, other than aging, have not been well identified.

OBJECTIVES

To evaluate whether a cumulative genetic risk identified in the recent genome-wide study is associated with penetrance of PD among p.G2019S mutation carriers.

METHODS

We included p.G2019S heterozygote carriers with European ancestry in three genetic cohorts in which the mutation carriers with and without PD were selectively recruited. We also included the carriers from two data sets: one from a case-control setting without selection of mutation carriers and the other from a population sampling. Associations between polygenic risk score constructed from 89 variants reported recently and PD were tested and meta-analyzed. We also explored the interaction of age and PRS.

RESULTS

After excluding eight homozygotes, 833 p.G2019S heterozygote carriers (439 PD and 394 unaffected) were analyzed. Polygenic risk score was associated with a higher penetrance of PD (odds ratio: 1.34; 95% confidence interval: [1.09, 1.64] per +1 standard deviation; P = 0.005). In addition, associations with polygenic risk score and penetrance were stronger in the younger participants (main effect: odds ratio 1.28 [1.04, 1.58] per +1 standard deviation; P = 0.022; interaction effect: odds ratio 0.78 [0.64, 0.94] per +1 standard deviation and + 10 years of age; P = 0.008).

CONCLUSIONS

Our results suggest that there is a genetic contribution for penetrance of PD among p.G2019S carriers. These results have important etiological consequences and potential impact on the selection of subjects for clinical trials. © 2020 International Parkinson and Movement Disorder Society.

Genomewide association study of Parkinson's disease clinical biomarkers in 12 longitudinal patients' cohorts

BACKGROUND

Several reports have identified different patterns of Parkinson's disease progression in individuals carrying missense variants in GBA or LRRK2 genes. The overall contribution of genetic factors to the severity and progression of Parkinson's disease, however, has not been well studied.

OBJECTIVES

To test the association between genetic variants and the clinical features of Parkinson's disease on a genomewide scale.

METHODS

We accumulated individual data from 12 longitudinal cohorts in a total of 4093 patients with 22,307 observations for a median of 3.81 years. Genomewide associations were evaluated for 25 cross-sectional and longitudinal phenotypes. Specific variants of interest, including 90 recently identified disease-risk variants, were also investigated post hoc for candidate associations with these phenotypes.

RESULTS

Two variants were genomewide significant. Rs382940(T>A), within the intron of SLC44A1, was associated with reaching Hoehn and Yahr stage 3 or higher faster (hazard ratio 2.04 [1.58-2.62]; P value = 3.46E-8). Rs61863020(G>A), an intergenic variant and expression quantitative trait loci for α-2A adrenergic receptor, was associated with a lower prevalence of insomnia at baseline (odds ratio 0.63 [0.52-0.75]; P value = 4.74E-8). In the targeted analysis, we found 9 associations between known Parkinson's risk variants and more severe motor/cognitive symptoms. Also, we replicated previous reports of GBA coding variants (rs2230288: p.E365K; rs75548401: p.T408M) being associated with greater motor and cognitive decline over time, and an APOE E4 tagging variant (rs429358) being associated with greater cognitive deficits in patients.

CONCLUSIONS

We identified novel genetic factors associated with heterogeneity of Parkinson's disease. The results can be used for validation or hypothesis tests regarding Parkinson's disease. © 2019 International Parkinson and Movement Disorder Society.

Identification of novel risk loci, causal insights, and heritable risk for Parkinson's disease: a meta-analysis of genome-wide association studies

BACKGROUND

Genome-wide association studies (GWAS) in Parkinson's disease have increased the scope of biological knowledge about the disease over the past decade. We aimed to use the largest aggregate of GWAS data to identify novel risk loci and gain further insight into the causes of Parkinson's disease.

METHODS

We did a meta-analysis of 17 datasets from Parkinson's disease GWAS available from European ancestry samples to nominate novel loci for disease risk. These datasets incorporated all available data. We then used these data to estimate heritable risk and develop predictive models of this heritability. We also used large gene expression and methylation resources to examine possible functional consequences as well as tissue, cell type, and biological pathway enrichments for the identified risk factors. Additionally, we examined shared genetic risk between Parkinson's disease and other phenotypes of interest via genetic correlations followed by Mendelian randomisation.

FINDINGS

Between Oct 1, 2017, and Aug 9, 2018, we analysed 7·8 million single nucleotide polymorphisms in 37 688 cases, 18 618 UK Biobank proxy-cases (ie, individuals who do not have Parkinson's disease but have a first degree relative that does), and 1·4 million controls. We identified 90 independent genome-wide significant risk signals across 78 genomic regions, including 38 novel independent risk signals in 37 loci. These 90 variants explained 16-36% of the heritable risk of Parkinson's disease depending on prevalence. Integrating methylation and expression data within a Mendelian randomisation framework identified putatively associated genes at 70 risk signals underlying GWAS loci for follow-up functional studies. Tissue-specific expression enrichment analyses suggested Parkinson's disease loci were heavily brain-enriched, with specific neuronal cell types being implicated from single cell data. We found significant genetic correlations with brain volumes (false discovery rate-adjusted p=0·0035 for intracranial volume, p=0·024 for putamen volume), smoking status (p=0·024), and educational attainment (p=0·038). Mendelian randomisation between cognitive performance and Parkinson's disease risk showed a robust association (p=8·00 × 10^-7).

INTERPRETATION

These data provide the most comprehensive survey of genetic risk within Parkinson's disease to date, to the best of our knowledge, by revealing many additional Parkinson's disease risk loci, providing a biological context for these risk factors, and showing that a considerable genetic component of this disease remains unidentified. These associations derived from European ancestry datasets will need to be followed-up with more diverse data.

FUNDING

The National Institute on Aging at the National Institutes of Health (USA), The Michael J Fox Foundation, and The Parkinson's Foundation (see appendix for full list of funding sources).

The Parkinson's Disease Mendelian Randomization Research Portal

BACKGROUND

Mendelian randomization is a method for exploring observational associations to find evidence of causality.

OBJECTIVE

To apply Mendelian randomization between risk factors/phenotypic traits (exposures) and PD in a large, unbiased manner, and to create a public resource for research.

METHODS

We used two-sample Mendelian randomization in which the summary statistics relating to single-nucleotide polymorphisms from 5,839 genome-wide association studies of exposures were used to assess causal relationships with PD. We selected the highest-quality exposure genome-wide association studies for this report (n = 401). For the disease outcome, summary statistics from the largest published PD genome-wide association studies were used. For each exposure, the causal effect on PD was assessed using the inverse variance weighted method, followed by a range of sensitivity analyses. We used a false discovery rate of 5% from the inverse variance weighted analysis to prioritize exposures of interest.

RESULTS

We observed evidence for causal associations between 12 exposures and risk of PD. Of these, nine were effects related to increasing adiposity and decreasing risk of PD. The remaining top three exposures that affected PD risk were tea drinking, time spent watching television, and forced vital capacity, but these may have been biased and were less convincing. Other exposures at nominal statistical significance included inverse effects of smoking and alcohol.

CONCLUSIONS

We present a new platform which offers Mendelian randomization analyses for a total of 5,839 genome-wide association studies versus the largest PD genome-wide association studies available (https://pdgenetics.shinyapps.io/MRportal/). Alongside, we report further evidence to support a causal role for adiposity on lowering the risk of PD. © 2019 The Authors. Movement Disorders published by Wiley Periodicals, Inc. on behalf of International Parkinson and Movement Disorder Society.

Genetic variability and potential effects on clinical trial outcomes: perspectives in Parkinson's disease

BACKGROUND

Classical randomisation of clinical trial patients creates a source of genetic variance that may be contributing to the high failure rate seen in neurodegenerative disease trials. Our objective was to quantify genetic difference between randomised trial arms and determine how imbalance can affect trial outcomes.

METHODS

5851 patients with Parkinson's disease of European ancestry data and two simulated virtual cohorts based on public data were used. Data were resampled at different sizes for 1000 iterations and randomly assigned to the two arms of a simulated trial. False-negative and false-positive rates were estimated using simulated clinical trials, and per cent difference in genetic risk score (GRS) and allele frequency was calculated to quantify variance between arms.

RESULTS

5851 patients with Parkinson's disease (mean (SD) age, 61.02 (12.61) years; 2095 women (35.81%)) as well as simulated patients from virtually created cohorts were used in the study. Approximately 90% of the iterations had at least one statistically significant difference in individual risk SNPs between each trial arm. Approximately 5%-6% of iterations had a statistically significant difference between trial arms in mean GRS. For significant iterations, the average per cent difference for mean GRS between trial arms was 130.87%, 95% CI 120.89 to 140.85 (n=200). Glucocerebrocidase (GBA) gene-only simulations see an average 18.86%, 95% CI 18.01 to 19.71 difference in GRS scores between trial arms (n=50). When adding a drug effect of -0.5 points in MDS-UPDRS per year at n=50, 33.9% of trials resulted in false negatives.

CONCLUSIONS

Our data support the hypothesis that within genetically unmatched clinical trials, genetic heterogeneity could confound true therapeutic effects as expected. Clinical trials should undergo pretrial genetic adjustment or, at the minimum, post-trial adjustment and analysis for failed trials.

Frequency of Loss of Function Variants in LRRK2 in Parkinson Disease

Importance

Pathogenic variants in LRRK2 are a relatively common genetic cause of Parkinson disease (PD). Currently, the molecular mechanism underlying disease is unknown, and gain and loss of function (LOF) models of pathogenesis have been postulated. LRRK2 variants are reported to result in enhanced phosphorylation of substrates and increased cell death. However, the double knockout of Lrrk2 and its homologue Lrrk1 results in neurodegeneration in a mouse model, suggesting that disease may occur by LOF. Because LRRK2 inhibitors are currently in development as potential disease-modifying treatments in PD, it is critical to determine whether LOF variants in LRRK2 increase or decrease the risk of PD.

Objective

To determine whether LRRK1 and LRRK2 LOF variants contribute to the risk of developing PD.

Design, Setting, and Participants

To determine the prevailing mechanism of LRRK2-mediated disease in human populations, next-generation sequencing data from a large case-control cohort (>23 000 individuals) was analyzed for LOF variants in LRRK1 and LRRK2. Data were generated at 5 different sites and 5 different data sets, including cases with clinically diagnosed PD and neurologically normal control individuals. Data were collected from 2012 through 2017.

Main Outcomes and Measures

Frequencies of LRRK1 and LRRK2 LOF variants present in the general population and compared between cases and controls.

Results

Among 11 095 cases with PD and 12 615 controls, LRRK1 LOF variants were identified in 0.205% of cases and 0.139% of controls (odds ratio, 1.48; SE, 0.571; 95% CI, 0.45-4.44; P = .49) and LRRK2 LOF variants were found in 0.117% of cases and 0.087% of controls (odds ratio, 1.48; SE, 0.431; 95% CI, 0.63-3.50; P = .36). All association tests suggested lack of association between LRRK1 or LRRK2 variants and PD. Further analysis of lymphoblastoid cell lines from several heterozygous LOF variant carriers found that, as expected, LRRK2 protein levels are reduced by approximately half compared with wild-type alleles.

Conclusions and Relevance

Together these findings indicate that haploinsufficiency of LRRK1 or LRRK2 is neither a cause of nor protective against PD. Furthermore, these results suggest that kinase inhibition or allele-specific targeting of mutant LRRK2 remain viable therapeutic strategies in PD.

SNCA and mTOR Pathway Single Nucleotide Polymorphisms Interact to Modulate the Age at Onset of Parkinson's Disease

BACKGROUND

Single nucleotide polymorphisms (SNPs) in the α-synuclein (SNCA) gene are associated with differential risk and age at onset (AAO) of both idiopathic and Leucine-rich repeat kinase 2 (LRRK2)-associated Parkinson's disease (PD). Yet potential combinatory or synergistic effects among several modulatory SNPs for PD risk or AAO remain largely underexplored.

OBJECTIVES

The mechanistic target of rapamycin (mTOR) signaling pathway is functionally impaired in PD. Here we explored whether SNPs in the mTOR pathway, alone or by epistatic interaction with known susceptibility factors, can modulate PD risk and AAO.

METHODS

Based on functional relevance, we selected a total of 64 SNPs mapping to a total of 57 genes from the mTOR pathway and genotyped a discovery series cohort encompassing 898 PD patients and 921 controls. As a replication series, we screened 4170 PD and 3014 controls available from the International Parkinson's Disease Genomics Consortium.

RESULTS

In the discovery series cohort, we found a 4-loci interaction involving STK11 rs8111699, FCHSD1 rs456998, GSK3B rs1732170, and SNCA rs356219, which was associated with an increased risk of PD (odds ratio = 2.59, P < .001). In addition, we also found a 3-loci epistatic combination of RPTOR rs11868112 and RPS6KA2 rs6456121 with SNCA rs356219, which was associated (odds ratio = 2.89; P < .0001) with differential AAO. The latter was further validated (odds ratio = 1.56; P = 0.046-0.047) in the International Parkinson's Disease Genomics Consortium cohort.

CONCLUSIONS

These findings indicate that genetic variability in the mTOR pathway contributes to SNCA effects in a nonlinear epistatic manner to modulate differential AAO in PD, unraveling the contribution of this cascade in the pathogenesis of the disease. © 2019 International Parkinson and Movement Disorder Society.

Author Correction: Genetic meta-analysis of diagnosed Alzheimer's disease identifies new risk loci and implicates Aβ, tau, immunity and lipid processing

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Blood Leukocyte DNA Methylation Predicts Risk of Future Myocardial Infarction and Coronary Heart Disease

BACKGROUND

DNA methylation is implicated in coronary heart disease (CHD), but current evidence is based on small, cross-sectional studies. We examined blood DNA methylation in relation to incident CHD across multiple prospective cohorts.

METHODS

Nine population-based cohorts from the United States and Europe profiled epigenome-wide blood leukocyte DNA methylation using the Illumina Infinium 450k microarray, and prospectively ascertained CHD events including coronary insufficiency/unstable angina, recognized myocardial infarction, coronary revascularization, and coronary death. Cohorts conducted race-specific analyses adjusted for age, sex, smoking, education, body mass index, blood cell type proportions, and technical variables. We conducted fixed-effect meta-analyses across cohorts.

RESULTS

Among 11 461 individuals (mean age 64 years, 67% women, 35% African American) free of CHD at baseline, 1895 developed CHD during a mean follow-up of 11.2 years. Methylation levels at 52 CpG (cytosine-phosphate-guanine) sites were associated with incident CHD or myocardial infarction (false discovery rate<0.05). These CpGs map to genes with key roles in calcium regulation (ATP2B2, CASR, GUCA1B, HPCAL1), and genes identified in genome- and epigenome-wide studies of serum calcium (CASR), serum calcium-related risk of CHD (CASR), coronary artery calcified plaque (PTPRN2), and kidney function (CDH23, HPCAL1), among others. Mendelian randomization analyses supported a causal effect of DNA methylation on incident CHD; these CpGs map to active regulatory regions proximal to long non-coding RNA transcripts.

CONCLUSION

Methylation of blood-derived DNA is associated with risk of future CHD across diverse populations and may serve as an informative tool for gaining further insight on the development of CHD.

Genetic risk of Parkinson disease and progression:: An analysis of 13 longitudinal cohorts

OBJECTIVE

To determine if any association between previously identified alleles that confer risk for Parkinson disease and variables measuring disease progression.

METHODS

We evaluated the association between 31 risk variants and variables measuring disease progression. A total of 23,423 visits by 4,307 patients of European ancestry from 13 longitudinal cohorts in Europe, North America, and Australia were analyzed.

RESULTS

We confirmed the importance of GBA on phenotypes. GBA variants were associated with the development of daytime sleepiness (p.N370S: hazard ratio [HR] 3.28 [1.69-6.34]) and possible REM sleep behavior (p.T408M: odds ratio 6.48 [2.04-20.60]). We also replicated previously reported associations of GBA variants with motor/cognitive declines. The other genotype-phenotype associations include an intergenic variant near LRRK2 and the faster development of motor symptom (Hoehn and Yahr scale 3.0 HR 1.33 [1.16-1.52] for the C allele of rs76904798) and an intronic variant in PMVK and the development of wearing-off effects (HR 1.66 [1.19-2.31] for the C allele of rs114138760). Age at onset was associated with TMEM175 variant p.M393T (-0.72 [-1.21 to -0.23] in years), the C allele of rs199347 (intronic region of GPNMB, 0.70 [0.27-1.14]), and G allele of rs1106180 (intronic region of CCDC62, 0.62 [0.21-1.03]).

CONCLUSIONS

This study provides evidence that alleles associated with Parkinson disease risk, in particular GBA variants, also contribute to the heterogeneity of multiple motor and nonmotor aspects. Accounting for genetic variability will be a useful factor in understanding disease course and in minimizing heterogeneity in clinical trials.

Parkinson's disease age at onset genome-wide association study: Defining heritability, genetic loci, and α-synuclein mechanisms

BACKGROUND

Increasing evidence supports an extensive and complex genetic contribution to PD. Previous genome-wide association studies (GWAS) have shed light on the genetic basis of risk for this disease. However, the genetic determinants of PD age at onset are largely unknown.

OBJECTIVES

To identify the genetic determinants of PD age at onset.

METHODS

Using genetic data of 28,568 PD cases, we performed a genome-wide association study based on PD age at onset.

RESULTS

We estimated that the heritability of PD age at onset attributed to common genetic variation was ∼0.11, lower than the overall heritability of risk for PD (∼0.27), likely, in part, because of the subjective nature of this measure. We found two genome-wide significant association signals, one at SNCA and the other a protein-coding variant in TMEM175, both of which are known PD risk loci and a Bonferroni-corrected significant effect at other known PD risk loci, GBA, INPP5F/BAG3, FAM47E/SCARB2, and MCCC1. Notably, SNCA, TMEM175, SCARB2, BAG3, and GBA have all been shown to be implicated in α-synuclein aggregation pathways. Remarkably, other well-established PD risk loci, such as GCH1 and MAPT, did not show a significant effect on age at onset of PD.

CONCLUSIONS

Overall, we have performed the largest age at onset of PD genome-wide association studies to date, and our results show that not all PD risk loci influence age at onset with significant differences between risk alleles for age at onset. This provides a compelling picture, both within the context of functional characterization of disease-linked genetic variability and in defining differences between risk alleles for age at onset, or frank risk for disease. © 2019 International Parkinson and Movement Disorder Society.