Large-scale rare variant burden testing in Parkinson's disease

Parkinson's disease has a large heritable component and genome-wide association studies have identified over 90 variants with disease-associated common variants, providing deeper insights into the disease biology. However, there have not been large-scale rare variant analyses for Parkinson's disease. To address this gap, we investigated the rare genetic component of Parkinson's disease at minor allele frequencies <1%, using whole genome and whole exome sequencing data from 7184 Parkinson's disease cases, 6701 proxy cases and 51 650 healthy controls from the Accelerating Medicines Partnership Parkinson's disease (AMP-PD) initiative, the National Institutes of Health, the UK Biobank and Genentech. We performed burden tests meta-analyses on small indels and single nucleotide protein-altering variants, prioritized based on their predicted functional impact. Our work identified several genes reaching exome-wide significance. Two of these genes, GBA1 and LRRK2, have variants that have been previously implicated as risk factors for Parkinson's disease, with some variants in LRRK2 resulting in monogenic forms of the disease. We identify potential novel risk associations for variants in B3GNT3, AUNIP, ADH5, TUBA1B, OR1G1, CAPN10 and TREML1 but were unable to replicate the observed associations across independent datasets. Of these, B3GNT3 and TREML1 could provide new evidence for the role of neuroinflammation in Parkinson's disease. To date, this is the largest analysis of rare genetic variants in Parkinson's disease.

Genome-Wide Analysis of Structural Variants in Parkinson Disease

OBJECTIVE

Identification of genetic risk factors for Parkinson disease (PD) has to date been primarily limited to the study of single nucleotide variants, which only represent a small fraction of the genetic variation in the human genome. Consequently, causal variants for most PD risk are not known. Here we focused on structural variants (SVs), which represent a major source of genetic variation in the human genome. We aimed to discover SVs associated with PD risk by performing the first large-scale characterization of SVs in PD.

METHODS

We leveraged a recently developed computational pipeline to detect and genotype SVs from 7,772 Illumina short-read whole genome sequencing samples. Using this set of SV variants, we performed a genome-wide association study using 2,585 cases and 2,779 controls and identified SVs associated with PD risk. Furthermore, to validate the presence of these variants, we generated a subset of matched whole-genome long-read sequencing data.

RESULTS

We genotyped and tested 3,154 common SVs, representing over 412 million nucleotides of previously uncatalogued genetic variation. Using long-read sequencing data, we validated the presence of three novel deletion SVs that are associated with risk of PD from our initial association analysis, including a 2 kb intronic deletion within the gene LRRN4.

INTERPRETATION

We identified three SVs associated with genetic risk of PD. This study represents the most comprehensive assessment of the contribution of SVs to the genetic risk of PD to date. ANN NEUROL 2023;93:1012-1022.

The Foundational Data Initiative for Parkinson Disease: Enabling efficient translation from genetic maps to mechanism

The Foundational Data Initiative for Parkinson Disease (FOUNDIN-PD) is an international collaboration producing fundamental resources for Parkinson disease (PD). FOUNDIN-PD generated a multi-layered molecular dataset in a cohort of induced pluripotent stem cell (iPSC) lines differentiated to dopaminergic (DA) neurons, a major affected cell type in PD. The lines were derived from the Parkinson's Progression Markers Initiative study, which included participants with PD carrying monogenic PD variants, variants with intermediate effects, and variants identified by genome-wide association studies and unaffected individuals. We generated genetic, epigenetic, regulatory, transcriptomic, and longitudinal cellular imaging data from iPSC-derived DA neurons to understand molecular relationships between disease-associated genetic variation and proximate molecular events. These data reveal that iPSC-derived DA neurons provide a valuable cellular context and foundational atlas for modeling PD genetic risk. We have integrated these data into a FOUNDIN-PD data browser as a resource for understanding the molecular pathogenesis of PD.

The IPDGC/GP2 Hackathon - an open science event for training in data science, genomics, and collaboration using Parkinson's disease data

Open science and collaboration are necessary to facilitate the advancement of Parkinson's disease (PD) research. Hackathons are collaborative events that bring together people with different skill sets and backgrounds to generate resources and creative solutions to problems. These events can be used as training and networking opportunities, thus we coordinated a virtual 3-day hackathon event, during which 49 early-career scientists from 12 countries built tools and pipelines with a focus on PD. Resources were created with the goal of helping scientists accelerate their own research by having access to the necessary code and tools. Each team was allocated one of nine different projects, each with a different goal. These included developing post-genome-wide association studies (GWAS) analysis pipelines, downstream analysis of genetic variation pipelines, and various visualization tools. Hackathons are a valuable approach to inspire creative thinking, supplement training in data science, and foster collaborative scientific relationships, which are foundational practices for early-career researchers. The resources generated can be used to accelerate research on the genetics of PD.

Classification of GBA1 Variants in Parkinson's Disease: The GBA1-PD Browser

BACKGROUND

GBA1 variants are among the most common genetic risk factors for Parkinson's disease (PD). GBA1 variants can be classified into three categories based on their role in Gaucher's disease (GD) or PD: severe, mild, and risk variant (for PD).

OBJECTIVE

This review aims to generate and share a comprehensive database for GBA1 variants reported in PD to support future research and clinical trials.

METHODS

We performed a literature search for all GBA1 variants that have been reported in PD. The data have been standardized and complemented with variant classification, odds ratio if available, and other data.

RESULTS

We found 371 GBA1 variants reported in PD: 22 mild, 84 severe, 3 risk variants, and 262 of unknown status. We created a browser containing up-to-date information on these variants (https://pdgenetics.shinyapps.io/GBA1Browser/).

CONCLUSIONS

The classification and browser presented in this work should inform and support basic, translational, and clinical research on GBA1-PD. © 2023 International Parkinson and Movement Disorder Society.

GALC variants affect galactosylceramidase enzymatic activity and risk of Parkinson’s disease

The association between glucocerebrosidase (GCase), encoded by GBA, and Parkinson’s disease highlights the role of the lysosome in Parkinson’s disease pathogenesis. Genome-wide association studies (GWAS) in Parkinson’s disease have revealed multiple associated loci, including the GALC locus on chromosome 14. GALC encodes the lysosomal enzyme galactosylceramidase (GalCase), which plays a pivotal role in the glycosphingolipid metabolism pathway. It is still unclear whether GALC is the gene driving the association in the chromosome 14 locus, and if so, by which mechanism.

We first aimed to examine whether variants in the GALC locus and across the genome are associated with GalCase activity. We performed a GWAS in two independent cohorts from a)Columbia University and b)the Parkinson’s Progression Markers Initiative study, followed by a meta-analysis with a total of 976 Parkinson’s disease patients and 478 controls with available data on GalCase activity. We further analyzed the effects of common GALC variants on expression and GalCase activity using genomic colocalization methods. Mendelian randomization was used to study whether GalCase activity may be causal in Parkinson’s disease. To study the role of rare GALC variants we analyzed sequencing data from 5,028 Parkinson’s disease patients and 5,422 controls. Additionally, we studied the functional impact of GALC knock-out on alpha-synuclein accumulation and on GCase activity in neuronal cell models and performed in silico structural analysis of common GALC variants associated with altered GalCase activity.

The top hit in Parkinson's disease GWAS in the GALC locus, rs979812, is associated with increased GalCase activity (b = 1.2; se = 0.06; p = 5.10E-95). No other variants outside the GALC locus were associated with GalCase activity. Colocalization analysis demonstrated that rs979812 was also associated with increased GalCase expression. Mendelian randomization suggested that increased GalCase activity may be causally associated with Parkinson’s disease (b = 0.025, se = 0.007, p = 0.0008). We did not find an association between rare GALC variants and Parkinson’s disease. GALC knockout using CRISPR-Cas9 did not lead to alpha-synuclein accumulation, further supporting that increased rather than reduced GalCase levels may be associated with Parkinson’s disease. The structural analysis demonstrated that the common variant p.I562T may lead to improper maturation of GalCase affecting its activity.

Our results nominate GALC as the gene associated with Parkinson’s disease in this locus and suggest that the association of variants in the GALC locus may be driven by their effect of increasing GalCase expression and activity. Whether altering GalCase activity could be considered as a therapeutic target should be further studied.

Analysis of Y chromosome haplogroups in Parkinson's disease

Parkinson's disease is a complex neurodegenerative disorder that is about 1.5 times more prevalent in males than females. Extensive work has been done to identify the genetic risk factors behind Parkinson's disease on autosomes and more recently on Chromosome X, but work remains to be done on the male-specific Y chromosome. In an effort to explore the role of the Y chromosome in Parkinson's disease, we analysed whole-genome sequencing data from the Accelerating Medicines Partnership-Parkinson's disease initiative (1466 cases and 1664 controls), genotype data from NeuroX (3491 cases and 3232 controls) and genotype data from UKBiobank (182 517 controls, 1892 cases and 3783 proxy cases), all consisting of male European ancestry samples. We classified sample Y chromosomes by haplogroup using three different tools for comparison (Snappy, Yhaplo and Y-LineageTracker) and meta-analysed this data to identify haplogroups associated with Parkinson's disease. This was followed up with a Y-chromosome association study to identify specific variants associated with disease. We also analysed blood-based RNASeq data obtained from the Accelerating Medicines Partnership-Parkinson's disease initiative (1020 samples) and RNASeq data obtained from the North American Brain Expression Consortium (171 samples) to identify Y-chromosome genes differentially expressed in cases, controls, specific haplogroups and specific tissues. RNASeq analyses suggest Y-chromosome gene expression differs between brain and blood tissues but does not differ significantly in cases, controls or specific haplogroups. Overall, we did not find any strong associations between Y-chromosome genetics and Parkinson's disease, suggesting the explanation for the increased prevalence in males may lie elsewhere.

Assessment of ANG variants in Parkinson's disease

Genetic risk factors are occasionally shared between different neurodegenerative diseases. Previous studies have linked ANG, a gene encoding angiogenin, to both Parkinson's disease (PD) and amyotrophic lateral sclerosis (ALS). Functional studies suggest ANG plays a neuroprotective role in both PD and amyotrophic lateral sclerosis by reducing cell death. We further explored the genetic association between ANG and PD by analyzing genotype data from the International Parkinson's Disease Genomics Consortium (14,671 cases and 17,667 controls) and whole genome sequencing data from the Accelerating Medicines Partnership - Parkinson's disease initiative (AMP-PD, https://amp-pd.org/) (1,647 cases and 1,050 controls). Our analysis did not replicate the findings of previous studies and identified no significant association between ANG variants and PD risk.

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.

The Parkinson's Disease DNA Variant Browser

BACKGROUND

Parkinson's disease (PD) is a genetically complex neurodegenerative disease with ~20 genes known to contain mutations that cause PD or atypical parkinsonism. Large-scale next-generation sequencing projects have revolutionized genomics research. Applying these data to PD, many genes have been reported to contain putative disease-causing mutations. In most instances, however, the results remain quite limited and rather preliminary. Our aim was to assist researchers on their search for PD-risk genes and variant candidates with an easily accessible and open summary-level genomic data browser for the PD research community.

METHODS

Sequencing and imputed genotype data were obtained from multiple sources and harmonized and aggregated.

RESULTS

In total we included a total of 102,127 participants, including 28,453 PD cases, 1650 proxy cases, and 72,024 controls.

CONCLUSIONS

We present here the Parkinson's Disease Sequencing Browser: a Shiny-based web application that presents comprehensive summary-level frequency data from multiple large-scale genotyping and sequencing projects https://pdgenetics.shinyapps.io/VariantBrowser/. Published © 2021 This article is a U.S. Government work and is in the public domain in the USA. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

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.

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.

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.