BackgroundSchizophrenia (SZ) typically manifests heterogeneous phenotypes involving positive, negative and cognitive symptoms. However, the underlying neural mechanisms of these symptoms keep unclear. Functional gradient is a fascinating measure to characterize continuous, hierarchical organization of brain.MethodsWe aimed to investigate whether reproducible disruptions of functional gradient existed in SZ compared to normal controls (NC), and these abnormalities were associated with severity of clinical and cognitive symptoms in SZ. All analyses were implemented in two independent large-sample multi-site datasets (discovery dataset, 400 SZ and 336 NC; replication dataset, 279 SZ and 262 NC). First, functional gradient across cerebral cortex was calculated in each subject. Second, vertex-wise comparisons of cortical gradient between SZ and NC groups were performed to identify abnormalities in SZ. Meanwhile, reproducible and robustness analyses were implemented to validate these abnormalities. Finally, regression analyses were performed using generalized additive models to link these abnormalities to severity of clinical and cognitive symptoms in SZ.ResultsWe found an abnormal gradient map in SZ in the discovery dataset, which was reproducible in the replication dataset. The abnormal gradient pattern was also robust when performing methodological alternatives and control analyses. Further, these reproducible abnormalities can reliably predict symptoms of clinical and cognitive domains across the two independent datasets.ConclusionThese findings demonstrated that alterations in functional gradient can provide a reliable signature of SZ, characterizing the heterogenous symptoms of clinical or cognitive domains, and may be further investigated to understand the neurobiological mechanisms of these symptoms.Impact StatementIn our study, using functional gradient measure and statistical learning technology and two independent multi-site case-control resting-state fMRI datasets (discovery dataset: 736 subjects; replication dataset: 541 subjects), we comprehensively investigated functional hierarchical organization in the cerebral cortex of SZ and its association with interindividual severity of symptoms. We found reproducible and robust abnormalities of functional gradient existed in SZ, which provided a reliable signature to characterize negative and general psychopathology symptoms, as well as cognitive deficits. Our findings can provide new insights to understand the neurobiological mechanisms of clinical and cognitive symptoms in SZ.