Hostname: page-component-586b7cd67f-2plfb Total loading time: 0 Render date: 2024-11-25T22:24:08.991Z Has data issue: false hasContentIssue false

Depression with atypical neurovegetative symptoms shares genetic predisposition with immuno-metabolic traits and alcohol consumption

Published online by Cambridge University Press:  06 July 2020

Isabella Badini
Affiliation:
Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, UK
Jonathan R.I. Coleman
Affiliation:
Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, UK National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King's College London, London, UK
Saskia P. Hagenaars
Affiliation:
Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, UK
Matthew Hotopf
Affiliation:
National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King's College London, London, UK Department of Psychological Medicine, Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, UK South London and Maudsley NHS Foundation Trust, London, UK
Gerome Breen
Affiliation:
Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, UK National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King's College London, London, UK
Cathryn M. Lewis
Affiliation:
Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, UK National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King's College London, London, UK Department of Medical and Molecular Genetics, King's College London, London, UK
Chiara Fabbri*
Affiliation:
Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, UK Department of Biomedical and Neuromotor Sciences, University of Bologna, Bologna, Italy
*
Author for correspondence: Chiara Fabbri, E-mail: [email protected], [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Background

Depression is a highly prevalent and heterogeneous disorder. This study aims to determine whether depression with atypical features shows different heritability and different degree of overlap with polygenic risk for psychiatric and immuno-metabolic traits than other depression subgroups.

Methods

Data included 30 069 European ancestry individuals from the UK Biobank who met criteria for lifetime major depression. Participants reporting both weight gain and hypersomnia were classified as ↑WS depression (N = 1854) and the others as non-↑WS depression (N = 28 215). Cases with non-↑WS depression were further classified as ↓WS depression (i.e. weight loss and insomnia; N = 10 142). Polygenic risk scores (PRS) for 22 traits were generated using genome-wide summary statistics (Bonferroni corrected p = 2.1 × 10−4). Single-nucleotide polymorphism (SNP)-based heritability of depression subgroups was estimated.

Results

↑WS depression had a higher polygenic risk for BMI [OR = 1.20 (1.15–1.26), p = 2.37 × 10−14] and C-reactive protein [OR = 1.11 (1.06–1.17), p = 8.86 × 10−06] v. non-↑WS depression and ↓WS depression. Leptin PRS was close to the significance threshold (p = 2.99 × 10−04), but the effect disappeared when considering GWAS summary statistics of leptin adjusted for BMI. PRS for daily alcohol use was inversely associated with ↑WS depression [OR = 0.88 (0.83–0.93), p = 1.04 × 10−05] v. non-↑WS depression. SNP-based heritability was not significantly different between ↑WS depression and ↓WS depression (14.3% and 12.2%, respectively).

Conclusions

↑WS depression shows evidence of distinct genetic predisposition to immune-metabolic traits and alcohol consumption. These genetic signals suggest that biological targets including immune-cardio-metabolic pathways may be relevant to therapies in individuals with ↑WS depression.

Type
Original Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s), 2020. Published by Cambridge University Press

Introduction

Depressive disorders are highly prevalent and a leading cause of global disability and are associated with premature mortality (James et al., Reference James, Abate, Abate, Abay, Abbafati, Abbasi and Murray2018; Kessler & Bromet, Reference Kessler and Bromet2013). Twin-based heritability of unipolar depression is estimated to be ~37% (Sullivan, Neale, & Kendler, Reference Sullivan, Neale and Kendler2000), with common single-nucleotide polymorphisms (SNPs) explaining ~9% of variation in depression liability (Howard et al., Reference Howard, Adams, Clarke, Hafferty, Gibson, Shirali and McIntosh2019; Wray et al., Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui2018). Efforts to identify genetic variants associated with the disease are hindered by the heterogeneity among depressed cases, who can vary greatly in symptom presentation and severity, clinical course and treatment response (Fried & Nesse, Reference Fried and Nesse2015). Clinical heterogeneity may also reflect different underlying biological and causal pathways. Increasing evidence suggests that depressive symptoms and subtypes are differentially associated with genetic risk factors that overlap with other disorders (Beijers, Wardenaar, van Loo, & Schoevers, Reference Beijers, Wardenaar, van Loo and Schoevers2019; Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Abdellaoui, Willemsen, Hottenga and Penninx2016, Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017). Investigating the clinical and genetic correlates of more homogenous subtypes may improve the understanding of specific aetiological mechanisms and the development of potential treatment targets. Similarly, identifying whether genetic risk for other disorders overlaps with certain characteristics of depression may help to elucidate biological mechanisms underlying common symptom presentations across multiple disorders.

The atypical subtype of major depression is specified in the DSM-5 by the presence of at least two of the following symptoms: hypersomnia, increased appetite and/or weight gain, leaden paralysis, and interpersonal rejection sensitivity (American Psychiatric Association, 2013). Reversed neurovegetative symptoms (hypersomnia, increased appetite and/or weight gain), in particular, have been found to be highly specific and predictive of clinically defined atypical MDD, as they identify patients with similar sociodemographic and clinical correlates as those reported for classification based on the full DSM criteria (Benazzi, Reference Benazzi2002). Classification of atypical depression based on reversed neurovegetative symptoms alone is more feasible in epidemiological studies, many of which have adopted this criterion instead of the full DSM-5 criteria (Lee, Ng, & Tsang, Reference Lee, Ng and Tsang2009; Matza, Revicki, Davidson, & Stewart, Reference Matza, Revicki, Davidson and Stewart2003).

Epidemiological studies have shown differences between atypical and non-atypical depression in sociodemographic factors, clinical features, lifestyle factors and comorbidities. Atypical depression has been associated with earlier age of onset, female gender, more severe and recurrent depressive episodes (Agosti & Stewart, Reference Agosti and Stewart2001; Blanco et al., Reference Blanco, Vesga-López, Stewart, Liu, Grant and Hasin2012; Brailean, Curtis, Davis, Dregan, & Hotopf, Reference Brailean, Curtis, Davis, Dregan and Hotopf2020). Findings from UK Biobank (UKB) showed higher rates of smoking, social isolation, loneliness, greater exposure to adverse life events and lower rates of moderate physical activity among atypical cases compared to non-atypical cases (Brailean et al., Reference Brailean, Curtis, Davis, Dregan and Hotopf2020). Atypical depression has also been associated with higher rates of bipolar disorder and psychiatric comorbidity such as anxiety disorders, binge eating disorder and substance abuse (Agosti & Stewart, Reference Agosti and Stewart2001; Blanco et al., Reference Blanco, Vesga-López, Stewart, Liu, Grant and Hasin2012; Brailean et al., Reference Brailean, Curtis, Davis, Dregan and Hotopf2020; Lee et al., Reference Lee, Ng and Tsang2009; Łojko, Buzuk, Owecki, Ruchała, & Rybakowski, Reference Łojko, Buzuk, Owecki, Ruchała and Rybakowski2015). Physical health comorbidities more strongly associated with atypical cases include higher body mass index (BMI), inflammation, metabolic syndrome and cardiovascular disease (CAD) (Brailean et al., Reference Brailean, Curtis, Davis, Dregan and Hotopf2020; Lasserre et al., Reference Lasserre, Glaus, Vandeleur, Marques-Vidal, Vaucher, Bastardot and Preisig2014; Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017). Specifically, evidence suggests stronger links between atypical features of increased appetite and/or weight and immuno-metabolic dysregulations, such as BMI, C-reactive protein (CRP) and leptin (Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017).

There is increasing evidence to suggest partially distinct genetic profiles among depressive subtypes. SNP-based heritability (h 2SNP) was found to vary across individual depressive symptoms (h 2SNP range from 6% to 9%), and patterns of SNP associations and genetic correlations differed across symptoms (Thorp et al., Reference Thorp, Marees, Ong, An, MacGregor and Derks2019). A study that classified cases into subtypes according to change in neurovegetative symptoms [i.e. no change, increased or decreased appetite and/or weight (A/W)] found similar SNP-heritability of 10–11% in all groups (Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017). Polygenic risk scoring analyses confirmed that the increased A/W subtype had higher polygenic risk for CRP, BMI, leptin and triglycerides levels than the decreased A/W subtype (Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Abdellaoui, Willemsen, Hottenga and Penninx2016, Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017). Overlapping genetic aetiology between depression and CAD was demonstrated (Hagenaars et al., Reference Hagenaars, Coleman, Choi, Gaspar, Adams, Howard and Lewis2019), but no studies investigated if depression with atypical features may have a greater genetic overlap with CAD compared to other depression subtypes.

In the present study, we examined the genetic overlap between depression with atypical features (increased weight and sleepiness: ↑WS) and a range of traits and disorders using polygenic risk scores (PRS). Based on previous findings, we hypothesised that ↑WS depression would show similar heritability to depression without ↑WS and higher polygenic risk for immune-cardio-metabolic, substance use and other psychiatric traits compared to depression without ↑WS.

Materials and methods

Sample

Individuals who met lifetime criteria for major depression were drawn from the UKB. UKB is a prospective population-based study of ~500 000 individuals recruited across the UK, aged between 40 and 69 at baseline (UK Biobank, 2019). A total of 157 387 participants completed an online Mental Health Questionnaire (MHQ) assessing self-reported psychiatric symptoms corresponding to clinical diagnostic criteria and self-reported professional diagnoses (Davis et al., Reference Davis, Coleman, Adams, Allen, Breen, Cullen and Hotopf2020). Genome-wide genetic data have been collected on all UKB participants (Bycroft et al., Reference Bycroft, Freeman, Petkova, Band, Elliott, Sharp and Marchini2018), as detailed below. All participants provided written informed consent and all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.

Measures

Lifetime psychiatric diagnoses were assessed in the MHQ using the Composite International Diagnostic Interview Short Form (CIDI-SF) (Kessler, Andrews, Mroczek, Ustun, & Wittchen, Reference Kessler, Andrews, Mroczek, Ustun and Wittchen1998). Criteria for lifetime major depressive episode were in accordance with DSM-V. The full CIDI is a validated measure of depression, demonstrated to have good concordance with direct clinical assessment (Haro et al., Reference Haro, Arbabzadeh-Bouchez, Brugha, De Girolamo, Guyer, Jin and Kessler2006).

Cases reporting both weight gain (data field 20 536) and hypersomnia (data field 20 534) were classified as ↑WS depression (N = 1854), and the remaining cases classified as depression without ↑WS (N = 28 215). From this group, depression with both weight loss and decreased sleep (data fields 20 533 and 20 535) was defined as ↓WS depression (N = 10 142). These definitions used the same coding conventions as those used in a previous study in the UKB (Brailean et al., Reference Brailean, Curtis, Davis, Dregan and Hotopf2020); however, we preferred to avoid the terms atypical and typical depression to avoid confusion with the standard nosological classification. We did not consider variations in appetite to distinguish the subgroups, since the available measure (data field 20 511) did not differentiate between hypophagia and hyperphagia.

Genotyping and quality control

Genetic data came from the full release of the UKB data (N = 488 377; Bycroft et al., Reference Bycroft, Freeman, Petkova, Band, Elliott, Sharp and Marchini2018). Genotyping was performed using two highly-overlapping arrays covering ~800 000 markers (UK Biobank Axiom Array Content Summary). Autosomal genotype data underwent centralised quality control to adjust for possible array effects, batch effects, plate effects and departures from Hardy–Weinberg equilibrium (HWE; Bycroft et al., Reference Bycroft, Freeman, Petkova, Band, Elliott, Sharp and Marchini2018). Variants for this analysis were limited to common variants (minor allele frequency >0.01) that were directly genotyped. SNPs were further excluded based on missingness (>0.02) and on HWE (p < 10–8). Individuals were removed for high levels of missingness (>0.05) or abnormal heterozygosity (as defined during centralised quality control), relatedness of up to third-degree kinship (KING r < 0.044; Manichaikul et al., Reference Manichaikul, Mychaleckyj, Rich, Daly, Sale and Chen2010), or phenotypic and genotypic gender discordance (phenotypic males with F X < 0.9, phenotypic females with F X > 0.6). Population structure within the UKB cohort was assessed using principal component analysis, with European ancestry defined by four-means clustering on the first two genetic principal components (Warren et al., Reference Warren, Evangelou, Cabrera, Gao, Ren and Mifsud2017). Among respondents to the MHQ, 95% were of European ancestry and therefore individuals from other ancestries were excluded from further analyses to maximise statistical power. After quality control, the final sample of respondents to the MHQ consisted of 126 522 individuals with genotype data.

Statistical analysis

Polygenic risk scores

PRS were calculated based on GWAS summary statistics for 22 traits reflecting the hypothesis formulated in the Introduction, including major psychiatric disorders, personality traits, substance use-related traits, cardio-metabolic traits and CRP (online Supplementary Table S1). There was no overlap between the samples included in these GWASs and the sample included in this study, except a very marginal overlap with the GWAS of anorexia nervosa [349/16 992 (2%) of cases included in Watson et al. (Reference Watson, Yilmaz, Thornton, Hübel, Coleman, Gaspar and Bulik2019)].

PRS were calculated using PRSice v.2 (Choi & O'Reilly, Reference Choi and O'Reilly2019; Euesden, Lewis, & O'Reilly, Reference Euesden, Lewis and O'Reilly2015). PRSice computes scores in an independent (target) sample by calculating the weighted sum of trait-associated alleles using summary data from GWAS discovery samples. SNPs in linkage disequilibrium [r 2 ⩾ 0.1 (250-kb window)] were removed using the clumping procedure. We used the default average option that calculates the ratio between the PRS and the number of alleles included in each individual and PRS were standardised (mean = 0, s.d. = 1). PRS were calculated at 11 p value thresholds P T (5 × 10−8, 1 × 10−5, 1 × 10−3, 0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 1) and the most predictive P T was selected. Logistic regression models were used to estimate associations between ↑WS depression v. non-↑WS depression and each PRS adjusting for covariates of six genetic ancestry principal components, centre and batch effects. We estimated the proportion of variance explained by PRS on the observed and liability scale (Lee, Goddard, Wray, & Visscher, Reference Lee, Goddard, Wray and Visscher2012), considering a range of possible values of prevalence among depressed cases (Levitan, Lesage, Parikh, Goering, & Kennedy, Reference Levitan, Lesage, Parikh, Goering and Kennedy1997; Łojko & Rybakowski, Reference Łojko and Rybakowski2017). For PRS of traits associated with ↑WS depression, we calculated: (1) the OR of ↑WS depression v. depression without ↑WS for each decile of the PRS, taking the first decile as reference and (2) if the effect was comparable when considering each of the symptoms separately (↑weight and ↑ sleep v. no increase in weight or sleep). We estimated differences between PRS results of different comparisons by comparing their estimates (E1 and E2) and SE (SE1 and SE2) using a Z-test $\left({Z = ( {E1-E2} ) /\sqrt {{( {SE1} ) }^2 + {( SE2) }^2} } \right)$. The code used for these analyses is available as online Supplementary material (code_used_for_analyses).

A Bonferroni correction was applied to account for multiple testing, providing a required significance level of p < 2.1 × 10−4 [the PRS of 17 traits were analysed at 11 P T, while for five traits there were no SNPs with p < 5 × 10−8 and 10 thresholds were tested (0.05/(11 × 17 + 10 × 5) = 2.1 × 10−4)]. The use of Bonferroni correction is conservative, since the different P T are highly correlated.

SNP-based heritability

Using genome-wide complex trait analysis software v.1.93.1beta (GCTA) (Yang, Lee, Goddard, & Visscher, Reference Yang, Lee, Goddard and Visscher2011), genetic relationship matrix-restricted maximum likelihood (GREML) methods (Lee, Wray, Goddard, & Visscher, Reference Lee, Wray, Goddard and Visscher2011) were used to estimate the variance in liability attributable to the additive effects of all SNPs (h 2SNP) for ↑WS depression, ↓WS depression and depression not falling in these two groups. A random subset of 8000 healthy controls were selected for the analyses among those who completed the MHQ. This number of controls was selected because it provided adequate power (at least 80%) to estimate heritability considering an expected heritability of 0.10 (Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017), using GCTA-GREML Power Calculator (Visscher et al., Reference Visscher, Hemani, Vinkhuyzen, Chen, Lee, Wray and Yang2014). The genetic relationship matrix (GRM) was adjusted for incomplete tagging of causal SNPs and we excluded related individuals using a GRM-cut-off of 0.05. We calculated the genetic correlation between depression subgroups using bivariate GREML with independent subsets of 8000 healthy controls for each depression subgroup.

We also calculated h 2SNP using GCTB (Genome-wide Complex Trait Bayesian analysis) Bayes S method, which estimates polygenicity from the data (i.e. the proportion of SNPs with non-zero effects, π). GCTB also calculates the relationship between effect size and MAF (S) which can be used to detect signatures of natural selection (Zeng et al., Reference Zeng, de Vlaming, Wu, Robinson, Lloyd-Jones, Yengo and Yang2018). We used the standard settings of 21 000 simulations with the first 1000 as burn-in and standard initial S and π values.

GCTA-GREML and GCTB Bayes S analyses were adjusted for the same covariates included in the PRS analysis, but we also adjusted for BMI because there was a significant difference in BMI between cases with ↑WS depression (30.55 ± 5.78), ↓WS depression (25.93 ± 4.41) and healthy controls (26.45 ± 4.14) (p = 1.51 × 10−221 and p = 1.25 × 10−154, respectively); heritability estimates would be confounded by this variable which was demonstrated to have a heritability of about 28% (Zeng et al., Reference Zeng, de Vlaming, Wu, Robinson, Lloyd-Jones, Yengo and Yang2018).

Possible differences among h 2SNP of depression subgroups were compared by a Z-test as explained for PRS results. h 2SNP estimates were transformed to a liability scale (Lee et al., Reference Lee, Goddard, Wray and Visscher2012) and we reported heritability considering a range of plausible prevalence values (Levitan et al., Reference Levitan, Lesage, Parikh, Goering and Kennedy1997; Lim et al., Reference Lim, Tam, Lu, Ho, Zhang and Ho2018; Łojko & Rybakowski, Reference Łojko and Rybakowski2017). The code used for these analyses is available as online Supplementary material (code_used_for_analyses).

Sensitivity analyses

Analyses were repeated (1) excluding cases with probable bipolar disorder or missing information for this variable [wider bipolar disorder definition, as described in Brailean et al. (Reference Brailean, Curtis, Davis, Dregan and Hotopf2020), n = 1747] and schizophrenia or missing information for this variable (n = 61); (2) comparing ↑WS depression with ↓WS depression (rather than non-↑WS depression); and (3) comparing both ↑WS depression and ↓WS depression with healthy controls who completed the MHQ (n = 64 604) (Fig. 1).

Fig. 1. Flow chart of the performed analyses. ↑WS, depression with increased weight and sleepiness; ↓WS depression, depression with decreased weight and insomnia.

Results

Sample characteristics

The total sample comprised of 30 069 participants who met criteria for lifetime major depressive episode, including 1854 cases classified as ↑WS depression and 28 215 classified as depression without ↑WS. Non-↑WS cases were further classified as ↓WS depression (N = 10 142). The ↑WS group was 75% female and had a mean age of 59.99 (s.d. = 7.09) years, compared to the non-↑WS group which was 68% female and had a mean age of 62.50 (7.53) years. The ↓WS depression group was 75% female and had a mean age of 62.55 (7.49). Further description of the clinical-demographic features of the subtypes can be found in a previous paper in UKB that used the same classification criteria (Brailean et al., Reference Brailean, Curtis, Davis, Dregan and Hotopf2020).

Polygenic risk analysis

Traits where the PRS were significantly associated with ↑WS depression compared to non-↑WS depression included BMI [OR = 1.20 (1.15–1.26), p = 2.37 × 10−14], CRP [OR = 1.11 (1.06–1.17), p = 8.86 × 10−06], daily alcohol use [OR = 0.88 (0.83–0.93), p = 1.04 × 10−05] and MDD [OR = 1.10 (1.05–1.15), p = 1.14 × 10−04], as shown in Table 1 (see online Supplementary Table S2 for all results). The OR for ↑WS depression was 2.03 (1.62–2.55), 1.36 (1.10–1.68), 0.62 (0.49–0.79) and 1.30 (1.04–1.61) in the 10th PRS decile v. the 1st PRS decile for BMI, CRP, daily alcohol use and MDD, respectively (Fig. 2). Results were consistent between different P T (online Supplementary Fig. S1). The same direction of effect was observed when comparing cases with weight gain v. cases without weight gain, with similar effect size for alcohol daily use, MDD, leptin and CRP PRS, while larger effect size for BMI PRS (z = 2.12, p = 0.03). These PRS showed the same direction of effect though not significant when considering depression with hypersomnia compared to depression without this symptom (online Supplementary Table S2C). The results were consistent when comparing ↑WS depression with ↓WS depression and ↑WS depression with healthy controls (online Supplementary Tables S3 and S4). Interestingly, the effect of BMI PRS was very close to the significance threshold for being inversely associated with ↓WS depression v. healthy controls (p = 1.24 × 10−04) while leptin, CRP and alcohol daily use PRS had no effect (online Supplementary Table S5).

Fig. 2. Decile plots showing the odds ratio of ↑WS depression (depression with atypical features) v. depression without ↑WS for BMI PRS (a), C-reactive protein (CRP) PRS (b), daily alcohol use PRS (c) and major depressive disorder (MDD) PRS (d). *p < 0.05; **p < 1 × 10−05.

Table 1. Polygenic risk scores associated with depression with atypical features (↑WS depression) compared with depression without atypical features

Results are shown for the P T threshold attaining the lowest p value.

R 2 estimated on the observed scale as the results are referred to case-only comparisons.

*Significant associations (p < 2.1 × 10−4).

The association between ↑WS v. depression without ↑WS and leptin PRS was close to the significance threshold [OR = 1.09 (1.04–1.14), p = 2.99 × 10−04] but disappeared when considering GWAS summary statistics of leptin adjusted for BMI. Nominal significant associations (p < 0.05) included PRS for type 2 diabetes, coronary artery disease, triglycerides, total HDL cholesterol and ischaemic stroke (Fig. 3; online Supplementary Table S2). The PRS for coronary artery disease reached statistical significance for association with both ↑WS depression and ↓WS depression when the comparator group was healthy controls, in addition to triglycerides PRS only for ↑WS depression v. controls. Case–healthy controls analyses showed expected associations with the PRS of psychiatric traits and contributed to clarifying the effect of alcohol use-related traits PRS: higher PRS for alcohol dependence but not daily alcohol use increased the risk of ↓WS depression, while ↑WS depression was still associated with lower daily alcohol use PRS (online Supplementary Tables S4–S5, Fig. 4). Sensitivity analysis revealed very similar results when excluding cases with probable bipolar disorder and schizophrenia or missing information for these variables (n = 196 and n = 12 among patients with ↑WS depression; n = 1551 and n = 49 among those with non-↑WS depression; online Supplementary Figs S2–S4).

Fig. 3. PRS odds ratio (OR) and 95% confidence intervals of ↑WS depression v. depression without ↑WS. Colors indicate different groups of traits (substance related disorders, major psychiatricdisorders, personality traits, immune metabolic traits). BP, bipolar disorder; SCZ, schizophrenia; ANX, anxiety disorders; PTSD, post-traumatic stress disorder; AN, anorexia nervosa; ALCDEP, alcohol dependence; ALCUSE, daily alcohol use; N_CIGARETTES, n cigarettes per day; CANN, cannabis use lifetime; EXTR, extraversion; NEU, neuroticism; DM2, type 2 diabetes mellitus; CAD, coronary artery disease; ISCH_STROKE, ischaemic stroke; TG, triglycerides; LDL, total LDL cholesterol; HDL, total HDL cholesterol; BMI, body max index; BMI_adjust_leptin, leptin adjusted for BMI; CRP, C-reactive protein.

Fig. 4. PRS odds ratio (OR) and 95% confidence intervals of ↑WS depression and ↓WS depression compared with healthy controls. BP, bipolar disorder; SCZ, schizophrenia; ANX, anxiety disorders; PTSD, posttraumatic stress disorder; AN, anorexia nervosa; ALCDEP, alcohol dependence; ALCUSE, daily alcohol use; N_CIGARETTES, n cigarettes per day; CANN, cannabis use lifetime; EXTR, extraversion; NEU, neuroticism; DM2, type 2 diabetes mellitus; CAD, coronary artery disease; ISCH_STROKE, ischaemic stroke; TG, triglycerides; LDL, total LDL cholesterol; HDL, total HDL cholesterol; BMI, body max index; BMI_adjust_leptin, leptin adjusted for BMI; CRP, C-reactive protein.

SNP-based heritability

SNP-based heritability (h 2SNP) on the liability scale for ↑WS depression was estimated to be 0.14 (95% CI 0.03–0.25) using GCTB-Bayes S and 0.21 (95% CI 0.07–0.35) using GCTA-GREML (Table 2). h 2SNP was 0.12 (95% CI 0.08–0.16) and 0.16 (95% CI 0.10–0.22) for ↓WS depression using the two methods, respectively. GCTA-GREML h 2SNP was not significantly different comparing ↑WS depression with ↓WS depression (z = 0.57, p = 0.56), depression with weight gain (z = 1.09, p = 0.27), depression with weight loss (z = 0.75, p = 0.46) or depression without ↑WS or ↓WS (z = 1.25, p = 0.21). GCTB-Bayes S estimates were not significantly different compared with GCTA-GREM estimates for any depression subgroup (Table 2). The estimated proportion of SNPs contributing to h 2SNP ranged from 2% to 4%; S estimation was negative for all traits but not significantly different from zero, except for depression with weight gain (z = 2.36, p = 0.018), which however did not survive Bonferroni correction for five tests (online Supplementary Fig. S5).

Table 2. SNP heritability (SNP-h2) of depression with weight gain and hypersomnia (↑WS depression), depression with weight loss and reduced sleep (↓WS depression), major depression not falling in any of these two groups (no ↑WS or ↓WS), depression with weight gain or weight loss

Heritability and s.e. in parenthesis are reported on the liability scale according to different possible values of prevalence (K) of each subtype in the population. The most credible values of population prevalence are in bold and a Z-test was performed to compare the estimate and s.e. of GCTA-GREML and GCTB Bayes S results.

Genetic correlations were 0.54 (s.e. = 0.14) for ↑WS depression and ↓WS depression, 1.05 (s.e. = 0.15) for ↑WS depression and depression without ↑WS or ↓WS, and 0.74 (s.e. = 0.14) between ↓WS depression and depression without ↑WS or ↓WS.

Discussion

This study examined the genetic overlap between depression with atypical features (↑WS) and a range of traits and disorders using PRS in 30 069 cases with major depression from the UKB. The findings showed that persons with higher BMI PRS, CRP PRS and MDD PRS are more likely to have ↑WS depression v. non-↑WS depression and ↓WS depression, while those with lower PRS for alcohol daily use were more likely to have non-↑WS depression rather than ↑WS depression. Associations with these PRS were consistent among different P T. The effect of BMI PRS and CRP PRS on the risk of ↑WS depression was similar when taking healthy controls as comparator group, while it was in the opposite direction or absent, respectively, when we compared ↓WS depression v. healthy controls. The analyses of individual symptoms (weight changes and sleep changes) showed that BMI, non-BMI adjusted leptin and CRP PRS had a higher effect size on depression with weight increase than sleep increase, but not the PRS of alcohol daily use and MDD; however, the direction of the effect on weight and sleep increase was the same for all these PRS, suggesting that the symptoms of weight gain and sleepiness in depression have shared rather than divergent genetics. The PRS of coronary artery disease and triglycerides was associated with ↑WS depression v. healthy controls, in line with nominal associations between the PRS of other cardio-metabolic traits (e.g. type 2 diabetes) and the risk of ↑WS depression compared to other depression subgroups as well as healthy controls. The PRS for alcohol dependence was associated with the risk of ↓WS depression compared to healthy controls but not for ↑WS depression compared to healthy controls, though the direction of the effect was the same. Interestingly, the PRS of alcohol daily use was similar between ↓WS depression and healthy controls, but significantly lower in ↑WS depression, suggesting that the pathogenesis of alcohol use disorders (AUD) may involve different mechanisms in ↑WS compared to ↓WS depression. A recent GWAS demonstrated that alcohol consumption and AUD show significant genetic differences, with the genetics of AUD being more closely related to other psychiatric disorders, and the genetics of alcohol consumption to that of some positive health outcomes, such as reduced risk of CAD, and lower BMI, in line with our findings (Kranzler et al., Reference Kranzler, Zhou, Kember, Vickers Smith, Justice, Damrauer and Gelernter2019).

Taken together, these results suggest partially distinct genetic pathways between depression with atypical and typical neurovegetative symptoms and that this divergence may be attributable to distinct genetic predisposition to immune-metabolic traits. Although there is no convincing evidence that these depression subgroups may respond differently to conventional antidepressant treatments, drugs acting on the specific biological mechanisms implicated in ↑WS depression may have clinical benefits. For example, peroxisome proliferator-activated receptor (PPAR)-γ agonists target insulin resistance and the related oxidative and pro-inflammatory changes (which are also involved in the pathogenesis of depressive symptoms); they were demonstrated to have antidepressant effects in patients with treatment-resistant bipolar depression and concomitant insulin resistance (Kemp et al., Reference Kemp, Schinagle, Gao, Conroy, Ganocy, Ismail-Beigi and Calabrese2014).

The SNP-based heritability was estimated to be similar between ↑WS depression and other subtypes and it did not significantly change after excluding cases with probable bipolar disorder and schizophrenia. GCTA-GREML h 2SNP was similar to a previous study of depression with weight/appetite gain [h 2SNP = 0.11 (s.e. = 0.03)] and weight/appetite loss [h 2SNP = 0.10 (s.e. = 0.02)] (z = 1.24, p = 0.22; z = 1.72, p = 0.09, respectively) (Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017). However, this study considered both weight and appetite changes in the definition of the subgroups, while we used only weight changes because of the lack of information on the direction of appetite changes in UKB. Genetic correlations suggested that ↑WS depression genetics is highly shared with depression without ↑WS or ↓WS, and the genetic correlation between these subgroups was significantly higher than the genetic correlation between ↑WS depression and ↓WS depression (z = 2.49, p = 0.013). In line with previous evidence, major depression subgroups did not show evidence of negative natural selection based on GCTB-Bayes S estimates (Zeng et al., Reference Zeng, de Vlaming, Wu, Robinson, Lloyd-Jones, Yengo and Yang2018).

Observational studies have shown differential pathophysiological correlates among subtypes of depression, with atypical depression more strongly associated with obesity and metabolic dysregulations than other depression groups (Brailean et al., Reference Brailean, Curtis, Davis, Dregan and Hotopf2020; Lasserre et al., Reference Lasserre, Glaus, Vandeleur, Marques-Vidal, Vaucher, Bastardot and Preisig2014; Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017). Consistent with previous research, the present findings suggest that the phenotypic association between depression with atypical features and cardio-metabolic traits may partly result from shared genetic and biological mechanisms (Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Abdellaoui, Willemsen, Hottenga and Penninx2016, Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017). For example, pleiotropy may occur when shared genetic variants influence obesity-related traits and ↑WS depression through common pathways such as inflammation and/or leptin system dysregulation (Guo & Lu, Reference Guo and Lu2014; Ring & Zeltser, Reference Ring and Zeltser2010). A previous study (Milaneschi et al., Reference Milaneschi, Lamers, Peyrot, Baune, Breen and Dehghan2017) reported that the genetic correlation between leptin and ↑WS depression was decreased but not absent when considering leptin adjusted for BMI, while we found no association between PRS for BMI-adjusted leptin levels and ↑WS depression or depression with weight increase. Leptin reduces food intake and obesity is associated with increased leptin levels, through the induction of leptin resistance (Myers, Leibel, Seeley, & Schwartz, Reference Myers, Leibel, Seeley and Schwartz2010). The adjustment of leptin levels for BMI eliminates an important source in leptin inter-individual variability, and separates the genetic factors regulating BMI from those that regulate only leptin production. Our results indicate that the genetic factors predisposing to increased BMI and increased leptin levels in ↑WS depression are highly correlated, in line with the hypothesis that increased leptin levels are generally a consequence of increased BMI rather than a cause (at least in common forms of obesity), although different pathogenetic mechanisms may lead to similar phenotypes (Myers et al., Reference Myers, Leibel, Seeley and Schwartz2010). This discrepancy may also relate to the different criteria used to define depression subgroups and sample-specific characteristics.

The current findings also showed that PRS for alcohol daily use was inversely associated with risk for ↑WS depression both in case only and case–control analyses, while it showed no difference between ↓WS depression and healthy controls. In case–control comparisons, the PRS of alcohol dependence was associated with ↓WS depression and at the nominal level with ↑WS depression. Epidemiological studies have shown increased alcohol consumption and/or risk of AUD among persons with atypical depression (Blanco et al., Reference Blanco, Vesga-López, Stewart, Liu, Grant and Hasin2012; Brailean et al., Reference Brailean, Curtis, Davis, Dregan and Hotopf2020). There are no previous studies which examined the genetic overlap between alcohol use and atypical depression, although a positive genetic correlation between MDD and alcohol dependence has been reported (Andersen et al., Reference Andersen, Pietrzak, Kranzler, Ma, Zhou, Liu and Han2017; Kranzler et al., Reference Kranzler, Zhou, Kember, Vickers Smith, Justice, Damrauer and Gelernter2019). Our results suggest that genetic risk for daily alcohol use is lower in ↑WS depression v. other groups, in contrast with epidemiological observations of higher risk of AUD in atypical cases. However, as previously noted, the genetics of alcohol daily use only partially overlaps with the genetics of alcohol dependence. In addition, the increased rate of alcohol use among cases with ↑WS depression may result from secondary or environmental risk factors rather than from a genetic aetiology. In line with this hypothesis, depression with atypical features in UKB show longer, more severe and recurrent disease episodes, increased risk of comorbidities, and more lifetime deprivation and adversity (Brailean et al., Reference Brailean, Curtis, Davis, Dregan and Hotopf2020), and these factors may be responsible for increased risk of AUD (Yang et al., Reference Yang, Tao, He, Liu, Wang and Zhang2018).

Several limitations should be considered when interpreting the results from the present study. First, the UKB is not representative of the general population, with respondents more likely to be older, female, healthier, of a higher socioeconomic background and better educated (Fry et al., Reference Fry, Littlejohns, Sudlow, Doherty, Adamska, Sprosen and Allen2017). Compared to the overall UKB cohort, respondents of the MHQ have on average higher educational and occupational attainment, are less likely to smoke or report longstanding illness/disability (Davis & Hotopf, Reference Davis and Hotopf2019). This is likely to impact on prevalence and severity of depression within the sample, with persons with severe depression less likely to have completed the MHQ. Second, measures for depression were based on recall of the worst episode of depression and assessed using self-reported symptoms, rather than clinical diagnoses. This raises the possibility that self-reported symptoms may be affected by recall bias or other medical conditions, and we are not able to exclude that participants had other depressive episodes with different neurovegetative symptoms compared to the one reported in the MHQ. Third, our measurement of depression with atypical features did not capture all DSM-5 symptoms but only the neurovegetative component. Our results were in line with previous studies and contributed to characterise the genetics of these symptoms, but future research using the full atypical spectrum should be considered and dimensional classifications should be considered as an alternative to binary categorisation. Fourth, due to the small numbers of individuals of non-European groups in UKB, and their systematic under-representation in GWAS studies used to derive PRS, our sample is limited to individuals of European descent, limiting the generalisability of our findings. Last, we had no statistical power to perform a genome-wide association study of ↑WS depression or Mendelian randomisation to assess the causal relationship between immune-metabolic traits and ↑WS depression.

In conclusion, this study showed specific genetic overlap of ↑WS and ↓WS depression with immune-metabolic and alcohol-related traits. Our findings suggest that depression with typical and atypical neurovegetative symptoms may represent relatively homogenous subtypes characterised by partially distinct genetic liabilities. Understanding the shared genetic aetiology between immune-metabolic traits and depression subtypes could play an important role in the prevention/treatment of depressive episodes and the development of tailored treatments.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291720002342.

Acknowledgements

The authors acknowledge the use of the research computing facility at King's College London, Rosalind (https://rosalind.kcl.ac.uk), which is delivered in partnership with the National Institute for Health Research (NIHR) Biomedical Research Centres at South London & Maudsley and Guy's & St. Thomas' NHS Foundation Trusts, and part-funded by capital equipment grants from the Maudsley Charity (award 980) and Guy's & St. Thomas' Charity (TR130505). UK Biobank data were accessed under application 18177 ‘Multi-trait GWAS analyses in the UK Biobank’. We thank the CHARGE Inflammation Working Group (CIWG), the Psychiatric Genomics Consortium (PGC), the International Cannabis Consortium (ICC), the Genetics of Personality Consortium (GPC), the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium, the Global Lipids Genetics Consortium (GLGC), The ADIPOGen Consortium, The AGEN-BMI Working Group, The CARDIOGRAMplusC4D Consortium, The CKDGen Consortium, The ICBP, The MAGIC Investigators, The MuTHER Consortium, The MIGen Consortium, The PAGE Consortium, The ReproGen Consortium, The GENIE Consortium, The International Endogene Consortium, AFGen Consortium, Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium, International Genomics of Blood Pressure (iGEN-BP) Consortium, INVENT Consortium, STARNET, COMPASS Consortium, EPIC-CVD Consortium, EPIC-InterAct Consortium, International Stroke Genetics Consortium (ISGC), METASTROKE Consortium, Neurology Working Group of the CHARGE Consortium, NINDS Stroke Genetics Network (SiGN), UK Young Lacunar DNA Study, MEGASTROKE Consortium and all the authors that contributed to generate the summary statistics used for the analyses of this study.

Financial support

This research was supported by the UK Medical Research Council (MR/N015746 and MR/S0151132). This paper represents independent research part-funded by the National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King's College London. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care. Chiara Fabbri is supported by Fondazione Umberto Veronesi (https://www.fondazioneveronesi.it).

Conflict of interest

Cathryn Lewis is a member of the Scientific Advisory Board of Myriad Neurosciences. Matthew Hotopf is the principal investigator of RADAR-CNS (funded by European Commission, Janssen, Biogen, UCB, MSD and Lundbeck). The other authors declare no potential conflict of interest.

References

Agosti, V., & Stewart, J. W. (2001). Atypical and non-atypical subtypes of depression: Comparison of social functioning, symptoms, course of illness, co-morbidity and demographic features. Journal of Affective Disorders, 65(1), 7579. https://doi.org/10.1016/S0165-0327(00)00251-2CrossRefGoogle ScholarPubMed
American Psychiatric Association. (2013). Diagnostic and statistical manual of mental disorders (5th Edn). Washington DC: American Psychiatric Association Publishing. https://doi.org/10.1176/appi.books.9780890425596.Google Scholar
Andersen, A. M., Pietrzak, R. H., Kranzler, H. R., Ma, L., Zhou, H., Liu, X., … Han, S. (2017). Polygenic scores for major depressive disorder and risk of alcohol dependence. JAMA Psychiatry, 74(11), 1153. https://doi.org/10.1001/jamapsychiatry.2017.2269CrossRefGoogle ScholarPubMed
Beijers, L., Wardenaar, K. J., van Loo, H. M., & Schoevers, R. A. (2019). Data-driven biological subtypes of depression: Systematic review of biological approaches to depression subtyping. Molecular Psychiatry, 24(6), 888900. https://doi.org/10.1038/s41380-019-0385-5CrossRefGoogle ScholarPubMed
Benazzi, F. (2002). Can only reversed vegetative symptoms define atypical depression? European Archives of Psychiatry and Clinical Neuroscience, 252(6), 288293. https://doi.org/10.1007/s00406-002-0395-0CrossRefGoogle ScholarPubMed
Blanco, C., Vesga-López, O., Stewart, J. W., Liu, S.-M., Grant, B. F., & Hasin, D. S. (2012). Epidemiology of major depression with atypical features: Results from the National Epidemiologic Survey on Alcohol and Related Conditions (NESARC). The Journal of Clinical Psychiatry, 73(02), 224232. https://doi.org/10.4088/JCP.10m06227CrossRefGoogle ScholarPubMed
Brailean, A., Curtis, J., Davis, K., Dregan, A., & Hotopf, M. (2020). Characteristics, comorbidities, and correlates of atypical depression: Evidence from the UK Biobank Mental Health Survey. Psychological Medicine, 50(7), 11291138. https://doi.org/10.1017/S0033291719001004CrossRefGoogle ScholarPubMed
Bycroft, C., Freeman, C., Petkova, D., Band, G., Elliott, L. T., Sharp, K., … Marchini, J. (2018). The UK Biobank resource with deep phenotyping and genomic data. Nature, 562(7726), 203209. https://doi.org/10.1038/s41586-018-0579-zCrossRefGoogle ScholarPubMed
Choi, S. W., & O'Reilly, P. F. (2019). PRSice-2: Polygenic Risk Score software for biobank-scale data. GigaScience, 8(7), giz082. https://doi.org/10.1093/gigascience/giz082.CrossRefGoogle ScholarPubMed
Davis, K. A. S., Coleman, J. R. I., Adams, M., Allen, N., Breen, G., Cullen, B., … Hotopf, M. (2020). Mental health in UK Biobank – development, implementation and results from an online questionnaire completed by 157 366 participants: A reanalysis. BJPsych Open, 6(2), e18. https://doi.org/10.1192/bjo.2019.100.CrossRefGoogle ScholarPubMed
Davis, K., & Hotopf, M. (2019). Mental health phenotyping in UK Biobank. Progress in Neurology and Psychiatry, 23(1), 47. https://doi.org/10.1002/pnp.522CrossRefGoogle Scholar
Euesden, J., Lewis, C. M., & O'Reilly, P. F. (2015). PRSice: Polygenic Risk Score software. Bioinformatics, 31(9), 14661468. https://doi.org/10.1093/bioinformatics/btu848CrossRefGoogle ScholarPubMed
Fried, E. I., & Nesse, R. M. (2015). Depression is not a consistent syndrome: An investigation of unique symptom patterns in the STAR*D study. Journal of Affective Disorders, 172, 96102. https://doi.org/10.1016/j.jad.2014.10.010CrossRefGoogle Scholar
Fry, A., Littlejohns, T. J., Sudlow, C., Doherty, N., Adamska, L., Sprosen, T., … Allen, N. E. (2017). Comparison of sociodemographic and health-related characteristics of UK Biobank participants with those of the general population. American Journal of Epidemiology, 186(9), 10261034. https://doi.org/10.1093/aje/kwx246CrossRefGoogle ScholarPubMed
Guo, M., & Lu, X.-Y. (2014). Leptin receptor deficiency confers resistance to behavioral effects of fluoxetine and desipramine via separable substrates. Translational Psychiatry, 4(12), e486. https://doi.org/10.1038/tp.2014.126CrossRefGoogle ScholarPubMed
Hagenaars, S., Coleman, J., Choi, S., Gaspar, H., Adams, M., Howard, D., … Lewis, C. (2019). Genetic comorbidity between major depression and cardio-metabolic disease, stratified by age at onset of major depression [Preprint]. bioRxiv. https://doi.org/10.1101/645077Google Scholar
Haro, J. M., Arbabzadeh-Bouchez, S., Brugha, T. S., De Girolamo, G., Guyer, M. E., Jin, R., … Kessler, R. C. (2006). Concordance of the Composite International Diagnostic Interview Version 3.0 (CIDI 3.0) with standardized clinical assessments in the WHO World Mental Health Surveys. International Journal of Methods in Psychiatric Research, 15(4), 167180. https://doi.org/10.1002/mpr.196CrossRefGoogle ScholarPubMed
Howard, D. M., Adams, M. J., Clarke, T.-K., Hafferty, J. D., Gibson, J., Shirali, M., … McIntosh, A. M. (2019). Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nature Neuroscience, 22(3), 343352. https://doi.org/10.1038/s41593-018-0326-7CrossRefGoogle ScholarPubMed
James, S. L., Abate, D., Abate, K. H., Abay, S. M., Abbafati, C., Abbasi, N., … Murray, C. J. L. (2018). Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: A systematic analysis for the Global Burden of Disease Study 2017. The Lancet, 392(10159), 17891858. https://doi.org/10.1016/S0140-6736(18)32279-7CrossRefGoogle Scholar
Kemp, D. E., Schinagle, M., Gao, K., Conroy, C., Ganocy, S. J., Ismail-Beigi, F., & Calabrese, J. R. (2014). PPAR-γ agonism as a modulator of mood: Proof-of-concept for pioglitazone in bipolar depression. CNS Drugs, 28(6), 571581. https://doi.org/10.1007/s40263-014-0158-2CrossRefGoogle ScholarPubMed
Kessler, R. C., Andrews, G., Mroczek, D., Ustun, B., & Wittchen, H.-U. (1998). The World Health Organization Composite International Diagnostic Interview short-form (CIDI-SF). International Journal of Methods in Psychiatric Research, 7(4), 171185. https://doi.org/10.1002/mpr.47CrossRefGoogle Scholar
Kessler, R. C., & Bromet, E. J. (2013). The epidemiology of depression across cultures. Annual Review of Public Health, 34(1), 119138. https://doi.org/10.1146/annurev-publhealth-031912-114409CrossRefGoogle ScholarPubMed
Kilpeläinen, T. O., Carli, J. F. M., Skowronski, A. A., Sun, Q., Kriebel, J., Feitosa, M. F., … Loos, R. J. F. (2016). Genome-wide meta-analysis uncovers novel loci influencing circulating leptin levels. Nature Communications, 7, 10494. https://doi.org/10.1038/ncomms10494CrossRefGoogle ScholarPubMed
Kranzler, H. R., Zhou, H., Kember, R. L., Vickers Smith, R., Justice, A. C., Damrauer, S., … Gelernter, J. (2019). Genome-wide association study of alcohol consumption and use disorder in 274424 individuals from multiple populations. Nature Communications, 10(1), 1499. https://doi.org/10.1038/s41467-019-09480-8.CrossRefGoogle Scholar
Lasserre, A. M., Glaus, J., Vandeleur, C. L., Marques-Vidal, P., Vaucher, J., Bastardot, F., … Preisig, M. (2014). Depression with atypical features and increase in obesity, body mass index, waist circumference, and fat mass: A prospective, population-based study. JAMA Psychiatry, 71(8), 880. https://doi.org/10.1001/jamapsychiatry.2014.411CrossRefGoogle Scholar
Lee, S. H., Goddard, M. E., Wray, N. R., & Visscher, P. M. (2012). A better coefficient of determination for genetic profile analysis: A better coefficient of determination. Genetic Epidemiology, 36(3), 214224. https://doi.org/10.1002/gepi.21614CrossRefGoogle ScholarPubMed
Lee, S., Ng, K. L., & Tsang, A. (2009). Prevalence and correlates of depression with atypical symptoms in Hong Kong. Australian & New Zealand Journal of Psychiatry, 43(12), 11471154. https://doi.org/10.3109/00048670903279895CrossRefGoogle ScholarPubMed
Lee, S. H., Wray, N. R., Goddard, M. E., & Visscher, P. M. (2011). Estimating missing heritability for disease from genome-wide association studies. The American Journal of Human Genetics, 88(3), 294305. https://doi.org/10.1016/j.ajhg.2011.02.002CrossRefGoogle ScholarPubMed
Levitan, R. D., Lesage, A., Parikh, S. V., Goering, P., & Kennedy, S. H. (1997). Reversed neurovegetative symptoms of depression: A community study of Ontario. The American Journal of Psychiatry, 154(7), 934940. https://doi.org/10.1176/ajp.154.7.934Google ScholarPubMed
Ligthart, S., Vaez, A., Võsa, U., Stathopoulou, M. G., de Vries, P. S., Prins, B. P., … Alizadeh, B. Z. (2018). Genome analyses of >200000 individuals identify 58 loci for chronic inflammation and highlight pathways that link inflammation and complex disorders. American Journal of Human Genetics, 103(5), 691706. https://doi.org/10.1016/j.ajhg.2018.09.009CrossRefGoogle Scholar
Lim, G. Y., Tam, W. W., Lu, Y., Ho, C. S., Zhang, M. W., & Ho, R. C. (2018). Prevalence of depression in the community from 30 countries between 1994 and 2014. Scientific Reports, 8(1), 2861. https://doi.org/10.1038/s41598-018-21243-xCrossRefGoogle ScholarPubMed
Locke, A. E., Kahali, B., Berndt, S. I., Justice, A. E., Pers, T. H., Day, F. R., … Speliotes, E. K. (2015). Genetic studies of body mass index yield new insights for obesity biology. Nature, 518(7538), 197206. https://doi.org/10.1038/nature14177CrossRefGoogle ScholarPubMed
Łojko, D., Buzuk, G., Owecki, M., Ruchała, M., & Rybakowski, J. K. (2015). Atypical features in depression: Association with obesity and bipolar disorder. Journal of Affective Disorders, 185, 7680. https://doi.org/10.1016/j.jad.2015.06.020CrossRefGoogle ScholarPubMed
Łojko, D., & Rybakowski, J. (2017). Atypical depression: Current perspectives. Neuropsychiatric Disease and Treatment, Volume, 13, 24472456. https://doi.org/10.2147/NDT.S147317CrossRefGoogle ScholarPubMed
Manichaikul, A., Mychaleckyj, J. C., Rich, S. S., Daly, K., Sale, M., & Chen, W.-M. (2010). Robust relationship inference in genome-wide association studies. Bioinformatics (Oxford, England), 26(22), 28672873. https://doi.org/10.1093/bioinformatics/btq559Google ScholarPubMed
Matza, L. S., Revicki, D. A., Davidson, J. R., & Stewart, J. W. (2003). Depression with atypical features in the National Comorbidity Survey: Classification, description, and consequences. Archives of General Psychiatry, 60(8), 817. https://doi.org/10.1001/archpsyc.60.8.817CrossRefGoogle ScholarPubMed
Milaneschi, Y., Lamers, F., Peyrot, W. J., Abdellaoui, A., Willemsen, G., Hottenga, , … Penninx, B. W. J. H. (2016). Polygenic dissection of major depression clinical heterogeneity. Molecular Psychiatry, 21(4), 516522. https://doi.org/10.1038/mp.2015.86CrossRefGoogle ScholarPubMed
Milaneschi, Y., Lamers, F., Peyrot, W. J., Baune, B. T., Breen, G., Dehghan, A., … CHARGE Inflammation Working Group and the Major Depressive Disorder Working Group of the Psychiatric Genomics Consortium. (2017). Genetic association of major depression with atypical features and obesity-related immunometabolic dysregulations. JAMA Psychiatry, 74(12), 12141225. https://doi.org/10.1001/jamapsychiatry.2017.3016CrossRefGoogle ScholarPubMed
Myers, M. G., Leibel, R. L., Seeley, R. J., & Schwartz, M. W. (2010). Obesity and leptin resistance: Distinguishing cause from effect. Trends in Endocrinology & Metabolism, 21(11), 643651. https://doi.org/10.1016/j.tem.2010.08.002CrossRefGoogle ScholarPubMed
Ring, L. E., & Zeltser, L. M. (2010). Disruption of hypothalamic leptin signaling in mice leads to early-onset obesity, but physiological adaptations in mature animals stabilize adiposity levels. Journal of Clinical Investigation, 120(8), 29312941. https://doi.org/10.1172/JCI41985CrossRefGoogle ScholarPubMed
Schumann, G., Liu, C., O'Reilly, P., Gao, H., Song, P., Xu, B., … Elliott, P. (2016). KLB Is associated with alcohol drinking, and its gene product β-Klotho is necessary for FGF21 regulation of alcohol preference. Proceedings of the National Academy of Sciences of the USA, 113(50), 1437214377. https://doi.org/10.1073/pnas.1611243113CrossRefGoogle ScholarPubMed
Sullivan, P. F., Neale, M. C., & Kendler, K. S. (2000). Genetic epidemiology of major depression: Review and meta-analysis. The American Journal of Psychiatry, 157(10), 15521562. https://doi.org/10.1176/appi.ajp.157.10.1552CrossRefGoogle ScholarPubMed
Thorp, J. G., Marees, A. T., Ong, J.-S., An, J., MacGregor, S., & Derks, E. M. (2019). Genetic heterogeneity in self-reported depressive symptoms identified through genetic analyses of the PHQ-9. Psychological Medicine, 112. https://doi.org/10.1017/S0033291719002526Google ScholarPubMed
UK Biobank (2019). UK Biobank: Improving the health of future generations. https://www.ukbiobank.ac.ukGoogle Scholar
Visscher, P. M., Hemani, G., Vinkhuyzen, A. A. E., Chen, G.-B., Lee, S. H., Wray, N. R., … Yang, J. (2014). Statistical power to detect genetic (co)variance of complex traits using SNP data in unrelated samples. PLoS Genetics, 10(4), e1004269. https://doi.org/10.1371/journal.pgen.1004269CrossRefGoogle ScholarPubMed
Warren, H. R., Evangelou, E., Cabrera, C. P., Gao, H., Ren, M., Mifsud, B., … UK Biobank CardioMetabolic Consortium BP working group. (2017). Genome-wide association analysis identifies novel blood pressure loci and offers biological insights into cardiovascular risk. Nature Genetics, 49(3), 403415. https://doi.org/10.1038/ng.3768CrossRefGoogle ScholarPubMed
Watson, H. J., Yilmaz, Z., Thornton, L. M., Hübel, C., Coleman, J. R. I., Gaspar, H. A., … Bulik, C. M. (2019). Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa. Nature Genetics, 51(8), 12071214. https://doi.org/10.1038/s41588-019-0439-2CrossRefGoogle ScholarPubMed
Wray, N. R., Ripke, S., Mattheisen, M., Trzaskowski, M., Byrne, E. M., Abdellaoui, A., … Major Depressive Disorder Working Group of the Psychiatric Genomics Consortium. (2018). Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nature Genetics, 50(5), 668681. https://doi.org/10.1038/s41588-018-0090-3CrossRefGoogle ScholarPubMed
Yang, J., Lee, S. H., Goddard, M. E., & Visscher, P. M. (2011). GCTA: A tool for genome-wide complex trait analysis. American Journal of Human Genetics, 88(1), 7682. https://doi.org/10.1016/j.ajhg.2010.11.011CrossRefGoogle ScholarPubMed
Yang, P., Tao, R., He, C., Liu, S., Wang, Y., & Zhang, X. (2018). The risk factors of the alcohol use disorders – through review of its comorbidities. Frontiers in Neuroscience, 12, 303. https://doi.org/10.3389/fnins.2018.00303.CrossRefGoogle ScholarPubMed
Zeng, J., de Vlaming, R., Wu, Y., Robinson, M. R., Lloyd-Jones, L. R., Yengo, L., … Yang, J. (2018). Signatures of negative selection in the genetic architecture of human complex traits. Nature Genetics, 50(5), 746753. https://doi.org/10.1038/s41588-018-0101-4CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Flow chart of the performed analyses. ↑WS, depression with increased weight and sleepiness; ↓WS depression, depression with decreased weight and insomnia.

Figure 1

Fig. 2. Decile plots showing the odds ratio of ↑WS depression (depression with atypical features) v. depression without ↑WS for BMI PRS (a), C-reactive protein (CRP) PRS (b), daily alcohol use PRS (c) and major depressive disorder (MDD) PRS (d). *p < 0.05; **p < 1 × 10−05.

Figure 2

Table 1. Polygenic risk scores associated with depression with atypical features (↑WS depression) compared with depression without atypical features

Figure 3

Fig. 3. PRS odds ratio (OR) and 95% confidence intervals of ↑WS depression v. depression without ↑WS. Colors indicate different groups of traits (substance related disorders, major psychiatricdisorders, personality traits, immune metabolic traits). BP, bipolar disorder; SCZ, schizophrenia; ANX, anxiety disorders; PTSD, post-traumatic stress disorder; AN, anorexia nervosa; ALCDEP, alcohol dependence; ALCUSE, daily alcohol use; N_CIGARETTES, n cigarettes per day; CANN, cannabis use lifetime; EXTR, extraversion; NEU, neuroticism; DM2, type 2 diabetes mellitus; CAD, coronary artery disease; ISCH_STROKE, ischaemic stroke; TG, triglycerides; LDL, total LDL cholesterol; HDL, total HDL cholesterol; BMI, body max index; BMI_adjust_leptin, leptin adjusted for BMI; CRP, C-reactive protein.

Figure 4

Fig. 4. PRS odds ratio (OR) and 95% confidence intervals of ↑WS depression and ↓WS depression compared with healthy controls. BP, bipolar disorder; SCZ, schizophrenia; ANX, anxiety disorders; PTSD, posttraumatic stress disorder; AN, anorexia nervosa; ALCDEP, alcohol dependence; ALCUSE, daily alcohol use; N_CIGARETTES, n cigarettes per day; CANN, cannabis use lifetime; EXTR, extraversion; NEU, neuroticism; DM2, type 2 diabetes mellitus; CAD, coronary artery disease; ISCH_STROKE, ischaemic stroke; TG, triglycerides; LDL, total LDL cholesterol; HDL, total HDL cholesterol; BMI, body max index; BMI_adjust_leptin, leptin adjusted for BMI; CRP, C-reactive protein.

Figure 5

Table 2. SNP heritability (SNP-h2) of depression with weight gain and hypersomnia (↑WS depression), depression with weight loss and reduced sleep (↓WS depression), major depression not falling in any of these two groups (no ↑WS or ↓WS), depression with weight gain or weight loss

Supplementary material: PDF

Badini et al. Supplementary Materials

Badini et al. Supplementary Materials

Download Badini et al. Supplementary Materials(PDF)
PDF 4.8 MB