BackgroundParkinson's disease (PD) is a heterogeneous condition. Cluster analysis based on cortical thickness has been used to define distinct patterns of brain atrophy in PD. However, the potential of other neuroimaging modalities, such as white matter (WM) fractional anisotropy (FA), which has also been demonstrated to be altered in PD, has not been investigated. ObjectiveWe aim to characterize PD subtypes using a multimodal clustering approach based on cortical and subcortical gray matter (GM) volumes and FA measures. MethodsWe included T1-weighted and diffusion-weighted MRI data from 62 PD patients and 33 healthy controls. We extracted mean GM volumes from 48 cortical and 17 subcortical regions using FSL-VBM, and the mean FA from 20 WM tracts using Tract-Based Spatial Statistics (TBSS). Hierarchical cluster analysis was performed with the PD sample using Ward's linkage method. Whole-brain voxel-wise intergroup comparisons of VBM and TBSS data were also performed using FSL. Neuropsychological and demographic statistical analyses were conducted using IBM SPSS Statistics 25.0. ResultsWe identified three PD subtypes, with prominent differences in GM patterns and little WM involvement. One group (n = 15) with widespread cortical and subcortical GM volume and WM FA reductions and pronounced cognitive deficits; a second group (n = 21) with only cortical atrophy limited to frontal and temporal regions and more specific neuropsychological impairment, and a third group (n = 26) without detectable atrophy or cognition impairment. ConclusionMultimodal MRI data allows classifying PD patients into groups according to GM and WM patterns, which in turn are associated with the cognitive profile.