
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39223129
52121
10.1038/s41467-024-52121-y
Article
Investigating the shared genetic architecture between depression and subcortical volumes
Liu Mengge 1
Wang Lu 2
Zhang Yujie 1
Dong Haoyang 1
Wang Caihong 13
Chen Yayuan 1
Qian Qian 1
Zhang Nannan 1
Wang Shaoying 1
Zhao Guoshu 1
Zhang Zhihui 1
Lei Minghuan 1
http://orcid.org/0000-0001-8949-303X
Wang Sijia wangsijia@tmu.edu.cn

1
http://orcid.org/0009-0003-2960-3024
Zhao Qiyu zhaoqy1216@163.com

1
http://orcid.org/0000-0002-3570-4222
Liu Feng fengliu@tmu.edu.cn

1
1 https://ror.org/003sav965 grid.412645.0 0000 0004 1757 9434 Department of Radiology and Tianjin Key Laboratory of Functional Imaging & Tianjin Institute of Radiology, Tianjin Medical University General Hospital, Tianjin, China
2 https://ror.org/003sav965 grid.412645.0 0000 0004 1757 9434 Department of Geriatrics and Tianjin Geriatrics Institute, Tianjin Medical University General Hospital, Tianjin, China
3 https://ror.org/056swr059 grid.412633.1 Department of Magnetic Resonance Imaging, the First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
2 9 2024
2 9 2024
2024
15 764722 1 2024
27 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Depression, a widespread and highly heritable mental health condition, profoundly affects millions of individuals worldwide. Neuroimaging studies have consistently revealed volumetric abnormalities in subcortical structures associated with depression. However, the genetic underpinnings shared between depression and subcortical volumes remain inadequately understood. Here, we investigate the extent of polygenic overlap using the bivariate causal mixture model (MiXeR), leveraging summary statistics from the largest genome-wide association studies for depression (N = 674,452) and 14 subcortical volumetric phenotypes (N = 33,224). Additionally, we identify shared genomic loci through conditional/conjunctional FDR analyses. MiXeR shows that subcortical volumetric traits share a substantial proportion of genetic variants with depression, with 44 distinct shared loci identified by subsequent conjunctional FDR analysis. These shared loci are predominantly located in intronic regions (58.7%) and non-coding RNA intronic regions (25.4%). The 269 protein-coding genes mapped by these shared loci exhibit specific developmental trajectories, with the expression level of 55 genes linked to both depression and subcortical volumes, and 30 genes linked to cognitive abilities and behavioral symptoms. These findings highlight a shared genetic architecture between depression and subcortical volumetric phenotypes, enriching our understanding of the neurobiological underpinnings of depression.

Depression affects millions of people worldwide. Here, the authors show a substantial polygenic overlap between depression and brain subcortical volumes, identifying 44 shared loci.

Subject terms

Disease genetics
Genome-wide association studies
Molecular neuroscience
Depression
Natural Science Foundation of China (82072001) Beijing-Tianjin-Hebei Basic Research Collaboration Project (J230040) Tianjin Natural Science Foundation (23JCZXJC00120) Tianjin Key Medical Discipline (Specialty) Construction Project (TJYXZDXK-001A)issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Depression is a complex and pervasive mental health disorder that affects millions of people worldwide1. It is characterized by a range of behavioral symptoms, including persistent feelings of sadness, hopelessness, worthlessness, and a lack of interest or pleasure in daily activities. These symptoms are often accompanied by concentration difficulties, changes in appetite, sleep disturbances—either insomnia or hypersomnia—and suicide attempts2–4. Additionally, cognitive impairments, which occur in 85–94% of clinical cases5, affect various domains such as attention, executive functions, memory, and processing speed6,7. The diverse nature of depression extends beyond these aspects, as researchers have explored its connection with neurobiology. Subcortical brain structures, such as the amygdala and hippocampus, play a crucial role in regulating emotions, memory, and cognitive functions, establishing their significance in the neurobiology of depression.

Neuroimaging studies have progressively revealed volumetric changes within subcortical structures in depression. For instance, the amygdala, pivotal in emotional processing, often exhibits alterations in volume under depressive conditions8–10. These changes may contribute to increased emotional reactivity and the heightened emotional experiences commonly observed in individuals with depression11. Similarly, the hippocampus, crucial for memory formation and stress regulation, frequently shows a reduction in volume among those with depression12–14. This hippocampal atrophy is associated with cognitive deficits and difficulties in stress response regulation in individuals experiencing depression15. The observed alterations in the nucleus accumbens, a key region involved in reward processing, suggest a broader impact of depression on neural circuits related to emotional and reward-related processes16. Additionally, other subcortical regions, such as the thalamus, caudate, and putamen, also exhibit structural changes in depression17–19. These structural alterations highlight the influence of depression on the physical dimensions of critical brain regions, providing insight into the neurobiological underpinnings of the disorder.

Recent advances in genetic research, particularly through genome-wide association study (GWAS), have provided invaluable insights into the genetic factors underlying depression and subcortical volumetric traits. In a meta-GWAS involving over 1.3 million individuals, 243 independent risk loci for depression were identified20. Concurrently, ongoing investigations into the genetic architecture of subcortical volumes have contributed significantly to our understanding of the complex interactions between genetic factors and these traits21–23. Given the high heritability and polygenic nature observed in both depression24 and subcortical volumetric traits21, it is likely that the neuroimaging findings could be attributed to a shared genetic basis influencing these conditions. Several studies exploring this genetic overlap, using approaches such as linkage disequilibrium score regression (LDSC) and polygenic risk score (PRS), have uncovered weak or absent genome-wide genetic correlations25–27. However, these approaches might underestimate genetic overlap, particularly in scenarios involving a balanced mixture of genetic variants with both concordant and discordant effects on the phenotypes. Recently, Frei et al. introduced the bivariate causal mixture model (MiXeR)28, a method that evaluates genetic overlap by considering all variants, regardless of the direction of effect. Furthermore, conditional and conjunctional false discovery rate (cond/conjFDR) analyses, which utilize an empirical Bayesian statistical framework, enhance the power of both GWASs in discovering novel risk loci and identifying shared genetic loci for the two traits29–31. These innovative approaches have been successfully applied in various studies, illuminating the genetic relationships between complex traits32,33.

In this work, we aim to identify the shared genetic basis between depression and subcortical volumes. Specifically, the MiXeR method is employed to investigate the shared polygenic architecture. Following this, the cond/conjFDR approaches are employed to identify both novel loci for each trait and the genetic loci shared between these traits. Furthermore, functional annotation and gene mapping of the shared loci are conducted. The expression patterns of these shared genes are then examined across various developmental stages, and their associations with depression and subcortical volumes, as well as their potential impacts on cognitive abilities and behavioral symptoms, are explored. By integrating advanced statistical models with recent large-cohort GWASs, the power to discover shared genomic loci in depression and subcortical brain structures is increased, thereby better revealing the shared genetic basis between depression and subcortical volumes. A systematic flow of the study design is shown in Supplementary Fig. 1.

Results

Quantification of genetic overlap with MiXeR

Using univariate MiXeR34, we estimated the heritability, polygenicity, and discoverability for depression and each subcortical volumetric phenotype. All 14 subcortical phenotypes showed a good model fit, as indicated by positive Akaike information criterion (AIC) values. As shown in Fig. 1a and Supplementary Data 1, univariate MiXeR analysis revealed that depression exhibited lower single nucleotide polymorphism (SNP) heritability (0.051), higher polygenicity (9400 trait-influencing variants), and lower discoverability (8.3E−6) compared with all subcortical volumetric phenotypes. In contrast, these volumetric phenotypes demonstrated SNP heritability ranging from 0.181 to 0.295, polygenicity ranging from 939 to 3207 trait-influencing variants, and discoverability ranging from 1.2E−4 to 3.5E−4. Additionally, due to the lower discoverability, the estimated sample size required to reach 90% heritability was more than ten times larger for depression than for subcortical volumetric phenotypes (65.8 million vs. 2.5–6.6 million, Supplementary Fig. 2, Supplementary Data 1).Fig. 1 Genetic overlap between depression and subcortical volumes using MiXeR.

a Heritability, polygenicity, and discoverability for each phenotype estimated by univariate MiXeR. b Venn diagrams showing MiXeR results of the estimated number of variants shared between depression and subcortical volumes. The numbers of shared and trait-specific variants are shown in thousands with standard errors. The size of the circles reflects the degree of polygenicity. c The points indicate the proportions of shared variants with concordant effect directions. Source data are provided as a Source Data file. L left, R right.

Applying bivariate MiXeR28, we found considerable quantities of shared variants in each pair of depression-volumetric phenotypes. Of the 14 subcortical structures, eight demonstrated a good model fit for MiXeR analyses with depression (both AICbest vs. min and AICbest vs. max were positive; Supplementary Data 2). As shown in Fig. 1b and Supplementary Data 2, subcortical volumetric traits shared 23.4% to 83.3% of their trait-influencing variants with depression. Conversely, depression shared 4.1% to 17.9% of its trait-influencing variants with subcortical structures. The Dice coefficient, which is the ratio of shared variants to the total number of variants, ranged from 7.4% to 28.5%, indicating varying degrees of genetic overlap between depression and the various subcortical traits. Additionally, the allelic effects of the variants shared between depression and subcortical volumes exhibited divergent directions, with the proportion of shared variants with concordant effects ranging from 19.7% to 73.6% (Fig. 1c, Supplementary Data 2). Within the shared component, the bilateral caudate, bilateral putamen, and right pallidum demonstrated positive correlations with depression, with the largest coefficient observed in the left caudate (rho = 0.634). In contrast, other subcortical structures displayed negative correlations with depression, with the largest coefficient observed in the left hippocampus (rho = -0.784), as shown in Supplementary Data 2. Details of the MiXeR results are provided in Supplementary Figs. 3–16.

Identification of novel and shared genetic loci (cond/conjFDR)

To improve the discovery of genetic variants associated with depression and subcortical volumes, we applied condFDR analysis29–31 to examine the associations of depression with subcortical volumes, and vice versa. By conditioning on subcortical volumetric phenotypes, we identified 94 distinct loci associated with depression at condFDR <0.01, including 7 novel loci (Supplementary Data 3, Supplementary Fig. 17). The candidate SNPs in these novel loci were mapped to 14 genes (Supplementary Data 4). Enrichment analysis was conducted on these genes to explore potential biological pathways, however, due to the small number of genes, no significant results were obtained. Additionally, when conditioning on depression, a total of 410 loci associated with the volume of the thalamus (left N = 22, right N = 28), caudate (left N = 43, right N = 41), putamen (left N = 44, right N = 43), pallidum (left N = 20, right N = 35), hippocampus (left N = 25, right N = 33), amygdala (left N = 19, right N = 12), and accumbens (left N = 25, right N = 20) were identified (Supplementary Data 5, Supplementary Fig. 18). Combining loci across all subcortical volumetric traits resulted in a total of 233 distinct loci, among which 13 were newly identified (Supplementary Data 6). The candidate SNPs in these novel loci were mapped to 57 genes (Supplementary Data 7), and subsequent enrichment analysis revealed associations with several gene ontology (GO) biological process terms, including cerebral cortex GABAergic interneuron development (FDR q = 4.0E−2) and pyrimidine nucleoside biosynthetic process (FDR q = 4.0E−2).

The conditional quantile‒quantile (Q‒Q) plots demonstrated substantial genetic overlap for each depression-volume pair (Supplementary Figs. 19–32). Thus, conjFDR analysis was subsequently performed to identify genomic loci jointly associated with both traits. As shown in Fig. 2a, b and Supplementary Data 8, conjFDR unveiled 77 loci shared between depression and subcortical volumetric phenotypes: the thalamus (left N = 9, right N = 13), caudate (left N = 8, right N = 4), putamen (left N = 3, right N = 6), pallidum (left N = 2, right N = 6), hippocampus (left N = 3, right N = 9), amygdala (left N = 4, right N = 5), and accumbens (left N = 2, right N = 3). Among them, the most significant locus was chr18: 50358109-51055069 for depression and bilateral putamen volume (lead SNP: rs4632195, left: conjFDR = 1.0E−8; right: conjFDR = 4.3E−9; Fig. 2c). Of the 77 shared loci, a total of 51 genomic loci were jointly associated with depression and at least two subcortical volumetric phenotypes, resulting in 44 distinct shared loci across all depression-volume pairs. Within these shared loci, 7 were novel for depression, 17 were novel for subcortical volumetric phenotypes, and 2 were novel for both (Supplementary Data 9).Fig. 2 Results of loci shared between depression and subcortical volumes.

a Manhattan plot showing the associations between depression and subcortical volumes identified by conjFDR. Each point represents a SNP and the point with a larger size and black border indicates a lead SNP. The y-axis represents the -log10-transformed conjFDR values for each SNP, and the x-axis reflects the chromosomal position. The dashed line represents the threshold of conjFDR = 0.05, and the color of the dots represents different phenotypes. b Brain map showing the number of loci of each subcortical structure shared with depression. A color transition from white to red indicates an increasing number of shared loci. c LocusZoom plots of a genomic locus (chr18: 50358109-51055069, lead SNP: rs4632195) showing the regional SNP associations (two-sided -log10(p)) for depression and bilateral putamen volumes from the original GWAS summary data d. MRI slice plot of bilateral putamen volumes. Source data are provided as a Source Data file. Chr chromosome, L left, R right.

The primary locations of the candidate SNPs in the distinct shared loci were within intronic regions (58.7%) and non-coding RNA intronic regions (25.4%), as illustrated in Fig. 3a and detailed in Supplementary Data 10. Additionally, 18 SNPs with non-synonymous mutations in exons were identified. Approximately 3.7% of the candidate SNPs were indicated as pathogenic, with combined annotation-dependent depletion (CADD) scores exceeding 12.37. About 10.3% of the candidate SNPs exhibited a higher likelihood of regulatory functionality, as indicated by a RegulomeDB level of less than 3. Moreover, 90.1% of the candidate SNPs were located in regions with a predominantly open chromatin configuration, characterized by a chromatin state score of less than 8. Among the 44 lead SNPs, five (rs11047209, rs1976052, rs2246556, rs4593875, rs62618693) had CADD larger than 12.37, with rs62618693, a non-synonymous mutation in QSER1, having the highest CADD value of 31 (Supplementary Data 9). Furthermore, the candidate SNPs in 44 distinct shared loci were mapped to 269 protein-coding genes (Supplementary Data 11), and enrichment analysis revealed 42 GO biological process terms (FDR q < 0.05), mainly including the regulation of neuron projection development (FDR q = 3.4E−3), regulation of axonogenesis (FDR q = 1.1E−2), and central nervous system development (FDR q = 2.8E−2); see Supplementary Data 12 and Fig. 3b for details.Fig. 3 Results of functional annotation and gene expression analyses for the shared loci.

a Distribution of functional consequences of candidate SNPs in the genomic loci shared between depression and 14 subcortical volumes. b The bubble plot shows the enriched GO terms of biological processes for protein-coding genes. The x-axis indicates -log10-transformed values of FDR-adjusted two-sided p-values, and the y-axis shows each GO term. The color of the dots indicates the magnitude of the -log10(q) values, and the size of the dots represents the number of genes enriched in a given term. c Developmental trajectories of average gene expression in whole brain postmortem-bulk RNA data (PsychEncode) for all genes mapped by shared candidate SNPs (conjFDR < 0.05) between depression and seven subcortical volumetric phenotypes (left and right subcortical phenotypes are combined through averaging). The shaded region indicates 95% confidence intervals, and the line indicates the usual timing of birth. The x-axis represents the nine developmental windows, and the y-axis represents the log2-transformed mean expression value. d The plot shows all significant expression-trait associations (two-sided q < 0.05, FDR corrected) identified by S-MultiXcan for both depression and at least one subcortical volume. The x-axis provides the names of 55 genes and the y-axis lists the names of depression and 14 subcortical structures. The magnitudes of the -log10-transformed q values are coded as different colors from white to red. Source data are provided as a Source Data file. L left, ncRNA non-coding RNA, PCW postconceptional week, PY postnatal year, R right.

To quantify the differences and similarities in how shared lead SNPs identified by conjFDR influence different subcortical volumetric phenotypes, we utilized Cochran’s Q test to detect heterogeneity in 44 shared lead SNPs across 14 subcortical volumetric phenotypes, and for each volumetric phenotype pair. We found that 30 shared lead SNPs showed heterogeneity across all subcortical phenotypes (p < 0.05), whereas 14 lead SNPs showed no evidence of effect heterogeneity (Supplementary Data 13). Additionally, the heterogeneity for each pair of subcortical phenotypes is shown in Supplementary Data 14 and Supplementary Fig. 33. Specifically, we observed no heterogeneity in shared lead SNPs between bilateral subcortical structures, except for the bilateral accumbens, which exhibited heterogeneity in one locus. Among the other subcortical phenotype pairs, the pair with the highest number of heterogeneous SNPs was the left caudate and the right hippocampus, with 20 out of 44 SNPs showing heterogeneity. In contrast, the pair with the lowest number of heterogeneous SNPs was the left hippocampus and the left amygdala, with only 2 out of 44 SNPs showing heterogeneity.

Developmental trajectories of shared genes

Using spatiotemporal brain gene expression data from PsychENCODE35, we investigated the developmental expression trajectories of genes shared between depression and each subcortical structure. As shown in Fig. 3c, genes shared between depression and the accumbens, pallidum and thalamus showed similar trajectories with peak expression during the mid-gestation period (windows 3-4, post-conceptional weeks [PCW] 16-22), followed by a slight decline in expression, and then entering a plateau phase after window 6 (postnatal years [PY] 0.5–2.6). The expression of genes in the hippocampus, amygdala, and caudate displayed a peak expression level during the early prenatal period (window 1, PCW 8-9), followed by a continuous decline until reaching the lowest point in PY 0.5–2.6. Regarding the putamen, its expression peaked at PCW 8-9 and then decreased steadily.

Expression-trait associations

For each depression-volume pair, we investigated whether the expression of shared genes in the brain was associated with both depression and the corresponding subcortical volumes using the S-MultiXcan tool36. Out of these shared genes, we found 137 significant associations between gene expression (involving 71 genes) and volumetric phenotypes. Among these 71 genes, 55 genes were also found to be associated with depression (Supplementary Data 15, 16, Fig. 3d). For example, the expression levels of LRRC37A2, SPPL2C, ARL17A, FMNL1, PLEKHM1, ARHGAP27, and CRHR1 in the brain were associated with the volumes of the bilateral thalamus, bilateral pallidum, left accumbens, and right hippocampus, as well as with depression. Moreover, the expression of NPEPPS, GINM1, PRMT6, EFCAB13, NUP43, LATS1, ULBP2, ACYP1, and AREL1 showed correlations with the bilateral thalamic volume and depression. The MSH6 was correlated with the bilateral hippocampal volume and depression. The most significant expression-depression association was found for CELSR3, which was also linked to the volume of the right pallidum.

Cognitive and behavioral associations

To investigate whether the genes shared between depression and subcortical volumetric phenotypes exerted further influence on cognitive abilities or behavioral symptoms associated with depression, GWAS summary statistics for 8 cognitive abilities and 15 behavioral symptoms were acquired (Supplementary Data 17, 18). The associations between the expression of the 55 shared genes (i.e., those associated with depression and at least one subcortical volume) and cognitive abilities, as well as behavioral symptoms, were also analyzed using S-MultiXcan. Out of these 55 shared genes, we found 30 significant associations between gene expression and behavioral symptoms, as well as 59 significant associations with cognition, involving a total of 30 genes (Fig. 4a, Supplementary Data 19). Notably, the expression of five genes in the brain, including ARL17A, CRHR1, GPX1, SPPL2C, and LRRC37A2, had significant associations with more than five cognitive or behavioral traits, indicating their widespread influence on cognition and behavior. The most significant expression-cognition association was found for LRRC37A2 with reaction time (RT), which was also correlated with depression and six subcortical structures.Fig. 4 Associations of the expression of shared genes in the brain with cognitive abilities and behavioral symptoms.

a Bubble plot showing significant expression-trait associations (two-sided q < 0.05, FDR corrected). The x-axis shows the shared genes, and the y-axis shows cognitive abilities and behavioral symptoms. The size and color of the bubbles demonstrate the significance level. The gray bubbles represent non-significant associations. b The bar plots show the significance level (uncorrected two-sided p-values) of the expression of LRRC37A2 in five subcortical structures with the VNR and RT. c The bar plots show the significance level (uncorrected two-sided p-values) of PLEKHM1 expression in five subcortical structures associated with feelings of worthlessness during the worst period of depression and VNR. Source data are provided as a Source Data file. ConcWorst difficulty concentrating during worst depression, DSST digit symbol substitution task, Inadequacy recent feelings of inadequacy, Interest recent lack of interest or pleasure in doing things, Movement recent changes in speed/amount of moving or speaking, PM pairs matching, RT reaction time, SleepIssues trouble falling or staying asleep, or sleeping too much, SleepMuch sleeping too much, SleepTrouble trouble falling asleep, TiredRecent recent feelings of tiredness or low energy, TiredWorst feelings of tiredness during the worst episode of depression, TMT-B trail making test B, VNR verbal–numerical reasoning, WorthlessWorst feelings of worthlessness during worst period of depression.

Additionally, we observed that the gene expression levels of the same gene across different subcortical regions were associated with distinct cognitive abilities or behavioral symptoms (Fig. 4b, c). For instance, the expression of LRRC37A2 in the amygdala exhibited the strongest association with verbal–numerical reasoning (VNR), whereas its expression in the accumbens did not correlate with VNR. Conversely, the expression of LRRC37A2 in the accumbens was predominantly associated with RT (Fig. 4b). Furthermore, the expression of PLEKHM1 in the hippocampus was most significantly associated with feelings of worthlessness during the worst period of depression, while its expression in the caudate was primarily associated with VNR (Fig. 4c).

Control analyses

Despite previous studies finding weak or no genome-wide genetic correlation between depression and subcortical volumes25,26, we employed LDSC analyses37,38 to investigate the correlation using updated GWAS summary statistics, in order to confirm or refute the relationship. We found that the absolute genetic correlations ranged from 0.002 for left pallidum volume to 0.093 for right caudate volume, with nominal statistical significance only observed for bilateral caudate volume (left p = 0.015, right p = 0.004; see details in Supplementary Data 20).

To validate the shared loci identified by conjFDR analysis, further analysis was undertaken using local analysis of [co]variant association (LAVA)39. In general, LAVA analyses supported the findings of the conjFDR analysis. Among the 77 shared loci, LAVA identified 56 genetic loci showing nominally significant local heritability in both depression and volumetric phenotypes. Within these 56 loci, 33 exhibited a significant local genetic correlation (FDR q < 0.05), with 20 displaying a positive correlation (ranging from 0.63 to 1.00) and 13 showing a negative correlation (ranging from -0.54 to -1.00), as shown in Supplementary Data 8. The most significant positive local genetic correlation was observed at the same locus (chr18: 50358109-51055069, lead SNP: rs4632195) between depression and the left putamen (r = 0.64, p = 6.0E−4, conjFDR = 1.0E−8), as well as the right putamen (r = 0.63, p = 4.5E−4, conjFDR = 4.3E−9). The most significant negative local genetic correlation was observed between depression and right thalamic volume at chr4: 3247007-3285389 (lead SNP: rs3135170, r = −1.00, p = 1.8E-3, conjFDR = 0.044).

Discussion

Here, we explored the shared genetic architecture between depression and subcortical volumetric phenotypes. Bivariate MiXeR analysis revealed that subcortical volumetric traits shared a substantial proportion of trait-influencing variants with depression, with the left amygdala showing the highest Dice coefficient of 28.5%. The identification of shared loci, encompassing 44 distinct loci associated with both depression (including 7 novel loci) and subcortical volumes (including 17 novel loci), through conjFDR analysis enriches our understanding of their intricate genetic connections. Functional annotation, gene mapping, and gene expression analyses provide vital information about the molecular mechanisms behind the identified shared loci. Unraveling these mechanisms will deepen our comprehension of how genetic factors contribute to the co-occurrence of depression and alterations in subcortical brain volumes.

In previous studies utilizing LDSC25,26, it was demonstrated that there are weak or absent genome-wide genetic correlations between depression and subcortical brain volumes, even when applying a lenient statistical threshold (i.e., uncorrected p < 0.05). For example, one study found no significant genetic correlation between depression and subcortical volumes25, and another study found that depression was weakly positively correlated with hippocampal volume26. In the present study, we validated these findings using updated GWAS summary statistics, revealing nominal significance exclusively in the genetic correlation between depression and bilateral caudate volume (Supplementary Data 20). Notably, genome-wide genetic correlations can only capture consistent allelic effect directions among shared variants and do not consider mixed patterns of directional effects. In contrast, MiXeR analysis indicated that all subcortical volumes shared a significant proportion of variants (ranging from 23.4% to 83.3%) with depression. These shared variants exhibit divergent directions, as evidenced by the proportion of shared variants with concordant effects ranging from 19.7% to 73.6%. This balanced mixture of directional effects aligns with the nonsignificant genome-wide genetic correlation between depression and subcortical volumetric phenotypes, suggesting the potential involvement of both agonistic and antagonistic overlapping molecular genetic mechanisms in its pathophysiology.

The condFDR analysis, focused on uncovering novel genetic variants associated with depression and subcortical volumetric phenotypes, yielded compelling results. Under the stringent condFDR threshold (< 0.01), conditioning on subcortical volumetric phenotypes revealed 94 distinct loci linked to depression, including 7 novel loci. Conversely, conditioning on depression identified 233 distinct loci associated with subcortical volumetric phenotypes, including 13 novel loci. Among the lead SNPs of these 13 novel loci, rs636926 is a nonsynonymous exonic variant within TMEM171 (Supplementary Data 6). This variant results in a change in the corresponding amino acid, potentially directly impacting the structure or function of the protein. Subsequent enrichment analysis linked 57 mapped genes from these 13 novel loci predominantly to the development of GABAergic interneurons in the cerebral cortex and pyrimidine nucleoside biosynthetic process. GABAergic interneurons are inhibitory neurons in the nervous system that play a vital role in the functional maturation of the cortex and hippocampal neural circuitry and activity40,41. Several studies have indicated that alterations in GABAergic interneurons contribute to clinical features in neuropsychiatric disorders such as autism spectrum disorder42 and schizophrenia43. Dysregulation of pyrimidine nucleoside biosynthetic process may influence neurogenesis, synaptic plasticity, and neurotransmitter signaling within subcortical structures, potentially leading to changes in the brain’s structural architecture. By exploring the biological functions of these genes, biological processes or pathways related to the functioning of subcortical brain regions can be revealed, providing valuable insights for further investigation.

The Q‒Q plots demonstrate a significant genetic overlap for each pair of depression and subcortical volume, as illustrated in Supplementary Figs. 19–32. Additionally, the conjFDR analysis identified 77 loci shared between depression and subcortical volumetric phenotypes. Within the shared loci, there are both concordant effects (42 out of 77 SNPs) and discordant effects (35 out of 77 SNPs) in the top lead SNPs (Supplementary Data 8), providing additional insights into the lack of a genome-wide genetic correlation revealed by the LDSC25. Notably, there were 44 distinct shared loci across all depression-volume pairs, indicating unique associations that contribute to the shared genetic architecture between depression and subcortical brain structures. The most significant shared locus, chr18: 50358109-51055069 (lead SNP: rs4632195), between depression and bilateral putamen volumes (left: conjFDR = 1.0E−8; right: conjFDR = 4.3E−9; Supplementary Data 8), showed positive local genetic correlations (left: rho = 0.644, right: rho = 0.627). This locus is located within the intronic region of the DCC gene, which encodes a netrin receptor involved in axon guidance and migration. Previous studies have linked SNPs within DCC to putamen volume44 and suggested correlations with depression45,46 and schizophrenia47. Moreover, we observed that rs62618693 (CADD = 31), a non-synonymous mutation within the exons of the QSER1, was jointly associated with bilateral amygdala volume and depression, and the locus (chr11: 32623621-32956492) with rs62618693 as the lead SNP was a novel locus for subcortical volume.

Among the 44 distinct shared loci, two are novel for both depression and subcortical volumes. One novel locus (chr2: 48024876-48024876, lead SNP: rs3136329, conjFDR = 3.5E−2, Supplementary Data 9) is located in FBXO11, which is ubiquitously expressed in the brain (reads per kilobase of transcript per million mapped reads [RPKM] 21.6) and has been found to be associated with schizophrenia48. Another novel locus (chr10: 98011878-98162702, lead SNP: rs11188695, conjFDR = 1.6E−2, Supplementary Data 9) can be mapped to TLL2 by expression quantitative trait loci (eQTL) mapping. TLL2 encodes an astacin-like zinc-dependent metalloprotease and is associated with stress sensitivity49. Additionally, the lead SNPs in other shared loci implicate genes like CDK14, NRD1, DYRK1A, KLF11, ARHGEF10, TWF1P1, and LRRC71 in the genetic relationship between depression and subcortical volumetric phenotypes (Supplementary Data 11). Furthermore, the candidate SNPs in the loci shared between depression and subcortical volumetric traits were mapped to 269 protein-coding genes. Enrichment analysis of these genes elucidated their involvement in essential biological processes, such as the regulation of neuron projection development, regulation of axonogenesis, and central nervous system development. These findings provide valuable insights into potential molecular mechanisms and pathways underlying the connection between depression and subcortical brain structures.

By employing S-MultiXcan to explore expression-trait associations, our study revealed 55 shared genes that not only exhibited correlations with volumetric phenotypes but also demonstrated associations with depression. For example, the expression of LRRC37A2, SPPL2C, ARL17A, FMNL1, PLEKHM1, ARHGAP27, and CRHR1 in the brain is linked to the bilateral thalamus, bilateral pallidum, left accumbens, right hippocampus and depression, suggesting shared genetic underpinnings between subcortical brain characteristics and depression. Moreover, we found 89 significant associations between the expression of shared genes and cognitive abilities as well as behavioral symptoms, involving a total of 30 genes. Among these, five genes—ARL17A, CRHR1, GPX1, SPPL2C, and LRRC37A2—had significant associations with more than five cognitive or behavioral traits. Our study also revealed diverse implications for these genes. For instance, PLEKHM1 variants are associated with neuroticism50 and global cortical volume51. LRRC37A2 and ARHGAP27 have been linked to Parkinson’s disease52,53. Additionally, CRHR1 encodes a G protein-coupled receptor implicated in the regulation of the hypothalamic-pituitary-adrenal axis, which is associated with stress-related psychopathology54–56. CRHR1 has been identified as a genetic risk factor for affective disorders and is linked to metabolic activity in the anterior hippocampus and amygdala57. ARL17A is involved in intracellular protein and vesicle-mediated transport and is a risk gene for progressive supranuclear palsy58. GPX1, a critical component of the antioxidant system involved in the major antioxidant molecule glutathione, may affect cognition and emotional behavior through mechanisms involving oxidative stress-induced neuronal damage, synaptic dysfunction, and neurotransmitter dysregulation, and is also related to the risk of depression59,60. These complex relationships provide a comprehensive understanding of the interactions among gene expression, subcortical volumetric phenotypes, depression, and cognitive abilities as well as behavioral symptoms, offering valuable perspectives into potential contributing mechanisms.

Our findings should be considered for the following limitations. First, the analysis was restricted to individuals of European ancestry due to the absence of adequately powered GWAS on other ancestries. Future research in large non-European populations is crucial to enhance the external validity and generalizability of our results. Second, our results are confined to common variants on autosomal chromosomes, while rare variants also play an important role in depression61,62. Therefore, future studies on rare variants are warranted. Third, the absence of subtype-specific data restricts our ability to thoroughly investigate the shared genetic mechanisms underlying various depression subtypes and subcortical brain structures, and our findings may not fully elucidate the associations between shared genes and the diverse clinical symptoms observed across different depression subtypes. This emphasizes the need for future research efforts to incorporate subtype-specific data to enhance our understanding of this condition. Fourth, both the MiXeR and condFDR/conjFDR approaches are limited to bivariate analyses, preventing the determination of whether the identified genetic overlap is influenced by additional factors. Finally, the observation of negative AIC values in the bivariate MiXeR analysis indicates that these results should be interpreted with caution. Future studies with larger sample sizes and more comprehensive data may be necessary to validate our findings.

In conclusion, our study employed innovative methods, including MiXeR and cond/conjFDR analyses, to explore the genetic overlap between depression and subcortical brain volumes. By unraveling the shared genetic architecture, we identified several loci associated with both depression and specific volumetric traits, and uncovered novel loci that provide insights into the molecular mechanisms underlying this complex relationship. The functional annotation highlighted pathways related to neuron projection development, regulation of axonogenesis, and central nervous system development. These findings enhance our understanding of the shared genetic basis of depression and subcortical brain structures, contributing valuable knowledge to the intersection of psychiatric genetics and neurobiology.

Methods

GWAS data

The GWAS summary statistics for depression were derived from a recent meta-analysis involving 371,184 individuals diagnosed with depression and 978,703 controls20. Participants from the UK Biobank were excluded to avoid sample overlap, and samples from 23andMe were omitted due to access restrictions, resulting in a GWAS subsample consisting of 166,773 cases and 507,679 controls. Additionally, summary statistics for 14 subcortical volumetric phenotypes, encompassing bilateral volumes of the thalamus, caudate, putamen, pallidum, hippocampus, amygdala and accumbens, were acquired from the UK Biobank with a sample size of 33,224 individuals23. Notably, both the depression and volumetric trait datasets consisted of individuals of European ancestry. For further details regarding the GWAS summary statistics, please refer to the Supplementary Methods.

Univariate and bivariate MiXeR analyses

The MiXeR (v1.3, https://github.com/precimed/mixer) approach was used to assess the polygenic overlap between depression and subcortical volumetric traits28,34. Univariate MiXeR analyses were first conducted to estimate heritability, polygenicity, and discoverability for each phenotype. Heritability was calculated as the cumulative effect of all trait-influencing variants, which are common variants exhibiting a non-zero additive effect on the trait. Polygenicity was defined as the number of trait-influencing variants accounting for 90% of the heritability for each phenotype. Traits with higher polygenicity have a larger number of genetic variants affecting them than those with lower polygenicity. Discoverability, measured as the variance of additive genetic effects, indicated that traits with higher discoverability exhibit more detectable genetic signals. MiXeR’s power analysis was further conducted to determine the proportions of phenotypic variance explained by causal SNPs reaching genome-wide significance at current sample sizes and to estimate the sample sizes necessary to explain larger portions of SNP heritability. Traits that are less genetically discoverable require larger sample sizes to explain a given proportion of SNP heritability.

Bivariate MiXeR analyses28 were then performed to model the estimated additive genetic effects as a mixture of four Gaussian components. These components represent variants that do not influence either trait, variants affecting only one of the traits, and variants influencing both traits. MiXeR simultaneously estimated the proportion of shared trait-influencing variants with concordant effects (i.e., the same effect directions) on both traits. To evaluate the polygenic overlap, MiXeR calculated a Dice coefficient, which is the ratio of shared variants to the total number of variants, providing a measure of the similarity in trait-influencing variants between depression and volume traits. The results for MiXeR are presented in a Venn diagram illustrating unique and shared polygenic components across the traits.

The AIC was employed to assess model fitting performance. In the univariate MiXeR analysis, positive AIC values indicate that the GWAS summary statistics have sufficient power to justify the use of the MiXeR model compared to the LDSC model. In the bivariate MiXeR, a positive AIC value indicates adequate discrimination between the best model fit and the model assuming minimal or maximum possible polygenic overlap28. Notably, SNPs within the major histocompatibility complex (MHC) region (chr6: 25000000-35000000) were excluded from the MiXeR analyses to avoid potential bias from extreme complex LD in this region. For a detailed description of MiXeR, please see the Supplementary Methods.

Cond/conjFDR analyses

CondFDR and conjFDR analyses were employed by “pleiofdr” software (https://github.com/precimed/pleiofdr) to identify novel loci for each trait and specific shared loci for each depression-volume pair, respectively29–31. To enhance the discovery of genetic variants for each individual trait, the test statistics of the primary phenotype were re-ranked by the condFDR statistical framework based on the strength of their association with a secondary phenotype. Additionally, conditional Q‒Q plots were generated to visually depict cross-trait enrichment patterns. These plots present p-value distributions across all SNPs for the primary phenotype and within SNP strata determined by their association significance with the secondary trait (p < 0.1, p < 0.01, and p < 0.001). In these plots, enrichment is observed as an increased degree of leftward deflection from the expected line as the association significance increases in the secondary phenotype. To ensure the model fit, both the MHC and 8p23.1 (chr8: 7200000-12500000) regions were excluded from the conditional Q‒Q plot63,64. Subsequently, conjFDR analysis was performed to identify genomic loci jointly associated with both traits, with the conjFDR value defined as the maximum of the two mutual condFDR values. The conjFDR approach estimates a posterior probability, indicating the likelihood that a SNP is null for either one or both phenotypes, given that the p-values for both phenotypes are lower than the observed p-values. The significance thresholds were set at 0.01 for condFDR and 0.05 for conjFDR, as recommended in previous studies33,65–67. For more details, please see the Supplementary Methods.

Genomic locus definition

For the cond/conjFDR results, independent genomic loci were defined using the “sumstats.py” script (https://github.com/precimed/python_convert), which implements the same logic as in functional mapping and annotation (FUMA, https://fuma.ctglab.nl/)68. Specifically, independent significant SNPs were those that achieved a condFDR <0.01 or conjFDR <0.05 in their respective analyses and were also required to be independent from each other with an r2 < 0.6. Among these significant independent SNPs, lead SNPs were defined as those that were independent of each other with an r2 < 0.1. Candidate SNPs were defined as all SNPs with an r2 ≥ 0.6 with one of the independent significant SNPs, which determined the borders of each locus. If two loci were within 250 kb of each other, they were merged into a single genomic locus. All LD information was calculated from the 1000 Genomes European reference panel. Subsequently, independent loci across all depression-volume pairs were combined into distinct loci. In cases where the genetic positions of two or more loci from different depression-volume pairs overlapped, they were merged into one distinct locus with the union of their genomic boundaries. The lead SNP with the minimum FDR value was designated the final lead SNP in the merged locus.

To determine the novelty of independent loci identified by cond/conjFDR, several steps were taken. First, we verified whether any candidate SNPs within the identified loci had been reported to be significantly associated with depression or subcortical volumetric traits in the NHGRI-EBI GWAS Catalog69. Next, the identified loci were compared with significant loci (p < 1.0E−6) reported in the original GWAS20,23 and findings from previous cond/conjFDR studies33,70–81 or other GWASs44,73,82,83 (Supplementary Data 21, 22). A locus was considered novel if all candidate SNPs had not been previously reported in the NHGRI-EBI GWAS Catalog and if it extended beyond 500 kb from the borders of any reported loci in the original GWAS summary, previous cond/conjFDR studies, or other GWAS studies. Notably, if the start and end positions of the locus were not provided, the locus border was defined as 1000 kb away from the lead SNP.

Functional annotation

The FUMA online platform was used to perform functional annotation for candidate SNPs within specific genomic loci. These SNPs were annotated using CADD scores, RegulomeDB scores, and 15-core chromatin states. CADD scores, derived from over 60 genomic features, predict the likelihood of a variant being pathogenic, with scores above 12.37 suggesting a deleterious effect84. RegulomeDB scores assess the regulatory functionality of SNPs, inferred from the eQTL and chromatin markers—a lower score indicates a greater likelihood of regulatory functionality85. The 15-core chromatin state classification, based on ChromHMM predictions and utilizing multiple chromatin markers across 127 epigenomes, denotes various functional states of the genome, including promoters, enhancers, and transcriptional activity regions. Categories 1–7 indicate open chromatin states, which are associated with transcriptional regulation and gene expression86,87. Moreover, candidate SNPs were mapped to protein-coding genes using any of three strategies: (1) positional mapping within a 10-kb window, (2) eQTL mapping by associating cis-eQTL SNPs with genes whose expression varies with allelic variation at the SNP level (GTEx v8 brain), and (3) chromatin interaction mapping to establish connections through three-dimensional DNA‒DNA interactions between the genomic regions of SNPs and proximal or distal genes. Furthermore, an assessment of the enrichment of GO biological processes for all mapped genes was conducted using g:Profiler (https://biit.cs.ut.ee/gprofiler)88, and the Benjamini‒Hochberg FDR method was applied to correct for multiple comparisons, with a significance threshold of q < 0.05.

Heterogeneity of shared loci across subcortical volumetric phenotypes

To assess the differences and similarities in how shared genetic variants affect different subcortical volumetric phenotypes, the effect sizes and standard errors of 44 lead SNPs for 14 subcortical volumetric phenotypes were extracted from the original GWAS summary statistics. Cochran’s Q test, a statistical method designed to detect heterogeneity, was then applied to these statistical values. Specifically, two types of heterogeneity tests were performed: one assessing all 14 subcortical phenotypes together and the other conducting pairwise comparisons between each pair of the 14 subcortical phenotypes. A p-value of less than 0.05 from Cochran’s Q test was considered to indicate significant heterogeneity among the subcortical volumetric phenotypes for a given SNP.

Spatiotemporal gene expression trajectory analysis

Spatiotemporal brain expression trajectory analysis was performed to explore how the expression of genes shared between depression and subcortical volumetric phenotypes varies across different developmental stages, aiming to elucidate the potential roles these genes play in the pathophysiology of depression and associated structural changes in the brain. Processed mRNA-seq data from PsychENCODE35, comprising 607 tissue samples from 41 postmortem brains ranging from 8 PCW to 40 PY, were utilized to identify the expression patterns of these shared genes across lifespan windows. Following the methodology of the PsychENCODE study35, we divided the data into nine developmental windows: 8-9 PCW, 12-13 PCW, 16-17 PCW, 19-22 PCW, 35 PCW to 4 months, 0.5-2.6 PY, 2.8-10.7 PY, 13-19 PY, and 21-40 PY. Gene expression levels, represented in RPKM, were log2-transformed and centered to the mean expression level for each sample89. The expression values for each volumetric phenotype were then calculated as the mean values across all shared genes mapped by their shared loci with depression. Finally, smoothed LOESS curves were plotted to visualize the expression trajectories.

Expression-trait association analysis

To evaluate whether the brain’s expression levels of shared genes, identified through conjFDR analyses of each depression-volume pair, are significantly correlated with these traits (i.e., depression and corresponding volumetric phenotypes), we employed the “S-MultiXcan” tool (https://github.com/hakyimlab/MetaXcan)36. This method integrates GWAS summary data with pre-computed eQTL prediction models of all the brain tissues available from the GTEx v8, including the amygdala, anterior cingulate cortex (BA24), caudate, cerebellar hemisphere, cerebellum, cortex, frontal cortex (BA9), hippocampus, hypothalamus, nucleus accumbens, putamen, cervical spinal cord C1, and substantia nigra. Genes lacking significant brain eQTL associations were excluded from the analysis, and a significant expression-trait association was determined with an FDR corrected q < 0.05.

Cognitive and behavioral associations

To explore whether the expression of genes shared between depression and subcortical brain regions influences specific cognitive or behavioral functions, we utilized GWAS summary data for 8 cognitive abilities90 and 15 behavioral symptoms. The cognitive abilities included digit symbol substitution task (DSST), VNR, pairs matching (PM), RT, matrix pattern recognition (Matrix), tower rearranging (Tower), trail-making test B (TMT-B), and a general factor of intelligence (g-factor), as determined by a genomic structural equation model. These cognitive ability domains cover reasoning ability, executive function, processing speed, and memory, offering a comprehensive view of neurocognitive performance (for more details, see Supplementary Data 17). Additionally, the GWAS summary for behavioral symptoms included appetite change, sleep disturbance, and difficulty concentrating, among others (for more details, see Supplementary Data 18). Using the “S-MultiXcan” tool, we performed association analyses between the expression of shared genes—those whose expression was linked to both depression and at least one subcortical volume—and these cognitive abilities and behavioral symptoms, with significance determined by an FDR-corrected q < 0.05.

Control analyses

Genome-wide genetic correlation

The LDSC analysis was performed using “ldsc” software (v1.0.1, https://github.com/bulik/ldsc)37,38 to calculate genome-wide genetic correlations for each depression-volume pair based on the summary statistics. All analyses were exclusively performed on HapMap3 variants, excluding the MHC and 8p23.1 regions, and utilized pre-calculated LD scores from the 1000 Genomes European data.

Local genetic correlation

LAVA (v0.1.0, https://github.com/josefin-werme/LAVA)39 was employed to validate the shared loci identified by conjFDR analysis. As the detection of valid and interpretable local genetic correlations depends on the presence of sufficient local genetic signals, univariate tests in the LAVA approach served as a filtering step for bivariate local genetic correlation analyses. Bivariate local genetic correlation analyses were then exclusively performed for pairs of depression and volumetric traits with a nominally significant univariate local genetic signal.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Peer Review File

Description of additional supplementary files

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Supplementary Data 4

Supplementary Data 5

Supplementary Data 6

Supplementary Data 7

Supplementary Data 8

Supplementary Data 9

Supplementary Data 10

Supplementary Data 11

Supplementary Data 12

Supplementary Data 13

Supplementary Data 14

Supplementary Data 15

Supplementary Data 16

Supplementary Data 17

Supplementary Data 18

Supplementary Data 19

Supplementary Data 20

Supplementary Data 21

Supplementary Data 22

Reporting Summary

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52121-y.

Acknowledgements

This work was funded by the Natural Science Foundation of China (82072001, 82102318), the Beijing-Tianjin-Hebei Basic Research Collaboration Project (J230040), the Tianjin Natural Science Foundation (23JCZXJC00120), and the Tianjin Key Medical Discipline (Specialty) Construction Project (TJYXZDXK-001A) to F.L.

Author contributions

F.L. and M.Liu designed the study and wrote the manuscript. M.Liu, L.W. and Y.Z. analyzed data. F.L., Q.Z. and Sijia Wang supervised this work. H.D. and C.W. acquired and processed the data. Y.C., Q.Q. and N.Z. advised on the MiXeR and cond/conjFDR analyses. Shaoying Wang and G.Z. contributed to the expression-trait analyses. Z.Z. and M.Lei advised on the expression trajectory analyses. All authors critically reviewed the manuscript.

Peer review

Peer review information

Nature Communications thanks Yao Wu, and the other, anonymous, reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

All datasets used in this study are publicly accessible. Specifically, the GWAS summary statistics for the 14 subcortical volumes based on UK biobank participants are available at https://open.win.ox.ac.uk/ukbiobank/big40/, and GWAS summary statistics for depression can be downloaded from https://ipsych.dk/en/research/downloads/. The processed mRNA-seq data were obtained from http://development.psychencode.org/. Additionally, the corresponding accessible links for GWAS summary statistics for 8 cognitive abilities and 15 behavioral symptoms are provided in Supplementary Data 17, 18. Source data are provided with this paper.

Code availability

The custom code that supports the findings of this study is available at 10.5281/zenodo.12912958 (ref. 91).

Competing interests

The authors declare no competing interests.

Ethics

Our research complies with all relevant ethical regulations. All GWASs used in the current study were approved by the ethics committees of the original study, and informed consent was obtained from all participants.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Mengge Liu, Lu Wang, Yujie Zhang.
==== Refs
References

1. Disease GBD, Injury I, Prevalence C. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017 Lancet 2018 392 1789 1858 10.1016/S0140-6736(18)32279-7 30496104
Disease GBD, Injury I, Prevalence C. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 392, 1789–1858 (2018).30496104 10.1016/S0140-6736(18)32279-7
2. Kaplan KA Harvey AG Hypersomnia across mood disorders: a review and synthesis Sleep. Med Rev. 2009 13 275 285 10.1016/j.smrv.2008.09.001 19269201
Kaplan, K. A. & Harvey, A. G. Hypersomnia across mood disorders: a review and synthesis. Sleep. Med Rev. 13, 275–285 (2009).19269201 10.1016/j.smrv.2008.09.001
3. Baglioni C Insomnia as a predictor of depression: a meta-analytic evaluation of longitudinal epidemiological studies J. Affect. Disord. 2011 135 10 19 10.1016/j.jad.2011.01.011 21300408
Baglioni, C. et al. Insomnia as a predictor of depression: a meta-analytic evaluation of longitudinal epidemiological studies. J. Affect. Disord. 135, 10–19 (2011).21300408 10.1016/j.jad.2011.01.011
4. Maxwell MA Cole DA Weight change and appetite disturbance as symptoms of adolescent depression: toward an integrative biopsychosocial model Clin. Psychol. Rev. 2009 29 260 273 10.1016/j.cpr.2009.01.007 19250729
Maxwell, M. A. & Cole, D. A. Weight change and appetite disturbance as symptoms of adolescent depression: toward an integrative biopsychosocial model. Clin. Psychol. Rev. 29, 260–273 (2009).19250729 10.1016/j.cpr.2009.01.007
5. Conradi HJ Ormel J de Jonge P Presence of individual (residual) symptoms during depressive episodes and periods of remission: a 3-year prospective study Psychol. Med 2011 41 1165 1174 10.1017/S0033291710001911 20932356
Conradi, H. J., Ormel, J. & de Jonge, P. Presence of individual (residual) symptoms during depressive episodes and periods of remission: a 3-year prospective study. Psychol. Med 41, 1165–1174 (2011).20932356 10.1017/S0033291710001911
6. Millan MJ Cognitive dysfunction in psychiatric disorders: characteristics, causes and the quest for improved therapy Nat. Rev. Drug Discov. 2012 11 141 168 10.1038/nrd3628 22293568
Millan, M. J. et al. Cognitive dysfunction in psychiatric disorders: characteristics, causes and the quest for improved therapy. Nat. Rev. Drug Discov. 11, 141–168 (2012).22293568 10.1038/nrd3628
7. Semkovska M Cognitive function following a major depressive episode: A systematic review and meta-analysis Lancet Psychiatry 2019 6 851 861 10.1016/S2215-0366(19)30291-3 31422920
Semkovska, M. et al. Cognitive function following a major depressive episode: A systematic review and meta-analysis. Lancet Psychiatry 6, 851–861 (2019).31422920 10.1016/S2215-0366(19)30291-3
8. Hamilton JP Siemer M Gotlib IH Amygdala volume in major depressive disorder: A meta-analysis of magnetic resonance imaging studies Mol. Psychiatry 2008 13 993 1000 10.1038/mp.2008.57 18504424
Hamilton, J. P., Siemer, M. & Gotlib, I. H. Amygdala volume in major depressive disorder: A meta-analysis of magnetic resonance imaging studies. Mol. Psychiatry 13, 993–1000 (2008).18504424 10.1038/mp.2008.57
9. Saleh K Impact of family history and depression on amygdala volume Psychiatry Res 2012 203 24 30 10.1016/j.pscychresns.2011.10.004 22867951
Saleh, K. et al. Impact of family history and depression on amygdala volume. Psychiatry Res 203, 24–30 (2012).22867951 10.1016/j.pscychresns.2011.10.004
10. Sandu AL Amygdala and regional volumes in treatment-resistant versus nontreatment-resistant depression patients Depress. Anxiety 2017 34 1065 1071 10.1002/da.22675 28792656
Sandu, A. L. et al. Amygdala and regional volumes in treatment-resistant versus nontreatment-resistant depression patients. Depress. Anxiety 34, 1065–1071 (2017).28792656 10.1002/da.22675
11. Bylsma LM Morris BH Rottenberg J A meta-analysis of emotional reactivity in major depressive disorder Clin. Psychol. Rev. 2008 28 676 691 10.1016/j.cpr.2007.10.001 18006196
Bylsma, L. M., Morris, B. H. & Rottenberg, J. A meta-analysis of emotional reactivity in major depressive disorder. Clin. Psychol. Rev. 28, 676–691 (2008).18006196 10.1016/j.cpr.2007.10.001
12. Schmaal L Subcortical brain alterations in major depressive disorder: findings from the ENIGMA Major Depressive Disorder working group Mol. Psychiatry 2016 21 806 812 10.1038/mp.2015.69 26122586
Schmaal, L. et al. Subcortical brain alterations in major depressive disorder: findings from the ENIGMA Major Depressive Disorder working group. Mol. Psychiatry 21, 806–812 (2016).26122586 10.1038/mp.2015.69
13. Videbech P Ravnkilde B Hippocampal volume and depression: A meta-analysis of MRI studies Am. J. Psychiatry 2004 161 1957 1966 10.1176/appi.ajp.161.11.1957 15514393
Videbech, P. & Ravnkilde, B. Hippocampal volume and depression: A meta-analysis of MRI studies. Am. J. Psychiatry 161, 1957–1966 (2004).15514393 10.1176/appi.ajp.161.11.1957
14. Whittle S Structural brain development and depression onset during adolescence: a prospective longitudinal study Am. J. Psychiatry 2014 171 564 571 10.1176/appi.ajp.2013.13070920 24577365
Whittle, S. et al. Structural brain development and depression onset during adolescence: a prospective longitudinal study. Am. J. Psychiatry 171, 564–571 (2014).24577365 10.1176/appi.ajp.2013.13070920
15. Hickie I Reduced hippocampal volumes and memory loss in patients with early- and late-onset depression Br. J. Psychiatry 2005 186 197 202 10.1192/bjp.186.3.197 15738499
Hickie, I. et al. Reduced hippocampal volumes and memory loss in patients with early- and late-onset depression. Br. J. Psychiatry 186, 197–202 (2005).15738499 10.1192/bjp.186.3.197
16. Auerbach RP Reward-related neural circuitry in depressed and anxious adolescents: A Human Connectome Project J. Am. Acad. Child Adolesc. Psychiatry 2022 61 308 320 10.1016/j.jaac.2021.04.014 33965516
Auerbach, R. P. et al. Reward-related neural circuitry in depressed and anxious adolescents: A Human Connectome Project. J. Am. Acad. Child Adolesc. Psychiatry 61, 308–320 (2022).33965516 10.1016/j.jaac.2021.04.014
17. Nugent AC Davis RM Zarate CA Jr. Drevets WC Reduced thalamic volumes in major depressive disorder Psychiatry Res 2013 213 179 185 10.1016/j.pscychresns.2013.05.004 23850106
Nugent, A. C., Davis, R. M., Zarate, C. A. Jr. & Drevets, W. C. Reduced thalamic volumes in major depressive disorder. Psychiatry Res 213, 179–185 (2013).23850106 10.1016/j.pscychresns.2013.05.004
18. Talati A Putamen structure and function in familial risk for depression: A multimodal imaging study Biol. Psychiatry 2022 92 932 941 10.1016/j.biopsych.2022.06.035 36038379
Talati, A. et al. Putamen structure and function in familial risk for depression: A multimodal imaging study. Biol. Psychiatry 92, 932–941 (2022).36038379 10.1016/j.biopsych.2022.06.035
19. Koolschijn PC van Haren NE Lensvelt-Mulders GJ Hulshoff Pol HE Kahn RS Brain volume abnormalities in major depressive disorder: a meta-analysis of magnetic resonance imaging studies Hum. Brain Mapp. 2009 30 3719 3735 10.1002/hbm.20801 19441021
Koolschijn, P. C., van Haren, N. E., Lensvelt-Mulders, G. J., Hulshoff Pol, H. E. & Kahn, R. S. Brain volume abnormalities in major depressive disorder: a meta-analysis of magnetic resonance imaging studies. Hum. Brain Mapp. 30, 3719–3735 (2009).19441021 10.1002/hbm.20801
20. Als TD Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses Nat. Med 2023 29 1832 1844 10.1038/s41591-023-02352-1 37464041
Als, T. D. et al. Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses. Nat. Med 29, 1832–1844 (2023).37464041 10.1038/s41591-023-02352-1
21. Zhao B Genome-wide association analysis of 19,629 individuals identifies variants influencing regional brain volumes and refines their genetic co-architecture with cognitive and mental health traits Nat. Genet 2019 51 1637 1644 10.1038/s41588-019-0516-6 31676860
Zhao, B. et al. Genome-wide association analysis of 19,629 individuals identifies variants influencing regional brain volumes and refines their genetic co-architecture with cognitive and mental health traits. Nat. Genet 51, 1637–1644 (2019).31676860 10.1038/s41588-019-0516-6
22. Elliott LT Genome-wide association studies of brain imaging phenotypes in UK Biobank Nature 2018 562 210 216 10.1038/s41586-018-0571-7 30305740
Elliott, L. T. et al. Genome-wide association studies of brain imaging phenotypes in UK Biobank. Nature 562, 210–216 (2018).30305740 10.1038/s41586-018-0571-7
23. Smith SM An expanded set of genome-wide association studies of brain imaging phenotypes in UK Biobank Nat. Neurosci. 2021 24 737 745 10.1038/s41593-021-00826-4 33875891
Smith, S. M. et al. An expanded set of genome-wide association studies of brain imaging phenotypes in UK Biobank. Nat. Neurosci. 24, 737–745 (2021).33875891 10.1038/s41593-021-00826-4
24. Sullivan PF Neale MC Kendler KS Genetic epidemiology of major depression: review and meta-analysis Am. J. Psychiatry 2000 157 1552 1562 10.1176/appi.ajp.157.10.1552 11007705
Sullivan, P. F., Neale, M. C. & Kendler, K. S. Genetic epidemiology of major depression: review and meta-analysis. Am. J. Psychiatry 157, 1552–1562 (2000).11007705 10.1176/appi.ajp.157.10.1552
25. Ohi K Genetic correlations between subcortical brain volumes and psychiatric disorders Br. J. Psychiatry 2020 216 280 283 10.1192/bjp.2019.277 32000869
Ohi, K. et al. Genetic correlations between subcortical brain volumes and psychiatric disorders. Br. J. Psychiatry 216, 280–283 (2020).32000869 10.1192/bjp.2019.277
26. Wigmore EM Do regional brain volumes and major depressive disorder share genetic architecture? A study of Generation Scotland (n = 19762), UK Biobank (n = 24 048) and the English Longitudinal Study of Ageing (n = 5766) Transl. Psychiatry 2017 7 e1205 10.1038/tp.2017.148 28809859
Wigmore, E. M. et al. Do regional brain volumes and major depressive disorder share genetic architecture? A study of Generation Scotland (n = 19762), UK Biobank (n = 24 048) and the English Longitudinal Study of Ageing (n = 5766). Transl. Psychiatry 7, e1205 (2017).28809859 10.1038/tp.2017.148
27. Cullen, H. et al. Common genetic variation important in early subcortical brain development. medRxiv, (2022).
28. Frei O Bivariate causal mixture model quantifies polygenic overlap between complex traits beyond genetic correlation Nat. Commun. 2019 10 2417 10.1038/s41467-019-10310-0 31160569
Frei, O. et al. Bivariate causal mixture model quantifies polygenic overlap between complex traits beyond genetic correlation. Nat. Commun. 10, 2417 (2019).31160569 10.1038/s41467-019-10310-0
29. Andreassen OA Thompson WK Dale AM Boosting the power of schizophrenia genetics by leveraging new statistical tools Schizophr. Bull. 2014 40 13 17 10.1093/schbul/sbt168 24319118
Andreassen, O. A., Thompson, W. K. & Dale, A. M. Boosting the power of schizophrenia genetics by leveraging new statistical tools. Schizophr. Bull. 40, 13–17 (2014).24319118 10.1093/schbul/sbt168
30. Andreassen OA Improved detection of common variants associated with schizophrenia by leveraging pleiotropy with cardiovascular-disease risk factors Am. J. Hum. Genet 2013 92 197 209 10.1016/j.ajhg.2013.01.001 23375658
Andreassen, O. A. et al. Improved detection of common variants associated with schizophrenia by leveraging pleiotropy with cardiovascular-disease risk factors. Am. J. Hum. Genet 92, 197–209 (2013).23375658 10.1016/j.ajhg.2013.01.001
31. Andreassen OA Improved detection of common variants associated with schizophrenia and bipolar disorder using pleiotropy-informed conditional false discovery rate PLoS Genet 2013 9 e1003455 10.1371/journal.pgen.1003455 23637625
Andreassen, O. A. et al. Improved detection of common variants associated with schizophrenia and bipolar disorder using pleiotropy-informed conditional false discovery rate. PLoS Genet 9, e1003455 (2013).23637625 10.1371/journal.pgen.1003455
32. Karch CM Selective genetic overlap between amyotrophic lateral sclerosis and diseases of the frontotemporal dementia spectrum JAMA Neurol. 2018 75 860 875 10.1001/jamaneurol.2018.0372 29630712
Karch, C. M. et al. Selective genetic overlap between amyotrophic lateral sclerosis and diseases of the frontotemporal dementia spectrum. JAMA Neurol. 75, 860–875 (2018).29630712 10.1001/jamaneurol.2018.0372
33. Bahrami S Genetic loci shared between major depression and intelligence with mixed directions of effect Nat. Hum. Behav. 2021 5 795 801 10.1038/s41562-020-01031-2 33462475
Bahrami, S. et al. Genetic loci shared between major depression and intelligence with mixed directions of effect. Nat. Hum. Behav. 5, 795–801 (2021).33462475 10.1038/s41562-020-01031-2
34. Holland D Beyond SNP heritability: Polygenicity and discoverability of phenotypes estimated with a univariate Gaussian mixture model PLoS Genet. 2020 16 e1008612 10.1371/journal.pgen.1008612 32427991
Holland, D. et al. Beyond SNP heritability: Polygenicity and discoverability of phenotypes estimated with a univariate Gaussian mixture model. PLoS Genet. 16, e1008612 (2020).32427991 10.1371/journal.pgen.1008612
35. Li M Integrative functional genomic analysis of human brain development and neuropsychiatric risks Science 2018 362 eaat7615 10.1126/science.aat7615 30545854
Li, M. et al. Integrative functional genomic analysis of human brain development and neuropsychiatric risks. Science 362, eaat7615 (2018).30545854 10.1126/science.aat7615
36. Barbeira AN Integrating predicted transcriptome from multiple tissues improves association detection PLoS Genet. 2019 15 e1007889 10.1371/journal.pgen.1007889 30668570
Barbeira, A. N. et al. Integrating predicted transcriptome from multiple tissues improves association detection. PLoS Genet. 15, e1007889 (2019).30668570 10.1371/journal.pgen.1007889
37. Bulik-Sullivan BK LD Score regression distinguishes confounding from polygenicity in genome-wide association studies Nat. Genet. 2015 47 291 295 10.1038/ng.3211 25642630
Bulik-Sullivan, B. K. et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat. Genet. 47, 291–295 (2015).25642630 10.1038/ng.3211
38. Bulik-Sullivan B An atlas of genetic correlations across human diseases and traits Nat. Genet. 2015 47 1236 1241 10.1038/ng.3406 26414676
Bulik-Sullivan, B. et al. An atlas of genetic correlations across human diseases and traits. Nat. Genet. 47, 1236–1241 (2015).26414676 10.1038/ng.3406
39. Werme J van der Sluis S Posthuma D de Leeuw CA An integrated framework for local genetic correlation analysis Nat. Genet. 2022 54 274 282 10.1038/s41588-022-01017-y 35288712
Werme, J., van der Sluis, S., Posthuma, D. & de Leeuw, C. A. An integrated framework for local genetic correlation analysis. Nat. Genet. 54, 274–282 (2022).35288712 10.1038/s41588-022-01017-y
40. Le Magueresse C Monyer H GABAergic interneurons shape the functional maturation of the cortex Neuron 2013 77 388 405 10.1016/j.neuron.2013.01.011 23395369
Le Magueresse, C. & Monyer, H. GABAergic interneurons shape the functional maturation of the cortex. Neuron 77, 388–405 (2013).23395369 10.1016/j.neuron.2013.01.011
41. Kelsom C Lu W Development and specification of GABAergic cortical interneurons Cell Biosci. 2013 3 19 10.1186/2045-3701-3-19 23618463
Kelsom, C. & Lu, W. Development and specification of GABAergic cortical interneurons. Cell Biosci. 3, 19 (2013).23618463 10.1186/2045-3701-3-19
42. Coghlan S GABA system dysfunction in autism and related disorders: from synapse to symptoms Neurosci. Biobehav. Rev. 2012 36 2044 2055 10.1016/j.neubiorev.2012.07.005 22841562
Coghlan, S. et al. GABA system dysfunction in autism and related disorders: from synapse to symptoms. Neurosci. Biobehav. Rev. 36, 2044–2055 (2012).22841562 10.1016/j.neubiorev.2012.07.005
43. Xu MY Wong AHC GABAergic inhibitory neurons as therapeutic targets for cognitive impairment in schizophrenia Acta Pharm. Sin. 2018 39 733 753 10.1038/aps.2017.172
Xu, M. Y. & Wong, A. H. C. GABAergic inhibitory neurons as therapeutic targets for cognitive impairment in schizophrenia. Acta Pharm. Sin. 39, 733–753 (2018).10.1038/aps.2017.172
44. Hibar DP Common genetic variants influence human subcortical brain structures Nature 2015 520 224 229 10.1038/nature14101 25607358
Hibar, D. P. et al. Common genetic variants influence human subcortical brain structures. Nature 520, 224–229 (2015).25607358 10.1038/nature14101
45. Torres-Berrío A DCC confers susceptibility to depression-like behaviors in humans and mice and is regulated by miR-218 Biol. Psychiatry 2017 81 306 315 10.1016/j.biopsych.2016.08.017 27773352
Torres-Berrío, A. et al. DCC confers susceptibility to depression-like behaviors in humans and mice and is regulated by miR-218. Biol. Psychiatry 81, 306–315 (2017).27773352 10.1016/j.biopsych.2016.08.017
46. Torres-Berrío A Hernandez G Nestler EJ Flores C The Netrin-1/DCC Guidance Cue pathway as a molecular target in depression: translational evidence Biol. Psychiatry 2020 88 611 624 10.1016/j.biopsych.2020.04.025 32593422
Torres-Berrío, A., Hernandez, G., Nestler, E. J. & Flores, C. The Netrin-1/DCC Guidance Cue pathway as a molecular target in depression: translational evidence. Biol. Psychiatry 88, 611–624 (2020).32593422 10.1016/j.biopsych.2020.04.025
47. Grant A Fathalli F Rouleau G Joober R Flores C Association between schizophrenia and genetic variation in DCC: a case-control study Schizophr. Res. 2012 137 26 31 10.1016/j.schres.2012.02.023 22418395
Grant, A., Fathalli, F., Rouleau, G., Joober, R. & Flores, C. Association between schizophrenia and genetic variation in DCC: a case-control study. Schizophr. Res. 137, 26–31 (2012).22418395 10.1016/j.schres.2012.02.023
48. Aberg KA A comprehensive family-based replication study of schizophrenia genes JAMA Psychiatry 2013 70 573 581 10.1001/jamapsychiatry.2013.288 23894747
Aberg, K. A. et al. A comprehensive family-based replication study of schizophrenia genes. JAMA Psychiatry 70, 573–581 (2013).23894747 10.1001/jamapsychiatry.2013.288
49. Arnau-Soler A Adams MJ Scotland G Consortium MDDWGotPG, Hayward C, Thomson PA. Genome-wide interaction study of a proxy for stress-sensitivity and its prediction of major depressive disorder PloS One 2018 13 e0209160 10.1371/journal.pone.0209160 30571770
Arnau-Soler, A., Adams, M. J. & Scotland, G. Consortium MDDWGotPG, Hayward C, Thomson PA. Genome-wide interaction study of a proxy for stress-sensitivity and its prediction of major depressive disorder. PloS One 13, e0209160 (2018).30571770 10.1371/journal.pone.0209160
50. Wendt FR Multivariate genome-wide analysis of education, socioeconomic status and brain phenome Nat. Hum. Behav. 2021 5 482 496 10.1038/s41562-020-00980-y 33349686
Wendt, F. R. et al. Multivariate genome-wide analysis of education, socioeconomic status and brain phenome. Nat. Hum. Behav. 5, 482–496 (2021).33349686 10.1038/s41562-020-00980-y
51. Hofer E Genetic correlations and genome-wide associations of cortical structure in general population samples of 22,824 adults Nat. Commun. 2020 11 4796 10.1038/s41467-020-18367-y 32963231
Hofer, E. et al. Genetic correlations and genome-wide associations of cortical structure in general population samples of 22,824 adults. Nat. Commun. 11, 4796 (2020).32963231 10.1038/s41467-020-18367-y
52. Tian Y Shared genetics and comorbid genes of amyotrophic lateral sclerosis and Parkinson’s disease Mov. Disord. 2023 38 1813 1821 10.1002/mds.29572 37534731
Tian, Y. et al. Shared genetics and comorbid genes of amyotrophic lateral sclerosis and Parkinson’s disease. Mov. Disord. 38, 1813–1821 (2023).37534731 10.1002/mds.29572
53. Yao S A transcriptome-wide association study identifies susceptibility genes for Parkinson’s disease NPJ Parkinsons Dis. 2021 7 79 10.1038/s41531-021-00221-7 34504106
Yao, S. et al. A transcriptome-wide association study identifies susceptibility genes for Parkinson’s disease. NPJ Parkinsons Dis. 7, 79 (2021).34504106 10.1038/s41531-021-00221-7
54. Claes SJ Corticotropin-releasing hormone (CRH) in psychiatry: From stress to psychopathology Ann. Med 2004 36 50 61 10.1080/07853890310017044 15000347
Claes, S. J. Corticotropin-releasing hormone (CRH) in psychiatry: From stress to psychopathology. Ann. Med 36, 50–61 (2004).15000347 10.1080/07853890310017044
55. Tyrka AR Interaction of childhood maltreatment with the corticotropin-releasing hormone receptor gene: effects on hypothalamic-pituitary-adrenal axis reactivity Biol. Psychiatry 2009 66 681 685 10.1016/j.biopsych.2009.05.012 19596121
Tyrka, A. R. et al. Interaction of childhood maltreatment with the corticotropin-releasing hormone receptor gene: effects on hypothalamic-pituitary-adrenal axis reactivity. Biol. Psychiatry 66, 681–685 (2009).19596121 10.1016/j.biopsych.2009.05.012
56. Smoller JW The corticotropin-releasing hormone gene and behavioral inhibition in children at risk for panic disorder Biol. Psychiatry 2005 57 1485 1492 10.1016/j.biopsych.2005.02.018 15953484
Smoller, J. W. et al. The corticotropin-releasing hormone gene and behavioral inhibition in children at risk for panic disorder. Biol. Psychiatry 57, 1485–1492 (2005).15953484 10.1016/j.biopsych.2005.02.018
57. Rogers J CRHR1 genotypes, neural circuits and the diathesis for anxiety and depression Mol. Psychiatry 2013 18 700 707 10.1038/mp.2012.152 23147386
Rogers, J. et al. CRHR1 genotypes, neural circuits and the diathesis for anxiety and depression. Mol. Psychiatry 18, 700–707 (2013).23147386 10.1038/mp.2012.152
58. Allen M Gene expression, methylation and neuropathology correlations at progressive supranuclear palsy risk loci Acta Neuropathol. 2016 132 197 211 10.1007/s00401-016-1576-7 27115769
Allen, M. et al. Gene expression, methylation and neuropathology correlations at progressive supranuclear palsy risk loci. Acta Neuropathol. 132, 197–211 (2016).27115769 10.1007/s00401-016-1576-7
59. Johnson LA The impact of GPX1 on the association of groundwater selenium and depression: a Project FRONTIER study BMC Psychiatry 2013 13 7 10.1186/1471-244X-13-7 23289525
Johnson, L. A. et al. The impact of GPX1 on the association of groundwater selenium and depression: a Project FRONTIER study. BMC Psychiatry 13, 7 (2013).23289525 10.1186/1471-244X-13-7
60. Li X Transcriptome-wide association study identifies new susceptibility genes and pathways for depression Transl. Psychiatry 2021 11 306 10.1038/s41398-021-01411-w 34021117
Li, X. et al. Transcriptome-wide association study identifies new susceptibility genes and pathways for depression. Transl. Psychiatry 11, 306 (2021).34021117 10.1038/s41398-021-01411-w
61. Amin N Exome-sequencing in a large population-based study reveals a rare Asn396Ser variant in the LIPG gene associated with depressive symptoms Mol. Psychiatry 2017 22 634 10.1038/mp.2016.141 27502476
Amin, N. et al. Exome-sequencing in a large population-based study reveals a rare Asn396Ser variant in the LIPG gene associated with depressive symptoms. Mol. Psychiatry 22, 634 (2017).27502476 10.1038/mp.2016.141
62. Cheng S Exome-wide screening identifies novel rare risk variants for major depression disorder Mol. Psychiatry 2022 27 3069 3074 10.1038/s41380-022-01536-4 35365804
Cheng, S. et al. Exome-wide screening identifies novel rare risk variants for major depression disorder. Mol. Psychiatry 27, 3069–3074 (2022).35365804 10.1038/s41380-022-01536-4
63. Smeland OB Genome-wide association analysis of Parkinson’s disease and Schizophrenia reveals shared genetic architecture and identifies novel risk Loci Biol. Psychiatry 2021 89 227 235 10.1016/j.biopsych.2020.01.026 32201043
Smeland, O. B. et al. Genome-wide association analysis of Parkinson’s disease and Schizophrenia reveals shared genetic architecture and identifies novel risk Loci. Biol. Psychiatry 89, 227–235 (2021).32201043 10.1016/j.biopsych.2020.01.026
64. Hope S Bidirectional genetic overlap between autism spectrum disorder and cognitive traits Transl. Psychiatry 2023 13 295 10.1038/s41398-023-02563-7 37709755
Hope, S. et al. Bidirectional genetic overlap between autism spectrum disorder and cognitive traits. Transl. Psychiatry 13, 295 (2023).37709755 10.1038/s41398-023-02563-7
65. Ma DR Shared genetic architecture between Parkinson’s disease and brain structural phenotypes Mov. Disord. 2023 38 2258 2268 10.1002/mds.29598 37990409
Ma, D. R. et al. Shared genetic architecture between Parkinson’s disease and brain structural phenotypes. Mov. Disord. 38, 2258–2268 (2023).37990409 10.1002/mds.29598
66. O’Connell KS Characterizing the genetic overlap between psychiatric disorders and sleep-related phenotypes Biol. Psychiatry 2021 90 621 631 10.1016/j.biopsych.2021.07.007 34482950
O’Connell, K. S. et al. Characterizing the genetic overlap between psychiatric disorders and sleep-related phenotypes. Biol. Psychiatry 90, 621–631 (2021).34482950 10.1016/j.biopsych.2021.07.007
67. Hindley G Charting the landscape of genetic overlap between mental disorders and related traits beyond genetic correlation Am. J. Psychiatry 2022 179 833 843 10.1176/appi.ajp.21101051 36069018
Hindley, G. et al. Charting the landscape of genetic overlap between mental disorders and related traits beyond genetic correlation. Am. J. Psychiatry 179, 833–843 (2022).36069018 10.1176/appi.ajp.21101051
68. Watanabe K Taskesen E Van Bochoven A Posthuma D Functional mapping and annotation of genetic associations with FUMA Nat. Commun. 2017 8 1826 10.1038/s41467-017-01261-5 29184056
Watanabe, K., Taskesen, E., Van Bochoven, A. & Posthuma, D. Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 8, 1826 (2017).29184056 10.1038/s41467-017-01261-5
69. Buniello A The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019 Nucleic Acids Res. 2019 47 D1005 D1012 10.1093/nar/gky1120 30445434
Buniello, A. et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 47, D1005–D1012 (2019).30445434 10.1093/nar/gky1120
70. Bahrami S Dissecting the shared genetic basis of migraine and mental disorders using novel statistical tools Brain 2022 145 142 153 10.1093/brain/awab267 34273149
Bahrami, S. et al. Dissecting the shared genetic basis of migraine and mental disorders using novel statistical tools. Brain 145, 142–153 (2022).34273149 10.1093/brain/awab267
71. Bahrami S Shared genetic Loci Between Body Mass Index and major psychiatric disorders: a genome-wide association study JAMA Psychiatry 2020 77 503 512 10.1001/jamapsychiatry.2019.4188 31913414
Bahrami, S. et al. Shared genetic Loci Between Body Mass Index and major psychiatric disorders: a genome-wide association study. JAMA Psychiatry 77, 503–512 (2020).31913414 10.1001/jamapsychiatry.2019.4188
72. Bang L Genome-wide analysis of anorexia nervosa and major psychiatric disorders and related traits reveals genetic overlap and identifies novel risk loci for anorexia nervosa Transl. Psychiatry 2023 13 291 10.1038/s41398-023-02585-1 37658054
Bang, L. et al. Genome-wide analysis of anorexia nervosa and major psychiatric disorders and related traits reveals genetic overlap and identifies novel risk loci for anorexia nervosa. Transl. Psychiatry 13, 291 (2023).37658054 10.1038/s41398-023-02585-1
73. Elvsåshagen T The genetic architecture of the human thalamus and its overlap with ten common brain disorders Nat. Commun. 2021 12 2909 10.1038/s41467-021-23175-z 34006833
Elvsåshagen, T. et al. The genetic architecture of the human thalamus and its overlap with ten common brain disorders. Nat. Commun. 12, 2909 (2021).34006833 10.1038/s41467-021-23175-z
74. Hindley G The shared genetic basis of mood instability and psychiatric disorders: A cross-trait genome-wide association analysis Am. J. Med Genet. B Neuropsychiatr. Genet 2022 189 207 218 10.1002/ajmg.b.32907 35841185
Hindley, G. et al. The shared genetic basis of mood instability and psychiatric disorders: A cross-trait genome-wide association analysis. Am. J. Med Genet. B Neuropsychiatr. Genet 189, 207–218 (2022).35841185 10.1002/ajmg.b.32907
75. Jung K Leveraging genetic overlap between irritability and psychiatric disorders to identify genetic variants of major psychiatric disorders Exp. Mol. Med. 2023 55 1193 1202 10.1038/s12276-023-01005-0 37258574
Jung, K. et al. Leveraging genetic overlap between irritability and psychiatric disorders to identify genetic variants of major psychiatric disorders. Exp. Mol. Med. 55, 1193–1202 (2023).37258574 10.1038/s12276-023-01005-0
76. Karadag N Identification of novel genomic risk loci shared between common epilepsies and psychiatric disorders Brain 2023 146 3392 3403 10.1093/brain/awad038 36757824
Karadag, N. et al. Identification of novel genomic risk loci shared between common epilepsies and psychiatric disorders. Brain 146, 3392–3403 (2023).36757824 10.1093/brain/awad038
77. Li Z Li D Chen X Characterizing the polygenic overlaps of bipolar disorder subtypes with schizophrenia and major depressive disorder J. Affect. Disord. 2022 309 242 251 10.1016/j.jad.2022.04.097 35487438
Li, Z., Li, D. & Chen, X. Characterizing the polygenic overlaps of bipolar disorder subtypes with schizophrenia and major depressive disorder. J. Affect. Disord. 309, 242–251 (2022).35487438 10.1016/j.jad.2022.04.097
78. Parker N Psychiatric disorders and brain white matter exhibit genetic overlap implicating developmental and neural cell biology Mol. Psychiatry 2023 28 4924 4932 10.1038/s41380-023-02264-z 37759039
Parker, N. et al. Psychiatric disorders and brain white matter exhibit genetic overlap implicating developmental and neural cell biology. Mol. Psychiatry 28, 4924–4932 (2023).37759039 10.1038/s41380-023-02264-z
79. Torgersen K Shared genetic loci between depression and cardiometabolic traits PLoS Genet. 2022 18 e1010161 10.1371/journal.pgen.1010161 35560157
Torgersen, K. et al. Shared genetic loci between depression and cardiometabolic traits. PLoS Genet. 18, e1010161 (2022).35560157 10.1371/journal.pgen.1010161
80. Zheng H Identify novel, shared and disorder-specific genetic architecture of major depressive disorder, insomnia and chronic pain J. Psychiatr. Res. 2022 155 511 517 10.1016/j.jpsychires.2022.09.036 36191519
Zheng, H. et al. Identify novel, shared and disorder-specific genetic architecture of major depressive disorder, insomnia and chronic pain. J. Psychiatr. Res. 155, 511–517 (2022).36191519 10.1016/j.jpsychires.2022.09.036
81. Cheng W Shared genetic architecture between schizophrenia and subcortical brain volumes implicates early neurodevelopmental processes and brain development in childhood Mol. Psychiatry 2022 27 5167 5176 10.1038/s41380-022-01751-z 36100668
Cheng, W. et al. Shared genetic architecture between schizophrenia and subcortical brain volumes implicates early neurodevelopmental processes and brain development in childhood. Mol. Psychiatry 27, 5167–5176 (2022).36100668 10.1038/s41380-022-01751-z
82. Satizabal CL Genetic architecture of subcortical brain structures in 38,851 individuals Nat. Genet. 2019 51 1624 1636 10.1038/s41588-019-0511-y 31636452
Satizabal, C. L. et al. Genetic architecture of subcortical brain structures in 38,851 individuals. Nat. Genet. 51, 1624–1636 (2019).31636452 10.1038/s41588-019-0511-y
83. Bahrami S., et al. Unveiling the genetic landscape of Basal Ganglia: Implications for common brain disorders. medRxiv, 2023.2007. 2026.23293206 (2023).
84. Kircher M A general framework for estimating the relative pathogenicity of human genetic variants Nat. Genet. 2014 46 310 315 10.1038/ng.2892 24487276
Kircher, M. et al. A general framework for estimating the relative pathogenicity of human genetic variants. Nat. Genet. 46, 310–315 (2014).24487276 10.1038/ng.2892
85. Boyle AP Annotation of functional variation in personal genomes using RegulomeDB Genome Res. 2012 22 1790 1797 10.1101/gr.137323.112 22955989
Boyle, A. P. et al. Annotation of functional variation in personal genomes using RegulomeDB. Genome Res. 22, 1790–1797 (2012).22955989 10.1101/gr.137323.112
86. Roadmap EC Integrative analysis of 111 reference human epigenomes Nature 2015 518 317 330 10.1038/nature14248 25693563
Roadmap, E. C. et al. Integrative analysis of 111 reference human epigenomes. Nature 518, 317–330 (2015).25693563 10.1038/nature14248
87. Ernst J Kellis M ChromHMM: automating chromatin-state discovery and characterization Nat. Methods 2012 9 215 216 10.1038/nmeth.1906 22373907
Ernst, J. & Kellis, M. ChromHMM: automating chromatin-state discovery and characterization. Nat. Methods 9, 215–216 (2012).22373907 10.1038/nmeth.1906
88. Kolberg L g:Profiler-interoperable web service for functional enrichment analysis and gene identifier mapping (2023 update) Nucleic Acids Res. 2023 51 W207 w212 10.1093/nar/gkad347 37144459
Kolberg, L. et al. g:Profiler-interoperable web service for functional enrichment analysis and gene identifier mapping (2023 update). Nucleic Acids Res. 51, W207–w212 (2023).37144459 10.1093/nar/gkad347
89. Lee PH Genomic relationships, novel loci, and pleiotropic mechanisms across eight psychiatric disorders Cell 2019 179 1469 1482.e1411 10.1016/j.cell.2019.11.020 31835028
Lee, P. H. et al. Genomic relationships, novel loci, and pleiotropic mechanisms across eight psychiatric disorders. Cell 179, 1469–1482.e1411 (2019).31835028 10.1016/j.cell.2019.11.020
90. de la Fuente J Davies G Grotzinger AD Tucker-Drob EM Deary IJ A general dimension of genetic sharing across diverse cognitive traits inferred from molecular data Nat. Hum. Behav. 2021 5 49 58 10.1038/s41562-020-00936-2 32895543
de la Fuente, J., Davies, G., Grotzinger, A. D., Tucker-Drob, E. M. & Deary, I. J. A general dimension of genetic sharing across diverse cognitive traits inferred from molecular data. Nat. Hum. Behav. 5, 49–58 (2021).32895543 10.1038/s41562-020-00936-2
91. Liu M. Uncovering the shared genetic architecture between depression and subcortical volumes. Zenodo 10.5281/zenodo.12912958 (2024).
