Elucidating individual aberrance is a critical first step toward precision medicine for heterogeneous disorders such as depression. The neuropathology of depression is related to abnormal inter-regional structural covariance indicating a brain maturational disruption. However, most studies focus on group-level structural covariance aberrance and ignore the interindividual heterogeneity. For that reason, we aimed to identify individualized structural covariance aberrance with the help of individualized differential structural covariance network (IDSCN) analysis.

T1-weighted anatomical images of 195 first-episode untreated patients with depression and matched healthy controls (n = 78) were acquired. We obtained IDSCN for each patient and identified subtypes of depression based on shared differential edges.

As a result, patients with depression demonstrated tremendous heterogeneity in the distribution of differential structural covariance edges. Despite this heterogeneity, altered edges within subcortical-cerebellum network were often shared by most of the patients. Two robust neuroanatomical subtypes were identified. Specifically, patients in subtype 1 often shared decreased motor network-related edges. Patients in subtype 2 often shared decreased subcortical-cerebellum network-related edges. Functional annotation further revealed that differential edges in subtype 2 were mainly implicated in reward/motivation-related functional terms.

In conclusion, we investigated individualized differential structural covariance and identified that decreased edges within subcortical-cerebellum network are often shared by patients with depression. The identified two subtypes provide new insights into taxonomy and facilitate potential clues to precision diagnosis and treatment of depression.

As a neuroprogressive illness, depression is accompanied by brain structural abnormality that extends to many brain regions. However, the progressive structural alteration pattern remains unknown.

To elaborate the progressive structural alteration of depression according to illness duration, we recruited 195 never-treated first-episode patients with depression and 130 healthy controls (HCs) undergoing T1-weighted MRI scans. Voxel-based morphometry method was adopted to measure gray matter volume (GMV) for each participant. Patients were first divided into three stages according to the length of illness duration, then we explored stage-specific GMV alterations and the causal effect relationship between them using causal structural covariance network (CaSCN) analysis.

Overall, patients with depression presented stage-specific GMV alterations compared with HCs. Regions including the hippocampus, the thalamus and the ventral medial prefrontal cortex (vmPFC) presented GMV alteration at onset of illness. Then as the illness advanced, others regions began to present GMV alterations. These results suggested that GMV alteration originated from the hippocampus, the thalamus and vmPFC then expanded to other brain regions. The results of CaSCN analysis revealed that the hippocampus and the vmPFC corporately exerted causal effect on regions such as nucleus accumbens, the precuneus and the cerebellum. In addition, GMV alteration in the hippocampus was also potentially causally related to that in the dorsolateral frontal gyrus.

Consistent with the neuroprogressive hypothesis, our results reveal progressive morphological alteration originating from the vmPFC and the hippocampus and further elucidate possible details about disease progression of depression.