Parkinson's disease (PD) is a neurodegenerative disorder characterized by a wide array of motor and non-motor symptoms. It remains unclear whether neurodegeneration in discrete loci gives rise to discrete symptoms, or whether network-wide atrophy gives rise to the unique behavioural and clinical profile associated with PD. Here we apply a data-driven strategy to isolate largescale, multivariate associations between distributed atrophy patterns and clinical phenotypes in PD. In a sample of N = 229 de novo PD patients, we estimate disease-related atrophy using deformation based morphometry (DBM) of T1 weighted MR images. Using partial least squares (PLS), we identify a network of subcortical and cortical regions whose collective atrophy is associated with a clinical phenotype encompassing motor and non-motor features. Despite the relatively early stage of the disease in the sample, the atrophy pattern encompassed lower brainstem, substantia nigra, basal ganglia and cortical areas, consistent with the Braak hypothesis. In addition, individual variation in this putative atrophy network predicted longitudinal clinical progression in both motor and non-motor symptoms. Altogether, these results demonstrate a pleiotropic mapping between neurodegeneration and the clinical manifestations of PD, and that this mapping can be detected even in de novo patients.