White matter tract integrity is reduced in major depression and in individuals with genetic liability for depression ==================================================================================================================== * David Nothdurfter * Philippe Jawinski * Sebastian Markett ## 0. Abstract Major depression has been linked to alterations of the brain’s white matter architecture. Major depression has also a considerable genetic component. It is, however, unclear if genetic liability for depression relates directly to the suggested alterations at the brain level. We re-evaluated white matter fiber tract alterations in major depression and aimed to relate polygenic risk for depression in unaffected individuals to these alterations. We conducted a population-based retrospective cohort study with N=17183 participants (n=287 currently depressed, n=5536 formerly depressed, n= 13360 healthy) aged 44 to 82 from the UK Biobank imaging cohort, collected between 2015 and 2020. Depression status was assessed based on questionnaire data. We derived 27 major white matter tracts from diffusion MRI and indexed their integrity according to mean fractional anisotropy. Genetic liability for depression was quantified by polygenic scores based on three recent case-control GWAS. White matter integrity was globally reduced in depression. The reduction was most pronounced on thalamic and intracortical fiber tracts, where the reduction was observed irrespective of current symptom status. Non-depressed individuals showed already a reduction in relevant fiber tracts with increasing polygenic risk, particularly in thalamic radiations. Subsequent in-silico simulations and non-parametric tests confirmed that the polygenic association was stronger than expected given the low-level statistical dependencies between target and discovery samples. White matter integrity on thalamic and association tracts is reduced in depression irrespective of current disease status and relates to the underlying genetic risk factors in unaffected individuals. The observed statistical links between genes, brain, and disorder adhere to the definition of an intermediate phenotype at the brain level. ## 1. Introduction Depression is among the most prevalent psychiatric disorders in Europe and one of the leading causes of disability around the globe1,2. Recent genome-wide association studies (GWAS) have confirmed that depression is considerably heritable and highly polygenic with many individual genetic variants of very small effect sizes each3–6. At the time of writing, about 3.2% of the risk towards major depression can be explained by polygenic scores, i.e., weighted sum scores that aggregate information of multiple trait-associated alleles from GWAS4. Little is known, however, about how this genetic risk might translate into depression. In the present study, we aim at identifying brain phenotypes that may mediate the effects of genetic variation on disease risk. Neuroimaging studies suggest alterations of the integrity of white matter pathways in people suffering from depression7. While the findings of lower global white matter integrity in people suffering from depression are robust8, it remains open which specific brain pathways are altered. Previous research has yielded mixed results, also due to different MRI protocols or limited power in studies with low sample size9. Meta-analyses have found most consistently alterations in the genu of the corpus callosum but have produced mixed results regarding other fiber tracts8–11. Even larger sample size may be needed to overcome limitations and to resolve inconsistencies regarding implicated fiber tracts. Another reason for inconsistent findings might be heterogeneity in case definition when not distinguishing between currently depressed individuals and individuals with a life-time history of depression. It has been suggested, for instance, that individuals suffering from acute depression exhibit more pronounced alterations of connectivity patterns9. The goals of the present study are twofold. First, we seek to re-evaluate alterations in white matter integrity. By distinguishing between life-time as well as currently depressed individuals in a large imaging sample (UK Biobank with 19,183 participants) we seek to establish both global and local connectivity markers of depression. Second, we seek to link depression-related white matter alterations to the genetic risk for depression in the healthy population. Studies on the relationship between polygenic risk for depression and brain connectivity have yielded divergent results thus far12,13 but were also not based on the most recent GWAS findings on depression. Establishing neural markers for depression (goal 1) and linking these markers to polygenic liability (goal 2) are major steps towards a deeper understanding of the etiology of depression14. ## 2. Methods ### 2.1. Data resource All participants were drawn from the UK Biobank (UK Biobank, 2020a). Access to the UK Biobank was granted under project 42032 (PI: Sebastian Markett). We focused on a subset of *N* = 19,183 participants with neuroimaging (T1- and diffusion-weighted) and with complete genotyping data (50.27% men; mean age *M* = 64.08 years (*SD* = 7.52, range: 45 – 82). ### 2.2 Exclusion criteria We excluded participants with poor quality in the structural T1-weighted image, participants with sex-aneuploidy and whose reported sex did not match their genetic sex, and participants with a history of neurological disorders such as dementia or psychiatric disorders other than depression, social phobia or anxiety, panic attacks, or generalized anxiety disorder. Details on the exclusion procedure are described in the supplementary methods. ### 2.3 Depression status Depression status was assessed based on participants’ responses to a touchscreen questionnaire presented during the imaging visit15. We assigned participants into three groups: *depressed without current symptoms, currently depressed* and *healthy controls*. Details on the assignment to the groups are described in the supplementary methods. In total, depression status was obtained for 19,183 subjects (13360 healthy, 5536 formerly depressed and 287 currently depressed individuals). For descriptive statistics see supplementary table T1. ### 2.4 Group matching Because measures of white matter integrity are influenced by sex, age and image quality, we performed an additional analysis in a subset of the sample in which the three groups were matched for these variables (via propensity score matching; see supplementary methods). The matched sub-sample included 8653 participants (see figure 1 and supplementary tables T1-T4). ![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/11/14/2022.11.09.22282111/F1.medium.gif) [Figure 1](http://medrxiv.org/content/early/2022/11/14/2022.11.09.22282111/F1) Figure 1 Overview of the sub-samples included for the different analyses. Each depicted person-icon represents 1,000 participants. Smaller icons represent proportionally lower numbers of participants (unmatched sub-samples: depressed with current symptoms 287 participants, depressed without current symptoms 5536 participants, no depression 13360 participants, unclear depression status 15876 participants; matched sub-samples: depressed with current symptoms 286 participants, depressed without current symptoms 4185 participants, no depression 4182 participants; no depression (remaining unmatched) 8092 participants). ### 2.5 MRI analyses MRI data were acquired on Siemens Skyra 3T scanners at three different imaging facilities across the UK. We used the preprocessed imaging-derived phenotypes for DTI data, generated by an imaging-processing pipeline run by the UK Biobank16. In brief, tractography was perfomed in FSL’s probtrackx and 27 major tracts were automatically reconstructed from the data and characterized by mean fractional anisotropy (FA). More details are given in the supplementary methods. Because fiber tracts share a substantial amount of variance17, we computed summary FA measures (via principal component analyses) for groups of fiber tracts (all fiber tracts, association/ commissural fibers, thalamic radiations and projection fibers), following suggestions from previous work18 (see supplementary methods for more details). Multiple regression analyses were then performed to investigate differences in fractional anisotropy between the three groups, individually for the three comparisons (currently vs. formerly depressed, currently depressed vs. healthy, formerly depressed vs. healthy). The four PCs (summary FA measures) were each set as the dependent variable, while depression status served as the independent variable. Also, multiple regressions were performed with FA of individual fiber tracts as dependent variable (applying the Benjamini-Hochberg procedure with a false discovery rate of 5%19). ### 2.6. Polygenic risk score analyses Polygenic risk scores were based on three genome-wide association studies (GWAS)4–6. Details on the calculation of the polygenic risk sores are described in the supplementary methods. To investigate if higher polygenic risk for depression is associated with lower white matter integrity on depression-related fiber tracts, we analyzed the remaining healthy participants that were not used for the previous analyses with the matched subgroups. To ensure full equivalence between this analysis and the group comparison, we used PCs from the group comparison and projected the unmatched healthy individuals’ data onto these components. The resulting factor scores as well as FA values of individual fiber tracts were then correlated with participants’ polygenic risk scores for depression. Multiple regressions were performed with factor scores as dependent variable and each polygenic risk score as independent variables. In addition to the above-mentioned covariates, genotyping array and 20 genetic principal components were included. As the three GWAS discovery samples showed varying sample overlap (6.19, 44.74 and 59.58%) with the full UKB cohort from which we selected a subsample for the present study, we performed additional analyses to ensure that the overlap does not confound our conclusions. We benchmarked our results in comparison to a null model that was based on a series of simulations in which we explored how random variables with varying intercorrelations with depression status are related to the polygenic risk scores. The rationale behind these analyses was to quantify the impact of sample overlap between GWAS discovery and target samples under different boundary conditions. As part of this benchmarking, each empirical PGS association was compared to the distribution of PGS associations with 10,000 random synthetic variables which were similarly intercorrelated with depression status as white matter integrity. We refer to the supplementary material for more details. ### 2.7. Genetic correlations We conducted GWAS using *PLINK2* for the FA phenotypes that differed between groups. All participants with available genetic and imaging data were included for this step. We then performed linkage disequilibrium score regression to calculate SNP-based heritabilities20. Finally, genetic correlations were performed between the depression phenotypes and the FA phenotypes. ### 2.8. Covariates and Outlier exclusion Highest degree of education, intake of antidepressants, sex, age, age2, image quality and the location of the assessment center were entered as covariates (see supplementary methods for more details).For all analyses, we excluded participants with outliers in the PC or FA value of interest by applying the Tukey criterion. On average, 1.27% (*SD=*0.76) of the participants were excluded in the different analyses due to outliers. Figure 1 provides an overview of the different sub-samples of the UK Biobank used for the various analyses. ### 2.9 Data and Code availability statement All analyses were performed in *R* version 3.6.3. All data were taken from the UK Biobank. Data access is limited to registered researchers via [https://www.ukbiobank.ac.uk/](https://www.ukbiobank.ac.uk/). We will publish analysis code for the null model for PGS testing as well as code to reproduce the figures on the open science framework upon final acceptance of the manuscript. ## 3. Results ### 3.1. Differences in global FA components We first investigated group differences in global FA as quantified by the first principal component (PC) across all fiber tracts and across all fiber tract groups. Supplementary tables T5 and T6 depict loadings of individual fiber tracts on these PCs. Figure 2 shows the comparison of the four fiber groups’ PCs for the subgroup of matched participants. Multiple regression analyses revealed that depressed participants without current symptoms showed significantly lower values on the PC derived from all fiber tracts (*b* = −0.184, *p* = 0.008), association/commissural fibers and thalamic radiations than healthy participants (*b* = −0.219, *p* = .018; *b* = −0.132, *p* < .001). Between healthy participants and depressed participants with current symptoms, the PCs derived from association/commissural fibers and thalamic radiations differed significantly (*b* = −0.38, *p* = .019; *b*= −0.252, *p* = .023). No differences in any of the four PCs could be observed between the two depression groups. The results were similar to the analyses including the full sample, except for a lack of significant difference in thalamic radiations between currently depressed and healthy participants (see tables T7 to T9). ![Figure 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/11/14/2022.11.09.22282111/F2.medium.gif) [Figure 2](http://medrxiv.org/content/early/2022/11/14/2022.11.09.22282111/F2) Figure 2 Overall integrity of all white matter fiber tracts (A) was reduced in individuals with a history of depression but no current symptoms (D-) relative to healthy participants (H). For association fiber tracts (B) and thalamic radiations (C), the reduction was additionally visible for the smaller group of individuals with current symptoms (D+). No such difference was observed for projection fibers (D). The bar graphs in the middle column give the factor scores of the fiber tracts in the matched subgroups (caliper 0.05). The red lines mark significant group differences. The heatmaps on the right list association strengths (standardized beta coefficients) between fiber tract integrity and polygenic scores based on the GWAS by Howard et al. (H), Wray et al. (W), and Baselmans et al. (B) across different thresholds. Significant coefficients are highlighted according to the corresponding t-statistic. ### 3.2. Differences in individual fiber tracts Eleven individual fiber tracts showed decreased FA in depressed participants without current symptoms in at least the full sample or the sub-sample including matched participants, compared to healthy participants (figure 2). Between healthy participants and depressed participants with current symptoms significant differences in FA were found for five fiber tracts (figure 2). None of the fiber tracts differed between the two depression groups. ### 3.3. Polygenic risk scores Figure 2 depicts the associations between PGS and fiber integrity in healthy participants of the fiber tract groups that were linked to depression in the group comparison (all fiber tracts, commissural/association fibers, thalamic radiations). With increasing genetic liability for depression (higher PGS scores), healthy individuals showed reduced integrity of commissural/association fiber tracts and of the thalamic radiations, suggesting a subclinical alteration of white matter integrity. Of note, this relationship was observed for all three operationalizations of polygenic risk and extended across several PGS thresholds. Figure 3 depicts the association of PGS and the integrity of the individual fiber tracts for which we observed group differences between healthy participants and one of the two depression groups. Every fiber tract except the left inferior longitudinal fasciculus showed significant association with at least one of the three PGS, most consistently for the anterior thalamic radiations and the forceps minor. ![Figure 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/11/14/2022.11.09.22282111/F3.medium.gif) [Figure 3](http://medrxiv.org/content/early/2022/11/14/2022.11.09.22282111/F3) Figure 3 Eleven fiber tracts (top panel) showed group differences individually (middle panel). We present standardized beta coefficients for group differences between healthy individuals and individuals with a history of depression (no current symptoms) and individuals with current symptoms, based on matched groups (m) and all participants (a). Only FDR-significant coefficients are given and color-coded by the corresponding t-value. Polygenic score associations are depicted in the bottom panel (different discovery GWAS in columns, different thresholds in rows). Only significant associations are highlighted and color coded by the corresponding t-value. ### 3.4 Genetic correlations We sought to corroborate previous heritability estimates of depression through LD score regression: Depending on the summary statistics from the original GWAS, heritability estimates were 9.6%6, 6.3%4, and 3.1%5. SNP-based heritability for the global FA components was considerable, ranging between 29-30% (supplementary table T21). Genetic correlations between depression and global FA components were all in the expected directions but not statistically significant (supplementary table T21). SNP-based heritability for individual fiber tracts ranged between 19 and 33%. Genetic correlations with depression were all in the expected direction but did not survive the correction for multiple comparisons (see supplementary table T22). ### 3.5 Validation analysis To investigate the influence of partly overlapping participants between the present and the GWAS discovery samples, we simulated synthetic random data for which we systematically manipulated the degree of intercorrelation with actual depression status. We then ran partial correlation analyses identical to our actual analysis of white matter fiber statistics (see methods and supplementary methods). As expected, we observed a spurious increase in polygenic association with increasing intercorrelations between depression status and random data when including patients and controls (Figure 4A). When limiting the analysis to healthy controls, however, the partial correlations where close to zero (< |1e-15|), irrespective of the degree of intercorrelation between the synthetic data and depression status (figure 4B). While the simulation results suggest that partially overlapping discovery and target samples do not affect the validity of PGS associations in a scenario where GWAS summary statistics are obtained by case-control comparisons and the target sample includes only controls, we still decided to complement the simulations with a further validation step: For each fiber group, we established a null model of 10,000 random variables whose intercorrelation with depression status matched the empirical correlations between fiber integrity and depression status. If the here reported PGS associations with fiber integrity were simply a result of the partially overlapping samples, we would except similar PGS associations with the random data. This, however, was not the case: For association/commissural tracts and for the thalamic radiations, we confirmed significant reductions in fiber integrity in healthy participants with higher genetic liability for depression (figure 4C). Only for the global measure of white matter integrity of all fiber tracts, the comparison to the null model did not confirm the descriptive PGS association, suggesting that the observed correlation is not stronger than the inherent interdependencies resulting from sample overlap. ![Figure 4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/11/14/2022.11.09.22282111/F4.medium.gif) [Figure 4](http://medrxiv.org/content/early/2022/11/14/2022.11.09.22282111/F4) Figure 4 Results from the validation analyses: PGS associations with synthetic data are only substantial when the target sample includes cases and controls, and increase as a function of similarity between the variable of interest and disease status (A). When only healthy participants are included in the target sample, the PGS associations with random data is practically zero (B). The different lines reflect PGS on different thresholds, with more yellow hues reflecting less liberal thresholds. The highlighted empirical PGS associations in panel C are significantly larger (i.e. more negative) than PGS associations with comparable random data (p<.05, non-parametric test with 10,000 random variables). ## 4. Discussion The present study aimed at establishing structural connectivity markers for major depression, and linking depression-related white matter alterations to the genetic risk for depression. We found lower white matter fiber tract integrity in depression using data from the UK Biobank, extending previous work with an earlier data release18. Our results support previous findings on reduced integrity of interhemispheric fiber tracts, particularly the forceps minor8–11,21, fiber tracts belonging to the Papez circuit which plays a role in negative affect22, and the medial forebrain bundle which includes the cingulum bundle and the uncinate fasiculus23 (Bracht et al., 2015). We further observed that similar reductions in white matter integrity present already subclinically in healthy individuals with higher genetic liability for depression. Mental disorders are considerably heritable24 and contemporary GWAS results can be used to quantify an individual’s genetic liability for a disorder in the form of polygenic scores25. Linking PGS to brain structure and function bears the potential to elucidate the neural mechanisms by which genetic risk translates into the disorder, to assess interactions with environmental factors, and to characterize sub-clinical expression in unaffected populations. So far, PGS-brain relationships have mostly been investigated in the context of schizophrenia and autism1 but not all attempts have been successful31,32. PGS for depression and related phenotypes have also been linked to brain structure and function33–36, but not all studies have distinguished between cases and controls, or have linked the PGS-brain associations to case/control differences. Particularly the latter, however, is crucial to establish links between genes, brain, and behavior, and to establish intermediate phenotypes for the disorder. To the best of our knowledge, our current results are the first to confirm neurostructural differences in individuals with depression and link the same neurostructural variables to genetic liability for depression in unaffected individuals. Our results thus provide major implications for structural connectivity as an intermediate phenotype for depression. To qualify as an intermediate phenotype (or endophenotype) of depression, an imaging-derived phenotype needs to be related to the disorder, needs to be heritable, and needs to be relatable to the genetic underpinnings of the disorder in unaffected individuals14. We confirmed these criteria for white matter integrity: Table 1 summarizes our findings regarding these criteria for different groups of fiber tracts. Using LD-score regression, we confirmed a SNP-based heritability for white matter tract integrity ranging from 19 and 33%. Using polygenic scores, we confirmed that genetic liability for depression is not only associated with the disorder itself, but also with white matter tract integrity, which is a key postulate for an endophenotype37,38. Even though total effect sizes were low, the significant reductions in white matter integrity with increasing genetic risk load in healthy individuals indicate that white matter integrity of implicated fiber tracts fulfill the endophenotype criteria. The small effect sizes are a likely consequence of the high polygenicity of the genetic underpinnings of depression and the limited prediction accuracy of polygenic scores to date.39 View this table: [Table 1](http://medrxiv.org/content/early/2022/11/14/2022.11.09.22282111/T1) Table 1 Summary of the main results with respect to the endophenotype concept. White matter integrity of association and commissural fiber tracts is heritable, reduced in depression independent of current symptoms, and subclinically reduced in healthy individuals who are genetically at risk. A well-known pitfall in PGS analyses is the spurious increase in association strength when the GWAS discovery sample and the target sample overlap40. Even partial sample overlap can dramatically inflate observed associations. Partially overlapping samples, however, should be less of a concern when the GWAS summary statistics reflect a group comparison (such as individuals with depression vs. healthy controls) and the target variable of interest is only assessed in individuals belonging to either group (such as healthy participants only). We corroborated this conjecture through a series of simulations. And we established through a simple non-parametric test that the observed relationship between polygenic liability for depression and white matter integrity of intracortical and thalamic fiber tracts is substantial, even in the face of overlapping samples. Given the large sample size requirements of discovery GWAS39, sample overlap will remain a concern, at least until the predictive accuracy of PGS will improve to the point that meaningful associations can already be observed in small target samples41. The here presented non-parametric test does not depend on prior knowledge of sample overlap, does not require access to individual-level data from the discovery cohort, and is easy to implement, particularly in psychiatric research contexts where summary statistics often result from case/control comparisons. Contrary to the polygenic score analyses, genetic correlations did not reveal significant associations between FA phenotypes and depression. This may be explained by the lower statistical power of summary-based (genetic correlations via LD score regression) relative to individual-level association analysis (polygenic score analysis). Findings from genetic correlations may indicate that the total genetic overlap between FA phenotypes and depression is modest, in line with previous findings linking psychiatric diseases to brain structure42. The lack of significant genetic correlations between FA and depression needs further investigation. Another general limitation to the present study is the group definition which was not based on structured clinical interviews but on self-report data. Even though this assessment led to similar prevalence rates of depression as other studies15, it is still possible that depression in the present study may have been over-diagnosed. In summary, we report structural connectivity differences between healthy and depressed individuals. These alterations are present in currently depressed as well as formerly depressed individuals, are heritable, and seem to be related to higher genetic liability for depression in healthy individuals. Our results thus provide indications for an involvement of structural connectivity in the genetic etiology of depression. Identifying molecular substrates of brain connectivity may thus reflect a promising strategy to inform individualized treatment43–45 and improved diagnostics in the near future46,47 ## Supporting information Supplementary Materials [[supplements/282111_file02.pdf]](pending:yes) ## Data Availability All data produced are available online at ukbiobank.ac.uk. * Received November 9, 2022. * Revision received November 9, 2022. * Accepted November 14, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## 5. References 1. 1.1. Samet J Mathers CD, Loncar D. Projections of Global Mortality and Burden of Disease from 2002 to 2030. Samet J, ed. PLoS Med. 2006;3(11):e442. doi:10.1371/journal.pmed.0030442 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.0030442&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17132052&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 2. 2.Wittchen HU, Jacobi F, Rehm J, et al. The size and burden of mental disorders and other disorders of the brain in Europe 2010. Eur Neuropsychopharmacol J Eur Coll Neuropsychopharmacol. 2011;21(9):655–679. doi:10.1016/j.euroneuro.2011.07.018 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.euroneuro.2011.07.018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21896369&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 3. 3.Sullivan PF, Neale MC, Kendler KS. Genetic Epidemiology of Major Depression: Review and Meta-Analysis. Am J Psychiatry. 2000;157(10):1552–1562. doi:10.1176/appi.ajp.157.10.1552 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.157.10.1552&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11007705&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000089633900002&link_type=ISI) 4. 4.Howard DM, Adams MJ, Clarke TK, et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nat Neurosci. 2019;22(3):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%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 5. 5.Baselmans BML, Jansen R, Ip HF, et al. Multivariate genome-wide analyses of the well-being spectrum. Nat Genet. 2019;51(3):445. doi:10.1038/s41588-018-0320-8 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0320-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30643256&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 6. 6.Wray NR, Ripke S, Mattheisen M, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genet. 2018;50(5):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%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 7. 7.Zhuo C, Li G, Lin X, et al. The rise and fall of MRI studies in major depressive disorder. Transl Psychiatry. 2019;9(1):335. doi:10.1038/s41398-019-0680-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41398-019-0680-6&link_type=DOI) 8. 8.Liao Y, Huang X, Wu Q, et al. Is depression a disconnection syndrome? Meta-analysis of diffusion tensor imaging studies in patients with MDD. J Psychiatry Neurosci. 2013;38(1):49–56. doi:10.1503/jpn.110180 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianBuIjtzOjU6InJlc2lkIjtzOjc6IjM4LzEvNDkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8xMS8xNC8yMDIyLjExLjA5LjIyMjgyMTExLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 9. 9.van Velzen LS, Kelly S, Isaev D, et al. White matter disturbances in major depressive disorder: a coordinated analysis across 20 international cohorts in the ENIGMA MDD working group. Mol Psychiatry. 2020;25(7):1511–1525. doi:10.1038/s41380-019-0477-2 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41380-019-0477-2&link_type=DOI) 10. 10.Wise T, Radua J, Nortje G, Cleare AJ, Young AH, Arnone D. Voxel-Based Meta-Analytical Evidence of Structural Disconnectivity in Major Depression and Bipolar Disorder. Biol Psychiatry. 2016;79(4):293–302. doi:10.1016/j.biopsych.2015.03.004 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2015.03.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25891219&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 11. 11.Chen G, Hu X, Li L, et al. Disorganization of white matter architecture in major depressive disorder: a meta-analysis of diffusion tensor imaging with tract-based spatial statistics. Sci Rep. 2016;6(1):21825. doi:10.1038/srep21825 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/srep21825&link_type=DOI) 12. 12.Reus LM, Shen X, Gibson J, et al. Association of polygenic risk for major psychiatric illness with subcortical volumes and white matter integrity in UK Biobank. Sci Rep. 2017;7(1):42140. doi:10.1038/srep42140 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/srep42140&link_type=DOI) 13. 13.Whalley HC, Sprooten E, Hackett S, et al. Polygenic Risk and White Matter Integrity in Individuals at High Risk of Mood Disorder. Biol Psychiatry. 2013;74(4):280–286. doi:10.1016/j.biopsych.2013.01.027 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2013.01.027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23453289&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000323156100009&link_type=ISI) 14. 14.Gottesman II, Gould TD. The Endophenotype Concept in Psychiatry: Etymology and Strategic Intentions. Am J Psychiatry. 2003;160(4):636–645. doi:10.1176/appi.ajp.160.4.636 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.160.4.636&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12668349&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000182096300008&link_type=ISI) 15. 15.Smith DJ, Nicholl BI, Cullen B, et al. Prevalence and Characteristics of Probable Major Depression and Bipolar Disorder within UK Biobank: Cross-Sectional Study of 172,751 Participants. Potash JB, ed. PLoS ONE. 2013;8(11):e75362. doi:10.1371/journal.pone.0075362 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0075362&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24282498&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 16. 16.Alfaro-Almagro F, Jenkinson M, Bangerter NK, et al. Image processing and Quality Control for the first 10,000 brain imaging datasets from UK Biobank. NeuroImage. 2018;166:400–424. doi:10.1016/j.neuroimage.2017.10.034 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2017.10.034&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29079522&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 17. 17.Cox SR, Ritchie SJ, Tucker-Drob EM, et al. Ageing and brain white matter structure in 3,513 UK Biobank participants. Nat Commun. 2016;7(1):13629. doi:10.1038/ncomms13629 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms13629&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27976682&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 18. 18.Shen X, Reus LM, Cox SR, et al. Subcortical volume and white matter integrity abnormalities in major depressive disorder: findings from UK Biobank imaging data. Sci Rep. 2017;7(1). doi:10.1038/s41598-017-05507-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-017-05507-6&link_type=DOI) 19. 19.Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300. doi:10.1111/j.2517-6161.1995.tb02031.x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.2517-6161.1995.tb02031.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24443148&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 20. 20.Bulik-Sullivan BK, Loh PR, Finucane HK, et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet. 2015;47(3):291–295. doi:10.1038/ng.3211 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3211&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 21. 21.Deng F, Wang Y, Huang H, et al. Abnormal segments of right uncinate fasciculus and left anterior thalamic radiation in major and bipolar depression. Prog Neuropsychopharmacol Biol Psychiatry. 2018;81:340–349. doi:10.1016/j.pnpbp.2017.09.006 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pnpbp.2017.09.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28912043&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 22. 22.Papez J. A proposed mechanism of emotion. 1937 [classical article]. J Neuropsychiatry Clin Neurosci. 1995;7(1):103–112. doi:10.1176/jnp.7.1.103 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/jnp.7.1.103&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7711480&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QD72700019&link_type=ISI) 23. 23.Bracht T, Linden D, Keedwell P. A review of white matter microstructure alterations of pathways of the reward circuit in depression. J Affect Disord. 2015;187:45–53. doi:10.1016/j.jad.2015.06.041 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jad.2015.06.041&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26318270&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 24. 24.Pettersson E, Lichtenstein P, Larsson H, et al. Genetic influences on eight psychiatric disorders based on family data of 4 408 646 full and half-siblings, and genetic data of 333 748 cases and controls. Psychol Med. 2019;49(07):1166–1173. doi:10.1017/S0033291718002039 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291718002039&link_type=DOI) 25. 25.Murray GK, Lin T, Austin J, McGrath JJ, Hickie IB, Wray NR. Could Polygenic Risk Scores Be Useful in Psychiatry?: A Review. JAMA Psychiatry. 2021;78(2):210. doi:10.1001/jamapsychiatry.2020.3042 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamapsychiatry.2020.3042&link_type=DOI) 26. 26.Grama S, Willcocks I, Hubert JJ, et al. Polygenic risk for schizophrenia and subcortical brain anatomy in the UK Biobank cohort. Transl Psychiatry. 2020;10(1):309. doi:10.1038/s41398-020-00940-0 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41398-020-00940-0&link_type=DOI) 27. 27.Ranlund S, Rosa MJ, de Jong S, et al. Associations between polygenic risk scores for four psychiatric illnesses and brain structure using multivariate pattern recognition. NeuroImage Clin. 2018;20:1026–1036. doi:10.1016/j.nicl.2018.10.008 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.nicl.2018.10.008&link_type=DOI) 28. 28.Sha Z, Schijven D, Francks C. Patterns of brain asymmetry associated with polygenic risks for autism and schizophrenia implicate language and executive functions but not brain masculinization. Mol Psychiatry. Published online July 1, 2021. doi:10.1038/s41380-021-01204-z [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41380-021-01204-z&link_type=DOI) 29. 29.Stauffer EM, Bethlehem RAI, Warrier V, et al. Grey and white matter microstructure is associated with polygenic risk for schizophrenia. Mol Psychiatry. Published online August 30, 2021. doi:10.1038/s41380-021-01260-5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41380-021-01260-5&link_type=DOI) 30. 30.Zhu X, Ward J, Cullen B, et al. Polygenic Risk for Schizophrenia, Brain Structure, and Environmental Risk in UK Biobank. Schizophr Bull Open. 2021;2(1):sgab042. doi:10.1093/schizbullopen/sgab042 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/schizbullopen/sgab042&link_type=DOI) 31. 31.Alemany S, Blok E, Jansen PR, Muetzel RL, White T. Brain morphology, autistic traits, and polygenic risk for autism: A p OPULATION-BASED neuroimaging study. Autism Res. 2021;14(10):2085–2099. doi:10.1002/aur.2576 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/aur.2576&link_type=DOI) 32. 32.van der Merwe C, Passchier R, Mufford M, Ramesar R, Dalvie S, Stein DJ. Polygenic risk for schizophrenia and associated brain structural changes: A systematic review. Compr Psychiatry. 2019;88:77–82. doi:10.1016/j.comppsych.2018.11.014 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.comppsych.2018.11.014&link_type=DOI) 33. 33.Schmitt S, Meller T, Stein F, et al. Effects of polygenic risk for major mental disorders and cross-disorder on cortical complexity. Psychol Med. Published online April 8, 2021:1–12. doi:10.1017/S0033291721001082 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291721001082&link_type=DOI) 34. 34.Shen X, Howard DM, Adams MJ, et al. A phenome-wide association and Mendelian Randomisation study of polygenic risk for depression in UK Biobank. Nat Commun. 2020;11(1):2301. doi:10.1038/s41467-020-16022-0 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-16022-0&link_type=DOI) 35. 35.Ward J, Lyall LM, Bethlehem RAI, et al. Novel genome-wide associations for anhedonia, genetic correlation with psychiatric disorders, and polygenic association with brain structure. Transl Psychiatry. 2019;9(1):327. doi:10.1038/s41398-019-0635-y [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41398-019-0635-y&link_type=DOI) 36. 36.Zhang X, Yan H, Yu H, et al. Childhood urbanicity interacts with polygenic risk for depression to affect stress-related medial prefrontal function. Transl Psychiatry. 2021;11(1):522. doi:10.1038/s41398-021-01650-x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41398-021-01650-x&link_type=DOI) 37. 37.Kendler KS, Neale MC. Endophenotype: a conceptual analysis. Mol Psychiatry. 2010;15(8):789–797. doi:10.1038/mp.2010.8 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2010.8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20142819&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280854800006&link_type=ISI) 38. 38.Meyer-Lindenberg A, Weinberger DR. Intermediate phenotypes and genetic mechanisms of psychiatric disorders. Nat Rev Neurosci. 2006;7(10):818–827. doi:10.1038/nrn1993 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn1993&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16988657&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241191800014&link_type=ISI) 39. 39.Zhang Y, Qi G, Park JH, Chatterjee N. Estimation of complex effect-size distributions using summary-level statistics from genome-wide association studies across 32 complex traits. Nat Genet. 2018;50(9):1318–1326. doi:10.1038/s41588-018-0193-x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0193-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30104760&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 40. 40.Wray NR, Yang J, Hayes BJ, Price AL, Goddard ME, Visscher PM. Pitfalls of predicting complex traits from SNPs. Nat Rev Genet. 2013;14(7):507–515. doi:10.1038/nrg3457 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg3457&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23774735&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 41. 41.Plomin R, von Stumm S. The new genetics of intelligence. Nat Rev Genet. 2018;19:148–159. doi:10.1038/nrg.2017.104 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg.2017.104&link_type=DOI) 42. 42.van der Meer D, Shadrin AA, O’Connell K, et al. Boosting Schizophrenia Genetics by Utilizing Genetic Overlap With Brain Morphology. Biol Psychiatry. 2022;92(4):291–298. doi:10.1016/j.biopsych.2021.12.007 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2021.12.007&link_type=DOI) 43. 43.Drobisz D, Damborská A. Deep brain stimulation targets for treating depression. Behav Brain Res. 2019;359:266–273. doi:10.1016/j.bbr.2018.11.004 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbr.2018.11.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30414974&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 44. 44.Riva-Posse P, Choi KS, Holtzheimer PE, et al. A connectomic approach for subcallosal cingulate deep brain stimulation surgery: prospective targeting in treatment-resistant depression. Mol Psychiatry. 2018;23(4):843–849. doi:10.1038/mp.2017.59 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2017.59&link_type=DOI) 45. 45.Korgaonkar MS, Rekshan W, Gordon E, et al. Magnetic Resonance Imaging Measures of Brain Structure to Predict Antidepressant Treatment Outcome in Major Depressive Disorder. EBioMedicine. 2015;2(1):37–45. doi:10.1016/j.ebiom.2014.12.002 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ebiom.2014.12.002&link_type=DOI) 46. 46.Gong Q, He Y. Depression, Neuroimaging and Connectomics: A Selective Overview. Biol Psychiatry. 2015;77(3):223–235. doi:10.1016/j.biopsych.2014.08.009 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2014.08.009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25444171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom) 47. 47.Korgaonkar MS, Cooper NJ, Williams LM, Grieve SM. Mapping inter-regional connectivity of the entire cortex to characterize major depressive disorder: a whole-brain diffusion tensor imaging tractography study. NeuroReport. 2012;23(9):566–571. doi:10.1097/WNR.0b013e3283546264 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/WNR.0b013e3283546264&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22562047&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F14%2F2022.11.09.22282111.atom)