Rare-variant burden across lysosomal genes implicates sialylation and ganglioside metabolism in Parkinson's disease

Lysosomal dysfunction is central to Parkinsons disease pathogenesis, with GBA1 as the strongest established genetic risk factor. Numerous other genes involved in lysosomal sphingolipid, glycosphingolipid and ceramide metabolism have been proposed as contributors to Parkinsons disease, underscoring the need for comprehensive genetic analyses across these pathways. We analysed rare variants (minor allele frequency < 0.01) across 36 lysosomal genes (excluding GBA1) in 8,267 individuals with Parkinsons disease and 68,208 controls, including a subset of 793 early-onset Parkinsons disease ([&le;]50 years) cases. Targeted sequencing was performed in four cohorts at McGill University (3,456 Parkinsons disease patients and 2,664 controls) and results were combined with whole-genome sequencing data from the UK Biobank (2,848 cases, 62,451 controls), and from the Accelerating Medicines Partnership - Parkinsons Disease (1,963 cases, 3,093 controls). We analysed the association of rare variants in these genes with Parkinsons disease using Sequence Kernel Association Test-Optimal (SKAT-O) across variant classes (all rare variants, nonsynonymous, loss-of-function and predicted damaging variants with a Combined Annotation Dependent Depletion (CADD) score >20), with meta-analysis across cohorts. We additionally performed per-domain analyses for variants in gene segments encoding functional domains. False discovery rate correction was applied. Meta-analysis identified a significant association between rare variants in ST3GAL3 and Parkinsons disease (Pfdr=0.04). Several additional lysosomal genes showed nominal associations (P<0.05), including HGSNAT, ASAH1, CTSD, HEXA, ST3GAL4 and SGPP1. Domain-based analyses identified a strong enrichment of nonsynonymous variants within the beta-acetyl-hexosaminidase-like domain of HEXA (P = 8.0 x 10), although this signal did not survive correction for multiple testing (Pfdr=0.154). In early-onset Parkinsons disease, domain-based analyses revealed significant associations in NAGLU (Pfdr=7.3x10) and ST3GAL5 (Pfdr=0.03). Together, these results provide genetic evidence that rare variants across multiple lysosomal pathways, particularly those related to sialylation, ganglioside metabolism, ceramide biology, and lysosomal proteolysis, may contribute to Parkinsons disease susceptibility beyond GBA1, highlighting biologically coherent pathways for future replication and functional investigation.