You're viewing our updated article page. If you need more time to adjust, you can return to the old layout.

ORIGINAL RESEARCH article

Front. Neurol., 04 April 2025

Sec. Applied Neuroimaging

Volume 16 - 2025 | https://doi.org/10.3389/fneur.2025.1477827

Association of free-water imaging data for the cholinergic nucleus with the motor function and subtypes in Parkinson’s disease

  • 1. Department of Neurology, Affiliated Zhejiang Hospital Zhejiang University School of Medicine, Hangzhou, China

  • 2. Department of Neurology, The Second Affiliated Hospital Zhejiang University School of Medicine, Hangzhou, China

  • 3. Department of Radiology, Affiliated Zhejiang Hospital Zhejiang University School of Medicine, Hangzhou, China

Article metrics

View details

1,2k

Views

339

Downloads

Abstract

Background:

Despite the importance of clinical heterogeneity in Parkinson’s disease (PD), its underlying pathophysiology remains unclear.

Objective:

This study aimed to distinguish the association of free-water (FW) imaging data for the cholinergic nuclei with the motor subtypes of PD.

Methods:

The study included 150 cases of idiopathic PD from the Parkinson’s Progression Markers Initiative cohort. FW imaging, including FW-corrected diffusion tensor imaging, was used to extract structural metrics from cholinergic nucleus 4 (Ch4) in the basal forebrain and the pedunculopontine nucleus. The motor subtypes were classified as tremor-dominant (TD, n = 99) and non-tremor-dominant (non-TD, n = 51). Statistical analyses were performed at baseline and the 4-year follow-up.

Results:

At baseline, FW value for Ch4 (FW-Ch4) was correlated with the tremor subscore, while FW-corrected fractional anisotropy in Ch4 (FA-t-Ch4) was negatively correlated with the rigidity subscore. However, the TD and non-TD groups showed no differences in cholinergic FW imaging data. Among the 84 patients who were followed-up, 36.36% (20/55) in the TD group and 34.48% (10/29) in the non-TD group showed a subtype shift after 4 years. Multivariate binary logistic regression analysis showed that the normalized FW value for Ch4 (nFW-Ch4) was a predictor of subtype at the 4-year follow-up (p = 0.041). In the TD subgroup, both nFW-Ch4 (p = 0.015) and normalized FW-corrected mean diffusivity in Ch4 (MD-t-Ch4) (p = 0.013) predicted subtype stability. The area under the receiver operating characteristic curve values were 0.69 and 0.73, respectively.

Conclusion:

Tremor and rigidity subscores were correlated with Ch4 FW imaging data. Moreover, Ch4 FW imaging predicted the motor subtype at the 4-year follow-up, especially identifying potential postural instability and gait difficulty (PIGD) subtype converters from the TD group.

Introduction

Parkinson’s disease (PD) is a progressive neurodegenerative disease characterized by heterogeneous clinical and neuropathological manifestations. The most common manifestations include the classic triad of motor symptoms, namely, resting tremor, rigidity, and bradykinesia and postural instability. Establishment of different subcategories is an important priority in clinical research on PD (1), since these subtypes have implications for diagnosis, prognosis, and expected treatment responses (2, 3). Although different algorithms to determine PD motor subtypes have been proposed, the subtypes based on Jankovic’s classification system, which categorizes patients into tremor-dominant (TD), indeterminate (ID), and postural instability and gait difficulty (PIGD) subtypes on the basis of motor dominance, are most frequently studied in the literature (4–6) and more suitable for identifying non-motor abnormalities (6). Nevertheless, due to the instability of motor subtypes, the Parkinson’s Progression Markers Initiative (PPMI) consortium decided to categorize patients with PD into TD and non-TD (PIGD and ID) subtypes (7).

In addition to experiencing motor symptoms, patients with PD also show non-motor symptoms (NMS) such as sleep disturbances, autonomic dysfunction, cognitive decline, and psychiatric symptoms. NMS may even appear earlier than the classic motor features in the course of PD (8). NMS are valuable for the differential diagnosis of PD from atypical parkinsonisms (9) and may precede motor symptoms and influence patients’ quality of life (10). For example, mild cognitive impairment is associated with poor gait performance in patients with PD (11). Both motor and NMS vary among patients with PD and show complex interactions with each other.

Despite the substantial importance of clinical heterogeneity in PD, its underlying pathophysiology largely remains unclear. Numerous extrastriatal pathologies have been recently described as correlates for NMS. While the severity of dopaminergic depletion in different subtypes has been evaluated (2), the pathogenic changes in other non-dopamine systems are not well understood. Progressive degeneration of the cholinergic system has been widely implicated in PD (12). Besides its role in cognitive deterioration, this form of degeneration has been also proposed to be the functional link between cholinergic alterations and PD motor symptoms, sleep abnormalities, autonomic dysfunction, and altered olfactory function (12, 13). Cholinergic system-related NMS also affect PD motor classification and its stability (14, 15). Therefore, we speculated that changes in the cholinergic system, most importantly the nucleus basalis magnocellularis (NBM) and the pedunculopontine nucleus (PPN), may be associated with motor impairment and may influence the motor classification.

Diffusion tensor imaging (DTI) is a sensitive method for isolating white matter (WM) trajectories and exploring microstructural integrity in PD (16). DTI models each voxel as a single compartment. Free-water (FW) DTI uses an advanced 2-compartment diffusion model: one from the tissue and the other from free diffusing water (17); recent advancements in this technique have been shown to minimize the partial-volume effects of the cerebrospinal fluid and reflect neuroinflammation (18). In one study, FW imaging findings for the substantia nigra in different stages of PD were evaluated and were found to be a useful disease-state biomarker (16). FW imaging has also been shown to be informative for assessment of the cholinergic nucleus and cholinergic projection system (19–21), and has been shown to be a useful biomarker of cognitive symptoms in neurodegenerative disease. In the current study, we used cholinergic FW imaging as well as clinical data from the PPMI cohort of patients with idiopathic de novo PD to evaluate (1) the clinical and cholinergic FW imaging data in relation to the PD motor subtype at baseline and at the 4-year follow-up; (2) the correlation of cholinergic microstructural integrity and motor subscores; and (3) the ability of FW markers to detect changes in the PD phenotype after 4 years.

Materials and methods

Study sample

All the clinical and imaging data used in the study were downloaded from the PPMI database1 in August 2023. The PPMI is a comprehensive longitudinal, international, multicenter database consisting of clinical, genetic, neuroimaging and blood/cerebrospinal fluid (CSF) biomarker data of patients with early PD.

The inclusion criteria for the PD cohort in PPMI were as follows: idiopathic PD, early-stage disease classified as Hoehn–Yahr stage 1–2 at baseline, T1 imaging and DTI data available at baseline, and absence of anti-PD or anticholinergic medication. Each PPMI site received approval from its relevant institutional review board, and all participants provided written informed consent before participating in the project.

Clinical evaluation and subtype classification

All baseline clinical evaluations were performed without dopaminergic replacement therapy, and 4-year follow-up evaluations were performed in the off-medication state. The levodopa-equivalent daily dose (LEDD; mg/day) was calculated from the dosages of the prescribed dopaminergic drugs using the standard conversion table (22). Motor dysfunction and disease severity were assessed on the basis of the UPDRS part II and III scores and the modified Hoehn and Yahr (H-Y) stages, respectively. For each patient, we obtained the rigidity subscore (total score of item 3.3), tremor subscore (sum of the scores for items 2.10 and 3.15–3.18), and PIGD score (sum of the scores for items 2.12, 2.13, and 3.10–3.12) from the PPMI. The TD/PIGD phenotype of PD classification was based on the mean tremor subscore/mean PIGD score ratio. If the ratio was >1.15 or the PIGD score was 0 and tremor subscore was >0, then the participant was classified into the TD subtype; if the ratio was <0.9, then the participant was classified into the PIGD subtype; and if the ratio was 0.9–1.15 or if the tremor and PIGD scores were 0, then the participant was classified into the ID subtype (5). The PIGD and ID subtypes were together classified as the “non-tremor” type.

Additionally, evaluations using Montreal Cognitive Assessment (MoCA), Parkinson’s disease-Autonomic (SCOPA-AUT), Geriatric Depression Scale (Short Version), State–Trait Anxiety Inventory (STAI), Epworth Sleepiness Scale (ESS), and Rapid Eye Movement Scale (REM) were also performed at baseline and at the 4-year follow-up.

Dopamine transporter imaging

Dopamine transporter (DAT) imaging using I123 Ioflupane single-photon emission computed tomography (SPECT) was acquired at PPMI imaging centers in accordance with the PPMI imaging protocol, and sent to the Institute for Neurodegenerative Disorders for processing and calculation. SPECT raw projection data were imported to a HERMES (Hermes Medical Solutions, Skeppsbron 44, 111 30 Stockholm, Sweden) system for iterative reconstruction, then transferred to the PMOD (PMOD Technologies, Zurich, Switzerland) for subsequent processing. After attenuation correction, Gaussian filtering, and alignment and identification of the regions of interest, the striatal binding ratio was calculated as the ratio of the values for the target and reference regions (see text footnote 1). We acquired the data for the mean caudate, putamen, and striatum uptake relative to the uptake of the occipital area for analysis.

MRI protocol and processing

T1-weighted MRI and DTI were performed using Siemens 3 T scanners at multiple centers. The detailed MRI protocol is available in the PPMI database. We downloaded the baseline image data of the study participants. DTI data were preprocessed for denoising and removal of Gibbs artifacts and correction of head motion, eddy current-induced distortion, and intensity bias using the FMRIB Software Library (FSL). Multi-fiber directions within each voxel of the brain were estimated based on the ball-and-two-sticks model (23). With the processed images, fractional anisotropy (FA) and mean diffusivity (MD) images for each participant were obtained by using the “dtifit” function within FSL. FW, FW-corrected FA (FA-t), and FW-corrected MD (MD-t) values were derived from preprocessed DTI data using the two-compartment free-water elimination model (24). This model, which separates tissue-specific diffusion from the contribution of extracellular free water, was implemented using the Diffusion Imaging in Python (DIPY) package (version 1.4.0) (24, 25) with publicly available reference code2 (TensorImaging). The image processing procedure is summarized in Figure 1.

Figure 1

Figure 1

Overview of image processing. (A) Raw diffusion data; (B) Pre-processing diffusion data; (C) Free water (FW) related image, (C.1) FW image, (C.2) FW-corrected fractional anisotropy (FA) image, (C.3) FW-corrected mean diffusivity (MD) image; (D) Diffusion tensor imaging (DTI) related image, (D.1) FA image, (D.2) MD image; (E.1) PPN mask in standard space, (E.2) NBM mask in standard space; (F.1) Standard template where PPN mask are; (F.2) Standard template where NBM are; (G) Individual’s T1 image; (H) Individual’s FA image; (I.1) PPN mask in individual’s diffusion space, (I.2) NBM mask in individual’s diffusion mask. ANTs, Advanced Normalization Tools.

Regions of interest

The NBM mask was based on the cholinergic nucleus 4 (Ch4) region of a cytoarchitectonic map of the cholinergic basal forebrain in the Montreal Neurological Institute (MNI) T1 space, which was derived by combining the findings of histological and MRI assessments of a postmortem brain (26). For the PPN, the mask was based on a PPN-specific high-resolution atlas, which was obtained by combining a recently constructed MNI-space unbiased quantitative susceptibility mapping atlas, a myelin staining histology human brain atlas, and an implicit representation-based self-supervised image super-resolution technique (27).

Extraction of diffusivity metrics and volumetry for NBM and PPN

The mean FW, MD-t, and FA-t values were extracted from the NBM and the PPN. To ensure that these were constrained to gray matter (GM), the NBM was masked with the GM mask estimated from the segmentation of T1-weighted images (20). For the PPN, which had WM pathways from the brainstem running through it, the analysis was restricted to voxels with fractional anisotropy between 0.47 and 0.77, in accordance with previous studies (19, 20). To evaluate cell body damage in the cholinergic neurons, the NBM and PPN volumes in each individual’s diffusion space were calculated. The final NBM and PPN volumes were adjusted on the basis of the total intracranial volume [estimated using the Freesurfer’s SAMSEG segmentation (28)] to account for between-subject variability in head size.

Statistical analysis

Statistical analysis was performed using Statistical Package for Social Sciences (SPSS 23.0) software. Continuous variables were expressed as mean ± standard deviation and compared using independent-samples t-test (normally distributed data) or the exact Mann–Whitney U-test (non-normally distributed data). Categorical variables were expressed as frequencies and compared using the chi-square test. Spearman’s correlation analysis was performed to evaluate the relationships between cholinergic imaging data and motor subscores at baseline. Multivariable binary logistic regression analysis was used to investigate the association of the baseline variables in the TD and non-TD subtypes and the stability of the subtypes at the 4-year follow-up. Zero-mean normalization was used for cholinergic imaging data in regression analysis.

Results

Clinical characteristics at baseline

A total of 150 patients were eligible for the baseline clinical evaluation. The baseline demographic and clinical data in the TD and non-TD groups are summarized in Table 1. The two groups showed significant differences in the years of education; UPDRS II, PIGD, and ESS scores; and the tremor subscore (p < 0.05). Although the TD group tended to show lower values on DAT imaging, significant intergroup differences were obtained only for the mean caudate (t = 2.143, p = 0.034). In a cross-sectional comparison of the baseline FW, FA-t, MD-t, and voxel values on a priori regions of the Ch4 and PPN between the TD and non-TD groups, none of these values for the regions of interest showed intergroup differences (Table 1).

Table 1

Clinical variables Baseline 4-year follow-up
Missing TD group
(n = 99)
Non-TD group
(n = 51)
p value Missing TD group
(n = 45)
Non-TD group
(n = 39)
p value
Demographic features
Male (n, %) 0/0 62 (62.63%) 34 (66.67%) 0.625 0/0 27 (60.0%) 30 (76.92%) 0.098
Age at visit (years) 0/0 61.38 ± 8.79 60.58 ± 10.87 0.795 0/0 65.44 ± 8.56 63.31 ± 10.45 0.306
Age at onset (years) 1/2 59.04 ± 9.02 59.19 ± 11.03 0.936 0/0 59.58 ± 9.12 57.55 ± 10.42 0.346
Education (years) 0/0 15.93 ± 2.77 14.24 ± 3.03 0.002 0/0 15.22 ± 3.28 15.51 ± 2.79 0.666
PD duration (months) 0/0 6.98 ± 7.21 6.20 ± 6.07 0.702 0/0 55.38 ± 6.30 54.82 ± 8.41 0.730
LEDD - - - - 0/0 516.55 ± 319.24 536.74 ± 285.30 0.762
Clinical characteristics
Hoehn & Yahr 0/0 1.54 ± 0.50 1.65 ± 0.48 0.318 1/1 1.71 ± 0.51 2.03 ± 0.58 0.010
UPDRS II score 0/0 4.85 ± 3.66 6.65 ± 4.40 0.009 0/0 7.64 ± 4.56 9.38 ± 5.54 0.188
UPDRS III score 0/0 21.11 ± 9.18 19.49 ± 8.68 0.299 0/1 28.64 ± 10.86 29.08 ± 15.42 0.714
UPDRS total score 0/0 30.81 ± 13.10 32.08 ± 14.08 0.584 0/1 44.20 ± 15.54 46.82 ± 21.75 0.830
Rigidity subscore 0/0 4.00 ± 2.69 4.53 ± 2.91 0.268 0/0 5.36 ± 2.98 6.95 ± 3.52 0.028
Tremor subscore 0/0 6.43 ± 3.19 1.92 ± 2.05 <0.001 1/0 7.56 ± 3.08 2.10 ± 2.47 <0.001
PIGD subscore 0/0 0.75 ± 0.83 1.88 ± 1.36 <0.001 1/0 1.82 ± 1.57 4.31 ± 5.01 0.026
Non-motor symptoms
MoCA 0/0 27.62 ± 2.04 27.39 ± 2.23 0.538 0/1 26.98 ± 2.58 27.68 ± 2.53 0.213
ESS 0/0 5.41 ± 3.17 6.88 ± 3.32 0.009 0/0 7.11 ± 4.39 9.13 ± 4.71 0.046
GDS 0/0 2.08 ± 2.36 2.90 ± 2.83 0.061 0/0 2.18 ± 2.54 2.59 ± 2.11 0.426
SCOPA-AUT 2/1 9.00 ± 6.47 8.88 ± 4.65 0.907 0/0 11.18 ± 5.88 13.18 ± 6.55 0.144
STAI 0/0 65.24 ± 19.17 68.27 ± 15.45 0.330 0/0 63.47 ± 17.63 65.44 ± 14.99 0.586
REM 0/0 3.89 ± 2.51 4.33 ± 2.59 0.333 1/0 4.05 ± 2.64 5.87 ± 3.66 0.025
DAT imaging
Mean caudate score 0/0 1.94 ± 0.51 1.75 ± 0.51 0.034 8/5 1.55 ± 0.44 1.38 ± 0.42 0.114
Mean putamen score 0/0 0.81 ± 0.24 0.77 ± 0.29 0.441 8/5 0.57 ± 0.18 0.54 ± 0.21 0.216
Mean striatum score 0/0 1.37 ± 0.36 1.26 ± 0.38 0.079 8/5 1.06 ± 0.29 0.96 ± 0.29 0.161
Cholinergic imaging
Ch4 volume (mm3) 0/0 416.27 ± 72.07 424.01 ± 0.41 0.566 0/0 439.43 ± 64.12 422.56 ± 68.87 0.249
FA-t-Ch4 0/0 0.50 ± 0.60 0.50 ± 0.64 0.771 0/0 0.50 ± 0.56 0.50 ± 0.67 0.559
MD-t-Ch4 *10−3 0/0 0.63 ± 0.14 0.60 ± 0.10 0.096 0/0 0.65 ± 0.16 0.58 ± 0.09 0.029
FW-Ch4 0/0 0.37 ± 0.05 0.37 ± 0.04 0.105 0/0 0.38 ± 0.51 0.36 ± 0.33 0.042
PPN volume (mm3) 0/0 65.94 ± 22.50 71.15 ± 25.33 0.200 0/0 68.84 ± 23.05 68.34 ± 20.50 0.918
FA-t-PPN 0/0 0.89 ± 0.37 0.90 ± 0.34 0.367 0/0 0.89 ± 0.04 0.90 ± 0.27 0.344
MD-t-PPN*10−3 0/0 0.36 ± 0.04 0.35 ± 0.05 0.906 0/0 0.36 ± 0.04 0.36 ± 0.04 0.557
FW-PPN 0/0 0.36 ± 0.28 0.36 ± 0.04 0.539 0/0 0.35 ± 0.30 0.36 ± 0.02 0.520

Demographic and clinical data in Parkinson’s disease patients with different subtypes at baseline and the 4-year follow-up.

Results are expressed as mean ± standard deviation for continuous variables and as frequencies for categorical variables. Bold font indicates significant differences (p < 0.05).

TD, tremor-dominant subtype; LEDD, levodopa-equivalent daily dose; UPDRS, Unified Parkinson’s Disease Rating Scale; PIGD, postural instability and gait difficulty; MoCA, Montreal Cognitive Assessment; ESS, Epworth Sleepiness Scale; SCOPA-AUT, Parkinson’s disease-Autonomic; GDS, Geriatric Depression Scale; STAI, State–Trait Anxiety Inventory; REM, Rapid Eye Movement scale; DAT, dopamine transporter; Ch4, cholinergic nucleus 4; FA-t-Ch4, free-water-corrected fractional anisotropy in Ch4; MD-t-Ch4, free-water-corrected mean diffusivity in Ch4; FW-Ch4, free-water value for Ch4; PPN, pedunculopontine nucleus; FA-t-PPN, free-water-corrected fractional anisotropy in PPN; MD-t-PPN, free-water-corrected mean diffusivity in PPN; FW-PPN, free-water value for PPN.

Correlation of cholinergic FW imaging findings and motor symptoms

Analysis of the correlation between the findings of cholinergic FW imaging and clinical characteristics was performed in all PD patients (Figure 2A). In brief, age was negatively correlated with the volume of Ch4 and positively correlated with the FW value for Ch4 (FW-Ch4) and MD-t-Ch4. Years of education was negatively correlated with FA-t-Ch4 and positively correlated with MD-t-Ch4. FW-Ch4 was correlated with the UPDRS II and UPDRS total scores and the tremor subscore, while FA-t-Ch4 was negatively correlated with the rigidity subscore (Figures 2BD). Our findings also showed trends for correlations between the UPDRS III score and FA-t-Ch4 (r = −0.161, p = 0.050) and MD-t-Ch4 (r = 0.155, p = 0.057) and between the tremor subscore and MD-t-Ch4 (r = 0.155, p = 0.058). Nevertheless, no correlation was observed between any of the clinical signs and the results of PPN diffusion tensor water imaging.

Figure 2

Figure 2

Spearman’s correlation analysis of cholinergic free-water imaging data and clinical data. (A) The correlation coefficients are shown in the table. **indicates p < 0.01 and *indicates p < 0.05. (B) Correlation between rigidity score and FW-corrected fractional anisotropy of cholinergic nucleus 4 (FA-t-Ch4); (C) correlation between the tremor subscore and the free-water value for cholinergic nucleus 4 (FW-Ch4); (D) correlation between the tremor subscore and the FW-corrected mean diffusivity of cholinergic nucleus 4 (MD-t-Ch4).

Stability and fluctuation of PD subtypes

Data for the stability and fluctuation of PD subtypes after 4 years are presented in Figure 3. The 4-year follow-up assessments were conducted in 84 (56.00%) individuals, and the remaining 66 (44.00%) individuals were excluded because of the lack of off-stage follow-up data. Among the 84 PD patients, 55 (65.48%) and 29 (34.52%) patients were categorized into the TD and non-TD (indeterminate, 13; PIGD, 16) types, respectively, at baseline. At the 4-year follow-up, the number of patients categorized into the TD and non-TD types was 45 (53.57%) and 39 (46.43%; indeterminate, 14 and PIGD, 25), respectively. Thus, 36.36% (20/55) of the cases initially classified under the TD subtype changed to the non-TD subtype, while 34.48% (10/29) of the cases initially classified under the non-TD subtype changed to the TD subtype.

Figure 3

Figure 3

Fluctuations in Parkinson’s disease subtypes over the 4-year follow-up period.

Association of cholinergic FW imaging data with PD subtypes at the 4-year follow-up

The clinical characteristics of patients categorized by the 4-year follow-up data are shown in Table 1. The H-Y stage and ESS and REM scores indicated greater severity in the non-TD group. However, other characteristics, including DAT imaging findings, showed no significant differences. Additionally, MD-t-Ch4 and FW-Ch4 were significantly higher in the TD group.

Multivariate binary logistic regression analysis showed that normalized FW-Ch4 (nFW-Ch4) was a predictor of subtype at the 4-year follow-up. In the initial TD subgroup, both nFW-Ch4 and normalized MD-t-Ch4 (nMD-t-Ch4) could predict subtype stability separately; however, in the initial non-TD group, nMD-t-Ch4 was a predictor for subtype stability in model 3 but not in model 2 (Table 2). Receiver operating characteristic curve analysis showed that the initial nFW-Ch4 (OR: 0.304, 95% CI: 0.116–0.796, p = 0.015) and nMD-t-Ch4 (OR: 0.154, 95% CI: 0.035–0.672, p = 0.013) could adequately predict subtype stability at the 4-year follow-up in the initial TD group (Figure 4).

Table 2

Total sample Initial TD group Initial non-TD group
B OR 95% CI p value B OR 95% CI p value B OR 95% CI p value
Model 1 based on FW-Ch4
Age −0.05 0.95 0.90–1.01 0.080 −0.88 0.92 0.84–1.00 0.060 −0.05 0.95 0.86–1.04 0.273
Male 0.65 1.91 0.63–5.79 0.255 1.94 6.97 1.00–48.82 0.051 −0.09 0.92 0.12–6.88 0.932
REM 0.14 1.15 0.95–1.38 0.155 0.29 1.33 1.01–1.76 0.042 −0.03 0.97 0.69–1.35 0.850
ESS 0.02 1.02 0.90–1.17 0.729 −0.04 0.96 0.79–1.17 0.695 −0.02 0.98 0.79–1.21 0.856
H-Y 1.19 3.29 1.14–9.53 0.028 1.75 5.75 0.51–64.36 0.156 1.21 3.34 0.65–17.07 0.147
nFW-Ch4 −0.64 0.53 0.29–0.98 0.041 −1.19 0.30 0.12–0.80 0.015 0.04 1.04 0.24–4.42 0.963
Model 2 based on MD-t-Ch4
Age −0.05 0.95 0.90–1.01 0.102 −0.07 0.93 0.85–1.02 0.139 −0.08 0.93 0.84–1.02 0.129
Male 0.56 1.74 0.58–5.22 0.321 2.02 7.51 1.00–56.37 0.050 0.68 1.98 0.21–18.68 0.552
REM 0.14 1.15 0.95–1.38 0.151 0.28 1.33 0.99–1.78 0.063 −0.08 0.92 0.65–1.31 0.645
ESS 0.03 1.10 0.91–1.18 0.644 0.04 1.04 0.84–1.29 0.731 0.03 1.03 0.82–1.28 0.817
H-Y 1.04 0.95 0.90–1.01 0.049 1.01 2.75 0.23–33.58 0.428 1.14 3.13 0.68–14.55 0.145
nMD-t-Ch4 −0.60 0.55 0.29–1.03 0.060 −1.88 0.15 0.04–0.67 0.013 1.18 3.25 0.60–17.63 0.171
Model 3 based on FW-Ch4 & MD-t-Ch4
Age −0.05 0.95 0.90–1.01 0.089 −0.07 0.93 0.85–1.03 0.158 −0.07 0.93 0.82–1.06 0.274
Male 0.63 1.88 0.62–5.73 0.267 2.00 7.39 0.96–57.05 0.055 0.38 1.45 0.07–32.77 0.814
REM 0.14 1.15 0.95–1.38 0.153 0.28 1.32 0.98–1.79 0.068 −0.64 0.53 0.25–1.11 0.090
ESS 0.03 1.03 0.90–1.17 0.710 0.41 1.04 0.83–1.31 0.721 0.07 1.07 0.80–1.43 0.652
H-Y 1.16 3.18 1.08–9.36 0.036 0.99 2.70 0.21–33.97 0.443 0.45 1.56 0.97–2.61 0.067#
nFW -Ch4 −0.52 0.60 0.20–1.82 0.365 0.076 1.08 0.23–5.12 0.924 −0.75 0.47 0.22–1.02 0.055#
nMD-t-Ch4 −0.15 0.86 0.27–2.72 0.801 −1.95 0.14 0.02–1.20 0.073 0.76 2.17 1.07–4.41 0.032 #

Ch4 FW imaging data as a predictor of Parkinson’s disease subtype stability at the 4-year follow-up.

TD, tremor-dominant subtype; ESS, Epworth Sleepiness Scale; REM, Rapid Eye Movement scale; Ch4, cholinergic nucleus 4; nFA-t-Ch4, normalized free-water-corrected fractional anisotropy in Ch4; nMD-t-Ch4:normalized free-water-corrected mean diffusivity in Ch4; FW, free-water value for Ch4; H-Y, Hoehn & Yahr.

# the value of H-Y, nFW -Ch4 and nMD-t-Ch4 were magnified by a factor of 10. Bold font indicates significant differences (p < 0.05).

Figure 4

Figure 4

Receiver operating characteristic curves of FW imaging data for cholinergic nucleus 4 in relation to subtype stability over 4 years in the initial tremor-dominant subtype group.

Discussion

In this study, we used FW imaging to investigate the microstructural changes of the two main cholinergic nuclei in different PD subtypes and determine their associations with motor function. Our findings revealed that FW-Ch4 was associated with the tremor subscore, while FA-t-Ch4 was negatively associated with the rigidity subscore. Moreover, nFW-Ch4 showed discriminative ability for predicting the stability of PD subtypes over 4 years. However, we did not observe any relationship between the volume of the cholinergic nuclei and motor symptoms, indicating that volume changes were not significant in the early stage and that the microstructural diffusion alterations measured by FW imaging may be more sensitive markers. These findings also indicated that in comparison with PPN, Ch4 is more associated with the motor score and motor subtypes.

The cholinergic nuclei NBM and PPN are at the epicenter of the cholinergic system, with most cholinergic projections in the brain originating from them and connecting with the brainstem, thalamus, and cortex (29). The cholinergic system has been associated with motor and NMS in PD (12, 13). Because of the interplay between cognitive and gait functions (30), previous studies mainly focused on the association of the cholinergic system with gait problems, freezing of gait, falls, and postural control (31, 32). DTI studies indicated that baseline increased axial diffusivity of the PPN detected with DTI was independently associated with the development of PIGD symptoms (33), and that FW in the cholinergic basal forebrain was correlated with a worse PIGD score, which may be mediated by attentional function (21). In addition to the cholinergic nuclei, the integrity of the cholinergic pathway also plays a role in motor symptoms (21). However, our study did not identify a significant relationship between the FW imaging data for the cholinergic nucleus and PIGD severity, and our findings could not characterize the role of the PPN. This may have occurred since the patients in the present study had early-stage PD. Thus, follow-up evaluation of the changes in cholinergic FW imaging may provide a better understanding of this aspect.

Surprisingly, we found that increased FW was positively correlated with the tremor subscore, which indicated that tremor severity was associated with microchanges in the Ch4 in the early de novo stage of PD. Tremor is one of the most prevalent symptoms in PD, but the pathophysiology of parkinsonian tremors remains poorly understood (34). Postmortem studies have demonstrated that PD patients with the TD subtype had less cell loss in the dopaminergic and nondopaminergic neural circuits than those without tremors (35). Moreover, the severity of parkinsonian tremors has been reported to be unrelated to the degree of dopaminergic denervation (36), consistent with our findings showing that the TD group had higher levels of dopamine in the caudate nucleus than the non-TD group. A growing body of evidence suggests that parkinsonian tremors may be related to degeneration of non-dopaminergic systems such as the serotoninergic and noradrenergic systems (34). Few studies have linked the cholinergic system to the severity of tremors. In fact elevated cholinergic activity is an important pathological marker in early PD, and anticholinergic therapy has shown remarkable effects on tremor symptoms. One study also indicated that cholinergic relevant functional reactivity is associated with the dopamine responsiveness of tremors in PD (37). In our study, an increased FW indicated degeneration or inflammation of Ch4, although the correlation between the FW-Ch4 alterations and acetylcholine release was not clear. Nevertheless, these findings provide the basis for further research on the role of cholinergic system damage in PD tremors.

Rigidity is another cardinal motor sign in PD, but it has been evaluated in only a relatively small number of clinical and experimental studies. An overnight sleep study showed higher and more symmetric upper limb rigidity indexes during wakefulness in patients with PD and REM sleep without atonia (RSWA) than in patients without RSWA and healthy controls (38). Although no correlation was observed between the rigidity subscore and the REM score in our sample, the findings could not clarify whether these scores are connected by a cholinergic mechanism.

Many studies on motor subtypes have been sourced from the PPMI database (2, 6, 7, 39), and the differences in other indicators, such as CSF biomarkers, among different subtypes were not shown in this study. In this study, the majority of the patients belonged to the TD subtype at baseline, while 35.71% of the individuals showed a shift in the motor subtypes at the fourth-year visit. Fluctuations in PD motor subtypes are common (7, 40, 41), especially during the first 2 years (42). These subtype shifts may be related to the progression of PIGD (40, 41), the effects of NMS (14), and drug interventions (41). In this regard, “true” TD patients are thought to show better subtype stability over time (40). The importance of non-motor manifestations as drivers of PD subtyping and prognosis has been proposed. The type composite with critical non-motor features (named as “diffuse malignant” type) demonstrates a more profound dopaminergic deficit, increased atrophy in PD brain networks, a more Alzheimer’s disease-like CSF profile, and faster progression of motor and cognitive deficits (2). The pathophysiology of these NMS, especially cognitive impairment and rapid eye movement sleep behavior disorder, are related to alterations of the cholinergic systems (12), indicating the role of the cholinergic system in the motor subtype. Meanwhile, as the disease progresses, resulting in a more severe dopaminergic deficit, the effect of the cholinergic system on subtype classification may increase. In this study, TD patients had higher putamen dopamine transporter availability than non-TD patients at baseline, but not at the 4-year follow-up. Moreover, the asymmetry in striatal denervation indices between the TD and PIGD subtypes has been reported to be the most prominent at the baseline, but this difference became non-significant after the 4-year follow-up (43). As the effects of dopamine reduce, the effects of the cholinergic system may become more apparent.

Patients showing PD with PIGD often face a more malignant prognosis than patients with TD. A model capable of identifying potential PIGD converters from the TD group would enable more personalized and potentially more effective therapeutic interventions. In our studies, the initial nFW-Ch4 and nMD-t-Ch4 showed predictive value for the 4-year motor subtype in the TD subgroup. FW indicates degeneration or inflammation, and MD is associated with alterations due to changes in membrane integrity. These changes may help us understand the pathological processes underlying the motor subtype in PD. Further studies in larger populations are required to validate these findings.

Limitations

This study had some limitations that require consideration. First, the sample size of this study was relatively small, and some cases were lost during follow-up. Second, we did not consider NMS in our subtyping algorithm and chose the traditional scoring method for grouping motor subtypes, which may have more widespread application to clinical care. However, PD motor subtypes determined by different established algorithms are inconsistent and unstable over time (40), and different categorizations may yield different outcomes. Third, to reduce the interference of medication, we chose to evaluate patients in the off state. However, in the PPMI study protocol, the off state only refers to overnight withdrawal of the patients’ anti-parkinsonian medications. Consequently, the drugs may not really have been metabolized to completion. Fourth, the majority of DTI data in the PPMI database were acquired using a single-shell protocol, which generally yields lower-quality FW analysis than multiple-shell DTI (17). Moreover, PPMI is a multicenter database, and we did not distinguish or correct the data from different centers. Fifth, we did not investigate the changes in cholinergic FW imaging data over time, which can precisely reflect the mechanisms underlying the disease-induced changes. Moreover, FW imaging only provided an indirect measure, and could not directly exhibit the neurotransmitter-related alterations in brain.

Conclusion

Our findings indicate that motor symptoms in PD were correlated with FW imaging data for the cholinergic nucleus and that nFW-Ch4 can be used to predict the transition of motor subtypes over a 4-year follow-up period, especially in patients initially classified into the TD group. Thus, our study provides valuable insights regarding the association of motor subtypes with the cholinergic system.

Statements

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. 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

SX: Conceptualization, Funding acquisition, Methodology, Data curation, Investigation, Writing – original draft. XS: Conceptualization, Writing – original draft, Formal analysis. MC: Investigation, Methodology, Writing – original draft. FF: Investigation, Writing – original draft. JT: Formal analysis, Visualization, Writing – original draft. BZ: Conceptualization, Funding acquisition, Supervision, Writing – review & editing. XL: Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by grants from the National Key Research and Development Program of China (No. 2022YFE0210100), the National Natural Science Foundation of China (No. 82271268), National Health Commission Scientific Research Fund -Zhejiang Provincial Major Health Science and Technology Plan (No. WKJ-ZJ-2218) and Program of Zhejiang Provincial TCM Sci-tech Plan (2024ZL235).

Acknowledgments

The authors are grateful to all PPMI colleagues and patients who participated in this study.

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.

References

  • 1.

    Sieber BA Landis S Koroshetz W Bateman R Siderowf A Galpern WR et al . Prioritized research recommendations from the National Institute of Neurological Disorders and Stroke Parkinson's disease 2014 conference. Ann Neurol. (2014) 76:46972. doi: 10.1002/ana.24261

  • 2.

    Fereshtehnejad SM Zeighami Y Dagher A Postuma RB . Clinical criteria for subtyping Parkinson's disease: biomarkers and longitudinal progression. Brain. (2017) 140:195976. doi: 10.1093/brain/awx118

  • 3.

    De Pablo-Fernandez E Lees AJ Holton JL Warner TT . Prognosis and Neuropathologic correlation of clinical subtypes of Parkinson disease. JAMA Neurol. (2019) 76:4709. doi: 10.1001/jamaneurol.2018.4377

  • 4.

    Thenganatt MA Jankovic J . Parkinson disease subtypes. JAMA Neurol. (2014) 71:499504. doi: 10.1001/jamaneurol.2013.6233

  • 5.

    Stebbins GT Goetz CG Burn DJ Jankovic J Khoo TK Tilley BC . 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. (2013) 28:66870. doi: 10.1002/mds.25383

  • 6.

    Choi SM Kim BC Cho BH Kang KW Choi KH Kim JT et al . Comparison of two motor subtype classifications in de novo Parkinson's disease. Parkinsonism Relat Disord. (2018) 54:748. doi: 10.1016/j.parkreldis.2018.04.021

  • 7.

    Simuni T Caspell-Garcia C Coffey C Lasch S Tanner C Marek K et al . How stable are Parkinson's disease subtypes in de novo patients: analysis of the PPMI cohort?Parkinsonism Relat Disord. (2016) 28:627. doi: 10.1016/j.parkreldis.2016.04.027

  • 8.

    Armstrong MJ Okun MS . Diagnosis and treatment of Parkinson disease: a review. JAMA. (2020) 323:54860. doi: 10.1001/jama.2019.22360

  • 9.

    Alster P Madetko-Alster N . Significance of dysautonomia in Parkinson's disease and atypical parkinsonisms. Neurol Neurochir Pol. (2024) 58:1479. doi: 10.5603/pjnns.98678

  • 10.

    Antonini A Barone P Marconi R Morgante L Zappulla S Pontieri FE et al . The progression of non-motor symptoms in Parkinson's disease and their contribution to motor disability and quality of life. J Neurol. (2012) 259:262131. doi: 10.1007/s00415-012-6557-8

  • 11.

    Kang SH Kim J Lee J Koh SB . Mild cognitive impairment is associated with poor gait performance in patients with Parkinson's disease. Front Aging Neurosci. (2022) 14:1003595. doi: 10.3389/fnagi.2022.1003595

  • 12.

    Pasquini J Brooks DJ Pavese N . The cholinergic brain in Parkinson's disease. Mov Disord Clin Pract. (2021) 8:101226. doi: 10.1002/mdc3.13319

  • 13.

    Perez-Lloret S Barrantes FJ . Deficits in cholinergic neurotransmission and their clinical correlates in Parkinson's disease. NPJ Parkinsons Dis. (2016) 2:16001. doi: 10.1038/npjparkd.2016.1

  • 14.

    Ren J Pan C Li Y Li L Hua P Xu L et al . Consistency and stability of motor subtype classifications in patients with de novo Parkinson's disease. Front Neurosci. (2021) 15:637896. doi: 10.3389/fnins.2021.637896

  • 15.

    Fereshtehnejad SM Romenets SR Anang JB Latreille V Gagnon JF Postuma RB . New clinical subtypes of Parkinson disease and their longitudinal progression: a prospective cohort comparison with other phenotypes. JAMA Neurol. (2015) 72:86373. doi: 10.1001/jamaneurol.2015.0703

  • 16.

    Mitchell T Lehericy S Chiu SY Strafella AP Stoessl AJ Vaillancourt DE . Emerging neuroimaging biomarkers across disease stage in Parkinson disease: a review. JAMA Neurol. (2021) 78:126272. doi: 10.1001/jamaneurol.2021.1312

  • 17.

    Nemmi F Levardon M Peran P . Brain-age estimation accuracy is significantly increased using multishell free-water reconstruction. Hum Brain Mapp. (2022) 43:236576. doi: 10.1002/hbm.25792

  • 18.

    Hoy AR Kecskemeti SR Alexander AL . Free water elimination diffusion tractography: a comparison with conventional and fluid-attenuated inversion recovery, diffusion tensor imaging acquisitions. J Magn Reson Imaging. (2015) 42:157281. doi: 10.1002/jmri.24925

  • 19.

    Ray NJ Lawson RA Martin SL Sigurdsson HP Wilson J Galna B et al . Free-water imaging of the cholinergic basal forebrain and pedunculopontine nucleus in Parkinson's disease. Brain. (2023) 146:105364. doi: 10.1093/brain/awac127

  • 20.

    Schumacher J Ray NJ Hamilton CA Bergamino M Donaghy PC Firbank M et al . Free water imaging of the cholinergic system in dementia with Lewy bodies and Alzheimer's disease. Alzheimers Dement. (2023) 19:454963. doi: 10.1002/alz.13034

  • 21.

    Wu C Wu H Zhou C Guan X Guo T Cao Z et al . Cholinergic basal forebrain system degeneration underlies postural instability/gait difficulty and attention impairment in Parkinson's disease. Eur J Neurol. (2024) 31:e16108. doi: 10.1111/ene.16108

  • 22.

    Jost ST Kaldenbach MA Antonini A Martinez-Martin P Timmermann L Odin P et al . Levodopa dose equivalency in Parkinson's disease: updated systematic review and proposals. Mov Disord. (2023) 38:123652. doi: 10.1002/mds.29410

  • 23.

    Jbabdi S Sotiropoulos SN Savio AM Graña M Behrens TE . Model-based analysis of multishell diffusion MR data for tractography: how to get over fitting problems. Magn Reson Med. (2012) 68:184655. doi: 10.1002/mrm.24204

  • 24.

    Pasternak O Sochen N Gur Y Intrator N Assaf Y . Free water elimination and mapping from diffusion MRI. Magn Reson Med. (2009) 62:71730. doi: 10.1002/mrm.22055

  • 25.

    Garyfallidis E Brett M Amirbekian B Rokem A van der Walt S Descoteaux M et al . Dipy, a library for the analysis of diffusion MRI data. Front Neuroinform. (2014) 8:8. doi: 10.3389/fninf.2014.00008

  • 26.

    Kilimann I Grothe M Heinsen H Alho EJ Grinberg L Amaro E Jr et al . Subregional basal forebrain atrophy in Alzheimer's disease: a multicenter study. J Alzheimers Dis. (2014) 40:687700. doi: 10.3233/JAD-132345

  • 27.

    Li J Guan X Wu Q He C Zhang W Lin X et al . Direct localization and delineation of human pedunculopontine nucleus based on a self-supervised magnetic resonance image super-resolution method. Hum Brain Mapp. (2023) 44:378194. doi: 10.1002/hbm.26311

  • 28.

    Puonti O Iglesias JE Van Leemput K . Fast and sequence-adaptive whole-brain segmentation using parametric Bayesian modeling. NeuroImage. (2016) 143:23549. doi: 10.1016/j.neuroimage.2016.09.011

  • 29.

    Selden NR Gitelman DR Salamon-Murayama N Parrish TB Mesulam MM . Trajectories of cholinergic pathways within the cerebral hemispheres of the human brain. Brain. (1998) 121:224957. doi: 10.1093/brain/121.12.2249

  • 30.

    Monaghan AS Gordon E Graham L Hughes E Peterson DS Morris R . Cognition and freezing of gait in Parkinson's disease: a systematic review and meta-analysis. Neurosci Biobehav Rev. (2023) 147:105068. doi: 10.1016/j.neubiorev.2023.105068

  • 31.

    Bohnen NI Kanel P Zhou Z Koeppe RA Frey KA Dauer WT et al . Cholinergic system changes of falls and freezing of gait in Parkinson's disease. Ann Neurol. (2019) 85:53849. doi: 10.1002/ana.25430

  • 32.

    Roytman S Paalanen R Griggs A David S Pongmala C Koeppe RA et al . Cholinergic system correlates of postural control changes in Parkinson's disease freezers. Brain. (2023) 146:324357. doi: 10.1093/brain/awad134

  • 33.

    Craig CE Jenkinson NJ Brittain JS Grothe MJ Rochester L Silverdale M et al . Pedunculopontine nucleus microstructure predicts postural and gait symptoms in Parkinson's disease. Mov Disord. (2020) 35:1199207. doi: 10.1002/mds.28051

  • 34.

    Dirkx MF Bologna M . The pathophysiology of Parkinson's disease tremor. J Neurol Sci. (2022) 435:120196. doi: 10.1016/j.jns.2022.120196

  • 35.

    Paulus W Jellinger K . The neuropathologic basis of different clinical subgroups of Parkinson's disease. J Neuropathol Exp Neurol. (1991) 50:74355. doi: 10.1097/00005072-199111000-00006

  • 36.

    Pirker W . Correlation of dopamine transporter imaging with parkinsonian motor handicap: how close is it?Mov Disord. (2003) 18:S4351. doi: 10.1002/mds.10579

  • 37.

    Wu JJ Zhou C Guo T Guan XJ Gao T Bai XQ et al . Cholinergic relevant functional reactivity is associated with dopamine responsiveness of tremor in Parkinson's disease. Brain Imaging Behav. (2022) 16:123445. doi: 10.1007/s11682-021-00610-9

  • 38.

    Linn-Evans ME Petrucci MN Huffmaster SLA Chung JW Tuite PJ Howell MJ et al . REM sleep without atonia is associated with increased rigidity in patients with mild to moderate Parkinson's disease. Clin Neurophysiol. (2020) 131:200816. doi: 10.1016/j.clinph.2020.04.017

  • 39.

    Jeong EH Lee JY Han SK Song YS . Motor subtypes and clinical characteristics in sporadic and genetic Parkinson's disease groups: analysis of the PPMI cohort. Front Neurol. (2023) 14:1276251. doi: 10.3389/fneur.2023.1276251

  • 40.

    von Coelln R Gruber-Baldini AL Reich SG Armstrong MJ Savitt JM Shulman LM . The inconsistency and instability of Parkinson's disease motor subtypes. Parkinsonism Relat Disord. (2021) 88:138. doi: 10.1016/j.parkreldis.2021.05.016

  • 41.

    Kohat AK Ng SYE Wong ASY Chia NSY Choi X Heng DL et al . Stability of MDS-UPDRS motor subtypes over three years in early Parkinson's disease. Front Neurol. (2021) 12:704906. doi: 10.3389/fneur.2021.704906

  • 42.

    Erro R Picillo M Amboni M Savastano R Scannapieco S Cuoco S et al . Comparing postural instability and gait disorder and akinetic-rigid subtyping of Parkinson disease and their stability over time. Eur J Neurol. (2019) 26:12128. doi: 10.1111/ene.13968

  • 43.

    Jeong EH Sunwoo MK Lee JY Han SK Hyung SW Song YS . Serial changes of I-123 FP-CIT SPECT binding asymmetry in Parkinson's disease: analysis of the PPMI data. Front Neurol. (2022) 13:976101. doi: 10.3389/fneur.2022.976101

Summary

Keywords

Parkinson’s disease, motor subtypes, free-water imaging, cholinergic nucleus 4, pedunculopontine nucleus

Citation

Xu S, Si X, Cai M, Fu F, Tian J, Zhang B and Liu X (2025) Association of free-water imaging data for the cholinergic nucleus with the motor function and subtypes in Parkinson’s disease. Front. Neurol. 16:1477827. doi: 10.3389/fneur.2025.1477827

Received

08 August 2024

Accepted

10 March 2025

Published

04 April 2025

Volume

16 - 2025

Edited by

Xinglong Yang, The First Affiliated Hospital of Kunming Medical University, China

Reviewed by

Piotr Alster, Medical University of Warsaw, Poland

Luoyu Wang, Hangzhou First People’s Hospital, China

Updates

Copyright

*Correspondence: Baorong Zhang, ; Xiaoli Liu

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics