AbstractDysexecutive Alzheimer’s disease (dAD) manifests as a progressive dysexecutive syndrome without prominent behavioral features, and previous studies suggest clinico-radiological heterogeneity within this syndrome. We uncovered this heterogeneity using unsupervised machine learning in 52 dAD patients with multimodal imaging and cognitive data. A spectral decomposition of covariance between FDG-PET images yielded six latent factors (“eigenbrains”) accounting for 48% of variance in patterns of hypometabolism. These eigenbrains differentially related to age at onset, clinical severity, and cognitive performance. A hierarchical clustering on the eigenvalues of these eigenbrains yielded four dAD subtypes, i.e. “left-dominant,” “right-dominant,” “bi-parietal-dominant,” and “heteromodal-diffuse.” Patterns of FDG-PET hypometabolism overlapped with those of tau-PET distribution and MRI neurodegeneration for each subtype, whereas patterns of amyloid deposition were similar across subtypes. Subtypes differed in age at onset and clinical severity where the heteromodal-diffuse exhibited a worse clinical picture, and the bi-parietal had a milder clinical presentation. We propose a conceptual framework of executive components based on the clinico-radiological associations observed in dAD. We demonstrate that patients with dAD, despite sharing core clinical features, are diagnosed with variability in their clinical and neuroimaging profiles. Our findings support the use of data-driven approaches to delineate brain–behavior relationships relevant to clinical practice and disease physiology.