Associative role of HLA-DRB1 as a protective factor for susceptibility and progression of Parkinson’s disease: a Chinese cross-sectional and longitudinal study

Background: Previous genome-wide association studies investigating the relationship between the HLA-DRB1 and the risk of Parkinson’s disease (PD) have shown limited racial diversity and have not explored clinical heterogeneity extensively.

Methods: The study consisted of three parts: a case–control study, a cross-sectional study, and a longitudinal cohort study. The case–control study included 477 PD patients and 477 healthy controls to explore the relationship between rs660895 and PD susceptibility. The cross-sectional study utilized baseline data from 429 PD patients to examine the correlation between rs660895 and PD features. The longitudinal study included 388 PD patients who completed a 3-year follow-up to investigate the effects of rs660895 on PD progression.

Results: In the case–control study, HLA-DRB1 rs660895-G allele was associated with a decreased risk of PD in allele model (adjusted OR=0.72, p = 0.003) and dominant model (AG + GG vs. AA: adjusted OR = 0.67, p = 0.003). In the cross-sectional analysis, there was no association between rs660895 and the onset age, motor phenotype, or initial motor symptoms. In the longitudinal analysis, PD patients with the G allele exhibited a slower progression of motor symptoms (MDS-UPDRS-III total score: β = −5.42, p < 0.001, interaction p_{time × genotype} < 0.001) and non-motor symptoms (NMSS score: β = −4.78, p = 0.030, interaction p_{time × genotype} < 0.001).

Conclusion: Our findings support HLA-DRB1 rs660895-G allele is a protective genetic factor for PD risk in Chinese population. Furthermore, we also provide new evidence for the protective effect of rs660895-G allele in PD progression.

1 Introduction

Parkinson’s disease (PD) is a multifactorial neurodegenerative disease characterized by complex pathogenesis. The pathophysiology of PD involves dopaminergic neuron death and α-synuclein accumulation (Lewy bodies) (Bloem et al., 2021). The study conducted by von Bernhardi et al. has revealed that enhanced microglial activation along with neuroinflammatory processes in the brain contributes to the progression of PD (von Bernhardi et al., 2015). Abnormalities in the peripheral immune system may drive onset and progression of PD by facilitating the immune cells infiltration and neuroinflammation in the central nervous system (Tansey and Romero-Ramos, 2019; Yan et al., 2021). Moreover, previous genome-wide association studies (GWAS) have identified numerous potential genetic loci associated with immune-inflammatory mechanisms, which may be strongly linked to susceptibility to PD (Satake et al., 2009; Simón-Sánchez et al., 2009; Nalls et al., 2019). This series of evidence from genetics and immunology further supports the involvement of immune-inflammatory mechanisms in PD (Simón-Sánchez et al., 2009; Herrero et al., 2015).

Human Leukocyte Antigen (HLA) gene, located on the major histocompatibility complex (MHC) on chromosome 6, encodes a group of proteins that bind peptides from foreign or self-antigens, such as α-synuclein. This binding allows T cells to recognize and subsequently coordinate the immune response. Previous genetic studies have shown that several single nucleotide polymorphisms (SNPs) in HLA region were associated with PD risk (Hamza et al., 2010; Nalls et al., 2011; Ahmed et al., 2012; Wissemann et al., 2013; Nalls et al., 2014). Particularly, the latest researches highlighted the strong correlation between PD and HLA-DR and HLA-DQ alleles within the HLA class II region (Naito et al., 2021). These specific HLAs may cause peptide-specific T cell responses, thereby involved in the pathogenesis of immune-inflammatory processes (Yu et al., 2021).

Due to geographical differences, there are extensive variations in HLA allele and haplotype frequencies among populations worldwide (Buhler and Sanchez-Mazas, 2011; Fernandez Vina et al., 2012). Previous studies have reported multiple loci in the HLA-DR region associated with PD. One of these loci is the SNP rs660895-G allele in the HLA-DRB1 region, which has been confirmed by multiple studies in the European population (Ahmed et al., 2012; Wissemann et al., 2013; Chuang et al., 2017). Nevertheless, very few studies on the link between HLA-DR gene polymorphisms and PD have been conducted with Asian population (Foo et al., 2020), and the results are not consistent with previous large-scale European GWASs (Ahmed et al., 2012). Additionally, there is currently a scarcity of clinical studies investigating the relationship between the genetic heterogeneity of the HLA-DRB1 gene and the clinical characteristics of PD, particularly in longitudinal studies. Given that, we focused on SNP rs660895 to verify the association of HLA-DRB1 with risk of PD within the Han Chinese population. Furthermore, we aimed to delve deeper into the impact of this loci on the clinical characteristics and the trajectory of disease progression among PD patients.

2 Methods

2.1 Study design and data extraction

This multicenter, retrospective observational study conducted in China comprised three phases:

a. Case–control study (1st phase): this sub-study assessed the association between HLA-DRB1 rs660895 polymorphism and PD risk in Chinese individuals.

b. Cross-sectional study (2nd phase): this sub-study investigated the association between clinical characteristics (e.g., onset age, motor phenotypes, and initial motor symptoms) and HLA-DRB1 rs660895 genotypes.

c. Longitudinal cohort study (3rd phase): this component examined the associations of HLA-DRB1 rs660895 with changes of both motor and non-motor symptom over 3 years.

This study conformed to the code of Ethics of the World Medical Association Declaration of Helsinki and was approved by the Ethics Committee of Fujian Medical University Union Hospital (approval number: 2023KY177) and the Ethics Committee of Xuanwu Hospital of Capital Medical University (approval number: 2020060). All enrolled participants furnished written informed consent, which encompassed written informed consent for genetic studies.

Data for this study were obtained from an existing prospective observational cohort study, including consecutive individuals with PD from the Movement Disorders Clinic at Fujian Medical University Union Hospital (Fuzhou, China) and the Department of Neurology of Xuanwu Hospital of Capital Medical University (Beijing, China) between February 2017 and July 2023. All patients were diagnosed with idiopathic PD by two movement disorder neurologists and met the International Parkinson and Movement Disorder Society (MDS) Clinical Diagnostic Criteria for the clinical diagnosis of PD (Postuma et al., 2015). To minimize the effects of other interference factors, subjects with active or chronic inflammatory, autoimmune diseases, malignancy, immunodeficiency and any other neurodegenerative disorders were excluded from the study.

In the initial phase (the case–control study), alongside 477 PD patients whom were successfully typed for HLA-DRB1 rs660895 and met the inclusion criteria, we also included 477 healthy individuals. These controls, genotyped for SNP rs660895 and devoid of any family history of PD, were randomly selected from 20 community units in Beijing and 5 in Fuzhou between April 2019 and September 2021. The next phase, a cross-sectional cohort study, focused on 429 of the 477 PD patients who had thorough baseline assessment data. These patients, who had complete baseline records and were tracked over a three-year span with annual clinical check-ups, progressed to the longitudinal cohort study. When considering consistent data availability throughout the follow-up, 388 (representing 90.4%) of the initial 429 PD patients were engaged in this longitudinal phase. All research subjects were of Han nationality. A comprehensive overview of the study’s methodology and patient inclusion is provided in Figure 1.

Figure 1

Figure 1. Study flow chart and analysis plan.

2.2 Clinical assessment

For each PD patient, a comprehensive clinical evaluation was conducted at baseline and during each annual follow-up visit over a period of 3 years. The data collected included medical history, medication usage, and family history. All participants underwent general medical and neurological examinations. Patients with age at onset ≤ 50 were classified as early-onset PD (EOPD), and age at onset ≥50 years as late-onset PD (LOPD) (Schrag and Schott, 2006). The levodopa equivalent daily dose (LEDD) was calculated based on published recommendations (Tomlinson et al., 2010). The Movement Disorders Society Unified Parkinson’s Disease Rating Scale (MDS-UPDRS) was used to assess Parkinson’s disease-related signs and symptoms (Goetz et al., 2008). The stage of the disease was categorized according to the modified Hoehn and Yahr stage (H&Y) (Goetz et al., 2004). The off-state Parkinsonian motor symptom progression was evaluated using Part 3 of the MDS-UPDRS. Motor phenotypes were determined based on the MDS-UPDRS score, classifying patients into tremor-dominant phenotype (TD), postural instability gait disorder phenotype (PIGD), or indeterminate phenotype. The ratio of tremor score to PIGD score was used to define patients with TD (ratio ≥ 1.15), patients with indeterminate phenotype (0.9 < ratio < 1.15), and patients with PIGD (ratio ≤ 0.9) (Stebbins et al., 2013). Additionally, cardinal motor features including tremor (items 20 and 21), rigidity (item 22), bradykinesia (items 23–26 and 31), and postural and gait abnormality (items 27–30) were calculated using the MDS-UPDRS for a more detailed assessment (Eisinger et al., 2020). The severity of overall non-motor dysfunction was assessed using the Non-motor Symptoms Scale (NMSS) (Chaudhuri et al., 2007). Global cognitive function was evaluated using the Montreal Cognitive Assessment (MoCA) (Nasreddine et al., 2005), and rapid eye movement sleep behavior disorder (RBD) was screened using the Rapid Eye Movement Sleep Behavior Disorder Screening Questionnaire (RBDSQ) (Stiasny-Kolster et al., 2007). Olfactory impairment was measured using a 12-item Sniffin’ Sticks test (SS-12) (Çomoğlu et al., 2015).

2.3 DNA extraction and genotyping

Genomic DNA was extracted from peripheral venous blood samples of all subjects using the Wizard Genomic DNA purification Kit (Promega, USA). Genotyping of the target SNP was performed using a matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometer on the MassARRAY^Ⓡ Analyzer 4 platform (Sequenom Inc., San Diego, CA, USA). The target SNP call rate was 98.9%. Individuals with discordant gender information were removed. The genotype frequency distribution in the control group was consistent with Hardy–Weinberg Equilibrium (HWE) (p = 0.152, p > 0.05) (Supplementary Table S1). SNP genotyping and data analysis were conducted by CapitalBio Technology Inc., Beijing, China.

2.4 Statistical methods

The Fisher exact test was utilized to test for Hardy–Weinberg equilibrium (HWE). For the case–control analysis of this study, allelic and genotypic frequencies were analyzed using logistic regression analysis to compare PD patients and controls, with adjustments made for age and gender. Comparison of baseline cross-sectional parameters was performed by Student’s t test, the Mann–Whitney rank sum test, or the Chi-square test, as appropriate.

For the longitudinal cohort analysis, generalized linear mixed models (GLMMs) were utilized to evaluate the association between the SNP and changes in clinical characteristics over time. The fixed effects included genotype, time, and the interaction between genotype and time. Subjects were included as random effects. The models were presented as unadjusted model, minimally adjusted model (adjusted for age and gender), and fully adjusted model (adjusted for age, gender, onset age, disease course, and LEDDs). Furthermore, GLMMs with the variance components structure (including the time × genotype interaction term) were used to test the genotype effect on the trajectory of disease progression. p values less than 0.05 were considered statistically significant.

Statistical analysis was performed using IBM SPSS Statistics 24.0 (SPSS, Inc., Chicago, IL, USA), and figures were created using GraphPad Prism 7.0 (GraphPad Software, San Diego, CA, USA).

3 Results

3.1 Association of HLA-DRB1 rs660895 and PD risk in the case–control study

In this study, 477 PD patients (male: 235, female: 242; mean age: 64.1 ± 9.6 years) and 477 healthy controls (male: 207, female: 270, mean age: 64.0 ± 5.7 years) were included. The demographic and clinical characteristics of the study participants in the case–control study are presented in Table 1. The two groups in the case–control study were similar in terms of age at diagnosis, gender, and education level.

Table 1

Table 1. Demographic and clinical characteristics of participants in the case–control study.

In our study, the frequency of the minor G-allele of rs660895 was 0.255 in healthy controls, closed to the 0.240 MAF in 1000 Genomes of East Asian populations (Auton et al., 2015). As shown in Table 2, the frequency of the G allele was significantly lower in individuals with PD compared to the control group (crude OR = 0.72, p = 0.003), indicating that the rs660895-G allele was strongly associated with a reduced risk of PD. This association remained statistically significant even after adjusting for age and gender (adjusted OR = 0.72, p = 0.003). Furthermore, we assessed the association between rs660895 and PD using three different genetic models (co-dominant, dominant and recessive). The results showed significant differences between cases and controls in both the co-dominant model (AG vs. AA: crude OR = 0.67, p = 0.005) and dominant model (AG + GG vs. AA: crude OR = 0.67, p = 0.003). After adjusting for age and gender, the association remained significant in the co-dominant model (AG vs. AA: adjusted OR = 0.68, p = 0.007) and the dominant model (AG vs. AA: adjusted OR = 0.67, p = 0.003). These findings suggest the HLA-DRB1 rs660895-G allele is associated with the lower risk of PD in the Han Chinese population.

Table 2

Table 2. Allele and genotype frequencies in HLA-DRB1 rs660895 in cases and controls.

In addition, considering that epigenetic modifications are influenced by geographical traits, we conducted a comparative analysis of samples derived from Fuzhou and Beijing populations. The analysis revealed no significant differences in demographics and genotype composition between two groups (Supplementary Table S2). And after adjusting the model based on different centers, consistent protection against the risk of PD was observed with the rs660895-G allele (Supplementary Table S3 for detail).

3.2 Association of HLA-DRB1 rs660895 and clinical profiles of PD in the cross-sectional study

Under the dominant model of rs660895, the baseline cross-sectional study included 429 PD patients who were divided into two HLA-DRB1 genotype groups: the G allele carriers (AG + GG) and the G allele non-carriers (AA). No significant difference in the demographic characteristics between groups. Additional analysis of PD motor characteristics revealed no significant association between rs660895-G allele carrier status and age of onset (p = 0.352), motor phenotype (p = 0.743), or initial motor symptoms (p = 0.293). The detailed data be shown in the Table 3.

Table 3

Table 3. Demographic and clinical characteristics of PD patients under the dominant model for HLA-DRB1 rs660895 in the cross-sectional cohort study (n = 429).

3.3 Association of HLA-DRB1 rs660895 and disease progression of PD in the longitudinal cohort study

We further explored the relationship between the rs660895-G allele and disease progression in patients with PD. Under the dominant model, the study included a total of 388 patients with complete data. This comprised of 152 G allele carriers (AG + GG) and 236 G allele non-carriers (AA). There were no significant differences between the two trajectories in demographic and clinical variables at baseline. The retention rate was also similar between G allele carriers (89.9%) and G allele non-carriers (90.8%) (p = 0.776). The participants lost to follow-up had similar baseline characteristics to those who completed the study. Supplementary Table S4 provides a summary of the demographic and clinical parameters of the longitudinal analyzed population at baseline in detail.

3.3.1 Effect of HLA-DRB1 rs660895 on the progression of motor symptoms in PD patients

The results of our study indicated that SNP rs660895-G allele showed a strong influence on the progression of motor symptoms over time. As shown in the three-year trajectories of motor symptoms progression (Figures 2A–E), patients with the G allele had lower MDS-UPDRS-III total score, rigidity subscore, bradykinesia subscore, and axial subscore at the 3rd year follow-up in the unadjusted model (MDS-UPDRS-III total score: p < 0.001; rigidity subscore: p = 0.006; bradykinesia subscore: p = 0.001; axial subscore: p < 0.001), except for tremor subscore (p = 0.090). These associations remained significant even in the minimally adjusted model and the fully adjusted model (Table 4). We further assessed the main effect of time, genotype, and the interaction of time and genotype, showing that the effect of genotype on overall motor symptoms, rigidity, bradykinesia, and axial impairment significantly varied by the time of assessment. Longitudinal analysis showed genotype effect with a significant genotype effect and genotype-by-time interaction on overall motor symptoms and axial impairment (MDS-UPDRS-III total scores: p _genotype = 0.032 and p _{time×genotype} < 0.001 in model 2, p _genotype = 0.047 and p _{time×genotype} < 0.001 in model 3; Axial subscore: p _genotype = 0.040 and p _{time×genotype} < 0.001 in model 1, p _genotype = 0.021 and p _{time×genotype} < 0.001 in model 2, p _genotype = 0.042 and p _{time×genotype} < 0.001 in model 3). And the genotype-by-time interaction for rigidity and bradykinesia was significant, although there was no significant genotype effect (Rigidity subscores: p _{time×genotype} = 0.011 in model 1, p _{time×genotype} = 0.011 in model 2, p _{time×genotype} = 0.023 in model 3; Bradykinesia subscore: p _{time×genotype} < 0.001 in model 1, p _{time×genotype} < 0.001 in model 2, p _{time×genotype} = 0.002 in model 3). Detailed information about the fixed-effects statistical analyses on motor symptoms is in Supplementary Table S5.

Figure 2

Figure 2. The best-fit trajectories of motor symptom scores and non-motor symptom scores over 3 years in PD patients according to the G allele of HLA-DRB1 rs660895 (in the fully adjusted model). (A) MDS-UPDRS-III total score; (B) tremor subscore; (C) rigidity subscore; (D) bradykinesia subscore; (E) axial score; (F) NMSS score; (G) MoCA score; (H) RBD-SQ score; (I) SS-12 score.

Table 4

Table 4. Changes in motor symptoms and non-motor symptoms scores according to HLA-DRB1 rs660895 genotypes under the dominant model.

3.3.2 Effect of HLA-DRB1 rs660895 on the progression of non-motor symptoms in PD patients

As shown in the three-year trajectories of over-all non-motor symptoms progression (Figure 2F), patients with the G allele had lower NMSS score at the 3rd year follow-up (p = 0.0140 in model 1; p = 0.019 in model 2; p = 0.030 in model 3) (Table 4). Longitudinal analysis showed that NMSS score increased with time in two groups (p _time < 0.001 in all models). No significant group effect (p _genotype > 0.05 in all models), while a significant genotype-by-time interaction effect was found (p _{time×genotype} < 0.001 in all models).

Given the aforementioned results demonstrating the obvious effect of the G allele on the progression of overall non-motor symptoms, subsequent longitudinal analysis was conducted to assess the impact of the G allele status on cognitive impairment, RBD, and olfactory impairment. Our findings revealed significant time effects and time-genotype interactions on the above three aspects, but no significant group effects. Additionally, we observed three-year trajectories of the MoCA and SS-12 scores, which showed that patients with the G allele had better MoCA score (p = 0.012 in model 1; p = 0.005 in model 2; p = 0.018 in model 3) and SS-12 score (p = 0.004 in model 1; p = 0.004 in model 2; p = 0.007 in model 3) (Figures 2G–I; Table 4). It is worth noting that, although there was an interaction effect of genotype and time on RBD, the differences in the RBD-SQ score between different genotype groups did not reach statistical significance at the 3-year follow-up (p > 0.05 in all models) (Figure 2H; Table 4). Detailed information about the fixed-effects statistical analyses on non-motor symptoms is in Supplementary Table S6.

4 Discussion

This study replicated and verified an inverse relationship between HLA-DRB1 rs660895 and susceptibility to PD in Chinese Han population. Importantly, through subsequent gene-based analysis in a 3-year PD cohort, we discovered that the rs660895-G allele may has a protective effect on the longitudinal changes in both motor and non-motor symptoms of PD. This finding demonstrates the potential importance of HLA gene in the long-term progression of PD.

The highly polymorphic HLA gene plays a crucial role in acquired immunity. Several studies have explored the relationship between the HLA gene and PD in recent years (Supplementary Table S7). HLA gene, especially HLA-DRB1, widely regarded relevant to the susceptibility of PD in Caucasian population, although specifics loci differ (Hamza et al., 2010; Guo et al., 2011; Nalls et al., 2011; Ahmed et al., 2012; Wissemann et al., 2013; Nalls et al., 2014; Chuang et al., 2017; Hollenbach et al., 2019; Chang et al., 2020;Naito et al., 2021; Yu et al., 2021; Le Guen et al., 2023). The latest research once again confirmed that specific HLA-DRB1 variants are associated with a reduced risk of PD and suggested that this protective effect is primarily caused by specific amino acid polymorphisms present in most HLA-DRB1*04 subtypes (Naito et al., 2021; Yu et al., 2021; Le Guen et al., 2023). In Asians, the association between HLA and PD in GWASs is much less. Although Naito et al. identified HLA-DRB1 locus (rs504594) associated with HLA-DRB1*04 in East Asians (Naito et al., 2021), it is important to note that this finding was not replicated in two other PD-related GWAS in Asia (Foo et al., 2020; Pan et al., 2023), which did not identify any HLA alleles. Thus, in this study, we focused on the Chinese Han population and conducted a comprehensive analysis to investigate the relationship between PD and rs660896, which has a strong correlation with HLA-DRB1*04. Our findings suggest that the G allele of significantly decreases the risk of PD in different genetic inheritance models, including allele, co-dominant, and dominant models. Particularly in the dominant model, the risk of PD was found to be approximately 33–36% lower compared to the healthy control group. These results confirm the strong correlation between rs660895 and the risk of PD. While this finding appears to be inconsistent with those reported in previous studies in Asian (Foo et al., 2020; Naito et al., 2021), it is consistent with previous studies conducted in Caucasian populations (Ahmed et al., 2012; Chuang et al., 2017).

Moreover, PD is a disorder with a highly variable clinical phenotype. Understanding genetic variants that modify disease presentation and progression is crucial for gaining insights into disease mechanisms and potential therapeutic targets. Unfortunately, subsequent cross-sectional analysis grouped based on the dominant inheritance model has not shown that possessing the rs660896 G allele affects the onset age, motor phenotypes, and initial motor symptoms of PD. Meanwhile, certain studies have identified a temporal longitudinal connection between α-synuclein-specific T cell reactivity and PD (Lindestam Arlehamn et al., 2020). Thus, considering the positive effect of rs660895 on reducing the risk of PD in this study, we further extended our observations to evaluate whether the rs660895 also affects symptoms progression in our 3-year PD cohort. Notably, our longitudinal study found that HLA-DRB1 rs660895 showed evidence of strong protection against the progression of PD. During the 3-year follow-up period, both G allele carriers and non-carriers generally exhibited a decline in terms of overall motor symptoms (including rigidity, bradykinesia, and axial impairment), overall non-motor symptoms, cognitive impairment, RBD, and olfactory impairment. Nevertheless, PD patients carrying the protective G allele exhibited an obviously slower rate of both motor and non-motor symptoms progression, particularly after the 2nd year of follow-up. This suggests that the beneficial effect of the G allele is continuing, and appears greater for a long-term impact on PD progression over time. This novel evidence highlights the previously unexplored correlation between PD progression and HLA genotypes.

Although the causal role of HLA molecules in the pathogenesis of PD is not fully understood, recent studies have shown that α-synuclein peptides can induce differential T cell reactivity associated with HLA alleles (Sulzer et al., 2017; Lindestam Arlehamn et al., 2020). Yu et al. further conducted trans-ethnic fine-mapping of the MHC region and discovered that HLA-DRB1*04 subtypes (HLA-DRB1 alleles with His13) may have a protective role in the development of PD by reducing the binding affinity to α-synuclein epitopes (Yu et al., 2021). Thus, we hypothesize that the polymorphism of rs660895, which has been previously linked to a strong association with HLA-DRB1*04, could potentially impede T cell activation by altering the binding affinity between HLA-DRB1 molecules and α-synuclein. However, further research is needed to explore the precise amino acid positions and the epitope presentation to α-synuclein peptides.

This study had several noteworthy strengths, limitations and uncertainties. One key strength of our study is that, despite having a smaller sample size compared to previous large-scale GWAS studies, we were able to observe the role of rs660895 in the progression of PD through a longitudinal cohort with an excellent follow-up rate and detailed clinical information. This is a rare opportunity to identify potential protective trends in HLA-DRB1 with disease progression, providing additional clinical support toward the hypothesis that neuroinflammation mediates the progression of PD. Another unique feature of this study is that it analyzed a specific Asian population (i.e., ethnic Han Chinese), a previously under-studied group. Given the limited number of PD-related GWAS studies conducted on Asians, our research contributes to expanding the applicability of previous findings. Specifically, we investigate whether the HLA-DRB1 rs660895 polymorphism has a protective effect on PD in Asian populations.

However, there are limitations to consider in the interpretation of our findings. First, our study only focused on a single SNP, rs660895, located in HLA-DRB1 gene. Although we have discovered an association of rs660895 with the risk and progression of PD, it is still unknown whether we have identified the functional variant. Additionally, due to limitations in follow-up duration and patient compliance issues with genetic testing, we were unable to test other newly discovered potential HLA-related loci and impute and separate the HLA-DRB1*04:05 with other DRB1 variants. It may have influenced our evaluation of the protective strength of rs660895 to some extent, as it is likely that rs660895 is only in linkage disequilibrium with the causal SNP. Hence, future studies should consider conducting a more comprehensive linkage analysis to address this issue. Second, despite many clinical variables, there were no adjustments for multiple comparisons in the longitudinal analyses. It is worth noting that the severity of non-motor symptoms (secondary indicators) is closely associated with motor symptoms (main indicators) (Santos-García et al., 2020). Thus, the results of our longitudinal exploratory study are hypothesis generating and focused on magnitudes of differences rather than statistical significance on disease progression. Additional high quality longitudinal studies are needed to further evaluate this protective effect.

In conclusion, this study comprehensively investigates the relationship between HLA-DRB1 rs660895 and PD in the Han Chinese population. The findings clearly demonstrate that the rs660895-G allele is associated with a decreased risk of PD. Furthermore, the study provides additional evidence suggesting that this allele may contribute to slowing the progression of both motor and non-motor symptoms of PD. Despite some limitations, the findings of this study enhance our understanding of the genetic factors underlying PD in the Chinese population and offer a new gene target for therapeutic development.

Data availability statement

The data analyzed in this study is subject to the following licenses/restrictions: the data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions. Requests to access these datasets should be directed to dW5pb25xeXllOEBmam11LmVkdS5jbg==.

Ethics statement

The studies involving humans were approved by Ethics Committee of Fujian Medical University Union Hospital Ethics Committee of Xuanwu Hospital of Capital Medical University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

RH: Data curation, Formal analysis, Investigation, Writing – original draft. YZ: Investigation, Writing – original draft. CW: Data curation, Resources, Writing – review & editing. LC: Data curation, Investigation, Writing – review & editing. GC: Formal analysis, Methodology, Writing – review & editing. YC: Data curation, Writing – review & editing. YW: Data curation, Writing – review & editing. QY: Conceptualization, Funding acquisition, Supervision, Writing – review & editing, Resources. XC: Conceptualization, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by grants from the Project of Industry-university-research Collaborative Innovation in Fujian Province Universities (2022Y4004), the National Natural Science Foundation of China (81870995), the Natural Science Foundation of Fujian Province (2022J01727), and the Key Clinical Specialty Discipline Construction Program of China (21281003).

Acknowledgments

We thank the participating subjects for their involvement in this study as well as the staff of the Department of Neurology at Fujian Medical University Union Hospital. Our thanks also extend to CW and his team at Xuanwu Hospital of Capital Medical University for ethics and data collection.

Conflict of interest

The 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.2024.1361492/full#supplementary-material

References

Ahmed, I., Tamouza, R., Delord, M., Krishnamoorthy, R., Tzourio, C., Mulot, C., et al. (2012). Association between Parkinson's disease and the HLA-DRB1 locus. Mov. Disord. 27, 1104–1110. doi: 10.1002/mds.25035

PubMed Abstract | Crossref Full Text | Google Scholar

Auton, A., Brooks, L. D., Durbin, R. M., Garrison, E. P., Kang, H. M., Korbel, J. O., et al. (2015). A global reference for human genetic variation. Nature 526, 68–74. doi: 10.1038/nature15393

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

PubMed Abstract | Crossref Full Text | Google Scholar

Buhler, S., and Sanchez-Mazas, A. (2011). HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events. PLoS One 6:e14643. doi: 10.1371/journal.pone.0014643

PubMed Abstract | Crossref Full Text | Google Scholar

Chang, K. H., Wu, Y. R., Chen, Y. C., Fung, H. C., and Chen, C. M. (2020). Association of genetic variants within HLA-DR region with Parkinson's disease in Taiwan. Neurobiol. Aging 87, 140.e13–140.e18. doi: 10.1016/j.neurobiolaging.2019.11.002

PubMed Abstract | Crossref Full Text | Google Scholar

Chaudhuri, K. R., Martinez-Martin, P., Brown, R. G., Sethi, K., Stocchi, F., Odin, P., et al. (2007). The metric properties of a novel non-motor symptoms scale for Parkinson's disease: results from an international pilot study. Mov. Disord. 22, 1901–1911. doi: 10.1002/mds.21596

PubMed Abstract | Crossref Full Text | Google Scholar

Chuang, Y. H., Lee, P. C., Vlaar, T., Mulot, C., Loriot, M. A., Hansen, J., et al. (2017). Pooled analysis of the HLA-DRB1 by smoking interaction in Parkinson disease. Ann. Neurol. 82, 655–664. doi: 10.1002/ana.25065

PubMed Abstract | Crossref Full Text | Google Scholar

Çomoğlu, Ş., Orhan, K. S., Kocaman, S., Çelik, M., Keleş, N., and Değer, K. (2015). Olfactory function assessment of blind subjects using the Sniffin' sticks test. Otolaryngol. Head Neck Surg. 153, 286–290. doi: 10.1177/0194599815583975

PubMed Abstract | Crossref Full Text | Google Scholar

Eisinger, R. S., Martinez-Ramirez, D., Ramirez-Zamora, A., Hess, C. W., Almeida, L., Okun, M. S., et al. (2020). Parkinson's disease motor subtype changes during 20 years of follow-up. Parkinsonism Relat. Disord. 76, 104–107. doi: 10.1016/j.parkreldis.2019.05.024

PubMed Abstract | Crossref Full Text | Google Scholar

Fernandez Vina, M. A., Hollenbach, J. A., Lyke, K. E., Sztein, M. B., Maiers, M., Klitz, W., et al. (2012). Tracking human migrations by the analysis of the distribution of HLA alleles, lineages and haplotypes in closed and open populations. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 367, 820–829. doi: 10.1098/rstb.2011.0320

PubMed Abstract | Crossref Full Text | Google Scholar

Foo, J. N., Chew, E. G. Y., Chung, S. J., Peng, R., Blauwendraat, C., Nalls, M. A., et al. (2020). Identification of risk loci for Parkinson disease in Asians and comparison of risk between Asians and Europeans: a genome-wide association study. JAMA Neurol. 77, 746–754. doi: 10.1001/jamaneurol.2020.0428

PubMed Abstract | Crossref Full Text | Google Scholar

Goetz, C. G., Poewe, W., Rascol, O., Sampaio, C., Stebbins, G. T., Counsell, C., et al. (2004). Movement Disorder Society task force report on the Hoehn and Yahr staging scale: status and recommendations. Mov. Disord. 19, 1020–1028. doi: 10.1002/mds.20213

PubMed Abstract | Crossref Full Text | Google Scholar

Goetz, C. G., Tilley, B. C., Shaftman, S. R., Stebbins, G. T., Fahn, S., Martinez-Martin, P., et al. (2008). Movement Disorder Society-sponsored revision of the unified Parkinson's disease rating scale (MDS-UPDRS): scale presentation and clinimetric testing results. Mov. Disord. 23, 2129–2170. doi: 10.1002/mds.22340

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, Y., Deng, X., Zheng, W., Xu, H., Song, Z., Liang, H., et al. (2011). HLA rs3129882 variant in Chinese Han patients with late-onset sporadic Parkinson disease. Neurosci. Lett. 501, 185–187. doi: 10.1016/j.neulet.2011.05.245

PubMed Abstract | Crossref Full Text | Google Scholar

Hamza, T. H., Zabetian, C. P., Tenesa, A., Laederach, A., Montimurro, J., Yearout, D., et al. (2010). Common genetic variation in the HLA region is associated with late-onset sporadic Parkinson's disease. Nat. Genet. 42, 781–785. doi: 10.1038/ng.642

PubMed Abstract | Crossref Full Text | Google Scholar

Herrero, M. T., Estrada, C., Maatouk, L., and Vyas, S. (2015). Inflammation in Parkinson's disease: role of glucocorticoids. Front. Neuroanat. 9:32. doi: 10.3389/fnana.2015.00032

PubMed Abstract | Crossref Full Text | Google Scholar

Hollenbach, J. A., Norman, P. J., Creary, L. E., Damotte, V., Montero-Martin, G., Caillier, S., et al. (2019). A specific amino acid motif of HLA-DRB1 mediates risk and interacts with smoking history in Parkinson's disease. Proc. Natl. Acad. Sci. USA 116, 7419–7424. doi: 10.1073/pnas.1821778116

PubMed Abstract | Crossref Full Text | Google Scholar

Le Guen, Y., Luo, G., Ambati, A., Damotte, V., Jansen, I., Yu, E., et al. (2023). Multiancestry analysis of the HLA locus in Alzheimer's and Parkinson's diseases uncovers a shared adaptive immune response mediated by HLA-DRB1*04 subtypes. Proc. Natl. Acad. Sci. USA 120:e2302720120. doi: 10.1073/pnas.2302720120

PubMed Abstract | Crossref Full Text | Google Scholar

Lindestam Arlehamn, C. S., Dhanwani, R., Pham, J., Kuan, R., Frazier, A., Rezende Dutra, J., et al. (2020). α-Synuclein-specific T cell reactivity is associated with preclinical and early Parkinson's disease. Nat. Commun. 11:1875. doi: 10.1038/s41467-020-15626-w

PubMed Abstract | Crossref Full Text | Google Scholar

Naito, T., Satake, W., Ogawa, K., Suzuki, K., Hirata, J., Foo, J. N., et al. (2021). Trans-ethnic fine-mapping of the major histocompatibility complex region linked to Parkinson's disease. Mov. Disord. 36, 1805–1814. doi: 10.1002/mds.28583

PubMed Abstract | Crossref Full Text | Google Scholar

Nalls, M. A., Blauwendraat, C., Vallerga, C. L., Heilbron, K., Bandres-Ciga, S., Chang, D., et al. (2019). Identification of novel risk loci, causal insights, and heritable risk for Parkinson's disease: a meta-analysis of genome-wide association studies. Lancet Neurol. 18, 1091–1102. doi: 10.1016/S1474-4422(19)30320-5

PubMed Abstract | Crossref Full Text | Google Scholar

Nalls, M. A., Pankratz, N., Lill, C. M., Do, C. B., Hernandez, D. G., Saad, M., et al. (2014). Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson's disease. Nat. Genet. 46, 989–993. doi: 10.1038/ng.3043

PubMed Abstract | Crossref Full Text | Google Scholar

Nalls, M. A., Plagnol, V., Hernandez, D. G., Sharma, M., Sheerin, U. M., Saad, M., et al. (2011). Imputation of sequence variants for identification of genetic risks for Parkinson's disease: a meta-analysis of genome-wide association studies. Lancet 377, 641–649. doi: 10.1016/S0140-6736(10)62345-8

PubMed Abstract | Crossref Full Text | Google Scholar

Nasreddine, Z. S., Phillips, N. A., Bédirian, V., Charbonneau, S., Whitehead, V., Collin, I., et al. (2005). The Montreal cognitive assessment, MoCA: a brief screening tool for mild cognitive impairment. J. Am. Geriatr. Soc. 53, 695–699. doi: 10.1111/j.1532-5415.2005.53221.x

PubMed Abstract | Crossref Full Text | Google Scholar

Pan, H., Liu, Z., Ma, J., Li, Y., Zhao, Y., Zhou, X., et al. (2023). Genome-wide association study using whole-genome sequencing identifies risk loci for Parkinson's disease in Chinese population. NPJ Parkinsons Dis. 9:22. doi: 10.1038/s41531-023-00456-6

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

Santos-García, D., de Deus, F. T., Suárez Castro, E., Aneiros Díaz, A., McAfee, D., Catalán, M. J., et al. (2020). Non-motor symptom burden is strongly correlated to motor complications in patients with Parkinson's disease. Eur. J. Neurol. 27, 1210–1223. doi: 10.1111/ene.14221

PubMed Abstract | Crossref Full Text | Google Scholar

Satake, W., Nakabayashi, Y., Mizuta, I., Hirota, Y., Ito, C., Kubo, M., et al. (2009). Genome-wide association study identifies common variants at four loci as genetic risk factors for Parkinson's disease. Nat. Genet. 41, 1303–1307. doi: 10.1038/ng.485

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

PubMed Abstract | Crossref Full Text | Google Scholar

Simón-Sánchez, J., Schulte, C., Bras, J. M., Sharma, M., Gibbs, J. R., Berg, D., et al. (2009). Genome-wide association study reveals genetic risk underlying Parkinson's disease. Nat. Genet. 41, 1308–1312. doi: 10.1038/ng.487

PubMed Abstract | Crossref Full Text | Google Scholar

Stebbins, G. T., Goetz, C. G., Burn, D. J., Jankovic, J., Khoo, T. K., and Tilley, B. C. (2013). How to identify tremor dominant and postural instability/gait difficulty groups with the movement disorder society unified Parkinson's disease rating scale: comparison with the unified Parkinson's disease rating scale. Mov. Disord. 28, 668–670. doi: 10.1002/mds.25383

PubMed Abstract | Crossref Full Text | Google Scholar

Stiasny-Kolster, K., Mayer, G., Schäfer, S., Möller, J. C., Heinzel-Gutenbrunner, M., and Oertel, W. H. (2007). The REM sleep behavior disorder screening questionnaire—a new diagnostic instrument. Mov. Disord. 22, 2386–2393. doi: 10.1002/mds.21740

PubMed Abstract | Crossref Full Text | Google Scholar

Sulzer, D., Alcalay, R. N., Garretti, F., Cote, L., Kanter, E., Agin-Liebes, J., et al. (2017). T cells from patients with Parkinson's disease recognize α-synuclein peptides. Nature 546, 656–661. doi: 10.1038/nature22815

PubMed Abstract | Crossref Full Text | Google Scholar

Tansey, M. G., and Romero-Ramos, M. (2019). Immune system responses in Parkinson's disease: early and dynamic. Eur. J. Neurosci. 49, 364–383. doi: 10.1111/ejn.14290

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

von Bernhardi, R., Eugenín-von Bernhardi, L., and Eugenín, J. (2015). Microglial cell dysregulation in brain aging and neurodegeneration. Front. Aging Neurosci. 7:124. doi: 10.3389/fnagi.2015.00124

PubMed Abstract | Crossref Full Text | Google Scholar

Wissemann, W. T., Hill-Burns, E. M., Zabetian, C. P., Factor, S. A., Patsopoulos, N., Hoglund, B., et al. (2013). Association of Parkinson disease with structural and regulatory variants in the HLA region. Am. J. Hum. Genet. 93, 984–993. doi: 10.1016/j.ajhg.2013.10.009

PubMed Abstract | Crossref Full Text | Google Scholar

Yan, Z., Yang, W., Wei, H., Dean, M. N., Standaert, D. G., Cutter, G. R., et al. (2021). Dysregulation of the adaptive immune system in patients with early-stage Parkinson disease. Neurol. Neuroimmunol. Neuroinflamm. 8:1036. doi: 10.1212/NXI.0000000000001036

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, E., Ambati, A., Andersen, M. S., Krohn, L., Estiar, M. A., Saini, P., et al. (2021). Fine mapping of the HLA locus in Parkinson's disease in Europeans. NPJ Parkinsons Dis. 7:84. doi: 10.1038/s41531-021-00231-5

PubMed Abstract | Crossref Full Text | Google Scholar