Mingrong Zuo
Zhihao Wang
Wenhao Li
Siliang Chen
Yunbo Yuan
Yuan Yang
Yanhui Liu- Department of Neurosurgery, West China Hospital, Sichuan University, Chengdu, China
Background: Postpartum depression (PPD) is a type of depressive episode related to parents after childbirth, which causes a variety of symptoms not only for parents but also affects the development of children. The causal relationship between potential risk factors and PPD remains comprehensively elucidated.
Methods: Linkage disequilibrium score regression (LDSC) analysis was conducted to screen the heritability of each instrumental variant (IV) and to calculate the genetic correlations between effective causal factors and PPD. To search for the causal effect of multiple potential risk factors on the incidence of PPD, random effects of the inverse variance weighted (IVW) method were applied. Sensitivity analyses, including weighted median, MR-Egger regression, Cochrane’s Q test, and MR Pleiotropy Residual Sum and Outlier (MR-PRESSO), were performed to detect potential Mendelian randomization (MR) assumption violations. Multivariable MR (MVMR) was conducted to control potential multicollinearity.
Results: A total of 40 potential risk factors were investigated in this study. LDSC regression analysis reported a significant genetic correlation of potential traits with PPD. MR analysis showed that higher body mass index (BMI) (Benjamini and Hochberg (BH) corrected p = 0.05), major depression (MD) (BH corrected p = 5.04E-19), and schizophrenia (SCZ) (BH corrected p = 1.64E-05) were associated with the increased risk of PPD, whereas increased age at first birth (BH corrected p = 2.11E-04), older age at first sexual intercourse (BH corrected p = 3.02E-15), increased average total household income before tax (BH corrected p = 4.57E-02), and increased years of schooling (BH corrected p = 1.47E-11) led to a decreased probability of PPD. MVMR analysis suggested that MD (p = 3.25E-08) and older age at first birth (p = 8.18E-04) were still associated with an increased risk of PPD.
Conclusion: In our MR study, we found multiple risk factors, including MD and younger age at first birth, to be deleterious causal risk factors for PPD.
1 Introduction
Postpartum depression (PPD), also termed postnatal depression, is a type of depressive episode related to childbirth with a variety of symptoms of mood changes, negative attitudes to life, and mental and psychiatric changes, with few severe cases showing suicidality (1). The diagnostic criteria of PPD comprise five or more symptoms, such as depressed mood, loss of previous interest, weight or appetite alteration, somnipathy, fatigue, attention deficit, and feeling worthless or guilty (2). PPD occurs in approximately 13% of women, which may cause severe consequences for mothers and their children (3). Moreover, the reported overall incidence of paternal PPD between 3 months and 12 months after birth was about 10.4% (4), which suggested that both new parents may suffer from this common psychiatric disorder. Apart from causing mental disorders or even suicides of parents, the risks of PPD on children have been well-investigated, and it suggested that children with parents experiencing PPD are more likely to suffer from depressive disorders (5). Regarding the risk factors associated with PPD, a variety of research found that multiple factors, including a history of mental illness, hormonal changes, and social factors, were related to PPD (1, 6). The causal relationship between potential risk factors and PPD is still unelucidated, and a comprehensive exploration of the causal effects of factors on PPD would be helpful for the intervention of PPD.
Mendelian randomization (MR) is a useful way to infer a causal relationship, using genetic variants as instrumental variants (IVs) to estimate the causal assumption between exposure and outcome (7). Due to the random allocation of genetic variants, MR is less susceptible to other confounding factors and reverse causation that can impede causal inference in conventional observational studies (8). Multivariable MR (MVMR) takes pleiotropy among multiple exposures into account, which is an extension of univariable MR and is important for pleiotropic pathways (9). Thus, accumulating research has taken advantage of MR to gain insights into the causes of multiple diseases (10). Recently, MR analysis initially demonstrated a potential causal relationship between opioid use and the risk of PPD. In turn, PPD was also associated with a higher risk of opioid and non-opioid analgesic use (11). Another recent study investigated the causal association between PPD and cerebrovascular diseases and cognitive impairment, and the result recommended that cognitive impairment was a significant outcome induced by PPD (12). In addition, modifiable risk factors, such as body mass index (BMI) and glucose, are usually considered exposures that may be associated with diseases in MR analysis (13). However, the causal association of other potential risk factors, particularly modifiable risk factors with PPD, still remains unclear.
In the present study, we focused on the association of various potential risk factors and modifiable factors with PPD. Based on the summarized data of the genome-wide association study (GWAS), we investigated the genetic correlations between 40 potential risk factors and PPD by linkage disequilibrium score regression (LDSC) analysis and then performed a two-sample MR analysis.
2 Materials and methods
2.1 Study design
A procedure of MR analysis was established to investigate a causal relationship between potential risk factors and PPD. Figure 1 shows the procedures of our study. The principles of selecting potential risk factors were as follows: First, we sought out a review of the PPD (14), which mentioned the risk factors of PPD, including history of mental illness (psychiatric disorders), hormonal changes (sex hormones), and social factors (socioeconomic Factors). Another study reported the association between thyroid function and PPD (15). Second, we investigated the modifiable risk factors frequently used as exposures in MR analysis, including BMI, blood pressure, glucose and lipids, and diet hobbies such as smoking, drinking alcohol, tea, and coffee. Third, we added some risk factors that may be involved in the pathogenesis of depression, including disorder of sleep (16) and inflammatory processes (inflammatory biomarkers) (17). In total, there were 40 potential risk factors.
Figure 1. Design of the study. Procedures of MR are shown in this flow chart. Three principles of MR were listed: the IVs must be associated with exposure (I) and not be associated with confounders (II) and outcome (III). MR, mendelian randomization; GWAS, genome-wide association studies; SNP, single nucleotide polymorphism; IV, instrumental variant; LDSC, linkage disequilibrium score regression.
Three major assumptions were utilized. First, there existed a direct correlation between the IVs of potential exposures and the PPD (Relevance). Second, the confounders cannot confound the IVs of exposures (Exchangeability). Finally, the IVs have no direct connection to the outcome except through exposure (the exclusion restriction).
2.2 Sources of exposures and PPD in GWAS data
We obtained data on exposures from the GWAS data. The exposure factors of this study are as follows: BMI, total testosterone, bioavailable testosterone, sex hormone binding globulin, estradiol, coffee intake, tea intake, average total household income before tax, age at menarche1 (18–20), Schizophrenia (SCZ) (21), autism spectrum disorder (22), bipolar disorder (23), major depression (MD) (24), attention-deficit/hyperactivity disorder (25), anxiety (26), diastolic blood pressure, pulse pressure, systolic blood pressure (27), glycated hemoglobin, fasting glucose, fasting insulin, 2-h glucose (28), triglycerides, low-density lipoprotein cholesterol-c, high-density lipoprotein cholesterol-c, apolipoprotein A-I, apolipoprotein B (29), free thyroxine4, thyrotropin (TSH), incTSH/hypothyroidism, decTSH/hyperthyroidism (30), 25 hydroxyvitamin D (31), C-reactive protein (32), smoking initiation (SI), cigarettes per day, alcohol consumption (33), years of schooling (33), age at first sexual intercourse, age at first birth (34), and insomnia (35). Next, the GWAS data of PPD were obtained from the European cohort: the FinnGen study round 8 (36). The FinnGen study is dedicated to combining genome information with digital healthcare data (37). The criterion in FinnGen was based on the 10th edition of the International Classification of Diseases criteria, in which participants with delivery history (O15) diagnosed with F32, F33, or F530 were identified as PPD. A total of 13,657 PPD cases and 236,178 controls were enrolled in the dataset. GWAS summary data of PPD and all risk factors are presented in Table 1. It is worth noting that the ethnicity of the population of attention deficit/hyperactivity disorder was European, North American, and Chinese, and the ethnicity of other GWAS data was European. Given that the diagnostic criteria of the PPD included a history of delivery, we speculated that the gender of this study was restricted to biological females. We elected GWAS of sex hormones, free thyroxine, and TSH of female-only participants as the rest of the studies did not provide gender-specified data.
Table 1. Description of GWAS statistics included in the present study.
2.3 Linkage disequilibrium score regression
To screen the heritability of each trait and the genetic correlations between effective causal factors and PPD, LDSC regression was conducted by the regression slope using GWAS summary data (38). The European Ancestry 1,000 Genomes LD reference panel was used as a reference.
2.4 Genetic instruments selection
To select the qualified IVs, only single nucleotide polymorphisms (SNPs) associated with each exposure at a genome-wide threshold (p < 5E-08) were elected. A less strict threshold of 5E-06 was applied if less than four IVs were involved in the inverse variance weighted (IVW) analysis. SNPs with linkage disequilibrium were filtered based on the European ancestry 1,000 Genomes LD reference panel, with r2 > 0.001 on the clump window of 10,000 kb to keep the IVs independent from each other. The F-statistic of each SNP was calculated, and any SNP with a low F-statistic (< 10) was removed to avoid weak instrument bias. SNPs significantly associated with outcome (p < 5E-08) were excluded as they violated the third principle of MR assumption. SNPs that passed these precedes and existed in the GWAS of outcome were harmonized to avoid a mismatch of alleles based on allele frequencies for palindromes (19).
2.5 Mendelian randomization analysis
We used IVW to analyze the causal relationship of 40 potential risk factors in the PPD (7). In addition, weighted median, MR-Egger regression, heterogeneity test, Cochrane’s Q test, and MR-PRESSO were utilized to assess the robustness of IVW results. We appraised the heterogeneity using Cochran’s Q statistics and I2 statistics (39), and if heterogeneity exists (p < 0.05), the results should coincide with results estimated by IVW (39). In the pleiotropy test, intercepts calculated by MR-PRESSO were used to assess the horizontal pleiotropy of valid IVs; p > 0.05 meant a lower probability of horizontal pleiotropy existence (40). Outliners were detected and removed by MR-PRESSO. The validation of causal inference between exposures and PPD depended on the same tendency of IVW and weighted median and MR-Egger analyses, besides no horizontal pleiotropic effect existing. To control potential multicollinearity, multivariable MVMR was conducted with significant results of the IVW method as exposures. For MVMR, the IVs must be associated with at least one of the exposures, independent of all confounders of any exposure and independent of the outcome (41). Least absolute shrinkage and selection operator regression were utilized to avoid potential bias caused by multicollinearity.
2.6 Statistical analysis
All the statistical analyses were performed using R-4.2.32 with R packages. The R packages included the two-sample MR package (19), MR-PRESSO package (40), and LDSC package3 (38, 42, 43). To control type I error, BH (Benjamini and Hochberg) correction was employed in both LDSC and IVW analyses. The BH-adjusted p < 0.05 was considered strong evidence, while the standard p < 0.05 was identified as suggestive evidence. Scatter plots and forest plots were used for the visualization of MR results.
2.7 Ethical approval
The data from public sources of our study had been granted ethical approval by their own institutional review boards and thus it did not require ethical approval from our institutional review board.
3 Results
3.1 Linkage disequilibrium score regression
We performed LDSC to analyze genetic correlations between potential causal factors and PPD. The heritability of SNPs for each exposure factor with a range of 2 to 36% is listed in Supplementary Table S1. Among 40 factors, half of them were associated with PPD even following BH correction, as listed in Table 2.
Table 2. Genetic correlation between PPD and potential risk factors.
3.2 IVs selection
SNPs were selected for each pair of exposure and outcome, and detailed information on all SNPs was listed in the Supplementary Table S2. We defined the qualified IVs as only SNPs associated with each exposure at a genome-wide threshold (p < 5E-08). Due to the lack of IVs at the genome-wide significance with anxiety, autism spectrum disorder, and estradiol as exposure, the threshold of 5E-06 was applied in the three analyses. The total F-statistic of all pairs was above 10, indicating strong IVs.
3.3 IVW analysis of potential risk factors for PPD
Our results of MR analyses indicated that among 40 potential risk factors, eight of them had a significant causal relationship with the PPD. The functionalities of SNPs, such as higher BMI (odds ratio (OR) = 1.13 (1.03, 1.24), BH corrected p = 0.05), MD (OR = 2.05 (1.76, 2.39), BH corrected p = 5.04E-19), SCZ (OR = 1.13 (1.08, 1.19), BH corrected p = 1.64E-05), and earlier SI (OR = 1.33 (1.14, 1.54), BH corrected p = 1.33E-03), were associated with the increased risk of PPD. On the contrary, we also found that increased age at first birth (OR = 0.86 (0.81, 0.93), BH corrected p = 2.11E-04), older age at first sexual intercourse (OR = 0.47 (0.4, 0.57), BH corrected p = 3.02E-15), increased average total household income before tax (OR = 0.62 (0.44, 0.88), BH corrected p = 4.57E-02), and increased years of schooling (OR = 0.56 (0.48, 0.66), BH corrected p = 1.47E-11) led to a decreased probability of PPD. While no significant causal relationship was found between the other 32 exposures and PPD. Complete results are available in Figure 2 and Supplementary Table S3.
Figure 2. Forest plots of causal effect estimates based on Mendelian analysis. The significant results of causal inference of eight potential risk factors on PPD in random effect IVW analysis were plotted. SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval.
3.4 Sensitivity analysis of MR
The scatter plots of effective exposures with causal effects on PPD are presented in Figure 3. Apart from SI showing a discrepant correlation with PPD between MR Egger and IVW analysis and weighted median analysis (Figure 3), the rest of the effective causal effectors all presented a clear causal relationship with a similar tendency, which was in line with each corresponding forest plot (Figures 2, 3). Heterogeneity was detected in six causal factors identified by IVW, which are as follows: age at first sexual intercourse, average total household income before tax, BMI, SCZ, SI, and years of schooling. The pleiotropy test showed no horizontal pleiotropy in all analyses (p > 0.05). MR-PRESSO analysis suggested outliners were presented in the following exposures: average total household income before tax, BMI, SCZ, and years of schooling, deleting those that did not affect the results. The complete results of the sensitivity analysis are available in Supplementary Table S4.
Figure 3. Scatter plots of causal effect inference for effective causal factors on PPD, respectively. SNP, single nucleotide polymorphism; BMI, body mass index.
3.5 Multivariable Mendelian randomization
We conducted an MVMR analysis with eight potential causal factors identified by IVW analysis. IVs were selected from shared SNPs of all exposure GWAS (Supplementary Table S5). None of these exposures were excluded after the least absolute shrinkage and selection operator regression. Taking all eight factors into account, genetically predicted MD was still associated with an increased risk of PPD (OR = 1.63 (1.37, 1.94), p = 3.25E-08). In contrast, older age at first birth was associated with a lower risk of PPD (OR = 0.81 (0.72, 0.92), p = 8.18E-04). The rest of the factors were not significant. Detailed information is listed in Table 3.
Table 3. MVMR analysis of eight significant causal factors.
4 Discussion
In the present MR study, we found that higher BMI, MD, SCZ, and earlier SI increased the risk of PPD. On the contrary, the risk factors of increased age at first birth, older age at first sexual intercourse, increased income, and higher education level led to a decreased probability of PPD. In sensitivity analysis, contradictory results of SI were reported, which should be interpreted with caution. The other 32 risk factors did not show a statistically causal correlation with PPD. Adjustment for all eight significant traits in the MVMR model further found that MD and older age at first birth were still associated with an increased risk of PPD. Thus, our study shed fresh light on multiple causal factors associated with PPD.
Younger maternal age was reported to be associated with an increased risk of PPD in women (44). There was an MR analysis illustrating that the older age of first birth was a protective factor for mental disorders (45). Consistently, we found a negative causal relationship between age at first birth and PPD; thus, we suggested that relatively late-to-be parents may be less likely to develop PPD. Early sexual intercourse was associated with mental health outcomes (46). A novel MR study established a risky causal relationship between early sexual intercourse and MD (47). In our study, we initially brought forward that the experience of early sexual intercourse exerted a deleterious effect on PPD. As one of the social determinants, higher education level is identified as a protective factor for PPD (48). Our study also suggested that longer years of schooling were negatively associated with the risks of PPD, which was readily understandable as the higher the level of education, the higher the capability to withstand health and wealth risks. Whether income level or urbanization is associated with depression remains controversial (49). Here, we found that a lower average total household income before tax was associated with an increased risk of PPD. This conclusion needs further investigation using genetic biobanks of middle- and low-income countries. Based on these results, clinicians and parents should take advantage of these viable socioeconomic factors to prevent PPD, such as avoiding earlier first birth and immature sexual intercourse and acquiring a higher education level and a higher income.
Regarding the negative causalities of PPD, we found that modifiable risk factors, both higher BMI and earlier SI, significantly increased the likelihood of PPD. Recently, MR analysis has already found a promoting role of BMI in the risk of MD. In addition, it suggested that depression increases the genetic susceptibility to high BMI (50), which means a positive feedback between BMI and depression and needs the clinician’s close attention and intervention. Combined with our findings, further clinical trials should investigate the effect of preventing PPD by getting rid of excess weight. Non-smoking has been considered a protective causal factor in the incidence of MD (51). This was consistent with our MR results, which indicated a direct effect of earlier SI on an increased risk for PPD. As the risk factors above are relatively easy to modify, we postulate that reducing BMI and postponing SI may decrease the risk of PPD. Although insomnia was associated with an increased risk of perinatal depression (52), we did not conclude a causal relationship between insomnia and PPD. An observational study reported that 11% of women and 8% of men had PPD after the first birth, and it was partly attributed to postpartum sleep deprivation (53). In contrast, our finding suggested no causal effect of insomnia on PPD. Although insomnia is prevalently comorbid with depression, and treating insomnia is beneficial to the mood of patients (54, 55), further research is warranted to identify the role of insomnia in the pathogenesis of PPD.
Finally, we reported that a history of MD increased the risk of PPD, which was in line with previous findings (56, 57). In addition, there was evidence of a positive causal effect between SCZ and PPD in our study. An observational study suggested that PPD and postpartum psychosis were related to a higher risk of SCZ (58). Our MR analysis further strengthened the robustness of the positive causal connection between mental disorders, such as SCZ and MD, and the risk of PPD. More than half of women with PPD have a comorbidity of bipolar disorder, which usually being misdiagnosed or neglected (59). In our study, we did not conclude a causal relationship between these two diseases. We speculated that the negative result may be related to limited cases, distinct diagnostic criteria, and sequencing depth for bipolar disorder. Anyway, future studies should try to illustrate the association between bipolar disorder and PPD. Subsequently, we conducted an MVMR analysis to control potential multicollinearity. We included all eight significant traits from the IVW analysis, of which the MR Egger analysis showed that SI had an opposite direction compared to IVW and weighted median analysis, which should be interpreted with caution. In addition, BMI did not reach a significance following BH correcting. Importantly, compared to univariable MR, the effect of MD and older age at first birth on the increased risks of PPD were still observed. Thus, our study brought forward the most important causal factors of PPD, which may help future studies to better prevent and treat PPD.
Our study has limitations. First, the data sources are mainly from participants with European ancestry; thus, our result inevitably may not adapt to different ethnic groups, and the non-European participants in the GWAS of attention-deficit/hyperactivity disorder might lead to biases in results due to variations in allele frequencies and genetic associations among different populations and ancestry-specific effects. Second, the causal relationship of our findings does not conduct external verification as we only obtained one GWAS data of PPD. Therefore, we utilize more approaches to promote the robustness of our results, such as using the MR egger and weighted median analysis. In addition, one of the limitations of the methodology of MR is the potential confounding of the genetic variants and the outcome (Violation of IV condition two: exchangeability). Other than that, another bias of MR analysis that is uneasily corrected is the selection bias (10). In all, we should be prudent with our findings of MR because of the limitations above.
In summary, our MR analysis provided suggestive evidence of the protective effect of older age at first birth and older age at first sexual intercourse, higher education, and higher income on the risk of PPD. We also suggested that higher BMI and a history of MD and SCZ were hazardous causal factors for PPD. In addition, MD and older age at first birth were independently related to the increased risks of PPD. These novel shreds of evidence may help to guide the prevention and intervention strategies for PPD.
Data availability statement
Publicly available datasets were analyzed in this study. See Table 1 for detailed information. Further inquiries can be directed to the corresponding author.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
MZ: Conceptualization, Data curation, Investigation, Methodology, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. ZW: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. WL: Data curation, Methodology, Validation, Visualization, Writing – review & editing. SC: Formal analysis, Investigation, Validation, Writing – review & editing. YunY: Data curation, Formal analysis, Investigation, Project administration, Validation, Writing – review & editing. YuaY: Funding acquisition, Supervision, Validation, Writing – review & editing. QM: Data curation, Supervision, Validation, Writing – review & editing. YL: Conceptualization, Funding acquisition, Investigation, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Sichuan Science and Technology Program (2023YFQ0002 to YL) and the Sichuan Science and Technology Program (2023YFG0127 to YuaY).
Acknowledgments
The authors sincerely thank all the participants and investigators for constructing and publishing the summarized GWAS data, the developers of the R software, and all the related algorithms and websites.
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2023.1275834/full#supplementary-material
Footnotes
References
1. Field, T. Postpartum depression effects on early interactions, parenting, and safety practices: a review. Infant Behav Dev. (2010) 33:1–6. doi: 10.1016/j.infbeh.2009.10.005
2. American Psychiatric Association D, Association AP. Diagnostic and statistical manual of mental disorders: DSM-5. Washington, DC: American Psychiatric Association (2013).
3. Musters, C, McDonald, E, and Jones, I. Management of postnatal depression. BMJ. (2008) 337:a736. doi: 10.1136/bmj.a736
4. Lancet, T. Postnatal depression: fathers have it too. Lancet. (2010) 375:1846. doi: 10.1016/S0140-6736(10)60853-7
5. Kassier, S, and Madlala, S. Antenatal and postpartum depression: effects on infant and young child health and feeding practices. South Afr J Clin Nutr. (2018) 26:17–22. doi: 10.1080/16070658.2013.11734435
6. Howard, LM, Molyneaux, E, Dennis, C-L, Rochat, T, Stein, A, and Milgrom, J. Non-psychotic mental disorders in the perinatal period. Lancet. (2014) 384:1775–88. doi: 10.1016/S0140-6736(14)61276-9
7. Burgess, S, Butterworth, A, and Thompson, SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. (2013) 37:658–65. doi: 10.1002/gepi.21758
8. Davey Smith, G, and Ebrahim, S. “Mendelian randomization”: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. (2003) 32:1–22. doi: 10.1093/ije/dyg070
9. Zuber, V, Colijn, JM, Klaver, C, and Burgess, S. Selecting likely causal risk factors from high-throughput experiments using multivariable Mendelian randomization. Nat Commun. (2020) 11:29. doi: 10.1038/s41467-019-13870-3
10. Sanderson, E, Glymour, MM, Holmes, MV, Kang, H, Morrison, J, Munafò, MR, et al. Mendelian randomization. Nat Rev Methods Primers. (2022) 2:6. doi: 10.1038/s43586-021-00092-5
11. Jiang, Y, Wei, D, and Xie, Y. Causal effects of opioids on postpartum depression: a bidirectional, two-sample Mendelian randomization study. Front Psych. (2023) 14:1043854. doi: 10.3389/fpsyt.2023.1043854
12. Li, J, Li, J, Shen, L, Wang, H, Zheng, T, Hui, Y, et al. Investigating the causal association of postpartum depression with cerebrovascular diseases and cognitive impairment: a Mendelian randomization study. Front Psych. (2023) 14:1196055. doi: 10.3389/fpsyt.2023.1196055
13. Saunders, CN, Cornish, AJ, Kinnersley, B, Law, PJ, Claus, EB, Il’yasova, D, et al. Lack of association between modifiable exposures and glioma risk: a Mendelian randomization analysis. Neuro-Oncology. (2020) 22:207–15. doi: 10.1093/neuonc/noz209
14. Stewart, DE, and Vigod, S. Postpartum depression. N Engl J Med. (2016) 375:2177–86. doi: 10.1056/NEJMcp1607649
15. Minaldi, E, D’Andrea, S, Castellini, C, Martorella, A, Francavilla, F, Francavilla, S, et al. Thyroid autoimmunity and risk of post-partum depression: a systematic review and meta-analysis of longitudinal studies. J Endocrinol Investig. (2020) 43:271–7. doi: 10.1007/s40618-019-01120-8
16. Lawson, A, Murphy, KE, Sloan, E, Uleryk, E, and Dalfen, A. The relationship between sleep and postpartum mental disorders: a systematic review. J Affect Disord. (2015) 176:65–77. doi: 10.1016/j.jad.2015.01.017
17. Liu, H, Zhang, Y, Gao, Y, and Zhang, Z. Elevated levels of Hs-CRP and IL-6 after delivery are associated with depression during the 6 months post partum. Psychiatry Res. (2016) 243:43–8. doi: 10.1016/j.psychres.2016.02.022
18. Elsworth, B, Lyon, M, Alexander, T, Liu, Y, Matthews, P, Hallett, J, et al. The MRC IEU OpenGWAS data infrastructure. BioRxiv. (2020) doi: 10.1101/2020.08.10.244293
19. Hemani, G, Zheng, J, Elsworth, B, Wade, KH, Haberland, V, Baird, D, et al. The MR-base platform supports systematic causal inference across the human phenome. elife. (2018) 7:e34408. doi: 10.7554/eLife.34408
20. Lyon, MS, Andrews, SJ, Elsworth, B, Gaunt, TR, Hemani, G, and Marcora, E. The variant call format provides efficient and robust storage of GWAS summary statistics. Genome Biol. (2021) 22:32. doi: 10.1186/s13059-020-02248-0
21. Trubetskoy, V, Pardiñas, AF, Qi, T, Panagiotaropoulou, G, Awasthi, S, Bigdeli, TB, et al. Mapping genomic loci implicates genes and synaptic biology in schizophrenia. Nature. (2022) 604:502–8. doi: 10.1038/s41586-022-04434-5
22. Autism Spectrum Disorder Working Group of the Psychiatric Genomics Consortium, BUPGEN, Major Depressive Disorder Working Group of the Psychiatric Genomics Consortium, 23andMe Research Team Grove, J, Ripke, S, Als, TD, Mattheisen, M, Walters, RK, et al. Identification of common genetic risk variants for autism spectrum disorder. Nat Genet. (2019) 51:431–44. doi: 10.1038/s41588-019-0344-8
23. Hou, L, Bergen, SE, Akula, N, Song, J, Hultman, CM, Landén, M, et al. Genome-wide association study of 40,000 individuals identifies two novel loci associated with bipolar disorder. Hum Mol Genet. (2016) 25:3383–94. doi: 10.1093/hmg/ddw181
24. Howard, DM, Adams, MJ, Clarke, T-K, 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:343–52. doi: 10.1038/s41593-018-0326-7
25. ADHD Working Group of the Psychiatric Genomics Consortium (PGC), Early Lifecourse & Genetic Epidemiology (EAGLE) Consortium, 23andMe Research Team Demontis, D, Walters, RK, Martin, J, Mattheisen, M, Als, TD, et al. Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat Genet. (2019) 51:63–75. doi: 10.1038/s41588-018-0269-7
26. Otowa, T, Hek, K, Lee, M, Byrne, EM, Mirza, SS, Nivard, MG, et al. Meta-analysis of genome-wide association studies of anxiety disorders. Mol Psychiatry. (2016) 21:1391–9. doi: 10.1038/mp.2015.197
27. Surendran, P, Feofanova, EV, Lahrouchi, N, Ntalla, I, Karthikeyan, S, Cook, J, et al. Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals. Nat Genet. (2020) 52:1314–32. doi: 10.1038/s41588-020-00713-x
28. Chen, J, Spracklen, CN, Marenne, G, Varshney, A, Corbin, LJ, Luan, J’, et al. The trans-ancestral genomic architecture of glycemic traits. Nat Genet. (2021) 53:840–60. doi: 10.1038/s41588-021-00852-9
29. Richardson, TG, Sanderson, E, Palmer, TM, Ala-Korpela, M, Ference, BA, Davey Smith, G, et al. Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: a multivariable Mendelian randomisation analysis. PLoS Med. (2020) 17:e1003062. doi: 10.1371/journal.pmed.1003062
30. Teumer, A, Chaker, L, Groeneweg, S, Li, Y, di Munno, C, Barbieri, C, et al. Genome-wide analyses identify a role for SLC17A4 and AADAT in thyroid hormone regulation. Nat Commun. (2018) 9:4455. doi: 10.1038/s41467-018-06356-1
31. Revez, JA, Lin, T, Qiao, Z, Xue, A, Holtz, Y, Zhu, Z, et al. Genome-wide association study identifies 143 loci associated with 25 hydroxyvitamin D concentration. Nat Commun. (2020) 11:1–12. doi: 10.1038/s41467-020-15421-7
32. Said, S, Pazoki, R, Karhunen, V, Võsa, U, Ligthart, S, Bodinier, B, et al. Genetic analysis of over half a million people characterises C-reactive protein loci. Nat Commun. (2022) 13:2198. doi: 10.1038/s41467-022-29650-5
33. 23andMe Research Team, HUNT All-In Psychiatry Liu, M, Jiang, Y, Wedow, R, Li, Y, Brazel, DM, et al. Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat Genet. (2019) 51:237–44. doi: 10.1038/s41588-018-0307-5
34. Mills, MC, Tropf, FC, Brazel, DM, van Zuydam, N, Vaez, A, eQTLGen Consortium, et al. Identification of 371 genetic variants for age at first sex and birth linked to externalising behaviour. Nat Hum Behav. (2021) 5:1717–30. doi: 10.1038/s41562-021-01135-3
35. HUNT All In Sleep Lane, JM, Jones, SE, Dashti, HS, Wood, AR, Aragam, KG, et al. Biological and clinical insights from genetics of insomnia symptoms. Nat Genet. (2019) 51:387–93. doi: 10.1038/s41588-019-0361-7
36. Kurki, MI, Karjalainen, J, Palta, P, Sipilä, TP, Kristiansson, K, Donner, KM, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. (2023) 613:508–18. doi: 10.1038/s41586-022-05473-8
37. Kurki, MI, Karjalainen, J, Palta, P, Sipilä, TP, Kristiansson, K, Donner, K, et al. FinnGen: unique genetic insights from combining isolated population and national health register data. MedRxiv. (2022) doi: 10.1101/2022.03.03.22271360
38. Bulik-Sullivan, BK, Loh, P-R, Finucane, HK, Ripke, S, and Yang, J. Consortium SWGotPG, et al. LD score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet. (2015) 47:291–5. doi: 10.1038/ng.3211
39. Huedo-Medina, TB, Sánchez-Meca, J, Marín-Martínez, F, and Botella, J. Assessing heterogeneity in meta-analysis: Q statistic or I2 index? Psychol Methods. (2006) 11:193–206. doi: 10.1037/1082-989X.11.2.193
40. Verbanck, M, Chen, C-Y, Neale, B, and Do, R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
41. Sanderson, E. Multivariable Mendelian randomization and mediation. Cold Spring Harb Perspect Med. (2020) 11:a038984. doi: 10.1101/cshperspect.a038984
42. ReproGen Consortium, Psychiatric Genomics Consortium, Genetic Consortium for Anorexia Nervosa of the Wellcome Trust Case Control Consortium 3 Bulik-Sullivan, B, Finucane, HK, Anttila, V, Gusev, A, Day, FR, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. (2015) 47:1236–41. doi: 10.1038/ng.3406
43. Grotzinger, AD, Rhemtulla, M, de Vlaming, R, Ritchie, SJ, Mallard, TT, Hill, WD, et al. Genomic structural equation modelling provides insights into the multivariate genetic architecture of complex traits. Nat Hum Behav. (2019) 3:513–25. doi: 10.1038/s41562-019-0566-x
44. Agnafors, S, Bladh, M, Svedin, CG, and Sydsjö, G. Mental health in young mothers, single mothers and their children. BMC Psychiatry. (2019) 19:1–7. doi: 10.1186/s12888-019-2082-y
45. Cai, J, Wei, Z, Chen, M, He, L, Wang, H, Li, M, et al. Socioeconomic status, individual behaviors and risk for mental disorders: a Mendelian randomization study. Eur Psychiatry. (2022) 65:e28. doi: 10.1192/j.eurpsy.2022.18
46. Sabia, JJ. Does early adolescent sex cause depressive symptoms? J Policy Anal Manage. (2006) 25:803–25. doi: 10.1002/pam.20209
47. Lu, Z, Sun, Y, Liao, Y, Kang, Z, Feng, X, Zhao, G, et al. Identifying causal associations between early sexual intercourse or number of sexual partners and major depressive disorders: a bidirectional two-sample Mendelian randomization analysis. J Affect Disord. (2023) 333:121–9. doi: 10.1016/j.jad.2023.04.079
48. De Oliveira, G, Cianelli, R, Gattamorta, K, Kowalski, N, and Peragallo, N. Social determinants of depression among Hispanic women. J Am Psychiatr Nurses Assoc. (2017) 23:28–36. doi: 10.1177/1078390316669230
49. Sampson, L, Ettman, CK, and Galea, S. Urbanization, urbanicity, and depression: a review of the recent global literature. Curr Opin Psychiatry. (2020) 33:233–44. doi: 10.1097/YCO.0000000000000588
50. Mulugeta, A, Zhou, A, Vimaleswaran, KS, Dickson, C, and Hyppönen, E. Depression increases the genetic susceptibility to high body mass index: evidence from UK biobank. Depress Anxiety. (2019) 36:1154–62. doi: 10.1002/da.22963
51. Wan, B, Wu, Y, Ma, N, Zhou, Z, and Lu, W. Four modifiable factors that mediate the effect of educational time on major depressive disorder risk: a network Mendelian randomization study. PLoS One. (2023) 18:e0288034. doi: 10.1371/journal.pone.0288034
52. Yang, Q, Borges, MC, Sanderson, E, Magnus, MC, Kilpi, F, Collings, PJ, et al. Associations between insomnia and pregnancy and perinatal outcomes: evidence from mendelian randomization and multivariable regression analyses. PLoS Med. (2022) 19:e1004090. doi: 10.1371/journal.pmed.1004090
53. Parfitt, Y, and Ayers, S. Transition to parenthood and mental health in first-time parents. Infant Ment Health J. (2014) 35:263–73. doi: 10.1002/imhj.21443
54. Gebara, MA, Siripong, N, DiNapoli, EA, Maree, RD, Germain, A, Reynolds, CF, et al. Effect of insomnia treatments on depression: a systematic review and meta-analysis. Depress Anxiety. (2018) 35:717–31. doi: 10.1002/da.22776
55. Sharma, V, Sharkey, KM, Palagini, L, Mazmanian, D, and Thomson, M. Preventing recurrence of postpartum depression by regulating sleep. Expert Rev Neurother. (2023) 23:1–9. doi: 10.1080/14737175.2023.2237194
56. Lancaster, CA, Gold, KJ, Flynn, HA, Yoo, H, Marcus, SM, and Davis, MM. Risk factors for depressive symptoms during pregnancy: a systematic review. Am J Obstet Gynecol. (2010) 202:5–14. doi: 10.1016/j.ajog.2009.09.007
57. Wisner, KL, Sit, DK, McShea, MC, Rizzo, DM, Zoretich, RA, Hughes, CL, et al. Onset timing, thoughts of self-harm, and diagnoses in postpartum women with screen-positive depression findings. JAMA Psychiatry. (2013) 70:490–8. doi: 10.1001/jamapsychiatry.2013.87
58. Chen, M-H, Pan, T-L, Bai, Y-M, Huang, K-L, Tsai, S-J, Su, T-P, et al. Postpartum depression and psychosis and subsequent severe mental illnesses in mothers and neurodevelopmental disorders in children: a nationwide study. J Clin Psychiatry. (2021) 82:20m35520. doi: 10.4088/JCP.20m13735
59. Sharma, V, and Khan, M. Identification of bipolar disorder in women with postpartum depression. Bipolar Disord. (2010) 12:335–40. doi: 10.1111/j.1399-5618.2010.00809.x
Glossary
Keywords: postpartum depression, causal factors, mental disorder, Mendelian randomization, linkage disequilibrium score regression
Citation: Zuo M, Wang Z, Li W, Chen S, Yuan Y, Yang Y, Mao Q and Liu Y (2023) Causal effects of potential risk factors on postpartum depression: a Mendelian randomization study. Front. Psychiatry. 14:1275834. doi: 10.3389/fpsyt.2023.1275834
Edited by:
Sami Hamdan, Academic College Tel Aviv-Yaffo, IsraelReviewed by:
Weiqiu Cheng, University of Michigan, United StatesJun Yang, Zhejiang University, China
Copyright © 2023 Zuo, Wang, Li, Chen, Yuan, Yang, Mao and Liu. 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: Yanhui Liu, liuyh@scu.edu.cn
†These authors have contributed equally to this work and share first authorship