Abstract

Depression is one of the most common mental health disorders and one of the top causes of disability throughout the world. The present study sought to identify putative causal associations between depression and hundreds of complex human traits through a genome-wide screening of genetic data and a hypothesis-free approach. We leveraged genome-wide association studies summary statistics for depression and 1504 complex traits and investigated potential causal relationships using the latent causal variable method. We identified 559 traits genetically correlated with depression risk at FDR < 5%. Of these, 46 were putative causal genetic determinants of depression, including lifestyle factors, diseases of the nervous system, respiratory disorders, diseases of the musculoskeletal system, traits related to the health of the gastrointestinal system, obesity, vitamin D levels and the use of prescription medications, among others. No phenotypes were identified as potential outcomes of depression. Our results suggest that genetic liability to multiple complex traits may contribute to a higher risk for depression. In particular, we show a putative causal genetic effect of pain, obesity and inflammation on depression. These findings provide novel insights into the potential causal determinants of depression and should be interpreted as testable hypotheses for future studies to confirm, which may facilitate the design of new prevention strategies to reduce depression’s burden.

Introduction

Depression is one of the most common mental health disorders and one of the top causes of disability throughout the world (1). Its 12-month prevalence is estimated to be 4.4% worldwide, affecting 5.1% of females and 3.6% of males (2–4). Depression is characterized by a persistent low mood or anhedonia (i.e. loss of pleasure and interest) and other symptoms such as fatigue or loss of energy, significant weight loss and sleep disturbance (5).

Epidemiological studies have identified factors that contribute to the development of depression. For instance, there is evidence that social isolation, alcohol dependence, psychological stress and medical conditions, such as cancer, are risk factors for depression (6–9). Furthermore, several hypotheses propose mechanisms to explain the pathogenesis of depression, suggesting that genes, psychosocial stress and inflammation may play a role in the development of depression (10–12).

Depression is a heritable trait, with twin and family studies estimating its heritability around 38% (13). In addition, genome-wide association studies (GWAS) have identified 223 independent genetic variants related to depression (14) and have pointed out genetic correlations between depression and other complex traits such as anxiety, autism, ADHD, schizophrenia, bipolar disorder, coronary artery disease and body fat (15–18).

Pleiotropic effects, which can be either vertical or horizontal, explain genetic correlations (19). Horizontal pleiotropy refers to genetic variants that have an independent and direct effect on both traits (i.e. owing to a shared biological pathway). In contrast, vertical pleiotropy can be understood as a causal cascade in which the effect of a genetic variant on one trait is mediated by its effect on another trait (Fig. 1) (20).

Illustration of horizontal and vertical pleiotropy concepts. (A) Horizontal pleiotropy: genetic variants independently influence the risk of traits A and B. (B) Vertical pleiotropy: genetic variants contribute to the risk of trait B through trait A, showing a potential causal relationship. rG = genetic correlation.
Figure 1

Illustration of horizontal and vertical pleiotropy concepts. (A) Horizontal pleiotropy: genetic variants independently influence the risk of traits A and B. (B) Vertical pleiotropy: genetic variants contribute to the risk of trait B through trait A, showing a potential causal relationship. rG = genetic correlation.

Mendelian randomization methods, which employ genetic variants as instrumental variables, are commonly used in genetic epidemiological studies to infer potential causal relationships between a pair of traits. One of the essential assumptions of traditional Mendelian randomization methods is that vertical pleiotropy drives the association between both traits (i.e. genetic variants are only associated with the outcome through the exposure) (21,22). If this assumption is violated through horizontal pleiotropic effects, results can be biased, increasing the risk for false-positive findings (23). Therefore, the latent causal variable (LCV) method was developed as an alternative to Mendelian randomization. LCV aims to assess potential causal associations between complex traits while accounting for potential horizontal pleiotropy. Advantages of the LCV method include that it is less susceptible to confounding by horizontal pleiotropy, it leverages genome-wide data and it is robust to sample overlap (24).

Despite the extensive efforts made to identify potential risk factors and outcomes of depression in previous observational and Mendelian randomization studies (25,26), there is still a need to uncover disease-modifying traits to screen individuals at a higher risk and enhance the design of prevention strategies. Here, we leverage a vast collection of GWAS summary statistics (N = 1504) and a larger sample size for depression than previous studies to perform LCV analyses and infer putative causal relationships between depression and other complex phenotypes. Our results should be interpreted as a set of testable hypotheses for future epidemiological and interventional studies.

Results

We identified 559 traits genetically correlated with depression risk at False Discovery Rate (FDR) < 5% (Supplementary Material, File S1). Of those traits, 37 revealed a putative causal genetic association with depression [genetic causality proportion (GCP); |GCP| > 0.60], while 9 showed evidence for limited partial genetic causality (|GCP| < 0.60). No phenotypes were identified as outcomes of depression. Putative risk factors of depression included physical health conditions, such as respiratory disorders, diseases of the musculoskeletal system, diseases of the nervous system and phenotypes related to the health of the gastrointestinal system, among others (Fig. 2). Similarly, lifestyle factors and the use of several medications increased the likelihood of depression (Table 1).

Putative causal associations for depression. Each dot represents a trait with a significant genetic correlation with depression. The x-axis shows GCP values, while the y-axis shows the GCP-related Z-scores as a measurement of statistical significance. The horizontal red dashed lines represent the statistical significance threshold (FDR < 5%). The vertical gray dashed lines divide traits that are potential causal determinants of depression (on the left) and traits that are putative consequences of depression (on the right). Results are shown separately for traits (A) with a positive genetic correlation with depression and (B) with a negative genetic correlation with depression.
Figure 2

Putative causal associations for depression. Each dot represents a trait with a significant genetic correlation with depression. The x-axis shows GCP values, while the y-axis shows the GCP-related Z-scores as a measurement of statistical significance. The horizontal red dashed lines represent the statistical significance threshold (FDR < 5%). The vertical gray dashed lines divide traits that are potential causal determinants of depression (on the left) and traits that are putative consequences of depression (on the right). Results are shown separately for traits (A) with a positive genetic correlation with depression and (B) with a negative genetic correlation with depression.

Table 1

Traits with an inferred potential causal relationship with depression

PhenotypeGCPGCP P-valuerG
Impingement syndrome of the shoulder-0.921.69E−060.29
Taking drugs or alcohol for anxiety more than once-0.901.42E−220.81
Other diseases of anus and rectum (ICD10)-0.887.18E−160.57
Medication: diazepam-0.851.24E−100.57
Umbilical hernia (ICD10)-0.835.87E−110.21
Other and unspecified injuries of head (ICD10)-0.832.52E−020.30
Self-reported H. pylori-0.822.30E−080.38
Chest pain felt during physical activity-0.794.73E−060.52
Medication: beclomethasone-0.782.15E−060.30
Self-reported psoriasis-0.772.15E−060.19
Leg pain when walking uphill or hurrying-0.751.08E−040.37
Diaphragmatic hernia (ICD10)-0.751.06E−040.41
Never been injured or injured someone else through drinking alcohol-0.751.86E−04-0.36
Age of stopping smoking-0.742.11E−040.35
Vitamin D (male)-0.731.69E−06-0.17
Medication: ranitidine-0.732.12E−030.55
Breathing problems did not improve/stopped away from workplace or on holiday-0.731.37E−370.23
Triglycerides-0.727.97E−090.15
Medication: fluoxetine-0.721.98E−020.82
Self-reported ulcerative colitis-0.713.58E−050.19
Obesity-0.705.39E−030.32
Polyneuropathies and other disorders of the peripheral nervous system-0.702.94E−030.25
Diseases of the ear and mastoid process-0.703.85E−030.32
Rarely/never worked with materials containing asbestos-0.707.69E−03-0.43
Fruit consumers-0.699.30E−04-0.32
Osteoporosis without pathological fracture (ICD10)-0.672.48E−020.35
Ever addicted to any substance or behavior-0.672.73E−020.54
Triglycerides (male)-0.666.12E−040.15
Vitamin D (both sexes)-0.663.49E−05-0.15
Medication: dihydrocodeine-0.664.19E−020.49
Other diseases of intestine (ICD10)-0.654.92E−020.40
Sleeping too much-0.652.23E−020.37
Age hay fever or allergic rhinitis diagnosed by doctor-0.647.69E−030.22
Medication: gabapentin-0.633.64E−020.32
Mental and behavioral disorders owing to alcohol (ICD10)-0.612.15E−020.23
Numeric addition test-0.614.92E−02-0.26
Malignant neoplasm of bronchus and lung (ICD10)-0.602.86E−170.29
PhenotypeGCPGCP P-valuerG
Impingement syndrome of the shoulder-0.921.69E−060.29
Taking drugs or alcohol for anxiety more than once-0.901.42E−220.81
Other diseases of anus and rectum (ICD10)-0.887.18E−160.57
Medication: diazepam-0.851.24E−100.57
Umbilical hernia (ICD10)-0.835.87E−110.21
Other and unspecified injuries of head (ICD10)-0.832.52E−020.30
Self-reported H. pylori-0.822.30E−080.38
Chest pain felt during physical activity-0.794.73E−060.52
Medication: beclomethasone-0.782.15E−060.30
Self-reported psoriasis-0.772.15E−060.19
Leg pain when walking uphill or hurrying-0.751.08E−040.37
Diaphragmatic hernia (ICD10)-0.751.06E−040.41
Never been injured or injured someone else through drinking alcohol-0.751.86E−04-0.36
Age of stopping smoking-0.742.11E−040.35
Vitamin D (male)-0.731.69E−06-0.17
Medication: ranitidine-0.732.12E−030.55
Breathing problems did not improve/stopped away from workplace or on holiday-0.731.37E−370.23
Triglycerides-0.727.97E−090.15
Medication: fluoxetine-0.721.98E−020.82
Self-reported ulcerative colitis-0.713.58E−050.19
Obesity-0.705.39E−030.32
Polyneuropathies and other disorders of the peripheral nervous system-0.702.94E−030.25
Diseases of the ear and mastoid process-0.703.85E−030.32
Rarely/never worked with materials containing asbestos-0.707.69E−03-0.43
Fruit consumers-0.699.30E−04-0.32
Osteoporosis without pathological fracture (ICD10)-0.672.48E−020.35
Ever addicted to any substance or behavior-0.672.73E−020.54
Triglycerides (male)-0.666.12E−040.15
Vitamin D (both sexes)-0.663.49E−05-0.15
Medication: dihydrocodeine-0.664.19E−020.49
Other diseases of intestine (ICD10)-0.654.92E−020.40
Sleeping too much-0.652.23E−020.37
Age hay fever or allergic rhinitis diagnosed by doctor-0.647.69E−030.22
Medication: gabapentin-0.633.64E−020.32
Mental and behavioral disorders owing to alcohol (ICD10)-0.612.15E−020.23
Numeric addition test-0.614.92E−02-0.26
Malignant neoplasm of bronchus and lung (ICD10)-0.602.86E−170.29

This table shows traits with a significant (FDR < 5%) and strong GCP (|GCP| >  0.60) with depression. Owing to space restrictions, all nominally significant genetic correlations (i.e. P-value <  0.05 before multiple testing correction) for depression are shown in Supplementary Material, File S1. Phenotype, phenotype with a potential causal association with depression; GCP P-value, GCP P-value after multiple testing correction; rG, genetic correlation.

Table 1

Traits with an inferred potential causal relationship with depression

PhenotypeGCPGCP P-valuerG
Impingement syndrome of the shoulder-0.921.69E−060.29
Taking drugs or alcohol for anxiety more than once-0.901.42E−220.81
Other diseases of anus and rectum (ICD10)-0.887.18E−160.57
Medication: diazepam-0.851.24E−100.57
Umbilical hernia (ICD10)-0.835.87E−110.21
Other and unspecified injuries of head (ICD10)-0.832.52E−020.30
Self-reported H. pylori-0.822.30E−080.38
Chest pain felt during physical activity-0.794.73E−060.52
Medication: beclomethasone-0.782.15E−060.30
Self-reported psoriasis-0.772.15E−060.19
Leg pain when walking uphill or hurrying-0.751.08E−040.37
Diaphragmatic hernia (ICD10)-0.751.06E−040.41
Never been injured or injured someone else through drinking alcohol-0.751.86E−04-0.36
Age of stopping smoking-0.742.11E−040.35
Vitamin D (male)-0.731.69E−06-0.17
Medication: ranitidine-0.732.12E−030.55
Breathing problems did not improve/stopped away from workplace or on holiday-0.731.37E−370.23
Triglycerides-0.727.97E−090.15
Medication: fluoxetine-0.721.98E−020.82
Self-reported ulcerative colitis-0.713.58E−050.19
Obesity-0.705.39E−030.32
Polyneuropathies and other disorders of the peripheral nervous system-0.702.94E−030.25
Diseases of the ear and mastoid process-0.703.85E−030.32
Rarely/never worked with materials containing asbestos-0.707.69E−03-0.43
Fruit consumers-0.699.30E−04-0.32
Osteoporosis without pathological fracture (ICD10)-0.672.48E−020.35
Ever addicted to any substance or behavior-0.672.73E−020.54
Triglycerides (male)-0.666.12E−040.15
Vitamin D (both sexes)-0.663.49E−05-0.15
Medication: dihydrocodeine-0.664.19E−020.49
Other diseases of intestine (ICD10)-0.654.92E−020.40
Sleeping too much-0.652.23E−020.37
Age hay fever or allergic rhinitis diagnosed by doctor-0.647.69E−030.22
Medication: gabapentin-0.633.64E−020.32
Mental and behavioral disorders owing to alcohol (ICD10)-0.612.15E−020.23
Numeric addition test-0.614.92E−02-0.26
Malignant neoplasm of bronchus and lung (ICD10)-0.602.86E−170.29
PhenotypeGCPGCP P-valuerG
Impingement syndrome of the shoulder-0.921.69E−060.29
Taking drugs or alcohol for anxiety more than once-0.901.42E−220.81
Other diseases of anus and rectum (ICD10)-0.887.18E−160.57
Medication: diazepam-0.851.24E−100.57
Umbilical hernia (ICD10)-0.835.87E−110.21
Other and unspecified injuries of head (ICD10)-0.832.52E−020.30
Self-reported H. pylori-0.822.30E−080.38
Chest pain felt during physical activity-0.794.73E−060.52
Medication: beclomethasone-0.782.15E−060.30
Self-reported psoriasis-0.772.15E−060.19
Leg pain when walking uphill or hurrying-0.751.08E−040.37
Diaphragmatic hernia (ICD10)-0.751.06E−040.41
Never been injured or injured someone else through drinking alcohol-0.751.86E−04-0.36
Age of stopping smoking-0.742.11E−040.35
Vitamin D (male)-0.731.69E−06-0.17
Medication: ranitidine-0.732.12E−030.55
Breathing problems did not improve/stopped away from workplace or on holiday-0.731.37E−370.23
Triglycerides-0.727.97E−090.15
Medication: fluoxetine-0.721.98E−020.82
Self-reported ulcerative colitis-0.713.58E−050.19
Obesity-0.705.39E−030.32
Polyneuropathies and other disorders of the peripheral nervous system-0.702.94E−030.25
Diseases of the ear and mastoid process-0.703.85E−030.32
Rarely/never worked with materials containing asbestos-0.707.69E−03-0.43
Fruit consumers-0.699.30E−04-0.32
Osteoporosis without pathological fracture (ICD10)-0.672.48E−020.35
Ever addicted to any substance or behavior-0.672.73E−020.54
Triglycerides (male)-0.666.12E−040.15
Vitamin D (both sexes)-0.663.49E−05-0.15
Medication: dihydrocodeine-0.664.19E−020.49
Other diseases of intestine (ICD10)-0.654.92E−020.40
Sleeping too much-0.652.23E−020.37
Age hay fever or allergic rhinitis diagnosed by doctor-0.647.69E−030.22
Medication: gabapentin-0.633.64E−020.32
Mental and behavioral disorders owing to alcohol (ICD10)-0.612.15E−020.23
Numeric addition test-0.614.92E−02-0.26
Malignant neoplasm of bronchus and lung (ICD10)-0.602.86E−170.29

This table shows traits with a significant (FDR < 5%) and strong GCP (|GCP| >  0.60) with depression. Owing to space restrictions, all nominally significant genetic correlations (i.e. P-value <  0.05 before multiple testing correction) for depression are shown in Supplementary Material, File S1. Phenotype, phenotype with a potential causal association with depression; GCP P-value, GCP P-value after multiple testing correction; rG, genetic correlation.

Physical health conditions inferred to increase risk of depression included malignant neoplasm of bronchus and lung (ICD10) (GCP = −0.60, rG = 0.29), self-reported ulcerative colitis (GCP = −0.83, rG = 0.19), high triglyceride levels (GCP = −0.72, rG = 0.15), self-reported susceptibility to Helicobacter pylori infection (GCP = −0.82, rG = 0.38), self-reported psoriasis (GCP = −0.77, rG = 0.19), polyneuropathies and disorders of the peripheral nervous system (GCP = −0.70, rG = 0.25) and obesity (GCP = −0.70, rG = 0.32), among others. In contrast, high vitamin D levels were inferred to potentially decrease the risk of depression (GCP = −0.66, rG = −0.15).

Among lifestyle-related phenotypes, we observed that being a fruit consumer (GCP = −0.69, rG = −0.32) or rarely/never worked with materials containing asbestos (GCP = −0.70, rG = −0.43) decreased the risk for depression. In addition, potential vertical pleiotropic effects were observed between phenotypes related to alcohol consumption and depression. These included substances taken for anxiety: drugs or alcohol (more than once) (GCP = −0.90, rG = 0.81) and mental and behavioral disorders owing to use of alcohol (ICD10) (GCP = −0.61, rG = 0.23), increasing the likelihood of depression. Consistently, never being injured or injuring someone else through drinking alcohol (GCP = −0.75, rG = −0.36) decreased the risk of depression.

The use of several medications potentially increased the risk of depression. These medications included gabapentin (GCP = −0.63, rG = 0.32), which is usually prescribed for epilepsy and peripheral neuropathy; dihydrocodeine (GCP = −0.66, rG = 0.49), an opioid analgesic indicated for moderate to severe pain; beclometasone (GCP = −0.78, rG = 0.30), a corticosteroid prescribed for several conditions including asthma; and diazepam (GCP = −0.85, rG = 0.57), which is commonly prescribed for anxiety-related traits. In addition, medications prescribed for gastrointestinal ulcers, such as gaviscon liquid (GCP = −0.52, rG = 0.37) and ranitidine (GCP = −0.73, rG = 0.55), potentially increased the likelihood of depression.

Sensitivity analysis

As a sensitivity analysis, we performed LCV analyses to further explore potential causal associations identified in this study using publicly available GWAS summary statistics with larger sample sizes than those available in CTG-VL (Table 2). Overall, we did not observe any direct potential causal effects between depression and substance use–related phenotypes, such as cigarettes per day, smoking cessation, smoking initiation and drinks per week. Similarly, despite the large genetic overlap between phenotypes (rG = 0.91), we did not identify potential causal effects between lifetime anxiety and depression. Nonetheless, we replicated the causal association in which high vitamin D levels potentially decrease the risk for depression (GCP = −0.59, rG = −0.12). Consistent with our previous results, our sensitivity analysis failed to identify evidence for a causal association between BMI and depression.

Table 2

Latent casual variable results for depression and other complex traits in sensitivity analysis

PhenotypeGCPGCP P-valuerG
Vitamin D-0.591.00E−300-0.12
Cigarettes per day0.060.53-0.23
Smoking initiation-0.040.54-0.31
Drinks per week0.110.61-0.06
Smoking cessation-0.070.62-0.24
BMI0.260.870.11
PhenotypeGCPGCP P-valuerG
Vitamin D-0.591.00E−300-0.12
Cigarettes per day0.060.53-0.23
Smoking initiation-0.040.54-0.31
Drinks per week0.110.61-0.06
Smoking cessation-0.070.62-0.24
BMI0.260.870.11

This table shows results for the LCV method between depression and other complex traits.

Table 2

Latent casual variable results for depression and other complex traits in sensitivity analysis

PhenotypeGCPGCP P-valuerG
Vitamin D-0.591.00E−300-0.12
Cigarettes per day0.060.53-0.23
Smoking initiation-0.040.54-0.31
Drinks per week0.110.61-0.06
Smoking cessation-0.070.62-0.24
BMI0.260.870.11
PhenotypeGCPGCP P-valuerG
Vitamin D-0.591.00E−300-0.12
Cigarettes per day0.060.53-0.23
Smoking initiation-0.040.54-0.31
Drinks per week0.110.61-0.06
Smoking cessation-0.070.62-0.24
BMI0.260.870.11

This table shows results for the LCV method between depression and other complex traits.

Discussion

In this work, we conducted a phenome-wide screening and LCV analyses between depression and 1504 other complex traits to shed light on the putative causal architecture of depression. Previous genetic epidemiological studies have proposed different risk factors for depression using traditional Mendelian randomization methods (25,27,28). For instance, inflammatory gastrointestinal diseases, anxiety disorders and asthma have been associated with a higher risk for depression (27,29). Here, we sought to assess the robustness to the modeling of horizontal pleiotropy of findings from previous studies and provide novel insights into the potential causal determinants and outcomes of depression by uncovering a set of testable hypotheses for future studies.

In the present study, traits such as difficulty stopping worrying during the worst period of anxiety and taking drugs or alcohol for anxiety more than once potentially increased the risk of depression. In particular, for the first phenotype, half of its genetic correlation with depression was explained by potential vertical pleiotropic effects, while for the latter, almost all of the genetic overlap is consistent with vertical pleiotropic effects. These symptoms are observed in anxiety disorders, which are highly comorbid with depression (30,31). For instance, previous studies suggested that both disorders have a substantial genetic overlap (15,32). In the present study, we did not identify potential vertical pleiotropic effects between lifetime anxiety and depression. Although it can be challenging to establish causal associations between co-occurring disorders, it could be the case here that indirect anxiety-related phenotypes serve as a proxy for depression. In addition, it has been suggested that taking drugs or alcohol to overcome depression and anxiety may exacerbate their symptoms and influence the development of other mental disorders, perhaps as a consequence of intoxication or withdrawal (7,33).

Mental and behavioral disorders owing to alcohol (ICD10) were a potential causal determinant for depression in which 60% of the genetic overlap with depression is explained by vertical pleiotropic effects. In addition, in the present study, never being injured or injuring someone else through drinking alcohol, which indirectly suggests a lower alcohol consumption, decreased the risk of depression (three-quarters of its genetic correlation with depression was explained by potential vertical pleiotropic effects). Previous studies have observed that depressive episodes were prevalent after alcohol drinking, withdrawal and intoxication (7,33–35), which suggests that a higher alcohol intake is associated with an increase in the incidence of depression (36). However, our sensitivity analysis did not show any potential causal effects between drinks per week and depression. Therefore, another plausible explanation for the indirect effect we originally observed of alcohol on depression could be a common genetic factor or a shared environment between heavy alcohol drinkers and depression patients (37). Regarding the relationship between smoking and depression, our sensitivity analysis did not show any potential causal effects between smoking phenotypes (cigarettes per day, smoking initiation and smoking cessation) and depression.

It has been shown that depression has a higher incidence and prevalence among lung cancer patients (9,38–40). A likely explanation for this is that lung cancer patients face a relatively low survival rate of around 21%, which has a detrimental psychological effect (41,42). In some cases, patients may be subjected to high financial stress related to treatment costs, which has been associated with an increased risk of developing depressive symptoms (41,42). Our results show malignant neoplasm of bronchus and lung (ICD10) and lung cancer and mesothelioma, potentially increasing the likelihood of depression, perhaps as a consequence of a detrimental psychological effect on patients upon diagnosis. Both phenotypes revealed a moderate genetic correlation with depression. In particular, for the first phenotype, three-fifths of the genetic overlap was explained by potential vertical pleiotropic effects, while for the latter, 46% of the genetic overlap is consistent with potential vertical pleiotropic effects.

There is accumulating observational and genetic evidence showing that atypical features of depression, such as increased appetite and/or weight during an active depressive episode, are associated with obesity-related traits (43). In the present study, obesity influenced a higher risk of depression (70% of its genetic correlation with depression is consistent with vertical pleiotropic effects), which may be explained by psychological and inflammatory factors. From a psychological perspective, obesity can affect self-esteem and body image perception. It could, in some instances, lead to depressive symptoms, particularly for individuals in communities where their body image is compared with cultural beauty standards, which are often defined by a low body mass index (44–47). On the other side, inflammation has also been suggested as a putative mechanism underlying depression (48). For instance, patients with depression are known to show higher levels of inflammatory biomarkers in peripheral blood (10), and it is well established that obesity leads to a chronic inflammatory state (49) which in turn increases the production of pro-inflammatory cytokines and c-reactive protein (50). Therefore, our results align with the hypothesis that an inflammatory state may contribute to a higher risk of depression, perhaps as a consequence of obesity. Nonetheless, further research is required to pinpoint the potential molecular underpinnings underlying the association between depression and inflammation.

Previous studies have reported an association between high triglycerides, low HDL cholesterol levels and high risk for depression (51,52), while some other studies found an inverse association with depression (53). Our results support the hypothesis in which high HDL cholesterol levels could be protective for depression (Supplementary Material, File S1), while high triglyceride levels might increase the likelihood of depression (54). In particular, for the first phenotype, half of its genetic overlap with depression was explained by potential vertical pleiotropic effects, while for the latter, around 70% of the genetic overlap is consistent with vertical pleiotropic effects. This finding is most likely explained by obesity-related effects, which are well known to lower HDL cholesterol and increase triglyceride levels (55,56).

The association between chronic pain and depression has been investigated before (57,58). For instance, it has been reported that neuropathic pain and depression share common pathogenic and inflammatory pathways, and neuroimaging findings suggest that common brain areas are responsible for regulating emotions and pain experiences (59). Chronic pain is associated with higher rates of depression, and comorbid chronic pain and depression are associated with poorer outcomes, including higher rates of suicidality and unsatisfactory antidepressant response both on the phenotypic (60) and genetic levels (61). Pain related to chronic musculoskeletal disorders has been associated with a high psychological distress level (62). In the present study, musculoskeletal disorders, such as meniscus derangement, impingement syndrome of the shoulder and osteoporosis without pathological fracture (ICD10), were potential causal determinants of depression, as were disorders of the peripheral nervous system. In particular, for the first phenotype, two-fifths of its genetic overlap with depression was explained by potential vertical pleiotropic effects, while for the last two phenotypes, 66 and 70% of the genetic overlap are explained by potential vertical pleiotropic effects, respectively. In addition, almost all of the genetic overlap for the impingement syndrome of the shoulder phenotype with depression is consistent with potential vertical pleiotropic effects. Consistently, medications usually prescribed for easing chronic pain symptoms, such as gabapentin and dihydrocodeine, which can be used as proxies for neuropathic pain and moderate to severe pain, respectively, potentially increased the risk of depression. These findings suggest that identifying and addressing pain as a potential risk factor for depression could be of great importance.

We observed putative causal associations in which gastrointestinal disorders increased the risk for depression. These included susceptibility to H. pylori infection, ulcerative colitis and the use of gaviscon liquid and ranitidine, which can be interpreted as proxies for peptic ulcers (63,64). In particular, for the previously mentioned traits, 82, 71, 52 and 73% of their genetic overlap with depression was explained by potential vertical pleiotropic effects, respectively. Previous studies have proposed a mechanism underlying the relationship between H. pylori infection and depression, suggesting that the infection plays a role in immunity and serotonin receptors upregulation, which may increase the risk of psychiatric disorders, including depression (65). Furthermore, evidence from a previous Mendelian randomization study supported that depression could potentially increase the risk of developing peptic ulcer disease. However, combined gastrointestinal disorders (i.e. peptic ulcer and inflammatory bowel diseases) showed potential causal associations with depression in both directions (66). Our results suggest that gastrointestinal disorders could be a putative causal determinant of depression; this association is most likely explained by immune-inflammatory pathways (12,67) and the detrimental psychological effect of physical illness on vulnerable individuals (68).

Previous studies have investigated the potential supportive role of vitamin D supplementation (69,70), suggesting that vitamin D has homeostatic and immunomodulatory effects (71). Furthermore, a Mendelian randomization study found a putative causal indirect effect of depression decreasing vitamin D concentrations, which is most likely explained by behaviors that are well known to lead to reduced production of vitamin D, such as less outdoor and physical activities decreasing sun exposure (72). Our results are consistent with previous findings in which vitamin D may have a protective direct effect on depression, and our results suggest that two-thirds of the genetic correlation between vitamin D and depression was explained by vertical pleiotropic effects. Nonetheless, we cannot rule out the possibility of sun exposure (i.e. through outdoor activities) underlying this association.

Previous studies have observed that a low socioeconomic status (SES) is associated with a higher prevalence of depression (73,74). Individuals with lower SES are more likely to work in an environment exposed to hazardous substances and involuntarily tend to consume less healthy food, as its price is more accessible compared with nutritious meals (75–77). In the present study, traits such as being a fruit consumer and rarely/never worked with materials containing asbestos influenced a potentially lower risk for depression. Around 70% of the genetic overlap of each of the aforementioned phenotypes with depression is consistent with potential vertical pleiotropic effects. These findings add up to previous reports of a positive association between higher fruit consumption and better mental health (78,79), suggesting that fruit, which is rich in minerals and vitamins, may have a role in modulating neurotransmitter receptors associated with depressive symptoms (80). In addition, fruit consumption may protect the body from oxidative stress, which has been suggested to contribute to depression through inflammation and neurodegeneration (81–83).

Similarly, previous studies suggest that depression is closely related to occupational environments (84). In fact, depression is prevalent among individuals working with asbestos (85). Our results are consistent with observations in previous studies, suggesting a putative protective effect of healthy food choices and non-hazardous work environments on depression risk.

Epidemiological studies commonly use randomized controlled trials (RCTs) as the standard approach to identify risk factors of several diseases. However, RCTs have shortcomings as it could be expensive to conduct large RCTs, or it may be unethical to randomly allocate participants to a well-known harmful exposure or give a placebo to seriously ill patients (86,87). In addition, observational studies, such as case–control or cohort studies, have limited capacity to identify causal associations owing to potential confounding variables and biases (88). As an alternative to this, causal inference methods used in statistical genetics are based on genetic data to provide a different approach in the light of the limitations of interventional and observational studies (21). As previous studies have noted, we note that potential causal associations identified with the LCV method should be used as testable hypotheses for future studies to confirm. Therefore, we stress the importance of triangulating results from several studies with different study designs before establishing a causal relationship between complex traits.

We recognize that there are potential limitations to our study. First, most GWAS summary statistics used in the present study were retrieved from the UK Biobank cohort, which primarily includes participants of European ancestry. Therefore, results should not be generalized to other populations until findings are confirmed using data for different populations. Second, owing to its nature, the LCV method assesses the predominant causal pathway between two phenotypes and assumes no bidirectional causality (24). However, it is unclear whether bidirectional associations would still be captured by the LCV. In the future, this assumption should be tested in a systematic manner. Null findings in the present study may be explained owing to a lack of power because there is truly no causal relationship or because the relationship is bidirectional. Null findings in the present study must be addressed with caution.

Third, it is important to consider the complexity of genotype–phenotype associations, in which a GWAS may capture other traits beyond the one under study, potentially reflecting other risk factors. For instance, it is likely that the relationship between fruit consumption and depression is influenced by SES-related variables, and the association between lung cancer and depression is most likely explained by a detrimental psychological effect upon diagnosis and financial stress. In addition, in the present study, we identified a potential causal association between depression and obesity but not between depression and BMI. This discrepancy has been observed in previous studies (49) and is explained by the difference in phenotype definition. For example, obesity refers to the dichotomization of BMI based on the International Classification of Diseases, whereas BMI represents BMI as a continuous measurement which in turn does not entirely reflect obesity as a disease. Another example for this is that, as suggested in previous studies (89), phenotypes related to the use of medications were interpreted as a proxy for the trait they are most commonly prescribed for.

Lastly, although the LCV method offers increased statistical power owing to the use of genetic information across the whole genome, it is tied to the statistical power of the original GWAS included in the analysis. Thus, the ability to infer potential causal associations for phenotypes with an underpowered GWAS is limited. This might explain the absence of evidence for potential outcomes of depression, particularly the absence of a statistically significant association for suicide attempt as a potential outcome of depression (Supplementary Material, File S1), which comes in contrast with documented evidence in which suicidal behaviors are frequent among patients with depressive symptoms (90). This could be owing to the low statistical power of the suicide attempt GWAS. Nonetheless, previous studies report that the association between suicidal behaviors and depression could be confounded by traits such as hopelessness (91), which adds up to evidence suggesting that the association between suicide and depression is likely explained by horizontal pleiotropic effects (92). In addition, the etiology of individual suicide-related phenotypes, such as suicide ideation, suicide attempt and death by suicide, might differ and thus cannot all be captured by the suicide attempt GWAS only (93).

In summary, we provide evidence for potential causal genetic effects between depression and 1504 complex phenotypes. Our findings uncovered respiratory and gastrointestinal disorders potentially increasing depression risk. Also, we identified that obesity and pain might increase the likelihood of depression. Our analyses provide evidence for causal relationships observed in previous studies and provide novel testable hypotheses, increasing our understanding of depression and contributing to the design of new prevention and treatment strategies to reduce depression’s burden.

Materials and Methods

Depression dataset

We leveraged publicly available GWAS summary statistics for depression from the Psychiatric Genomics Consortium (PGC). Briefly, an inverse-variance weighted meta-analysis was performed on samples from the UK Biobank, the PGC and 23andMe. After excluding the data from 23andMe, the total sample size was 500 199 individuals of which 170 756 were cases and 329 443 were controls. Full details for these GWAS summary statistics are available at their corresponding publication (17). The summary statistics for depression can be freely accessed via the PGC’s website (https://www.med.unc.edu/pgc/download-results/).

Other datasets

This study used a compilation of 1504 phenotypes available at the Complex Traits Genetics Virtual Lab (CTG-VL) web-based tool (https://genoma.io) (94), which provides free access to public GWAS summary statistics and post-GWAS analyses to improve research collaboration and reproducibility (94). Most of the GWAS summary statistics available in CTG-VL were retrieved from the second wave of GWAS results from the UK Biobank released by Neale’s Lab (www.nealelab.is/uk-biobank/) (95), and the rest were retrieved from GWAS consortia. Therefore, most of these GWAS summary statistics represent individuals of European ancestry which in turn prevents putative biases owing to differences between genetic ancestries (94). Phenotypes in the CTG-VL include objective laboratory measurements, self-reported traits and consortia meta-analysis. UK Biobank GWAS were adjusted for age, age-squared, inferred sex (i.e. sex is inferred from the sex chromosomes in the genome), age * inferred sex, age-squared * inferred sex and 20 genetic ancestry principal components (94,95). CTG-VL platform included only GWAS with a nominally significant heritability (i.e. heritability estimate with a P-value < 0.05 before multiple testing correction) derived from LD-score regression to ensure that analyses such as LCV and genetic correlations can be performed (94).

Latent causal variable

The LCV method determines whether a genetic correlation between traits can be explained by a putative causal relationship rather than horizontal pleiotropic effects. The method fits a latent variable L that mediates the genetic correlation between two traits (trait A and trait B), assuming that L is causal for both traits (24). In particular, L represents the degree to which a one-way potential causal relationship explains a genetic correlation (i.e. the genetic effects of trait A are proportional on B but not vice versa). Furthermore, the GCP parameter, which is estimated by comparing the correlation between L and trait A and the correlation between L and trait B (24), indicates whether the genetic correlation is more likely explained by horizontal or vertical pleiotropy.

Under the assumption of no bi-directional causality, a GCP value equal to zero suggests the genetic correlation is explained by horizontal pleiotropic effects, implying the absence of potential causal genetic effects. In contrast, a |GCP| = 1 indicates the detection of vertical pleiotropic effects and full genetic causality. Furthermore, |GCP| values between 0 and 1 indicate partial genetic causality (24). A |GCP| > 0.6 is considered to be robust, which indicates a lower likelihood of false-positive findings, and interventions targeting trait A are likely to affect trait B owing to potential causal genetic effects (24). Therefore, to interpret LCV results, it is important to take into consideration three factors: the magnitude of the genetic correlation, the value of the GCP estimate and the direction of the GCP estimate, which can be positive or negative. The GCP value is not a measure of the magnitude of the potential causal effects but indicates the proportion of a genetic correlation that can be explained by potential causal effects. For instance, associations with a small genetic correlation (|rG| < 0.30) and a large |GCP| estimate close to 1 suggest that almost all of the genetic overlap between two traits, even though it is small, is explained by potential vertical pleiotropic effects. In addition, in the present study, a negative GCP between depression and trait A suggests potential causal genetic effects from trait A on depression, while a positive GCP value between depression and trait A would imply potential causal genetic effects from depression on trait A.

In the present study, we estimated the GCP between depression and 1504 complex traits. We used the phenome-wide analysis pipeline publicly available in CTG-VL as described in previous studies (96,97). The phenome-wide analysis pipeline was created using the original R script for the LCV method made available by the original authors in a GitHub repository (https://github.com/lukejoconnor/LCV) (24). Consistency of alleles and variants for all GWAS summary statistics was ensured using munge_sumstats.py, which is available in the LD score regression software. Also, HapMap SNPs were extracted using the list of SNPs w_hm3_noMHC.snplist (https://github.com/bulik/ldsc/wiki). In this study, to conduct the analyses, we uploaded GWAS summary statistics for depression to CTG-VL. Then, we used the phenome-wide analysis pipeline to perform LD score regression (98) and LCV analyses (24) to estimate genetic correlations and potential causal relationships, respectively. Finally, we created causal architecture plots to visualize the results. Full details on using and interpreting results from the publicly available phenome-wide analysis pipeline are available elsewhere (96,97). The LCV was applied to all phenotypes that showed a genetic correlation with depression at Benjamini-Hochberg FDR < 5%. Similarly, we applied FDR < 5% on all GCP estimates obtained from the LCV method.

Sensitivity analysis

As a sensitivity analysis, we performed LCV analyses to further explore potential causal associations identified in this study using publicly available GWAS summary statistics with larger sample sizes. These included cigarettes per day (N = 403 928) (99), smoking cessation (N = 820 192) (99), smoking initiation (N = 1 359 002) (99), drinks per week (N = 1 039 210) (99), vitamin D (N = 417 580) (72), lifetime anxiety (N = 83 566) (100) and BMI (N = ∼700 000) (101).

Data availability statement

Summary-level data used in the present study are publicly available at the Complex Traits Genomics Virtual Lab (https://genoma.io/) platform. Summary-level data for depression are publicly available at the PGC website (https://www.med.unc.edu/pgc/shared-methods/open-source-philosophy/).

Funding

L.M.G.-M. is supported by UQ Research Training Scholarship from The University of Queensland (UQ).

Conflict of Interest statement. G.C.-P. contributed to this study while employed at The University of Queensland. He is now an employee of 23andMe Inc., and he may hold stock or stock options. All other authors declare no conflicts of interest.

References

1.

GBD 2017
and
Disease and Injury Incidence and Prevalence Collaborators
(
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
.
Lancet
,
392
,
1789
1858
.

2.

Depression, WHO.
(
2017
)
Other Common Mental Disorders: Global Health Estimates
.
World Health Organization
,
Geneva
, pp.
1
24
.

3.

Gutiérrez-Rojas
,
L.
,
Porras-Segovia
,
A.
,
Dunne
,
H.
,
Andrade-González
,
N.
and
Cervilla
,
J.A.
(
2020
)
Prevalence and correlates of major depressive disorder: a systematic review
.
Braz. J. Psychiatry
,
42
,
657
672
.

4.

Bromet
,
E.
,
Andrade
,
L.H.
,
Hwang
,
I.
,
Sampson
,
N.A.
,
Alonso
,
J.
,
De Girolamo
,
G.
,
De Graaf
,
R.
,
Demyttenaere
,
K.
,
Hu
,
C.
,
Iwata
,
N.
et al. (
2011
)
Cross-national epidemiology of DSM-IV major depressive episode
.
BMC Med.
,
9
,
1
16
.

5.

American Psychiatric Association
(
2013
)
Diagnostic and Statistical Manual of Mental Disorders (DSM-5®)
.
American Psychiatric Pub
, United States of America, p.
2013
.

6.

Matthews
,
T.
,
Danese
,
A.
,
Wertz
,
J.
,
Odgers
,
C.L.
,
Ambler
,
A.
,
Moffitt
,
T.E.
and
Arseneault
,
L.
(
2016
)
Social isolation, loneliness and depression in young adulthood: a behavioural genetic analysis
.
Soc. Psychiatry Psychiatr. Epidemiol.
,
51
,
339
348
.

7.

Kuria
,
M.W.
,
Ndetei
,
D.M.
,
Obot
,
I.S.
,
Khasakhala
,
L.I.
,
Bagaka
,
B.M.
,
Mbugua
,
M.N.
and
Kamau
,
J.
(
2012
)
The association between alcohol dependence and depression before and after treatment for alcohol dependence
.
ISRN Psychiatry
,
2012
,
482802
.

8.

Yang
,
L.
,
Zhao
,
Y.
,
Wang
,
Y.
,
Liu
,
L.
,
Zhang
,
X.
,
Li
,
B.
and
Cui
,
R.
(
2015
)
The effects of psychological stress on depression
.
Curr. Neuropharmacol.
,
13
,
494
504
.

9.

Hung
,
M.-S.
,
Chen
,
I.-C.
,
Lee
,
C.-P.
,
Huang
,
R.-J.
,
Chen
,
P.-C.
,
Tsai
,
Y.-H.
and
Yang
,
Y.-H.
(
2017
)
Incidence and risk factors of depression after diagnosis of lung cancer: a nationwide population-based study
.
Medicine
,
96
,
e6864
.

10.

Miller
,
A.H.
,
Maletic
,
V.
and
Raison
,
C.L.
(
2009
)
Inflammation and its discontents: the role of cytokines in the pathophysiology of major depression
.
Biol. Psychiatry
,
65
,
732
741
.

11.

Hasler
,
G.
(
2010
)
Pathophysiology of depression: Do we have any solid evidence of interest to clinicians?
World Psychiatry
,
9
,
155
161
.

12.

Miller
,
A.H.
and
Raison
,
C.L.
(
2016
)
The role of inflammation in depression: from evolutionary imperative to modern treatment target
.
Nat. Rev. Immunol.
,
16
,
22
34
.

13.

Kendler
,
K.S.
,
Gatz
,
M.
,
Gardner
,
C.O.
and
Pedersen
,
N.L.
(
2006
)
A Swedish national twin study of lifetime major depression
.
Am. J. Psychiatry
,
163
,
109
114
.

14.

Levey
,
D.F.
,
Stein
,
M.B.
,
Wendt
,
F.R.
,
Pathak
,
G.A.
,
Zhou
,
H.
,
Aslan
,
M.
,
Quaden
,
R.
,
Harrington
,
K.M.
,
Nuñez
,
Y.Z.
,
Overstreet
,
C.
et al. (
2021
)
Bi-ancestral depression GWAS in the Million Veteran Program and meta-analysis in >1.2 million individuals highlight new therapeutic directions
.
Nat. Neurosci.
,
24
,
954
963
.

15.

Thorp
,
J.G.
,
Campos
,
A.I.
,
Grotzinger
,
A.D.
,
Gerring
,
Z.F.
,
An
,
J.
,
Ong
,
J.-S.
,
Wang
,
W.
,
23andMe Research Team
,
Shringarpure
,
S.
,
Byrne
,
E.M.
et al. (
2021
)
Symptom-level modelling unravels the shared genetic architecture of anxiety and depression
.
Nat. Hum. Behav
.,
5
,
1432
1442
.

16.

Cross-Disorder Group of the Psychiatric Genomics Consortium
Electronic address: [email protected] and Cross-Disorder Group of the Psychiatric Genomics Consortium
(
2019
)
Genomic relationships, novel loci, and pleiotropic mechanisms across eight psychiatric disorders
.
Cell
,
179
,
1469
1482.e11
.

17.

Howard
,
D.M.
,
Adams
,
M.J.
,
Clarke
,
T.-K.
,
Hafferty
,
J.D.
,
Gibson
,
J.
,
Shirali
,
M.
,
Coleman
,
J.R.I.
,
Hagenaars
,
S.P.
,
Ward
,
J.
,
Wigmore
,
E.M.
et al. (
2019
)
Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions
.
Nat. Neurosci.
,
22
,
343
352
.

18.

García-Marín
,
L.M.
,
Campos
,
A.I.
,
Cuéllar-Partida
,
G.
,
Medland
,
S.E.
,
Kollins
,
S.H.
and
Rentería
,
M.E.
(
2021
)
Large-scale genetic investigation reveals genetic liability to multiple complex traits influencing a higher risk of ADHD
.
Sci. Rep.
,
11
,
11
.

19.

van
Rheenen
,
W.
,
Peyrot
,
W.J.
,
Schork
,
A.J.
,
Lee
,
S.H.
and
Wray
,
N.R.
(
2019
)
Genetic correlations of polygenic disease traits: from theory to practice
.
Nat. Rev. Genet.
,
20
,
567
581
.

20.

García-Marín
,
L.M.
,
Campos
,
A.I.
,
Martin
,
N.G.
,
Cuéllar-Partida
,
G.
and
Rentería
,
M.E.
(
2021
)
Phenome-wide analysis highlights putative causal relationships between self-reported migraine and other complex traits
.
J. Headache Pain
,
22
,
66
.

21.

Pingault
,
J.-B.
,
O’Reilly
,
P.F.
,
Schoeler
,
T.
,
Ploubidis
,
G.B.
,
Rijsdijk
,
F.
and
Dudbridge
,
F.
(
2018
)
Using genetic data to strengthen causal inference in observational research
.
Nat. Rev. Genet.
,
19
,
566
580
.

22.

Davey Smith
,
G.
and
Hemani
,
G.
(
2014
)
Mendelian randomization: genetic anchors for causal inference in epidemiological studies
.
Hum. Mol. Genet.
,
23
,
R89
R98
.

23.

Verbanck
,
M.
,
Chen
,
C.-Y.
,
Neale
,
B.
and
Do
,
R.
(
2018
)
Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases
.
Nat. Genet.
,
50
,
693
698
.

24.

O’Connor
,
L.J.
and
Price
,
A.L.
(
2018
)
Distinguishing genetic correlation from causation across 52 diseases and complex traits
.
Nat. Genet.
,
50
,
1728
1734
.

25.

Choi
,
K.W.
,
Stein
,
M.B.
,
Nishimi
,
K.M.
,
Ge
,
T.
,
Coleman
,
J.R.I.
,
Chen
,
C.-Y.
,
Ratanatharathorn
,
A.
,
Zheutlin
,
A.B.
,
Dunn
,
E.C.
,
23andMe Research Team
et al. (
2020
)
An exposure-wide and Mendelian randomization approach to identifying modifiable factors for the prevention of depression
.
Am. J. Psychiatry
,
177
,
944
954
.

26.

Meng
,
X.
,
Brunet
,
A.
,
Turecki
,
G.
,
Liu
,
A.
,
D’Arcy
,
C.
and
Caron
,
J.
(
2017
)
Risk factor modifications and depression incidence: a 4-year longitudinal Canadian cohort of the Montreal Catchment Area Study
.
BMJ Open
,
7
,
e015156
.

27.

Mulugeta
,
A.
,
Zhou
,
A.
,
King
,
C.
and
Hyppönen
,
E.
(
2020
)
Association between major depressive disorder and multiple disease outcomes: a phenome-wide Mendelian randomisation study in the UK Biobank
.
Mol. Psychiatry
,
25
,
1469
1476
.

28.

Shen
,
X.
,
Howard
,
D.M.
,
Adams
,
M.J.
,
Hill
,
W.D.
,
Clarke
,
T.-K.
,
Major Depressive Disorder Working Group of the Psychiatric Genomics Consortium
,
Deary
,
I.J.
,
Whalley
,
H.C.
and
McIntosh
,
A.M.
(
2020
)
A phenome-wide association and Mendelian randomisation study of polygenic risk for depression in UK Biobank
.
Nat. Commun.
,
11
,
2301
.

29.

Shen
,
X.
,
Howard
,
D.
,
Adams
,
M.
,
Deary
,
I.
,
Whalley
,
H.
and
McIntosh
,
A.
(
2019
)
52 a phenome-wide association and Mendelian randomisation study of polygenic risk for depression in UK Biobank
.
Eur. Neuropsychopharmacol.
, (
2019
),
29
,
S88
.

30.

Kalin
,
N.H.
(
2020
)
The critical relationship between anxiety and depression
.
Am. J. Psychiatry
,
177
,
365
367
.

31.

Lamers
,
F.
,
van
Oppen
,
P.
,
Comijs
,
H.C.
,
Smit
,
J.H.
,
Spinhoven
,
P.
,
van
Balkom
,
A.J.L.M.
,
Nolen
,
W.A.
,
Zitman
,
F.G.
,
Beekman
,
A.T.F.
and
Penninx
,
B.W.J.H.
(
2011
)
Comorbidity patterns of anxiety and depressive disorders in a large cohort study: the Netherlands Study of Depression and Anxiety (NESDA)
.
J. Clin. Psychiatry
,
72
,
341
348
.

32.

Taporoski
,
T.P.
,
Negrão
,
A.B.
,
Horimoto
,
A.R.V.R.
,
Duarte
,
N.E.
,
Alvim
,
R.O.
,
de
Oliveira
,
C.M.
,
Krieger
,
J.E.
,
Von Schantz
,
M.
,
Vallada
,
H.
and
Pereira
,
A.C.
(
2015
)
Shared genetic factors of anxiety and depression symptoms in a Brazilian family-based cohort, the Baependi Heart Study
.
PLoS One
,
10
,
e0144255
.

33.

Davidson
,
K.M.
(
1995
)
Diagnosis of depression in alcohol dependence: changes in prevalence with drinking status
.
Br. J. Psychiatry
,
166
,
199
204
.

34.

Bonin
,
M.F.
,
McCreary
,
D.R.
and
Sadava
,
S.W.
(
2000
)
Problem drinking behavior in two community-based samples of adults: influence of gender, coping, loneliness, and depression
.
Psychol. Addict. Behav.
,
14
,
151
161
.

35.

Hämäläinen
,
J.
,
Kaprio
,
J.
,
Isometsä
,
E.
,
Heikkinen
,
M.
,
Poikolainen
,
K.
,
Lindeman
,
S.
and
Aro
,
H.
(
2001
)
Cigarette smoking, alcohol intoxication and major depressive episode in a representative population sample
.
J. Epidemiol. Community Health
,
55
,
573
576
.

36.

Gea
,
A.
,
Beunza
,
J.J.
,
Estruch
,
R.
,
Sánchez-Villegas
,
A.
,
Salas-Salvadó
,
J.
,
Buil-Cosiales
,
P.
,
Gómez-Gracia
,
E.
,
Covas
,
M.-I.
,
Corella
,
D.
,
Fiol
,
M.
et al. (
2013
)
Alcohol intake, wine consumption and the development of depression: the PREDIMED study
.
BMC Med.
,
11
,
192
.

37.

Kendler
,
K.S.
,
Heath
,
A.C.
,
Neale
,
M.C.
,
Kessler
,
R.C.
and
Eaves
,
L.J.
(
1993
)
Alcoholism and major depression in women. A twin study of the causes of comorbidity
.
Arch. Gen. Psychiatry
,
50
,
690
698
.

38.

Kovacevic
,
T.
,
Zaric
,
B.
,
Stanic
,
J.
,
Sakic
,
B.
and
Perin
,
B.
(
2014
)
Anxiety and depression in lung cancer patients—Are there any relations to clinico-pathological characteristics?
Eur. Respir. J.
,
44
,
44
.

39.

Yan
,
X.-R.
,
Chen
,
X.
and
Zhang
,
P.
(
2019
)
Prevalence and risk factors of depression in patients with lung cancer: protocol for a systematic review and meta-analysis
.
BMJ Open
,
9
,
e028994
.

40.

Hopwood
,
P.
and
Stephens
,
R.J.
(
2000
)
Depression in patients with lung cancer: prevalence and risk factors derived from quality-of-life data
.
J. Clin. Oncol.
,
18
,
893
903
.

41.

Lu
,
T.
,
Yang
,
X.
,
Huang
,
Y.
,
Zhao
,
M.
,
Li
,
M.
,
Ma
,
K.
,
Yin
,
J.
,
Zhan
,
C.
and
Wang
,
Q.
(
2019
)
Trends in the incidence, treatment, and survival of patients with lung cancer in the last four decades
.
Cancer Manag. Res.
,
11
,
943
953
.

42.

Sharp
,
L.
,
Carsin
,
A.-E.
and
Timmons
,
A.
(
2013
)
Associations between cancer-related financial stress and strain and psychological well-being among individuals living with cancer
.
Psycho-Oncology
,
22
,
745
755
.

43.

Milaneschi
,
Y.
,
Lamers
,
F.
,
Peyrot
,
W.J.
,
Baune
,
B.T.
,
Breen
,
G.
,
Dehghan
,
A.
,
Forstner
,
A.J.
,
Grabe
,
H.J.
,
Homuth
,
G.
,
Kan
,
C.
et al. (
2017
)
Genetic association of major depression with atypical features and obesity-related immunometabolic dysregulations
.
JAMA Psychiatry
,
74
,
1214
1225
.

44.

Weinberger
,
N.-A.
,
Kersting
,
A.
,
Riedel-Heller
,
S.G.
and
Luck-Sikorski
,
C.
(
2018
)
The relationship between weight status and depressive symptoms in a population sample with obesity: the mediating role of appearance evaluation
.
Obes. Facts
,
11
,
514
523
.

45.

Sarwer
,
D.B.
and
Polonsky
,
H.M.
(
2016
)
The psychosocial burden of obesity
.
Endocrinol. Metab. Clin. N. Am.
,
45
,
677
688
.

46.

Weinberger
,
N.-A.
,
Kersting
,
A.
,
Riedel-Heller
,
S.G.
and
Luck-Sikorski
,
C.
(
2016
)
Body dissatisfaction in individuals with obesity compared to normal-weight individuals: a systematic review and meta-analysis
.
Obes. Facts
,
9
,
424
441
.

47.

Grogan
,
S.
(
1999
)
Body image: understanding body dissatisfaction in men, women, and children
.
Choice (Middletown)
,
37
,
37–0628–37–0628
.

48.

Akhondzadeh
,
S.
(
2019
)
Depression and inflammation: Is there any role for biomarkers?
Avicenna J. Med. Biotechnol.
,
11
,
207
.

49.

García-Marín
,
L.M.
,
Campos
,
A.I.
,
Kho
,
P.-F.
,
Martin
,
N.G.
,
Cuéllar-Partida
,
G.
and
Rentería
,
M.E.
(
2021
)
Phenome-wide screening of GWAS data reveals the complex causal architecture of obesity
.
Hum. Genet.
,
140
,
1253
1265
.

50.

Ellulu
,
M.S.
,
Patimah
,
I.
,
Khaza’ai
,
H.
,
Rahmat
,
A.
and
Abed
,
Y.
(
2017
)
Obesity and inflammation: the linking mechanism and the complications
.
Arch. Med. Sci.
,
13
,
851
863
.

51.

Lehto
,
S.M.
,
Hintikka
,
J.
,
Niskanen
,
L.
,
Tolmunen
,
T.
,
Koivumaa-Honkanen
,
H.
,
Honkalampi
,
K.
and
Viinamäki
,
H.
(
2008
)
Low HDL cholesterol associates with major depression in a sample with a 7-year history of depressive symptoms
.
Prog. Neuro-Psychopharmacol. Biol. Psychiatry
,
32
,
1557
1561
.

52.

So
,
H.-C.
,
Chau
,
C.K.-L.
,
Cheng
,
Y.-Y.
and
Sham
,
P.C.
(
2020
)
Causal relationships between blood lipids and depression phenotypes: a Mendelian randomisation analysis
.
Psychol. Med.
,
14
,
1
13
.

53.

Shin
,
J.Y.
,
Suls
,
J.
and
Martin
,
R.
(
2008
)
Are cholesterol and depression inversely related? A meta-analysis of the association between two cardiac risk factors
.
Ann. Behav. Med.
,
36
,
33
43
.

54.

Parekh
,
A.
,
Smeeth
,
D.
,
Milner
,
Y.
and
Thure
,
S.
(
2017
)
The role of lipid biomarkers in major depression
.
Healthcare (Basel)
,
5
,
5
.

55.

Wang
,
H.
and
Peng
,
D.-Q.
(
2011
)
New insights into the mechanism of low high-density lipoprotein cholesterol in obesity
.
Lipids Health Dis.
,
10
,
176
.

56.

van
Reedt Dortland
,
A.K.B.
,
Giltay
,
E.J.
,
van
Veen
,
T.
,
van
Pelt
,
J.
,
Zitman
,
F.G.
and
Penninx
,
B.W.J.H.
(
2010
)
Associations between serum lipids and major depressive disorder: results from the Netherlands Study of Depression and Anxiety (NESDA)
.
J. Clin. Psychiatry
,
71
,
729
736
.

57.

Sheng
,
J.
,
Liu
,
S.
,
Wang
,
Y.
,
Cui
,
R.
and
Zhang
,
X.
(
2017
)
The link between depression and chronic pain: neural mechanisms in the brain
.
Neural Plast.
,
2017
,
9724371
.

58.

Campos
,
A.C.P.
,
Antunes
,
G.F.
,
Matsumoto
,
M.
,
Pagano
,
R.L.
and
Martinez
,
R.C.R.
(
2020
)
Neuroinflammation, pain and depression: an overview of the main findings
.
Front. Psychol.
,
11
,
1825
.

59.

Torta
,
R.
,
Ieraci
,
V.
and
Zizzi
,
F.
(
2017
)
A review of the emotional aspects of neuropathic pain: from comorbidity to co-pathogenesis
.
Pain Therapy
,
6
,
11
17
.

60.

Roughan
,
W.H.
,
Campos
,
A.I.
,
García-Marín
,
L.M.
,
Cuéllar-Partida
,
G.
,
Lupton
,
M.K.
,
Hickie
,
I.B.
,
Medland
,
S.E.
,
Wray
,
N.R.
,
Byrne
,
E.M.
,
Ngo
,
T.T.
et al. (
2021
)
Comorbid chronic pain and depression: shared risk factors and differential antidepressant effectiveness
.
Front. Psychiatry
,
0
,
1
13
.

61.

Campos
,
A.I.
,
Ngo
,
T.T.
,
Medland
,
S.E.
,
Wray
,
N.R.
,
Hickie
,
I.B.
,
Byrne
,
E.M.
,
Martin
,
N.G.
and
Rentería
,
M.E.
(
2021
)
Genetic risk for chronic pain is associated with lower antidepressant effectiveness: converging evidence for a depression subtype
.
Aust. N. Z. J. Psychiatry
,
48674211031491
.

62.

Crofford
,
L.J.
(
2015
)
Psychological aspects of chronic musculoskeletal pain
.
Best Pract. Res. Clin. Rheumatol.
,
29
,
147
155
.

63.

De Ruigh
,
A.
,
Roman
,
S.
,
Chen
,
J.
,
Pandolfino
,
J.E.
and
Kahrilas
,
P.J.
(
2014
)
Gaviscon double action liquid (antacid & alginate) is more effective than antacid in controlling post-prandial oesophageal acid exposure in GERD patients: a double-blind crossover study
.
Aliment. Pharmacol. Ther.
,
40
,
531
537
.

64.

Wagner
,
J.A.
and
Colombo
,
J.M.
(
2020
)
Medicine and media: the ranitidine debate
.
Clin. Transl. Sci.
,
13
,
649
.

65.

Al Quraan
,
A.M.
,
Beriwal
,
N.
,
Sangay
,
P.
and
Namgyal
,
T.
(
2019
)
The psychotic impact of Helicobacter pylori gastritis and functional dyspepsia on depression: a systematic review
.
Cureus
,
11
,
e5956
.

66.

Wu
,
Y.
,
Murray
,
G.K.
,
Byrne
,
E.M.
,
Sidorenko
,
J.
,
Visscher
,
P.M.
and
Wray
,
N.R.
(
2021
)
GWAS of peptic ulcer disease implicates Helicobacter pylori infection, other gastrointestinal disorders and depression
.
Nat. Commun.
,
12
,
1146
.

67.

Martin-Subero
,
M.
,
Anderson
,
G.
,
Kanchanatawan
,
B.
,
Berk
,
M.
and
Maes
,
M.
(
2016
)
Comorbidity between depression and inflammatory bowel disease explained by immune-inflammatory, oxidative, and nitrosative stress; tryptophan catabolite; and gut-brain pathways
.
CNS Spectr.
,
21
,
184
198
.

68.

Goodwin
,
G.M.
(
2006
)
Depression and associated physical diseases and symptoms
.
Dialogues Clin. Neurosci.
,
8
,
259
265
.

69.

Głąbska
,
D.
,
Kołota
,
A.
,
Lachowicz
,
K.
,
Skolmowska
,
D.
,
Stachoń
,
M.
and
Guzek
,
D.
(
2021
)
The influence of vitamin D intake and status on mental health in children: a systematic review
.
Nutrients
,
13
,
1
24
.

70.

Spedding
,
S.
(
2014
)
Vitamin D and depression: a systematic review and meta-analysis comparing studies with and without biological flaws
.
Nutrients
,
6
,
1501
1518
.

71.

Menon
,
V.
,
Kar
,
S.K.
,
Suthar
,
N.
and
Nebhinani
,
N.
(
2020
)
Vitamin D and depression: a critical appraisal of the evidence and future directions
.
Indian J. Psychol. Med.
,
42
,
11
21
.

72.

Revez
,
J.A.
,
Lin
,
T.
,
Qiao
,
Z.
,
Xue
,
A.
,
Holtz
,
Y.
,
Zhu
,
Z.
,
Zeng
,
J.
,
Wang
,
H.
,
Sidorenko
,
J.
,
Kemper
,
K.E.
et al. (
2020
)
Genome-wide association study identifies 143 loci associated with 25 hydroxyvitamin D concentration
.
Nat. Commun.
,
11
,
1647
.

73.

Gavin
,
A.R.
,
Walton
,
E.
,
Chae
,
D.H.
,
Alegria
,
M.
,
Jackson
,
J.S.
and
Takeuchi
,
D.
(
2010
)
The associations between socio-economic status and major depressive disorder among Blacks, Latinos, Asians and non-Hispanic Whites: findings from the Collaborative Psychiatric Epidemiology Studies
.
Psychol. Med.
,
40
,
51
61
.

74.

Freeman
,
A.
,
Tyrovolas
,
S.
,
Koyanagi
,
A.
,
Chatterji
,
S.
,
Leonardi
,
M.
,
Ayuso-Mateos
,
J.L.
,
Tobiasz-Adamczyk
,
B.
,
Koskinen
,
S.
,
Rummel-Kluge
,
C.
and
Haro
,
J.M.
(
2016
)
The role of socio-economic status in depression: results from the COURAGE (aging survey in Europe)
.
BMC Public Health
,
16
,
1098
.

75.

Clougherty
,
J.E.
,
Souza
,
K.
and
Cullen
,
M.R.
(
2010
)
Work and its role in shaping the social gradient in health
.
Ann. N. Y. Acad. Sci.
,
1186
,
102
124
.

76.

Burgard
,
S.A.
and
Lin
,
K.Y.
(
2013
)
Bad jobs, bad health? How work and working conditions contribute to health disparities
.
Am. Behav. Sci.
,
57
,
1105
1127
.

77.

Alkerwi
,
A.
,
‘A, Vernier
,
C.
,
Sauvageot
,
N.
,
Crichton
,
G.E.
and
Elias
,
M.F.
(
2015
)
Demographic and socioeconomic disparity in nutrition: application of a novel Correlated Component Regression approach
.
BMJ Open
,
5
,
e006814
.

78.

Głąbska
,
D.
,
Guzek
,
D.
,
Groele
,
B.
and
Gutkowska
,
K.
(
2020
)
Fruit and vegetable intake and mental health in adults: a systematic review
.
Nutrients
,
12
,
1
34
.

79.

Dharmayani
,
P.N.A.
,
Juergens
,
M.
,
Allman-Farinelli
,
M.
and
Mihrshahi
,
S.
(
2021
)
Association between fruit and vegetable consumption and depression symptoms in young people and adults aged 15–45: a systematic review of cohort studies
.
Int. J. Environ. Res. Public Health
,
18
,
780
.

80.

Cornish
,
S.
and
Mehl-Madrona
,
L.
(
2008
)
The role of vitamins and minerals in psychiatry
.
Integr. Med. Insights
,
3
,
33
42
.

81.

Irshad
,
M.
and
Chaudhuri
,
P.S.
(
2002
)
Oxidant-antioxidant system: role and significance in human body
.
Indian J. Exp. Biol.
,
40
,
1233
1239
.

82.

Maes
,
M.
,
Yirmyia
,
R.
,
Noraberg
,
J.
,
Brene
,
S.
,
Hibbeln
,
J.
,
Perini
,
G.
,
Kubera
,
M.
,
Bob
,
P.
,
Lerer
,
B.
and
Maj
,
M.
(
2009
)
The inflammatory & neurodegenerative (I&ND) hypothesis of depression: leads for future research and new drug developments in depression
.
Metab. Brain Dis.
,
24
,
27
53
.

83.

Maes
,
M.
,
Galecki
,
P.
,
Chang
,
Y.S.
and
Berk
,
M.
(
2011
)
A review on the oxidative and nitrosative stress (O&NS) pathways in major depression and their possible contribution to the (neuro)degenerative processes in that illness
.
Prog. Neuro-Psychopharmacol. Biol. Psychiatry
,
35
,
676
692
.

84.

Woo
,
J.-M.
and
Postolache
,
T.T.
(
2008
)
The impact of work environment on mood disorders and suicide: evidence and implications
.
Int. J. Disabil. Hum. Dev.
,
7
,
185
200
.

85.

Mounchetrou Njoya
,
I.
,
Paris
,
C.
,
Dinet
,
J.
,
Luc
,
A.
,
Lighezzolo-Alnot
,
J.
,
Pairon
,
J.-C.
and
Thaon
,
I.
(
2017
)
Anxious and depressive symptoms in the French Asbestos-Related Diseases Cohort: risk factors and self-perception of risk
.
Eur. J. Pub. Health
,
27
,
359
366
.

86.

Deaton
,
A.
and
Cartwright
,
N.
(
2018
)
Understanding and misunderstanding randomized controlled trials
.
Soc. Sci. Med.
,
210
,
2
21
.

87.

Kovesdy
,
C.P.
and
Kalantar-Zadeh
,
K.
(
2012
)
Observational studies versus randomized controlled trials: avenues to causal inference in nephrology
.
Adv. Chronic Kidney Dis.
,
19
,
11
18
.

88.

Boyko
,
E.J.
(
2013
)
Observational research — opportunities and limitations
.
J. Diabetes Complicat.
,
27
,
642
648
.

89.

Wu
,
Y.
,
Byrne
,
E.M.
,
Zheng
,
Z.
,
Kemper
,
K.E.
,
Yengo
,
L.
,
Mallett
,
A.J.
,
Yang
,
J.
,
Visscher
,
P.M.
and
Wray
,
N.R.
(
2019
)
Genome-wide association study of medication-use and associated disease in the UK Biobank
.
Nat. Commun.
,
10
,
1891
.

90.

Orsolini
,
L.
,
Latini
,
R.
,
Pompili
,
M.
,
Serafini
,
G.
,
Volpe
,
U.
,
Vellante
,
F.
,
Fornaro
,
M.
,
Valchera
,
A.
,
Tomasetti
,
C.
,
Fraticelli
,
S.
et al. (
2020
)
Understanding the complex of suicide in depression: from research to clinics
.
Psychiatry Investig.
,
17
,
207
221
.

91.

Zhang
,
J.
and
Li
,
Z.
(
2013
)
The association between depression and suicide when hopelessness is controlled for
.
Compr. Psychiatry
,
54
,
790
796
.

92.

Mullins
,
N.
,
Perroud
,
N.
,
Uher
,
R.
,
Butler
,
A.W.
,
Cohen-Woods
,
S.
,
Rivera
,
M.
,
Malki
,
K.
,
Euesden
,
J.
,
Power
,
R.A.
,
Tansey
,
K.E.
et al. (
2014
)
Genetic relationships between suicide attempts, suicidal ideation and major psychiatric disorders: a genome-wide association and polygenic scoring study
.
Am. J. Med. Genet. B Neuropsychiatr. Genet.
,
165B
,
428
437
.

93.

Albanese
,
B.J.
,
Macatee
,
R.J.
,
Stanley
,
I.H.
,
Bauer
,
B.W.
,
Capron
,
D.W.
,
Bernat
,
E.
,
Joiner
,
T.E.
and
Schmidt
,
N.B.
(
2019
)
Differentiating suicide attempts and suicidal ideation using neural markers of emotion regulation
.
J. Affect. Disord.
,
257
,
536
550
.

94.

Cuellar-Partida
,
G.
,
Lundberg
,
M.
,
Kho
,
P.F.
,
D’Urso
,
S.
,
Gutierrez-Mondragon
,
L.F.
and
Hwang
,
L.-D
.
Complex-Traits Genetics Virtual Lab: a community-driven web platform for post-GWAS analyses
. https://doi.org/10.1101/518027.

95.

Neale’s Lab
(
2018
)
GWAS Results
. (
accessed Aug 25, 2021
).

96.

García-Marín
,
L.M.
,
Campos
,
A.I.
,
Martin
,
N.G.
,
Cuéllar-Partida
,
G.
and
Rentería
,
M.E.
(
2021
)
Inference of causal relationships between sleep-related traits and 1,527 phenotypes using genetic data
.
Sleep
,
44
,
1
13
.

97.

Haworth
,
S.
,
Kho
,
P.F.
,
Holgerson
,
P.L.
,
Hwang
,
L.-D.
,
Timpson
,
N.J.
,
Rentería
,
M.E.
,
Johansson
,
I.
and
Cuellar-Partida
,
G.
(
2020
)
Assessment and visualization of phenome-wide causal relationships using genetic data: an application to dental caries and periodontitis
.
Eur. J. Hum. Genet.
,
29
,
300
308
.

98.

Bulik-Sullivan
,
B.K.
,
Loh
,
P.-R.
,
Finucane
,
H.K.
,
Ripke
,
S.
,
Yang
,
J.
,
Patterson
,
N.
,
Daly
,
M.J.
,
Price
,
A.L.
and
Neale
,
B.M.
(
2015
)
LD score regression distinguishes confounding from polygenicity in genome-wide association studies
.
Nat. Genet.
,
47
,
291
295
.

99.

Liu
,
M.
,
Jiang
,
Y.
,
Wedow
,
R.
,
Li
,
Y.
,
Brazel
,
D.M.
,
Chen
,
F.
,
Datta
,
G.
,
Davila-Velderrain
,
J.
,
McGuire
,
D.
,
Tian
,
C.
et al. (
2019
)
Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use
.
Nat. Genet.
,
51
,
237
244
.

100.

Purves
,
K.L.
,
Coleman
,
J.R.I.
,
Meier
,
S.M.
,
Rayner
,
C.
,
Davis
,
K.A.S.
,
Cheesman
,
R.
,
Bækvad-Hansen
,
M.
,
Børglum
,
A.D.
,
Wan
,
C.S.
,
Jürgen
,
D.J.
et al. (
2020
)
A major role for common genetic variation in anxiety disorders
.
Mol. Psychiatry
,
25
,
25
.

101.

Yengo
,
L.
,
Sidorenko
,
J.
,
Kemper
,
K.E.
,
Zheng
,
Z.
,
Wood
,
A.R.
,
Weedon
,
M.N.
,
Frayling
,
T.M.
,
Hirschhorn
,
J.
,
Yang
,
J.
,
Visscher
,
P.M.
et al. (
2018
)
Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of European ancestry
.
Hum. Mol. Genet.
,
27
,
3641
.

Author notes

Present address: 23andMe, Inc., Sunnyvale, CA, USA.

Asma M. Aman and Luis M. García-Marín authors contributed equally

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model)

Supplementary data