Abstract

Background. LRRK2 mutations have emerged as the most prevalent and potentially treatable determinants of Parkinson’s disease (PD). Peculiar geographic distribution of these mutations has triggered an interest in genotyping PD cohorts of different ethnic backgrounds for LRRK. Objective. Here, we report on the results of LRRK2 screening in the first Central Asian PD cohort. Methods. 246 PD patients were consecutively recruited by movement disorder specialists from four medical centers in Kazakhstan, and clinicodemographic data and genomic DNA from blood were systematically obtained and shipped to the Institute of Neurology University College London together with DNAs from 200 healthy controls. The cohort was genotyped for five LRRK2 mutations (p.Gly2019Ser, p.Arg1441His, p.Tyr1699Cys, p.Ile2020Thr, and p.Asn1437His) and three East Asian disease-associated variants (p.Gly2385Arg, p.Ala419Val, and p.Arg1628Pro) via Kompetitive allele-specific polymerase chain reaction assay analysis. Results. None of the study subjects carried LRRK2 mutations, whereas the following Asian variants were found with insignificant odds ratios (OR): p.Gly2385Arg (1.2%, minor allele frequency (MAF) 0.007, OR 1.25, ), p.Ala419Val (3.7%, MAF 0.02, OR 1.5, ), and p.Arg1628Pro was found only in 1% of controls. p.Gly2385Arg was positive in a big family with PD and tremor, although with incomplete segregation. One early-onset PD subject was homozygous for p.Ala419Val who developed fast progression and severe dyskinesias. p.Ala419Val was associated with early-onset PD. Conclusions. We showed that East Asian LRRK variants could be found in Central Asian populations but their pathogenicity remains to be elucidated in larger PD cohorts.

1. Introduction

The important recent achievement in Parkinson’s disease (PD) research has been the identification of several causative and risk genes with their putative functions. Leucine-rich repeat kinase 2 (LRRK2) gene has been distinguishable from other known PD genes with a number of functional and epidemiologic features. LRRK2 is a highly conserved and widely expressed gene that encodes a unique multifunctional and multidomain protein, named dardarin [1]. Dardarin is a complex chain of 2 527 polypeptides containing two distinct enzymes, namely protein kinase and guanosine triphosphatase (GTPase), as well as multiple protein interaction domains [2]. These domains might interact with each other and other cell signaling proteins, thus playing a putatively key role in cellular function [35].

Intriguing is the fact that almost every LRRK2 domain is susceptible to PD-associated mutations resulting largely in idiopathic PD- (iPD-) like phenotype and pleomorphic neuronal pathology [2, 4]. Amongst the number of known LRRK2 mutations, p.Gly2019Ser mutation has emerged as an important determinant of familial autosomal dominant PD and iPD in North African Arabs, Ashkenazi Jews, and to a lesser degree in European and North American populations [6].

Interestingly, it appears that some LRRK2 mutations and disease-associated variants are specific to particular ethnic groups, most likely due to common founder effects [7]. This is evidently applicable to p.Gly2019Ser mutation, which is common in PD patients from the Western hemisphere and has not yet been reported in the big PD cohorts from East Asia. Similarly, several PD risk variants of LRRK2 including p.Gly2385Arg, p.Ala419Val, and p.Arg1628Pro have only been reported in East Asian populations and have been absent in the Western PD cohorts [8].

PD genetics has been largely unexplored in several world regions, including Central Asia. Here, we investigate 8 LRRK2 mutations and East Asian risk variants in an interesting population residing between Europe and Asia, in the first cohort of PD patients and healthy controls from Kazakhstan.

2. Study Methodology

2.1. Study Subjects

A total of 246 PD patients were consecutively recruited, with no regard to nationality, during 14 months from the National Center for Neurosurgery in Nur-Sultan city, movement disorders clinics in Almaty city, and a regional hospital in Shymkent city in Kazakhstan. The diagnoses of clinically established and clinically probable PD were made on the basis of the agreement between two movement disorder specialists according to the Movement Disorders Society (MDS) PD criteria [9]. Both iPD and familial PD cases with their available first-degree relatives were included in the study. Clinicodemographic characteristics of the cohort were uploaded to University College London (UCL) Research Data Capturing Database (Redcap) online secure database and its summary is given in Table 1. There were 21 patients (8.5%, 21/246) with a family history of PD and/or tremor, and 31 (12.6%, 31/246) patients with the onset of PD before 40 years of age. The male to female ratio was 0.95 : 1. The mean age of PD onset was 55.06 ± 11.15 (range 14–77), and the mean age at the last examination was 61.7 ± 10.3 (range 28–83). Self-reported nationalities in 72.8% of the cohort were Kazakh, 20.7% were Russian, and the remaining nationalities were Uyghur (2.8%), Tatar (2%), Korean (1.3%), and Tajik (0.4%). Genomic DNAs of age- and gender-matched 200 unrelated control subjects were obtained from the research-ready database of neurologically healthy Kazakhs from the National Center for Biotechnology, Nur-Sultan (NCB).

This study was approved by the Research Ethics Committee of NCB (4/29.08.2017) and the Institute of Neurology University College London (IoN UCL) (07/Q0512/26). Written informed consent for participation in the study was obtained from each subject. All personal information was hidden with unique study identifiers.

2.2. Genetic Analysis

Whole venous blood was collected from the subjects at the recruitment centers and sent to NCB for DNA extraction using the standardized laboratory protocols. DNAs were then shipped to IoN UCL for genetic analysis. At IoN, DNAs were checked for quality and concentrations using NanoDrop Spectrophotometer (Thermo Scientific, Waltham, MA, USA). Samples were uniformly diluted to 25 ng/μl, and 25 μl volume of DNA from each sample was transferred to 96-well plates.

On the basis of the literature review, 8 LRRK2 single-nucleotide polymorphisms (SNPs) of interest were selected for the analysis. Five of them are commonly reported LRRK2 mutations: (c.6055G>A) p.Gly2019Ser (rs34637584), (c.4322G>A) p.Arg1441His (rs34995376), (c.5096A>G) p.Tyr1699Cys (rs35801418), (c.6059T>C) p.Ile2020Thr (rs35870237), and (c.4309A>C) p.Asn1437His (rs74163686); and three of them are the East Asian-specific PD-associated variants: (c.7153G>A) p.Gly2385Arg (rs34778348), (c.1256C>T) p.Ala419Val (rs34594498), and (c.4883G>C) p.Arg1628Pro (rs33949390) (Figure 1). The SNPs and surrounding 50 base-pairs were annotated in Ensembl genome browser (see S1 in the supplementary material for comprehensive analysis). The design of primers, adaptation of assays, and genetic analysis were performed in LGC Genomics, England. No positive and negative controls were available for the selected LRRK2 SNPs. LRRK2 genotyping was done using Kompetitive Allele-Specific Polymerase chain reaction (PCR) genotyping assay (KASP™, LGC Genomics. Herts, UK), a method that enables biallelic scoring of SNPs and insertion and deletions at specific loci through competitive allele-specific PCR. Genotyping followed the LGC Genomics protocol. SNP viewer Software (version 1.99, Hoddesdon, UK) was used to visualize the genotyping results (https://www.biosearchtech.com/support/tools/genotyping-software/snpviewer). Familial and young-onset PD cases with positive LRRK2 substitutions were sent for exome sequencing (WES) to Macrogen, South Korea. In addition, when LRRK2 mutations were identified in a proband, Sanger sequencing was performed for all available family members. Thus, we enrolled 13 additional living relatives for the p.Gly2385Arg and p.Ala419Val mutations from two families.

2.3. Statistical Analysis

Statistical analysis was performed using IBM SPSPS version 21 (Chicago, USA). Genotype frequency distributions were tested for conformity to Hardy–Weinberg equilibrium (HWE). T-test and chi-squared tests were used at the level of significance 0.05. Odds ratios (OR) were calculated and presented with 95% confidence interval (CI) values.

3. Results

3.1. KASP Coverage

1.6–3.2% of samples were uncalled in KASP assay analysis and this was within the expected values [10]. The number of uncalled samples for each LRRK2 SNP is shown in Supplementary Materials (see S2 and S4 for comprehensive analysis).

3.2. SPNs with Negative Findings

The following most pathogenic LRRK2 mutations p.Gly2019Ser, p.Arg1441His, p.Tyr1699Cys, p.Ile2020Thr, and p.Asn1437His were not found in our cohort of PD patients and controls. All of the 246 PD subjects and 200 controls were homozygous for wild-type alleles of these SNPs (Table 2). The allelic frequencies for the aforementioned SNPs were in Hardy–Weinberg equilibrium ().

3.3. SNPs with Positive Findings
3.3.1. p.Gly2385Arg (c.7153G>A)

Monoallelic p.Gly2385Arg variant (rs34778348) was found in three PD cases (1.2%, n = 3/239) and two controls (1%, n = 2/199) (Table 3). The first positive case was a young Kazakh male with an autosomal dominant family history of PD. He developed PD at the age of 38 years, and the onset symptoms were depression, anxiety, hyposmia, and left-hand tremor. He had an affected mother, who had developed PD at the age of 55 years and deceased at the age of 58 years. Interestingly, his mother had three female siblings with upper-limb tremor but no signs of bradykinesia (Figure 2). He had fast progression in the disease course and in two years from the disease onset developed difficulty in rising up from a chair, freezing episodes, urinary frequency, and gait abnormalities. His off-stage MDS UPDRS motor scale score was 79 with Hoehn–Yahr stage 3. There was a good response to levodopa. WES in the proband did not reveal any other known PD genes. p.Gly2385Arg variant was tested in six of his healthy relatives and in those with tremor (n = 4) by Sanger sequencing. The variant was present in one of the two unaffected siblings of the proband, unaffected 20-year-old daughter of the proband, and in only one out of the four relatives with tremor (Figure 2).

The second PD case was a 71-year-old Kazakh female with sporadic PD onset at the age of 66 years. She had a mild and slowly progressive disease course. Her off-stage MDS UPDRS motor score was 16 and Hoehn–Yahr stage 1. Dopamine agonists effectively controlled her motor symptoms.

The third PD case was a 65-year-old Kazakh female with sporadic PD with the onset at the age of 61 years. She had also a mild disease course and was not on levodopa. Her off-stage MDS UPDRS motor score was 14 and Hoehn–Yahr stage 2.

p.Gly2385Arg positive healthy controls were 50-year-old and 54-year-old male and female subjects. The allelic and genotypic frequencies for p.Gly2385Arg were in Hardy–Weinberg equilibrium () and did not statistically differ between cases and controls (OR 1.25, 95% C.I.: 0.2071–7.5688, ) (Tables 2 and 4).

3.3.2. p.Ala419Val (c.1256C>T)

LRRK2 p.Ala419Val variant (rs34594498) was positive in 9 PD cases (3.7%, n = 9/242) and 5 controls (2.5%, n = 5/199), giving the OR of 1.5 (95% C.I.: 0.4941–4.5463, ) (Table 4). WES in the probands with young-onset and familial PD did not reveal any other known PD genes. One positive case had homozygous substitution in c.1256C>T (T/T). This was a 52-year-old Kazakh patient with the onset of sporadic PD at the age of 48 years. The patient expressed akinetic-rigid PD with a good response to levodopa but early and severe motor complications. He reached HY stage 3 in four years from the onset of motor symptoms. Two unaffected children of the proband were heterozygous for LRRK2 c.1256C>T substitution, whereas one unaffected sibling of the proband did not have LRRK2 c.1256C>T substitution on Sanger sequencing segregation analysis (Figure 3) (see Figure S3 in the supplementary material for comprehensive image analysis).

Five out of the 9 LRRK2 p.Ala419Val carriers developed PD before the age of 50 years, the youngest manifestation being at the age of 26 years. The mean age at onset for the p.Ala419Val carriers was 48.3 ± 12.6 (range 26–69) and this did not significantly differ from the noncarriers (48.3 ± 12.6 v54.7 ± 12.3, ). The mean age at examination was 57.4 ± 12.9 (range 32–82), and mean disease duration was 9.1 ± 6.1 (range 2–20) (Table 5). The mean HYS score for the positive cases was 2.5 ± 0.5.

Interestingly, self-reported nationalities in three out of the nine positive cases were Russian, one case was half Kazakh and half Russian, and the remaining 5 cases were Kazakhs. Two of the Russian patients had an autosomal dominant history of PD with onset after the age of 50 years. All of the LRRK2 p.Ala419Val-positive cases had a good response to levodopa. The allele frequencies for LRRK2 p.Ala419Val in PD cases deviated from Hardy–Weinberg equilibrium () (Table 2).

3.3.3. p.Arg1628Pro (c.4883G>C)

PD cases were negative for p.Arg1628Pro (rs33949390) variant, whereas two control subjects were found to be positive (1%, n = 2/199). The variant was in Hardy–Weinberg equilibrium () in controls.

4. Discussion

This study screened cohorts of PD patients and controls from Kazakhstan for five LRRK2 mutations and three Asian disease-associated variants. To date, over 80 LRRK2 disease-causing and disease-associated variants have been described in literature since the gene was discovered in 2004. However, only eight of them have been acknowledged as PD-causing mutations including p.Gly2019Ser, p.Arg1441His, p. Arg1441Cys, p.Arg1441Gly, p.Tyr1699Cys, p.Ile2020Thr, p. Asn1437His, and p. Ile2012Thr. All of these mutations affect the catalytic core of the LRRK2 enzyme [11]. Among these mutations, p.Gly2019Ser is the most common followed by the substitution of arginine (Arg) by glycine (Gly), cysteine (Cys), or histidine (His) in the position 1441 of LRRK2 gene [11]. p.Gly2019Ser mutation has a global prevalence of 1% in patients with iPD and around 4% in familial PD. It is noteworthy that p.Gly2019Ser has been predominantly reported in the North African population where it is responsible for 30–42% of familial and 30–34% of sporadic PD cases. Its prevalence shows high figures in Ashkenazi Jews (28% of familial PD and 10% in iPD) and among European as well as North American populations (6% and 3%, respectively) [12]. Conversely, p.Gly2019Ser mutation has not been reported in Asians (<0.1%) [13, 14]. p.Gly2019Ser was reported in Russian PD cohorts (%), possibly in subjects of Ashkenazi Jewish origin [15].

The currently known most deleterious LRRK2 mutations associated with PD were not found in our cohort. The absence of these mutations in our cohort is similar to other Asian studies [6] This is probably due to the founder effects of these mutations, which seem to be specific to Western populations.

4.1. p.Gly2385Arg

LRRK2 p.Gly2385Arg substitution is located in the WD-40 domain, a toroidal beta-propeller structure responsible for protein-protein interactions [16]. There have been 5 independent case-control studies reporting this variant [6, 13, 1618] and 14 studies screening their cohorts for p.Gly2385Arg [14, 19].

Originally, p.Gly2385Arg variant was identified in a small Taiwanese PD family, in a proband and his affected father [20]. Later, a number of studies in Chinese, Taiwanese, Korean, and Japanese populations have reported p.Gly2385Arg variant significantly more among PD patients than controls, with minor allele frequency (MAF) up to 0.4 in patients (Table 6). Thus, p.Gly2385Arg variant was attributed to a risk factor for sporadic and familial PD. The population attributable risk for p.Gly2385Arg in Han-Chinese ethnicity was estimated to be around 4% [18], and the variant probably originated around 4800 years ago from one ancestor [16]. p.Gly2385Arg variant had a tendency for equal distribution across genders and age groups [6]. Regarding the clinical presentation, p.Gly2385Arg carriers expressed typical PD and homozygous cases were not clinically different from heterozygotes and non-carriers [13]. This variant seems not to influence the age of PD onset, as the mean ages of onset in p.Gly2385Arg carriers and non-carriers have not been consistently reported to significantly differ between these groups [6, 13, 17]. Reports on the association of p.Gly2385Arg carrier status with a family history of PD have also been inconsistent [16, 17].

In regards to non-Chinese populations, p.Gly2385Arg was found in 1.2% (2/166) of PD patients and 0.6% (2/306) of controls in Malay/Indian ethnicity from Singapore [8] The frequency of p.Gly2385Arg in cases and controls was not significantly different in this Singaporean study, and considerably lower than the frequency of 8–10% reported in Chinese PD subjects [6, 8]. A recent study on Malaysian PD subjects has found a significant association between p.Gly2385Arg and increased risk of PD [14] The variant was absent in 405 Iranian subjects [22] and positive in only one individual of Northern European origin among 14,002 screened Caucasian subjects [25].

Evidence from functional studies suggests that p.Gly2385Arg substitution leads to the replacement of hydrophobic glycine with the hydrophilic arginine and increases the net positive charge on the 40WD domain of LRRK2. Both LRRK2 proteins with wild-type and p.Gly2385Arg variant localized to the cytoplasm forming aggregates, but the intensity of apoptosis is higher in p.Gly2385Arg variant under oxidative stress conditions [18].

The frequency of p.Gly2385Arg variant in our study (1.2% patients and 1% controls) was almost similar to non-Chinese Singaporean subjects. If MAF for Gly2385Arg in Chinese and Japanese PD populations was between 0.05 and 0.4 (Table 6), MAF in our study was 0.007. This might suggest that p.Gly2385Arg could be found in non-Chinese Asians but in considerably less proportion. Due to a small amount of non-Chinese subjects screened for p.Gly2385Arg and its equal distribution between patients and healthy controls, currently, it is difficult to ascertain the role of this LRRK2 variant in the risk of PD among Central Asian populations. Contemporary data suggest that p.Gly2385Arg could be a risk factor for PD only in selected Asian races.

On the other hand, among three patients positive to p.Gly2385Arg in our study, we had an interesting familial PD case, where proband and his deceased mother had PD, whereas maternal siblings and one maternal cousin of the proband had unilateral asymmetric UL tremor. Although all other known PD genes have been excluded by WES in the proband, p.Gly2385Arg did not completely segregate in the family, being positive in some unaffected family members and negative in some affected (Figure 2). To date, the variant has been shown to segregate with PD in only one small Taiwanese family with affected proband, affected father of the proband, and unaffected sibling. We showed the segregation of p.Gly2385Arg, although incomplete, in a larger family presenting not only with PD but tremor. The incomplete segregation in our family could be due to reduced penetrance or other unknown genetic factors.

4.2. p.Ala419Val

There have been nine reports describing p.Ala419Val variant [26]. The variant resides in the LRRK2 Armadillo domain and is predicted deleterious by online prediction tools with high conservation in the vertebrates [26].

Initially, the association between this variant and PD was described by Ross et al. [23] where p.Ala419Val was tested in 2,338 Asian subjects from Japan, Korea, and Taiwan (1,376 PD cases and 962) in a large-scale multicenter study. The study results revealed the OR of 2.27 (95% CI: 1.35–3.83, ). Several studies before and after the reported association of p.Ala419Val with PD have found either no carriers of this variant in large cohorts of PD patients and controls or insignificant OR (Table 7), thus considering the variant as putatively nonpathogenic population-specific SNP.

The MAF for p.Ala419Val in PD patients has been 0.002–0.018 in Chinese, and 0.026–0.029 in Japanese and Korean populations in studies reporting positive association [23, 31]. Interestingly, the replication studies in the Chinese and Taiwanese ethnicities have yielded inconsistent results. While p.Ala419Val was negative or positive with insignificant OR in some studies, several studies and their meta-analysis reported a significant association between this variant and predominantly early-onset PD (Table 7, [26]). This has been explained by possible natural sampling variation and population heterogeneity [30] on the one hand. On the other hand, Li et al. [26] argued that the discrepancy is likely to result from different mean ages at onset (AAO) of PD patients in these studies. Thus, while the mean AAO in p.Ala419Val-positive reports on Chinese ethnicity was <55 years [23, 26, 29], negative reports on the same population had AAO above 60 years [28, 30, 31, 33]. Provided the fact that p.Ala419Val might have a strong association with early-onset Chinese PD, the likelihood of yielding positive findings could be higher in young-onset PD cohorts.

The MAF for p.Ala419Val in PD patients in our study was 0.02, which is higher than in Chinese and Taiwanese populations and closer to Japanese and Korean (Table 7). Taking into account the young mean AAO in our Kazakhstani PD cohort (55.06 ± 11.15) and reference to Li et al. [26], one might explain the high MAF for p.Ala419Val. Moreover, the AAO in p.Ala419Val-positive PD patients was remarkably young (48.3 ± 12.6) in our study. However, the high frequency (MAF 0.012) of p.Ala419Val in healthy Kazakhstani controls, who also had a mean age of below 55 years, and insignificant OR do not allow us to ascertain the pathogenicity of p.Ala419Val in Kazakhstani PD.

We have found a previously unreported homozygous carrier of p.Ala419Val variant with PD onset before 50 years and relatively aggressive disease course. If this LRRK2 variant is nonpathogenic and not rare in our population, the likelihood of p.Ala419Val homozygous carriers would increase and this might result in the HWE deviation.

p.Ala419Val variant seems to be present not only in Kazakhs but also in self-reported Russian patients with late-onset or familial PD. Considering the reported specificity of p.Ala419Val to Asian populations, we could speculate that these self-reported Russian subjects in our study could have a mixed ethnic background, particularly with Tatars, an Asian population with some Russian phenotypic features. Alternatively, the variant could also be present in the Russian population.

4.3. p.Arg1628Pro

Regarding p.Arg1628Pro, data from a meta-analysis, including 19 studies with a total of 9,927 PD patients and 8,602 controls, suggest that the variant is significantly associated with the risk of PD in East Asian populations [34]. We failed to find this variant in our PD patients but it was present in controls. This might suggest that p.Arg1628Pro could be a common benign polymorphism in Kazakhstani population.

We have to acknowledge the limitations in our study due to the relatively small sample size, the nonhomogeneous ethnic composition of the PD cohort, as only 72.8% of the PD cohort were Kazakhs. In addition, controls were not perfectly matched to cases by age, gender, and ethnicity.

5. Conclusions

The negative findings on common LRRK2 PD causing mutations, and the presence of LRRK2 Asian-specific variants in our PD patients, although at insignificant level compared with controls, suggest that further large-cohort genetic studies are required in Central Asia to ascertain the pathogenicity of LRRK2 Asian-specific variants in the Central Asian PD population.

Data Availability

VCF files from exome sequencing used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors have no conflicts of interest.

Acknowledgments

Cases were collected as part of the SYNaPS Study Group collaboration funded by The Wellcome Trust and Strategic Award (Synaptopathies) Funding (WT093205MA and WT104033AIA). This research was conducted as part of the Queen Square Genomics group at University College London, supported by the National Institute for Health Research University College London Hospitals Biomedical Research Centre. This research was funded by the Medical Research Council (MRC) (MR/S01165X/1, MR/S005021/1, and G0601943).

Supplementary Materials

S1: the 8 LRRK2 SNPs and surrounding 50 base-pairs annotated in Ensembl genome browser S2: calculations of allelic and genotypic frequencies, odds ratios, and Hardy–Weinberg equilibrium S3: chromatogram for Ala419Val homozygous variant with family segregation S4: the results of KASP analysis. (Supplementary Materials)