BackgroundBipolar disorder I (BD-I) is a heterogeneous disorder with a high prevalence of comorbid anxiety. The aim of this study was to investigate whether anxiety and mania symptoms define distinct subgroups within BD-I and to explore potential differences in functional network characteristics between these subgroups. MethodsSubgroups were identified using scores from clinical anxiety and mania scales. After dimension reduction of these scores, data-driven clustering analysis with cross-validation was employed to reveal the existence of subgroups. Resting-state functional magnetic resonance imaging (rs-fMRI) scans were pre-processed using fMRIPrep. After parcellation and network construction, global and regional graph theoretical measures were calculated per subgroup. ResultsClustering results revealed that, based on anxiety symptomatology, subjects fell into two distinct subgroups, whereas mania symptoms divided subjects into four unique subgroups. These subgroups varied notably on several symptom scales. Network assortativity was significantly associated with anxiety subgroups. Post-hoc pairwise comparisons did not reveal significant global functional network differences between the anxiety subgroups or between mania subgroups. Regional network differences between clinical subgroups were especially apparent for strength and degree in the temporal and frontal lobes. LimitationsSmall sample size of some subgroups is a limitation of this study as is the categorical rather than continuous representation of anxiety and mania symptoms. ConclusionsBD-I populations may be stratified into robust subgroups based on anxiety and mania symptoms, showing differences in functional network connectivity. Our findings highlight new avenues of research for investigating heterogeneity in psychiatric populations.