Abstract To elucidate the brain‐wide information interactions that vary and contribute to individual differences in schizophrenia (SCZ), an information‐resolved method is employed to construct individual synergistic and redundant interaction matrices based on regional pairwise BOLD time‐series from 538 SCZ and 540 normal controls (NC). This analysis reveals a stable pattern of regionally‐specific synergy dysfunction in SCZ. Furthermore, a hierarchical Bayesian model is applied to deconstruct the patterns of whole‐brain synergy dysfunction into three latent factors that explain symptom heterogeneity in SCZ. Factor 1 exhibits a significant positive correlation with Positive and Negative Syndrome Scale (PANSS) positive scores, while factor 3 demonstrates significant negative correlations with PANSS negative and general scores. By integrating the neuroimaging data with normative gene expression information, this study identifies that each of these three factors corresponded to a subset of the SCZ risk gene set. Finally, by combining data from NeuroSynth and open molecular imaging sources, along with a spatially heterogeneous mean‐field model, this study delineates three SCZ synergy factors corresponding to distinct symptom profiles and implicating unique cognitive, neurodynamic, and neurobiological mechanisms.

Mapping Brain Synergy Dysfunction in Schizophrenia: Understanding Individual Differences and Underlying Molecular Mechanisms