Skip to main content

ORIGINAL RESEARCH article

Front. Aging Neurosci., 11 November 2022
Sec. Parkinson’s Disease and Aging-related Movement Disorders
Volume 14 - 2022 | https://doi.org/10.3389/fnagi.2022.1040293

Subtyping of early-onset Parkinson’s disease using cluster analysis: A large cohort study

Zhou Zhou1,2† Xiaoxia Zhou1† Yaqin Xiang1 Yuwen Zhao1 Hongxu Pan1 Juan Wu1 Qian Xu1 Yase Chen1 Qiying Sun2,3 Xinyin Wu4 Jianping Zhu5 Xuehong Wu5 Jianhua Li6 Xinxiang Yan1,3 Jifeng Guo1,3,7 Beisha Tang1,3,7 Lifang Lei8* Zhenhua Liu1,3,7* for Parkinson’s Disease Movement Disorders Multicenter Database Collaborative Network in China (PD-MDCNC)
  • 1Department of Neurology, Xiangya Hospital, Central South University, Changsha, China
  • 2Department of Geriatrics, Xiangya Hospital, Central South University, Changsha, China
  • 3National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha, China
  • 4Department of Epidemiology and Health Statistics, Xiangya School of Public Health, Central South University, Changsha, China
  • 5Hunan KeY Health Technology Co., Ltd., Changsha, China
  • 6Hunan Creator Information Technology Co., Ltd., Changsha, China
  • 7Key Laboratory of Hunan Province in Neurodegenerative Disorders, Central South University, Changsha, China
  • 8Department of Neurology, The Third Xiangya Hospital, Central South University, Changsha, China

Background: Increasing evidence suggests that early-onset Parkinson’s disease (EOPD) is heterogeneous in its clinical presentation and progression. Defining subtypes of EOPD is needed to better understand underlying mechanisms, predict disease course, and eventually design more efficient personalized management strategies.

Objective: To identify clinical subtypes of EOPD, assess the clinical characteristics of each EOPD subtype, and compare the progression between EOPD subtypes.

Materials and methods: A total of 1,217 patients were enrolled from a large EOPD cohort of the Parkinson’s Disease & Movement Disorders Multicenter Database and Collaborative Network in China (PD-MDCNC) between January 2017 and September 2021. A comprehensive spectrum of motor and non-motor features were assessed at baseline. Cluster analysis was performed using data on demographics, motor symptoms and signs, and other non-motor manifestations. In 454 out of total patients were reassessed after a mean follow-up time of 1.5 years to compare progression between different subtypes.

Results: Three subtypes were defined: mild motor and non-motor dysfunction/slow progression, intermediate and severe motor and non-motor dysfunction/malignant. Compared to patients with mild subtype, patients with the severe subtype were more likely to have rapid eye movement sleep behavior disorder, wearing-off, and dyskinesia, after adjusting for age and disease duration at baseline, and showed a more rapid progression in Unified Parkinson’s Disease Rating Scale (UPDRS) total score (P = 0.002), UPDRS part II (P = 0.014), and III (P = 0.001) scores, Hoehn and Yahr stage (P = 0.001), and Parkinson’s disease questionnaire-39 item version score (P = 0.012) at prospective follow-up.

Conclusion: We identified three different clinical subtypes (mild, intermediate, and severe) using cluster analysis in a large EOPD cohort for the first time, which is important for tailoring therapy to individuals with EOPD.

Introduction

Parkinson’s disease (PD) is the second most common neurodegenerative movement disorder characterized by typical motor symptoms and many less visible non-motor symptoms (NMSs; Poewe et al., 2017; Bloem et al., 2021). It is characterized pathologically by dopaminergic neuronal loss in the substantia nigra pars compacta and intracellular inclusions containing α-synuclein aggregates (Armstrong and Okun, 2020; Liu et al., 2022; Zhou et al., 2022). PD is a progressive and complex neurological disorder with heterogeneous symptomatology (Jafari et al., 2020; Bloem et al., 2021). Although PD is an age-related disease that typically appears after the age of 65 years, the age of onset for approximately 10% of affected individuals is younger than 50 years, referring to early-onset PD (EOPD; Zhang et al., 2005; Schrag and Schott, 2006; Poewe et al., 2017; Niemann and Jankovic, 2019). Patients with EOPD have relatively high clinical heterogeneity and a longer disease course and typically develop motor fluctuations and dyskinesias earlier, which vary dramatically in its clinical manifestations and prognosis (Pagano et al., 2016; Schirinzi et al., 2020). Thus, EOPD requires more personalized treatment and long-term management.

Several previous studies have used cluster analysis to define clinical PD subtypes based on motor severity, motor complications, some non-motor features, and demographic characteristics (Fereshtehnejad et al., 2015; Lawton et al., 2018; De Pablo-Fernández et al., 2019; Belvisi et al., 2021; Brendel et al., 2021; Mestre et al., 2021). Growing evidences have shown that there are distinct subtypes of PD with diverging trends of progression (Qian and Huang, 2019; Hendricks and Khasawneh, 2021; Mestre et al., 2021). However, all previous cluster analyses were limited to patients with PD, and an EOPD cluster analysis was not available. With some patients with EOPD following a relatively benign course and others progressing rapidly to disability, subtyping EOPD is required. It is essential to perform cluster analysis based on deep phenotyping, followed by prospective validation of subtypes. Defining different subcategories of EOPD is key to better understand its underlying disease mechanisms, predict its disease course, and subsequently design more efficient personalized management strategies (Qian and Huang, 2019).

This study aimed to (1) identify clinical EOPD subtypes using cluster analysis based on a comprehensive baseline dataset, (2) assess the clinical characteristics of each EOPD subtype, and (3) compare disease progression between different EOPD subtypes.

Materials and methods

Participants

Participants were enrolled from a large EOPD cohort of the Parkinson’s Disease & Movement Disorders Multicenter Database and Collaborative Network in China (PD-MDCNC) between January 2017 and September 2021. The clinical diagnosis of PD was confirmed by at least two neurological specialists according to the Movement Disorder Society Clinical Diagnostic Criteria for Parkinson’s Disease (Postuma et al., 2015), including diagnoses of either clinically established or probable PD. The exclusion criteria were as follows: (1) familial history, (2) missing data ≥ 10%, and (3) diagnosis of other causes of parkinsonism on baseline or follow-up assessments. The clinical data of all participants were stored in the PD-MDCNC.1 Written informed consent was obtained from all the participants. This study was approved by the Ethics Committee of Xiangya Hospital and was conducted in accordance with the ethical guidelines of the Declaration of Helsinki.

Clinical assessment

All participants underwent comprehensive and standardized clinical assessments. Before conducting the assessment, all researchers were trained to ensure an equal understanding of the scales used and the methods and phrasing for clinical data collection. Demographic information and clinical characteristics were collected, as described in our previous study (Zhao et al., 2020; Zhou et al., 2022). Clinical examinations of motor symptoms were performed based on the Unified Parkinson’s Disease Rating Scale (UPDRS) and Hoehn and Yahr (H&Y) stages, defining motor subtypes as either tremor dominant, postural instability/gait difficulty or indeterminate group. Motor complications, such as dyskinesia and wearing-off were diagnosed by clinicians, and the severities of dyskinesia and wearing-off were evaluated by UPDRS part IV-A and 9-item End-of-dose Wearing-off Questionnaire (WOQ-9), respectively. Freezing of gait (FOG) was evaluated using New Freezing of Gait Questionnaire (NFOGQ).

In addition to motor symptoms, we evaluated a broad range of NMSs based on the NMSS, the Scale for Outcomes in Parkinson’s Disease for Autonomic Dysfunction (SCOPA-AUT), Mini-Mental State Examination (MMSE), Rapid Eye Movement Sleep Behavior Disorder Questionnaire-Hong Kong (RBDQ-HK), Epworth Sleepiness Scale (ESS), Parkinson’s Disease Sleep Scale (PDSS), Hyposmia Rating Scale (HRS), Functional Constipation Diagnostic Criteria Rome III, and Hamilton Depression Scale (HAMD-17). Quality of life was assessed using the Parkinson’s disease questionnaire-39 item version (PDQ-39). Details regarding the clinical scales were provided in our previous study (Zhou et al., 2022). Patients with illiteracy, primary education, and above junior education were identified as having cognitive impairment when the MMSE scores were below 17, 20, and 24 points, respectively. Hyposmia was defined as a total HRS score less than 22.5. Rapid eye movement sleep behavior disorder (RBD) was defined as a total RBDQ-HK scale score no less than 18. Excessive daytime sleepiness (EDS) was defined as a total ESS score higher than 10. Depression was defined as a total HAMD-17 score higher than 7. The levodopa equivalent daily dose (LEDD) was calculated based on a commonly used method (Tomlinson et al., 2010).

After a mean follow-up period of 1.5 years, the same patients were reassessed on the same variables as the baseline.

Database

Our team established the PD-MDCNC, a comprehensive yet flexible, user friendly, secure, and easily accessible database. And we launched the Chinese Early-Onset Parkinson’s Disease Registry (CEOPDR), which is a large and longitudinal study designed to assess clinical features, genetic architecture, imaging, and biologic markers of EOPD progression in China. All study data will be integrated in the CEOPDR study database through the PD-MDCNC.

Data preprocessing

We conducted standardization to eliminate the influence of various dimensions by scaling the variables to zero mean and unit variance. Correlation analysis was performed between the two variables. When the calculated coefficient was greater than 0.98, only one feature was retained. Finally, the analysis did not exclude any factors. To address high-dimensional and multicollinearity problems, we used principal component analysis (PCA). PCA was performed using the Python software (version 3.6).

Cluster analysis

All data downloaded from the PD-MDCNC database were analyzed using R version 4.1.2.2 Agglomerative hierarchical clustering, K-means clustering, and spectral clustering analyses were synchronously performed. We computed the Calinski–Harabasz score to estimate the optimal clustering methods. Ultimately, agglomerative hierarchical clustering was performed because of the higher Calinski–Harabasz score and better-balanced data distribution (Supplementary Tables 1, 2). Visualization of the final hierarchical cluster solution was performed using Python software (version 3.6) (Supplementary Figure 1). A flowchart of data-driven clustering is shown in Supplementary Figure 2.

Statistical analyses

For the analysis of cross-sectional data, continuous variables were analyzed using one-way analysis of variance or non-parametric tests. Categorical variables were analyzed using the chi-squared test. Comparison of the baseline demography and clinical features between the three statistical clusters was also applied, adjusting for age and disease duration (continuous variables were analyzed by linear regression, and categorical variables were analyzed by logistic regression model).

We used general linear models (GLMs) for a comprehensive longitudinal comparison of the progression of the three subtypes. In each GLM, change of clinical characteristics was defined as the dependent variable. To reduce the regression toward the mean bias, the analysis was adjusted by the follow-up duration and baseline values of the clinical factors (Vickers and Altman, 2001; Fereshtehnejad et al., 2017). Statistical significance was defined as P < 0.05. All data were analyzed using the IBM SPSS Statistics version 23.0 (IBM Corp., Armonk, NY, USA).

Results

Overview

A total of 1,217 patients with EOPD were included in this study. The mean age was 50.54 ± 6.82 years, 53.66% were male patients, and the mean age at onset was 44.12 ± 5.54 years, with an average disease duration of 6.34 ± 5.22 years. The mean UPDRS part I, II, III, and total scores were 2.28 ± 1.98, 10.98 ± 6.56, 25.42 ± 15.73, and 40.58 ± 23.25, respectively. Among the entire study population, RBD and dementia were found in 322 (26.46%) and 103 (8.46%) patients, respectively, at baseline. Supplementary Table 3 summarizes the baseline clinical characteristics of the patients.

Cluster results and baseline characteristics in different clusters

The following 25 variables were included in the final clustering solution: age, age at onset, sex, duration, body mass index, LEDD, UPDRS total score, UPDRS part I–III scores, NMSS score, PDSS score, SCOPA-AUT score, PDQ-39 score, H&Y stage score, motor subtypes, dyskinesia, wearing-off, FOG, RBD, depression, EDS, dementia, hyposmia, and constipation (Supplementary Table 3). As illustrated in Supplementary Figure 1 (visualization of the final hierarchical cluster solution in the EOPD cohort), cluster analysis revealed three distinct clusters of patients with EOPD. Detailed characteristics of the three clusters are listed in Table 1. We also demonstrated the discriminative power of these features using the heatmap shown in Figure 1. We observed evident baseline differences in motor and non-motor manifestations among clusters, with clinically important effect sizes.

TABLE 1
www.frontiersin.org

Table 1. Comparison of the baseline demography and clinical features between the three clusters of EOPD cohort based on hierarchical clustering solution.

FIGURE 1
www.frontiersin.org

Figure 1. Heatmap of the three subtypes at baseline. The figure was depicted according to the mean values or percentages of subtypes. The red color represented a more severe deficit and the blue color referred to a less severe symptom. The darker the figure described, the larger difference were shown among the subtypes. Variables with P value <0.05 are shown. LEDD, Levodopa Equivalent Daily Dose; UPDRS, Unified Parkinson’s disease Rating Scale; H&Y, Hoehn and Yahr; FOG, freezing of gait; NMSS, Non-Motor Symptoms Scale; PDSS, Parkinson’s Disease Sleep Scale; SCOPA-AUT, the Scale for Outcomes in Parkinson’s Disease for Autonomic Dysfunction; PDQ-39, Parkinson’s disease questionnaire-39 item version; RBD, Rapid Eye Movement Sleep Behavior Disorder; EDS, excessive daytime sleepiness.

The first cluster of 533 patients (cluster I, termed mild motor and non-motor dysfunction based on baseline features) was characterized by a low frequency of RBD, wearing-off, and dyskinesia and mild motor and NMSs. Motor symptoms/signs were relatively mild, with the lowest mean UPDRS part II and, III and total scores (P < 0.05). The H&Y stage of the patients was relatively mild. RBD [109 (20.45%) patients], wearing-off [104 (19.51%) patients], and dyskinesia [69 (12.95%) patients] were uncommon. Autonomic symptoms were generally mild (5.75 ± 5.35 in the SCOPA-AUT scale, P = 0.013). Moreover, the average LEDD was lower in cluster I than in clusters II and III (P < 0.001).

At the other extreme, the third cluster of 172 patients (cluster III, termed severe motor and non-motor dysfunction based on baseline features) was characterized by a high frequency of RBD, wearing-off and dyskinesia and more severe motor and NMSs. Motor symptoms/signs were relatively severe, with the highest mean UPDRS part II and, III and total scores (P < 0.05). The H&Y stage score were also relatively worse. RBD [70 (40.70%) patients], wearing-off [81 (47.09%) patients], and dyskinesia [60 (34.88%) patients] were more common in cluster III than in clusters I and II (P < 0.05). Autonomic symptoms were the most severe (8.72 ± 6.95 in the SCOPA-AUT scale, P = 0.013). Moreover, the average LEDD was significantly higher in cluster III than in clusters I and II (P < 0.001).

The patients in cluster II (512 patients with the subtype of intermediate EOPD) had intermediate motor and NMSs between clusters I and III. The UPDRS part II, III, and total scores were intermediate. The H&Y stage of the patients was relatively moderate. RBD (143 [27.93%] patients), wearing-off (148 [28.91%] patients) and dyskinesia (103 [20.12%] patients) were moderately frequent.

Disease progression in different clusters

After a mean duration of 1.5 years, follow-up data were available for 454 patients (Table 2 and Figure 2). Patients in cluster III had a dramatically worse prognosis, with a more rapid progression in the UPDRS total (P = 0.009), UPDRS part II (P = 0.035), UPDRS part III (P = 0.004), H&Y stage (P = 0.006), and PDQ-39 (P = 0.037) scores. The intermediate cluster had a medium progression rate, which was slightly higher than that of cluster I.

TABLE 2
www.frontiersin.org

Table 2. Longitudinal changes in clinical motor and non-motor outcomes in three different clinical phenotypes of EOPD at follow-up.

FIGURE 2
www.frontiersin.org

Figure 2. Longitudinal changes in outcomes of interest in different phenotypes of EOPD at follow-up. In each general linear models, change of clinical characteristics was defined as the dependent variable, and was adjusted by the follow-up duration and baseline value of clinical factors. Mean follow-up duration in the entire population = 1.5 years. UPDRS, Unified Parkinson’s disease Rating Scale; H&Y, Hoehn and Yahr; PDQ-39, Parkinson’s disease questionnaire-39 item version.

Results from the GLM adjusted for baseline values and follow-up duration showed that compared to the cluster I, the cluster III subtypes had significantly greater progression in UPDRS total score (7.44 units more increase in compared to the mild subtype), UPDRS part II score (1.67 points more decline), UPDRS part III score (5.63 points more decline), H&Y stage score (0.31 points more decline), and PDQ-39 score (7.63 points more decline), which demonstrated the worst prognosis of all groups. Similar hierarchical progression was also observed in several non-motor features, namely, NMSS (1.92 points faster in cluster II, 5.00 points faster in cluster III), SCOPA-AUT (0.51 and 0.56 points faster), RBDQ-HK (0.25 and 2.69 points faster), ESS (0.82 and 1.44 points faster), HRS (−0.02 and −0.05 points faster), and HAMD (0.17 and 0.17 points faster) scores. Nevertheless, all showed no statistical significance.

However, the rate of progression was not statistically different between clusters I and II. As illustrated in Figure 2, the faster slope of progression in cluster III was most observed for the UPDRS total, UPDRS part III, and PDQ-39 scores. Based on this prognostic information, we updated the terminology of cluster III to severe motor and non-motor dysfunction/malignant and that of cluster I to mild motor and non-motor dysfunction/slow progression, leaving cluster II terminology unchanged (i.e., intermediate). The data distribution of all clusters and the prevalence of clinical symptoms at baseline and follow-up are shown in Supplementary Figures 3, 4, respectively.

Discussion

To the best of our knowledge, this is the first study to explore the classification of patients with EOPD with a longitudinal large-sample cohort in a comprehensive database on a broad spectrum of motor and non-motor characteristics. Our study found that the most critical determinants of EOPD subtype and prognosis were motor and some NMSs, especially the UPDRS total score, motor complications, RBD, and autonomic dysfunction. Three subtypes were identified: mild motor and non-motor dysfunction/slow progression (cluster I), intermediate (cluster II) and severe motor and non-motor dysfunction/malignant (cluster III). Identification of these subtypes at baseline helps predict prognosis.

Evidence suggests that early-onset Parkinson’s disease is clinically heterogeneous; however, cluster studies of clinical subtypes in patients with EOPD remain scarce. We compared different clustering solutions in our EOPD cohort and selected the best solution to compare the baseline differences and disease progression among the clusters. As expected, baseline differences in motor and NMSs were observed among the clusters. The mild subtype represents patients with EOPD who have mild motor and non-motor manifestations, and RBD and motor complications might be present but are milder than those present in cluster III. These patients had the most favorable disease course, with the least worsening of the UPDRS total score after 1.5 years. In contrast, the severe subtype had a high frequency of RBD, wearing-off, and dyskinesia at baseline. These patients also exhibited more severe motor and autonomic symptoms. This subtype showed the most rapid and malignant progression rate in terms of the UPDRS total, UPDRS part II, UPDRS part III, H&Y stage, and PDQ-39 total scores. Between these two extremes, the intermediate subtype was defined as having intermediate motor and NMSs. This subtype showed a moderate progression rate, which was slightly higher than that in cluster I. Our results further confirmed the clinical heterogeneity of EOPD. Defining EOPD subtypes contributes to a better understanding of the underlying mechanisms of EOPD, predicts the disease course of EOPD, and leads to tailored treatment strategies (Hendricks and Khasawneh, 2021).

A heterogeneous clinical presentation and prognosis of patients with PD is increasingly recognized; in a recent longitudinal study, the presence of RBD, cognitive impairment, and autonomic dysfunction were the best predictors of a diffuse/malignant phenotype of PD (Fereshtehnejad et al., 2015). In our study, patients with the severe subtype also showed a higher propensity for RBD and more worsened autonomic disturbance, in agreement with previous studies (Fereshtehnejad et al., 2015, 2017). Previous study (Fereshtehnejad et al., 2015, 2017; Udow et al., 2016; Jozwiak et al., 2017) also found that these NMSs often coincided, suggesting a common neurobiological factor (Udow et al., 2016). Reduced noradrenergic function in PD was associated with the presence of RBD and autonomic dysfunction and noradrenergic impairment may contribute to the high prevalence of these NMSs in PD (Sommerauer et al., 2018). In addition, compared to patients with mild and intermediate EOPD subtypes, patients with severe EOPD subtype had more severe motor symptoms and higher frequency of motor complications in our study, which are consistent with the results of other similar studies in PD (van Rooden et al., 2011; Belvisi et al., 2021), and severe motor and non-motor dysfunction/malignant subtype may indicate a relatively diffuse neurodegenerative process. Different from the results of other cohort studies, our study further highlights the importance of motor complications as drivers of EOPD subtyping and prognosis.

The mechanism underlying the EOPD subtype differences remains unclear. Studies have indicated that genetic, environmental factors, gene-environment interactions, or all of these factors may play an important role in the pathogenesis of EOPD. There is accumulating evidence indicating that different genotypes of EOPD may lead to distinct clinical characteristics and rate of progression. For example, EOPD with PRKN, PINK1, or DJ-1 mutations are characterized by good response to L-dopa treatment, dystonia and dyskinesia being relatively common, cognitive decline relatively uncommon and a slower deterioration of the disease when compared with idiopathic PD (Kasten et al., 2018; Lin et al., 2019; Zhao et al., 2020; Sun et al., 2021). Patients of EOPD with SNCA mutations show asymmetric onset, good responsiveness for L-dopa in initial time, early motor complications, rapid progression and worse cognitive impairment (Trinh et al., 2018; Chen et al., 2020; Zhao et al., 2020). Potential explanations of subtype differences include different PD-associated protein dysfunctions (Quinn et al., 2020; Mencke et al., 2021; Menozzi and Schapira, 2021), different accumulation rates of pathogenic alpha-synuclein (Lawton et al., 2022), different trajectories of pathology progression (De Pablo-Fernández et al., 2019) or possibly even different compensatory capacities in the neural circuits (Palermo et al., 2020; Wang et al., 2022). Future studies should focus on understanding the underlying disease pathophysiology that drives these different clinical clusters in EOPD and their subsequent progression.

This study has some limitations. Our cluster analysis solution classified EOPD based on clinical presentation only, and additional variables (i.e., genetics, neuroimaging markers, and other biomarkers) may be able to further refine the clusters. Secondly, our study has the heterogeneity of clinical data collection and quality from different investigators. Nevertheless, to retain the consistency and quality of the data as much as possible, all the researchers involved in the study received standardized and unified training. Thirdly, the follow-up data are incomplete, and the reasons of causing this are varied. For example, the follow-up time of some patients with EOPD has not been reached, and the impact of COVID-19 containment measures. Moreover, the prognostic value and applicability of our recommended EOPD subtyping method should be confirmed in other cohorts.

In summary, we developed three clinical EOPD subtypes based on the cluster analysis of a large EOPD cohort. Motor symptoms, motor complications, and NMSs were distinguished among the clusters, some of which were verified in follow-up data. These findings improve our understanding of EOPD heterogeneity. Exploring EOPD subtypes will shed light on precision medicine and develop more effective approaches for clinical trials and treatment strategies for patients with EOPD in the future.

Conclusion

We identified three different clinical subtypes (mild, intermediate, and severe) using cluster analysis in a large EOPD cohort for the first time, which is important for tailoring therapy to individuals with EOPD.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Xiangya Hospital of Central South University. The patients/participants provided their written informed consent to participate in this study.

Author contributions

ZZ and XZ: research project conception, organization, execution, statistical analysis design and execution of the study, patients’ enrollment and follow up, and writing of the first manuscript draft. YX, YZ, HP, and JW: research project execution and patients’ enrollment and follow up. QX, YC, QS, XY, JG, and BT: research project organization, execution, and manuscript review and critique. XYW, JZ, XHW, and JL: statistical analysis design, execution, review and critique of the study, and manuscript review and critique. LL and ZL: research project conception, organization, execution, and manuscript review and critique. All authors contributed to the article and approved the submitted version for publication.

Funding

This study was supported by the Hunan Innovative Province Construction Project (Grant Nos. 2019SK2335 and 2021SK1010), Natural Science Foundation of Hunan Province (Grant Nos. 2022JJ40843 and 2020JJ5929), National Key Research and Development Program of China (Grant Nos. 2021YFC2501204 and 2016YFC306000), National Natural Science Foundation of China (Grant No. 82001359), and Scientific Research Project of Hunan Provincial Health Commission (Grant No. 202203074637).

Acknowledgments

We thank all subjects for their participation, all investigators for their efforts and support in this study. We also thank Editage (www.editage.cn) for English language editing.

Conflict of interest

Authors JZ and XHW were employees of Hunan KeY Health Technology Co., Ltd. Author JL was employee of Hunan Creator Information Technology Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2022.1040293/full#supplementary-material

Footnotes

  1. ^ http://www.pd-mdcnc.com
  2. ^ https://www.r-project.org

References

Armstrong, M. J., and Okun, M. S. (2020). Diagnosis and treatment of parkinson disease: A review. JAMA 323, 548–560. doi: 10.1001/jama.2019.22360

PubMed Abstract | CrossRef Full Text | Google Scholar

Belvisi, D., Fabbrini, A., De Bartolo, M. I., Costanzo, M., Manzo, N., Fabbrini, G., et al. (2021). The pathophysiological correlates of Parkinson’s disease clinical subtypes. Mov. Disord. 36, 370–379. doi: 10.1002/mds.28321

PubMed Abstract | CrossRef Full Text | Google Scholar

Bloem, B. R., Okun, M. S., and Klein, C. (2021). Parkinson’s disease. Lancet 397, 2284–2303. doi: 10.1016/s0140-6736(21)00218-x

CrossRef Full Text | Google Scholar

Brendel, M., Su, C., Hou, Y., Henchcliffe, C., and Wang, F. (2021). Comprehensive subtyping of Parkinson’s disease patients with similarity fusion: A case study with BioFIND data. NPJ Parkinsons Dis. 7:83. doi: 10.1038/s41531-021-00228-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Gu, X., Ou, R., Zhang, L., Hou, Y., Liu, K., et al. (2020). Evaluating the role of SNCA, LRRK2, and GBA in Chinese patients with early-onset Parkinson’s disease. Mov. Disord. 35, 2046–2055. doi: 10.1002/mds.28191

PubMed Abstract | CrossRef Full Text | Google Scholar

De Pablo-Fernández, E., Lees, A. J., Holton, J. L., and Warner, T. T. (2019). Prognosis and neuropathologic correlation of clinical subtypes of Parkinson disease. JAMA Neurol. 76, 470–479. doi: 10.1001/jamaneurol.2018.4377

PubMed Abstract | CrossRef Full Text | Google Scholar

Fereshtehnejad, S. M., Romenets, S. R., Anang, J. B., Latreille, V., Gagnon, J. F., and Postuma, R. B. (2015). New clinical subtypes of parkinson disease and their longitudinal progression: A prospective cohort comparison with other phenotypes. JAMA Neurol. 72, 863–873. doi: 10.1001/jamaneurol.2015.0703

PubMed Abstract | CrossRef Full Text | Google Scholar

Fereshtehnejad, S. M., Zeighami, Y., Dagher, A., and Postuma, R. B. (2017). Clinical criteria for subtyping Parkinson’s disease: Biomarkers and longitudinal progression. Brain 140, 1959–1976. doi: 10.1093/brain/awx118

PubMed Abstract | CrossRef Full Text | Google Scholar

Hendricks, R. M., and Khasawneh, M. T. (2021). A systematic review of Parkinson’s disease cluster analysis research. Aging Dis. 12, 1567–1586. doi: 10.14336/ad.2021.0519

PubMed Abstract | CrossRef Full Text | Google Scholar

Jafari, Z., Kolb, B. E., and Mohajerani, M. H. (2020). Auditory dysfunction in Parkinson’s disease. Mov. Disord. 35, 537–550. doi: 10.1002/mds.28000

PubMed Abstract | CrossRef Full Text | Google Scholar

Jozwiak, N., Postuma, R. B., Montplaisir, J., Latreille, V., Panisset, M., Chouinard, S., et al. (2017). REM sleep behavior disorder and cognitive impairment in Parkinson’s disease. Sleep 40:zsx101. doi: 10.1093/sleep/zsx101

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasten, M., Hartmann, C., Hampf, J., Schaake, S., Westenberger, A., Vollstedt, E. J., et al. (2018). Genotype-phenotype relations for the Parkinson’s disease genes Parkin, PINK1, DJ1: MDSGene systematic review. Mov. Disord. 33, 730–741. doi: 10.1002/mds.27352

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawton, M., Ben-Shlomo, Y., May, M. T., Baig, F., Barber, T. R., Klein, J. C., et al. (2018). Developing and validating Parkinson’s disease subtypes and their motor and cognitive progression. J. Neurol. Neurosurg. Psychiatry 89, 1279–1287. doi: 10.1136/jnnp-2018-318337

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawton, M., Tan, M. M., Ben-Shlomo, Y., Baig, F., Barber, T., Klein, J. C., et al. (2022). Genetics of validated Parkinson’s disease subtypes in the Oxford discovery and tracking Parkinson’s cohorts. J. Neurol. Neurosurg. Psychiatry 93, 952–959. doi: 10.1136/jnnp-2021-327376

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, C. H., Chen, P. L., Tai, C. H., Lin, H. I., Chen, C. S., Chen, M. L., et al. (2019). A clinical and genetic study of early-onset and familial parkinsonism in taiwan: An integrated approach combining gene dosage analysis and next-generation sequencing. Mov. Disord. 34, 506–515. doi: 10.1002/mds.27633

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Yang, N., Dong, J., Tian, W., Chang, L., Ma, J., et al. (2022). Deficiency in endocannabinoid synthase DAGLB contributes to early onset Parkinsonism and murine nigral dopaminergic neuron dysfunction. Nat. Commun. 13:3490. doi: 10.1038/s41467-022-31168-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Mencke, P., Boussaad, I., Romano, C. D., Kitami, T., Linster, C. L., and Krüger, R. (2021). The role of DJ-1 in cellular metabolism and pathophysiological implications for Parkinson’s disease. Cells 10:347. doi: 10.3390/cells10020347

PubMed Abstract | CrossRef Full Text | Google Scholar

Menozzi, E., and Schapira, A. H. V. (2021). Exploring the genotype-phenotype correlation in GBA-Parkinson disease: Clinical aspects, biomarkers, and potential modifiers. Front. Neurol. 12:694764. doi: 10.3389/fneur.2021.694764

PubMed Abstract | CrossRef Full Text | Google Scholar

Mestre, T. A., Fereshtehnejad, S. M., Berg, D., Bohnen, N. I., Dujardin, K., Erro, R., et al. (2021). Parkinson’s disease subtypes: Critical appraisal and recommendations. J. Parkinsons Dis. 11, 395–404. doi: 10.3233/jpd-202472

PubMed Abstract | CrossRef Full Text | Google Scholar

Niemann, N., and Jankovic, J. (2019). Juvenile parkinsonism: Differential diagnosis, genetics, and treatment. Parkinsonism Relat. Disord. 67, 74–89. doi: 10.1016/j.parkreldis.2019.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Pagano, G., Ferrara, N., Brooks, D. J., and Pavese, N. (2016). Age at onset and Parkinson disease phenotype. Neurology 86, 1400–1407. doi: 10.1212/wnl.0000000000002461

PubMed Abstract | CrossRef Full Text | Google Scholar

Palermo, G., Giannoni, S., Frosini, D., Morganti, R., Volterrani, D., Bonuccelli, U., et al. (2020). Dopamine transporter, age, and motor complications in Parkinson’s disease: A clinical and single-photon emission computed tomography study. Mov. Disord. 35, 1028–1036. doi: 10.1002/mds.28008

PubMed Abstract | CrossRef Full Text | Google Scholar

Poewe, W., Seppi, K., Tanner, C. M., Halliday, G. M., Brundin, P., Volkmann, J., et al. (2017). Parkinson disease. Nat. Rev. Dis. Primers 3:17013. doi: 10.1038/nrdp.2017.13

PubMed Abstract | CrossRef Full Text | Google Scholar

Postuma, R. B., Berg, D., Stern, M., Poewe, W., Olanow, C. W., Oertel, W., et al. (2015). MDS clinical diagnostic criteria for Parkinson’s disease. Mov. Disord. 30, 1591–1601. doi: 10.1002/mds.26424

PubMed Abstract | CrossRef Full Text | Google Scholar

Qian, E., and Huang, Y. (2019). Subtyping of Parkinson’s disease–Where are we up to? Aging Dis. 10, 1130–1139. doi: 10.14336/ad.2019.0112

PubMed Abstract | CrossRef Full Text | Google Scholar

Quinn, P. M. J., Moreira, P. I., Ambrósio, A. F., and Alves, C. H. (2020). PINK1/PARKIN signalling in neurodegeneration and neuroinflammation. Acta Neuropathol. Commun. 8:189. doi: 10.1186/s40478-020-01062-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Schirinzi, T., Di Lazzaro, G., Sancesario, G. M., Summa, S., Petrucci, S., Colona, V. L., et al. (2020). Young-onset and late-onset Parkinson’s disease exhibit a different profile of fluid biomarkers and clinical features. Neurobiol. Aging 90, 119–124. doi: 10.1016/j.neurobiolaging.2020.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Schrag, A., and Schott, J. M. (2006). Epidemiological, clinical, and genetic characteristics of early-onset parkinsonism. Lancet Neurol. 5, 355–363. doi: 10.1016/s1474-4422(06)70411-2

CrossRef Full Text | Google Scholar

Sommerauer, M., Fedorova, T. D., Hansen, A. K., Knudsen, K., Otto, M., Jeppesen, J., et al. (2018). Evaluation of the noradrenergic system in Parkinson’s disease: An 11C-MeNER PET and neuromelanin MRI study. Brain 141, 496–504. doi: 10.1093/brain/awx348

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y. M., Yu, H. L., Zhou, X. Y., Xiong, W. X., Luo, S. S., Chen, C., et al. (2021). Disease progression in patients with parkin-related Parkinson’s disease in a longitudinal cohort. Mov. Disord. 36, 442–448. doi: 10.1002/mds.28349

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomlinson, C. L., Stowe, R., Patel, S., Rick, C., Gray, R., and Clarke, C. E. (2010). Systematic review of levodopa dose equivalency reporting in Parkinson’s disease. Mov. Disord. 25, 2649–2653. doi: 10.1002/mds.23429

PubMed Abstract | CrossRef Full Text | Google Scholar

Trinh, J., Zeldenrust, F. M. J., Huang, J., Kasten, M., Schaake, S., Petkovic, S., et al. (2018). Genotype-phenotype relations for the Parkinson’s disease genes SNCA, LRRK2, VPS35: MDSGene systematic review. Mov. Disord. 33, 1857–1870. doi: 10.1002/mds.27527

PubMed Abstract | CrossRef Full Text | Google Scholar

Udow, S. J., Robertson, A. D., MacIntosh, B. J., Espay, A. J., Rowe, J. B., Lang, A. E., et al. (2016). ‘Under pressure’: Is there a link between orthostatic hypotension and cognitive impairment in α-synucleinopathies? J. Neurol. Neurosurg. Psychiatry 87, 1311–1321. doi: 10.1136/jnnp-2016-314123

PubMed Abstract | CrossRef Full Text | Google Scholar

van Rooden, S. M., Colas, F., Martínez-Martín, P., Visser, M., Verbaan, D., Marinus, J., et al. (2011). Clinical subtypes of Parkinson’s disease. Mov. Disord. 26, 51–58. doi: 10.1002/mds.23346

PubMed Abstract | CrossRef Full Text | Google Scholar

Vickers, A. J., and Altman, D. G. (2001). Statistics notes: Analysing controlled trials with baseline and follow up measurements. BMJ 323, 1123–1124. doi: 10.1136/bmj.323.7321.1123

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Wu, P., Brown, P., Zhang, W., Liu, F., Han, Y., et al. (2022). Association of structural measurements of brain reserve with motor progression in patients with Parkinson disease. Neurology. doi: 10.1212/wnl.0000000000200814 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z. X., Roman, G. C., Hong, Z., Wu, C. B., Qu, Q. M., Huang, J. B., et al. (2005). Parkinson’s disease in China: Prevalence in Beijing, Xian, and Shanghai. Lancet 365, 595–597. doi: 10.1016/s0140-6736(05)17909-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Qin, L., Pan, H., Liu, Z., Jiang, L., He, Y., et al. (2020). The role of genetics in Parkinson’s disease: A large cohort study in Chinese mainland population. Brain 143, 2220–2234. doi: 10.1093/brain/awaa167

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Liu, Z., Zhou, X., Xiang, Y., Zhou, Z., Zhao, Y., et al. (2022). The Chinese Parkinson’s disease registry (CPDR): Study design and baseline patient characteristics. Mov. Disord. 37, 1335–1345. doi: 10.1002/mds.29037

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: early-onset Parkinson’s disease, heterogeneous, subtype, PD-MDCNC, cluster analysis

Citation: Zhou Z, Zhou X, Xiang Y, Zhao Y, Pan H, Wu J, Xu Q, Chen Y, Sun Q, Wu X, Zhu J, Wu X, Li J, Yan X, Guo J, Tang B, Lei L and Liu Z (2022) Subtyping of early-onset Parkinson’s disease using cluster analysis: A large cohort study. Front. Aging Neurosci. 14:1040293. doi: 10.3389/fnagi.2022.1040293

Received: 09 September 2022; Accepted: 27 October 2022;
Published: 11 November 2022.

Edited by:

Woon-Man Kung, Chinese Culture University, Taiwan

Reviewed by:

Robert Fekete, New York Medical College, United States
Yafu Yin, Shanghai Jiao Tong University School of Medicine, China

Copyright © 2022 Zhou, Zhou, Xiang, Zhao, Pan, Wu, Xu, Chen, Sun, Wu, Zhu, Wu, Li, Yan, Guo, Tang, Lei and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhenhua Liu, liuzhenhua@csu.edu.cn; Lifang Lei, Lei_lifang@qq.com

These authors have contributed equally to this work and share first authorship

Download