Abstract

Background. Parkinson’s disease (PD) is an age-related neurodegenerative disease affected by multiple genetic and environmental factors. We performed a case-control study on candidate gene to scrutinize whether genetic variants in LRRK2, SNCA, and ITGA8 genes could be associated with sporadic PD in Chinese Han population. Methods. Five single-nucleotide polymorphisms (SNPs) of LRRK2 (rs1491942), SNCA (rs2301134, rs2301135, and rs356221), and ITGA8 (rs7077361) were selected and genotyped among 583 unrelated PD patients and 558 healthy controls. Results. Rs1491942 of LRRK2 gene had a significantly higher genotype frequency () and allelic G/C frequencies () in PD patients than controls. Rs2301135 of SNCA gene also showed an obvious difference in genotype frequency () and allelic G/C frequencies () between PD patients and controls. SNPs rs2301134 and rs356221 of SNCA gene and rs7077361 of ITGA8 gene lacked the significant association with the susceptibility of PD in Chinese Han population. Conclusions. Our study firstly expresses that rs1491942 of LRRK2 and rs2301135 of SNCA gene are substantially associated with sporadic Parkinson’s disease in Chinese Han population.

1. Introduction

Parkinson’s disease (PD), the second most common neurodegenerative disease after Alzheimer’s disease, consists of two major pathological hallmarks: loss of dopaminergic neurons and the presence of Lewy bodies (LB). The classic manifestations of PD are characterized by resting tremor, rigidity, bradykinesia, and impairment of postural reflexes. In addition, some untypical nonmotor features, such as sleep disturbances, mood disorders, autonomic dysfunction, sensory problems, and cognitive impairment, are highly concerned recently. Even when treated with effective therapies, PD is progressive and somehow leads to disability or even mortality.

Increasing evidence supports that complex factors contribute a lot to the susceptibility of PD, which includes genetic and environmental factors [13]. In the past decades, a large number of Genome-Wide Association Studies (GWAS), Candidate Gene Replication Study (CGRS), and subsequent meta-analysis studies have found that a number of potential genes and single-nucleotide polymorphisms (SNPs) associated with PD, including both risk variants and protective variants [4]. In addition, the previous candidate genetic studies provided conclusive evidence showing SNPs in LRRK2, SNCA, and ITGA8 genes significantly impact PD susceptibility and disease characteristics.

Several variations of LRRK2 gene were identified as risk factors for PD. For example, rs34778348 (G2385R, c.7153G>A) and rs33949390 (R1628P, c.4883G>C) were seen to associate with PD in Asian population [5, 6]. Another novel SNP within LRRK2, rs1491942, was found to be responsible for PD in Caucasian populations [7, 8]. However, it was never reported in Chinese Han population before.

SNCA, as the first pathogenic gene identified in PD, encodes α-synuclein, the primary component of LB, the pathological hallmark of PD. From then on, several SNPs of SNCA were highly considered as the genetic risk factors for sporadic PD. Located in the promoter region of SNCA, two SNPs (rs2301134 and rs2301135) with high allele frequency were reported in some studies as PD-related SNPs in European and Taiwanese cohorts [9, 10]. One SNP (rs356221) in the 3′UTR region of SNCA gene showed association with susceptibility to sporadic PD in Japanese and Taiwanese cohorts [10, 11]. All these three SNPs of SNCA (rs2301134, rs2301135, and rs356221) have not been investigated in Han population on the Mainland of China.

While ITGA8 (encoding integrin alpha 8, a type-I transmembrane protein) gene was firstly proved to connect with idiopathic PD in Caucasian population in Simón-Sánchez’s study [4], it was not featured as a PD relevant gene until Lill’s study revealed its potential association with PD [8]. Additional studies are needed to screen the potential pathogenic variants within this gene and assess the potential role of these variants in PD pathogenesis.

There are no study that explores the association of the three genes and their SNPs with Parkinson’s disease in Chinese Han population. Here, we perform the first SNP replication study on previously published SNPs within SNCA (rs356221, rs2301134, and rs2301135), LRRK2 (rs1491942), and ITGA8 (rs7077361) gene in Chinese Han population to explore the ethnic differences and recognize predictive factors for the diagnosis of PD.

2. Methods

2.1. Subjects

This study recruits 1136 cases in the Neurology Department of the First Affiliated Hospital of Xiamen University, which includes 583 Chinese Han sporadic PD patients and 553 matched healthy controls. PD diagnosis coincided well with the diagnostic criteria of UK Parkinson’s Disease Society Brain Bank [12]. Among all the PD patients, the mean age is 65.10 8.90 and the ratio of male to female patients is 320 : 263. The group of controls consists of healthy volunteers from the Medical Center of the First Affiliated Hospital of Xiamen University, the mean age of which is 65.37 9.03, and the ratio of male to female patients is 286 : 267 (Table 1). All subjects are Han population, and the two groups are matched for age, gender, ethnicity, and area of residence. Moreover, this study has gained approval of the local ethics committees, and all patients and controls signed informed consents.

2.2. Genetic Analysis

Venous blood specimens are collected directly from all PD patients and the healthy controls with ethylene diamine tetraacetic acid (EDTA) anticoagulant. Genomic DNA is extracted from the blood samples with the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) under the manufacturer’s instructions and stored at −20°C. The five SNPs are genotyped by Multiplex Snapshot technique (Applied Biosystems by Life Technologies, Foster City, CA, USA). The primers are designed for each SNP locus using Primer 5 and listed in Table 2. Multiplex PCR reactions are performed to amplify target regions containing the selected SNPs. All products are analyzed by the ABI PRISM 3730 DNA Sequence, of which the sequence analyses are conducted by DNA Sequencing Analysis software, GeneMapper4.0. To confirm the results, 10% patients and 10% controls are randomly selected for Sanger sequencing approaches. The concordance rate for replicate approaches was 100%.

2.3. Statistical Analysis

The statistical analyses are processed with SPSS, version 20.0 (IBM, Armonk, NY, USA). The clinical data are expressed as the means standard deviation (SD) for the continuous variables and as numbers (percentage) for the quantitative variables. Student’s t-test is used to compare the age variables between the patients and controls. The gender variables are assessed by the chi-square test. Differences in frequencies of the alleles and genotypes between cases and controls are tested for each SNP through Pearson’s chi-square test and Fisher’s exact test. The criterion for significance is set at based on two sides for all of the tests. The statistical power is calculated by Power and Sample Size Calculations version 3.1.2. The Hardy-Weinberg equilibrium (HWE) is tested by adopting the public statistics web tool (http://ihg.gsf.de/cgi-bin/hw/hwa1.pl). The Haploview program [13] is used for the calculation of linkage disequilibrium (LD) among the three SNPs in SNCA.

3. Results

Genotype and allele frequencies of each SNP of all 1136 subjects (583 patients and 553 healthy controls) are shown in Table 3. Among PD patients, 131 (22%) had an early age of onset (<50 years) and around 77% of the patients were LOPD (≧50 years). The age () and gender () show no statistical difference between the PD patients and the controls in our study. All the subjects were ethnic Hans.

Linkage disequilibrium between the SNCA SNPs rs2301134 and rs356221 is 0.235, while for rs356221 and rs2301135 it is = 0.066 and for rs2301134 and rs2301135 it is = 0.127. This shows weakly correlation with in rs2301134, rs2301135, and rs356221.

Single marker analysis showed a number of significantly statistical associations in our study. Two SNPs from SNCA and LRRK2 genes displayed values < 0.05 prior to correction, both of which were estimated ORs > 1.4. The best value SNPs from both LRRK2 and SNCA genes were further analyzed for age stratification.

Between all PD patients and controls, LRRK2 gene showed a significant difference in the genotype frequency of variant rs1491942 () and allelic G/C frequencies (, , and 95% CI: 1.55–2.30). In the three SNPs of SNCA gene, only variant rs2301135 met the statistics standard in genotype frequencies () and the allelic G/C frequencies (, OR 7.857, and 95% CI: 4.05–15.26). The ITGA8 variant rs7077361 failed to show significant difference in this group ().

In the subgroup of EOPD patients and the controls (age < 50 years), the differences were still obvious in the genotype frequencies of variant rs1491942 of LRRK2 gene () and allelic G/C frequencies (, OR 1.924, and 95% CI: 1.24–2.98). Three SNPs of SNCA gene () and one SNP of ITGA8 gene () failed to show significant difference in this group. Our statistic data of the LOPD patients and controls aged ≧ 50 years also surpassed the significance thresholds. The variant rs1491942 of LRRK2 gene showed a value of 2.538E-07 in genotype frequencies. And the allelic G/C frequencies have a value of 2.459E-08 while the OR is 1.874, and the 95% CI ranged from 1.50 to 2.34. The difference is obvious in the variant rs2301135 of SNCA gene in the genotype frequencies () and allelic G/C frequencies (, OR = 7.846, and 95% CI: 4.03–15.29) in the subgroup. All SNPs in our study met Hardy-Weinberg equilibrium except for the rs2301135 of SNCA gene, shown in Table 3.

4. Discussion

Since the first Genome-Wide Association Study on sporadic Parkinson’s disease was performed in 2005, a new era starts to gain attention in the genetic basis of Parkinson’s disease [14]. Advances in genotyping technology and meta-analysis have allowed researchers to rapidly identify common variants related to PD in different populations [4]. Though some of the previously nominated PD risk genes were firstly reported in familiar Parkinson’s disease (such as SNCA and LRRK2), both of them were successfully replicated in unrelated sporadic PD patients [15]. Additional associated studies and subsequent meta-analysis contributed a lot to identifying the unnoticed variants that can also drive PD risk, ITGA8 as an example [8].

LRRK2 gene was firstly featured as a PD-related gene in Zimprich’s study of families with autosomal-dominant, late-onset Parkinsonism in 2004 [16]. Variants in different domains of LRRK2 have been identified in both familial and sporadic PD in different populations [1719]. Rs1494942 in LRRK2 was previously found associated with PD in US and European series [8, 20]. Our result, being consistent with previous studies, suggests that polymorphism rs1491942 of LRRK2 is a risk loci of sporadic PD, and the variant carriers may share a similar pathomechanism in different populations. As it is reported, LRRK2 variant carriers share similar clinical and pathological features, including a wild range of onset ages, typical Parkinsonism presentation, and sensitivity to L-dopa therapy [21]. In our PD patients, the rs1491942 shows significantly higher frequencies in both EOPD and LOPD subgroups compared with matched controls. These suggest that rs1491942 does not influence the onset age but contributes to the pathogenesis of EOPD and LOPD in a similar way.

SNCA is the first identified causal gene in familial PD [22]. Our findings are partly consistent with previous reports of the association of polymorphisms in SNCA with the susceptibility of PD in US, Norway, and Italian studies [2325]. The linkage disequilibrium of three SNPs of SNCA gene showed that those SNPs are independent. Only one SNP near the promoter region (rs2301135) shows significant differences between PD patients and the controls in Chinese Han population. Concerning the age of onset, rs2301135 is more likely to associate with late-onset PD in our study, while another study in UK suggested that SNCA risk alleles for PD may associate with earlier onset of PD [26]. Given that genotype frequencies and allelic frequencies of rs2301135 of SNCA gene do not follow the Hardy-Weinberg equilibrium in our study, this suggests the possibility of inappropriate population stratification and selection or other confounding factors in our study. Therefore, these results should be interpreted carefully.

ITGA8 expressing in brain mediates cell-cell interactions and regulates neurite outgrowth of sensory and motor neurons [8]. ITGA8 gene was firstly shown associated with PD in Caucasian population in Simón-Sánchez’s study but failed to replicate in other studies of Greece, Irish, and Polish series [4, 19, 20]. ITGA8 variant rs7077361 showed no evidence of relation to PD in our population. In patients and controls, the observed MAFs of the SNP rs7077361 were similar to those reported in the 1000 genomes Southern Han Chinese (CHS) population. However, there is insufficient power to detect the association of rs7077361 with PD in the current sample size. The lack of association of the rs7077361 in Chinese Han population could be ascribed to the limited sample size and the rare existence of this SNP in Chinese population. Large-sample trials from multicenter are required to better understand the contribution of rs7077361 in ITGA8 to PD susceptibility.

In conclusion, the present study provides considerable evidence to support the significant influence of genetic variants on PD risk. It does not replicate the susceptibility of rs2301134 and rs356221 in SNCA and rs7077361 in ITGA8 for PD but confirms that single-nucleotide polymorphisms rs1491942 of LRRK2 and rs2301135 of SNCA gene are susceptible to sporadic PD in Chinese Han population. Certain variants are responsible for the incidence of this disease, while other variants may modify the onset age, which suggests that distinct aspects of PD have a specific genetic architecture. Further studies are required to enrich genetic architecture of PD.

Competing Interests

The authors declare that they have no competing interests.

Authors’ Contributions

Jie Fang and Kehui Yi contributed equally to this work as first authors.

Acknowledgments

The authors would like to thank all the patients and control subjects who participated in this study. The authors would like to show many thanks to Jingjing Yin, a teacher from New Oriental English training school (E-mail: [email protected]), for the substantial English editing. This work was supported by the Fujian Medical Technology Innovation Programs (no. 2014-CXB-33), Xiamen Important Joint Research Project of Major Diseases (no. 3502Z20149028), Science and Technology Program of Xiamen (no. 3502Z20154014), National Natural Science Foundation of China (no. 81400912), and Natural Science Foundation of Fujian Province (no. 2015J01544).