Targeted next-generation sequencing identifies novel variants in candidate genes for Parkinson’s disease in Black South African and Nigerian patients
BMC Medical Genetics volume 21, Article number: 23 (2020)
The prevalence of Parkinson’s disease (PD) is increasing in sub-Saharan Africa, but little is known about the genetics of PD in these populations. Due to their unique ancestry and diversity, sub-Saharan African populations have the potential to reveal novel insights into the pathobiology of PD. In this study, we aimed to characterise the genetic variation in known and novel PD genes in a group of Black South African and Nigerian patients.
We recruited 33 Black South African and 14 Nigerian PD patients, and screened them for sequence variants in 751 genes using an Ion AmpliSeq™ Neurological Research panel. We used bcftools to filter variants and annovar software for the annotation. Rare variants were prioritised using MetaLR and MetaSVM prediction scores. The effect of a variant on ATP13A2’s protein structure was investigated by molecular modelling.
We identified 14,655 rare variants with a minor allele frequency ≤ 0.01, which included 2448 missense variants. Notably, no common pathogenic mutations were identified in these patients. Also, none of the known PD-associated mutations were found highlighting the need for more studies in African populations. Altogether, 54 rare variants in 42 genes were considered deleterious and were prioritized, based on MetaLR and MetaSVM scores, for follow-up studies. Protein modelling showed that the S1004R variant in ATP13A2 possibly alters the conformation of the protein.
We identified several rare variants predicted to be deleterious in sub-Saharan Africa PD patients; however, further studies are required to determine the biological effects of these variants and their possible role in PD. Studies such as these are important to elucidate the genetic aetiology of this disorder in patients of African ancestry.
Parkinson’s disease (PD) is a debilitating neurodegenerative disorder that impairs patients’ motor skills, and speech coordination. It is one of the leading causes of disability and mortality among neurological disorders globally . The neuropathological hallmark of PD is the progressive loss of predominantly dopaminergic neurons of the substantia nigra pars compacta of the midbrain, which regulate voluntary movement. The diagnosis of this disorder is largely clinical using criteria such as the UK PD Society Brain Bank criteria (UKPDSBBC) to differentiate ageing related symptoms from PD . The pathobiology is yet to be fully elucidated, but environmental and genetic factors have been linked to PD aetiology [3, 4]. PD symptoms usually manifest in the same way in all patients, but the prevalence, incidence and risk factors may vary according to the geographical region . Estimates of PD prevalence in sub-Saharan Africa (SSA) vary widely across previous studies and range from 10 to 235/100,000 in urban populations [6, 7].
Genetics as an etiologic concept in PD has been well-established [8, 9]. Approximately 5–10% of PD patients have a familial form of the disease, which is due to highly penetrant, rare pathogenic mutations . For sporadic forms of this disorder, the genetics is complex as common genetic variants may act in concert with environmental factors [9,10,11]. The genetic discoveries have led to important hypotheses about the mechanisms underlying PD, which include dysfunction of the ubiquitin–proteasome system and mitochondrial dysfunction coupled with oxidative stress .
Most of the studies on the established PD genes or genes associated with PD including SNCA, LRRK2, PRKN, PINK1, PARK7, ATP13A2 and GBA, have been performed in European, North American, North African Arab or Asian populations [9, 13, 14]. In general, limited studies exist on the genetics of PD in the Black African populations . It has been suggested that the variants most commonly associated with PD are rare among South African PD patients [16, 17]. Similarly, a previous genetic study screened for mutations in LRRK2, PRKN and ATXN3 in 57 Nigerian PD patients but did not identify any pathogenic mutations . African populations have a diverse ancestry, and have more private alleles than any other population, suggesting that the genetic aetiology of PD in African populations could be unique .
Next-Generation Sequencing (NGS) provides a way to explore the genetic basis of diseases, and has resulted in the discovery of a large number of disease-associated mutations . In contrast to whole-genome or whole-exome sequencing [21, 22], targeted sequencing panels  focus the analysis on specific genes of interest. The Ion AmpliSeq™ Neurological Research Panel is a commercially available panel designed to screen genes linked to neurological disorders as well as genes involved in brain function. The primary goals of the present study were to use this panel to determine whether a common pathogenic mutation was present, and to characterise the genetic variation in known and novel PD genes, in a group of Black South African and Nigerian PD patients.
The study group consisted of 33 unrelated Black South African PD patients and 14 unrelated Nigerian PD patients. South African patients were primarily recruited at the Neurology Clinic of Tygerberg Academic Hospital, Cape Town, South Africa, and at the Neurology Clinic of Frere Hospital, East London, South Africa. Nigerian PD patients were recruited at the Neurology Clinic, Obafemi Awolowo University Teaching Hospitals Complex, Ile-Ife, Nigeria. All patients were confirmed to have PD by neurologists, based on the UKPDSBBC diagnostic criteria. All patients provided written informed consent to take part in the study and provided peripheral blood samples for genetic studies. This study was approved by the Health Research Ethics Committee of Stellenbosch University (HREC 2002/C059, N16/04/041 and S16/08/151), and the Ethics and Research Committee of Obafemi Awolowo University Teaching Hospitals (ERC/2015/08/15). Demographic information and clinical characteristics of the patients are provided in Table 1 and Additional file 1: Table S1.
Quality control and annotation of targeted next-generation sequencing (tNGS) data
The Ion AmpliSeq™ Neurological Research panel and the Ion AmpliSeq™ Library Kit 2.0 (Thermo Scientific, Waltham, Massachusetts, USA) were used for multiplex PCR amplification of 751 genes (Additional file 2: Table S2). The intronic regions incorporated as part of the exon targets are listed separately in Additional file 3: Table S3. Details on the library construction and next-generation-sequencing protocols are available in Additional file 4.
The flow space calibration, base calling, alignment with the reference genome (GRCh38–hg19), coverage analysis and variant calling were performed using standard parameters in the Ion Torrent Software Suite (ISS) version 5.4.0. Sequenced variants, including insertions and deletions (INDELs), splice site variants, single nucleotide variants (SNVs), multiple nucleotide variants (MNVs), as well as variants in the 3′ untranslated region (UTR3) and 5′ untranslated region (UTR5) were identified. The variant call format (VCF) files produced by the ISS were filtered using bcftools to ensure that:
FILTER = PASS (the ISS internal criteria for setting the PASS filter were met).
QUAL > 100 (the quality score exceeded 100).
FMT/AO ≥ 20 (there were at least 20 reads for the alternative allele).
FMT/DP ≥ 40 (there were at least 40 reads in total).
FMT/SAF ≥ 5 (there were at least 5 reads of the alternative allele in the forward direction).
FMT/SAR ≥ 5 (there were at least 5 reads of the alternative allele in the reverse direction).
These criteria ensured that the observation was made in both directions with at least five reads in either direction, the overall depth was at least 40 and the alternative allele depth was at least 20.
The sequencing data on the 47 samples were of good overall quality. We plotted a graph for the target region coverage using the bam files generated by the Ion Torrent Variant Caller (Additional file 5: Figure S1). This graph showed that 41 samples had at least 80% coverage of the target region at an average read depth of 40X, three samples had a coverage of 78–79%, and another three samples had a coverage of 61–76%.
The VCF files were merged into a single file and processed with the utility annovar (annovar.openbioinformatics.org/) to produce an annotation file for all the variants that passed the criteria above. All available annotations were included. These included conservation scores, allele frequencies and functional predictions (Additional file 6: Table S4). Perl (https://www.perl.org/) was used to extract data and analyses were carried out in R (R Core Team, 2018) . Variants were extracted as homozygous or heterozygous for the alternative allele. The quality scores for all the variants passing the filters were extracted and analysed in R. Using the bam files, bedtools was used to generate statistics on the coverage (depth of sequencing) for the regions in the Ampliseq capture panel as defined by the manufacturer’s bed file. We focused on variants that are rare in control populations as defined by a minor allele frequency (MAF) threshold of 0.01. We created global classifications of variants and generated a summary of variant types, to encode variants as synonymous, missense or frameshift in the variable amino acid class, insertion, deletion or substitution in the variable mutation type, as well as single or multiple base variants in the variable mutated base. We merged the variant summary (whether the variant was observed as a homozygous or heterozygous), with the annotation.
We used the prediction scores MetaLR and MetaSVM for selecting deleterious sequence variants (Additional file 6: Table S4). MetaLR and MetaSVM are themselves ensembles (composite models) of many other scores . Currently, these two have the best performance on curated data sets (training and test) of non-synonymous variants that contain both deleterious (protein-function altering) and benign variants. We therefore used these to prioritize the rare variants in our data. Both metrics were scaled as probabilities [0, 1] with scores close to 1 indicating certainty that the variant is deleterious. We used a score of > 0.8 as a cut-off for including the variant into our list of rare “pathogenic variants” as recommended by Liu et al . We also used 24 other variant scoring algorithms. We plotted the correlation matrix of all 26 scoring algorithm outputs used in our study (Additional file 7: Figure S2). We generated Radar plots (http://www.cmap.polytechnique.fr/~lepennec/R/Radar/RadarAndParallelPlots.html) for each variant likely to be deleterious to demonstrate the correlation among 17 different scoring algorithms. All scores were standardized to 0–1 scale with score 1 (furthest from the centre of the graph) indicating strongest evidence that the variant is deleterious.
Protein structures and modelling for functional prediction
To determine the consequences introduced by potential pathogenic variants on the protein structure, we selected a variant in ATP13A2, for this analysis. We extracted protein information from the Protein Data Bank (PDB) of the Research Collaboratory for Structural Bioinformatics (RCSB) (https://www.rcsb.org/) . The structure of ATP13A2 was modelled by submitting the 1180 amino acid UniProt accession ACQ9NQ11 to the Phyre2 server . Phyre2 is a suite of tools available on the web to predict and analyse protein structure, function and mutations. The predicted ATP13A2 structure conformed well with those of known P-type ATPase cation pumps [28,29,30,31]. This structure was used for additional modelling using the Maestro 11.4 suite of software (Schrödinger Inc., Cambridge, MA). The protein parameterization was carried out using the Maestro Preparation Wizard. The Ca2+ binding sites were identified based on the availability of coordinating glutamate, aspartate, asparagine and glutamine sidechains as found in the other P-type ATPase structures. The two Ca2+ ions were manually docked into the active sites and the structure's energy minimized. The S1004R mutation was generated using in silico mutagenesis with subsequent energy minimization. Based on these structures, the interaction network between the mutated site and Ca2+ was identified (Kenyon et al. unpublished results).
The mean age-at-onset (AAO) of PD in patients varied between the two study groups. It was 48 ± 8 years and 63 ± 13 years for the South African and Nigerian patients, respectively (Table 1). This may be because recruitment in South Africa was predominantly focussed on patients with earlier AAO (< 50 years). Two of the South African PD patients had a positive family history with both having an affected sibling and an affected parent.
Identification of sequence variants
We applied stringent threshold criteria for the filtering and annotation of the variants to exclude low quality variants. Altogether 25,917 sequence variants passed quality control. We then removed all variants with MAF > 0.01 in any of the sequencing databases used as reference databases for the study (Additional file 6: Table S4) and were left with 14,655 rare variants. These rare variants could be classified into 7934 intronic and 5695 exonic variants (Fig. 1; an interactive html-version of the figure is at BMC website). They included 198 UTR5 and 341 UTR3 variants, as well as 32 frameshift, 3175 synonymous and 2448 missense variants. Altogether 14,057 were SNVs and 598 MNVs. There were 261 insertions, 600 deletions and 13,794 substitutions.
In addition, we separately screened 16 PD genes on the panel (SNCA, LRRK2, PRKN, PINK1, PARK7, ATP13A2, EIF4G1, GIGYF2, PLA2G6, FBXO7, VPS35, MAPT, HTRA2, SPG11, GRN and DCTN1) for all sequence variants, and these results are shown in Additional file 8: Table S5).
Pathogenicity prediction of variants
To determine which rare variants are likely to be deleterious and could potentially contribute to the PD pathobiology in the study participants, we used MetaLR and MetaSVM. We focused on identifying rare (MAF ≤ 0.01) or novel (not seen in any of the databases listed in Additional file 6: Table S4) exonic variants predicted to be deleterious. The goal was to minimize the number of false positives by applying stringent filtering criteria. Altogether, 52 heterozygous, one hemizygous and one homozygous exonic rare (MAF ≤ 0.01) missense variants predicted to be deleterious were found in 42 genes (Tables 2 and 3). This included a heterozygous missense variant in one of the known PD genes, ATP13A2 (S1004R) which was validated by Sanger sequencing (data not shown). Radar plots demonstrating pathogenicity scores for each of these 54 rare variants are shown in Additional file 9: Figure S3.
We then analysed separately the two South African individuals (S43_059 and S94_069) with a positive family history of PD (one affected sibling and an affected parent). They each had three heterozygous variants (Table 3). Both individuals had a pyruvate carboxylase (PC) R732G variant (rs112948607). Additionally, individual S43_059 carried galactosylceramidase (GALC) T445S (rs34134328) and TSC complex subunit 2 (TSC2) S1092 L (rs148527903) variants; while S94_069 carried arylsulfatase A (ARSA) N442S (rs6151427) and prickle planar cell polarity protein 2 (PRICKLE2) Q274P (rs564701683) variants. All of the variants had a high certainty of being predicted to be deleterious (pathogenicity score > 0.8) (Table 2). The GALC variant was excluded based on its high MAF in GnomAD African controls (MAF = 0.016) but all of the other variants are rare (MAF < 0.01) (; Table 2) and are therefore potential candidates. Notably, the PC R732G variant that they both share was not found in any of the other patients screened. However, since there was only DNA available for one affected sibling for each of these patients, co-segregation analysis of the variants with disease could not be performed.
We attempted to prioritise one possible pathogenic variant per patient based on MAF (< 0.01), pathogenicity prediction scores (> 0.8) and evidence of prior association of the gene/protein with PD or Parkinsonism (Table 2; Additional file 10: Table S6). In some cases, the MAF of the variant in African controls in GnomAD was ≥0.01, similar to the frequency observed in the patients (Table 2), and those variants were therefore excluded. The prioritised variants are shown in bold and in green font in Table 3. In a few individuals, one variant could not be prioritised over others as more than one variant fulfilled these criteria.
Protein modelling for the S1004R variant in ATP13A2
When the S1004R variant was inserted into the ATP13A2 structure and energy minimized (Kenyon et al. unpublished results), we found that the peptide backbone around the cation binding site was displaced, changing the distance between the cation and a coordinating atom from 3.44 to 2.66 Å (Additional file 11: Figure S4). It could be postulated that the conformational change may alter the efficiency of the pump by interfering with the reaction cycle [28,29,30].
In this study, 47 Black South African and Nigerian PD patients were screened and 54 potentially deleterious sequence variants with MAF ≤ 0.01 in 42 different genes were identified. The 751-gene panel used in the current study contains only 16 of the 34 known PD genes, but it does have six genes (ATP13A2, LRRK2, PARK7, PINK1, PRKN, and SNCA) with strong prior evidence of being involved in PD pathobiology. We identified a rare sequence variant predicted to be deleterious in only one of these genes, ATP13A2. Notably, we did not identify any of the previously reported pathogenic PD mutations catalogued in the PDmutDB database (https://www.molgen.vib-ua.be/PDMutDB/database) in the SSA patients. One possible reason is that, as seen in previous genetic studies on SSA PD patients, common mutations such as LRRK2 G2019S may be a rare cause of PD in these populations .
Protein modelling analysis of the ATP13A2 S1004R variant, which was found in a South African patient (AAO of 39 years), revealed that the variant is potentially functionally important. An interaction between R1004 and the cation binding site was identified suggesting that the variant would interfere with the function of ATP13A2 as a pump of inorganic cations such as metal ions. A previous functional study demonstrated that increased expression of ATP13A2 supresses α-synuclein toxicity in neural cells and that ATP13A2 was likely to act as a Zn2+ pump . Thus, this variant could potentially contribute to PD however, wet-laboratory functional studies are necessary to prove that the variant is indeed pathogenic.
In the two patients with a possible Mendelian inheritance of PD, five heterozygous variants were identified. Co-segregation analysis of the variants with disease in these families was not possible due to a lack of DNA of the family members. However, none of the genes in which these variants were found has been linked to Mendelian forms of PD. In fact, mutations within the PC, PRICKLE2 and TSC2 genes have previously been associated with non-neurodegenerative diseases including diseases involved in energy deficiency, tumour formation and seizures [34,35,36]. Therefore, it is unlikely that mutations in these genes would contribute to a Mendelian inheritance pattern of PD in these patients. ARSA mutations, similar to GBA mutations, have been previously linked to lysosomal storage diseases (LSDs) . Although lysosomal mechanisms are increasingly being shown to be important in PD pathogenesis, the interplay between genetic mutations, lysosomal storage biology and PD is complex and require further elucidation to understand the underlying biology connecting lysosomal storage and PD. However, there is no evidence currently indicating that ARSA mutations cause familial forms of PD.
Limitations of our study include the fact that the sample size was small making it difficult to estimate the actual contribution of genetic factors to PD in the SSA populations. The belief among Black SSA populations that PD is caused by witchcraft and does not have a genetic link  and the notion that it is part of normal ageing, may have contributed to the difficulty in recruiting more patients for the study. Also, 18 of the 34 previously identified PD genes were not on the panel. In addition, annotation of sequence variants in terms of effect on the protein using bioinformatic tools remains problematic. We chose to use the MetaLR and MetaSVM algorithms that currently appear to perform best, but functional studies are needed for validation of these results. Future studies will involve using a custom-panel that captures all of the known PD genes; using whole-exome or whole-genome sequencing; screening of the patients for copy number variations especially in the PINK1 and PRKN genes; and recruitment of a large number of ethnic-matched controls to determine the frequency of prioritized variants in these populations.
Studies are urgently needed to characterise the genetic variation in the known as well as novel PD genes in the understudied SSA populations. Multi-national collaborations across Africa are essential to recruit the large numbers of patients and controls required. The current study provides a starting point to address this need and although it is acknowledged that the sample sizes used here are relatively small, the use of NGS technologies means that the full spectrum of sequence variation in 751 genes has been captured and is available for future studies. We identified several rare variants predicted to be deleterious and they provide new putative candidates for PD but further studies are required to assess their role in PD pathobiology. It is important to include SSA populations in PD genetic studies to ensure that they do not miss out on the potential benefits and opportunities promised by precision medicine .
Availability of data and materials
The bam files of the tNGS data have been deposited to the European Nucleotide Archive (ENA) and can be retrieved with an accession number PRJEB30330 from https://www.ebi.ac.uk/ena/browse/data-retrieval-rest.
Age at onset
ATPase cation transporting 13A2
Basic local alignment search tool
Hidden markov model
Insertion and deletion
Ion Torrent software suite
Leucine-rich repeat kinase 2
Minor allele frequency
Multiple nucleotide variant
Parkinsonism associated deglycase
Parkinson’s disease mutation database
PTEN induced putative kinase 1
Parkin RBR E3 ubiquitin protein ligase
Single nucleotide variant
UK PD society brain bank criteria
3′ untranslated region
5′ untranslated region
Variant call format
GBD 2015 Neurological Disorders Collaborator Group. Global, regional, and national burden of neurological disorders during 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet Neurol. 2017;16:877–97.
Jankovic J. Parkinson’s disease: clinical features and diagnosis. J Neurol Neurosurg Psychiatry. 2008;79:368–76.
Cannon JR, Greenamyre JT. Gene-environment interactions in Parkinson’s disease: specific evidence in humans and mammalian models. Neurobiol Dis. 2013;57. https://doi.org/10.1016/j.nbd.2012.06.025.
Chin-Chan M, Navarro-Yepes J, Quintanilla-Vega B. Environmental pollutants as risk factors for neurodegenerative disorders: Alzheimer and Parkinson diseases. Front Cell Neurosci. 2015;9. https://doi.org/10.3389/fncel.2015.00124.
Pringsheim T, Jette N, Frolkis A, Steeves TDL. The prevalence of Parkinson’s disease: a systematic review and meta-analysis. Mov Disord Off J Mov Disord Soc. 2014;29:1583–90.
Lekoubou A, Echouffo-Tcheugui JB, Kengne AP. Epidemiology of neurodegenerative diseases in sub-Saharan Africa: a systematic review. BMC Public Health. 2014;14:653.
Williams U, Bandmann O, Walker R. Parkinson’s disease in sub-Saharan Africa: a review of epidemiology, genetics and access to care. J Mov Disord. 2018;11:53–64.
Klein C, Westenberger A. Genetics of Parkinson’s disease. Cold Spring Harb Perspect Med. 2012;2. https://doi.org/10.1101/cshperspect.a008888.
Lill CM. Genetics of Parkinson’s disease. Mol Cell Probes. 2016;30:386–96.
Polito L, Greco A, Seripa D. Genetic profile, environmental exposure, and their interaction in Parkinson’s disease. Park Dis. 2016;2016. https://doi.org/10.1155/2016/6465793.
Zschiedrich K, König IR, Brüggemann N, Kock N, Kasten M, Leenders KL, et al. MDR1 variants and risk of Parkinson disease. Association with pesticide exposure? J Neurol. 2009;256:115–20.
Jain S, Wood NW, Healy DG. Molecular genetic pathways in Parkinson’s disease: a review. Clin Sci Lond Engl 1979. 2005;109:355–64.
Clarimón J, Kulisevsky J. Parkinson’s disease: from genetics to clinical practice. Curr Genomics. 2013;14:560–7.
Popescu C. Mechanisms implicated in Parkinson disease from genetic perspective. Med Clin Rev. 2016;2. https://doi.org/10.21767/2471-299X.1000028.
Blanckenberg J, Bardien S, Glanzmann B, Okubadejo NU, Carr JA. The prevalence and genetics of Parkinson’s disease in sub-Saharan Africans. J Neurol Sci. 2013;335:22–5.
Haylett WL, Keyser RJ, du Plessis MC, van der Merwe C, Blanckenberg J, Lombard D, et al. Mutations in the parkin gene are a minor cause of Parkinson’s disease in the south African population. Parkinsonism Relat Disord. 2012;18:89–92.
Blanckenberg J, Ntsapi C, Carr JA, Bardien S. EIF4G1 R1205H and VPS35 D620N mutations are rare in Parkinson’s disease from South Africa. Neurobiol Aging. 2014;35:445.e1–3.
Okubadejo N, Britton A, Crews C, Akinyemi R, Hardy J, Singleton A, et al. Analysis of Nigerians with apparently sporadic Parkinson disease for mutations in LRRK2, PRKN and ATXN3. PLoS One. 2008;3:e3421.
Gurdasani D, Carstensen T, Tekola-Ayele F, Pagani L, Tachmazidou I, Hatzikotoulas K, et al. The African genome variation project shapes medical genetics in Africa. Nature. 2015;517:327.
Olgiati S, Quadri M, Bonifati V. Genetics of movement disorders in the next-generation sequencing era. Mov Disord Off J Mov Disord Soc. 2016;31:458–70.
Farlow JL, Robak LA, Hetrick K, Bowling K, Boerwinkle E, Coban-Akdemir ZH, et al. Whole-exome sequencing in familial Parkinson disease. JAMA Neurol. 2016;73:68–75.
Steele JC, Guella I, Szu-Tu C, Lin MK, Thompson C, Evans DM, et al. Defining neurodegeneration on Guam by targeted genomic sequencing. Ann Neurol. 2015;77:458–68.
Gorostidi A, Martí-Massó JF, Bergareche A, Rodríguez-Oroz MC, López de Munain A, Ruiz-Martínez J. Genetic mutation analysis of Parkinson’s disease patients using multigene next-generation sequencing panels. Mol Diagn Ther. 2016;20:481–91.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018. (https://www.R-project.org/). https://www.blopig.com/blog/2013/07/citing-r-packages-in-your-thesispaperassignments/. Accessed 5 Oct 2018
Liu X, Wu C, Li C, Boerwinkle E. dbNSFP v3.0: a one-stop database of functional predictions and annotations for human nonsynonymous and splice-site SNVs. Hum Mutat. 2016;37:235–41.
Kouranov A, Xie L, de la Cruz J, Chen L, Westbrook J, Bourne PE, et al. The RCSB PDB information portal for structural genomics. Nucleic Acids Res. 2006;34(Database issue):D302–5.
Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJE. The Phyre2 web portal for protein modeling, prediction and analysis. Nat Protoc. 2015;10:845–58.
Toyoshima C. How Ca2+−ATPase pumps ions across the sarcoplasmic reticulum membrane. Biochim Biophys Acta. 2009;1793:941–6.
Toyoshima C, Cornelius F. New crystal structures of PII-type ATPases: excitement continues. Curr Opin Struct Biol. 2013;23:507–14.
Toyoshima C, Mizutani T. Crystal structure of the calcium pump with a bound ATP analogue. Nature. 2004;430:529–35.
van Veen S, Sørensen DM, Holemans T, Holen HW, Palmgren MG, Vangheluwe P. Cellular function and pathological role of ATP13A2 and related P-type transport ATPases in Parkinson’s disease and other neurological disorders. Front Mol Neurosci. 2014;7. https://doi.org/10.3389/fnmol.2014.00048.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes. bioRxiv. 2019:531210. https://doi.org/10.1101/531210.
Kong SMY, Chan BKK, Park J-S, Hill KJ, Aitken JB, Cottle L, et al. Parkinson’s disease-linked human PARK9/ATP13A2 maintains zinc homeostasis and promotes α-Synuclein externalization via exosomes. Hum Mol Genet. 2014;23:2816–33.
Monnot S, Serre V, Chadefaux-Vekemans B, Aupetit J, Romano S, De Lonlay P, et al. Structural insights on pathogenic effects of novel mutations causing pyruvate carboxylase deficiency. Hum Mutat. 2009;30:734–40.
Crino PB, Nathanson KL, Henske EP. The tuberous sclerosis complex. N Engl J Med. 2006;355:1345–56.
Sowers LP, Loo L, Wu Y, Campbell E, Ulrich JD, Wu S, et al. Disruption of the non-canonical Wnt gene PRICKLE2 leads to autism-like behaviors with evidence for hippocampal synaptic dysfunction. Mol Psychiatry. 2013;18:1077–89 Note: Erratum: Molec. Psychiat. 19: 742.
Filocamo M, Morrone A. Lysosomal storage disorders: molecular basis and laboratory testing. Hum Genomics. 2011;5:156–69.
Mokaya J, Gray WK, Carr J. Beliefs, knowledge and attitudes towards Parkinson’s disease among a Xhosa speaking black population in South Africa: a cross-sectional study. Parkinsonism Relat Disord. 2017;41:51–7.
Rotimi CN, Bentley AR, Doumatey AP, Chen G, Shriner D, Adeyemo A. The genomic landscape of African populations in health and disease. Hum Mol Genet. 2017;26:R225–36.
We thank the study participants and the recruitment staff in the Division of Neurology at Tygerberg Academic Hospital, Cape Town, and in the Division of Neurology, Department of Medicine, Faculty of Health Sciences, Walter Sisulu University, East London, in South Africa, as well as the Obafemi Awolowo University Teaching Hospitals Complex, Ile-Ife, in Nigeria. We thank Alexandra I. Beasley for assistance with Additional file 8: Table S5.
This work was supported by a grant (R21NS098862) from the National Institute of Neurological Disorders and Stroke, and Fogarty International Center, the National Institutes of Health, USA, Grant 106052 from the South African National Research Foundation, and a Self-Initiated Research Grant from the South African Medical Research Council to SB. HK was supported by the Faculty of Medicine and Health Sciences, Stellenbosch University, South Africa. JC received support from the South African National Research Foundation. SA was supported by the NRF Innovation Postdoctoral Fellowship from the South African National Research Foundation. GT and DLT were supported by the South African Tuberculosis Bioinformatics Initiative (SATBBI), a Strategic Health Innovation Partnership grant from the South African Medical Research Council and South African Department of Science and Technology. This work was also supported by an NIH grant (P50NS072187 to OAR) for the Mayo Clinic Morris K. Udall Center of Excellence in Parkinson’s Disease Research Lewy Body Dementia Association (LBDA) Research Center of Excellence, and an American Parkinson Disease Association (APDA) Center for Advanced Research grant to OAR.
The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
All patients provided written informed consent to take part in the study. This study was submitted to and approved by the Health Research Ethics Committee of Stellenbosch University in South Africa (Approval reference numbers HREC 2002/C059, N16/04/041 and S16/08/151), as well as the Ethics and Research Committee of Obafemi Awolowo University Teaching Hospitals in Nigeria (Approval reference number ERC/2015/08/15).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Individual clinical data on study participants.
Neurological Research panel information provided by Ion AmpliSeq™ (www.ampliseq.com).
Intron regions covered as part of the exon targets.
Protocols for library construction and Ion Torrent sequencing.
Target region coverage for the 47 samples.
Tools and databases used for the annotation of sequence variants.
Correlation between prediction and conservation scores.
Sequence variants found in the known PD genes.
Radar plots for 54 selected rare variants.
Candidate genes with a link to Parkinson’s disease or Parkinsonism.
Stereoscopic crystal structure models of ATP13A2 showing the position of S1004R.
About this article
Cite this article
Oluwole, O.G., Kuivaniemi, H., Abrahams, S. et al. Targeted next-generation sequencing identifies novel variants in candidate genes for Parkinson’s disease in Black South African and Nigerian patients. BMC Med Genet 21, 23 (2020). https://doi.org/10.1186/s12881-020-0953-1