Heritability Enrichment Implicates Microglia in Parkinson's Disease Pathogenesis.

OBJECTIVE: Understanding how different parts of the immune system contribute to pathogenesis in Parkinson's disease is a burning challenge with important therapeutic implications. We studied enrichment of common variant heritability for Parkinson's disease stratified by immune and brain cell types. METHODS: We used summary statistics from the most recent meta-analysis of genomewide association studies in Parkinson's disease and partitioned heritability using linkage disequilibrium score regression, stratified for specific cell types, as defined by open chromatin regions. We also validated enrichment results using a polygenic risk score approach and intersected disease-associated variants with epigenetic data and expression quantitative loci to nominate and explore a putative microglial locus. RESULTS: We found significant enrichment of Parkinson's disease risk heritability in open chromatin regions of microglia and monocytes. Genomic annotations overlapped substantially between these 2 cell types, and only the enrichment signal for microglia remained significant in a joint model. We present evidence suggesting P2RY12, a key microglial gene and target for the antithrombotic agent clopidogrel, as the likely driver of a significant Parkinson's disease association signal on chromosome 3. INTERPRETATION: Our results provide further support for the importance of immune mechanisms in Parkinson's disease pathogenesis, highlight microglial dysregulation as a contributing etiological factor, and nominate a targetable microglial gene candidate as a pathogenic player. Immune processes can be modulated by therapy, with potentially important clinical implications for future treatment in Parkinson's disease. ANN NEUROL 2021;89:942-951.

Common X-Chromosome Variants Are Associated with Parkinson Disease Risk.

OBJECTIVE: The objective of this study was to identify genetic variants on the X-chromosome associated with Parkinson disease (PD) risk. METHODS: We performed an X-chromosome-wide association study (XWAS) of PD risk by meta-analyzing results from sex-stratified analyses. To avoid spurious associations, we designed a specific harmonization pipeline for the X-chromosome and focused on a European ancestry sample. We included 11,142 cases, 280,164 controls, and 5,379 proxy cases, based on parental history of PD. Additionally, we tested the association of significant variants with (1) PD risk in an independent replication with 1,561 cases and 2,465 controls and (2) putamen volume in 33,360 individuals from the UK Biobank. RESULTS: In the discovery meta-analysis, we identified rs7066890 (odds ratio [OR] = 1.10, 95% confidence interval [CI] = 1.06-1.14, p = 2.2 × 10-9 ), intron of GPM6B, and rs28602900 (OR = 1.10, 95% CI = 1.07-1.14, p = 1.6 × 10-8 ) in a high gene density region including RPL10, ATP6A1, FAM50A, and PLXNA3. The rs28602900 association with PD was replicated (OR = 1.16, 95% CI = 1.03-1.30, p = 0.016) and shown to colocalize with a significant expression quantitative locus (eQTL) regulating RPL10 expression in the putamen and other brain tissues in the Genotype-Tissue Expression Project. Additionally, the rs28602900 locus was found to be associated with reduced brain putamen volume. No results reached genome-wide significance in the sex-stratified analyses. INTERPRETATION: We report the first XWAS of PD and identify 2 genome-wide significant loci. The rs28602900 association was replicated in an independent PD dataset and showed concordant effects in its association with putamen volume. Critically, rs26802900 is a significant eQTL of RPL10. These results support a role for ribosomal proteins in PD pathogenesis and show that the X-chromosome contributes to PD genetic risk. ANN NEUROL 2021;90:22-34.

Investigation of Autosomal Genetic Sex Differences in Parkinson's Disease.

OBJECTIVE: Parkinson's disease (PD) is a complex neurodegenerative disorder. Men are on average ~ 1.5 times more likely to develop PD compared to women with European ancestry. Over the years, genomewide association studies (GWAS) have identified numerous genetic risk factors for PD, however, it is unclear whether genetics contribute to disease etiology in a sex-specific manner. METHODS: In an effort to study sex-specific genetic factors associated with PD, we explored 2 large genetic datasets from the International Parkinson's Disease Genomics Consortium and the UK Biobank consisting of 13,020 male PD cases, 7,936 paternal proxy cases, 89,660 male controls, 7,947 female PD cases, 5,473 maternal proxy cases, and 90,662 female controls. We performed GWAS meta-analyses to identify distinct patterns of genetic risk contributing to disease in male versus female PD cases. RESULTS: In total, 19 genomewide significant regions were identified and no sex-specific effects were observed. A high genetic correlation between the male and female PD GWAS were identified (rg = 0.877) and heritability estimates were identical between male and female PD cases (~ 20%). INTERPRETATION: We did not detect any significant genetic differences between male or female PD cases. Our study does not support the notion that common genetic variation on the autosomes could explain the difference in prevalence of PD between males and females cases at least when considering the current sample size under study. Further studies are warranted to investigate the genetic architecture of PD explained by X and Y chromosomes and further evaluate environmental effects that could potentially contribute to PD etiology in male versus female patients. ANN NEUROL 2021;90:41-48.

Rapid eye movement sleep behaviour disorder: Past, present, and future.

This manuscript presents an overview of REM sleep behaviour disorder (RBD) with a special focus on European contributions. After an introduction examining the history of the disorder, we address the pathophysiological and clinical aspects, as well as the diagnostic issues. Further, implications of RBD diagnosis and biomarkers are discussed. Contributions of European researchers to this field are highlighted.

Targeted sequencing of Parkinson's disease loci genes highlights SYT11, FGF20 and other associations.

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

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.

Genetic, Structural, and Functional Evidence Link TMEM175 to Synucleinopathies.

OBJECTIVE: The TMEM175/GAK/DGKQ locus is the 3rd strongest risk locus in genome-wide association studies of Parkinson disease (PD). We aimed to identify the specific disease-associated variants in this locus, and their potential implications. METHODS: Full sequencing of TMEM175/GAK/DGKQ followed by genotyping of specific associated variants was performed in PD (n = 1,575) and rapid eye movement sleep behavior disorder (RBD) patients (n = 533) and in controls (n = 1,583). Adjusted regression models and a meta-analysis were performed. Association between variants and glucocerebrosidase (GCase) activity was analyzed in 715 individuals with available data. Homology modeling, molecular dynamics simulations, and lysosomal localization experiments were performed on TMEM175 variants to determine their potential effects on structure and function. RESULTS: Two coding variants, TMEM175 p.M393T (odds ratio [OR] = 1.37, p = 0.0003) and p.Q65P (OR = 0.72, p = 0.005), were associated with PD, and p.M393T was also associated with RBD (OR = 1.59, p = 0.001). TMEM175 p.M393T was associated with reduced GCase activity. Homology modeling and normal mode analysis demonstrated that TMEM175 p.M393T creates a polar side-chain in the hydrophobic core of the transmembrane, which could destabilize the domain and thus impair either its assembly, maturation, or trafficking. Molecular dynamics simulations demonstrated that the p.Q65P variant may increase stability and ion conductance of the transmembrane protein, and lysosomal localization was not affected by these variants. INTERPRETATION: Coding variants in TMEM175 are likely to be responsible for the association in the TMEM175/GAK/DGKQ locus, which could be mediated by affecting GCase activity. ANN NEUROL 2020;87:139-153.

Cytokines and Gaucher Biomarkers in Glucocerebrosidase Carriers with and Without Parkinson Disease.

BACKGROUND: Homozygous and compound heterozygous variants in glucocerebrosidase (GBA) can cause Gaucher disease (GD), whereas heterozygous variants increase the risk of developing Parkinson's disease (PD). GD patients display altered peripheral immune proteins. However, it is unknown if these are altered in GBA carriers with PD. OBJECTIVES: To determine whether plasma cytokines and immune biomarkers associated with GD are also altered in GBA carriers with or without PD. METHODS: Inflammatory cytokines and established GD biomarkers, ferritin, CD162, CCL18, and chitotriosidase (28 biomarkers) were measured in GBA pathogenic variant carriers with (n = 135) and without (n = 83) PD, and non-carriers with (n = 75) and without PD (n = 77). RESULTS: PD patients with biallelic pathogenic variants in GBA had elevated plasma levels of ferritin, CCL18, and MIP1α. These biomarkers were not elevated in heterozygous GBA carriers. CONCLUSION: GD plasma biomarkers are not promising candidates for stratifying the risk for PD in carriers of heterozygous GBA pathogenic variants. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Lack of Causal Effects or Genetic Correlation between Restless Legs Syndrome and Parkinson's Disease.

BACKGROUND: Epidemiological studies have reported an association between Parkinson's disease (PD) and restless legs syndrome. OBJECTIVES: We aimed to use genetic data to study whether these 2 disorders are causally linked or share genetic architecture. METHODS: We performed two-sample Mendelian randomization and linkage disequilibrium score regression using summary statistics from recent genome-wide meta-analyses of PD and restless legs syndrome. RESULTS: We found no evidence for a causal relationship between restless legs syndrome (as the exposure) and PD (as the outcome, inverse variance-weighted; b = -0.003, SE = 0.031, P = 0.916; F statistic = 217.5). Reverse Mendelian randomization also did not demonstrate any causal effect of PD on restless legs syndrome (inverse variance-weighted; b = -0.012, SE = 0.023, P = 0.592; F statistic = 191.7). Linkage disequilibrium score regression analysis demonstrated lack of genetic correlation between restless legs syndrome and PD (rg = -0.028, SE = 0.042, P = 0.507). CONCLUSIONS: There was no evidence for a causal relationship or genetic correlation between restless legs syndrome and PD. The associations observed in epidemiological studies could be attributed, in part, to confounding or nongenetic determinants. © 2021 International Parkinson and Movement Disorder Society.

Analysis of Heterozygous PRKN Variants and Copy-Number Variations in Parkinson's Disease.

BACKGROUND: Biallelic PRKN mutation carriers with Parkinson's disease (PD) typically have an earlier disease onset, slow disease progression, and, often, different neuropathology compared to sporadic PD patients. However, the role of heterozygous PRKN variants in the risk of PD is controversial. OBJECTIVES: Our aim was to examine the association between heterozygous PRKN variants, including single-nucleotide variants and copy-number variations (CNVs), and PD. METHODS: We fully sequenced PRKN in 2809 PD patients and 3629 healthy controls, including 1965 late-onset (63.97 ± 7.79 years, 63% men) and 553 early-onset PD patients (43.33 ± 6.59 years, 68% men). PRKN was sequenced using targeted next-generation sequencing with molecular inversion probes. CNVs were identified using a combination of multiplex ligation-dependent probe amplification and ExomeDepth. To examine whether rare heterozygous single-nucleotide variants and CNVs in PRKN are associated with PD risk and onset, we used optimized sequence kernel association tests and regression models. RESULTS: We did not find any associations between all types of PRKN variants and risk of PD. Pathogenic and likely-pathogenic heterozygous single-nucleotide variants and CNVs were less common among PD patients (1.52%) than among controls (1.8%, false discovery rate-corrected P = 0.55). No associations with age at onset and in stratified analyses were found. CONCLUSIONS: Heterozygous single-nucleotide variants and CNVs in PRKN are not associated with PD. Molecular inversion probes allow for rapid and cost-effective detection of all types of PRKN variants, which may be useful for pretrial screening and for clinical and basic science studies targeting specifically PRKN patients. © 2020 International Parkinson and Movement Disorder Society.

Comprehensive Analysis of Familial Parkinsonism Genes in Rapid-Eye-Movement Sleep Behavior Disorder.

BACKGROUND: There is only partial overlap in the genetic background of isolated rapid-eye-movement sleep behavior disorder (iRBD) and Parkinson's disease (PD). OBJECTIVE: To examine the role of autosomal dominant and recessive PD or atypical parkinsonism genes in the risk of iRBD. METHODS: Ten genes, comprising the recessive genes PRKN, DJ-1 (PARK7), PINK1, VPS13C, ATP13A2, FBXO7, and PLA2G6 and the dominant genes LRRK2, GCH1, and VPS35, were fully sequenced in 1039 iRBD patients and 1852 controls of European ancestry, followed by association tests. RESULTS: We found no association between rare heterozygous variants in the tested genes and risk of iRBD. Several homozygous and compound heterozygous carriers were identified, yet there was no overrepresentation in iRBD patients versus controls. CONCLUSION: Our results do not support a major role for variants in these genes in the risk of iRBD. © 2020 International Parkinson and Movement Disorder Society.

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.

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.

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.

SPTAN1 variants as a potential cause for autosomal recessive hereditary spastic paraplegia.

More than 80 known or suspected genes/loci have been reported to be involved in hereditary spastic paraplegia (HSP). Genetic and clinical overlap have been reported between HSP and other neurological condition, yet about 50% of HSP patients remain genetically undiagnosed. To identify novel genes involved in HSP, we performed a genetic analysis of 383 HSP patients from 289 families with HSP. Two patients with biallelic SPTAN1 variants were identified; one carried the c.2572G>T p.(Ala858Ser) and c.4283C>G p.(Ala1428Gly) variants, and the second also carried the c.2572G>T p.(Ala858Ser) variant, and an additional variant, c.6990G>C p.(Met2330Ile). In silico predictive and structural analyses suggested that these variants are likely to be deleterious. SPTAN1 was highly intolerant for functional variants (in the top 0.31% of intolerant genes) with much lower observed vs. expected number of loss-of-function variants (8 vs. 142.7, p < 5 × 10-15). Using public databases of animal models and previously published data, we have found previously described zebrafish, mouse, and rat animal models of SPTAN1 deficiency, all consistently showing axonal degeneration, fitting the pathological features of HSP in humans. This study expands the phenotype of SPTAN1 mutations, which at the heterozygous state, when occurred de novo, may cause early infantile epileptic encephalopathy-5 (EIEE5). Our results further suggest that SPTAN1 may cause autosomal recessive HSP, and that it should be included in genetic screening panels for genetically undiagnosed HSP patients.

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.

Transcript-aware analysis of rare predicted loss-of-function variants in the UK Biobank elucidate new isoform-trait associations.

A single gene can produce multiple transcripts with distinct molecular functions. Rare-variant association tests often aggregate all coding variants across individual genes, without accounting for the variants' presence or consequence in resulting transcript isoforms. To evaluate the utility of transcript-aware variant sets, rare predicted loss-of-function (pLOF) variants were aggregated for 17,035 protein-coding genes using 55,558 distinct transcript-specific variant sets. These sets were tested for their association with 728 circulating proteins and 188 quantitative phenotypes across 406,921 individuals in the UK Biobank. The transcript-specific approach resulted in larger estimated effects of pLOF variants decreasing serum cis-protein levels compared to the gene-based approach (pbinom ≤ 2x10-16). Additionally, 251 quantitative trait associations were identified as being significant using the transcript-specific approach but not the gene-based approach, including PCSK5 transcript ENST00000376752 and standing height (transcript-specific statistic, P = 1.3x10-16, effect = 0.7 SD decrease; gene-based statistic, P = 0.02, effect = 0.05 SD decrease) and LDLR transcript ENST00000252444 and apolipoprotein B (transcript-specific statistic, P = 5.7x10-20, effect = 1.0 SD increase; gene-based statistic, P = 3.0x10-4, effect = 0.2 SD increase). This approach demonstrates the importance of considering the effect of pLOFs on specific transcript isoforms when performing rare-variant association studies.

Genome-wide association study of glucocerebrosidase activity modifiers.

One of the most common genetic risk factors for Parkinson's disease (PD) are variants in GBA1, which encodes the lysosomal enzyme glucocerebrosidase (GCase). GCase deficiency has been associated with an increased PD risk, but not all individuals with low GCase activity are carriers of GBA1 mutations, suggesting other factors may be acting as modifiers. We aimed to discover common variants associated with GCase activity, as well as replicate previously reported associations, by performing a genome-wide association study using two independent cohorts: a Columbia University cohort consisting of 697 PD cases and 347 controls and the Parkinson's Progression Markers Initiative (PPMI) cohort consisting of 357 PD cases and 163 controls. As expected, GBA1 variants have the strongest association with decreased activity, led by p.N370S (beta = -4.36, se = 0.32, p = 5.05e-43). We also identify a novel association in the GAA locus (encoding for acid alpha-glucosidase, beta = -0.96, se = 0.17, p = 5.23e-09) that may be the result of an interaction between GCase and acid alpha-glucosidase based on various interaction analyses. Lastly, we show that several PD-risk loci are potentially associated with GCase activity. Further research will be needed to replicate and validate our findings and to uncover the functional connection between acid alpha-glucosidase and GCase.

Genome-wide association study stratified by MAPT haplotypes identifies potential novel loci in Parkinson's disease.

Objective: To identify genetic factors that may modify the effects of the MAPT locus in Parkinson's disease (PD). Methods: We used data from the International Parkinson's Disease Genomics Consortium (IPDGC) and the UK biobank (UKBB). We stratified the IPDGC cohort for carriers of the H1/H1 genotype (PD patients n=8,492 and controls n=6,765) and carriers of the H2 haplotype (with either H1/H2 or H2/H2 genotypes, patients n=4,779 and controls n=4,849) to perform genome-wide association studies (GWASs). Then, we performed replication analyses in the UKBB data. To study the association of rare variants in the new nominated genes, we performed burden analyses in two cohorts (Accelerating Medicines Partnership - Parkinson Disease and UKBB) with a total sample size PD patients n=2,943 and controls n=18,486. Results: We identified a novel locus associated with PD among MAPT H1/H1 carriers near EMP1 (rs56312722, OR=0.88, 95%CI= 0.84-0.92, p= 1.80E-08), and a novel locus associated with PD among MAPT H2 carriers near VANGL1 (rs11590278, OR=1.69 95%CI=1.40-2.03, p=2.72E-08). Similar analysis of the UKBB data did not replicate these results and rs11590278 near VANGL1 did have similar effect size and direction in carriers of H2 haplotype, albeit not statistically significant (OR= 1.32, 95%CI= 0.94-1.86, p=0.17). Rare EMP1 variants with high CADD scores were associated with PD in the MAPT H2 stratified analysis (p=9.46E-05), mainly driven by the p.V11G variant. Interpretation: We identified several loci potentially associated with PD stratified by MAPT haplotype and larger replication studies are required to confirm these associations.