Characterizing a complex CT-rich haplotype in intron 4 of SNCA using large-scale targeted amplicon long-read sequencing.

Parkinson's disease (PD) is a common neurodegenerative disorder with a significant risk proportion driven by genetics. While much progress has been made, most of the heritability remains unknown. This is in-part because previous genetic studies have focused on the contribution of single nucleotide variants. More complex forms of variation, such as structural variants and tandem repeats, are already associated with several synucleinopathies. However, because more sophisticated sequencing methods are usually required to detect these regions, little is understood regarding their contribution to PD. One example is a polymorphic CT-rich region in intron 4 of the SNCA gene. This haplotype has been suggested to be associated with risk of Lewy Body (LB) pathology in Alzheimer's Disease and SNCA gene expression, but is yet to be investigated in PD. Here, we attempt to resolve this CT-rich haplotype and investigate its role in PD. We performed targeted PacBio HiFi sequencing of the region in 1375 PD cases and 959 controls. We replicate the previously reported associations and a novel association between two PD risk SNVs (rs356182 and rs5019538) and haplotype 4, the largest haplotype. Through quantitative trait locus analyzes we identify a significant haplotype 4 association with alternative CAGE transcriptional start site usage, not leading to significant differential SNCA gene expression in post-mortem frontal cortex brain tissue. Therefore, disease association in this locus might not be biologically driven by this CT-rich repeat region. Our data demonstrates the complexity of this SNCA region and highlights that further follow up functional studies are warranted.

An ANXA11 P93S variant dysregulates TDP-43 and causes corticobasal syndrome.

INTRODUCTION: Variants of uncertain significance (VUS) surged with affordable genetic testing, posing challenges for determining pathogenicity. We examine the pathogenicity of a novel VUS P93S in Annexin A11 (ANXA11) - an amyotrophic lateral sclerosis/frontotemporal dementia-associated gene - in a corticobasal syndrome kindred. Established ANXA11 mutations cause ANXA11 aggregation, altered lysosomal-RNA granule co-trafficking, and transactive response DNA binding protein of 43 kDa (TDP-43) mis-localization. METHODS: We described the clinical presentation and explored the phenotypic diversity of ANXA11 variants. P93S's effect on ANXA11 function and TDP-43 biology was characterized in induced pluripotent stem cell-derived neurons alongside multiomic neuronal and microglial profiling. RESULTS: ANXA11 mutations were linked to corticobasal syndrome cases. P93S led to decreased lysosome colocalization, neuritic RNA, and nuclear TDP-43 with cryptic exon expression. Multiomic microglial signatures implicated immune dysregulation and interferon signaling pathways. DISCUSSION: This study establishes ANXA11 P93S pathogenicity, broadens the phenotypic spectrum of ANXA11 mutations, underscores neuronal and microglial dysfunction in ANXA11 pathophysiology, and demonstrates the potential of cellular models to determine variant pathogenicity. HIGHLIGHTS: ANXA11 P93S is a pathogenic variant. Corticobasal syndrome is part of the ANXA11 phenotypic spectrum. Hybridization chain reaction fluorescence in situ hybridization (HCR FISH) is a new tool for the detection of cryptic exons due to TDP-43-related loss of splicing regulation. Microglial ANXA11 and related immune pathways are important drivers of disease. Cellular models are powerful tools for adjudicating variants of uncertain significance.

A SINE-VNTR-Alu at the LRIG2 locus is associated with proximal and distal gene expression in CRISPR and population models.

SINE-VNTR-Alu (SVA) retrotransposons represent mobile regulatory elements that have the potential to influence the surrounding genome when they insert into a locus. Evolutionarily recent mobilisation has resulted in loci in the human genome where a given retrotransposon might be observed to be present or absent, termed a retrotransposon insertion polymorphism (RIP). We previously observed that an SVA RIP ~ 2 kb upstream of LRIG2 on chromosome 1, the 'LRIG2 SVA', was associated with differences in local gene expression and methylation, and that the two were correlated. Here, we have used CRISPR-mediated deletion of the LRIG2 SVA in a cell line model to validate that presence of the retrotransposon is directly affecting local expression and provide evidence that is suggestive of a modest role for the SVA in modulating nearby methylation. Additionally, in leveraging an available Hi-C dataset we observed that the LRIG2 SVA was also involved in long-range chromatin interactions with a cluster of genes ~ 300 kb away, and that expression of these genes was to varying degrees associated with dosage of the SVA in both CRISPR cell line and population models. Altogether, these data support a regulatory role for SVAs in the modulation of gene expression, with the latter potentially involving chromatin looping, consistent with the model that RIPs may contribute to interpersonal differences in transcriptional networks.

Human brain single nucleus cell type enrichments in neurodegenerative diseases.

Background: Single-cell RNA sequencing has opened a window into clarifying the complex underpinnings of disease, particularly in quantifying the relevance of tissue- and cell-type-specific gene expression. Methods: To identify the cell types and genes important to therapeutic target development across the neurodegenerative disease spectrum, we leveraged genome-wide association studies, recent single-cell sequencing data, and bulk expression studies in a diverse series of brain region tissues. Results: We were able to identify significant immune-related cell types in the brain across three major neurodegenerative diseases: Alzheimer's disease, amyotrophic lateral sclerosis, and Parkinson's disease. Subsequently, putative roles of 30 fine-mapped loci implicating seven genes in multiple neurodegenerative diseases and their pathogenesis were identified. Conclusions: We have helped refine the genetic regions and cell types effected across multiple neurodegenerative diseases, helping focus future translational research efforts.

An ANXA11 P93S variant dysregulates TDP-43 and causes corticobasal syndrome.

As genetic testing has become more accessible and affordable, variants of uncertain significance (VUS) are increasingly identified, and determining whether these variants play causal roles in disease is a major challenge. The known disease-associated Annexin A11 (ANXA11) mutations result in ANXA11 aggregation, alterations in lysosomal-RNA granule co-trafficking, and TDP-43 mis-localization and present as amyotrophic lateral sclerosis or frontotemporal dementia. We identified a novel VUS in ANXA11 (P93S) in a kindred with corticobasal syndrome and unique radiographic features that segregated with disease. We then queried neurodegenerative disorder clinic databases to identify the phenotypic spread of ANXA11 mutations. Multi-modal computational analysis of this variant was performed and the effect of this VUS on ANXA11 function and TDP-43 biology was characterized in iPSC-derived neurons. Single-cell sequencing and proteomic analysis of iPSC-derived neurons and microglia were used to determine the multiomic signature of this VUS. Mutations in ANXA11 were found in association with clinically diagnosed corticobasal syndrome, thereby establishing corticobasal syndrome as part of ANXA11 clinical spectrum. In iPSC-derived neurons expressing mutant ANXA11, we found decreased colocalization of lysosomes and decreased neuritic RNA as well as decreased nuclear TDP-43 and increased formation of cryptic exons compared to controls. Multiomic assessment of the P93S variant in iPSC-derived neurons and microglia indicates that the pathogenic omic signature in neurons is modest compared to microglia. Additionally, omic studies reveal that immune dysregulation and interferon signaling pathways in microglia are central to disease. Collectively, these findings identify a new pathogenic variant in ANXA11, expand the range of clinical syndromes caused by ANXA11 mutations, and implicate both neuronal and microglia dysfunction in ANXA11 pathophysiology. This work illustrates the potential for iPSC-derived cellular models to revolutionize the variant annotation process and provides a generalizable approach to determining causality of novel variants across genes.

Scalable Nanopore sequencing of human genomes provides a comprehensive view of haplotype-resolved variation and methylation.

Long-read sequencing technologies substantially overcome the limitations of short-reads but have not been considered as a feasible replacement for population-scale projects, being a combination of too expensive, not scalable enough or too error-prone. Here we develop an efficient and scalable wet lab and computational protocol, Napu, for Oxford Nanopore Technologies long-read sequencing that seeks to address those limitations. We applied our protocol to cell lines and brain tissue samples as part of a pilot project for the National Institutes of Health Center for Alzheimer's and Related Dementias. Using a single PromethION flow cell, we can detect single nucleotide polymorphisms with F1-score comparable to Illumina short-read sequencing. Small indel calling remains difficult within homopolymers and tandem repeats, but achieves good concordance to Illumina indel calls elsewhere. Further, we can discover structural variants with F1-score on par with state-of-the-art de novo assembly methods. Our protocol phases small and structural variants at megabase scales and produces highly accurate, haplotype-specific methylation calls.

Human brain single nucleus cell type enrichments in neurodegenerative diseases.

Single cell RNA sequencing has opened a window into clarifying the complex underpinnings of disease, particularly in quantifying the relevance of tissue- and cell-type-specific gene expression. To identify the cell types and genes important to therapeutic target development across the neurodegenerative disease spectrum, we leveraged genome-wide association studies, recent single cell sequencing data, and bulk expression studies in a diverse series of brain region tissues. We were able to identify significant immune-related cell types in the brain across three major neurodegenerative diseases: Alzheimer's Disease, Amyotrophic Lateral Sclerosis, and Parkinson's Diseases. Subsequently, we identified the major role of 30 fine-mapped loci implicating seven genes in multiple neurodegenerative diseases and their pathogenesis.

Inhibition of p38α MAPK restores neuronal p38γ MAPK and ameliorates synaptic degeneration in a mouse model of DLB/PD.

Alterations in the p38 mitogen-activated protein kinases (MAPKs) play an important role in the pathogenesis of dementia with Lewy bodies (DLB) and Parkinson's disease (PD). Activation of the p38α MAPK isoform and mislocalization of the p38γ MAPK isoform are associated with neuroinflammation and synaptic degeneration in DLB and PD. Therefore, we hypothesized that p38α might be associated with neuronal p38γ distribution and synaptic dysfunction in these diseases. To test this hypothesis, we treated in vitro cellular and in vivo mouse models of DLB/PD with SKF-86002, a compound that attenuates inflammation by inhibiting p38α/ß, and then investigated the effects of this compound on p38γ and neurodegenerative pathology. We found that inhibition of p38α reduced neuroinflammation and ameliorated synaptic, neurodegenerative, and motor behavioral deficits in transgenic mice overexpressing human α-synuclein. Moreover, treatment with SKF-86002 promoted the redistribution of p38γ to synapses and reduced the accumulation of α-synuclein in mice overexpressing human α-synuclein. Supporting the potential value of targeting p38 in DLB/PD, we found that SKF-86002 promoted the redistribution of p38γ in neurons differentiated from iPS cells derived from patients with familial PD (carrying the A53T α-synuclein mutation) and healthy controls. Treatment with SKF-86002 ameliorated α-synuclein-induced neurodegeneration in these neurons only when microglia were pretreated with this compound. However, direct treatment of neurons with SKF-86002 did not affect α-synuclein-induced neurotoxicity, suggesting that SKF-86002 treatment inhibits α-synuclein-induced neurotoxicity mediated by microglia. These findings provide a mechanistic connection between p38α and p38γ as well as a rationale for targeting this pathway in DLB/PD.

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.

LRRK2: Genetic mechanisms vs genetic subtypes.

In 2004, the identification of pathogenic variants in the LRRK2 gene across several families with autosomal dominant late-onset Parkinson's disease (PD) revolutionized our understanding of the role of genetics in PD. Previous beliefs that genetics in PD was limited to rare early-onset or familial forms of the disease were quickly dispelled. Currently, we recognize LRRK2 p.G2019S as the most common genetic cause of both sporadic and familial PD, with more than 100,000 affected carriers across the globe. The frequency of LRRK2 p.G2019S is also highly variable across populations, with some regions of Asian or Latin America reporting close to 0%, contrasting to Ashkenazi Jews or North African Berbers reporting up to 13% and 40%, respectively. Patients with LRRK2 pathogenic variants are clinically and pathologically heterogeneous, highlighting the age-related variable penetrance that also characterizes LRRK2-related disease. Indeed, the majority of patients with LRRK2-related disease are characterized by a relatively mild Parkinsonism with less motor symptoms with variable presence of α-synuclein and/or tau aggregates, with pathologic pleomorphism widely described. At a functional cellular level, it is likely that pathogenic variants mediate a toxic gain-of-function of the LRRK2 protein resulting in increased kinase activity perhaps in a cell-specific manner; by contrast, some LRRK2 variants appear to be protective reducing PD risk by decreasing the kinase activity. Therefore, employing this information to define appropriate patient populations for clinical trials of targeted kinase LRRK2 inhibition strategies is very promising and demonstrates a potential future application for PD using precision medicine.

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.

Scalable Nanopore sequencing of human genomes provides a comprehensive view of haplotype-resolved variation and methylation.

Long-read sequencing technologies substantially overcome the limitations of short-reads but to date have not been considered as feasible replacement at scale due to a combination of being too expensive, not scalable enough, or too error-prone. Here, we develop an efficient and scalable wet lab and computational protocol for Oxford Nanopore Technologies (ONT) long-read sequencing that seeks to provide a genuine alternative to short-reads for large-scale genomics projects. We applied our protocol to cell lines and brain tissue samples as part of a pilot project for the NIH Center for Alzheimer's and Related Dementias (CARD). Using a single PromethION flow cell, we can detect SNPs with F1-score better than Illumina short-read sequencing. Small indel calling remains to be difficult inside homopolymers and tandem repeats, but is comparable to Illumina calls elsewhere. Further, we can discover structural variants with F1-score comparable to state-of the-art methods involving Pacific Biosciences HiFi sequencing and trio information (but at a lower cost and greater throughput). Using ONT based phasing, we can then combine and phase small and structural variants at megabase scales. Our protocol also produces highly accurate, haplotype-specific methylation calls. Overall, this makes large-scale long-read sequencing projects feasible; the protocol is currently being used to sequence thousands of brain-based genomes as a part of the NIH CARD initiative. We provide the protocol and software as open-source integrated pipelines for generating phased variant calls and assemblies.

Association of a common genetic variant with Parkinson's disease is mediated by microglia.

Studies of multiple neurodegenerative disorders have identified many genetic variants that are associated with risk of disease throughout a lifetime. For example, Parkinson's disease (PD) risk is attributed in part to both coding mutations in the leucine-rich repeat kinase 2 (LRRK2) gene and to a common noncoding variation in the 5' region of the LRRK2 locus, as identified by genome-wide association studies (GWAS). However, the mechanisms linking GWAS variants to pathogenicity are largely unknown. Here, we found that the influence of PD-associated noncoding variation on LRRK2 expression is specifically propagated through microglia and not by other cell types that express LRRK2 in the human brain. We find microglia-specific regulatory chromatin regions that modulate the LRRK2 expression in human frontal cortex and substantia nigra and confirm these results in a human-induced pluripotent stem cell-derived microglia model. We showed, using a large-scale clustered regularly interspaced short palindromic repeats interference (CRISPRi) screen, that a regulatory DNA element containing the single-nucleotide variant rs6581593 influences the LRRK2 expression in microglia. Our study demonstrates that cell type should be considered when evaluating the role of noncoding variation in disease pathogenesis and sheds light on the mechanism underlying the association of the 5' region of LRRK2 with PD risk.

Coding and Noncoding Variation in LRRK2 and Parkinson's Disease Risk.

BACKGROUND: The leucine-rich repeat kinase 2 (LRRK2) gene harbors both rare highly damaging missense variants (eg, p.G2019S) and common noncoding variants (eg, rs76904798) with lower effect sizes that are associated with Parkinson's disease (PD) risk. OBJECTIVES: This study aimed to investigate in a large meta-analysis whether the LRRK2 Genome-Wide Association Study (GWAS) signal represented by rs76904798 is independently associated with PD risk from LRRK2 coding variation and whether complex linkage disequilibrium structures with p.G2019S and the 5' noncoding haplotype account for the association of LRRK2 coding variants. METHODS: We performed a meta-analysis using imputed genotypes from 17,838 patients, 13,404 proxy patients, and 173,639 healthy controls of European ancestry. We excluded carriers of p.G2019S and/or rs76904798 to clarify the role of LRRK2 coding variation in mediating disease risk and excluded carriers of relatively rare LRRK2 coding variants to assess the independence of rs76904798. We also investigated the co-inheritance of LRRK2 coding variants with p.G2019S, rs76904798, and p.N2081D. RESULTS: LRRK2 rs76904798 remained significantly associated with PD after excluding the carriers of relatively rare LRRK2 coding variants. LRRK2 p.R1514Q and p.N2081D were frequently co-inherited with rs76904798, and the allele distribution of p.S1647T significantly changed among patients after removing rs76904798 carriers. CONCLUSIONS: These data suggest that the LRRK2 coding variants previously related to PD (p.N551K, p.R1398H, p.M1646T, and p.N2081D) do not drive the 5' noncoding GWAS signal. These data, however, do not preclude the independent association of the haplotype p.N551K-p.R1398H and p.M1646T with altered disease risk. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson Movement Disorder Society. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Transcriptional signatures in iPSC-derived neurons are reproducible across labs when differentiation protocols are closely matched.

Reproducibility of expression patterns in iPSC-derived cells from different labs is an important first step in ensuring replication of biochemical or functional assays that are performed in different labs. Here we show that reproducible gene expression patterns from iPSCs and iPSC-derived neurons matured and collected at two separate laboratory locations can be achieved by closely matching protocols and reagents. While there are significant differences in gene expression between iPSCs and differentiated neurons, as well as between different donor lines of the same cell type, transcriptional changes that vary with laboratory sites are relatively small. These results suggest that making great efforts to match protocols, reagents and technical methods between labs may improve the reproducibility of iPSC-derived cell models.

THAP1 modulates oligodendrocyte maturation by regulating ECM degradation in lysosomes.

Mechanisms controlling myelination during central nervous system (CNS) maturation play a pivotal role in the development and refinement of CNS circuits. The transcription factor THAP1 is essential for timing the inception of myelination during CNS maturation through a cell-autonomous role in the oligodendrocyte lineage. Here, we demonstrate that THAP1 modulates the extracellular matrix (ECM) composition by regulating glycosaminoglycan (GAG) catabolism within oligodendrocyte progenitor cells (OPCs). Thap1-/- OPCs accumulate and secrete excess GAGs, inhibiting their maturation through an autoinhibitory mechanism. THAP1 controls GAG metabolism by binding to and regulating the GusB gene encoding ß-glucuronidase, a GAG-catabolic lysosomal enzyme. Applying GAG-degrading enzymes or overexpressing ß-glucuronidase rescues Thap1-/- OL maturation deficits in vitro and in vivo. Our studies establish lysosomal GAG catabolism within OPCs as a critical mechanism regulating oligodendrocyte development.

Generation of fourteen isogenic cell lines for Parkinson's disease-associated leucine-rich repeat kinase (LRRK2).

Mutations in leucine-rich repeat kinase 2 (LRRK2) are associated with inherited forms of Parkinson's disease (PD), causing disease by a gain of kinase function. Here, we describe a series of isogenic iPSC lines with any of five pathogenic mutations (N1437H, R1441C, Y1699C, G2019S and I2020T); two hypothesis testing mutations (GTP binding null, T1348N, and kinase dead, K1906M) and two LRRK2 knockouts. This resource could be used to assess effects of mutations on the function of endogenous LRRK2 and/or to study LRRK2 interactors and substrates in iPSC-derived cellular models.

Advancing Personalized Medicine in Common Forms of Parkinson's Disease through Genetics: Current Therapeutics and the Future of Individualized Management.

Parkinson's disease (PD) is a condition with heterogeneous clinical manifestations that vary in age at onset, rate of progression, disease course, severity, motor and non-motor symptoms, and a variable response to antiparkinsonian drugs. It is considered that there are multiple PD etiological subtypes, some of which could be predicted by genetics. The characterization and prediction of these distinct molecular entities provides a growing opportunity to use individualized management and personalized therapies. Dissecting the genetic architecture of PD is a critical step in identifying therapeutic targets, and genetics represents a step forward to sub-categorize and predict PD risk and progression. A better understanding and separation of genetic subtypes has immediate implications in clinical trial design by unraveling the different flavors of clinical presentation and development. Personalized medicine is a nascent area of research and represents a paramount challenge in the treatment and cure of PD. This manuscript summarizes the current state of precision medicine in the PD field and discusses how genetics has become the engine to gain insights into disease during our constant effort to develop potential etiological based interventions.

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.