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

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.


Introduction
Depressive disorders are highly prevalent and a leading cause of global disability and are associated with premature mortality (James et al., 2018;Kessler & Bromet, 2013). Twin-based heritability of unipolar depression is estimated to be ∼37% (Sullivan, Neale, & Kendler, 2000), with common single-nucleotide polymorphisms (SNPs) explaining ∼9% of variation in depression liability Wray et al., 2018). 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, 2015). 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, 2019;Milaneschi et al., 2016Milaneschi et al., , 2017. 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, 2002). 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, 2009;Matza, Revicki, Davidson, & Stewart, 2003).
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, 2001;Blanco et al., 2012;Brailean, Curtis, Davis, Dregan, & Hotopf, 2020). 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., 2020). 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, 2001;Blanco et al., 2012;Brailean et al., 2020;Lee et al., 2009;Łojko, Buzuk, Owecki, Ruchała, & Rybakowski, 2015). 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., 2020;Lasserre et al., 2014;Milaneschi et al., 2017). 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., 2017).
There is increasing evidence to suggest partially distinct genetic profiles among depressive subtypes. SNP-based heritability (h 2 SNP ) was found to vary across individual depressive symptoms (h 2 SNP range from 6% to 9%), and patterns of SNP associations and genetic correlations differed across symptoms (Thorp et al., 2019). 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., 2017). 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., 2016(Milaneschi et al., , 2017. Overlapping genetic aetiology between depression and CAD was demonstrated (Hagenaars et al., 2019), 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.

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 selfreported professional diagnoses . Genome-wide genetic data have been collected on all UKB participants (Bycroft et al., 2018), 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, 1998). 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., 2006).
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., 2020); 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., 2018). 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., 2018). 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 thirddegree kinship (KING r < 0.044; Manichaikul et al., 2010), 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., 2017). 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. (2019)].
PRS were calculated using PRSice v.2 (Choi & O'Reilly, 2019;Euesden, Lewis, & O'Reilly, 2015). 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, 2012), considering a range of possible values of prevalence among depressed cases (Levitan, Lesage, Parikh, Goering, & Kennedy, 1997;Łojko & Rybakowski, 2017). 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 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, 2011), genetic relationship matrix-restricted maximum likelihood (GREML) methods (Lee, Wray, Goddard, & Visscher, 2011) were used to estimate the variance in liability attributable to the additive effects of all SNPs (h 2 SNP ) 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., 2017), using GCTA-GREML Power Calculator (Visscher et al., 2014). 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 2 SNP 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., 2018). 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., 2018).
Possible differences among h 2 SNP of depression subgroups were compared by a Z-test as explained for PRS results. h 2 SNP estimates were transformed to a liability scale (Lee et al., 2012) and we reported heritability considering a range of plausible prevalence values (Levitan et al., 1997;Lim et al., 2018;Łojko & Rybakowski, 2017). 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. (2020), 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).

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., 2020).  (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). 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).

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., 2019).
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., 2014).
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 2 SNP was similar to a previous study of depression with weight/appetite gain [h 2 SNP = 0.11 (S.E. = 0.03)] and weight/appetite loss [h 2 SNP = 0.10 (S.E. = 0.02)] (z = 1.24, p = 0.22; z = 1.72, p = 0.09, respectively) (Milaneschi et al., 2017). 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., 2018).

732
Isabella Badini et al. (Milaneschi et al., 2017) 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, 2010). 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., 2010). 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., 2012;Brailean et al., 2020). 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., 2017;Kranzler et al., 2019). 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., 2020), and these factors may be responsible for increased risk of AUD .
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., 2017). 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, 2019). This is likely to impact on prevalence and severity of depression within the sample, with persons with severe 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.
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 selfreported 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 alcoholrelated 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. 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 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