The genetic drivers of juvenile, young, and early-onset Parkinson’s Disease in India

Shan V Andrews, Prashanth L Kukkle, Ramesh Menon, Thenral S Geetha, Vinay Goyal, Rukmini Mridula Kandadai, Hrishikesh Kumar, Rupam Borgohain, Adreesh Mukherjee, Pettarusp M Wadia, Ravi Yadav, Soaham Desai, Niraj Kumar, Deepika Joshi, Sakthivel Murugan, Atanu Biswas, Pramod K Pal, Merina Oliver, Sandhya Nair, Anbu Kayalvizhi, Praveena L Samson, Manjari Deshmukh, Akshi Bassi, Charugulla Sandeep, Nitin Mandloi, Oliver B Davis, Melissa A Roberts, Dara E Leto, Anastasia G Henry, Gilbert Di Paolo, Uday Muthane, Shymal K Das, Andrew S Peterson, Thomas Sandmann, Ravi Gupta, Vedam L Ramprasad, Parkinson Research Alliance of India (PRAI)
doi: https://doi.org/10.1101/2023.06.18.23291407

Abstract

Background Recent studies have advanced our understanding of the genetic drivers of Parkinson’s Disease (PD). Rare variants in more than 20 genes are considered causal for PD, and the latest PD GWAS study identified 90 independent risk loci. However, there remains a gap in our understanding of PD genetics outside of the European populations in which the vast majority of these studies were focused.

Objectives To identify genetic risk factors for PD in a South Asian population.

Methods 674 PD subjects predominantly with age of onset ≤ 50 years (encompassing juvenile, young, or early-onset PD) were recruited from 10 specialty movement disorder centers across India over a 2-year period. 1,376 control subjects were selected from the reference population GenomeAsia, Phase 2. We performed various case-only and case-control genetic analyses for PD diagnosis and age of onset.

Results A genome-wide significant signal for PD diagnosis was identified in the SNCA region, strongly colocalizing with SNCA region signal from European PD GWAS. PD cases with pathogenic mutations in PD genes exhibited, on average, lower PD polygenic risk scores than PD cases lacking any PD gene mutations. Gene burden studies of rare, predicted deleterious variants identified BSN, encoding the presynaptic protein Bassoon that has been previously associated with neurodegenerative disease.

Conclusions This study constitutes the largest genetic investigation of PD in a South Asian population to date. Future work should seek to expand sample numbers in this population to enable improved statistical power to detect PD genes in this understudied group.

Introduction

Parkinson’s disease (PD) is a progressive movement disorder characterized by dopaminergic neuron loss in the substantia nigra. Aging, environmental factors, and genetics contribute to PD risk and progression1. To date, common variants in 90 independent risk loci, discovered via genome-wide association studies (GWAS), and rare variants in more than 20 genes have been associated with PD2,3. However, most genetic studies of PD have been conducted in individuals of European ancestry4. Recent large-scale PD genetics studies in Latinos5 and East Asians6,7 have both confirmed the cross-ancestry presence of well-established European PD GWAS signals (i.e. in the SNCA region) and also have provided evidence for the existence of ancestry-specific genetic drivers of PD. However, South Asian populations remain particularly underrepresented in PD genetics studies4 and in human genetics studies generally, comprising ∼25% of the world population but only ∼1% of individuals assessed overall in GWAS (as of 2018)8. Previous genetics studies of PD in South Asians have been limited by small sample size, a lack of genome-wide interrogation, and/or case-only analyses9–17. The success of recent efforts to broaden PD genetic discovery beyond European populations5–7, and the high rates of founder effects and consanguinity in South Asians (themselves imparting distinct advantages for variant and gene discovery)18,19 render expanded genetics efforts in South Asians a promising opportunity to better understand the genetic basis of PD.

Focusing on individuals with young-(YOPD; typically defined as age of onset, or AoO, ≤ 40) or early-onset (EOPD; 40 < AoO ≤ 50) forms of PD, impart a greater probability of success, because younger age of onset is associated with a greater genetic contribution to disease3,20. These results are also relevant for late-onset PD (LOPD; defined herein as AoO > 50) wherever rarer variants associated with earlier-onset or more monogenic PD overlap signals detected via GWAS studies of common variants in LOPD. Several such “pleomorphic” risk loci/genes are known already: SNCA, GBA1, LRRK2, and VPS13C3,21. Previous genetic studies of early-onset PD have been conducted in a wide variety of ancestral populations, including European22–26, East Asian20,27–35, and South Asian12,16,17 groups. However, most of these early-onset studies were limited by the same factors that have hindered South Asian PD genetics studies generally: small sample sizes, analysis of only a small number of previously known PD genes, and/or lack of control populations to allow for genetic association studies.

In the present study, we address these limitations by conducting the largest genetic study of PD in a South Asian population to date, including genome-wide and case-control discovery approaches. Specifically, we jointly profiled 674 PD subjects predominantly with earlier onset PD from the Parkinson’s Research Alliance of India (PRAI) cohort via genotype array and whole exome sequencing (WES), expanding upon our previous study of 99 subjects from this cohort17. Combining these PD data with 1,376 ancestry-matched controls profiled via whole genome sequencing (WGS) and derived from a novel South Asian reference population18,19, we performed several studies quantifying variants and genes associated with PD diagnosis and age of onset, utilizing both common and rare variants.

Methods

674 subjects with PD were recruited through a network of 10 specialty movement disorder centers and/or neurology clinics located across India. PD diagnosis was made via modified UK Brain Bank criteria, and the study was approved by the Ethics Committee of Vikram Hospital (Bengaluru, India) Private Limited; additional details can be found in the previous clinical characterization of this cohort36. 1,376 ancestry-matched controls with available WGS data were selected from the GenomeAsia, Phase 2 (GAsPh2) reference cohort18,19.

Common variants were profiled with the South Asian Research Genotyping Array for Medicine (SARGAM)19, imputed with 6,461 samples from a GAsPh2 reference panel18,19, and subjected to stringent quality control (see Supporting Information Methods). Several steps were taken to identify and address genotype measurement error stemming from use of the first iteration of the SARGAM platform (see Supporting Information Methods). A GWAS of PD risk was performed for 4,499,703 common variants (study minor allele frequency or MAF > 5%) by comparing sporadic PD cases (N = 515) to matched controls (N = 1,363), adjusting for sex, age, and the first 10 principal components to account for population stratification. A second association analysis of PD AoO (N = 484 sporadic PD cases) was performed similarly. A significance level of P < 5 × 10−8 was used to define genome-wide significance.

WES data was generated for 576 subjects with PD. Variants were called with Sentieon’s GATK best practices germline pipeline and combined with results from WGS data for an additional 92 cases. Using American College of Medical Genetics (ACMG) criteria, we identified pathogenic, likely pathogenic, or variants of uncertain significance (VUS) in known PD genes in these 668 cases, as described previously (see Supporting Information Methods)17. For one identified variant in GBA1, p.Ser164Arg, we compared levels and activity37 of the corresponding enzyme β-Glucocerebrosidase (GCase) in homozygous variant knock-in (KI), GBA1 knock-out (KO), and wild-type (WT) A549 cells (see Supporting Information Methods).

We established an ancestry-normalized polygenic risk score (PRS) using variant-specific weights from a recent GWAS study of PD in Europeans2, and an adjustment procedure to apply this score in a South Asian ancestry sample, as described previously38. PRS distributions were compared between cases and controls and within PD cases defined by presence of a PD mutation, with a significance level of P < 0.05.

To perform gene burden studies of rare variants (MAF < 1%) between cases (WES, N = 633) and controls (WGS, N = 1363), we first accounted for case-control platform differences by implementing a previously described calibration procedure informed by synonymous variants (see Supporting Information Methods)39. Selected parameters from the calibration procedure were then used in gene burden testing of qualifying variants (QVs) defined based on their predicted functional consequences. Case-control analysis was performed using Fisher’s exact test with multiple testing correction applied across all tests (see Supporting Information Methods).

Data sharing

Summary statistics for PD diagnosis GWAS and PD age of onset GWAS are available on Zenodo (https://doi.org/10.5281/zenodo.8436983). Analysis code is also available on Zenodo (https://doi.org/10.5281/zenodo.10011733). Summary statistics for the gene burden analyses are available in the supplementary tables. Access to GAsPh2 data can be obtained through https://browser.genomeasia100k.org/#tid=access. Access to European PD GWAS summary statistics including the 23andme discovery dataset can be obtained through https://research.23andme.com/collaborate/. Chromatin interaction data from CNS cell types can be obtained through https://genome.ucsc.edu/s/nottalexi/glassLab_BrainCellTypes_hg19.

Results

Data and Analysis Overview

Genetic data was available on a total of 674 PD subjects, encompassing sporadic PD (N = 532; 78.9%) and familial PD (N = 142; 21.1%). Of these 674 subjects, genotype array data was generated on 659 subjects. WGS data was available for a pilot cohort of 92 subjects, as described previously17, and we generated WES data on 576 subjects. Control subjects were selected from available WGS from GAsPh2; a total of 1,376 subjects met criteria for ancestry matching (Fig. 1, Fig. S1). AoO was available for 651 PD subjects (Table S1), and used to assign them to diagnostic groups, as previously described36: 23 juvenile-onset PD (JOPD; AoO ≤ 20 years, 3.5%), 313 YOPD (20 < AoO ≤ 40, 48.1%) and 289 EOPD (40 < AoO ≤ 50, 44.4%) and 26 LOPD (AoO > 50, 4.0%). Case and control genetic data were combined to define genetic drivers of PD, encompassing common variant-(minor allele frequency, MAF > 5%) and/or rare variant (MAF < 1%) - based studies (Fig. 1).

Figure 1
Figure 1 Study overview.

A total of 674 PD patients, derived from the PRAI cohort, had some form of genetic data. Of this group, 659 patients had available SARGAM (genotyping array) data. NGS data consisted of 92 subjects with WGS data17 and 576 subjects with WES data, which were merged together to effectively create WES data on 668 samples (see Materials and Methods). A total of 1,376 control subjects were used after ancestry-matched to PD cases; these subjects were derived from WGS samples from GAsPh2. 5 major analyses were completed, each of which included common (MAF > 5%) and/or rare variants (MAF < 1%) and different genetic data types from cases and/or controls.

A genome-wide significant signal for PD diagnosis in the SNCA region mirrors signal in Europeans

After implementing a procedure to identify and correct genotyping errors from the novel SARGAM platform (see Supporting Information Methods, Figs. S2-S3), we conducted a GWAS of PD diagnosis using genetic data from sporadic PD samples (Ncases = 515) and ancestry matched controls (Ncontrols = 1,363; Fig. S1). Two loci reached genome-wide significance, including one containing the canonical PD gene SNCA (Fig. 2A). The lead SNCA region SNP in both the PRAI cohort and the latest GWAS of PD diagnosis in European samples2, rs356182 (pPRAI = 4.11E-11), is harbored in a neuronal enhancer element that physically interacts with the SNCA promoter. This enhancer-promoter loop is not observed in the other cell types (microglia and oligodendrocytes) assayed for chromatin interactions in the same dataset40 (Fig. 2B). The PRAI PD and European PD SNCA signals are highly colocalized (coloc41 posterior probability = 1; Fig. 2C) and show a consistent direction of effect (Fig. 2D). The other genome-wide significant locus harbored a lead SNP (rs59330234; p = 1.02E-8) in an intron of CCDC85A and has not been previously observed in European PD GWAS. Overall, there was very little concordance genome-wide between PRAI PD GWAS and European PD GWAS signal (excluding the SNCA region), even when considering loci reaching a suggestive (p < 1E–5) significance threshold in the PRAI PD GWAS (Fig. S4). In our (case-only) GWAS of AoO in sporadic PD samples (Ncases = 484), no genome-wide significant signals were observed (Figs. S2 and S5). Summary statistics for both the PD diagnosis and age of onset GWAS analyses are available (see Data Sharing).

Figure 2
Figure 2 PD Diagnosis GWAS results and concordance of SNCA region signals in South Asian and European populations

(A) Manhattan plot for PD diagnosis in the joint PRAI (PD cases) - GAsPh2 (controls) cohort for PD diagnosis. Lighter dashed line indicates suggestive significance (p = 1E-5) and darker dashed line indicates genome-wide significance (p = 5E-8). (B) Regional association plots in SNCA region for South Asian cohort in current study (top panel) and Europeans (International Parkinson’s Disease Genomics Consortium, “IPDGC”, second panel from top) from Nalls et al2. Points are colored by r2 value in study cohort and 1000 Genomes European samples, respectively. Also plotted are chromatin interaction loops identified in microglia, neurons, and oligodendrocytes40 in this region and overlapping a SNP reaching suggestive significance from the current study. (C) + (D) Comparison of PRAI-GAsPh2 cohort and SNCA region signals via p-value (C) and effect size and direction (via log odds ratio; D).

PRKN variants are a major driver of PD across diagnostic groups

We combined WGS data from a pilot cohort of 92 subjects (previously described17) with newly generated WES data for 576 additional individuals (Fig. S6) and assessed them for the presence of pathogenic, likely pathogenic, risk variants, or VUS in a list of 64 previously-associated PD genes by applying ACMG criteria (see Supporting Information Methods and Figs. S7-S8). Among the 649 subjects with available AoO information, pathogenic or likely pathogenic variants (Table S2) were more frequently observed in cases with an earlier age of onset, ranging from 47.6% of JOPD cases (10 of 21 cases), 10.2% of YOPD cases (32 of 313), 3.8% of EOPD cases (11 of 289), and none in LOPD cases (0 of 26; Fig. 3A). The total diagnostic yield, defined as the percentage of cases having a pathogenic or likely pathogenic variant in a PD gene, was 8.1%. Of the 33 J/Y/EOPD cases with a pathogenic variant, 22 cases harbored a PRKN mutation (Fig. 3B). Almost all pathogenic PRKN variants were homozygous deletions (20/22 cases) which were concentrated in the 3rd and 4th exons of PRKN (Fig. S9). No clinical features displayed strong differences between PRKN pathogenic variant carriers and non-carriers (Table S3). Pathogenic mutations were also observed in PINK1 (Ncases = 4), CHCHD2 and VPS13C (both Ncases = 2), and ATP13A2, PLA2G6, and PRRT2 (all Ncases = 1; Fig. 3B). Likely pathogenic mutations were observed in PLA2G6 (Ncases = 6), PRKN (Ncases = 6), CHCHD2 and MAPT (both Ncases = 2), and GCH1, PINK1, SYNJ1, and WDR45 (all Ncases = 1; Table S4). Next, we broadened considered variants to include risk variants and VUS, and tabulated occurrences of these variants specifically in a list of 21 genes previously described as causal for PD3. GBA1 mutations were the most frequent in the cohort (Fig. 3C), including variants observed across multiple individuals, all in heterozygous form (Table S4), e.g. p.Ser164Arg (Ncases = 10), p.Asp448His (Ncases = 6), p.Arg170His, p.Ser276Phe, and p.Ala495Pro (all Ncases = 3), and p.Arg434Cys (Ncases = 2). GBA1 p.Ser164Arg was observed in the cohort with a high frequency, but had not been functionally characterized before. We investigated the impact of this variant in a human cell culture model (see Supporting Information Methods). Homozygous p.Ser164Arg KI A549 cells exhibited significantly reduced GCase levels (Fig. 3D) and activity (Fig. 3E) compared to wild-type cells, at magnitudes akin to GBA1 KO cells. A list of all identified variants in J/Y/EOPD cases, across all 64 PD-related genes, including the counts of heterozygous and homozygous carriers in cases and controls (when available) is available as Table S4.

Figure 3
Figure 3 Diagnostic pipeline results

(A) Distribution of mutation type in known PD–related genes, by diagnostic group and familial/sporadic PD status. (B) Distribution of genes affected by pathogenic variants, by diagnostic group and familial/sporadic PD status. (C) Counts of all pathogenic, likely pathogenic, risk, and VUS variants identified in 21 genes previously described as monogenic for PD. (D) Representative Western blot and quantification of total cellular GCase levels from A549 GBA1-wild type, GBA1-KO, and three clonal GBA1-S164R variant cell lines. Bars represent mean GCase levels (normalized to beta-actin loading control within each sample) normalized to the mean WT value across all replicates, and points represent values from individual replicates. (Error bars: SEM [n = 3]; **** = p<0.0001, ANOVA with Dunnett’s multiple comparisons test). (E) Flow cytometry measurement of GCase activity in A549 GBA1-wild type, GBA1-KO, and three clonal GBA1-S164R variant cell lines. Bars represent mean GCase activity normalized to the mean WT value across all replicates, and points represent values from individual replicates. (Error bars: SEM [n = 3]; **** = p<0.0001, ANOVA with Dunnett’s multiple comparisons test).

PD gene pathogenic variant carriers have significantly reduced polygenic burden relative to non-carriers

We next leveraged our imputed SARGAM microarray data to construct a PRS for PD, using weights from the most recent PD GWAS study in European samples2, and a previously described method to project a PRS informed by weights of one ancestral population to another38. PD cases had a significantly higher PRS score than controls (Wilcoxon test, p-value = 1.94E-11; Fig. 4A). To further quantify this difference, we binned samples into deciles (using breakpoints defined in the controls, as previously described38) and calculated decile-specific odds ratios for disease, adjusted for age and sex. Membership in the tenth decile was significantly associated with increased odds of disease relative to a joint 5th and 6th decile comparison group38 (Odds Ratio = 1.61, p = 0.006; Fig. 4B).

Figure 4
Figure 4 Polygenic risk scores and their integration with diagnostic categories

(A) Distributions of ancestry-normalized PRSs in PD cases and controls. (B) Odds ratios for disease in PRS-defined deciles in the control group, relative to deciles 5 and 6. (C) Distribution of PRSs in controls and diagnostic pipeline-defined categories for PD cases.

Next, we compared the PRS between controls and several groups of PD cases defined by presence of a PD-related mutation (“VUS”, “Likely Pathogenic”, or “Pathogenic”) or lack thereof (“None”), and within those same PD case groups. The “None” (p = 5.72E-9) and “VUS” (p = 4.08E-8) groups exhibited significantly higher PRS than the control group, adjusted for age and sex. Interestingly, the “Pathogenic” group exhibited significantly lower PRS compared to the “None” group (p = 1.49E-3; Fig. 4C), adjusted for age and sex. No other comparisons achieved nominal (p < 0.05) statistical significance.

Gene burden analyses of rare, predicted deleterious variants demonstrate suggestive evidence linking BSN to PD

We next performed gene burden analysis with WES data from PD cases (N = 633) and WGS data from controls (N = 1363; Fig. S7). Using a calibration procedure informed by synonymous variants39, we derived coverage and variant quality filters to adjust for platform differences between cases and controls (see Supporting Information Methods and Fig. S10). After imposing these filters, and further restricting to variants with MAF < 1% across all ancestral populations, a total of 134,871 variants from the PD case WES data and 195,903 variants from the control WGS data could be evaluated. We selected QVs according to several criteria (see Supporting Information Methods)41. For example, the most stringent analysis considered 91 genes with predicted loss-of-function (LoF) variants. Additional, more permissive definitions of QVs, infomed by the in silico prediction tools Combined Annotation Depletion-Dependent (CADD)43 and Polymorphism Phenotyping v2 (PolyPhen2)43 were implemented, along with a minimum carrier count across cases and controls of 5 samples, resulting in different gene counts corresponding to stringency: 1) 2,854 genes carrying LoF + PolyPhen2 “probably damaging” variants, 2) 4,572 genes with LoF + PolyPhen2 “probably damaging” or “possibly damaging” variants, 3) 7,627 genes with variants with CADD score > 20, and 4) 9,209 genes with broadly “Functional Rare” variants (encompassing missense, inframe, frameshift, and nonsense variants) for analysis (Fig. 5A). After removing the LoF only tests (low gene count) and gene-QV criterion combinations potentially impacted by confounding (see Supporting Information Methods), a total of 21,412 burden tests were performed, with many genes tested multiple times under different QV criteria (Table S5). No tests surpassed a Bonferroni correction threshold (p = 2.33E-6), likely due to the limited statistical power afforded by the cohort (Fig. S11). However, several genes achieved marginal significance (p < 5E-4) across multiple QV criteria, demonstrating association with PD not entirely restricted to a single set of variants (Fig. 5B). For example, we observed suggestive associations for BSN under the “Functional Rare” QV criterion (p = 2.57E-5, the most significant result observed overall), but also the LoF + PolyPhen2 “probably damaging” or “possibly damaging” (p = 3.56E-4) and the CADD score > 20 (p = 3.88E-4) criteria. BSN encodes the Bassoon protein and is highly expressed in the brain (Gene-Tissue Expression (GTEx) Project v845, Fig. S12A). Interestingly, the variants driving the differential gene burden signal were localized towards the C-terminus (a similar region of the protein implicated in a recent genetics study of progressive supranuclear palsy-like syndrome, or PSP46) and were never observed in controls (Fig. S12B).

Figure 5
Figure 5 Gene burden analyses of rare variants

(A) Number of genes available for analysis as a function of the number of carriers (across cases and controls), for several different definitions of a qualifying variant (QV). Genes with carrier count ≥ 5 were included in gene burden testing. LoF-only analysis not performed due to low number of genes available that met the minimum carrier count threshold. (B) QQ-plot displaying results of gene burden studies across QV criteria. Methods and plots adapted from previous gene burden analysis of Alzheimer’s Disease42.

Discussion

We have conducted the largest genetics study of PD in a South Asian population to date and to our knowledge, a study focused on PD patients with an earlier-onset phenotype (AoO ≤ 50 years) and utilizing a diverse set of genetic data types (genotype array and WES in cases, WGS in controls) and analysis approaches (including common and rare variation, variant and gene-level inferences, and PRS). By comparing these PD patients with ancestry-matched controls we identified a genome-wide significant signal in the region of SNCA, a canonical PD gene, that recapitulates the European PD GWAS signal from the same region. Approximately 8% of PD patients harbor a pathogenic or likely pathogenic mutation in a previously known PD gene (predominantly, PRKN), and carriers of such pathogenic variants have a significantly reduced polygenic burden relative to PD cases without any PD gene mutations. Finally, through case-control comparison of rare variants at the gene level, we demonstrate suggestive evidence for an association of BSN and increased PD risk.

The PD diagnosis GWAS yielded two genome-wide significant signals. The first was an intronic SNP in CCDC85A, implicating this gene, which is highly expressed in the brain45, for the first time in PD. The second signal, in the SNCA region, was driven by the lead SNP rs356182, which matches that of one of the two independent SNCA region signals reported in a recent European PD GWAS2, and has been functionally characterized extensively47,48. Neuronal chromatin interaction data40 indicates that rs356182 variants impact SNCA gene expression in this cell type as shown previously47,49, and consistent with the proposed impact of this SNP on upregulating a SNCA 5’ untranslated region transcript isoform in frontal cortex47. This SNP has also been reported as the lead SNP in SNCA region signals in recent PD studies in Latino5 and Chinese7 populations, and in a multi-ancestry meta-analysis including individuals of European, East Asian, Latin American, and African ancestry50. Our results extend the relevance of rs356182 to an additional ancestral population, South Asians, and to earlier onset forms of PD. Beyond this SNCA region signal, no regions demonstrating suggestive signal for PD diagnosis (p < 1E-5) in our South Asian cohort showed evidence of association in European samples. Many of these suggestive signals may turn out to be false positives, but some might also correspond to common variation associated with PD risk that is unique to South Asians populations. Future, larger studies of PD in South Asians will be better powered to distinguish these possibilities.

The examination of known PD genes revealed strong rare variant drivers, and the diagnostic yield (∼8%) in our patient population matches that of similar studies in East Asian populations20,33,34 (7.9% to 11.6%). Homozygous deletions in PRKN comprised the majority of pathogenic variants identified in this cohort, consistent with previous studies demonstrating that homozygous or compound heterozygous PRKN variants are the major driver of young- or early- onset PD51–53. Furthermore, the identified PRKN deletions were concentrated in exons 3 and 4 of PRKN, consistent with a recent comprehensive assessment of PRKN mutation frequency54. We also detected a number of GBA1 variants in PD cases. However, all were classified via ACMG criteria as being “Risk” or “VUS” variants, and therefore their roles as causal genetic drivers remains to be elucidated. One such GBA1 variant, p.Ser164Arg, was observed in 10 (unique) PD cases out of 668 (1.5%). A549 cells carrying a homozygous p.Ser164Arg KI displayed a marked decrease in both GCase levels and activity, effectively phenocopying gene KO cells. This result, the first functional evidence for p.Ser164Arg presented to date, is in line with the homozygous occurrence of this variant in Gaucher disease (GD) patients55–57, whose GCase levels and/or activity are severely reduced. Data from gnomAD58 and a recent comprehensive investigation of GBA1 variant frequencies59 only report occurrence of p.Ser164Arg in South Asians. Along these lines, all previous reports of p.Ser164Arg in GD patients and in PD patients16 have come from Indian populations. In the latter instance, p.Ser164Arg was observed in 3 of 250 (1.2%) EOPD patients, a frequency similar to that observed in this cohort. Two additional GBA1 variants, p.Ser276Phe and p.Arg434Cys, observed in 3 and 2 PD cases respectively, also appear to be only present in South Asian individuals58,59. Given this evidence, future functional studies should ascertain if these variants similarly reduce GCase levels and activity as demonstrated here for p.Ser164Arg, and whether they pinpoint novel aspects of GBA1 biology.

We observed a modestly increased polygenic burden in PD cases relative to controls. When considering rare and common variants jointly, we observed this difference was driven by PD cases lacking a PD gene mutation (“None” group) and those carrying only a VUS. Recent functional studies of induced pluripotent stem cell (iPSC) derived dopaminergic neurons from EOPD patient samples without known PD gene mutations displayed increased α-synuclein accumulation and reduced lysosomal membrane protein abundance compared to control iPSC lines60. Our findings are consistent with these results in arguing for the existence of genetic contributions to EOPD beyond rare mutations in known PD risk genes. It was also observed in our study that PD cases with a pathogenic PD gene mutation had significantly reduced polygenic burden relative to PD cases lacking a PD gene mutation. Similar findings of reduced polygenic burden in mutation carrier cases relative to non-carrier cases have also been observed for phenotypes such as autism spectrum disorder61, breast and colorectal cancer, diabetes, osteoporosis, and short stature62. These observations are consistent with a liability threshold model of disease by which rare, pathogenic variant carriers require a smaller contribution from common variants to reach a disease diagnosis, compared to non-carriers63. However, because only a limited number of studies have examined the interplay between rare and common variant genetic drivers in PD specifically64,65, this is an important area of future work.

Gene-based burden analyses revealed a suggestive association for BSN, which has had limited connection to PD previously. As mentioned above, BSN encodes Bassoon, a large synaptic scaffolding protein that plays a critical role in assembling the presynaptic active zone, where synaptic vesicle exocytosis primarily occurs66. Bassoon is also a key regulator of presynaptic proteostasis, controlling degradative pathways that have been extensively implicated in PD, such as autophagy, notably in concert with the bona fide PD-linked protein Parkin67. Recent human genetic and preclinical studies have also directly implicated Bassoon in neurodegenerative diseases, by identifying mutations in individuals with PSP46, a primary tauopathy, and showing a role of Bassoon in the seeding and toxicity of tau pathology in mice68. Our findings suggest a broader involvement of Bassoon in neurodegenerative disease, and future work should look to evaluate the functional impact of the individual BSN variants driving the gene-level association we observed.

Our study has limitations. The total sample size for our study population consisted of roughly 2,000 cases and controls, limiting the statistical power for all variant- and gene-based discovery analyses. Moreover, given this sample size, and the rarity of early-onset PD cases generally, we were not able to derive a suitable validation cohort for our findings. In some instances, our findings are validated by previous studies in other cohorts or populations, such as the SNCA region signal for PD diagnosis or high frequency of PRKN deletions in our cohort. Other novel findings, such as the reduced polygenic burden observed in PD cases carrying pathogenic variants, or the suggestive association of BSN, require replication in additional South Asian PD cohorts and in PD cohorts more generally.

The Global Parkinson’s Genetics Program (GP2) was recently initiated to directly address the limited insights into PD disease etiology and limited potential for equitable personalized medicine due to the strong bias of European ancestry samples in previous PD genetics studies4,69. This work is an initial step in addressing this knowledge gap in a South Asian population. With this effort and parallel GP2-fueled PD genetics studies in India70, the knowledge of PD genetic risk factors in South Asians, and therefore understanding of PD genetics as whole, should increase substantially in the coming years.

Data Availability

Summary statistics for PD diagnosis GWAS and PD age of onset GWAS are available on Zenodo (https://doi.org/10.5281/zenodo.8436983). Analysis code is also available on Zenodo (https://doi.org/10.5281/zenodo.10011733). Summary statistics for the gene burden analyses are available in the supplementary tables. Access to GAsPh2 data can be obtained through https://browser.genomeasia100k.org/#tid=access. Access to European PD GWAS summary statistics including the 23andme discovery dataset can be obtained through https://research.23andme.com/collaborate/. Chromatin interaction data from CNS cell types can be obtained through https://genome.ucsc.edu/s/nottalexi/glassLab_BrainCellTypes_hg19.

Authors’ Roles

Design: SVA, PLK, RM, OBD, MAR, DEL, AGH, AP, TS, RG, VLR

Execution: All authors

Analysis: SVA, PLK, RM, TSG, MO, SN, AK, PLS, MD, AB, CS, NM, OBD, MAR, TS, RG

Writing: SVA, PLK, RM, OBD, MAR, GDP, TS, RG

Editing of final version of the manuscript: SVA, PLK, RM, OBD, MAR, DEL, AGH, GDP, TS, RG, VLR

Acknowledgements

We thank the study participants for their time and contribution of genetic and clinical information to perform this study. We thank J. Fah Sathirapongsasuti for his contribution to the data analysis and study design.

Footnotes

  • ↵† Author is deceased.

  • Relevant conflicts of interest/financial disclosures: SVA, OBD, MAR, DEL, AGH, GDP, and TS are full-time employees and shareholders of Denali Therapeutics. RM, TSG, SM, MO, SN, AK, PLS, MD, AB, CS, NM, ASP, RG, and VLR are current or former employees of MedGenome Labs.

  • Funding: No funding was received towards this work.

  • Updated per reviewer feedback and to include new data characterizing the GBA1 p.Ser164Arg variant.

References

  1. 1.↵
    Dorsey ER, Bloem BR. The Parkinson Pandemic—A Call to Action. JAMA Neurol. 2018;75(1):9– 10.
    OpenUrl
  2. 2.↵
    Nalls MA, Blauwendraat C, Vallerga CL, et al. Identification of novel risk loci, causal insights, and heritable risk for Parkinson’s disease: a meta-analysis of genome-wide association studies. Lancet Neurol. 2019;18(12):1091–1102.
    OpenUrlCrossRefPubMed
  3. 3.↵
    Blauwendraat C, Nalls MA, Singleton AB. The genetic architecture of Parkinson’s disease. Lancet Neurol. 2020;19(2):170–178.
    OpenUrlCrossRefPubMed
  4. 4.↵
    Schumacher-Schuh AF, Bieger A, Okunoye O, et al. Underrepresented Populations in Parkinson’s Genetics Research: Current Landscape and Future Directions. Mov. Disord. Off. J. Mov. Disord. Soc. 2022;37(8):1593–1604.
    OpenUrlCrossRef
  5. 5.↵
    Loesch DP, Horimoto ARVR, Heilbron K, et al. Characterizing the Genetic Architecture of Parkinson’s Disease in Latinos. Ann. Neurol. 2021;90(3):353–365.
    OpenUrlCrossRef
  6. 6.↵
    Foo JN, Chew EGY, Chung SJ, et al. Identification of Risk Loci for Parkinson Disease in Asians and Comparison of Risk Between Asians and Europeans: A Genome-Wide Association Study. JAMA Neurol. 2020;77(6):746–754.
    OpenUrl
  7. 7.↵
    Pan H, Liu Z, Ma J, et al. Genome-wide association study using whole-genome sequencing identifies risk loci for Parkinson’s disease in Chinese population. Npj Park. Dis. 2023;9(1):1–11.
    OpenUrl
  8. 8.↵
    Martin AR, Kanai M, Kamatani Y, et al. Clinical use of current polygenic risk scores may exacerbate health disparities. Nat. Genet. 2019;51(4):584–591.
    OpenUrlCrossRefPubMed
  9. 9.↵
    Punia S, Behari M, Govindappa ST, et al. Absence/rarity of commonly reported LRRK2 mutations in Indian Parkinson’s disease patients. Neurosci. Lett. 2006;409(2):83–88.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.
    Sanyal J, Sarkar B, Ojha S, et al. Absence of commonly reported leucine-rich repeat kinase 2 mutations in Eastern Indian Parkinson’s disease patients. Genet. Test. Mol. Biomark. 2010;14(5):691–694.
    OpenUrlCrossRef
  11. 11.
    Vijayan B, Gopala S, Kishore A. LRRK2 G2019S mutation does not contribute to Parkinson’s disease in South India. Neurol. India 2011;59(2):157–160.
    OpenUrlCrossRefPubMed
  12. 12.↵
    Padmaja MV, Jayaraman M, Srinivasan AV, et al. PARK2 gene mutations in early onset Parkinson’s disease patients of South India. Neurosci. Lett. 2012;523(2):145–147.
    OpenUrlCrossRefPubMed
  13. 13.
    Vishwanathan Padmaja M, Jayaraman M, Srinivasan AV, et al. The SNCA (A53T, A30P, E46K) and LRRK2 (G2019S) mutations are rare cause of Parkinson’s disease in South Indian patients. Parkinsonism Relat. Disord. 2012;18(6):801–802.
    OpenUrl
  14. 14.
    Yadav R, Kapoor S, Madhukar M, et al. Genetic analysis of the glucocerebrosidase gene in South Indian patients with Parkinson’s disease. Neurol. India 2018;66(6):1649–1654.
    OpenUrlCrossRef
  15. 15.
    Kishore A, Ashok Kumar Sreelatha A, Sturm M, et al. Understanding the role of genetic variability in LRRK2 in Indian population. Mov. Disord. Off. J. Mov. Disord. Soc. 2019;34(4):496–505.
    OpenUrl
  16. 16.↵
    Kumar S, Yadav N, Pandey S, et al. Novel and reported variants in Parkinson’s disease genes confer high disease burden among Indians. Parkinsonism Relat. Disord. 2020;78:46–52.
    OpenUrlCrossRef
  17. 17.↵
    Kukkle PL, Geetha TS, Chaudhary R, et al. Genome-Wide Polygenic Score Predicts Large Number of High Risk Individuals in Monogenic Undiagnosed Young Onset Parkinson’s Disease Patients from India. Adv. Biol. 2022;e2101326.
  18. 18.↵
    Wall JD, Stawiski EW, Ratan A, et al. The GenomeAsia 100K Project enables genetic discoveries across Asia. Nature 2019;576(7785):106–111.
  19. 19.↵
    Wall JD, Sathirapongsasuti JF, Gupta R, et al. South Asian medical cohorts reveal strong founder effects and high rates of homozygosity. Nat. Commun. 2023;14(1):3377.
    OpenUrl
  20. 20.↵
    Zhao Y, Qin L, Pan H, et al. The role of genetics in Parkinson’s disease: a large cohort study in Chinese mainland population. Brain J. Neurol. 2020;143(7):2220–2234.
    OpenUrlCrossRef
  21. 21.↵
    Singleton A, Hardy J. A generalizable hypothesis for the genetic architecture of disease: pleomorphic risk loci. Hum. Mol. Genet. 2011;20(R2):R158–162.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    Siitonen A, Nalls MA, Hernández D, et al. Genetics of early-onset Parkinson’s disease in Finland: exome sequencing and genome-wide association study. Neurobiol. Aging 2017;53:195.e7–195.e10.
    OpenUrlCrossRef
  23. 23.
    Gustavsson EK, Trinh J, McKenzie M, et al. Genetic Identification in Early Onset Parkinsonism among Norwegian Patients. Mov. Disord. Clin. Pract. 2017;4(4):499–508.
    OpenUrl
  24. 24.
    Trinh J, Lohmann K, Baumann H, et al. Utility and implications of exome sequencing in early-onset Parkinson’s disease. Mov. Disord. Off. J. Mov. Disord. Soc. 2019;34(1):133–137.
    OpenUrl
  25. 25.
    Cristina T-P, Pablo M, Teresa PM, et al. A genetic analysis of a Spanish population with early onset Parkinson’s disease. PloS One 2020;15(9):e0238098.
    OpenUrlCrossRef
  26. 26.↵
    Olszewska DA, McCarthy A, Soto-Beasley AI, et al. PARKIN, PINK1, and DJ1 analysis in early-onset Parkinson’s disease in Ireland. Ir. J. Med. Sci. 2022;191(2):901–907.
    OpenUrlCrossRef
  27. 27.↵
    Lin C-H, Chen P-L, Tai C-H, et al. A clinical and genetic study of early-onset and familial parkinsonism in taiwan: An integrated approach combining gene dosage analysis and next-generation sequencing. Mov. Disord. 2019;34(4):506–515.
    OpenUrlCrossRef
  28. 28.
    Youn J, Lee C, Oh E, et al. Genetic variants of PARK genes in Korean patients with early-onset Parkinson’s disease. Neurobiol. Aging 2019;75:224.e9–224.e15.
    OpenUrlCrossRef
  29. 29.
    Li N, Wang L, Zhang J, et al. Whole-exome sequencing in early-onset Parkinson’s disease among ethnic Chinese. Neurobiol. Aging 2020;90:150.e5–150.e11.
    OpenUrlCrossRef
  30. 30.
    Jiang Y, Yu M, Chen J, et al. Parkin is the most common causative gene in a cohort of mainland Chinese patients with sporadic early-onset Parkinson’s disease. Brain Behav. 2020;10(9):e01765.
    OpenUrlCrossRef
  31. 31.
    Chen H, Jin Y-H, Xue Y-Y, et al. Novel ATP13A2 and PINK1 variants identified in Chinese patients with Parkinson’s disease by whole-exome sequencing. Neurosci. Lett. 2020;733:135075.
  32. 32.
    Chen Y-P, Gu X-J, Song W, et al. Rare Variants Analysis of Lysosomal Related Genes in Early-Onset and Familial Parkinson’s Disease in a Chinese Cohort. J. Park. Dis. 2021;11(4):1845–1855.
    OpenUrl
  33. 33.↵
    Hua P, Zhao Y, Zeng Q, et al. Genetic Analysis of Patients With Early-Onset Parkinson’s Disease in Eastern China. Front. Aging Neurosci. 2022;14:849462.
  34. 34.↵
    Chen Y-P, Yu S-H, Zhang G-H, et al. The mutation spectrum of Parkinson-disease-related genes in early-onset Parkinson’s disease in ethnic Chinese. Eur. J. Neurol. 2022;29(11):3218–3228.
    OpenUrlCrossRef
  35. 35.↵
    Do MD, Tran TN, Luong AB, et al. Clinical and genetic analysis of Vietnamese patients diagnosed with early-onset Parkinson’s disease. Brain Behav. 2023;e2950.
  36. 36.↵
    Kukkle PL, Goyal V, Geetha TS, et al. Clinical Study of 668 Indian Subjects with Juvenile, Young, and Early Onset Parkinson’s Disease. Can. J. Neurol. Sci. J. Can. Sci. Neurol. 2022;49(1):93–101.
    OpenUrlCrossRef
  37. 37.↵
    Deen MC, Zhu Y, Gros C, et al. A versatile fluorescence-quenched substrate for quantitative measurement of glucocerebrosidase activity within live cells. Proc. Natl. Acad. Sci. U. S. A. 2022;119(29):e2200553119.
    OpenUrl
  38. 38.↵
    Wang M, Menon R, Mishra S, et al. Validation of a Genome-Wide Polygenic Score for Coronary Artery Disease in South Asians. J. Am. Coll. Cardiol. 2020;76(6):703–714.
    OpenUrlFREE Full Text
  39. 39.↵
    Guo MH, Plummer L, Chan Y-M, et al. Burden Testing of Rare Variants Identified through Exome Sequencing via Publicly Available Control Data. Am. J. Hum. Genet. 2018;103(4):522–534.
    OpenUrlCrossRefPubMed
  40. 40.↵
    Nott A, Holtman IR, Coufal NG, et al. Brain cell type-specific enhancer-promoter interactome maps and disease-risk association. Science 2019;366(6469):1134–1139.
  41. 41.↵
    Giambartolomei C, Vukcevic D, Schadt EE, et al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. 2014;10(5):e1004383.
    OpenUrlCrossRefPubMed
  42. 42.↵
    Holstege H, Hulsman M, Charbonnier C, et al. Exome sequencing identifies rare damaging variants in ATP8B4 and ABCA1 as risk factors for Alzheimer’s disease. Nat. Genet. 2022;54(12):1786– 1794.
    OpenUrlCrossRef
  43. 43.↵
    Rentzsch P, Witten D, Cooper GM, et al. CADD: predicting the deleteriousness of variants throughout the human genome. Nucleic Acids Res. 2019;47(D1):D886–D894.
    OpenUrlCrossRefPubMed
  44. 44.
    Adzhubei IA, Schmidt S, Peshkin L, et al. A method and server for predicting damaging missense mutations. Nat. Methods 2010;7(4):248–249.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    GTEx Consortium. The Genotype-Tissue Expression (GTEx) project. Nat. Genet. 2013;45(6):580– 585.
    OpenUrlCrossRefPubMed
  46. 46.↵
    Yabe I, Yaguchi H, Kato Y, et al. Mutations in bassoon in individuals with familial and sporadic progressive supranuclear palsy-like syndrome. Sci. Rep. 2018;8(1):819.
    OpenUrlCrossRef
  47. 47.↵
    Pihlstrøm L, Blauwendraat C, Cappelletti C, et al. A Comprehensive Analysis of SNCA−Related Genetic Risk in Sporadic Parkinson Disease. Ann. Neurol. 2018;84(1):117–129.
    OpenUrlCrossRef
  48. 48.↵
    Prahl JD, Pierce SE, van der Schans EJC, et al. The Parkinson’s disease variant rs356182 regulates neuronal differentiation independently from alpha-synuclein. Hum. Mol. Genet. 2023;32(1):1–14.
    OpenUrlCrossRef
  49. 49.↵
    Schilder BM, Raj T. Fine-mapping of Parkinson’s disease susceptibility loci identifies putative causal variants. Hum. Mol. Genet. 2022;31(6):888–900.
    OpenUrlCrossRef
  50. 50.↵
    Kim JJ, Vitale D, Otani DV, et al. Multi-ancestry genome-wide meta-analysis in Parkinson’s disease [Internet]. 2022;2022.08.04.22278432.[cited 2023 Feb 28] Available from: https://www.medrxiv.org/content/10.1101/2022.08.04.22278432v1
  51. 51.↵
    Alcalay RN, Caccappolo E, Mejia-Santana H, et al. Frequency of Known Mutations in Early-Onset Parkinson Disease: Implication for Genetic Counseling: The Consortium on Risk for Early Onset Parkinson Disease Study. Arch. Neurol. 2010;67(9):1116–1122.
    OpenUrlCrossRefPubMedWeb of Science
  52. 52.
    Bonifati V. Autosomal recessive parkinsonism. Parkinsonism Relat. Disord. 2012;18 Suppl 1:S4–6.
    OpenUrl
  53. 53.↵
    Kilarski LL, Pearson JP, Newsway V, et al. Systematic review and UK-based study of PARK2 (parkin), PINK1, PARK7 (DJ-1) and LRRK2 in early-onset Parkinson’s disease. Mov. Disord. Off. J. Mov. Disord. Soc. 2012;27(12):1522–1529.
    OpenUrl
  54. 54.↵
    Zhu W, Huang X, Yoon E, et al. Heterozygous PRKN mutations are common but do not increase the risk of Parkinson’s disease. Brain J. Neurol. 2022;145(6):2077–2091.
    OpenUrlCrossRef
  55. 55.↵
    Agarwal S, Lahiri K, Muranjan M, Solanki N. The face of lysosomal storage disorders in India: a need for early diagnosis. Indian J. Pediatr. 2015;82(6):525–529.
    OpenUrlCrossRef
  56. 56.
    Sheth J, Pancholi D, Mistri M, et al. Biochemical and molecular characterization of adult patients with type I Gaucher disease and carrier frequency analysis of Leu444Pro - a common Gaucher disease mutation in India. BMC Med. Genet. 2018;19(1):178.
    OpenUrlCrossRef
  57. 57.↵
    Thomas DC, Sharma S, Puri RD, et al. Lysosomal storage disorders: Novel and frequent pathogenic variants in a large cohort of Indian patients of Pompe, Fabry, Gaucher and Hurler disease. Clin. Biochem. 2021;89:14–37.
    OpenUrlCrossRef
  58. 58.↵
    Karczewski KJ, Francioli LC, Tiao G, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 2020;581(7809):434–443.
  59. 59.↵
    Parlar SC, Grenn FP, Kim JJ, et al. Classification of GBA1 Variants in Parkinson’s Disease: The GBA1-PD Browser. Mov. Disord. Off. J. Mov. Disord. Soc. 2023;
  60. 60.↵
    Laperle AH, Sances S, Yucer N, et al. iPSC modeling of young-onset Parkinson’s disease reveals a molecular signature of disease and novel therapeutic candidates. Nat. Med. 2020;26(2):289–299.
    OpenUrlCrossRef
  61. 61.↵
    Klei L, McClain LL, Mahjani B, et al. How rare and common risk variation jointly affect liability for autism spectrum disorder. Mol. Autism 2021;12(1):66.
    OpenUrlCrossRef
  62. 62.↵
    Lu T, Zhou S, Wu H, et al. Individuals with common diseases but with a low polygenic risk score could be prioritized for rare variant screening. Genet. Med. 2021;23(3):508–515.
    OpenUrl
  63. 63.↵
    Walsh R, Tadros R, Bezzina CR. When genetic burden reaches threshold. Eur. Heart J. 2020;41(39):3849–3855.
    OpenUrlCrossRef
  64. 64.↵
    Blauwendraat C, Reed X, Krohn L, et al. Genetic modifiers of risk and age at onset in GBA associated Parkinson’s disease and Lewy body dementia. Brain J. Neurol. 2020;143(1):234–248.
    OpenUrlCrossRef
  65. 65.↵
    Iwaki H, Blauwendraat C, Makarious MB, et al. Penetrance of Parkinson’s Disease in LRRK2 p.G2019S Carriers Is Modified by a Polygenic Risk Score. Mov. Disord. Off. J. Mov. Disord. Soc. 2020;35(5):774–780.
    OpenUrl
  66. 66.↵
    Gundelfinger ED, Reissner C, Garner CC. Role of Bassoon and Piccolo in Assembly and Molecular Organization of the Active Zone. Front. Synaptic Neurosci. 2015;7:19.
  67. 67.↵
    Gundelfinger ED, Karpova A, Pielot R, et al. Organization of Presynaptic Autophagy-Related Processes. Front. Synaptic Neurosci. 2022;14:829354.
  68. 68.↵
    Martinez P, Patel H, You Y, et al. Bassoon contributes to tau-seed propagation and neurotoxicity. Nat. Neurosci. 2022;25(12):1597–1607.
    OpenUrl
  69. 69.↵
    GP2: The Global Parkinson’s Genetics Program. Mov. Disord. 2021;36(4):842–851.
    OpenUrl
  70. 70.↵
    Rajan R, Divya KP, Kandadai RM, et al. Genetic Architecture of Parkinson’s Disease in the Indian Population: Harnessing Genetic Diversity to Address Critical Gaps in Parkinson’s Disease Research [Internet]. Front. Neurol. 2020;11[cited 2023 Feb 27] Available from: https://www.frontiersin.org/articles/10.3389/fneur.2020.00524
Back to top
Posted October 19, 2023.
Download PDF

Supplementary Material

Data/Code
Email
The genetic drivers of juvenile, young, and early-onset Parkinson’s Disease in India
Shan V Andrews, Prashanth L Kukkle, Ramesh Menon, Thenral S Geetha, Vinay Goyal, Rukmini Mridula Kandadai, Hrishikesh Kumar, Rupam Borgohain, Adreesh Mukherjee, Pettarusp M Wadia, Ravi Yadav, Soaham Desai, Niraj Kumar, Deepika Joshi, Sakthivel Murugan, Atanu Biswas, Pramod K Pal, Merina Oliver, Sandhya Nair, Anbu Kayalvizhi, Praveena L Samson, Manjari Deshmukh, Akshi Bassi, Charugulla Sandeep, Nitin Mandloi, Oliver B Davis, Melissa A Roberts, Dara E Leto, Anastasia G Henry, Gilbert Di Paolo, Uday Muthane, Shymal K Das, Andrew S Peterson, Thomas Sandmann, Ravi Gupta, Vedam L Ramprasad, Parkinson Research Alliance of India (PRAI)
medRxiv 2023.06.18.23291407; doi: https://doi.org/10.1101/2023.06.18.23291407
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area