A structural equation modeling approach using behavioral and neuroimaging markers in major depressive disorder

Major depressive disorder (MDD) has consistently proven to be a multifactorial and highly comorbid disease. Despite recent depression-related research demonstrating causalities between MDD-related factors and a small number of variables, including brain structural changes, a high-statistical power analysis of the various factors is yet to be conducted. We retrospectively analyzed data from 155 participants (84 healthy controls and 71 patients with MDD). We used magnetic resonance imaging and diffusion tensor imaging data, scales assessing childhood trauma, depression severity, cognitive dysfunction, impulsivity, and suicidal ideation. To simultaneously evaluate the causalities between multivariable, we implemented two types of MDD-specified structural equation models (SEM), the behavioral and neurobehavioral models. Behavioral SEM showed significant results in the MDD group: Comparative Fit Index [CFI] = 1.000, Root Mean Square Error of Approximation [RMSEA]) = 0.000), with a strong correlation in the scales for childhood trauma, depression severity, suicidal ideation, impulsivity, and cognitive dysfunction. Based on behavioral SEM, we established neurobehavioral models showing the best-fit in MDD, especially including the right cingulate cortex, central to the posterior corpus callosum, right putamen, pallidum, whole brainstem, and ventral diencephalon, including the thalamus (CFI >0.96, RMSEA <0.05). Our MDD-specific model revealed that the limbic-associated regions are strongly connected with childhood trauma rather than depression severity, and that they independently affect suicidal ideation and cognitive dysfunction. Furthermore, cognitive dysfunction could affect impulsivity.