Abstract
Major depressive disorder (MDD) is a clinically heterogeneous disorder. Its mechanism is still unknown. Although the altered intersubject variability in functional connectivity (IVFC) within gray-matter has been reported in MDD, the alterations to IVFC within white-matter (WM-IVFC) remain unknown. Based on the resting-state functional MRI data of discovery (145 MDD patients and 119 healthy controls [HCs]) and validation cohorts (54 MDD patients, and 78 HCs), we compared the WM-IVFC between the two groups. We further assessed the meta-analytic cognitive functions related to the alterations. The discriminant WM-IVFC values were used to classify MDD patients and predict clinical symptoms in patients. In combination with the Allen Human Brain Atlas, transcriptome-neuroimaging association analyses were further conducted to investigate gene expression profiles associated with WM-IVFC alterations in MDD, followed by a set of gene functional characteristic analyses. We found extensive WM-IVFC alterations in MDD compared to HCs, which were associated with multiple behavioral domains, including sensorimotor processes and higher-order functions. The discriminant WM-IVFC could not only effectively distinguish MDD patients from HCs with an area under curve ranging from 0.889 to 0.901 across three classifiers, but significantly predict depression severity (r = 0.575, p = 0.002) and suicide risk (r = 0.384, p = 0.040) in patients. Furthermore, the variability-related genes were enriched for synapse, neuronal system, and ion channel, and predominantly expressed in excitatory and inhibitory neurons. Our results obtained good reproducibility in the validation cohort. These findings revealed intersubject functional variability changes of brain WM in MDD and its linkage with gene expression profiles, providing potential implications for understanding the high clinical heterogeneity of MDD.