Major depressive disorder (MDD) is a prevalent and frequently disabling psychiatric disorder, associated with prolonged low mood, although many unique symptom combinations may lead to the same diagnosis (Gallo and Rabins, Reference Gallo and Rabins1999; Kennedy, Reference Kennedy2008). The high clinical heterogeneity suggests that MDD may group together a number of disease subtypes (Fava et al., Reference Fava, Uebelacker, Alpert, Nierenberg, Pava and Rosenbaum1997; Baumeister and Parker, Reference Baumeister and Parker2012; van Loo et al., Reference van Loo, Cai, Gruber, Li, de Jonge, Petukhova, Rose, Sampson, Schoevers, Wardenaar, Wilcox, Al-Hamzawi, Andrade, Bromet, Bunting, Fayyad, Florescu, Gureje, Hu, Huang, Levinson, Medina-Mora, Nakane, Posada-Villa, Scott, Xavier, Zarkov and Kessler2014). Differences in clinical presentation may reflect differences in aetiology and underlying neurobiological mechanisms. Identifying subtypes of MDD is therefore likely to be a crucial step toward more effective diagnosis and treatment of affected individuals, as well as providing better targets for genomic and neurobiological studies.
The heritability of MDD is estimated to be about 37% (Sullivan et al., Reference Sullivan, Neale and Kendler2000), determined by a large number of alleles each of small effect (Ripke et al., Reference Ripke, Wray, Lewis, Hamilton, Weissman, Breen, Byrne, Blackwood, Boomsma, Cichon, Heath, Holsboer, Lucae, Madden, Martin, McGuffin, Muglia, Noethen, Penninx, Pergadia, Potash, Rietschel, Lin, Müller-Myhsok, Shi, Steinberg, Grabe, Lichtenstein, Magnusson, Perlis, Preisig, Smoller, Stefansson, Uher, Kutalik, Tansey, Teumer, Viktorin, Barnes, Bettecken, Binder, Breuer, Castro, Churchill, Coryell, Craddock, Craig, Czamara, De Geus, Degenhardt, Farmer, Fava, Frank, Gainer, Gallagher, Gordon, Goryachev, Gross, Guipponi, Henders, Herms, Hickie, Hoefels, Hoogendijk, Hottenga, Iosifescu, Ising, Jones, Jones, Jung-Ying, Knowles, Kohane, Kohli, Korszun, Landen, Lawson, Lewis, Macintyre, Maier, Mattheisen, McGrath, McIntosh, McLean, Middeldorp, Middleton, Montgomery, Murphy, Nauck, Nolen, Nyholt, O'Donovan, Oskarsson, Pedersen, Scheftner, Schulz, Schulze, Shyn, Sigurdsson, Slager, Smit, Stefansson, Steffens, Thorgeirsson, Tozzi, Treutlein, Uhr, van den Oord, Van Grootheest, Völzke, Weilburg, Willemsen, Zitman, Neale, Daly, Levinson and Sullivan2013). MDD subtypes may differ in terms of which alleles (Flint and Kendler, Reference Flint and Kendler2014) and which biological mechanisms (Hasler et al., Reference Hasler, Drevets, Manji and Charney2004) are involved. For example, Milaneschi et al. (Reference Milaneschi, Lamers, Peyrot, Baune, Breen, Dehghan, Forstner, Grabe, Homuth, Kan, Lewis, Mullins, Nauck, Pistis, Preisig, Rivera, Rietschel, Streit, Strohmaier, Teumer, Van der Auwera, Wray, Boomsma and Penninx2017) demonstrated only a very small genetic overlap between typical MDD (associated with decreased appetite and weight) and an atypical subtype (associated with increased appetite and weight). Focusing more specifically on stratifying MDD by genetic factors, Howard et al. (Reference Howard, Clarke, Adams, Hafferty, Wigmore, Zeng, Hall, Gibson, Boutin, Hayward, Thomson, Porteous, Smith, Murray, Haley, Deary, Whalley and McIntosh2017) recently identified 10 variables related to genetic subgroups that might partly account for the heterogeneity of MDD. Further investigation of specific genetic subgroups may help to identify more clinically and empirically useful MDD subtypes.
Some genes that contribute to MDD overlap with those underlying schizophrenia (SCZ; Schulze et al., Reference Schulze, Akula, Breuer, Steele, Nalls, Singleton, Degenhardt, Nöthen, Cichon, Rietschel and McMahon2014) and the two disorders have a genetic correlation of r g = 0.43 (Cross-Disorder Group of the Psychiatric Genetics Consortium et al., Reference Lee, Ripke, Neale, Faraone, Purcell, Perlis, Mowry, Thapar, Goddard, Witte, Absher, Agartz, Akil, Amin, Andreassen, Anjorin, Anney, Anttila, Arking, Asherson, Azevedo, Backlund, Badner, Bailey, Banaschewski, Barchas, Barnes, Barrett, Bass, Battaglia, Bauer, Bayés, Bellivier, Bergen, Berrettini, Betancur, Bettecken, Biederman, Binder, Black, Blackwood, Bloss, Boehnke, Boomsma, Breen, Breuer, Bruggeman, Cormican, Buccola, Buitelaar, Bunney, Buxbaum, Byerley, Byrne, Caesar, Cahn, Cantor, Casas, Chakravarti, Chambert, Choudhury, Cichon, Cloninger, Collier, Cook, Coon, Cormand, Corvin, Coryell, Craig, Craig, Crosbie, Cuccaro, Curtis, Czamara, Datta, Dawson, Day, De Geus, Degenhardt, Djurovic, Donohoe, Doyle, Duan, Dudbridge, Duketis, Ebstein, Edenberg, Elia, Ennis, Etain, Fanous, Farmer, Ferrier, Flickinger, Fombonne, Foroud, Frank, Franke, Fraser, Freedman, Freimer, Freitag, Friedl, Frisén, Gallagher, Gejman, Georgieva, Gershon, Geschwind, Giegling, Gill, Gordon, Gordon-Smith, Green, Greenwood, Grice, Gross, Grozeva, Guan, Gurling, De Haan, Haines, Hakonarson, Hallmayer, Hamilton, Hamshere, Hansen, Hartmann, Hautzinger, Heath, Henders, Herms, Hickie, Hipolito, Hoefels, Holmans, Holsboer, Hoogendijk, Hottenga, Hultman, Hus, Ingason, Ising, Jamain, Jones, Jones, Jones, Tzeng, Kähler, Kahn, Kandaswamy, Keller, Kennedy, Kenny, Kent, Kim, Kirov, Klauck, Klei, Knowles, Kohli, Koller, Konte, Korszun, Krabbendam, Krasucki, Kuntsi, Kwan, Landén, Långström, Lathrop, Lawrence, Lawson, Leboyer, Ledbetter, Lee, Lencz, Lesch, Levinson, Lewis, Li, Lichtenstein, Lieberman, Lin, Linszen, Liu, Lohoff, Loo, Lord, Lowe, Lucae, MacIntyre, Madden, Maestrini, Magnusson, Mahon, Maier, Malhotra, Mane, Martin, Martin, Mattheisen, Matthews, Mattingsdal, McCarroll, McGhee, McGough, McGrath, McGuffin, McInnis, McIntosh, McKinney, McLean, McMahon, McMahon, McQuillin, Medeiros, Medland, Meier, Melle, Meng, Meyer, Middeldorp, Middleton, Milanova, Miranda, Monaco, Montgomery, Moran, Moreno-De-Luca, Morken, Morris, Morrow, Moskvina, Muglia, Mühleisen, Muir, Müller-Myhsok, Murtha, Myers, Myin-Germeys, Neale, Nelson, Nievergelt, Nikolov, Nimgaonkar, Nolen, Nöthen, Nurnberger, Nwulia, Nyholt, O'Dushlaine, Oades, Olincy, Oliveira, Olsen, Ophoff, Osby, Owen, Palotie, Parr, Paterson, Pato, Pato, Penninx, Pergadia, Pericak-Vance, Pickard, Pimm, Piven, Posthuma, Potash, Poustka, Propping, Puri, Quested, Quinn, Ramos-Quiroga, Rasmussen, Raychaudhuri, Rehnström, Reif, Ribasés, Rice, Rietschel, Roeder, Roeyers, Rossin, Rothenberger, Rouleau, Ruderfer, Rujescu, Sanders, Sanders, Santangelo, Sergeant, Schachar, Schalling, Schatzberg, Scheftner, Schellenberg, Scherer, Schork, Schulze, Schumacher, Schwarz, Scolnick, Scott, Shi, Shilling, Shyn, Silverman, Slager, Smalley, Smit, Smith, Sonuga-Barke, St Clair, State, Steffens, Steinhausen, Strauss, Strohmaier, Stroup, Sutcliffe, Szatmari, Szelinger, Thirumalai, Thompson, Todorov, Tozzi, Treutlein, Uhr, van den Oord, Van Grootheest, Van Os, Vicente, Vieland, Vincent, Visscher, Walsh, Wassink, Watson, Weissman, Werge, Wienker, Wijsman, Willemsen, Williams, Willsey, Witt, Xu, Young, Yu, Zammit, Zandi, Zhang, Zitman, Zöllner, Devlin, Kelsoe, Sklar, Daly, O'Donovan, Craddock, Sullivan, Smoller, Kendler and Wray2013). Individuals with MDD also at high polygenic risk of SCZ have been previously shown to differ in terms of clinical and behavioural phenotypes. Whalley et al. (Reference Whalley, Adams, Hall, Clarke, Fernandez-Pujals, Gibson, Wigmore, Hafferty, Hagenaars, Davies, Campbell, Hayward, Lawrie, Porteous, Deary and McIntosh2016) observed higher neuroticism and psychological distress among controls with a higher SCZ polygenic risk score (PRS), but not among high-SCZ-PRS MDD cases. Additionally, Power et al. (Reference Power, Tansey, Buttenschøn, Cohen-Woods, Bigdeli, Hall, Kutalik, Lee, Ripke, Steinberg, Teumer, Viktorin, Wray, Arolt, Baune, Boomsma, Børglum, Byrne, Castelao, Craddock, Craig, Dannlowski, Deary, Degenhardt, Forstner, Gordon, Grabe, Grove, Hamilton, Hayward, Heath, Hocking, Homuth, Hottenga, Kloiber, Krogh, Landén, Lang, Levinson, Lichtenstein, Lucae, MacIntyre, Madden, Magnusson, Martin, McIntosh, Middeldorp, Milaneschi, Montgomery, Mors, Müller-Myhsok, Nyholt, Oskarsson, Owen, Padmanabhan, Penninx, Pergadia, Porteous, Potash, Preisig, Rivera, Shi, Shyn, Sigurdsson, Smit, Smith, Stefansson, Stefansson, Strohmaier, Sullivan, Thomson, Thorgeirsson, Van der Auwera, Weissman, Breen and Lewis2017) found that MDD subtypes based on age of onset differed in terms of SCZ PRS, and early-onset cases have been found to show greater deficits in frontal (Jaworska et al., Reference Jaworska, MacMaster, Gaxiola, Cortese, Goodyear and Ramasubbu2014) and limbic (MacMaster et al., Reference MacMaster, Mirza, Szeszko, Kmiecik, Easter, Taormina, Lynch, Rose, Moore and Rosenberg2008; Clark et al., Reference Clark, Konduru, Kemp, Bray, Brown, Goodyear and Ramasubbu2018) regions. This may suggest a different aetiology and perhaps different neurobiological profile of MDD in the presence of high genetic risk for SCZ, i.e. an SCZ-risk-related MDD subtype. It seems that some of the genes that may contribute to both SCZ and MDD are involved in regulating synaptic function and the excitability of prefrontal neurons (Howard et al., Reference Howard, Adams, Clarke, Hafferty, Gibson, Shirali, Coleman, Hagenaars, Ward, Wigmore, Alloza, Shen, Barbu, Xu, Whalley, Marioni, Porteous, Davies, Deary, Hemani, Berger, Teismann, Rawal, Arolt, Baune, Dannlowski, Domschke, Tian, Hinds, Trzaskowski, Byrne, Ripke, Smith, Sullivan, Wray, Breen, Lewis and McIntosh2019). However, whether subtypes of MDD based on genetic loading for SCZ are distinct in terms of brain structure has yet to be determined.
MDD is typically associated with fronto-limbic deficits, including reduced cortical volume/thickness in prefrontal regions (Salvadore et al., Reference Salvadore, Nugent, Lemaitre, Luckenbaugh, Tinsley, Cannon, Neumeister, Zarate and Drevets2011; Rodriguez-Cano et al., Reference Rodríguez-Cano, Sarró, Monté, Maristany, Salvador, McKenna and Pomarol-Clotet2014; Schmaal et al., Reference Schmaal, Hibar, Sämann, Hall, Baune, Jahanshad, Cheung, van Erp, Bos, Ikram, Vernooij, Niessen, Tiemeier, Hofman, Wittfeld, Grabe, Janowitz, Bülow, Selonke, Völzke, Grotegerd, Dannlowski, Arolt, Opel, Heindel, Kugel, Hoehn, Czisch, Couvy-Duchesne, Rentería, Strike, Wright, Mills, de Zubicaray, McMahon, Medland, Martin, Gillespie, Goya-Maldonado, Gruber, Krämer, Hatton, Lagopoulos, Hickie, Frodl, Carballedo, Frey, van Velzen, Penninx, van Tol, van der Wee, Davey, Harrison, Mwangi, Cao, Soares, Veer, Walter, Schoepf, Zurowski, Konrad, Schramm, Normann, Schnell, Sacchet, Gotlib, MacQueen, Godlewska, Nickson, McIntosh, Papmeyer, Whalley, Hall, Sussmann, Li, Walter, Aftanas, Brack, Bokhan, Thompson and Veltman2017), and smaller subcortical limbic structures (Kim et al., Reference Kim, Hamilton and Gotlib2008; Lu et al., Reference Lu, Liang, Han, Mo, Li, Cheng, Xu, Shen, Tan, Zhao, Zhu and Sun2016; Schmaal et al., Reference Schmaal, Veltman, van Erp, Sämann, Frodl, Jahanshad, Loehrer, Tiemeier, Hofman, Niessen, Vernooij, Ikram, Wittfeld, Grabe, Block, Hegenscheid, Völzke, Hoehn, Czisch, Lagopoulos, Hatton, Hickie, Goya-Maldonado, Krämer, Gruber, Couvy-Duchesne, Rentería, Strike, Mills, de Zubicaray, McMahon, Medland, Martin, Gillespie, Wright, Hall, MacQueen, Frey, Carballedo, van Velzen, van Tol, van der Wee, Veer, Walter, Schnell, Schramm, Normann, Schoepf, Konrad, Zurowski, Nickson, McIntosh, Papmeyer, Whalley, Sussmann, Godlewska, Cowen, Fischer, Rose, Penninx, Thompson and Hibar2016). Reductions in integrity of connecting white matter tracts have also been observed (Cole et al., Reference Cole, Chaddock, Farmer, Aitchison, Simmons, McGuffin and Fu2012; Bessette et al., Reference Bessette, Nave, Caprihan and Stevens2014; Shen et al., Reference Shen, Reus, Cox, Adams, Liewald, Bastin, Smith, Deary, Whalley and McIntosh2017). Fronto-limbic impairments are also apparent in SCZ (Ross et al., Reference Ross, Margolis, Reading, Pletnikov and Coyle2006; Keshavan et al., Reference Keshavan, Tandon, Boutros and Nasrallah2008; Yao et al., Reference Yao, Lui, Liao, Du, Hu, Thomas and Gong2013), along with more widespread white matter deficits (Ellison-Wright and Bullmore, Reference Ellison-Wright and Bullmore2009) and greater cortical reductions in temporal regions (Wong and van Tol, Reference Wong and van Tol2003). Furthermore, some SCZ-like neuroimaging markers, such as reduced hippocampal, thalamic and prefrontal volumes, appear in individuals who are genetically predisposed to SCZ but remain healthy (Lawrie et al., Reference Lawrie, Whalley, Abukmeil, Kestelman, Donnelly, Miller, Best, Owens and Johnstone2001; Boos et al., Reference Boos, Aleman, Cahn, Hulshoff Pol and Kahn2007). As of yet, the influences of SCZ risk on such neurobiological measures in the presence of MDD are unclear.
Following on from previous work stratifying MDD traits by PRS for SCZ (Whalley et al., Reference Whalley, Adams, Hall, Clarke, Fernandez-Pujals, Gibson, Wigmore, Hafferty, Hagenaars, Davies, Campbell, Hayward, Lawrie, Porteous, Deary and McIntosh2016), we aimed to measure the neuroimaging associations of MDD as a function of SCZ PRS, testing for further evidence of MDD subtype heterogeneity related to SCZ risk. We assessed interactive effects between SCZ PRS and MDD case/control status on measures of cortical, subcortical and white matter structure. Based on previous evidence of neuroimaging abnormalities in healthy subjects at high genetic risk for SCZ (Lawrie et al., Reference Lawrie, Whalley, Abukmeil, Kestelman, Donnelly, Miller, Best, Owens and Johnstone2001; Boos et al., Reference Boos, Aleman, Cahn, Hulshoff Pol and Kahn2007), we expected negative correlations between SCZ PRS and brain measures among both cases and controls, with significant interactions between MDD and SCZ PRS indicating a specific subtype. Without any specific hypotheses regarding the neurobiological bases of differences in neuroticism and distress among high-SCZ-risk MDD cases, we took an exploratory approach, testing SCZ PRS by MDD interactions for structural measures covering the entire brain. We did, however, anticipate interactions in areas that normally show deficits in MDD, such as frontal cortical regions, limbic structures and connecting white matter tracts. Any significant interactions, indicating that MDD's associations with neuroimaging measures differed as a function of SCZ PRS, would provide evidence of an SCZ-risk-related MDD subtype.
Methods
Participants
Data were derived from 4944 participants in the large-scale prospective epidemiological study, UK Biobank (http://www/ukbiobank.ac.uk). UK Biobank includes data on 503 325 members of the general UK population, recruited between 2006 and 2010 (Allen et al., Reference Allen, Sudlow, Downey, Peakman, Danesh, Elliot, Gallacher, Green, Matthews, Pell, Sprosen and Collins2012). Participants originally provided information on a wide range of health, lifestyle, environment and other variables, and the majority have also provided genetic data. The project aims to submit approximately 100 000 of the participants to a full-body imaging protocol, which includes the acquisition of neuroimaging data in several modalities (Miller et al., Reference Miller, Alfaro-Almagro, Bangerter, Thomas, Yacoub, Xu, Bartsch, Jbabdi, Sotiropoulos, Andersson, Griffanti, Douaud, Okell, Weale, Dragonu, Garratt, Hudson, Collins, Jenkinson, Matthews and Smith2016). The present study focused on 4944 UK Biobank participants for whom both genetic and brain imaging data have already been acquired, who also had data on MDD and relevant covariates, described below. Participants who reported a diagnosis of schizophrenia were excluded from all analyses. Descriptive statistics for included participants are reported in Table 1.
Descriptive statistics for MDD cases and controls, and group differences
MDD, major depressive disorder; WBV, whole brain volume; SCZ, schizophrenia; PRS, polygenic risk score; DTI, diffusion tensor imaging.
Group differences were assessed using χ2 tests and two-sample t tests.
*, ** and *** represent significant results at p < 0.05, p < 0.01 and p < 0.001, respectively
Polygenic risk for schizophrenia
UK Biobank blood samples were genotyped using either the UK BiLEVE array or UK Biobank Axiom array, and quality controlled at the Wellcome Trust Centre for Human Genetics (Wain et al., Reference Wain, Shrine, Miller, Jackson, Ntalla, Soler Artigas, Billington, Kheirallah, Allen, Cook, Probert, Obeidat, Bossé, Hao, Postma, Paré, Ramasamy, Mägi, Mihailov, Reinmaa, Melén, O'Connell, Frangou, Delaneau, Freeman, Petkova, McCarthy, Sayers, Deloukas, Hubbard, Pavord, Hansell, Thomson, Zeggini, Morris, Marchini, Strachan, Tobin and Hall2015; Hagenaars et al., Reference Hagenaars, Harris, Davies, Hill, Liewald, Ritchie, Marioni, Fawns-Ritchie, Cullen, Malik, Worrall, Sudlow, Wardlaw, Gallacher, Pell, McIntosh, Smith, Gale and Deary2016). Further quality control measures included removal of participants based on missingness, relatedness, non-British ancestry and gender mismatch. Subsequent processing involved removal of SNPs with minor allele frequency <1%, and clump-based linkage disequilibrium pruning (r 2 < 0.25 within a 200 bp window). PRSs for SCZ were calculated for each participant using PRSice (Euesden et al., Reference Euesden, Lewis and O'Reilly2015) and genome-wide association study summary data (Schizophrenia Working Group of the Psychiatric Genetics Consortium, 2014), including SNPs selected according to the significance of their associations with the phenotype at thresholds of p < 0.01, p < 0.05, p < 0.10, p < 0.50 and all SNPs.
Major depressive disorder
As part of a web-based mental health questionnaire, UK Biobank participants completed the short form of the Composite International Diagnostic Interview (CIDI-SF; Kessler et al., Reference Kessler, Andrews, Mroczek, Ustun and Wittchen1998). Section A of the CIDI-SF focuses on MDD symptomatology, providing a score of between zero and eight. Participants with a score of five or more were classified as lifelong MDD cases; those with a score of four or less as controls. Participants who scored highly but did not report having experienced core symptoms of either low mood or anhedonia were excluded. CIDI-SF controls who had previously reported a past diagnosis of depression or showed signs of MDD on the 9-item Patient Health Questionnaire (PHQ-9; Kroenke et al., Reference Kroenke, Spitzer and Williams2001) were also excluded (N = 286). The resulting CIDI-assessed measure of MDD therefore included 1674 cases and 3270 controls (Table 1). Some of the earlier scanned subjects completed the web-based questionnaire up to 2.7 years later, but most completed it within a year of imaging data acquisition.
Brain imaging data
Full information on imaging data acquisition and processing are provided in supplementary materials and previous publications (Miller et al., Reference Miller, Alfaro-Almagro, Bangerter, Thomas, Yacoub, Xu, Bartsch, Jbabdi, Sotiropoulos, Andersson, Griffanti, Douaud, Okell, Weale, Dragonu, Garratt, Hudson, Collins, Jenkinson, Matthews and Smith2016; Alfaro-Almagro et al., Reference Alfaro-Almagro, Jenkinson, Bangerter, Andersson, Griffanti, Douaud, Sotiropoulos, Jbabdi, Hernandez-Fernandez, Vallee, Vidaurre, Webster, McCarthy, Rorden, Daducci, Alexander, Zhang, Dragonu, Matthews, Miller and Smith2018). To summarise, all structural and diffusion data were acquired on the same 3T scanner using the same protocol. Cortical metrics for 27 regions (N = 1877) were derived locally using FreeSurfer version 5.3 (Dale et al., Reference Dale, Fischl and Sereno1999; Fischl et al., Reference Fischl, Sereno and Dale1999, Reference Fischl, van der Kouwe, Destrieux, Halgren, Ségonne, Salat, Busa, Seidman, Goldstein, Kennedy, Caviness, Makris, Rosen and Dale2004; Desikan et al., Reference Desikan, Ségonne, Fischl, Quinn, Dickerson, Blacker, Buckner, Dale, Maguire, Hyman, Albert and Killiany2006). Volumes of seven subcortical structures (N = 4944) were derived by the UK Biobank imaging team using FIRST (Patenaude et al., Reference Patenaude, Smith, Kennedy and Jenkinson2011). Fractional anisotropy (FA) and mean diffusivity (MD) of 15 white matter tracts (N = 4420) were also derived by the UK Biobank imaging team using bedpostx and probtrackx (Behrens et al., Reference Behrens, Berg, Jbadi, Rushworth and Woolrich2007), and then AutoPtx (de Groot et al., Reference de Groot, Vernooij, Klein, Ikram, Vos, Smith, Niessen and Andersson2013).
At the time of beginning this study, subcortical volumes and white matter metrics were available for the first two releases of UK Biobank imaging data, while cortical metrics were available for only the first release. Following initial submission of our results for publication, cortical metrics became available for the second release as well, including for 794 MDD cases and 1604 controls who met the criteria for inclusion in this study. We have therefore used these additional data to replicate any significant findings discovered in the first release.
Additional covariates
In relation to SCZ PRS, genotyping array and the first 15 ancestry principal components (controlling for population stratification within the sample) were included as covariates in all analyses. We used as many as 15 principal components to control for the substantial genetic clustering recently found in UK Biobank (Abdellaoui et al., Reference Abdellaoui, Hugh-Jones, Kemper, Holtz, Nivard, Veul, Yengo, Zietsch, Frayling, Wray, Yang, Verweij and Visscher2018). Including more than 15 made negligible difference to our results. In relation to neuroimaging data, x, y and z coordinates of head position within the MRI scanner were included as covariates, along with a measure of whole-brain volume (including grey matter, white matter and ventricles), derived from the UK Biobank image processing pipeline (Alfaro-Almagro et al., Reference Alfaro-Almagro, Jenkinson, Bangerter, Andersson, Griffanti, Douaud, Sotiropoulos, Jbabdi, Hernandez-Fernandez, Vallee, Vidaurre, Webster, McCarthy, Rorden, Daducci, Alexander, Zhang, Dragonu, Matthews, Miller and Smith2018). Additionally, age in years at time of scanning, age2 and sex were included in all analyses.
Statistical analysis
Data were analysed in R version 3.2.3 (R Core Team, 2013). As additional quality control, outliers, defined as further than three standard deviations from the mean, were removed from all neuroimaging measures. Left and right measures were then combined bilaterally, as above; a general lack of hemispheric effects in the same UK Biobank subcortical and white matter phenotypes has been reported previously (Shen et al., Reference Shen, Reus, Cox, Adams, Liewald, Bastin, Smith, Deary, Whalley and McIntosh2017). We performed analyses of case-control differences in group size, sex and SCZ PRS, reported in Table 1. For each regional brain measure and each PRS threshold, we then ran linear regression models to estimate the strength of the interactive effect between SCZ PRS and MDD in all participants, controlling for the main effects of each, as well as effects of age, age2, genotyping array, the first 15 principal components of genetic ancestry, whole-brain volume and scanner head position coordinates. FDR correction was applied to p values across all interaction effects together.
As main effects of SCZ PRS among MDD cases would not demonstrate whether the effects were specific to cases, only interactions between MDD and SCZ PRS could provide indications of aetiological differences between MDD cases at high and low SCZ risk. However, where interactions were significant, we assessed main effects of SCZ PRS among MDD cases and controls separately, and tested case-control differences among those above and below mean SCZ PRS. Together, these analyses provided further information on how MDD and SCZ risk interact. We additionally tested whether significant interactions were influenced by a subgroup of MDD cases with both repeated k-mean and hierarchical cluster analyses, using package ‘fpc’ (Henning, Reference Henning2015) in addition to base packages. Further packages ‘ggplot2’ (Wickham and Chang, Reference Wickham and Chang2016) and ‘ellipse’ (Murdoch and Chow, Reference Murdoch and Chow2013) were also used to create figures. Results are reported in terms of standardised β coefficients, and p values and q values (FDR-corrected p values) <0.05 are considered significant.
Results
Descriptive statistics are reported for MDD cases and controls in Table 1. MDD cases were younger than controls (t 3252 = −11.15, p < 0.001), a higher proportion of cases were female (χ21 = 113.33, p < 0.001), and WBV was lower among cases (t 3460 = 5.11, p < 0.001). These variables were controlled in all subsequent analyses. Mean PRS was slightly higher among cases at thresholds 0.10 (t 3484 = 2.28, p = 0.02), 0.50 (t 3430 = 2.91, p < 0.01) and all SNPs (t 3427 = 3.01, p < 0.01). All SCZ PRS by MDD interactions tested for both MDD measures are summarised in Fig. 1; within each measure, fronto-limbic areas are presented first.
Significance of SCZ PRS by MDD interactions on all cortical, subcortical and white matter metrics. SCZ, schizophrenia; PRS, polygenic risk score; MDD, major depressive disorder; CT, cortical thickness; CSA, cortical surface area; CV, cortical volume; SCV, subcortical volume; FA, fractional anisotropy; MD, mean diffusivity. Results for fronto-limbic brain regions are presented first in red/darker grey, followed by all other brain regions in blue/lighter grey. Points above the lower dashed line represent results that were significant before correcting for multiple comparisons (p < 0.05). Points representing significant results after FDR correction (q < 0.05) are circled.
Cortical regions
Analyses of cortical metrics focused on a subset of 1877 subjects for whom these data were available. We first assessed interactive effects between SCZ PRS and MDD on regional mean cortical thickness. There were significant interactive effects across multiple PRS thresholds of at least β = 0.132 (p = 0.009) in rostral ACC, β = 0.101 (p = 0.042) in superior temporal gyrus and β = 0.120 (p = 0.037) in entorhinal cortex (Table 2). The interaction in rostral ACC was significant across all thresholds and the two strongest of these effects remained significant after FDR correction: for PRS threshold 0.50 (β = 0.191, q = 0.043) and for all SNPs (β = 0.191, q = 0.043). These interactions were driven by significant negative effects of β = −0.094 (p = 0.001) and β = −0.090 (p = 0.002) among controls, but positive effects of β = 0.103 (p = 0.017) and β = 0.109 (p = 0.013) among cases. These results mean that greater SCZ PRS was associated with thinner rostral ACC in controls, but not in MDD cases (Fig. 2). Paired t tests confirmed that the case-control difference was significant at PRS thresholds 0.05 (t 111 = 2.590, p = 0.011) and 0.50 (t 111 = 2.289, p = 0.024) among those of below-average SCZ PRS, but not at any threshold for those of above-average SCZ PRS.
Mean thickness by SCZ PRS and MDD group for the two significant SCZ PRS by MDD interactions in rostral anterior cingulate cortex. SCZ, schizophrenia; PRS, polygenic risk score; MDD, major depressive disorder; RACC, rostral anterior cingulate cortex. Lines represent least-squares regression lines (red/darker grey = MDD cases; blue/lighter grey = controls); shaded areas represent 95% confidence intervals. Plots illustrate the two significant interactions between SCZ PRS at threshold 0.5 (left) or including all SNPs (right), and MDD status on RACC mean thickness.
Interactive effects between SCZ PRS and MDD status on mean cortical thickness by region
SCZ, schizophrenia; PRS, polygenic risk score; MDD, major depressive disorder; PFC, prefrontal cortex; MFG, middle frontal gyrus; IFG, inferior frontal gyrus; OFC, orbitofrontal cortex; ACC, anterior cingulate cortex; STG, superior temporal gyrus; MTG, middle temporal gyrus; ITG, inferior temporal gyrus.
Reported results are for interactive effects between SCZ PRS and MDD status on regional mean cortical thickness, controlling covariates.
*,** and *** represent significant results at p < 0.05, p < 0.01 and p < 0.001, respectively, before correcting for multiple comparisons. Significant results after FDR correction are highlighted in bold.
Given the specific effect of SCZ PRS on rostral ACC mean thickness among MDD cases, we tested whether it was influenced by a subgroup of participants. Repeated k-means clustering of SCZ PRS and rostral ACC thickness data using the function ‘pamk’ identified three as the optimum (according to average silhouette width; also verified using hierarchical clustering) for each of the three thresholds showing significant interactions. However, as shown in online Supplementary Fig. S1, for each threshold, none of the three clusters clearly represented a high-SCZ-risk subgroup, and none was significantly distinct from the other two, with low Dunn indices (DIs) of 0.009.
As more cortical data became available later, we were able to retest the interaction between SCZ PRS and MDD on rostral ACC thickness in a replication sample. The interaction was weaker but significant for SCZ PRS thresholds 0.05 (β = 0.095, p = 0.024) and 0.10 (β = 0.083, p = 0.047). As shown in online Supplementary Fig. S2, the interaction was attenuated mainly by the controls of the replication sample, among whom the negative association between SCZ PRS and rostral ACC thickness was much closer to zero than in the discovery sample (β = −0.013, p = 0.595; β = −0.009, p = 0.724). The association among MDD cases of the replication sample was still positive, but also weaker than in the discovery sample and no longer significant (β = 0.065, p = 0.074; β = 0.057, p = 0.116).
For cortical surface area, there were significant interactions across multiple PRS thresholds of β = −0.091 (p = 0.014) to β = −0.105 (p = 0.005) in inferior frontal gyrus (IFG), again due to a stronger negative effect of SCZ PRS among MDD cases than controls, but these results failed to survive FDR correction. Full results are reported in online Supplementary Table S1. For cortical volume, there were significant interactions across multiple PRS thresholds of β = −0.088 (p = 0.020) and β = −0.090 (p = 0.017) for IFG, β = 0.091 (p = 0.029) and β = 0.094 (p = 0.026) for rostral ACC; and β = 0.098 (p = 0.035) to β = 0.110 (p = 0.017) for caudal ACC. However, following FDR correction, none of these interactive effects on regional cortical volumes remained significant (online Supplementary Table S2).
Subcortical structures
We next tested interactive effects on subcortical volumes, which were available for all 4944 subjects. Across all subcortical volumes and PRS thresholds, there was only one significant SCZ PRS by MDD interaction: a weak positive interactive effect of β = 0.047 (p = 0.016) on thalamic volume at PRS threshold 0.01 (online Supplementary Table S3). This was driven by weak and non-significant but opposite associations between SCZ PRS and thalamic volume in MDD cases (β = 0.030) and controls (β = −0.022), but the interaction did not remain significant after FDR correction.
White matter tracts
Finally, we explored interactions in white matter tracts for the 4420 subjects with DTI data. There were no significant interactive effects between SCZ PRS and MDD on the FA of any of the white matter tracts at any of the PRS thresholds (online Supplementary Table S4). For MD, there was a significant interactive effect of β = 0.057 (p = 0.045) at threshold 0.50 and β = 0.061 (p = 0.033) with all SNPs in the uncinate fasciculus (online Supplementary Table S5). Uncinate MD also showed weak and non-significant opposite effects of SCZ PRS in MDD cases (β = 0.037; β = 0.035) and controls (β = −0.024; β = −0.030). No results survived correction for multiple comparisons.
Discussion
Following on from findings that MDD differs behaviourally as a function of genetic risk for SCZ, we tested interactions between SCZ PRS and MDD status on a range of brain structure measures in a large sample of UK Biobank participants. Interactions were not significant for most brain measures, but results for rostral ACC mean thickness at two PRS thresholds remained significant after correction for multiple comparisons. In rostral ACC, greater SCZ PRS showed a weak association with thinner cortex among controls, but a weak association with thicker cortex among MDD cases. However, MDD cases and controls only completely diverged toward the lower end of the range in SCZ PRS, where cases still showed thinner rostral ACC. We were able to replicate this interaction in a second subsample, although the effect was attenuated. This significant result seems to indicate localised reserve or protection against the negative neurostructural effects of SCZ-associated genes among individuals with MDD, but alternative explanations are discussed below.
We focused on interactions between SCZ PRS and MDD, but also tested main effects of SCZ PRS within MDD case and control groups for significant interactions. For controls, results suggested that a higher SCZ PRS was associated with reductions in rostral ACC thickness, as expected. These results are consistent with the reductions in prefrontal volumes seen in SCZ (Ross et al., Reference Ross, Margolis, Reading, Pletnikov and Coyle2006; Keshavan et al., Reference Keshavan, Tandon, Boutros and Nasrallah2008), and therefore with previous evidence of schizophrenia-like neuromorphology among high-risk but healthy controls (Lawrie et al., Reference Lawrie, Whalley, Abukmeil, Kestelman, Donnelly, Miller, Best, Owens and Johnstone2001; Boos et al., Reference Boos, Aleman, Cahn, Hulshoff Pol and Kahn2007). Further, case-control t tests clearly showed deficits in rostral ACC thickness among MDD cases with a low SCZ PRS, compared to low-SCZ-risk controls. This is also expected and consistent with previous evidence of reduced cingulate grey matter (Salvadore et al., Reference Salvadore, Nugent, Lemaitre, Luckenbaugh, Tinsley, Cannon, Neumeister, Zarate and Drevets2011; Rodriguez-Cano et al., Reference Rodríguez-Cano, Sarró, Monté, Maristany, Salvador, McKenna and Pomarol-Clotet2014; Schmaal et al., Reference Schmaal, Hibar, Sämann, Hall, Baune, Jahanshad, Cheung, van Erp, Bos, Ikram, Vernooij, Niessen, Tiemeier, Hofman, Wittfeld, Grabe, Janowitz, Bülow, Selonke, Völzke, Grotegerd, Dannlowski, Arolt, Opel, Heindel, Kugel, Hoehn, Czisch, Couvy-Duchesne, Rentería, Strike, Wright, Mills, de Zubicaray, McMahon, Medland, Martin, Gillespie, Goya-Maldonado, Gruber, Krämer, Hatton, Lagopoulos, Hickie, Frodl, Carballedo, Frey, van Velzen, Penninx, van Tol, van der Wee, Davey, Harrison, Mwangi, Cao, Soares, Veer, Walter, Schoepf, Zurowski, Konrad, Schramm, Normann, Schnell, Sacchet, Gotlib, MacQueen, Godlewska, Nickson, McIntosh, Papmeyer, Whalley, Hall, Sussmann, Li, Walter, Aftanas, Brack, Bokhan, Thompson and Veltman2017) in MDD. However, others have observed increased thickness of rostral ACC in MDD (Ancelin et al., Reference Ancelin, Carrière, Artero, Maller, Meslin, Ritchie, Ryan and Chaudieu2019) and of other cortical regions in SCZ (Dukart et al., Reference Dukart, Smieskova, Harrisberger, Lenz, Schmidt, Walter, Huber, Riecher-Rössler, Simon, Lang, Fusar-Poli and Borgwardt2017), which does correspond to our findings among individuals at higher genetic risk of SCZ.
Increased rostral ACC thickness among individuals with both MDD and a high genetic predisposition to SCZ may initially seem counterintuitive, suggesting that two putative neuropathologies instead confer benefits when combined. However, while positive associations between SCZ PRS and brain measures among individuals with MDD also stand in contrast with most previous findings (as above), some previous findings are comparable. For example, Papmeyer et al. (Reference Papmeyer, Giles, Sussmann, Kielty, Stewart, Lawrie, Whalley and McIntosh2015) reported thickening of inferior frontal and precentral cortices among individuals with both MDD and familial risk of bipolar disorder. One possible explanation is that the subset of MDD cases with a higher SCZ PRS in this study included individuals who had not developed SCZ despite a high level of genetic predisposition. These individuals may have had a neurobiological advantage over others at high genetic risk who did develop SCZ, who were not included in the study. Resistance to aberrant neurodevelopmental influences could explain both why high-SCZ-risk MDD cases did not develop SCZ and why they were more comparable to controls than to low-SCZ-risk cases. Or, simply, the greater rostral ACC thickness observed among high-SCZ-risk MDD cases may have protected against more severe psychopathology. However, as only a small proportion of individuals develop SCZ – even among those with a higher PRS – their exclusion is unlikely to have had such an impact on our results.
Alternatively, it may be that the increase in rostral ACC thickness that we observed does not actually represent a benefit. Increased thickness and functional connectivity of rostral ACC have been associated with insomnia (Winkelman et al., Reference Winkelman, Plante, Schoerning, Benson, Buxton, O'Connor, Jensen, Renshaw and Gonenc2013) and delusions in SCZ (Schott et al., Reference Schott, Voss, Wagner, Wüstenberg, Düzel and Behr2015), respectively. Such symptoms may relate to a deficit in, for example, normal axonal pruning, as this could lead to increased cortical thickness. Some theories of SCZ associate the disorder with altered pruning throughout development (Feinberg, Reference Feinberg1982; Keshavan et al., Reference Keshavan, Anderson and Pettegrew1994; Glausier and Lewis, Reference Glausier and Lewis2013), although they typically attribute it to excessive pruning, rather than decreased pruning. However, disruption of normal neurodevelopmental mechanisms could lead to excessive pruning in some areas but decreased pruning in others. While this could in turn lead to greater cortical thickness, these increases might be associated with aberrant hyperactivity or reduced efficiency, contributing to psychiatric symptoms. Further results do suggest SCZ-related increases in thickness and volume of specific brain regions, such as superior parietal cortex and the amygdala (Spoletini et al., Reference Spoletini, Piras, Fagioli, Rubino, Martinotti, Siracusano, Caltagirone and Spalletta2011; van Haren et al., Reference van Haren, Schnack, Cahn, van den Heuvel, Lepage and Collins2011), supporting this account. Increases in ACC thickness have also been observed with ageing in some studies (Salat et al., Reference Salat, Buckner, Snyder, Greve, Desikan, Busa, Morris, Dale and Fischl2004; Abe et al., Reference Abe, Yamasue, Aoki, Suga, Yamada, Kasai, Masutani, Kato, Kato and Ohtomo2008). Whether our results relate to a benefit or impairment specifically among high-SCZ-risk MDD cases, they provide some evidence of an SCZ-risk-related MDD subtype distinct in terms of prefrontal morphology.
Considering that MDD cases and controls only completely diverged at the lower end of the SCZ PRS scale, it may be that case-control differences in rostral ACC thickness were not reversed at the higher end of the range in SCZ PRS, but merely attenuated. In this case, the results relate more directly to Whalley et al.'s (Reference Whalley, Adams, Hall, Clarke, Fernandez-Pujals, Gibson, Wigmore, Hafferty, Hagenaars, Davies, Campbell, Hayward, Lawrie, Porteous, Deary and McIntosh2016) findings of attenuated case-control differences in neuroticism and psychological distress at higher SCZ PRS. This rekindles the idea that both sets of findings represent a protective effect of genetic risk for SCZ among MDD cases, and although it remains unclear why SCZ risk might be protective, the correspondence between clinical and neuroimaging findings strengthens evidence for a possible SCZ-risk-related MDD subtype. Although our cluster analyses were not able to distinguish this subtype, it warrants further investigation.
For most other structural brain metrics assessed, however, we found no significant influence of SCZ PRS on MDD case-control differences. This may have been due to limitations of the study. We used data from UK Biobank, which recruited large numbers of participants from the general population, but our sample (and particularly the smaller subset with data on cortical metrics) may not have been large enough to detect subtle genetic effects. The generalizability of results from UK Biobank is also still potentially limited by selection bias (Allen et al., Reference Allen, Sudlow, Downey, Peakman, Danesh, Elliot, Gallacher, Green, Matthews, Pell, Sprosen and Collins2012), as the sample is healthier overall than the general population (Fry et al., Reference Fry, Littlejohns, Sudlow, Doherty, Adamska, Sprosen, Collins and Allen2017). This is a common limitation to most observational studies but, as Fry et al. argue, results are likely to be both robust and widely generalizable. However, a limitation of recruiting from the general population is that SCZ PRS may not have been particularly high for any of our subjects. We still found significant heterogeneity in MDD related to SCZ PRS despite this, which perhaps lends further weight to the conclusion drawn from our significant results for rostral ACC thickness. Investigation of similar effects in a sample including a group of individuals with a higher PRS for SCZ (e.g. relatives of patients) may reveal more about the potential SCZ-risk-related MDD subtype.
UK Biobank participants are also older than the general population, with an average age of over 60 years, although a range of age groups from middle to late adulthood are represented. Focusing on this age group did, however, make it unlikely that subjects were yet to be diagnosed with SCZ, meaning that our results do relate to the effects of genetic risk only, rather than to SCZ itself. Another possible limitation was our measure of MDD, as it was derived from a web-based questionnaire. However, the online CIDI-SF was the most reliable method that could be efficiently used in such a large sample, and provided clinically relevant information based on DSM-IV diagnostic criteria. Finally, there were a number of potential covariates that we were unable to include, as complete data were unavailable at the time of analysis, for example, subjects' use of alcohol, tobacco, other drugs and medications, age of disorder onset and disorder severity. This should be considered when interpreting the current findings.
In summary, we tested for the presence of a distinct subtype of MDD associated with high genetic risk for SCZ by exploring interactive effects between these two variables on structural measures across the brain. Results were generally not significant, but we did observe a significant interactive effect on cortical thickness in rostral ACC for two PRS thresholds, which we were also able to replicate in a second subsample. While the direction of this interaction was unexpected, our results indicate that differences between MDD cases and controls in rostral ACC thickness are significantly influenced by SCZ PRS. This finding, together with related previous findings, suggest subtle differences in MDD related to SCZ risk and rostral ACC, which could influence important factors such as development of the disorder and response to treatment.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S003329171900165X.
Acknowledgements
This study was supported by a Wellcome Trust Strategic Award, ‘Stratifying Resilience and Depression Longitudinally' (STRADL; ref. 104036/Z/14/Z) and conducted using the UK Biobank Resource (application no. 4844). We would like to thank UK Biobank participants and researchers for their contributions to the study. Part of the work was undertaken at The Centre for Cognitive Ageing and Cognitive Epidemiology (CCACE), part of the cross-council Lifelong Health and Wellbeing Initiative (MR/K026992/1); funding from the Medical Research Council (MRC) and Biotechnology and Biological Sciences Research Council (BBSRC) is gratefully acknowledged. Age UK also provided support for the work undertaken at CCACE (via the Disconnected Mind project). We would also like to acknowledge David Liewald of CCACE for his part in processing imaging data. MAH is supported by The Dr Mortimer and Theresa Sackler Foundation. SRC is supported by the MRC (grants MR/M013111/1 and MR/R024065/1) and the National Institutes of Health (R01AG054628). HCW is supported by a John, Margaret, Alfred and Stewart (JMAS) Sim fellowship from the Royal College of Physicians of Edinburgh (RCPE), and an Edinburgh Scientific Academic Track (ESAT) fellowship from the University of Edinburgh.
Conflict of interest
In the past three years, SML has received research grant support from Janssen, Lundbeck and Sunovion, as well as personal fees from Otsuka, Sunovion and Janssen. AMM has previously received grant support from Pfizer, Lilly and Janssen. None of these funding sources are connected to the current study. MAH is funded by the Dr Mortimer and Theresa Sackler Foundation. Remaining authors have no biomedical financial interests or potential conflicts of interest to declare.
Open access
