The Lancet NeurologyGenetic & Genomic MedicinePeer-reviewed journal article
Genome-wide Association Identifies Novel Etiological Insights Associated with Parkinsons Disease in African and African Admixed Populations
Mie RizigiD, Sara Bandrés‐CigaiD, Mary B. MakariousiD, Oluwadamilola Omolara OjoiD, Peter Wild CreaiD, Oladunni AbioduniD, Kristin S. LevineiD, Sani AbubakariD, Charles AchoruiD, Dan VitaleiD, Olaleye Akinmola AdenijiiD, Osigwe Paul AgabiiD, Mathew J. KoretskyiD, Uchechi AgulannaiD, Deborah A Hall, Rufus AkinyemiiD, Tao XieiD, Mohammed Wulgo AliiD, Ejaz A. ShamimiD, Ifeyinwa Ani-Osheku, Mahesh PadmanabaniD, Ohwotemu Arigbodi, David G. StandaertiD, Abiodun Hamzat BelloiD, Marissa N. DeaniD, , Francis Ojini
Abstract
BackgroundUnderstanding the genetic mechanisms underlying diseases in ancestrally diverse populations is a critical step towards the realization of the global application of precision medicine. The African and African admixed populations enable mapping of complex traits given their greater levels of genetic diversity, extensive population substructure, and distinct linkage disequilibrium patterns. MethodsHere we perform a comprehensive genome-wide assessment of Parkinsons disease (PD) in 197,918 individuals (1,488 cases; 196,430 controls) of African and African admixed ancestry, characterizing population-specific risk, differential haplotype structure and admixture, coding and structural genetic variation and polygenic risk profiling. FindingsWe identified a novel common risk factor for PD and age at onset at the GBA1 locus (risk, rs3115534-G; OR=1.58, 95% CI = 1.37 - 1.80, P=2.397E-14; age at onset, BETA =-2.004, SE =0.57, P = 0.0005), that was found to be rare in non-African/African admixed populations. Downstream short- and long-read whole genome sequencing analyses did not reveal any coding or structural variant underlying the GWAS signal. However, we identified that this signal mediates PD risk via expression quantitative trait locus (eQTL) mechanisms. While previously identified GBA1 associated disease risk variants are coding mutations, here we suggest a novel functional mechanism consistent with a trend in decreasing glucocerebrosidase activity levels. Given the high population frequency of the underlying signal and the phenotypic characteristics of the homozygous carriers, we hypothesize that this variant may not cause Gaucher disease. Additionally, the prevalence of Gauchers disease in Africa is low. InterpretationThe present study identifies a novel African-ancestry genetic risk factor in GBA1 as a major mechanistic basis of PD in the African and African admixed populations. This striking result contrasts to previous work in Northern European populations, both in terms of mechanism and attributable risk. This finding highlights the importance of understanding population-specific genetic risk in complex diseases, a particularly crucial point as the field moves toward precision medicine in PD clinical trials and while recognizing the need for equitable inclusion of ancestrally diverse groups in such trials. Given the distinctive genetics of these underrepresented populations, their inclusion represents a valuable step towards insights into novel genetic determinants underlying PD etiology. This opens new avenues towards RNA-based and other therapeutic strategies aimed at reducing lifetime risk. Research in Context Evidence Before this StudyOur current understanding of Parkinsons disease (PD) is disproportionately based on studying populations of European ancestry, leading to a significant gap in our knowledge about the genetics, clinical characteristics, and pathophysiology in underrepresented populations. This is particularly notable in individuals of African and African admixed ancestries. Over the last two decades, we have witnessed a revolution in the research area of complex genetic diseases. In the PD field, large-scale genome-wide association studies in the European, Asian, and Latin American populations have identified multiple risk loci associated with disease. These include 78 loci and 90 independent signals associated with PD risk in the European population, nine replicated loci and two novel population-specific signals in the Asian population, and a total of 11 novel loci recently nominated through multi-ancestry GWAS efforts. Nevertheless, the African and African admixed populations remain completely unexplored in the context of PD genetics. Added Value of this StudyTo address the lack of diversity in our research field, this study aimed to conduct the first genome-wide assessment of PD genetics in the African and African admixed populations. Here, we identified a genetic risk factor linked to PD etiology, dissected African-specific differences in risk and age at onset, characterized known genetic risk factors, and highlighted the utility of the African and African admixed risk haplotype substructure for future fine-mapping efforts. We identified a novel disease mechanism via expression changes consistent with decreased GBA1 activity levels. Future large scale single cell expression studies should investigate the neuronal populations in which expression differences are most prominent. This novel mechanism may hold promise for future efficient RNA-based therapeutic strategies such as antisense oligonucleotides or short interfering RNAs aimed at preventing and decreasing disease risk. We envisage that these data generated under the umbrella of the Global Parkinsons Genetics Program (GP2) will shed light on the molecular mechanisms involved in the disease process and might pave the way for future clinical trials and therapeutic interventions. This work represents a valuable resource in an underserved population, supporting pioneering research within GP2 and beyond. Deciphering causal and genetic risk factors in all these ancestries will help determine whether interventions, potential targets for disease modifying treatment, and prevention strategies that are being studied in the European populations are relevant to the African and African admixed populations. Implications of all the Available EvidenceWe nominate a novel signal impacting GBA1 as the major genetic risk factor for PD in the African and African admixed populations. The present study could inform future GBA1 clinical trials, improving patient stratification. In this regard, genetic testing can help to design trials likely to provide meaningful and actionable answers. It is our hope that these findings may ultimately have clinical utility for this underrepresented population.
- organism: human
- assay: GWAS summary statistics
- assay: whole-genome sequencing
- assay: long-read sequencing
- research focus: GWAS / association mapping
- research focus: fine-mapping
- research focus: polygenic risk
- research focus: eQTL / molecular QTL mapping
- approach: statistical model
- approach: experimental assay
In-depth summary
Claude Opus 5.5 · xhigh effort · 8 Nov 2026 · from the full text of version 1 (7 May 2023) · CC0 1.0
Intronic GBA1 variant rs3115534-G, common in African and African admixed ancestry but rare elsewhere, raises Parkinson's disease risk (OR 1.58) and lowers onset age.
Parkinson's disease (PD) genetics rests on European, Asian and Latin American cohorts, whereas African and African admixed populations, with high diversity and short linkage disequilibrium blocks, lacked any genome-wide study. This first such GWAS finds a population-specific GBA1 signal acting through a non-coding, likely expression-based route, unlike the coding GBA1 mutations known in Europeans, informing trial stratification and RNA-based therapy ideas.
Key results
Shownrs3115534-G, an intronic GBA1 variant, is the lead genome-wide significant PD risk signal in the African plus African admixed meta-analysis, and each risk allele also brings disease onset earlier.
Meta-analysis, 1,488 cases vs 196,430 controls: OR=1.58, 95% CI 1.37-1.80, P=2.397E-14 (conditional analysis supports one signal). African-only GWAS: OR=1.58, 95% CI 1.35-1.84, P=3.44E-09. Onset: BETA=-2.004, SE=0.57, P=0.0005 in 711 African cases; meta BETA=-3.06, SE=0.40, P=0.008 (about three years earlier per allele). · Abstract; GWAS reveals a novel genome-wide significant…
SuggestedThe risk allele is common in African-ancestry groups but very rare in European, East Asian, South Asian and Ashkenazi samples and tracks African ancestry. It lies on short haplotypes in African groups, which aids fine-mapping, while longer Esan and Yoruba haplotypes suggest an African founder origin.
Allele frequency 0.16 (gnomAD) and 0.21 (1000 Genomes African) vs 0.0015 European, 0.0005 East Asian, 0.0017 South Asian, 0.0009 Ashkenazi (gnomAD). Ancestry regression BETA=0.064, SE=0.024, P=0.01 (a second regression reports BETA=-0.001, P=0.011). Haplotype length 79.19 (90 SNPs) European vs 19.30 (29 SNPs) African and 15.15 (22 SNPs) admixed; lead SNP PP 71.4%. · GWAS reveals a novel genome-wide significant…; Genome-wide comparison of the GBA1 locus across…
ShownShort-read and long-read whole-genome sequencing (WGS) found no coding, structural or splice-altering GBA1 variant that could explain the association. This points to a non-coding mechanism.
Short-read WGS of 206 individuals (141 cases, 65 controls; 39 GG, 69 GT, 98 TT): 96.6% correlation with imputed genotypes and no coding differences by carrier status. Gauchian found no GBA1/GBAP1 rearrangement. ONT long reads (five GG cases, two GT, six TT controls) showed no structural variant. Splice tools predicted no effect. · Short-read Whole Genome Sequencing; Long-read Whole Genome Sequencing; Short- and long-read whole genome sequencing did…
SuggestedIn whole blood, the risk allele is an eQTL that raises GBA1 expression. Yet lymphoblastoid cells from the only available GG carrier showed the lowest glucocerebrosidase (GCase) activity, which the authors propose reflects non-functional transcript diversity.
Mak et al. 2021 eQTL (2,733 samples): P=9.99E-25, BETA=0.238, SE=0.022. GCase activity: GG 762.50 ± 273.50 U vs GT 2743.76 ± 1960.83 U (Welch p=0.00029) and vs TT 1879.94 ± 1010.84 U (p=0.00014). Only one GG donor line was available. · Expression quantitative trait locus analysis…; Discussion
ShownPolygenic risk scores (PRS) built from European data predict PD only modestly in both ancestries, with SNCA rs356182 as the largest contributor. Cases also carry more genome-wide homozygosity, hinting at undiscovered recessive variants.
PRS weighted by Nalls 2019: OR=1.43, 95% CI 1.26-1.61, P=4.37E-05 in African admixed (86 SNPs); OR=1.27, 95% CI 1.16-1.38, P=1.05E-05 in African (79 SNPs); compare OR 3.74 in the European training set. 16 of 91 known variants nominally significant. Runs-of-homozygosity segment count: OR=1.08, 95% CI 1.066-1.11, P=4.63E-16. · Characterization of PD known risk loci and…; Enrichment of homozygosity patterns in African…
Methods, data and reproducibility
Logistic-regression GWAS of imputed genotypes, run separately in African and African admixed groups (covariates: sex, age, ten principal components) and pooled with 23andMe statistics by fixed-effects METAL meta-analysis. The GBA1 hit was followed up with onset regression, fine-mapping, short- and long-read WGS, an eQTL lookup and lymphoblastoid GCase assays.
- African-ancestry individual-level genotypes from IPDGCAN and GP2 (mostly Nigerian) on NeuroBooster and NeuroChip arrays; GP2 data imputed via the TOPMed server · human · n = 997 PD cases; 1,294 controls
- African admixed GWAS: GP2 African-American and Afro-Caribbean genotypes combined with 23andMe African-American summary statistics · human · n = 467 PD cases; 195,120 controls
- 23andMe African admixed GWAS summary statistics, PD diagnosis self-reported · human · n = 288 cases; 193,985 controls
- Short-read whole-genome sequencing (Psomagen) of rs3115534 GG, GT and TT carriers · human · n = 206 (141 cases, 65 controls; 39 GG, 69 GT, 98 TT)
- Oxford Nanopore long-read whole-genome sequencing of DNA from blood or cell lines · human · n = five GG cases, two GT, six TT controls
- Patient-derived lymphoblastoid cell lines from the Coriell repository, used for fluorometric GCase activity assays · human
- Whole-blood eQTL summary statistics (Mak et al., 2021) from mostly African American and Indigenous American donors · human · n = 2,733 samples
- Published European, East Asian and Latin American PD GWAS summary statistics, plus GP2 European and Amerindian/indigenous genotypes for cross-ancestry comparison · human · n = GP2 European 13,186; Amerindian 230 cases, 182 controls
- Methods & tools
- NeuroBooster (v.1.0, Illumina) and NeuroChip genotyping arrays, GenoTools v1.0 (QC and machine-learning ancestry prediction), TOPMed Imputation Server (Minimac), PLINK 1.9 and 2.0, METAL fixed-effects meta-analysis, coloc Approximate Bayes Factor fine-mapping, LDlink, Broad Institute functional-equivalence and joint discovery WGS pipelines, Gauchian v1.0.2, STRipy v2.2 and ANNOVAR, ONT PromethION with Guppy v6.12, winnowmap v2.03, Sniffles2 v2.0.3, 4-methylumbelliferone fluorometric GCase activity assay
- Code
- https://github.com/GP2code/GP2-AFR-AAC-metaGWAShttps://github.com/GP2code/GenoToolshttps://github.com/GP2code
- Software & environment
- GenoTools v1.0; PLINK 1.9 and 2.0; Gauchian v1.0.2; STRipy v2.2; minKNOW v22.10.7; Guppy v6.12; winnowmap v2.03; Sniffles2 v2.0.3; short reads aligned to GRCh38DH; analyses run on Terra.
- Assay conditions & compute
- Analyses ran on the Terra platform and the NIH Biowulf Linux cluster; each ONT library ran 72 hours on a PromethION R9.4.1 flow cell. No CPU, GPU or memory details are given.
- Reproducibility gaps
- GP2 individual-level data are available only by application through AMP-PD. 23andMe summary statistics also require an application and are left out of the GP2 release 5 summary statistics. LCL donor counts per genotype are not stated.
Limitations (6)
- GCase activity comparisons rely on lymphoblastoid lines from only one homozygous risk-allele carrier, so the proposed loss-of-activity mechanism still needs corroboration. · stated by the authors
- The cohort was too small to uncover further susceptibility loci, and onset estimates in the African admixed group should be read cautiously because cases and GG carriers were few. · stated by the authors
- 23andMe data overlap with those used in the Kim et al. 2022 multi-ancestry meta-analysis, which biases the replication estimates for known PD risk loci. · stated by the authors
- No brain QTL data exist for non-European populations, and sparse African RNA sequencing plus GBA1/GBAP1 multi-mapping reads leave the expression-based mechanism unconfirmed. · stated by the authors
- The risk allele increases GBA1 expression in blood, yet activity was lowest in GG and GT lines exceeded TT lines. The activity pattern is therefore not allele-dose consistent, and the transcript-diversity explanation is not tested here. · noted in this summary
- S13 calls the admixture regression a positive association with African ancestry but reports BETA = -0.001, while S14 reports BETA = 0.064. The role of ancestry confounding in the signal needs clearer treatment. · noted in this summary
Read this if you work on Parkinson's disease genetics outside European cohorts, GBA1 biology or trial stratification, or want a worked example of following an ancestry-specific GWAS hit with long-read sequencing and eQTL data.
- Citations
- 167
- FWCI
- 6.79
- Field percentile
- top 1%
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Insights Into Parkinsons Disease Genetics in African Populations: Expanded GWAS Identifies Ancestry-Specific and Cross-Population Risk Loci
N. Okubadejo, O. O. Ojo, … Asli·medRxiv·Mar 3, 2026
- Large-scale genetic characterization of Parkinson's disease in the African and African admixed populations
Fulya Akçimen, Kimberly Paquette, … Kamalini Ghosh Galvelis·BrainPeer-reviewed·Oct 7, 2025
- The Global Landscape of Genetic Variation in Parkinson's disease: Multi-Ancestry Insights into Established Disease Genes and their Translational Relevance
Lara Mariah Lange, Zih‐Hua Fang, … the Global Parkinson’s Genetic Program (GP2)·medRxiv·Jul 8, 2025
- Genome-wide association analyses reveal susceptibility variants linked to Parkinson's disease in the South African population using inferred global and local ancestry
Kathryn Step, Thiago Peixoto Leal, … Soraya Bardien·medRxiv·Aug 1, 2025
- Polygenic Parkinson's disease genetic risk score as risk modifier of parkinsonism in Gaucher disease
Cornelis Blauwendraat, Nahid Tayebi, … Ellen Sidransky·Movement DisordersPeer-reviewed·Mar 3, 2023
- Genetic Modifiers of Age at Onset for Parkinson's Disease in European Populations
the Global Parkinson’s Genetic Program (GP2), Erin R. Voss·medRxiv·Oct 26, 2026