Mood disorders are characterized by great heterogeneity in clinical manifestation. Uncovering such heterogeneity using neuroimaging-based individual biomarkers, clinical behaviors, and genetic risks, might contribute to elucidating the etiology of these diseases and support precision medicine.

We recruited 174 drug-naïve and drug-free patients with major depressive disorder and bipolar disorder, as well as 404 healthy controls. T1 MRI imaging data, clinical symptoms, and neurocognitive assessments, and genetics were obtained and analyzed. We applied regional gray matter volumes (GMV) and quantile normative modeling to create maturation curves, and then calculated individual deviations to identify subtypes within the patients using hierarchical clustering. We compared the between-subtype differences in GMV deviations, clinical behaviors, cell-specific transcriptomic associations, and polygenic risk scores. We also validated the GMV deviations based subtyping analysis in a replication cohort.

Two subtypes emerged: subtype 1, characterized by increased GMV deviations in the frontal cortex, cognitive impairment, a higher genetic risk for Alzheimer's disease, and transcriptionally associated with Alzheimer's disease pathways, oligodendrocytes, and endothelial cells; and subtype 2, displaying globally decreased GMV deviations, more severe depressive symptoms, increased genetic vulnerability to major depressive disorder and transcriptionally related to microglia and inhibitory neurons. The distinct patterns of GMV deviations in the frontal, cingulate, and primary motor cortices between subtypes were shown to be replicable.

Our current results provide vital links between MRI-derived phenotypes, spatial transcriptome, genetic vulnerability, and clinical manifestation, and uncover the heterogeneity of mood disorders in biological and behavioral terms.