Convergent molecular, cellular, and neural signatures of major depressive disorder

Major depressive disorder emerges from the complex interactions of biological systems that span across genes and molecules through cells, circuits, networks, and behavior. Establishing how neurobiological processes coalesce to contribute to the onset and maintenance of depression requires a multi-scale approach, encompassing measures of brain structure and function as well as genetic and cell-specific genomic data. Here, we examined anatomical (cortical thickness) and functional (functional variability, global brain connectivity) correlates of depression and negative affect across three population-imaging datasets: UK Biobank, Genome Superstruct Project, and ENIGMA (combined N[≥]23,723). Integrative analyses incorporated measures of cortical gene expression, post-mortem patient transcriptional data, depression GWAS, and single-cell transcription. Neuroimaging correlates of depression and negative affect were consistent across the three independent datasets. Linking ex-vivo gene downregulation with in-vivo neuroimaging, we found that genomic correlates of depression-linked neuroimaging phenotypes tracked gene downregulation in post-mortem cortical tissue samples of patients with depression. Integrated analysis of single-cell and Allen Human Brain Atlas expression data implicated somatostatin interneurons and astrocytes as consistent cell associates of depression, through both in-vivo imaging and ex-vivo cortical gene dysregulation. Providing converging evidence for these observations, GWAS derived polygenic risk for depression was enriched for genes expressed in interneurons, but not glia. Underscoring the translational potential of multi-scale approaches, the genomic correlates of depression-linked brain function and structure were enriched for known and novel disorder relevant molecular pathways. These findings bridge across levels to connect specific genes, cell classes, and biological pathways to in-vivo imaging correlates of depression. Key FindingsO_LIMajor depressive disorder and negative affect are associated with replicable profiles of cortical anatomy and function across independent population-level neuroimaging datasets (combined N[≥]23,723). C_LIO_LISomatostatin interneurons are consistent spatial transcriptional associates of in-vivo depression-linked imaging phenotypes. C_LIO_LIIntegrative single-cell gene expression analysis associate somatostatin interneurons and astrocytes with both in-vivo depression-linked imaging and ex-vivo gene downregulation in independent MDD cortical tissue samples. C_LIO_LITranscriptional correlates of in-vivo depression imaging phenotypes selectively capture gene downregulation in post-mortem tissue samples from patients with depression, but not other psychiatric disorders. C_LIO_LIIndicating that some cell classes are preferentially sensitive to inherited disease liability, genome-wide risk for depression is enriched among interneurons, but not glia. C_LIO_LIGene associates of depression-linked anatomy and function identify specific neurotransmitter systems, molecular signaling pathways, and receptors, suggesting possible targets for pharmaceutical intervention. C_LI