Gene expression plays a critical role in the pathogenesis of Parkinson's disease (PD). How gene expression profiles are correlated with functional-metabolic architecture remains obscure. We enrolled 34 PD patients and 25 age-and-sex-matched healthy controls for simultaneous F-18-FDG-PET/functional MRI scanning during resting state. We investigated the functional gradients and the ratio of standard uptake value. Principal component analysis was used to further combine the functional gradients and glucose metabolism into functional-metabolic architecture. Using partial least squares (PLS) regression, we introduced the transcriptomic data from the Allen Institute of Brain Sciences to identify gene expression patterns underlying the affected functional-metabolic architecture in PD. Between-group comparisons revealed significantly higher gradient variation in the visual, somatomotor, dorsal attention, frontoparietal, default mode, and subcortical network (p(FDR) < .048) in PD. Increased FDG-uptake was found in the somatomotor and ventral attention network while decreased FDG-uptake was found in the visual network (p(FDR) < .008). Spatial correlation analysis showed consistently affected patterns of functional gradients and metabolism (p = 2.47 x 10(-8)). PLS analysis and gene ontological analyses further revealed that genes were mainly enriched for metabolic, catabolic, cellular response to ions, and regulation of DNA transcription and RNA biosynthesis. In conclusion, our study provided genetic pathological mechanism to explain imaging-defined brain functional-metabolic architecture of PD.