Genome-wide study of half a million individuals with major depression identifies 697 independent associations, infers causal neuronal subtypes and biological targets for novel pharmacotherapies ================================================================================================================================================================================================= * Andrew M McIntosh * Cathryn M Lewis * Mark J Adams for the Psychiatric Genomics Consortium Major Depressive Disorder Working Group ## Abstract In a genome-wide association study (GWAS) of 685,808 individuals with major depression (MD) and 4,364,225 controls from 29 countries and across diverse and admixed ancestries, we identify 697 independent associations at 636 genetic loci, 293 of which are novel. Using fine-mapping and functional genomic datasets, we find 308 high-confidence gene associations and enrichment of postsynaptic density and receptor clustering. Leveraging new single-cell gene expression data, we conducted a causal neural cell type enrichment analysis that implicated excitatory and inhibitory midbrain and forebrain neurons, peptidergic neurons, and medium spiny neurons in MD. Critically, our findings are enriched for the targets of antidepressants and provide potential antidepressant repurposing opportunities (e.g., pregabalin and modafinil). Polygenic scores (PGS) from European ancestries explained up to 5.7% of the variance in liability to MD in European samples and PGS trained using either European or multi-ancestry data significantly predicted case control status across all four diverse ancestries. These findings represent a major advance in our understanding of MD across global populations. We provide evidence that MD GWAS reveals known and novel biological targets that may be used to target and develop pharmacotherapies addressing the considerable unmet need for effective treatment. ## Introduction Major depression (MD) is a leading cause of worldwide disability and affects approximately 15% of the population during their lifetime. The peak age of onset is in early adulthood and the disorder is typically recurrent or chronic in nature, with persisting disability despite pharmacological and psychological therapies. Twin and family-based studies provide evidence of a significant genetic contribution to its etiology, with a heritability of approximately 37% (Polderman et al., 2015). Since 2013, genome-wide association studies (GWAS) have provided major insights into the polygenic nature of MD, its genetic risk factors and underlying mechanisms (Ripke et al., 2013; CONVERGE Consortium et al., 2015; Hyde et al., 2016; Howard et al., 2018; Wray et al., 2018; Howard et al., 2019; Levey et al., 2021; Als et al., 2022; Meng et al., 2024). The largest study conducted to date involved a meta- analysis of the Million Veteran Program, 23andMe, UK Biobank, FinnGen, and iPSYCH including 371K cases (Als et al., 2022) reporting 243 independent MD risk loci. Nevertheless, the molecular, cellular, and neurobiological mechanisms of MD remain largely unidentified, limiting the development of disease models and mechanism-informed drug treatments (Zhu, 2020). In the current study, we report results from the Psychiatric Genomics Consortium (PGC) Major Depressive Disorder Working Group’s largest GWAS of MD to date. We used strategies designed for analysis of multi-ancestry and admixed populations to implement the largest, most inclusive study of MD genetics to date. These results substantially extend previous GWAS findings, provide evidence of causal genes, cell types and tissues involved in MD, and demonstrate out-of-sample prediction across diverse ancestry groups. ## Results ### **GWAS associations** We meta-analyzed (see **Methods** and **Figure S1, Tables S1, S2, supplementary study information**) GWAS summary statistics from 109 ancestrally diverse cohort datasets with 685,808 cases and 4,364,225 controls (with power equivalent to a case-control study of 1,000,101 cases and 1,000,101 controls, with 23% in diverse/non-European ancestries, Table 1). For cohorts with diverse ancestries, associations were assessed using tools that explicitly model population structure, admixture, and relatedness (GENESIS). For a subset of cohorts with ancestrally diverse samples, we compared the sample size using the commonly used strategy of assigning individuals into ancestry groups followed by logistic regression (N=24,859) with our joint approach (N=47,642) and found a 92% sample size increase, leading to a final sample size of 160,611 cases and 1,001,890 controls (Supplementary Material) and the discovery of an additional 66 genome wide significant loci. Using conditional and joint GCTA-COJO (Yang et al., 2012) analysis with threshold P ≤ 5 × 10-8 within 10 Mb windows, we found 697 significant independent SNPs in 635 genomic regions, about half (293/635; 46%) of which were novel MD associations (Figure 1, Table S3). Of these, 27 were identified due to the inclusion of cohorts with ancestrally diverse samples. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/06/24/2024.04.29.24306535/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2024/06/24/2024.04.29.24306535/F1) Figure 1: Manhattan Plot of GWAS meta-analysis of 685,808 MD cases and 4,364,225 controls ***Legend****: Manhattan plot displaying the significance of each SNP’s association with MD across the genome (vertical axis). Chromosomal position is shown on the horizontal axis. The horizontal line at 7.3 (-log10(5 x 10-8)) depicts the statistical significance threshold*. View this table: [Table 1:](http://medrxiv.org/content/early/2024/06/24/2024.04.29.24306535/T1) Table 1: Details of diverse ancestry studies included in the current GWAS Legend to table 1: Table shows the breakdown of the studies included in the current meta- analysis broken down by ancestry group and admixture/multiple ancestries (analysed using GENESIS software). The complete dataset included 109 individual datasets from 93 individual studies, some of which (e.g. UK Biobank) contributed more than one non-overlapping ancestry specific dataset. Neff/2 represents the equal case and control sample sizes of an equivalently powered balanced study. View this table: [Table 2:](http://medrxiv.org/content/early/2024/06/24/2024.04.29.24306535/T2) Table 2: Significant drug target enrichments *Legend* *Table shows the top 16 most significantly enriched drugs based on capture of their targets within the gene-based associations of the current MD GWAS analysis. One topical preparation is not shown. The test for drug enrichment is not directional and may indicate compounds that confer risk of MD or exacerbate depressive symptoms, as well as those that ameliorate risk or depressive symptoms. Q value is false discovery rate, Benjamini–Yekutieli corrected (competitive analysis)*. We carried out a fixed-effects meta-analysis for samples of European ancestries (525,197 cases and 3,362,335 controls), for which we had a large single linkage disequilibrium reference dataset, to carry out downstream analyses including heritability, gene prioritization, enrichment, genetic correlation and Mendelian Randomization analyses To examine the consequences of MD phenotyping on the meta-analyses, we implemented genomic structural equation modelling (SEM) with a common factor meta-analysis of the European- ancestry summary statistics in Genomic SEM (Grotzinger et al., 2019) (**Figure S2**). Cohorts were first meta-analyzed based on how the MD phenotype was determined: clinical/interview, electronic health record [EHR], questionnaire, or self-report of MD diagnosis. The proportion of total effective sample size contributed by each phenotype definition was 4% clinical/interview, 54% EHR, 14% questionnaire, and 27% self-report. The different phenotype definitions of MD had strong genetic correlations (LD score *rg* from 0.78 to 0.88). We fit a common factor model in Genomic SEM and treated the clinical/interview phenotype as the primary phenotype by fixing its factor loading to 1 and its residual variance to 0. This factor model was consistent with the data (X2 = 4.49, P = 0.213), meaning that we could not reject null hypothesis that a single factor capturing all the variance of the primary phenotype explained the intercorrelations among the other depression phenotypes. Most MD phenotypes had strong positive loadings on the common factor (clinical/interview = 1.0 [reference], EHR = 0.92±0.04, questionnaire = 0.95±0.04) although the loading for self- reported diagnosis was lower (self-report loading = 0.85±04). One locus showed significant SNP heterogeneity between phenotyping definitions (rs12124523 intronic variant in *NEGR1*, common factor association P = 8.4 × 10-14, Q heterogeneity P = 2.9 × 10-10, *I2* = 0.71) with a stronger association found in self-reported depression studies (Self-report odds ratio (OR) = 1.081, CI = 1.065–1.098, other cohorts OR = 1.008, CI = 0.999–1.018). We found no evidence of heterogeneity at 569/570 loci, supporting the use of multiple phenotypes in genetic association studies of MD. SNP-based heritability was estimated in European ancestries using SBayesS (Zeng et al., 2021) to be 8.4% (s.e. 0.07%) on the liability scale (assuming lifetime MD risk of 15%) similar to prior estimates (Howard et al., 2018; Wray et al., 2018). SBayesS provides genetic architecture estimates of polygenicity of 6% and selection parameter of -0.54. Compared to previously reported estimates for 155 traits, MD has a relatively higher polygenicity, but its associated variants are under weaker negative selection (Zeng et al., 2021). ### Gene prioritization and pathway enrichment analysis We used a range of methods and functional genomic datasets to gain insight into the variants, genes, and pathways that may underlie MD-associated loci: SNP-based fine- mapping of MD-associated loci, integration of expression and protein quantitative trait loci (eQTL and pQTL) data, transcriptome- and proteome-wide association studies (TWAS and PWAS) with summary data-based Mendelian randomization (SMR), colocalization (COLOC) and TWAS-based fine-mapping (of eQTLs and pQTLs, in FOCUS) analyses. We also mapped associated loci to genes using gene-based association analysis in fastBAT, chromatin interaction datasets (HiC) and applied a novel gene prioritization package PsyOPS (**see Methods**). We undertook functionally informed SNP-based fine-mapping analyses, using the European ancestry GWAS findings, targeted at all autosomal GWAS loci excluding the Major Histocompatibility Complex (MHC) region. Twenty-four variants showed strong putative evidence of causality (posterior inclusion probability, PIP > 0.95, **Table S4**). Credible causal set sizes comprising ≤10 variants (cumulative PIP > 0.95) were identified at 224 loci (**Figure S3**) and 234/564 autosomal loci could be mapped to one or more genes (**Table S5**). Both eQTL and pQTL data were used to infer up- or down-regulated gene expression or protein levels associated with MD. Stringent criteria were used to identify high confidence associations with MD (**Methods**). MD genetic associations were found to correlate and colocalize with cis-regulated expression of 75 genes (**Table S6**) and cis-regulated levels of 10 proteins (**Table S7**). No gene was identified as high confidence by both eQTL and pQTL analyses. In total, across SNP-based fine-mapping, eQTL and pQTL analyses, 308 high-confidence associations were identified, with 14 eQTL genes and 1 pQTL gene also identified as high confidence by SNP-based fine-mapping. For example, SNP-based fine-mapping found all SNPs in one 95% credible set were within the cytochrome P450 gene *CYP7B1,* which was also inferred to have decreased expression in the dorsolateral prefrontal cortex of individuals with MD (TWAS p-value = 2.92x10-15, COLOC PP4 = 0.939, FOCUS PIP = 1). Positional mapping approaches were used to identify additional genes that may be involved in MD etiology, including identification of the nearest gene to lead MD variants, aggregating genetic associations across gene regions using fastBAT (**Table S8, Figure S4**), and linking associated loci to genes through Hi-C chromatin interactions using H-MAGMA (**Table S9**). Furthermore, the gene prioritization method PsyOPS was used to score genes based on prior information on mutational constraint, brain expression, and involvement in neurodevelopmental disorders (**Table S10**). Of the 18,737 genes assessed using fastBAT, 1,568 were associated with MD (*P* < 2.67 x 10-6) with the dopamine receptor D2 (*DRD2*) gene providing the strongest association (P = 9.39 × 10-47). *DRD2* was also identified as associated with MD by H-MAGMA with astrocyte cells providing the strongest evidence from the four brain tissue profiles analyzed (P = 2.72 × 10-15). An additional 1,033 genes were also identified as associated with MD based on three-dimensional chromatin data using H- MAGMA. While PsyOPS prioritized a neighboring gene, *NCAM1* (PsyOPS score = 0.402), *DRD2* had an equivalent score (0.399). Other genes with high PsyOPS prioritization scores were *PTPRT, SLC12A5, RFX3, ELAVL2, HCN1, KIF5A, and SHANK3*. ### Synaptic gene set enrichment We used the high-confidence gene list from SNP-based fine-mapping, TWAS, and PWAS to identify enriched synapse functions using SynGO (Koopmans et al., 2019). The 43 genes from the high-confidence gene list with SynGO annotations were compared against a background of 18,035 brain-expressed genes. We replicated earlier findings from Howard et al (2019), showing enrichment of neuron differentiation processes and postsynaptic membrane components. The current GWAS provided greatly increased specificity, implicating the cytosol, active zone membrane, calcium levels, vesicle cycle and presynaptic endocytosis. Postsynaptically, there was enrichment of synaptic specialization, density, and receptor clustering (**Table S11A+B**). ### Tissue and cell type enrichment analysis We conducted tissue and cell type enrichment analysis using published expression datasets including bulk RNA-sequencing data from human tissues (Bryois et al., 2020) and single-cell RNA-sequencing data from the adult mouse central and peripheral nervous system (Zeisel et al., 2018). Across human tissues, we found clearer enrichment patterns of MD SNP- heritability in neural tissues using the current GWAS association findings than those obtained from the previous two PGC MDD group analyses **(Figure S5)**. In the adult mouse central and peripheral nervous system, we found significant enrichment of MD SNP- heritability in 10 out of 39 cell types with two different methods (MAGMA and partitioned LD Score, see **Figure 2 and Figure S6**). We confirmed all the cell types identified in the previous GWAS (Wray et al., 2018) including both excitatory and inhibitory neurons, but implicate multiple additional inhibitory neuron categories and peptidergic neurons. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/06/24/2024.04.29.24306535/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2024/06/24/2024.04.29.24306535/F2) Figure 2: Broad brain cell category enrichment analysis**. Legend to** Figure 2**: Cell type enrichment analysis**. *20 categories of brain cell types are listed (from a total of 39 broad cell type categories tested) along the vertical axis, horizontal bar size represents the significance of the enrichment measured using MAGMA gene set enrichment test or partitioned LDSC. Bars drawn in salmon color represent enrichments that were significant using both methods, green – MAGMA only, blue partitioned LDSC only and purple when neither method showed significant enrichment. 19 broad categories not displayed were also not significant using either method. Columns represent the results of each test using the summary statistics from MDD2013 (Ripke et al., 2013), MDD2018 (Wray et al., 2018) and the current study*. Analysis at a more refined level of cell types further emphasized the enrichment in excitatory and inhibitory neurons in multiple brain regions (**Figure S6**, **Table S12**). Associated cell types using both methods included midbrain (mouse atlas reference: MEGLU7, MEGLU8, MEGLU10 and MEGLU11), amygdala (TEGLU22), hippocampal (CA1, TEGLU21), thalamic (DEGLU4) and cortical (TEGLU1, TEGLU4, TEGLU8, TEGLU8, TEGLU11, TEGLU13 and TEGLU20) excitatory neurons. We also found additional evidence for the involvement of D1/D2 midbrain and striatal medium spiny neurons (MSN2 by both methods and MSN1,3-4 by MAGMA only). ### Drug target enrichment analysis Using Drug Targetor, we searched for Anatomical Therapeutic Chemical (ATC) drug classes whose targets were enriched in the signals from the GWAS analysis (Gaspar et al., 2019). Drug Targetor harnesses drug bioactivity data, to prioritize drugs and targets for a given phenotype. Replicating an earlier analysis, we found the gene targets of antidepressants (ATC class N06A) are significantly enriched (**Figure S7**) in our association findings. Other drug classes that were significantly enriched included antipsychotics (N05A), which includes some medicines with antidepressant effects (**Table S13B**). The gene targets of *specific drugs* were also enriched in genetic associations with MD, although the analysis does not infer whether the effects of these agents were more likely to be congruent or opposed to the effects of genetic risk. The identified drugs provide possible repurposing opportunities and examples included several anti-cancer therapies, and the agents pregabalin (used in the management of pain and anxiety) and modafinil which is used to treat daytime sleepiness caused by narcolepsy (**Tables 2 and S13A**). ### Within and cross-trait prediction #### Polygenic score prediction in European ancestry samples Using the case-control cohorts in the meta-analysis, we conducted a leave-one-cohort-out GWAS meta-analysis for 43 European ancestry cohorts that had provided individual level data. Polygenic scores (PGS) were generated in the left-out target European samples generated using SNP weights for the multi-ancestry and the European ancestry meta- analyses derived using SBayesR (Lloyd-Jones et al., 2019). Other PGS methods including the standard p-value clumping and thresholding gave similar results (**Table S14**). Across all European cohorts, the variance explained on the liability scale (was 5.8% (s.e. 0.2%) (**Figure S8, S9)**, with an AUC statistic of 0.625 (**Figure S10**). Adding functional annotations into the algorithm to generate SNP weights for PGS (SBayesR) increased prediction accuracy by 0.1% (i.e., of 5.9%). The was more than 1.4 times greater than that reported in the PGC MDD 2018 analysis (Wray et al., 2018; Ni et al., 2021) (**Figure 3a and S8**). The OR for being a case per standard deviation (SD) increase in PGS was 1.57. The OR for being a case in the tenth compared to the first decile of polygenic scores was 4.92 (95% CI 4.57-5.29) (**Figure 3**), and the OR for the top vs bottom centiles was 11.8 (95% CI 8.4-15.2) (**Figure S11**). The non-linear shape of these decile and centile plots is expected under a polygenic architecture (Baselmans et al., 2021). Heterogeneity in the out-of-sample prediction results could be partly explained by the recorded ascertainment type **(Figure 3b and S9),** which we classified as ‘clinical’ (16 cohorts; ascertained from in- or out-patient settings, or EHR) or ‘community’ ascertained (32 cohorts; interviews or questionnaires self- reporting on lifetime depression). The difference in mean PGS between clinical vs community cases was 0.131 (s.e 0.012, P< 2x10-16) control sample SD units. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/06/24/2024.04.29.24306535/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2024/06/24/2024.04.29.24306535/F3) Figure 3: MD polygenic score prediction into European ancestry studies ***Legend: (A)*** *Comparison of liability R2 by input summary statistics by availability (full dataset with 23andMe versus public dataset without 23andMe, using P-value clumping + thresholding at P <= 0.05 [P+CT]), PGS method (P+CT versus SBayesR), and discovery dataset (previous Howard 2019 versus current MDD2024 SBayesR). The R2 are estimated across 43 cohorts with individual level data when each was left out of the discovery sample in turn. For the Discovery panel the R2 are estimated from the 20 cohorts with individual level data contributed to the PGC subsequent to the Howard 2019 study. The* r2 *was calculated using a lifetime prevalence of 0.15. **(C)** Odds ratio by decile, with reference to decile 1, for clinical and community ascertained studies (SBayesR). Bars reflecting the 95%CI are based on estimates from the logistic regression*. ### Cross-ancestry prediction of MD We used data from 9 diverse ancestry studies to assess PGS transferability (**Table S15**) using PGS derived from the clumping and thresholding approach. The PGS were derived from the multi-ancestry and the European meta-analysis excluding 23andMe because of data access restrictions (Neffective = 739,180 and 576,327, respectively) **(Figure 4, Table S16**). In the diverse ancestry studies, the r2, by the PGS based on the European ancestry training data ranged from ∼0.6-4.5%. The r2 values for prediction into European ancestry (excluding 23andMe) were 3.9% (s.e. 0.2%) using PT=0.05. Values were lowest in studies with participants of African descent, and in the largest African ancestry study, the Million Veteran Program (MVP), the PGS was not associated with MD r2=0.0018). Results using the multi-ancestry summary statistics showed only minor and non-significant differences from European-only PGS GWAS trained scores in all ancestry groups. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/06/24/2024.04.29.24306535/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2024/06/24/2024.04.29.24306535/F4) Figure 4: Polygenic prediction of MD status from European and multi-ancestry GWAS into ancestrally diverse non-European studies ***Legend:*** *Details of cohorts found in Table S2. The was calculated using a prevalence of 0.15 with the P+CT method. The error bars are confidence intervals calculated using bootstrap. The training data did not include 23andMe, because of access limitations. AFR: African ancestry; AMR: Hispanic and Latin American ethnicities; EAS: Easy Asian ancestries; EUR: European ancestries; SAS: South Asian ancestries*. ### Genetic correlation and phenome-wide Mendelian randomization analysis We estimated LD score regression genetic correlations using 3,229 sets of summary statistics obtained from the OpenGWAS catalog (**Table S17**). Among phenotypes that were not direct measures of depression, the largest genetic correlation effect sizes with MD were with neuroticism (*rg* = 0.70, P = 2.02 × 10-162) and subjective well-being (*rg* = -0.63, P = 1.27 × 10-26). Compared to Howard et al (2019), novel findings were restricted to phenotypes with smaller genetic correlations (|*rg*| <0.16) likely reflecting the greater power provided by current GWAS summary statistics. In a phenome wide association study (PheWAS), we first identified associations between MD PGS, health related traits and potentially modifiable environmental factors **(Table S18 & Fig S12)**. We then tested for evidence of possible causal associations between MD and MD-PGC associated traits using bidirectional two-sample Mendelian Randomization (MR) in MR-Base (149 phenotypes; Table S19) and using summary statistics from European- ancestry participants in UK Biobank (256 phenotypes; Table S20), leaving UK Biobank out of the European-ancestry MD GWAS. Lower salt usage, faster walking pace, and higher educational attainment were associated with a reduced liability to MD. Greater body mass index (BMI), bio-impedance trunk adiposity measures, and waist and hip circumference were all associated with an increased risk of MD (absolute effect sizes ranged from 0.031 to 0.109, pFDR ranged from 0.006 to 4.33×10-15). A full list of these and other potentially causal associations is shown in **Figure S13** and **Table S19A-C**. MR-Base analyses (**Table S20A- C)**, supported a causal role for BMI and measures of adiposity in conferring liability to MD (absolute effect sizes from 0.021 to 0.089, pFDR from 0.01 to 2.24×10-5). We found evidence of potential causal consequences of MD for a number of health behaviors (e.g., increased alcohol and salt intake, increased TV use, absolute effect sizes ranged from 0.023 to 0.207, pFDR ranged from 0.01 to 8.49×10-16) and known disease risk factors (e.g., higher triglycerides, C reactive protein (CRP), and gamma glutamyl transferase levels, lower vitamin D levels, higher diastolic blood pressure and waist circumference, absolute effect sizes ranged from 0.018 to 0.128, pFDR ranged 0.046 to 9.63×10-22) **(Figure S13; Table S19D-F)**. Lower gray matter, brain structure volumes, and fluid intelligence score were also implicated as causal consequences of MD (absolute effect sizes from 0.048 to 0.102, pFDR from 0.044 to 3.01×10-16, see **Figure S13, Table S20D-F)**. ## Discussion The current study represents the largest and most-inclusive GWAS of MD to date, providing a yield of 697 independent SNP associations located within 636 independent genetic loci and evidence that neuronal differentiation and receptor clustering are involved in the aetiology of the disorder. 286 high-confidence gene associations were identified (summarized in Table S21) in European ancestries, including *CYP7B1*, a cytochrome P450 enzyme involved in neurosteroid synthesis. Our results confirm and extend previous findings showing the enrichment of expression signals in excitatory and inhibitory neurons. Importantly, the increased power provided by this genetic analysis provided additional evidence of involvement of amygdala and hippocampal excitatory neurons including granule cells, and medium spiny neurons. The amygdala and hippocampus have been previously implicated from a wide range of human imaging (Whalen et al., 2002; Schmaal et al., 2016) and animal studies of depression (Hall et al., 2001; Jentsch et al., 2002; Warner-Schmidt and Duman, 2006) and medium spiny neurons have also been previously implicated in animal studies of reward and linked to depressive behaviors (Lammel et al., 2014; Soares-Cunha et al., 2020). The enrichment of expression signals in granule cells is of particular interest given the renewal of this cell type throughout adult life in the dentate gyrus (Boldrini et al., 2018), its role in stress resilience (Holland, 2012) and the increased hippocampal granule cell expansion associated with antidepressant treatment (Boldrini et al., 2009). Together, these findings underline the mechanistic insights provided by expansion of GWAS to over half a million depressed individuals. Lack of ancestral and global diversity remain a significant concern for GWAS, with 86% of studies conducted in participants of European ancestry (Fatumo et al., 2022). Our study included data from 160,611 cases and 1,001,890 controls of non-European diverse ancestries. Unlike most other multi-ancestry GWAS, we used a joint analysis approach and did not exclude individuals with mixed ancestry or ancestry not represented in reference sets. This is becoming ever more important as the number of people with mixed ancestry are increasing in countries such as the USA and the UK (Livingston, 2022). Overall, the additional ancestrally diverse participants helped identify 27 novel genetic associations and enabled for the first time to demonstrate significant genetic risk prediction across diverse ancestry groups. Using polygenic scores, the proportion of variation in liability to MD explained in European ancestry case control studies also showed a considerable increase from an *R2* of 3.2% in our previous analyses, to 5.8% using SbayesR. We also show a significant MD prediction in diverse non-European and admixed ancestries. The SNP-*h2* in the current study of 8.4% implies that approximately 69% of the additive genetic variance for MD associated with common SNPs across studies can now be accounted for by polygenic scores. This study provides the first evidence of limited transferability of MD PGS to multiple diverse ancestries, and further emphasizes the importance of conducting future GWAS studies across different global populations, especially in Africa where transferability is poorest. Whilst we did not find evidence for improved prediction based on multi-ancestry rather than European-only PGS, this may be due to the small proportion of participants *within* each individual ancestry group (23% of individuals of Non-European ancestries were divided across 4 major ancestry and admixed groups) relative to the European ancestry group alone. Causal inference analyses using Mendelian Randomization provide the strongest genetic evidence to date that increased adiposity may increase the risk of MD. We also found that a less favorable profile of many behavioral and biochemical disease risk factors, reduced brain volumes and decreased fluid intelligence may be causal consequences of MD. These findings provide further evidence that behavioral and pharmacological interventions to reduce adiposity may reduce the risk of MD, and that reducing MD risk may have favorable consequences for disease risk, brain health and cognition. Genome wide association signals for depression also showed enrichment for the targets of antidepressants, suggesting that they may also help to reveal other effective treatment targets and more effective interventions. Pregabalin is an effective treatment for neuropathic pain, focal seizures and generalized anxiety disorder (GAD) (Generoso et al., 2017). Studies have shown that it is effective for the treatment of GAD-related depressive symptoms (Dold et al., 2022), but randomized controlled trial (RCT) evidence supporting its use in MD is weak. While pregabalin is sometimes prescribed to augment other effective treatments in medication-resistant MD (Dold et al., 2022), there is now justification to consider its efficacy in MD within the context of an RCT. Together these findings highlight the value of ancestrally diverse genetic studies to prioritize the study of pathophysiological processes in MD. The clearer enrichment of antidepressant targets provides confidence that genetic association findings will be relevant to the deployment or repurposing of existing treatments. Critically, these findings suggest genetic associations will point to new drug targets and more effective therapies that may reduce the considerable disability caused by depression. ## Supporting information Full Authorship [[supplements/306535_file03.pdf]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors after publication. Data from 23andMe are available upon publication following application to 23andMe. ## Conflicts of interest Cathryn Lewis is a member of the SAB for Myriad Neuroscience and has received consultancy fees from UCB. ## Code and data availability Summary statistics are available for download from [https://pgc.unc.edu/for-researchers/download-results/](https://pgc.unc.edu/for-researchers/download-results/) and individual data by application to the Psychiatric Genomics Consortium Data Access Committee [https://pgc.unc.edu/for-researchers/data-access-committee/](https://pgc.unc.edu/for-researchers/data-access-committee/). Available summary statistics including 23andMe data require an approved application to 23andMe here: [https://research.23andme.com/dataset-access/](https://research.23andme.com/dataset-access/). Project code is available from [https://github.com/psychiatric-genomics-consortium/mdd-wave3-meta](https://github.com/psychiatric-genomics-consortium/mdd-wave3-meta). ## Acknowledgements We would like to thank the participants and investigators from all studies and the research participants and employees of 23andMe for making this meta-analysis possible. Major funding for the PGC is from the US National Institutes of Health (MH124873, MH124871). Statistical analyses were carried out on the NL Genetic Cluster Computer ([http://www.geneticcluster.org/](http://www.geneticcluster.org/)) hosted by SURFsara. The iPSYCH team acknowledges funding from the Lundbeck Foundation (grants R102- A9118 and R155-2014-1724), the Stanley Medical Research Institute, the Novo Nordisk Foundation for supporting the Danish National Biobank resource, and the GenomeDK HPC facility. This research has been conducted using the UK Biobank Resource (application 4844) and data from dbGaP (accession phs000021, phs000196, phs000187) and including data from: the Molecular Genetics of Schizophrenia Collaboration (Pablo Gejman, Northwestern University), the NINDS CIDR:NGRC Parkinson’s Disease Study, and the SNP Association Analysis of Melanoma: Case–Control and Outcomes Investigation (supported by FNIH GAIN study, [CA093459](http://medrxiv.org/lookup/external-ref?link\_type=GEN&access\_num=CA093459&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom), [CA097007](http://medrxiv.org/lookup/external-ref?link\_type=GEN&access\_num=CA097007&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom), [ES011740](http://medrxiv.org/lookup/external-ref?link\_type=GEN&access_num=ES011740&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom), and [CA133996](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=CA133996&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom)). Individual study funding and other acknowledgements are provided in the supplementary study information. This paper represents independent research partly-funded by the NIHR Maudsley Biomedical Research Centre and Maudsley NHS Foundation Trust and King’s College London; the views expressed are those of the authors and not necessarily those of the NIHR or the Department of Health and Social Care. The current work was also supported by the Wellcome Trust (220857/Z/20/Z) and the European Union under the Horizon 2020 research and innovation programme (No 847776 and 948561). ## Footnotes * * Full list of authors in supplementary information * The authors have been added * Received April 29, 2024. * Revision received June 22, 2024. * Accepted June 24, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. 1.Als, T.D., Kurki, M., Grove, J., Voloudakis, G., Therrien, K., Tasanko, E., Nielsen, T.T., Naamanka, J., Veerapen, K., Levey, D., et al. (2022). Identification of 64 new risk loci for major depression, refinement of the genetic architecture and risk prediction of recurrence and comorbidities. medRxiv, 2022.2008.2024.22279149. doi:10.1101/2022.08.24.22279149. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMi4wOC4yNC4yMjI3OTE0OXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDYvMjQvMjAyNC4wNC4yOS4yNDMwNjUzNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 2. 2.Baselmans, B.M.L., Yengo, L., van Rheenen, W., and Wray, N.R. (2021). Risk in Relatives, Heritability, SNP-Based Heritability, and Genetic Correlations in Psychiatric Disorders: A Review. Biol Psychiatry 89, 11–19. doi:10.1016/j.biopsych.2020.05.034. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2020.05.034&link_type=DOI) 3. 3.Boldrini, M., Fulmore, C.A., Tartt, A.N., Simeon, L.R., Pavlova, I., Poposka, V., Rosoklija, G.B., Stankov, A., Arango, V., Dwork, A.J., et al. (2018). Human Hippocampal Neurogenesis Persists throughout Aging. Cell Stem Cell 22, 589–599 e585. doi:10.1016/j.stem.2018.03.015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.stem.2018.03.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29625071&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 4. 4.Boldrini, M., Underwood, M.D., Hen, R., Rosoklija, G.B., Dwork, A.J., John Mann, J., and Arango, V. (2009). Antidepressants increase neural progenitor cells in the human hippocampus. Neuropsychopharmacology 34, 2376–2389. doi:10.1038/npp.2009.75. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2009.75&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19606083&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269760800003&link_type=ISI) 5. 5.Bryois, J., Skene, N.G., Hansen, T.F., Kogelman, L.J.A., Watson, H.J., Liu, Z., Eating Disorders Working Group of the Psychiatric Genomics, C., International Headache Genetics, C., andMe Research, T., Brueggeman, L., et al. (2020). Genetic identification of cell types underlying brain complex traits yields insights into the etiology of Parkinson’s disease. Nat Genet 52, 482-493. doi:10.1038/s41588-020-0610-9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-020-0610-9&link_type=DOI) 6. 6.CONVERGE Consortium, CARDIoGRAM Consortium, and Gerad Consortium (2015). Sparse whole-genome sequencing identifies two loci for major depressive disorder. Nature 523, 588–591. doi:10.1038/nature14659. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature14659&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26176920&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 7. 7.Dold, M., Bartova, L., Fugger, G., Mitschek, M.M., Fabbri, C., Serretti, A., Mendlewicz, J., Souery, D., Zohar, J., Montgomery, S., and Kasper, S. (2022). Pregabalin augmentation of antidepressants in major depression - results from a European multicenter study. J Affect Disord 296, 485–492. doi:10.1016/j.jad.2021.09.063. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jad.2021.09.063&link_type=DOI) 8. 8.Fatumo, S., Chikowore, T., Choudhury, A., Ayub, M., Martin, A.R., and Kuchenbaecker, K. (2022). A roadmap to increase diversity in genomic studies. Nature medicine 28, 243–250. doi:10.1038/s41591-021-01672-4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-021-01672-4&link_type=DOI) 9. 9.Gaspar, H.A., Hubel, C., and Breen, G. (2019). Drug Targetor: a web interface to investigate the human druggome for over 500 phenotypes. Bioinformatics 35, 2515–2517. doi:10.1093/bioinformatics/bty982. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bty982&link_type=DOI) 10. 10.Generoso, M.B., Trevizol, A.P., Kasper, S., Cho, H.J., Cordeiro, Q., and Shiozawa, P. (2017). Pregabalin for generalized anxiety disorder: an updated systematic review and meta- analysis. Int Clin Psychopharmacol 32, 49–55. doi:10.1097/YIC.0000000000000147. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/YIC.0000000000000147&link_type=DOI) 11. 11.Grotzinger, A.D., Rhemtulla, M., de Vlaming, R., Ritchie, S.J., Mallard, T.T., Hill, W.D., Ip, H.F., Marioni, R.E., McIntosh, A.M., Deary, I.J., et al. (2019). Genomic structural equation modelling provides insights into the multivariate genetic architecture of complex traits. Nature Human Behaviour 3, 513–525. doi:10.1038/s41562-019-0566-x. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41562-019-0566-x&link_type=DOI) 12. 12.Hall, J., Thomas, K.L., and Everitt, B.J. (2001). Fear memory retrieval induces CREB phosphorylation and Fos expression within the amygdala. The European journal of neuroscience 13, 1453–1458. doi:10.1046/j.0953-816x.2001.01531.x. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.0953-816x.2001.01531.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11298807&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000168275000021&link_type=ISI) 13. 13.Holland, S. (2012). Furthering the sceptical case against virtue ethics in nursing ethics. Nurs Philos 13, 266–275. doi:10.1111/j.1466-769X.2012.00541.x. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1466-769X.2012.00541.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22950730&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 14. 14.Howard, D.M., Adams, M.J., Clarke, T.-K., Hafferty, J.D., Gibson, J., Shirali, M., Coleman, J.R.I., Hagenaars, S.P., Ward, J., Wigmore, E.M., et al. (2019). Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nature Neuroscience 22, 343–352. doi:10.1038/s41593-018-0326-7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41593-018-0326-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30718901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 15. 15.Howard, D.M., Adams, M.J., Shirali, M., Clarke, T.-K., Marioni, R.E., Davies, G., Coleman, J.R.I., Alloza, C., Shen, X., Barbu, M.C., et al. (2018). Genome-wide association study of depression phenotypes in UK Biobank identifies variants in excitatory synaptic pathways. Nat Commun. 16. 16.Hyde, C.L., Nagle, M.W., Tian, C., Chen, X., Paciga, S.A., Wendland, J.R., Tung, J.Y., Hinds, D.A., Perlis, R.H., and Winslow, A.R. (2016). Identification of 15 genetic loci associated with risk of major depression in individuals of European descent. Nat Genet 48, 1031–1036. doi:10.1038/ng.3623. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3623&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27479909&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 17. 17.Jentsch, J.D., Olausson, P., Nestler, E.J., and Taylor, J.R. (2002). Stimulation of protein kinase a activity in the rat amygdala enhances reward-related learning. Biol Psychiatry 52, 111–118. doi:10.1016/s0006-3223(02)01358-6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0006-3223(02)01358-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12114002&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176777200005&link_type=ISI) 18. 18.Koopmans, F., van Nierop, P., Andres-Alonso, M., Byrnes, A., Cijsouw, T., Coba, M.P., Cornelisse, L.N., Farrell, R.J., Goldschmidt, H.L., Howrigan, D.P., et al. (2019). SynGO: An Evidence-Based, Expert-Curated Knowledge Base for the Synapse. Neuron 103, 217–234 e214. doi:10.1016/j.neuron.2019.05.002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2019.05.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31171447&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 19. 19.Lammel, S., Lim, B.K., and Malenka, R.C. (2014). Reward and aversion in a heterogeneous midbrain dopamine system. Neuropharmacology 76 *Pt B*, 351-359. doi:10.1016/j.neuropharm.2013.03.019. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuropharm.2013.03.019&link_type=DOI) 20. 20.Levey, D.F., Stein, M.B., Wendt, F.R., Pathak, G.A., Zhou, H., Aslan, M., Quaden, R., Harrington, K.M., Nunez, Y.Z., Overstreet, C., et al. (2021). Bi-ancestral depression GWAS in the Million Veteran Program and meta-analysis in >1.2 million individuals highlight new therapeutic directions. Nat Neurosci 24, 954–963. doi:10.1038/s41593-021-00860-2. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41593-021-00860-2&link_type=DOI) 21. 21.Livingston, G. (2022). The rise of multiracial and multiethnic babies in the U.S. [https://www.pewresearch.org/short-reads/2017/06/06/the-rise-of-multiracial-and-multiethnic-babies-in-the-u-s/](https://www.pewresearch.org/short-reads/2017/06/06/the-rise-of-multiracial-and-multiethnic-babies-in-the-u-s/). 22. 22.Lloyd-Jones, L.R., Zeng, J., Sidorenko, J., Yengo, L., Moser, G., Kemper, K.E., Wang, H., Zheng, Z., Magi, R., Esko, T., et al. (2019). Improved polygenic prediction by Bayesian multiple regression on summary statistics. Nat Commun. Nat Commun. 23. 23.Meng, X., Navoly, G., Giannakopoulou, O., Levey, D.F., Koller, D., Pathak, G.A., Koen, N., Lin, K., Adams, M.J., Renteria, M.E., et al. (2024). Multi-ancestry genome-wide association study of major depression aids locus discovery, fine mapping, gene prioritization and causal inference. Nat Genet 56, 222–233. doi:10.1038/s41588-023-01596-4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-023-01596-4&link_type=DOI) 24. 24.Ni, G., Zeng, J., Revez, J.A., Wang, Y., Zheng, Z., Ge, T., Restuadi, R., Kiewa, J., Nyholt, D.R., Coleman, J.R.I., et al. (2021). A Comparison of Ten Polygenic Score Methods for Psychiatric Disorders Applied Across Multiple Cohorts. Biol Psychiatry 90, 611–620. doi:10.1016/j.biopsych.2021.04.018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2021.04.018&link_type=DOI) 25. 25.Polderman, T.J., Benyamin, B., de Leeuw, C.A., Sullivan, P.F., van Bochoven, A., Visscher, P.M., and Posthuma, D. (2015). Meta-analysis of the heritability of human traits based on fifty years of twin studies. Nat Genet 47, 702–709. doi:10.1038/ng.3285. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3285&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25985137&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 26. 26.Prive, F., Aschard, H., Carmi, S., Folkersen, L., Hoggart, C., O’Reilly, P.F., and Vilhjalmsson, B.J. (2022). Portability of 245 polygenic scores when derived from the UK Biobank and applied to 9 ancestry groups from the same cohort. Am J Hum Genet 109, 373. doi:10.1016/j.ajhg.2022.01.007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2022.01.007&link_type=DOI) 27. 27.Ripke, S., Wray, N.R., Lewis, C.M., Hamilton, S.P., Weissman, M.M., Breen, G., Byrne, E.M., Blackwood, D.H., Boomsma, D.I., Cichon, S., et al. (2013). A mega-analysis of genome-wide association studies for major depressive disorder. Molecular psychiatry 18, 497–511. doi:10.1038/mp.2012.21. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2012.21&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=AMBIGUOUS (2&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 28. 28.Schmaal, L., Veltman, D.J., van Erp, T.G., Samann, P.G., Frodl, T., Jahanshad, N., Loehrer, E., Tiemeier, H., Hofman, A., Niessen, W.J., et al. (2016). Subcortical brain alterations in major depressive disorder: findings from the ENIGMA Major Depressive Disorder working group. Mol Psychiatry 21, 806–812. doi:10.1038/mp.2015.69. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2015.69&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26122586&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 29. 29.Soares-Cunha, C., de Vasconcelos, N.A.P., Coimbra, B., Domingues, A.V., Silva, J.M., Loureiro-Campos, E., Gaspar, R., Sotiropoulos, I., Sousa, N., and Rodrigues, A.J. (2020). Nucleus accumbens medium spiny neurons subtypes signal both reward and aversion. Mol Psychiatry 25, 3241–3255. doi:10.1038/s41380-019-0484-3. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41380-019-0484-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31462765&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 30. 30.Warner-Schmidt, J.L., and Duman, R.S. (2006). Hippocampal neurogenesis: opposing effects of stress and antidepressant treatment. Hippocampus 16, 239–249. doi:10.1002/hipo.20156. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hipo.20156&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16425236&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000235910900006&link_type=ISI) 31. 31.Whalen, P.J., Shin, L.M., Somerville, L.H., McLean, A.A., and Kim, H. (2002). Functional neuroimaging studies of the amygdala in depression. Semin Clin Neuropsychiatry 7, 234–242. doi:10.1053/scnp.2002.35219. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/scnp.2002.35219&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12382206&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 32. 32.Wray, N.R., Ripke, S., Mattheisen, M., Trzaskowski, M., Byrne, E.M., Abdellaoui, A., Adams, M.J., Agerbo, E., Air, T.M., Andlauer, T.F., et al. (2018). Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nature genetics 50, 668–681. doi:10.1038/s41588-018-0090-3. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0090-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29700475&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F24%2F2024.04.29.24306535.atom) 33. 33.Yang, J., Ferreira, T., Morris, A.P., Medland, S.E., Madden, P.A.F., Heath, A.C., Martin, N.G., Montgomery, G.W., Weedon, M.N., Loos, R.J., et al. (2012). Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nature genetics. Nat Genet. 34. 34.Zeisel, A., Hochgerner, H., Lönnerberg, P., Johnsson, A., Memic, F., van der Zwan, J., Häring, M., Braun, E., Borm, L.E., La Manno, G., et al. (2018). Molecular Architecture of the Mouse Nervous System. Cell. Cell. 35. 35.Zeng, J., Xue, A., Jiang, L., Lloyd-Jones, L.R., Wu, Y., Wang, H., Zheng, Z., Yengo, L., Kemper, K.E., Goddard, M.E., et al. (2021). Widespread signatures of natural selection across human complex traits and functional genomic categories. Nat Commun. Nat Commun. 36. 36.Zhu, T. (2020). Challenges of Psychiatry Drug Development and the Role of Human Pharmacology Models in Early Development-A Drug Developer’s Perspective. Front Psychiatry 11, 562660. doi:10.3389/fpsyt.2020.562660. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fpsyt.2020.562660&link_type=DOI)