Background/Aims: The topological organization of brain functional networks is impaired in Parkinson's disease (PD). However, the altered patterns of functional network hubs in different subtypes of PD are not completely understood. Methods: 3T resting-state functional MRI and voxel-based graph-theory analysis were employed to systematically investigate the intrinsic functional connectivity patterns of whole-brain networks. We enrolled 31 patients with PD (12 tremor dominant [TD] and 19 with postural instability/gait difficulty [PIGD]) and 22 matched healthy controls. Whole-brain voxel-wise functional networks were constructed by measuring the temporal correlations of each pair of brain voxels. Functional connectivity strength was calculated to explore the brain network hubs. Results: We found that both the TD and PIGD subtypes had comprehensive disrupted regions. These mainly involved the basal ganglia, cerebellum, superior temporal gyrus, pre- and postcentral gyri, inferior frontal gyrus, middle temporal gyrus, lingual gyrus, insula, and parahippocampal gyrus. Furthermore, the PIGD subgroup had more disrupted hubs in the cerebellum than the TD subgroup. These disruptions of hub connectivity were not correlated with the HY stage or disease duration. Conclusion: Our results emphasize the subtype-specific PD-related degeneration of brain hubs, providing novel insights into the pathophysiological mechanisms of connectivity dysfunction in different PD subgroups.