Twenty-two PD-TD and 21 PD-PIGD patients were recruited from the Movement Disorders Clinic at the University Medical Center Hamburg-Eppendorf between 2021 and 2024, along with 21 age- and sex-matched healthy controls (HC). All participants met the following inclusion criteria: (1) age above 18 years, and (2) written informed consent to participate in the study. Specific inclusion criteria for the patients included diagnosis of PD in accordance with the United Kingdom Parkinson's Disease Society Brain Bank (UKPDSBB) (Gibb and Lees, 1988) by a board-certified neurologist specialized in movement disorders.

Exclusion criteria applied uniformly across all participants encompassed: (1) coexistence of other neurologic disorders, including secondary or atypical PD, (2) history of deep brain stimulation or other significant head surgery, (3) electronic or metal implants, such as cardiac pacemakers and aneurysm clips, (4) pregnancy, (5) intake of other centrally acting medication such as antidepressants, (6) incompatibility with MRI safety standards, and (7) inability to provide informed consent.

Standardized motor symptom evaluation of PD patients was performed based on the Movement MDS-UPDRS parts II and III in a single session (Goetz et al., 2008). The examinations were carried out by movement disorder-trained medical personnel and videos were later reviewed by a board-certified neurologist who is specialized in movement disorders as quality control. Both the examiner and the reviewing physician were blinded to the MRI data. The ratio of MDS-UPDRS subscores for tremor severity to those for PIGD was computed and used to determine clinical subtypes of PD: PD-TD and PD-PIGD (Stebbins et al., 2013). Moreover, we focused on six key symptom domains derived from MDS-UPDRS Part III: rest tremor (items 3.17–3.18), postural tremor (item 3.15), kinetic tremor (item 3.16), limb bradykinesia (items 3.4–3.8), PIGD (items 3.9–3.12), and the total motor score (MDS-UPDRS-III), which reflects the sum of all Part III items.

All patients were evaluated in the “OFF” state i.e. not having taken medication for at least 12 h prior to the experiment. To assess the influence of the daily intake of antiparkinsonian drugs the levodopa equivalent daily dose (LEDD) was calculated for each patient (Jost et al., 2023). Additionally, the handedness of participants was determined by the Edinburgh Handedness Inventory (Oldfield, 1971). The staging of PD functional disability was further identified through the HY staging (Hoehn and Yahr, 1967). This study has been approved by the local ethics committee of Hamburg (2020-10281_BO-ff). All participants provided written informed consent, and the study was performed in accordance with the ethical standards as laid down in the Declaration of Helsinki.

All statistical analyses were conducted using R (version 4.3.2) and Python (version 3.8.12). A significance threshold of p < 0.05 was applied throughout. Where applicable, p-values were adjusted for multiple comparisons using the false discovery rate (FDR) correction (Benjamini and Hochberg, 1995). For classification algorithms and models associating volume with symptom severity, leave-one-out cross-validation (LOOCV) was used to estimate p-value stability, and median p-values across iterations were reported.

Two patients were excluded due to the insufficient quality of the acquired MRI images for processing. Additionally, one patient and two HC were excluded due to coincidental abnormal MRI findings. Table 1 summarizes the demographic information of the final study cohort. The HC group had a mean age of 63 ± 9 years (58% females); the PD-PIGD cohort 68 ± 8 years (female ratio: 39%), and the PD-TD cohort 61 ± 9 years (female ratio: 41%). No significant age difference was observed between HC and PD groups (t-test, p = 0.838), nor was there a difference in sex distribution (Fisher's exact test, p = 0.266). However, a significant age difference was found between PD motor subtypes, with PD-TD patients being younger than PD-PIGD patients (p = 0.010); we therefore consistently included age as a covariate in all analyses and interpreted age-sensitive findings with caution. In addition, we performed an age-matched sensitivity analysis for PD-TD vs. PD-PIGD single-ROI comparisons, which showed preserved effect directions and no contradictory findings relative to the main analysis (see Supplementary Table S1). No significant differences were observed between PD-TD and PD-PIGD groups in sex distribution (p = 1.000) or HY stage (p = 0.496). LEDD values did not differ significantly between PD-TD and PD-PIGD groups (Mann–Whitney U = −92, p = 0.321), suggesting that medication dosage is not a major confounding factor for the observed group differences.

All automated and manual segmentations were successfully completed across participants without major artifacts or processing failures. The extracted regional volumes served as the basis for subsequent statistical analyses. A representative example of the anatomical segmentation results is shown in Figure 1, including both cortical and subcortical regions from the BTC and CTC networks. Levene's tests revealed no significant violations of variance homogeneity across cohorts for any of the volumetric measures (all p > 0.05). Similarly, Shapiro–Wilk tests indicated that most ROI volumes conformed to normality (p > 0.05), with only mild deviations in certain smaller regions such as lobule IX, which upon closer inspection were deemed acceptable for parametric analysis. Detailed plots are provided in the Supplementary Figure S1.

We observed significant anatomical asymmetry in a number of regions. Rightward asymmetry was observed in the cerebellum (HC: AI_anatomical = 0.006; PD-TD: 0.009; PD-PIGD: 0.005; all FDR p < 0.001), Crus I (HC: AI_anatomical= 0.011; PD-TD: 0.017; PD-PIGD: 0.022; all FDR p < 0.01), Crus II (HC: AI_anatomical= 0.019; PD-TD: 0.031; PD-PIGD: 0.019; all FDR p < 0.01), and VLp (PD-TD: AI_anatomical= 0.036, FDR p < 0.001), with similar patterns across groups. In contrast, leftward asymmetry was found in the DN (PD-PIGD: −0.084, FDR p < 0.001), RN (PD-TD: −0.044, FDR p < 0.01), thalamus (PD-TD: −0.026, FDR p < 0.001), and VLa (PD-TD: −0.018, FDR p < 0.001). The remaining ROIs showed no significant asymmetry (all FDR p > 0.05).

The ANCOVA analyses (controlled for age, sex, HY stage, and handedness) showed that these asymmetry patterns were not significantly different between cohorts (all FDR p > 0.05), suggesting cohort-independent anatomical asymmetry. Furthermore, no significant effects of handedness were observed across regions (all FDR p > 0.05).

One-sample t-tests revealed that none of the clinical asymmetry indices (AI_clinical) significantly deviated from zero (all FDR p > 0.05), indicating the absence of systematic lateralization of volumes across the motor networks due to symptom expression. Furthermore, no significant alignment was observed between anatomical and clinical asymmetry, as confirmed by Kendall correlation analyses across all ROIs (all FDR p > 0.05). Group comparisons between PD subtypes revealed no significant differences in clinical asymmetry across any region (all FDR p > 0.05). Importantly, no significant associations were found between anatomical and clinical asymmetry (all FDR p > 0.05), or between asymmetry indices and handedness (all FDR p > 0.05). Detailed plots are provided in the Supplementary Figures S2, S3.

Therefore, for all subsequent analyses involving patient data, MAS ROIs were included based on clinical lateralization; while for HC, the average of both hemispheres was used. Anatomical asymmetry and handedness were not considered further, as neither was found to systematically influence the results.

No significant group differences in MAS volume were observed between PD patients and HC in any region (all FDR p > 0.05). These results remained consistent after outlier removal. As shown in Figure 2, volume distributions overlapped substantially across groups, reinforcing the absence of robust group-level effects.

No multicollinearity was detected among the ROI volumes, as verified by pairwise correlation analysis (see Supplementary Figure S4). Ridge regression and PCA yielded highly comparable classification performance across all tasks, with similar mean AUCs (ridge: 0.63 ± 0.13, 95 % CI = 0.62–0.64; PCA: 0.64 ± 0.13, 95 % CI = 0.63–0.65). Given this parity in predictive accuracy, we prioritized ridge regression for further interpretation due to its enhanced model transparency and more consistent feature importance estimates across splits and iterations.

Using z-scored MAS volumes, classification of PD vs. HC reached an average AUC of 0.88 and accuracy of 0.86 among the best-performing iterations, while the overall mean across 500 iterations was 0.63 ± 0.13 (95 % CI = 0.62–0.64). Feature-wise AUCs and importance scores are shown in Figure 3.

Feature importance analysis revealed that the strongest positive contributors to classification included the thalamus, VLa, VLp, and cerebellar lobules I–V, VIIIa, and X. These regions showed positive weights in the model, suggesting that preserved or increased volume in these areas supports PD classification. Conversely, regions such as M1, DN, RN, cerebellum (as a whole), and lobules VI, Crus I-II, and IX contributed negatively, indicating that their volume reduction is informative of PD. To complement importance values, we also computed individual-region AUCs. While some regions had small model weights, they still demonstrated strong standalone discriminative power—most notably lobule VIIIa (AUC = 0.83), VLp (AUC = 0.79), I–II (AUC = 0.79), and thalamus (AUC = 0.75). Several negatively weighted regions like M1 (AUC = 0.71) and VLa (AUC = 0.83) also achieved high AUCs individually, reinforcing their relevance in PD-related volume reduction.

Covariates showed moderate predictive value. Age was weakly informative (AUC = 0.75, importance near zero). Sex showed a small but consistent directional contribution (AUC = 0.58), with male sex slightly associated with PD classification.

To perform motor subtype classification, the ridge model was trained to identify PD-TD participants using z-scored MAS volumes, with positive contributions indicating regional preservation predictive of PD-TD and negative contributions indicating reduced volumes associated with PD-TD (Figure 4). The model achieved a mean AUC of 0.66 ± 0.15 (95 % CI = 0.65–0.68) across 500 iterations, with top-performing runs reaching an AUC of 0.95 and an accuracy of 0.72. PCA produced comparable results (mean AUC = 0.67 ± 0.14, 95 % CI = 0.65–0.68; top AUC = 0.92), confirming the robustness of the multivariate classification approach. The strongest positive contributors for PD-TD included lobules I–II, III, VIIb, VIIIa, and VLp. In contrast, M1, DN, thalamus, VLa, and SN–STN showed the most negative importance scores. Individual ROI classifiers (green dotted line in Figure 4) revealed the highest AUCs (≥0.83) for VLp, SN–STN, cerebellum, and lobules VI and VIIIa. HY stage also contributed positively (importance = 0.10), while age and sex showed negligible effects. LOOCV results remained consistent (AUC = 0.67, balanced accuracy = 50%).

To illustrate classification performance at the model level, receiver operating characteristic (ROC) curves summarizing sensitivity–specificity trade-offs across iterations are shown in Figure 5 for both PD vs. HC and PD subtype classification.

We examined how z-scored MAS volumes of the chosen regions relate to motor symptom severity using both multi-ROI and single-ROI linear regression models, adjusting for age, sex, and HY stage. In the multi-ROI models, total tremor was significantly associated with the volumes of the SN–STN, VLp, and lobule VIIb. Rest tremor was also associated with VLp, and postural tremor with M1, lobules I-II and IX. MDS-UPDRS-III total scores were significantly associated with volumes of the thalamus and the DN, while thalamus volume also showed a distinct association with bradykinesia. PIGD severity, in turn, was linked to reduced volume in Lobule IV.

The single-ROI models supported many of these findings and highlighted additional associations. For total tremor, significant contributions were again found for VLp and thalamus. Rest tremor remained associated with VLp, while kinetic tremor showed links to both VLp and lobule III of the cerebellum. PIGD severity was associated with VLa, and MDS-UPDRS-III score was associated with thalamus volume. Details of the significant results, including standardized beta coefficients, 95% confidence intervals, and LOOCV p-values, are summarized in Tables 2, 3 for multi-ROI and single-ROI models respectively.

This study establishes a new perspective on early-stage PD by showing that, in the present cohort, its structural substrate is distributed, not regional. Using an integrative multivariate framework that combines volumetric measures across BTC and CTC networks, we identified coherent network-wide structural alterations that differentiate PD from HC and further separate PD-TD from PD-PIGD motor subtypes. Unlike most prior volumetric studies that focused on isolated anatomical structures (Banwinkler et al., 2022; Charroud and Turella, 2021; Mak et al., 2015; Sterling et al., 2017; Torres-Parga et al., 2025), this approach uncovers distributed, multiregional volumetric patterns, emphasizing that in this sample, early PD-related alterations manifest at the network level rather than as pronounced focal atrophy (Caligiore et al., 2016).

On the regional level, these distributed signatures point to a continuum between relative volume preservation and early volume reduction, suggesting an imbalance between compensatory and degenerative processes within BTC–CTC pathways in PD. By shifting the focus from isolated regional atrophy to network-wide volumetry, our findings propose sensitive MRI biomarkers that reflect early disease structural network alterations and provide a basis for distinguishing PD subtypes.

The absence of significant regional atrophy in standard cohort comparisons contrasts with the distributed volumetric differences detected by our multivariate analysis. This apparent discrepancy may reflect an early, masked phase of PD, during which large-scale network alterations emerge before overt tissue loss becomes detectable with conventional morphometry. At this stage, compensatory remodeling within interconnected motor circuits may temporarily preserve function, while age-related atrophy patterns partially obscure subtle disease-specific effects (Johansson et al., 2025; Zhu et al., 2025). It has to be considered though that given the modest sample size and early disease stage, the null univariate findings in our cohort may reflect limited power to detect subtle regional effects rather than their true absence.

Consistent hemispheric asymmetries across cohorts, but without group differences or alignment with lateralization, suggest that motor volume asymmetry reflects general anatomical variability rather than a disease-specific marker. This pattern, together with our cohort being predominantly in the HY stage 2 (PD-TD: 20/22, PD-PIGD: 15/18), supports the notion that adaptive and degenerative processes coexist in early-stage PD. Several studies reflect the same interpretation. In the large-scale Parkinson's Progression Markers Initiative (PPMI) cohort, Zeighami et al. (2015) reported only subtle gray-matter alterations in de novo PD (n = 232 PD, 117 HC), largely overlapping with healthy-aging trajectories. Longitudinal data further demonstrate that macroscopic atrophy in basal ganglia and cerebellum progresses slowly compared with earlier functional and network-level disruptions (Zarkali et al., 2024). Similarly, aging cohorts show comparably small annual subcortical volume declines (n = 653, >10 years MRI; Fujita et al., 2023), suggesting that early-stage PD changes remain concealed within normative variance. Studies in newly diagnosed patients (Mak et al., 2015; Sterling et al., 2017) also report null volumetric alterations despite altered connectivity, with longitudinal evidence indicating that structural changes typically become detectable only after more than 5 years of disease progression.

Together, these observations highlight that early-stage PD represents a masked phase of anatomical reorganization, in which compensatory preservation and emerging degeneration coexist within distributed networks. Detecting this phase requires analytical approaches sensitive to interregional covariance, such as the multivariate modeling employed here, which can reveal coherent, subthreshold alterations invisible to conventional volumetric statistics.

Our multivariate classification results demonstrate that volumetric configurations of BTC and CTC regions can robustly differentiate PD from HC and further separate PD motor subtypes, even in the absence of regional atrophy. Using ridge regression, the model achieved high AUC values in the best-performing iterations, which confirmed that distributed volumetric patterns carry diagnostic information. However, the overall mean AUCs across all iterations were more modest, suggesting that the highest model's performance may be optimistic in some cases. Comparable multivariate approaches in large cohorts reported similar or slightly lower AUCs than in our analysis (Solana-Lavalle and Rosas-Romero, 2021; Ya et al., 2022). Although PCA yielded similar accuracy, ridge regression provided greater interpretability and stability in small-sample data (Schulz et al., 2020).

Feature-weight analyses revealed two opposing gradients across motor circuits. Thalamic nuclei (VLa, VLp) and cerebellar lobules I–V, VIIIa, and X were the strongest positive contributors to PD classification, consistent with early compensatory reorganization within the CTC network (Xiong et al., 2023). Preserved volume in these regions may reflect adaptive mechanisms that maintain network efficiency despite dopaminergic loss. Supporting this interpretation, functional connectivity studies have shown increased VLp–VLa coupling with motor and prefrontal cortices in PD, suggesting early network-level adaptations that precede overt atrophy (Chen Z. et al., 2023). Similarly, lobules VIIIa and X have been implicated in motor symptom modulation via altered connectivity with cortical motor areas (Xu et al., 2019). Moreover, multivariate structural connectivity analyses including cerebellar and thalamic nodes have demonstrated high discriminative power in PD (Owens-Walton et al., 2019), reinforcing the relevance of the same regions identified here. Individual ROI analyses further confirmed their importance, with lobule VIIIa, VLp, and thalamus achieving high AUCs, underscoring their standalone diagnostic potential within the distributed motor network.

In contrast, negatively weighted regions, including M1, DN, RN, Crus I–II, and lobule IX, showed reduced volumes that contributed to PD classification as markers of early neurodegeneration. These findings converge with prior evidence that subcortical and cerebellar atrophy plays a critical role in PD-related motor dysfunction (Camacho et al., 2024; Ya et al., 2022). DN atrophy, observed even in early disease stages (based on HY score), has been linked to disrupted cerebellar connectivity and impaired motor-cognitive integration (Liu et al., 2013; Pietracupa et al., 2024). Similarly, RN alterations, though variable across studies, have been associated with motor severity and progression (Camlidag et al., 2014), while reductions in Crus II and lobule IX volumes have been reported in both motor and non-motor phenotypes (Ma et al., 2018; Ya et al., 2022). Functional imaging further supports these structural findings, showing that diminished DN–Crus I–II coupling disrupts integration between motor and cognitive cerebellar loops (Liu et al., 2013). Cortical involvement followed a similar pattern, with M1 thinning and gray-matter loss reflecting degeneration along the corticospinal tract (Wang et al., 2024). Finally, posterior cerebellar atrophy has been shown to improve PD classification accuracy in multivariate frameworks (Zhao et al., 2025), emphasizing its diagnostic value. Collectively, these negatively weighted features delineate a coherent subcortical–cerebellar degeneration pattern that complements the compensatory CTC preservation identified in our positive-weight analysis.

In the classification of PD subtypes, our ridge regression model indicated a robust separation between PD-TD and PD-PIGD based on distributed volumetric features. Feature-weight analyses revealed that preserved volume in cerebellar lobules I–II, III, VIIb–VIIIa, and VLp, as well as reduced volume in M1, DN, thalamus, VLa, and SN–STN, both contributed to PD-TD classification, indicating that distinct patterns of preservation and loss within CTC and BTC circuits jointly characterize this subtype. These patterns align with prior neuroimaging and electrophysiological studies showing subtype-specific circuit organization. Chen Y. et al. (2023) and Chen Z. et al. (2023) reported preserved anterior cerebellar lobules and enhanced cerebello-thalamic connectivity in PD-TD, while PD-PIGD showed broader thalamocortical disruption. Similarly, Park et al. (2025) found selective atrophy in lobule VIIIa but preserved thalamic volume in PD-TD, underscoring the discriminative role of CTC network regions for subtype separation.

Subcortical degeneration also contributed to subtype classification. Functional connectivity studies have shown weakened DN–VLp coupling in PD-TD (Chen Z. et al., 2023), and volumetric reductions in M1 have been linked to disrupted cortical motor integration (Wang et al., 2023). Additionally, volumetric differentiation of thalamic nuclei across subtypes has been reported in de novo PD patients, with PD-TD showing reductions in VLa and related intralaminar nuclei (Park et al., 2020). This supports our findings where VLa volume contributed negatively to PD-TD classification. Moreover, histopathological and imaging evidence suggests that SN degeneration follows distinct topographies across subtypes: Boonstra et al. (2022) reported preferential medial SN neuron loss projecting to the VLa and STN, while Poston et al. (2020) linked SN volume reductions to bradykinesia and rigidity but not tremor. These findings underscore the anatomical specificity of dopaminergic degeneration in PD-TD and may explain the negative weighting of SN–STN and VLa in our model. Together, our classification findings reflect a biologically plausible dissociation: cerebellar preservation in PD-TD may reflect compensatory plasticity, while atrophy in the SN–STN–VLa–M1 axis captures a core pathophysiological substrate of PD-TD.

This distributed dissociation across CTC and BTC networks reinforces that motor subtypes emerge from distinct patterns of structural imbalance rather than uniform atrophy. Network-level analyses further support this view: graph-theoretical studies demonstrate early disruptions in global topology (Fang et al., 2017) and improved classification performance when using multiregional, network-level features instead of isolated regional measures (Gratton et al., 2019). Consistent with this systems-level organization, HY stage contributed moderate predictive value (high importance, low AUC), in line reflecting its role as a coarse clinical correlate that complements network-derived structural metrics (Kaplan et al., 2022).

The regression models, controlling for age, sex, and HY stage, revealed that distinct motor-network nodes are associated with specific symptom severities. Tremor severity was positively associated with SN–STN volume and negatively with VLp and lobule VIIb volumes, mirroring the evidence that both BTC and CTC networks are involved in tremor generation (Dirkx and Bologna, 2022; Helmich, 2018), and specifically that lobule VIIb atrophy correlates with tremor severity (Piccinin et al., 2017; Sadeghi et al., 2023). Furthermore, functional imaging has shown that decreased connectivity between the DN and VLp is linked to increased tremor severity in PD-TD, reinforcing the role of VLp in tremor mechanism (Chen Z. et al., 2023; Sadeghi et al., 2024). Rest tremor was specifically linked to VLp volume, while postural tremor engaged the M1 and lobules I–II and IX, consistent with studies showing M1 and anterior cerebellar activation during postural stabilization (Wang et al., 2023; Zeng et al., 2019). Moreover, PD-TD patients show increased M1–cerebellum functional connectivity through VLp, potentially reflecting compensatory engagement of motor networks (Bu et al., 2024; Chen Z. et al., 2023; Wang et al., 2023). This pattern aligns with evidence that CTC interactions—particularly through VLp—are more functionally intact or hyperactive in PD-TD relative to PD-PIGD (Sadeghi et al., 2024).

Bradykinesia and total MDS-UPDRS-III scores were associated with thalamus and DN volumes, supporting mechanistic models where altered basal ganglia–thalamic output to M1 underlies bradykinesia (Berardelli et al., 2001). Ultra-high field MRI confirms that SN degeneration is selectively associated with bradykinesia-rigidity but not tremor (Poston et al., 2020). Furthermore, a recent longitudinal MRI study demonstrated that multivariate gray matter volume metrics predicted progression in bradykinesia and gait symptoms over 5–10 years, underscoring their prognostic value (Vijayakumari et al., 2023).

PIGD severity correlated with lobule IV atrophy in both multi- and single-ROI models. Lobule IV, part of the spinocerebellum, maintains strong functional connections to vestibulocerebellar regions involved in balance and postural stabilization (Snijders et al., 2016; Wang et al., 2023). Consistent with our findings, PD-PIGD patients exhibit reduced white matter integrity in cerebellar tracts such as the middle cerebellar peduncle and diminished thalamic connectivity relative to PD-TD (Bu et al., 2024).

Together, these symptom-specific volumetric profiles highlight that Parkinsonian motor features emerge from network-wide structural alterations rather than localized atrophy (Caligiore et al., 2016), and suggest that targeted assessment of these network nodes could inform personalized prognosis and circuit-based interventions (Vijiaratnam et al., 2021). Our findings may have implications for clinical decision support. Although the average classification performance was overall modest, the use of multivariate volumetric patterns demonstrates proof of concept for clinical applications. The distinct volumetric signatures between PD-TD and PD-PIGD could inform stratification in clinical trials, addressing the heterogeneity that often confounds therapeutic studies (Greenland et al., 2019).

Several limitations should be considered when interpreting our findings. First, our cross-sectional design limits causal inferences about the relationship between volumetric changes and symptom progression. Longitudinal studies are needed to establish whether the identified volumetric signatures represent stable trait markers or evolve with disease progression. Second, the modest sample size, particularly for PD motor subtype comparisons in single ROI analysis, may have limited statistical power to detect subtle regional volumetric effects, and thus null univariate findings may not be interpreted as evidence for the absence of focal structural changes. This may explain the absence of significant regional atrophy in traditional group comparisons despite robust multivariate classification performance. Third, replication of these findings in larger and independent cohorts will be essential to confirm the generalizability and robustness of the observed network-level volumetric patterns. Fourth, our focus on motor network regions, while hypothesis-driven, may have overlooked contributions from other brain areas. Non-motor symptoms, which are prevalent even in early PD, were not consistently assessed, limiting our understanding of the full spectrum of brain-behavior relationships in our cohort. Finally, cognitive status and mood data were not available for analysis in this cohort, which limits the ability to examine their potential influence on volumetric measures and motor phenotype. These factors should be considered in future studies.

Despite these limitations, our multivariate approach successfully identified network-level volumetric signatures that robustly classified PD patients and revealed association with severity of motor symptoms, suggesting that these patterns capture biologically meaningful disease-related changes.

To conclude, this study demonstrates that multivariate modeling across PD motor network reveals meaningful, disease-relevant alterations even when regional volumetric differences may not be detectable in early-stage patient using traditional univariate analyses. By applying regularized classification techniques to ROI volumes within the BTC-CTC networks, we achieved robust discrimination of PD from HC and identified distinct volumetric profiles predictive of PD motor subtypes. These findings highlight the value of network-level structural biomarkers—particularly within the CTC circuit—as early and subtype-sensitive indicators of PD. Importantly, the convergence of our anatomical and symptom-association results highlights a structure–function relationship within these motor networks. Future studies are encouraged to expand on these findings by incorporating longitudinal data and multimodal imaging to track individualized disease progression and ultimately enable personalized, network-based diagnostics in PD.