- 1Department of Psychiatry, First Hospital/First Clinical Medical College of Shanxi Medical University, Taiyuan, China
- 2School of Systems Biology, George Mason University, Fairfax, VA, United States
- 3Institute of Systems Medicine, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, China
- 4Suzhou Institute of Systems Medicine, Suzhou, China
- 5Research Centre for Medical Genetics, Moscow, Russia
- 6Department of Psychiatry, The Affiliated Brain Hospital of Nanjing Medical University, Nanjing, China
- 7Institute of Neuropsychiatry, The Affiliated Brain Hospital of Nanjing Medical University, Nanjing, China
Objectives: Deciphering the genetic relationships between major depressive disorder (MDD) and atopic diseases (asthma, hay fever, and eczema) may facilitate understanding of their biological mechanisms as well as the development of novel treatment regimens. Here we tested the genetic correlation between MDD and atopic diseases by linkage disequilibrium score regression.
Methods: A polygenic overlap analysis was performed to estimate shared genetic variations between the two diseases. Causal relationships between MDD and atopic diseases were investigated using two-sample bidirectional Mendelian randomization analysis. Genomic loci shared between MDD and atopic diseases were identified using cross-trait meta-analysis. Putative functional genes were evaluated by fine-mapping of transcriptome-wide associations.
Results: The polygenic analysis revealed approximately 15.8 thousand variants causally influencing MDD and 0.9 thousand variants influencing atopic diseases. Among these variants, approximately 0.8 thousand were shared between the two diseases. Mendelian randomization analysis indicates that genetic liability to MDD has a causal effect on atopic diseases (b = 0.22, p = 1.76 × 10-6), while genetic liability to atopic diseases confers a weak causal effect on MDD (b = 0.05, p = 7.57 × 10-3). Cross-trait meta-analyses of MDD and atopic diseases identified 18 shared genomic loci. Both fine-mapping of transcriptome-wide associations and analysis of existing literature suggest the estrogen receptor β-encoding gene ESR2 as one of the potential risk factors for both MDD and atopic diseases.
Conclusion: Our findings reveal shared genetic liability and causal links between MDD and atopic diseases, which shed light on the phenotypic relationship between MDD and atopic diseases.
Introduction
Mental disorders confer a heavy burden on society (1). Major depressive disorder (MDD), the most prevalent mental disorder accompanied by considerable morbidity, mortality, and risk of suicide, is characterized by persistent low mood (2). MDD and depressive symptoms have close associations with certain physical conditions. Generally speaking, long-term depression adds to the risk for somatic illness, and, vice versa, chronic somatic diseases are frequently accompanied by depression (3). When comorbid with other ailments, for example, atopic diseases (ADs), MDD produces worse clinical outcomes and incurs higher healthcare costs.
ADs are driven by the dysfunction of the immune system. Three kinds of common ADs, namely, asthma, hay fever (allergic rhinitis), and eczema (atopic dermatitis), may coexist in the same individuals (4). Asthma, a chronic airway disease that is common worldwide, is characterized by coughing, wheezing, shortness of breath, and/or chest tightness due to increased airway reactivity, inflammation, and/or mucus production. In 2015, asthma affected 358 million people globally and caused about 400,000 deaths (5). Allergic rhinitis is an inflammatory disease characterized by nasal congestion, rhinorrhea, sneezing, and/or nasal itching. Allergic rhinitis is one of the most common diseases in adults (20%~30%), and the most common chronic disease in children (up to 40%) in the United States (6). Eczema is an inflammatory skin disease that is caused by a dysfunction of a skin barrier followed by aberrant inflammation/immune responses; this disease is affecting 5% of the population worldwide (7). Together, symptoms of ADs significantly impair quality of life and impose a heavy cost on society. Common comorbidities of MDD with ADs have been documented previously (8–12). Specifically, allergic rhinitis has been shown to have a positive association with MDD (odds ratio: 1.24) (8). In patients with asthma, the hazard ratio of MDD increases by 35%, and MDD patients show about 25% increased hazard ratio for being affected by asthma (9). Atopic eczema is also associated with an increased incidence of new depression (hazard ratio: 1.14) (10).
Although previous studies have detected associations between MDD and ADs, several key questions remain pending: 1) to what extent may the two conditions share genetic components? 2) Are the phenotypic associations mediated by genetic variations? 3) What molecular and cellular mechanisms underline these associations?
Genetic relationships between two traits are commonly quantified by genetic correlation coefficients. The sign of the correlation coefficient indicates directions of the shared genetic effects. When dealing with mixtures of effect directions across shared genetic variants, genetic correlation analyses may be underpowered (13). A polygenic overlap was recently proposed to measure the fraction of genetic variants causally associated with both traits over the total number of causal variants across a pair of traits involved (13).
Mendelian randomization (MR) is an analytic framework that utilizes genetic variants as instrumental variables to test for causative association between an exposure and an outcome (14). Recently, a general type of SMR (GSMR) had been developed by leveraging power from multiple genetic variants to account for linkage disequilibrium (LD) between the variants (15).
Recently, Zhu et al. reported a causal effect of MDD on asthma and identified 10 loci shared by asthma and MDD by cross-trait meta-analysis (16). The GWAS dataset for MDD, however, did not include the 23andMe samples. We set on taking this line of investigation further, by both utilizing a larger MDD dataset and including two other ADs related to asthma, namely, allergic rhinitis and atopic dermatitis. Asthma, allergic rhinitis, and atopic dermatitis genetically correlate with each other and are often comorbid (17). The genetic liability to MDD may confer a causal effect on all of these ADs. Dissection of this shared genetic liability may deliver novel insights into the pathophysiology of both MDD and ADs.
Methods
GWAS Summary Datasets and Quality Control
This study relied on both de-identified publicly available summary-level GWAS data and the pre-approval 23andMe dataset. The resultant MDD dataset included 135,458 cases and 344,901 healthy controls (18), and the AD dataset included 96,794 cases and 145,775 healthy controls (19). For the inclusion of each dataset, both bi-allelic SNPs and imputation INFO above 0.80 were required. Each SNP was compared between the two datasets, and SNPs with conflicting alleles were excluded. If an SNP was mapped to opposite strands in the two datasets, alleles of this SNP in the second dataset were flipped, and the effect direction was reversed.
Genetic Correlation and Polygenic Overlap Analysis
GWAS summary results were utilized to analyze the genetic correlation of MDD with ADs by LD score regression software (LDSC, v1.0.1) (20, 21). A polygenic overlap was analyzed by MiXeR v1.2 using default parameters (13). Using GWAS summary statistics, MiXeR quantifies the polygenic overlap irrespective of the genetic correlation between traits. Based on the univariate causal mixture model (22), MiXeR builds four bivariate normal distributions, with two causal components for variants specific to each trait, one causal component for variants affecting both traits, and a null component for variants with no effect on either trait. The likelihood function of the observed signed test statistics (GWAS Z-scores) is produced from the prior distribution of genetic effects, incorporating effects of the LD structure, sample size, minor allele frequency (MAF), cryptic relationships, and sample overlap. The summary statistics are used to estimate the parameters of the mixture model by optimization of the likelihood function. The number of causal variants reported by the software is 22.6% of the total estimated variants, which account for 90% of SNP heritability for each trait.
MR Analysis
Bidirectional causal associations between MDD and ADs were inferred using GSMR v1.0.9 (15). Instrumental variants were selected based on default p ≤ 5×10-8. It is well accepted that pleiotropy is a potential source of bias and an inflated estimation in an MR analysis (23). In GSMR, the HEIDI-outlier statistical approach allows the detection and elimination of genetic instruments with apparent pleiotropic effects on both risk factors and disease (15, 24). It was suggested that genetic correlation may confound Mendelian randomization estimates (25). To examine this possibility, we performed a latent causal variable model (LCV) analysis between MDD and ADs (26). The LCV framework utilizes the genetic causality proportion (GCP) to quantify the partial causality of trait 1 on trait 2. The GCP ranges from 0 (no partial genetic causality) to 1 (full genetic causality). A high value of GCP indicates a causal effect of interventions targeting trait 1 on trait 2.
Cross-Trait Meta-Analysis
A cross-trait meta-analysis of the MDD and the ADs was executed by the subset-based fixed-effect method ASSET v2.4.0, which permits the characterization of each SNP with respect to its pattern of effects on multiple phenotypes (27). For each assessed variant, this type of analysis returns a p-value for the best subset containing the studies contributing to the overall association signal. The meta-analysis pools the effect of a given SNP across K studies, weighting the effects by the size of the respective study. After subset-based meta-analysis, SNP-related findings were considered statistically significant, if two-tailed p values were lower than 5 × 10-8. In the meta-analysis results, functional annotation and gene-mapping of variants and identifying LD-independent genomic regions were performed on a FUMA platform (28). Firstly, independent significant SNPs (IndSigSNPs) were identified based on their p-value being genome-wide significant (p ≤ 5.0 × 10−8) and independent of each other (r2 < 0.6). Secondly, lead SNPs were identified as a subset of the independent significant SNPs that were in LD with each other at r2 < 0.1 within a 250-kb window. The gene-based association for the meta-analysis of MDD and ADs was conducted using MAGMA (29).
To ensure that sample overlap did inflate estimates of genetic overlap between MDD and ADs, λmeta statistics, which use effect size concordance to detect sample overlap or heterogeneity, were calculated (30). Under the null hypothesis, λmeta equals 1 when the pair of cohorts are completely independent. When there are overlapping samples, λmeta is less than 1.
Fine-Mapping of TWAS Associations
To prioritize putatively causal genes, fine-mapping of causal gene sets (FOCUS v0.6.10) (31) to the meta-analysis result of MDD and ADs was performed in four relevant tissues, including the brain, whole blood, lung, and skin. Using FOCUS, predicted expression correlations were modeled and posterior inclusion probabilities (PIP) are assigned to genes within each transcriptome-wide association study (TWAS) region in the relevant tissue types. A multi-tissue eQTL reference weight database from the software was used as eQTL weights, while LD information from LDSC was used as a reference. Multiple-testing correction was used to account for all gene–tissue pairs based on Benjamini–Hochberg adjusted TWAS p-values (FDR < 0.05).
Knowledge-Based Analysis
GWAS results, including meta-analysis, were obtained for depression (major depressive disorder and depressive symptoms) and for ADs from the GWAS Catalog database (access date: April 17, 2020) (32). We explore whether the genes shared by MDD and ADs have been identified in previous genome-wide association studies. Protein–protein interaction analysis was conducted using STRING v11 (33). Enrichment of the 27 genes in the GWAS catalog reported genes was analyzed using FUMA (28).
All the statistical analyses were conducted in R 3.6.1 or Python 3.7 environment. A detailed description of the methods is provided in the Supplementary File.
Results
Genetic Correlation and Polygenic Overlap Analysis
MDD displayed a significant genetic correlation with ADs (r = 0.18, s.e. = 0.03, p = 1.04 × 10-9). The LD score intercept did not deviate from zero (0.017). The polygenic analysis highlighted approximately 15.8 thousand variants causally influencing MDD and 0.9 thousand variants influencing ADs. Among these variants, approximately 0.8 thousand variants were shared between the two diseases (Figure 1A). MDD has much larger numbers of causal variants than ADs, indicating a higher polygenic property of MDD.
Figure 1 Shared causal variants and causal associations between MDD and ADs. (A) Venn diagrams of unique and shared polygenic components at the causal level, showing a polygenic overlap between MDD and ADs. The numbers indicate the estimated quantity of causal variants (in thousands) per component, explaining 90% of SNP heritability in each phenotype. The size of the circles reflects the degree of polygenicity. (B) Causal associations between MDD and ADs. The lines denote effect sizes (B). The left panel denotes the causal effect of MDD on ADs. The left panel denotes the causal effect of ADs on MDD.
MR Analysis
Mendelian randomization analysis indicated that genetic liability to MDD has a causal effect on ADs (b = 0.22, s.e. = 0.05, OR = 1.25, 95%CI: 1.13–1.37, p = 1.76 × 10-6), with 45 independent instrumental variants being involved. The genetic liability of ADs conferred a causal effect on MDD (b = 0.05, s.e. = 0.02, OR = 1.05, 95%CI: 1.01–1.09, p = 7.57 × 10-3), with 115 independent instrumental variants being involved (Figure 1B). The LCV analysis showed that GCP was 0.49 (0.32), supporting a causal effect of genetic liability to MDD on ADs.
Cross-Trait Meta-Analysis
The cross-trait meta-analysis of MDD and ADs revealed the involvement of 103 loci, 470 significant independent SNPs (IndSigSNPs), and 141 lead SNPs, including 44 pleiotropic IndSigSNPs located in 18 loci (associated with both traits) (Figure 2A, Table 1 and Supplementary Tables 1, 2). The 14q23 locus is shown in Figure 2B. A total of 82 pleiotropic protein-coding genes were identified, including 27 protein-coding genes implicated by the pleiotropic IndSigSNPs and another 55 protein-coding genes implicated by SNPs tagged by IndSigSNPs (Supplementary Table 3). The gene-based association for the meta-analysis of MDD and ADs identified a total of 273 significant genes at the threshold of 2.70 × 10-6 (Bonferroni correction, 0.05/18,545) (Supplementary Table 4). Compared with SNP-based analysis, an additional 63 genes were identified by the gene-based analysis, including DRD2. The λmeta value was at 1.18 for datasets between MDD and ADs, indicating no significant overlap between MDD and AD GWAS samples. Quantile–quantile (QQ) plots to display the observed meta-analysis statistics versus the expected statistics under the null model of no associations in the -log10(p) scale are shown in Supplementary Figure 1.
Figure 2 Cross-trait meta-analysis of MDD and ADs. (A) Manhattan plot of meta-analysis of MDD with ADs. The x-axis is the chromosomal position of SNPs, and the y-axis is the significance of the SNPs (-log10P). Protein-coding genes containing or adjacent to independent significant SNPs shared by two traits were annotated. PCDHA1_6: PCDHA1, PCDHA2, PCDHA3, PCDHA4, PCDHA5, and PCDHA6. (B) The 14q23 locus containing the ESR2 gene. Each SNP is colored based on the highest r2 to one of the independent significant SNPs. (C) Enrichment of the 27 protein-coding genes in GWAS catalog gene sets.
Fine-Mapping of TWAS Associations
To prioritize putatively causal genes, we used the fine-mapping of TWAS associations. A total of 126 gene–tissue pairs were identified between the 82 genes and the four tissues, with 36 genes being associated with two or more tissues (Supplementary Table 5). A total of 31 gene–tissue pairs were in the credible sets. Fifteen genes associated with three or more tissues are listed in Table 2. However, most genes in Table 2 had low PIP. Of note, the ESR2 gene was associated with three tissues (skin, lung, and blood) with relatively high posterior probability (Figure 3).
Figure 3 Transcriptome-wide association study of the meta-analysis of MDD and ADs. (A) skin; (B) blood; (C) lung; (D) brain. Within each panel, the top part is the transcriptome-wide association signal indicating strength of the predicted expression association with trait, and the bottom part is the induced correlation of the predicted expression.
Knowledge-Based Analysis
A total of 23 out of the 27 pleiotropic protein-coding genes have been identified in previous GWASs on depression or ADs (Supplementary Table 6). Among these 23 genes were 16 genome-wide risk genes for depression, including BEND4, DENND1A, ESR2, L3MBTL2, NEGR1, PCDHA1, PCDHA2, PCDHA3, PCDHA4, PCDHA5, PCDHA6, RBX1, SLC30A9, SORCS3, SYNE2, and TMEM161B, and 8 genome-wide risk genes for ADs, including BOLL, DIAPH1, GLB1, MOB4, NEK6, OVOL1, RANGAP1, and RBX1. Enrichment of the 27 genes in the GWAS catalog-reported genes revealed that these genes were enriched in several mental disorders and basophil neutrophil counts, as well as neutrophil counts, supporting the involvement of these genes in neurodevelopmental conditions and atopic diseases (Figure 2C and Supplementary Table 7).
PPI analysis showed that a majority of the 82 genes are interconnected, forming one large network and several small networks (Supplementary Figure 2). Schematics of ESR2 gene interactions with depression and ADs are shown in Figure 4.
Discussion
In this study, we detected a significant genetic correlation between MDD and ADs (r = 0.18), at a level comparable to that for a previously reported correlation of MDD and autism spectrum disorder (r = 0.16) (34). Our results indicate a much higher polygenicity of MDD when compared to ADs, with substantial polygenic overlap between these conditions identified. Nearly 90% of causal variants influencing the risk of ADs may also affect MDD. Cumulative evidence supports a close relationship between these two conditions in the context of underlining genetics.
More importantly, causal relationships between MDD and ADs were discovered. In particular, a major causal effect of genetic liability to depression on ADs was detected. Although liability to ADs also exerts a statistically significant causal effect on MDD, the size of this effect is relatively small (b = 0.05). Previous studies already showed the possible influence exerted by MDD on ADs. For instance, patients with MDD show elevated levels of non-esterified fatty acids in plasma (35); other studies showed that fatty acids may contribute to the development of atopic diseases such as hay fever and asthma (36). Elevated serum interferon levels may contribute to eczema and also are commonly detected in MDD (37). Moreover, MDD has been shown to stimulate the production of cytokines (38), including IL-13 and IL-6, both of which are also strongly involved in asthma pathogenesis (39). Our findings are consistent with these previous studies and partially explain the previously reported comorbidity of MDD and ADs (8–12), while adding novel insights into underlying pathogenetic mechanisms. Notably, one previous study reported that depression may lead to asthma rather than the opposite (40). The causal effect of ADs on MDD should be further evaluated in additional datasets.
Shared genetic liability between MDD and ADs offers the possibility of employing polygenic risk scores (PRS) for evaluating allergic risks in MDD patients and the risk of developing depression in AD patients. This strategy may lead to an improvement in the clinical management of these conditions. Shared biological markers of MDD and ADs are far from being well studied. The cross-trait analysis revealed that MDD and ADs share 18 loci and a panel of protein-coding genes. The majority of these pleiotropic protein-coding genes have been previously implicated either in depression or in ADs, with a genome-wide significance level. For example, the RBX1 gene was reported as a significant contributor to both depression (41) and ADs (42). To shed new light on the genetic susceptibility of ADs and MDD, we have concentrated on the estrogen receptor β encoding gene ESR2 for further discussion.
Estrogen is capable of modulating neurotransmitter turnover to enhance the levels of serotonin and noradrenaline and participates in the regulation of serotonin receptor amounts and function (43). Accumulating evidence indicates the involvement of estrogen signaling in depression (44). In females, estrogen fluctuations are associated with depressed mood (45), and the beneficial effects of estrogen-containing hormone treatments were reported (46, 47). The gene for estrogen receptor β, ESR2, has been previously identified as a genome-wide significant gene contributing to MDD (18, 48). As the levels of estrogen are easily modulated by pharmacological means, the association between ESR2 and MDD may inform the development of personalized treatment modalities for this condition. Notably, model studies in neonatal rats treated with antidepressant clomipramine uncovered both the changes in the levels of estrogen receptors on the surface of brain cells and the neurochemical changes that resemble human depression (49). The role of estrogens in the development of ADs is noticeable as well. Women have a higher prevalence of asthma and display its greater severity than men (50). Estrogen receptors are found on numerous immune-regulatory cells, with estrogen-dependent responses favoring the shift toward allergy. In particular, estrogens promote allergic response by stimulating Th2 polarization, boosting class switching of B cells to IgE production, and prompting mast cell and basophil degranulation (51). ESR2 and its product, estrogen receptor β, have been suggested as potential targets for asthma treatment (52). There is also accumulating evidence supporting estrogens’ role in hay fever and eczema (53, 54). In particular, there is a correlation between the mean number of ER-β-positive cells in the nasal mucosa and seasonal allergy symptoms (55).
This study identified ESR2 as a novel genome-wide significant contributor to ADs, providing strong support for the involvement of the estrogen pathway in ADs. Fine-mapping of TWAS had assigned the posterior probability for causality for ESR2 in the skin, blood PBMCs, and lung tissue at 0.998, 0.256, and 0.243, respectively. Although the fine-mapping of TWAS hits did not support the involvement of ESR2 in the brain, analysis of existing literature points at its role in neurodevelopment and mental disorders. Together, our findings highlight ESR2 as a critical gene for both MDD and ADs and point to its relevance at the therapy target.
The presented study has several strengths. First, we utilized the largest combination of available datasets as a study backbone. Furthermore, to avoid potential population heterogeneity across the studies, we limited our analysis to individuals of European ancestry. Lastly, the genetic relationship between MDD and ADs was explored systemically by employing multiple analytic frameworks.
However, several limitations should also be noted. The datasets employed in this study only contained data of three subtypes of ADs. Further studies using more datasets covering other subtypes of ADs are warranted to evaluate the associations between MDD and ADs. In TWAS, the gene expression levels are imputed from weighted linear combinations of SNPs and, therefore, may report noise. As our analysis was limited to a genetic component of each trait, hence, the presented results should be interpreted cautiously, with the understanding that human traits arise from a complex web of interactions of various psycho-social-environmental factors.
In summary, our findings reveal shared genetic liability and causal links between MDD and atopic diseases, which may underline the phenotypic relationship between MDD and ADs. Presented results may have implications both for the therapy and for the management of MDD and ADs.
Data Availability Statement
Publicly available datasets were analyzed in this study. These data can be found here: https://www.ebi.ac.uk/gwas.
Author Contributions
FZ contributed to the study design and data analysis. HC, SL, and AB contributed to drafting and revising the work. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (81471364).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank members of the Psychiatric Genomics Consortium, the UK Biobank, 23andMe, and other teams, who generously shared the GWAS data.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.665160/full#supplementary-material
References
1. DALYs GBD, Collaborators H. Global, Regional, and National Disability-Adjusted Life-Years (DALYs) for 359 Diseases and Injuries and Healthy Life Expectancy (HALE) for 195 Countries and Territories, 1990-2017: A Systematic Analysis for the Global Burden of Disease Study 2017. Lancet (2018) 392(10159):1859–922. doi: 10.1016/S0140-6736(18)32335-3
2. Ferrari AJ, Charlson FJ, Norman RE, Patten SB, Freedman G, Murray CJ, et al. Burden of Depressive Disorders by Country, Sex, Age, and Year: Findings From the Global Burden of Disease Study 2010. PloS Med (2013) 10(11):e1001547. doi: 10.1371/journal.pmed.1001547
3. Veerman JL, Dowrick C, Ayuso-Mateos JL, Dunn G, Barendregt JJ. Population Prevalence of Depression and Mean Beck Depression Inventory Score. Br J Psychiatry (2009) 195(6):516–9. doi: 10.1192/bjp.bp.109.066191
4. Pinart M, Benet M, Annesi-Maesano I, von Berg A, Berdel D, Carlsen KC, et al. Comorbidity of Eczema, Rhinitis, and Asthma in IgE-Sensitised and Non-IgE-Sensitised Children in MeDALL: A Population-Based Cohort Study. Lancet Respir Med (2014) 2(2):131–40. doi: 10.1016/S2213-2600(13)70277-7
5. Disease GBD, Injury I, Prevalence C. Global, Regional, and National Incidence, Prevalence, and Years Lived With Disability for 310 Diseases and Injuries, 1990-2015: A Systematic Analysis for the Global Burden of Disease Study 2015. Lancet (2016) 388(10053):1545–602. doi: 10.1016/S0140-6736(16)31678-6
6. Mattos JL, Woodard CR, Payne SC. Trends in Common Rhinologic Illnesses: Analysis of U.S. Healthcare Surveys 1995-2007. Int Forum Allergy Rhinol (2011) 1(1):3–12. doi: 10.1002/alr.20003
7. Barbarot S, Auziere S, Gadkari A, Girolomoni G, Puig L, Simpson EL, et al. Epidemiology of Atopic Dermatitis in Adults: Results From an International Survey. Allergy (2018) 73(6):1284–93. doi: 10.1111/all.13401
8. Misery L, Taieb C, Schollhammer M, Bertolus S, Coulibaly E, Feton-Danou N, et al. Psychological Consequences of the Most Common Dermatoses: Data From the Objectifs Peau Study. Acta Derm Venereol (2020) 100(13):adv00175. doi: 10.2340/00015555-3531
9. Choi HG, Kim JH, Park JY, Hwang YI, Jang SH, Jung KS. Association Between Asthma and Depression: A National Cohort Study. J Allergy Clin Immunol Pract (2019) 7(4):1239–45.e1231. doi: 10.1016/j.jaip.2018.10.046
10. Schonmann Y, Mansfield KE, Hayes JF, Abuabara K, Roberts A, Smeeth L, et al. Atopic Eczema in Adulthood and Risk of Depression and Anxiety: A Population-Based Cohort Study. J Allergy Clin Immunol Pract (2020) 8(1):248–57.e216. doi: 10.1016/j.jaip.2019.08.030
11. Cheng CM, Hsu JW, Huang KL, Bai YM, Su TP, Li CT, et al. Risk of Developing Major Depressive Disorder and Anxiety Disorders Among Adolescents and Adults With Atopic Dermatitis: A Nationwide Longitudinal Study. J Affect Disord (2015) 178:60–5. doi: 10.1016/j.jad.2015.02.025
12. Wan J, Takeshita J, Shin DB, Gelfand JM. Mental Health Impairment Among Children With Atopic Dermatitis: A United States Population-Based Cross-Sectional Study of the 2013-2017 National Health Interview Survey. J Am Acad Dermatol (2020) 82(6):1368–75. doi: 10.1016/j.jaad.2019.10.019
13. Frei O, Holland D, Smeland OB, Shadrin AA, Fan CC, Maeland S, et al. Bivariate Causal Mixture Model Quantifies Polygenic Overlap Between Complex Traits Beyond Genetic Correlation. Nat Commun (2019) 10(1):2417. doi: 10.1038/s41467-019-10310-0
14. Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian Randomization: Using Genes as Instruments for Making Causal Inferences in Epidemiology. Stat Med (2008) 27(8):1133–63. doi: 10.1002/sim.3034
15. Zhu Z, Zheng Z, Zhang F, Wu Y, Trzaskowski M, Maier R, et al. Causal Associations Between Risk Factors and Common Diseases Inferred From GWAS Summary Data. Nat Commun (2018) 9(1):224. doi: 10.1038/s41467-017-02317-2
16. Zhu Z, Zhu X, Liu CL, Shi H, Shen S, Yang Y, et al. Shared Genetics of Asthma and Mental Health Disorders: A Large-Scale Genome-Wide Cross-Trait Analysis. Eur Respir J (2019) 54(6):1901507. doi: 10.1183/13993003.01507-2019
17. Zhu Z, Lee PH, Chaffin MD, Chung W, Loh PR, Lu Q, et al. A Genome-Wide Cross-Trait Analysis From UK Biobank Highlights the Shared Genetic Architecture of Asthma and Allergic Diseases. Nat Genet (2018) 50(6):857–64. doi: 10.1038/s41588-018-0121-0
18. Wray NR, Ripke S, Mattheisen M, Trzaskowski M, Byrne EM, Abdellaoui A, et al. Genome-Wide Association Analyses Identify 44 Risk Variants and Refine the Genetic Architecture of Major Depression. Nat Genet (2018) 50(5):668–81. doi: 10.1038/s41588-018-0090-3
19. Ferreira MA, Vonk JM, Baurecht H, Marenholz I, Tian C, Hoffman JD, et al. Shared Genetic Origin of Asthma, Hay Fever and Eczema Elucidates Allergic Disease Biology. Nat Genet (2017) 49(12):1752–7. doi: 10.1038/ng.3985
20. Bulik-Sullivan BK, Loh PR, Finucane HK, Ripke S, Yang J, Schizophrenia Working Group of the Psychiatric Genomics C, et al. LD Score Regression Distinguishes Confounding From Polygenicity in Genome-Wide Association Studies. Nat Genet (2015) 47(3):291–5. doi: 10.1038/ng.3211
21. Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh PR, et al. An Atlas of Genetic Correlations Across Human Diseases and Traits. Nat Genet (2015) 47(11):1236–41. doi: 10.1038/ng.3406
22. Holland D, Frei O, Desikan R, Fan CC, Shadrin AA, Smeland OB, et al. Beyond SNP Heritability: Polygenicity and Discoverability of Phenotypes Estimated With a Univariate Gaussian Mixture Model. PloS Genet (2020) 16(5):e1008612. doi: 10.1371/journal.pgen.1008612
23. Solovieff N, Cotsapas C, Lee PH, Purcell SM, Smoller JW. Pleiotropy in Complex Traits: Challenges and Strategies. Nat Rev Genet (2013) 14(7):483–95. doi: 10.1038/nrg3461
24. Ong JS, MacGregor S. Implementing MR-PRESSO and GCTA-GSMR for Pleiotropy Assessment in Mendelian Randomization Studies From a Practitioner’s Perspective. Genet Epidemiol (2019) 43(6):609–16. doi: 10.1002/gepi.22207
25. Verbanck M, Chen CY, Neale B, Do R. Detection of Widespread Horizontal Pleiotropy in Causal Relationships Inferred From Mendelian Randomization Between Complex Traits and Diseases. Nat Genet (2018) 50(5):693–8. doi: 10.1038/s41588-018-0099-7
26. O’Connor LJ, Price AL. Distinguishing Genetic Correlation From Causation Across 52 Diseases and Complex Traits. Nat Genet (2018) 50(12):1728–34. doi: 10.1038/s41588-018-0255-0
27. Bhattacharjee S, Rajaraman P, Jacobs KB, Wheeler WA, Melin BS, Hartge P, et al. A Subset-Based Approach Improves Power and Interpretation for the Combined Analysis of Genetic Association Studies of Heterogeneous Traits. Am J Hum Genet (2012) 90(5):821–35. doi: 10.1016/j.ajhg.2012.03.015
28. Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional Mapping and Annotation of Genetic Associations With FUMA. Nat Commun (2017) 8(1):1826. doi: 10.1038/s41467-017-01261-5
29. de Leeuw CA, Mooij JM, Heskes T, Posthuma D. MAGMA: Generalized Gene-Set Analysis of GWAS Data. PloS Comput Biol (2015) 11(4):e1004219. doi: 10.1371/journal.pcbi.1004219
30. Chen GB, Lee SH, Zhu ZX, Benyamin B, Robinson MR. EigenGWAS: Finding Loci Under Selection Through Genome-Wide Association Studies of Eigenvectors in Structured Populations. Hered (Edinb) (2016) 117(1):51–61. doi: 10.1038/hdy.2016.25
31. Mancuso N, Freund MK, Johnson R, Shi H, Kichaev G, Gusev A, et al. Probabilistic Fine-Mapping of Transcriptome-Wide Association Studies. Nat Genet (2019) 51(4):675–82. doi: 10.1038/s41588-019-0367-1
32. Buniello A, MacArthur JAL, Cerezo M, Harris LW, Hayhurst J, Malangone C, et al. The NHGRI-EBI GWAS Catalog of Published Genome-Wide Association Studies, Targeted Arrays and Summary Statistics 2019. Nucleic Acids Res (2019) 47(D1):D1005–12. doi: 10.1093/nar/gky1120
33. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING V11: Protein-Protein Association Networks With Increased Coverage, Supporting Functional Discovery in Genome-Wide Experimental Datasets. Nucleic Acids Res (2019) 47(D1):D607–13. doi: 10.1093/nar/gky1131
34. Brainstorm C, Anttila V, Bulik-Sullivan B, Finucane HK, Walters RK, Bras J, et al. Analysis of Shared Heritability in Common Disorders of the Brain. Science (2018) 360(6395):eaap8757. doi: 10.1126/science.aap8757
35. Badawy AA. Tryptophan: The Key to Boosting Brain Serotonin Synthesis in Depressive Illness. J Psychopharmacol (2013) 27(10):878–93. doi: 10.1177/0269881113499209
36. Kompauer I, Demmelmair H, Koletzko B, Bolte G, Linseisen J, Heinrich J. Association of Fatty Acids in Serum Phospholipids With Hay Fever, Specific and Total Immunoglobulin E. Br J Nutr (2005) 93(4):529–35. doi: 10.1079/bjn20041387
37. Berger L, Descamps V, Marck Y, Dehen L, Grossin M, Crickx B, et al. Alpha Interferon-Induced Eczema in Atopic Patients Infected by Hepatitis C Virus: 4 Case Reports. Ann Dermatol Venereol (2000) 127(1):51–5.
38. Wright CE, Strike PC, Brydon L, Steptoe A. Acute Inflammation and Negative Mood: Mediation by Cytokine Activation. Brain Behav Immun (2005) 19(4):345–50. doi: 10.1016/j.bbi.2004.10.003
39. Bullone M, Lavoie JP. The Equine Asthma Model of Airway Remodeling: From a Veterinary to a Human Perspective. Cell Tissue Res (2020) 380(2):223–36. doi: 10.1007/s00441-019-03117-4
40. Gao YH, Zhao HS, Zhang FR, Gao Y, Shen P, Chen RC, et al. The Relationship Between Depression and Asthma: A Meta-Analysis of Prospective Studies. PloS One (2015) 10(7):e0132424. doi: 10.1371/journal.pone.0132424
41. Nagel M, Jansen PR, Stringer S, Watanabe K, de Leeuw CA, Bryois J, et al. Meta-Analysis of Genome-Wide Association Studies for Neuroticism in 449,484 Individuals Identifies Novel Genetic Loci and Pathways. Nat Genet (2018) 50(7):920–7. doi: 10.1038/s41588-018-0151-7
42. Johansson A, Rask-Andersen M, Karlsson T, Ek WE. Genome-Wide Association Analysis of 350 000 Caucasians From the UK Biobank Identifies Novel Loci for Asthma, Hay Fever and Eczema. Hum Mol Genet (2019) 28(23):4022–41. doi: 10.1093/hmg/ddz175
43. Summer BE, Fink G. Estrogen Increases the Density of 5-Hydroxytryptamine(2A) Receptors in Cerebral Cortex and Nucleus Accumbens in the Female Rat. J Steroid Biochem Mol Biol (1995) 54(1-2):15–20. doi: 10.1016/0960-0760(95)00075-b
44. Ancelin ML, Scali J, Ritchie K. Hormonal Therapy and Depression: Are We Overlooking an Important Therapeutic Alternative? J Psychosom Res (2007) 62(4):473–85. doi: 10.1016/j.jpsychores.2006.12.019
45. Soares CN, Zitek B. Reproductive Hormone Sensitivity and Risk for Depression Across the Female Life Cycle: A Continuum of Vulnerability? J Psychiatry Neurosci (2008) 33(4):331–43.
46. Zweifel JE, O’Brien WH. A Meta-Analysis of the Effect of Hormone Replacement Therapy Upon Depressed Mood. Psychoneuroendocrinology (1997) 22(3):189–212. doi: 10.1016/s0306-4530(96)00034-0
47. Gleason CE, Dowling NM, Wharton W, Manson JE, Miller VM, Atwood CS, et al. Effects of Hormone Therapy on Cognition and Mood in Recently Postmenopausal Women: Findings From the Randomized, Controlled KEEPS-Cognitive and Affective Study. PloS Med (2015) 12(6):e1001833. doi: 10.1371/journal.pmed.1001833. discussion e1001833.
48. Howard DM, Adams MJ, Clarke TK, Hafferty JD, Gibson J, Shirali M, et al. Genome-Wide Meta-Analysis of Depression Identifies 102 Independent Variants and Highlights the Importance of the Prefrontal Brain Regions. Nat Neurosci (2019) 22(3):343–52. doi: 10.1038/s41593-018-0326-7
49. Limon-Morales O, Arteaga-Silva M, Rojas-Castaneda JC, Molina-Jimenez T, Guadarrama-Cruz GV, Cerbon M, et al. Neonatal Treatment With Clomipramine Modifies the Expression of Estrogen Receptors in Brain Areas of Male Adult Rats. Brain Res (2019) 1724:146443. doi: 10.1016/j.brainres.2019.146443
50. Yung JA, Fuseini H, Newcomb DC. Hormones, Sex, and Asthma. Ann Allergy Asthma Immunol (2018) 120(5):488–94. doi: 10.1016/j.anai.2018.01.016
51. Fan Z, Che H, Yang S, Chen C. Estrogen and Estrogen Receptor Signaling Promotes Allergic Immune Responses: Effects on Immune Cells, Cytokines, and Inflammatory Factors Involved in Allergy. Allergol Immunopathol (Madr) (2019) 47(5):506–12. doi: 10.1016/j.aller.2019.03.001
52. Wang Y, Chen YJ, Xiang C, Jiang GW, Xu YD, Yin LM, et al. Discovery of Potential Asthma Targets Based on the Clinical Efficacy of Traditional Chinese Medicine Formulas. J Ethnopharmacol (2020) 252:112635. doi: 10.1016/j.jep.2020.112635
53. Kanda N, Hoashi T, Saeki H. The Roles of Sex Hormones in the Course of Atopic Dermatitis. Int J Mol Sci (2019) 20(19):4660. doi: 10.3390/ijms20194660
54. Klis K, Wronka I. Association of Estrogen-Related Traits With Allergic Rhinitis. Adv Exp Med Biol (2017) 968:71–8. doi: 10.1007/5584_2016_190
Keywords: major depressive disorder, Mendelian randomization, meta-analyses, asthma, atopic diseases
Citation: Cao H, Li S, Baranova A and Zhang F (2021) Shared Genetic Liability Between Major Depressive Disorder and Atopic Diseases. Front. Immunol. 12:665160. doi: 10.3389/fimmu.2021.665160
Received: 07 February 2021; Accepted: 16 August 2021;
Published: 08 September 2021.
Edited by:
Lam Cheung Tsoi, University of Michigan, United StatesReviewed by:
Cunyou Zhao, Southern Medical University, ChinaMona-Rita Yacoub, San Raffaele Hospital (IRCCS), Italy
Copyright © 2021 Cao, Li, Baranova and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fuquan Zhang, zhangfq@njmu.edu.cn
†These authors have contributed equally to this work