Abstract

Introduction. Genetic factors play an important role in Parkinson’s disease (PD) risk. However, the genetic contribution to progression in Chinese PD patients has rarely been studied. This study investigated genetic associations with progression based on 30 PD risk loci common in a longitudinal cohort of Chinese PD patients and the Parkinson’s Progression Markers Initiative (PPMI) cohort. Methods. PD patients from the true world (TW) Chinese PD longitudinal cohort and the PPMI cohort with demographic information and assessment scales were assessed. A panel containing 30 PD risk single nucleotide polymorphisms was tested. Progression rates of each scale were derived from random-effect slope values of mixed-effects regression models. Progression rates of multiple assessments were combined by using principal component analysis (PCA) to derive scores for composite, motor, and nonmotor progression. The association of genetic polymorphism and separate scales or PCA progression was analysed via linear regression. Results. In the Chinese PD cohort, MAOB rs1799836 was associated with progression based on the Montreal Cognitive Assessment, the top 3 principal components (PCs) of nonmotor PCA and PC1 of the composite PCA. In the PPMI cohort, both MDS-Unified Parkinson’s Disease Rating Scale II and motor PC1 progression were associated with RIT2 rs12456492. The PARK16 haplotype was associated with Geriatric Depression Scale and the State-Trait Anxiety Inventory for Adults progression, and the SNCA haplotype was associated with the Hoehn-Yahr staging progression and motor PC1 progression. Ethnicity-stratified analysis showed that the association between MAOB rs1799836 and PD progression may be specific to Asian or Chinese patients. Conclusion. MAOB rs1799836 was associated with the progression of nonmotor symptoms, especially cognitive impairment, and the composite progression of motor and nonmotor symptoms within our Chinese PD cohort. The RIT2 rs12456492 and SNCA haplotypes were associated with motor function decline, and the PARK16 haplotype was associated with progression in mood in the PPMI cohort.

1. Introduction

Parkinson’s disease (PD) is the second most common neurodegenerative disease in the world and is characterized by motor disabilities, including tremor, bradykinesia, rigidity, and gait disorder, as well as nonmotor symptoms, including hyposmia, REM sleep behaviour disorder (RBD), cognitive impairment, and others. However, the progression of these different clinical features varies among individuals, and the mechanisms underlying this heterogeneity are not well understood.

There are a number of genes that increase the risk for PD in a largely Mendelian fashion, including SNCA, LRRK2, VPS35, Parkin, PINK1, DJ1, and so on [1]. Genome-wide association studies (GWASs) of PD have identified 90 independent loci associated with disease risk, with an overall heritability estimated to be between 22% and 27% [24]. Several studies have explored certain gene mutations, such as those in LRRK2, Parkin, SNCA, GBA, and others, with distinct clinical presentations and motor or cognitive impairment progression [59]. Nevertheless, these studies were either based on cross-sectional PD cohorts or involved limited genetic variants. Although a recent GWAS revealed that APOE4 was associated with cognitive progression and ATP8B2 was associated with motor progression in longitudinal PD cohorts [7], a longitudinal genome-wide survival study revealed that RIMS2, TMEM108, WWOX, GBA, and APOE were involved in cognitive impairment progression in PD patients [10]. However, the genetic contribution to progression in Chinese PD patients has rarely been studied. Thus, the current study investigated genetic associations with progression in a longitudinal PD cohort based on 30 PD risk loci common in Chinese patients. To consider the inconsistent frequency of PD risk loci between different races, we also explored the associations between genetics and PD progression in a Parkinson’s Progression Markers Initiative (PPMI) cohort in which the Caucasian population was the majority.

2. Methods

2.1. Patients and Grouping

The true world longitudinal Chinese PD (TW) cohort included outpatients from November 2014 to September 2017 from the Department of Neurology, Ruijin Hospital, affiliated to Shanghai Jiao Tong University School of Medicine. These patients were enrolled according to the following criteria: (1) a diagnosis of PD according to Movement Disorder Society Clinical Diagnostic Criteria [11]; (2) available clinical assessment scores; (3) multiple visiting records in the Health Information System; and (4) previous genetic testing related to PD. All the diagnoses of PD patients were made by the senior movement disorder specialist. The exclusion criteria were secondary parkinsonism caused by other neurological conditions, such as vascular, inflammatory, toxic, and drug-induced parkinsonism, or parkinsonism plus syndrome including progressive supranuclear palsy, cortico-basal degeneration, and multiple system atrophy. This study was approved by the local ethics committees of Shanghai Ruijin Hospital. Written informed consent was obtained from all TW participants. As indicated in Figure S1, 50 patients were included for further analyses.

For the PPMI cohort, approval was obtained to download and analyze the publicly accessible whole-genome sequence (WGS) and clinical data [12]. Patients whose longitudinal follow-up evaluations were not consistent with a diagnosis of PD were excluded.

2.2. Clinical Assessments

The collection of demographic and clinical data included age, sex, education, and family history of PD. Overlapping scales shared by the two cohorts were collected to perform genotype-phenotype association analyses (Table S1); these scales included Hoehn-Yahr (H-Y) staging, MDS-Unified Parkinson’s Disease Rating Scale (MDS-UPDRS), Montreal Cognitive Assessment (MoCA), Scales for Outcomes in Parkinson’s Disease—Autonomic Dysfunction (SCOPA-AUT), Epworth Sleepiness Scale Score (ESS), Hamilton Anxiety Rating Scale (HAMA) in the TW cohort or State-Trait Anxiety Inventory for Adults (STAI) in the PPMI cohort, Rapid Eye Movement Sleep Behavior Disorder Hong Kong (RBD-HK) in the TW cohort or Rapid Eye Movement Sleep Behavior Disorder Questionnaire (RBD-Q) in the PPMI cohort, and Hamilton Rating Scale for Depression (HAMD)-17 in the TW cohort or Geriatric Depression Scale (GDS) in the PPMI cohort. In the PPMI cohort, we used motor assessments conducted in the “off” medication state.

2.3. Genetic Testing and Haplotype Block Construction

A panel containing 30 PD-related single nucleotide polymorphisms (SNPs) was designed and tested in the TW PD patients. The genotyping methods were previously reported [13]. Peripheral blood samples were collected, and DNA was extracted from leukocytes by using the phenol–chloroform isopropyl alcohol method. Primers were designed by using Primer 5 (version 5.00, PREMIER Biosoft International). After purification of polymerase chain reaction products by both phosphorylase (FastAP, Applied Biosystems) and exonuclease I (EXO I, Applied Biosystems), the extension reaction was performed by using the SNaPshot Multiplex kit from ABI. The extended products were further purified by using phosphorylase (FastAP, Applied Biosystems) and sampled by using ABI3730xl (Applied Biosystems). The results of SNP typing were analyzed by using Genemap 4.0 (Applied Biosystems). Population structure was determined by using EIGENSOFT v6.1.4 [14]. We excluded rs2230288 from our analysis because all the patients in the TW cohort had the G/G genotype at this locus.

For the PPMI cohort, variants of interest were extracted from whole-genome sequencing (WGS) data (project 118, hg38 aligned January 2021 VCFs) using BCFtools [15].

For the construction of haplotype blocks, two steps were followed: (1) a list of variants was pruned based on linkage disequilibrium (LD) among 29 SNPs in the CHB, CHS, and CEU populations of 1 KG, and those pruned variants on the same chromosome were grouped as LD blocks in the PPMI and TW cohorts (Figure S2), and (2) haplotype alleles were constructed through PHASE version 2.1 [16, 17]. Considering previous studies reporting LD of LRRK2, the LRRK2 haplotype was also assigned in the analysis in addition to the two SNPs [6]. Thus, a total of 3 LD blocks were created, namely, the PARK16 haplotype, SNCA haplotype, and LRRK2 haplotype blocks, and the frequencies of the haplotype alleles are listed in Table S2. Finally, 21 SNPs and 3 LD blocks were included to perform the genotype-phenotype association analyses (Figure S3). All the genetic information collected from the PPMI Cohort Genetic Database (N = 1830) and the TW cohort (N = 415) was used to construct haplotypes (Figure S1). The LDlinkR package in R was employed to construct the haplotypes.

2.4. Progression and PCA Progression

Progression rates were derived from random-effect slope values of mixed-effects regression models for two populations [18]. The disease status of PD patients is always assessed by multiple clinical scales, and significant correlations existed between scales in the TW cohort (Table S3). Thus, the collinearity between these scales implied redundancy and made it difficult to clarify genotype-phenotype associations. To further investigate the specific influence of gene loci on the progression of motor and nonmotor PD symptoms, we performed principal component analysis (PCA) on the motor or nonmotor symptom progression rate after zero centring and scaling to have unit variance after a combination of TW and PPMI PD populations [18]. We found that the first principal components (PCs) of motor and nonmotor progression were moderately correlated (r = −0.435, ). Therefore, the third PCA was conducted to combine all the motor and nonmotor measures to create a composite progression score to illustrate the composite progression of PD. PC1 from this cross-domain PCA accounted for 37.61% of the joint variance. The correlations of polymorphisms versus nonmotor PC1-PC3, motor PC1-PC2, and composite PC1 were considered (Table S4), and the contributions of variables to the PCs are illustrated in Table S5 and Figure S4.

2.5. Ethnicity-Stratified Genotype-Phenotype Correlation

The distinct genotype-phenotype correlation might be caused by population race since the genetic background varied between the PPMI and Chinese PD cohorts mentioned above. To illustrate whether there were ethnicity-specific associations, we employed an ethnicity-stratified association analysis. Since the PPMI PD cohort is a multi-ethnic cohort, among which there were 346 Caucasian people (346/377, 92%), 7 Asian people (7/377, 1.8%), 6 black people, and 18 individuals of other races, the TW cohort and the PPMI cohort were merged and then secondarily grouped into four other cohorts: the PPMI Caucasian, PPMI Asian, TW, and All Asian (merging the PPMI Asian and TW cohorts) cohorts. Genotype-phenotype associations were compared across those 4 groups. Considering the limited sample size of the PPMI Asian cohort, both mild () and strict () association settings were used to explore differences and commonalities with pooled raw progression and PCA progression. Those polymorphisms were defined as mono-ethnic specific genetic features, which were merely correlated with raw progression or PCA progression in mono-ethnic patients.

2.6. Statistical Analysis

R (version 4.1.0) with R Studio (version 1.4.1717) was used to perform the statistical analysis. The Factoextra package was used to extract and visualize the PCA results, and other visualization was performed by using ggplot2. A linear mixed model was constructed by using the lme4 package.

All continuous data are presented as the median (IQR) and were compared by using the Mann‒Whitney test. Categorical data are presented as N (%) and were compared by using the chi-square test. A value < 0.05 (two-tailed) was considered statistically significant. The association between the SNPs or haplotypes and the progression of clinical assessments was performed using a linear regression model and partial correlation with adjustment for sex, age onset, education, and disease duration at baseline under an additive minor genetic model. Bonferroni correction was used for multiple comparisons. A value < 0.002381 (two-tailed) was considered statistically significant after Bonferroni correction.

3. Results

Table 1 presents the demographic information of the TW PD and PPMI PD cohorts. The TW cohort included 50 PD patients with an average follow-up of 1.33 years, and the PPMI cohort included 413 PD patients with an average follow-up of 4.35 years. The TW and PPMI PD cohorts were comparable regarding disease duration (Table 1), while the TW patients were significantly younger and had a lower education level than the PPMI patients (Table 1). Despite 30 common variants, genetic background discrepancies were detected in the two cohorts (Figure S5), and the genotype frequency of most variants was significantly different between the two cohorts (Table S6). PCA progression rates were comparable in the two cohorts; thus, we combined the two cohorts to perform PCA.

3.1. Genotype-Phenotype Association in the TW Cohort

A linear mixed-effect model was used to calculate patient-specific disease progression in the two cohorts separately. An adjusted linear regression model with additive genetic assumption was used to assess the association between disease progression rates and gene polymorphisms (21 SNPs and 3 LD blocks). Age at baseline, sex, disease duration in years, and education level were selected as covariates. In the TW cohort, progression measured by the MoCA total score was significantly associated with B-type monoamine oxidation (MAO-B) rs1799836. One additional C allele led to 0.967 (0.495–1.439) more of the corresponding progression rate () (Table 2). LD SNPs were grouped into LD blocks (noted PARK16 haplotype, SNCA haplotype, and LRRK2 haplotype). No multiple-testing significant correlation between LD blocks and progression of separated scales was found.

For PCA-based progression rates, we found that MAO-B rs1799836 was associated with PC1 (), PC2 (), and PC3 () of nonmotor progression PCA and PC1 () of composite progression PCA (Table 3). No multiple-testing significant or mild correlation was detected between LD blocks and PCs of the three types of PCA.

3.2. Genotype-Phenotype Association in the PPMI Cohort

In the PPMI cohort, both progressions measured by MDS-UPDRSII () and PC1 of motor progression PCA () were associated with SYT4.RIT2 rs12456492, and one more copy of the G allele led to 0.286 (0.104–0.468) and −0.333 (−0.536−0.131) increases in the corresponding progression rate, respectively (Tables 2 and 3). The PARK16 haplotype was associated with progression measured by GDS scores and STAI scores, while the SNCA haplotype was associated with progression based on H-Y staging and PC1 of motor progression PCA (Table 4).

3.3. Ethnicity-Stratified Genotype-Phenotype Association

The significant ethnicity-stratified genotype-phenotype association () is listed in Table S6. As indicated in Table 5 and Figure 1, we obtained 4 Asian- and 1 Caucasian-specific strict-progression-correlated polymorphisms (PARK16 haplotype, MAO-B, NUCKS1, and PRKN; SYT4.RIT2), and 1 shared polymorphism among Asians and Caucasians (SNCA haplotype) based on strict setting association ().

Since most Asians were from the TW cohort and all Caucasians were from the PPMI cohort, ethnicity-specific polymorphisms are likely biased by cohort. Thus, we investigated the genotype-phenotype association within the PPMI Asian cohort with a sample size of 7 PD patients. Considering the limited sample size, we set the statistical significance level and found 10 unique genetic polymorphisms in the PPMI Asian cohort. Among them, the association between MAO-B and MoCA again was significant, as was the association between SYT4.RIT2 and H-Y staging (Table 5). In addition, when all the associations in the PPMI Asian cohort were sorted by the negative logarithm of the P value, the former was the 6th significantly associated gene, while the latter was the 11th (Table S7).

4. Discussion

In our study, we explored the association of 30 PD risk loci common in Chinese patients with PD progression assessed by comprehensive scales in a TW PD cohort, which has rarely been reported in previous studies. In our Chinese cohort, we found that MAO-B rs1799836 was associated with PD progression, which may be specific to Asian or Chinese populations. Furthermore, we also performed the analysis of the association between haplotype and PD progression, as most previous studies performed the analysis based on SNPs. In addition to using the traditional method of analyzing the association of SNPs and raw scale progression in PD separately, a new method of analyzing clinical progression in PD was used in our clinical cohorts by combining multiple assessments in a data-driven PCA to derive scores of composite, motor, and nonmotor progression, combining multiple measures to identify latent components that explain the most variability in the data and reflect disease progression more accurately.

Two main findings were observed in our Chinese cohort. First, we found that MAO-B rs1799836 was associated with a change in MoCA scores in our TW cohort. Previous studies found that MAO-B rs1799836 was associated with the risk of PD in Chinese and Asian populations, whereas this association was insignificant in the Caucasian population [1921]. A randomized, double-blind, placebo-controlled prospective study found that rasagiline, a selective monoamine oxidase type-B inhibitor, exerts beneficial effects on attention and executive functions in nondemented but cognitively impaired PD patients [22]. Another study found that Aβ42 levels were also decreased in carriers of the A allele in the MAO-B rs1799836 polymorphism [23]. These studies indicated that the MAO-B gene may be associated with cognition in PD patients. To further improve the phenotypic measure, we used the PCA data reduction method and found that MAO-B rs1799836 was associated with PC1-3 within nonmotor progression PCA and PC1 of composite progression PCA, indicating that MAO-B rs1799836 may be associated with nonmotor progression and composite progression of PD. One possible mechanism is that MAO-B plays a role in degrading norepinephrine (NE) in addition to dopamine (DA), as NE involves many processes, such as cognitive functions, sleep, pain modulation, mood, autonomic functions, blood‒brain barrier permeability, neuroprotection, and immunological mechanisms in the brain, all of which have been impaired in PD and related to nonmotor symptoms [24, 25].

The finding that MAO-B rs1799836 was associated with PD progression in our Chinese cohort was seldom or never reported in previous studies. As the cohorts used in the previous studies included very few Asian patients, MAO-B rs1799836 might be the Chinese- or Asian-specific SNP associated with PD progression. Thus, we explored the associations of the 30 PD risk loci common in Chinese patients and PD progression in the PPIM cohort, which included mostly Caucasian patients, and found that RIT2 rs12456492 was associated with the change in MDS-UPDRS II and dim1 of motor progression PCA, while no significant association was found between MAO-B rs1799836 and PD progression. The protein encoded by the human RIT2 gene binds to the product of human calmodulin 1 (CALM1), and the latter binds to α-synuclein and microtubule-associated protein tau, suggesting that RIT2 may share the same pathway as SNCA and MAPT in the pathogenesis of PD [2628]. The association of RIT2 rs12456492 with PD risk among the Caucasian population has been investigated and found to be significant [29]. The relevant results in Asian populations in different studies have been inconsistent, but a recent meta-analysis of Asian populations demonstrated that RIT2 rs12456492 may increase PD risk in Asian individuals [30], while RIT2 rs12456492 was not associated with disease progression in our Chinese PD cohort. In a previous GWAS meta-analysis, a suggestive association of PD motor subtype with RIT2 rs12456492 was detected, implicating disease progression. In addition to single SNPs, the association of the SNCA haplotype with a change in H-Y staging was found in the PPMI cohort. SNCA codes for the protein α-synuclein, the most important pathologic protein in PD, and is frequently reported to be associated with alpha-synuclein levels and motor function decline [5, 3133]. PARK16 polymorphisms, including NUCKS1 and RAB29, are associated with decreased PD risk in previous Asian and Caucasian studies [34]. However, the association of the PARK16 haplotype and progression in PD has rarely been reported in previous studies.

The inconsistency between our TW cohort and the PPMI cohort may be caused by the different patient ethnicities. Thus, we performed a study within the PPMI Asian population and found that MAO-B rs1799836 was significantly associated with cognitive decline at a level of . However, we failed to find significance after correction, as only 7 patients were included in this analysis. We also assessed only the Caucasian population of the PPMI cohort and found that MAO-B rs1799836 was not associated with PD progression even at a level of . In addition, we performed the same analysis within all Asian patients in the 2 cohorts and found that MAO-B rs1799836 was associated with PD progression, while RIT2 rs12456492 was not after Bonferroni correction. Thus, we suspected that the association between MAO-B rs1799836 and PD progression may be specific in Asian or Chinese populations. However, this analysis was underpowered due to the small sample size, especially in Asia. Thus, a longitudinal study with a large sample size of Asian or Chinese PD patients is needed.

Our study has several limitations. The sample size of our TW study was small. Second, follow-up was limited. The limited follow-up and small sample size may lead to underpower for variants with smaller effects on progression. A third limitation is that we only detected 30 SNPs rather than GWAS. Thus, a larger and longer study of Chinese PD patients is needed.

5. Conclusion

MAO-B rs1799836 was associated with PD progression in nonmotor symptom, particularly cognitive impairment, and composite progression of motor and nonmotor symptoms in a Chinese PD cohort. The RIT2 rs12456492 and SNCA haplotypes were associated with motor function decline, and the PARK16 haplotype was associated with progression in mood in the PPMI cohort. The association between MAO-B rs1799836 and PD progression may be specific to Asian or Chinese patients.

Abbreviations

ESS:Epworth sleepiness scale score
GDS:Geriatric depression scale
GWAS:Genome-wide association study
HAMA:Hamilton anxiety rating scale
HAMD-17:Hamilton rating scale for depression
LD:Linkage disequilibrium
MDS-UPDRS:MDS-unified parkinson’s disease rating scale
MoCA:Montreal cognitive assessment
PC:Principal component
PCA:Principal component analysis
PD:Parkinson’s disease
PPMI:Parkinson’s progression markers initiative
RBD:REM sleep behaviour disorder
RBD-HK:Rapid eye movement sleep behavior disorder Hong Kong
RBD-Q:Rapid eye movement sleep behavior disorder questionnaire
SCOPA-AUT:Scales for outcomes in Parkinson’s disease—autonomic dysfunction
SNP:Single nucleotide polymorphism
STAI:State-trait anxiety inventory for adults.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

CSS was involved in data acquisition and drafting of the manuscript; WLY did data analysis and interpretation and revised the manuscript; LG did data acquisition and data analysis and interpretation; DJJ, HP, and LJ were involved in data acquisition and genetic testing; LY and RK were responsible for conception and design of the study and research project organization; CZL did conception and design of the study and manuscript revision; CSD was concerned with conception and design of the study, diagnosis and recruitment of participants, and manuscript revision. Shi-Shuang Cui and Ling-Yu Wu contributed equally.

Acknowledgments

This work was supported by grants from the National Natural Science Foundation of China (Grant Nos. 81430022, 81771374, and 81971187), Foundation of Shanghai Municipal Science and Technology Major Project (Grant No. 2018SHZDZX05), Shanghai Municipal Education Commission (Grant No. 2017NKX001), and Shanghai Sailing Program (Grant No. 19YF1429900).

Supplementary Materials

Table S1. Clinical assessments. Table S2. Constructed haplotype with PHASE. Table S3. Correlation between scale-wise progression. Table S4. Correlation between raw PD progression and PCA. Table S5. PC variance. Table S6. Genotype distribution of 30 SNPs. Table S7. Diversity of gene-PD-progression correlation among diverse populations. Figure S1. Flowchart of the study. Figure S2. Chrome-wise pairwise LD. A: PARK16 haplotype; B: SNCA haplotype; C: chr6snps; D: chr12snps. Figure S3. Flowchart of SNP preprocessing. Figure S4. Heatmap of the contribution of raw scale progressions to PCs derived from three types of PCA. A: contribution of motor scales of PC from motor PCA; B: contribution of nonmotor scales to PC from nonmotor PCA; C: contribution of all scales to PC from composite PCA. Figure S5. PCA of progression rate and 30 SNPs. A. PCA of progression rate; B. PCA of 30 SNPs. (Supplementary Materials)