Major depressive disorder recurrence and medication status shape brain network topology

Introduction: Major depressive disorder (MDD) is a highly prevalent and disabling psychiatric disorder. Human neuroimaging studies increasingly frame its neurobiological substrate in terms of alterations of large-scale brain network organization. Resting-state fMRI findings broadly align with this view, yet remaining highly heterogeneous, reflecting both clinical and analytic variability. Predominant approaches largely characterize functional organization in terms of pairwise relationships between regions, which may limit the ability to capture more global features of brain organization associated with depression and its course. Here we use novel methods from Topological Data Analysis to provide new perspectives on these brain-behavior associations Objectives: The main objective of this study was to determine whether whole-brain topological descriptors of functional connectivity capture alterations associated with MDD. In particular, we examined whether these features vary along clinically relevant dimensions of heterogeneity, focusing on illness course (single-episode vs. recurrent MDD) and current antidepressant medication status. To this end, we analyzed the area under the curve (AUC) of the zeroth and first Betti numbers (B and B), which index global network integration and higher-order cycle structure across scales, respectively. Methods: Neuroimaging (resting state fMRI) and phenotypic data were obtained from 1,490 participants (776 individuals with MDD and 714 never-depressed controls) included in the REST-meta-MDD project of the DIRECT consortium. Functional connectivity matrices were computed using Pearson correlations between regions defined by the Power-264 atlas and harmonized across sites using CovBat. Topological Data Analysis was applied to characterize whole-brain network organization across connectivity thresholds using B and B curves and their areas under the curve (AUC). B associated cycle counts and their coarse anatomical configurations were further examined. Group differences and effects of illness course and medication status were assessed using multiple linear and generalized linear regression models. Results: Compared with controls, individuals with MDD showed significantly higher B AUC, indicating altered global integration profiles across connectivity thresholds. B AUC did not differ by illness course, illness duration, or symptom severity. In contrast, B AUC showed a graded increase across illness courses, with the highest values observed in recurrent depression. This effect was driven by an increased number of one-dimensional cycles rather than greater cycle persistence. Medication status further modulated B alterations, with elevated values observed in unmedicated single-episode and medicated recurrent MDD. Anatomical decomposition revealed that higher-order alterations were primarily driven by inter-network configurations spanning multiple large-scale functional systems. Discussion:These findings indicate that MDD is associated with altered multiscale functional organization, combining less compact global integration with an increased prevalence of cross-network hole-defining cycles, with these features varying systematically with illness course and current medication, suggesting that topology-informed measures capture clinically relevant variability in brain organization in depression.