Skip to main content

DNA methylation in cord blood in association with prenatal depressive symptoms

Abstract

Background

Prenatal symptoms of depression (PND) and anxiety affect up to every third pregnancy. Children of mothers with mental health problems are at higher risk of developmental problems, possibly through epigenetic mechanisms together with other factors such as genetic and environmental. We investigated DNA methylation in cord blood in relation to PND, taking into consideration a history of depression, co-morbidity with anxiety and selective serotonin reuptake inhibitors (SSRI) use, and stratified by sex of the child. Mothers (N = 373) prospectively filled out web-based questionnaires regarding mood symptoms and SSRI use throughout pregnancy. Cord blood was collected at birth and DNA methylation was measured using Illumina MethylationEPIC array at 850 000 CpG sites throughout the genome. Differentially methylated regions were identified using Kruskal–Wallis test, and Benjamini-Hochberg adjusted p-values < 0.05 were considered significant.

Results

No differential DNA methylation was associated with PND alone; however, differential DNA methylation was observed in children exposed to comorbid PND with anxiety symptoms compared with healthy controls in ABCF1 (log twofold change − 0.2), but not after stratification by sex of the child. DNA methylation in children exposed to PND without SSRI treatment and healthy controls both differed in comparison with SSRI exposed children at several sites and regions, among which hypomethylation was observed in CpGs in the promoter region of CRBN (log2 fold change − 0.57), involved in brain development, and hypermethylation in MDFIC (log2 fold change 0.45), associated with the glucocorticoid stress response.

Conclusion

Although it is not possible to assess if these methylation differences are due to SSRI treatment itself or to more severe depression, our findings add on to existing knowledge that there might be different biological consequences for the child depending on whether maternal PND was treated with SSRIs or not.

Introduction

Common complications during pregnancy regard poor mental health. Maternal prenatal symptoms of depression (PND) prevalence range from 4 to 29%, depending on context and assessment methods [1], with the potential to affect the development of the child [2]. In a previous study from our group, PND (compared with postpartum depression) was particularly influential on child behavioral problems [3]. Furthermore, children exposed to maternal depression during pregnancy are reported to be at increased risk of altered brain development [4] as well as cognitive difficulties [5, 6]. A sex specific impact from exposure to PND has also been suggested, such as larger effect on cognitive development in boys and on emotional development in girls [5].

Anxiety is a common mental health problem in pregnant women. High rates of anxiety symptoms have been reported among pregnant women, varying between 4.4 and 39% [7]. Anxiety is often comorbid with depression and depression in the peripartum period is sometimes viewed as comprising much of the anxiety phenotype [5]. Depressive symptoms in women with a history of depression [8, 9] or in women with co-morbid antenatal symptoms of anxiety [10] are thus of special interest when considering eventual differences in pathogenesis and consequences. Finally, treatment with selective serotonin reuptake inhibitor (SSRI) has been shown to be associated with biological correlates such as altered cortico-releasing hormone levels and differential gene expression in the placenta [11, 12].

The link between maternal mental health problems and child development is complex and includes genetic factors, some with possibly pleiotropic effects, as well as environmental and family factors [13], which could impact on fetal development through many biological pathways. Greater levels of stress hormones and inflammation have been reported in cord blood of preterm neonates exposed to poor maternal mental health [14]. Apart from direct impact on fetal biological systems, alterations resulting from maternal depression during pregnancy are thought to induce epigenetic changes in developmental pathways. A hypothesis of epigenetic programming due to maternal behavior has evolved and includes research on experimental animal models [15]. There is also a growing body of human studies investigating the impact of prenatal exposure to maternal depressed mood in infants [16,17,18].

The focus of these epigenetic studies has been on DNA methylation, which is a chemical signature that can be alongside DNA replication, through transfer of methylation groups by DNA methyltransferase, and influence gene expression, thus representing a candidate mechanism for long-term effects of environmental exposure on phenotype [19]. The role of epigenetics in neural development has been mainly explored in epigenome-wide association studies (EWAS). In fetal brain samples, associations have been demonstrated between DNA methylation and brain development, these specific sites have also been seen to be schizophrenia risk loci [20]. The use of the more accessible surrogate tissue such as blood for these investigations has also been evaluated [21, 22]. Common environmental exposures could cause both epigenetic modifications in umbilical cord blood and brain, affecting the child neural development, for example, paracetamol usage during pregnancy is associated with attention deficit/hyperactivity disorder (ADHD) (reviewed by Bauer et al. [23]) and long-term exposure to paracetamol has shown altered DNA methylation in cord blood in children who later develop ADHD [24].

Studies on how fetal exposure to depressive symptoms alone, or with comorbidities, with or without SSRI treatment, affects DNA methylation have shown varied results, largely due to different methodologies and cohorts [16, 25,26,27]. Of special interest is the difference in cell composition between cord blood and whole blood from adults [28, 29] and previous studies on the association between PND and DNA methylation in cord blood have, however, not taken this into account [16, 17]. The sex of the child may also moderate the epigenetic influence by PND, since there are differences between sexes in DNA methylation patterns [30].

Altered DNA methylation in cord blood has been observed in children exposed to anxiety and also exposure to SSRIs has been associated with hypomethylation in ZNF575, a gene whose functions are largely unknown [16]. However, this study did not consider comorbidity of PND with anxiety nor did it compare exposure to SSRI with exposure to untreated PND [16]. A review by Viuff and colleagues assessed studies investigating the effect of SSRIs on cord blood DNA methylation and showed the need of new investigations using wide-range DNA methylation analysis, validation of exposure to SSRI, and analyses combining PND with other mental illness [31].

The aim of the present study was therefore to investigate the possible association between PND and DNA methylation in cord blood. Secondary aims, unique to this investigation, were to reassess the association among different subtypes of PND and to investigate the effect of sex. Specifically, subtypes were those with a history of depression or co-morbidity with symptoms of anxiety in comparison with PND alone, and also those treated with SSRIs in comparison with untreated PND and healthy controls.

Results

Demographics

Characteristics of the participating mothers are presented in Table 1. A university level education was less frequent among mothers with PND. Comorbid anxiety was reported in over 57 percent of the women with PND. The children of mothers with and without PND had not significantly different birth weight. There were no significant difference in the proportion of mothers born in Scandinavia between any of the analyzed study groups. Children born by women with PND had slightly shorter gestational length.

Table 1 Characteristics of participating women by maternal pregnancy depressive symptoms, N = 373

Cord blood DNA methylation in relation to maternal depression symptoms during pregnancy

After adjusting for potential confounding factors and controlling for multiple testing, in all group analyses, the number of DM CpGs between all study groups are presented in Table 2 and the data are displayed as quantile–quantile (QQ)-plots in Additional file 1: Fig. 1. In our main analysis, no CpG sites in cord blood were associated with maternal PND in comparison with healthy controls (HC). Further, there were no DM CpGs when comparing PND with or without history of depression to each other or HC. On the other hand, PND comorbid with anxiety showed a total of two DM CpGs in comparison with HC: at cg05385119 located just upstream of the ATP Binding Cassette Subfamily F Member 1 gene (ABCF1) (log twofold change − 0.21) and at cg00017362 (log twofold change 0.18) located upstream of Homo sapiens integrator complex subunit 10 gene (INTS10) (Table 3, Additional file 2: Fig. 2). Stratification by sex of the child did not reveal any DM between study groups (Table 2).

Table 2 Overview of the number of differentially methylated CpG sites in different group comparisons
Fig. 1
figure1

Flowchart describing selection of mothers for DNA methylation analysis from the BASIC study cohort

Additionally, PND with or without SSRI treatment was tested (Table 3, Additional file 2: Fig. 2). A total of 14 CpGs were DM among PND cases without SSRIs compared with those on SSRI treatment, while HC showed 14 CpG sites to be DM compared PND on SSRI treatment. Nine of the CpG sites identified coincided in both comparisons, meaning they were DM in both PND without SSRI and HC compared with PND on SSRIs, Fig. 2. Pathway enrichment analysis did not show any association between the DM genes and known pathways. Stratified analyses based on the sex of the child showed, among girls, 131 DM CpGs sites in PND without SSRIs compared with PND treated with SSRI, while HC showed 109 CpG sites to be DM compared with PND treated with SSRIs. In boys, PND without SSRIs compared with PND treated with SSRIs showed 12 CpGs to be DM, while 195 CpG sites were DM healthy controls in comparison with PND on SSRIs.

Table 3 Differentially methylated CpGs
Fig. 2
figure2

Venn diagram displaying genes with CpGs differently methylated in cord blood from children with mothers experiencing prenatal depressive symptoms (PND) without selective serotonin receptor inhibitor (SSRI) treatment compared with PND with SSRI treatment (blue) and healthy controls compared with PND with SSRI treatment (grey). Arrow pointing upwards indicates higher DNA methylation in the PND without SSRI group (blue) or HC (grey) in comparison with PND with SSRI treatment while arrows pointing downwards indicate lower DNA methylation

In addition to specific DM CpGs, we also identified one differentially methylated region (DMR) when comparing HC with PND on SSRI treatment. The DMR included two CpGs located on chromosome 11 but in a unknown gene. Stratified analyses based on the sex of the child showed, among girls, one DMR in PND without SSRIs compared with PND treated with SSRI including two CpGs in the C10of11 gene. The exact same DMR was also identified in HC compared with PND treated with SSRIs. In boys, one DMR was identified in HC compared with PND treated with SSRIs, including two CpGs in the Family with sequence similarity 46 (FAM46B) gene.

We identified no DM CpGs within the genes of interest (NR3C1, COL7A1, ZNF575, MEG3, SLC6A4, HSD11B2, FKBP5, BDNF, IGF2, PLAGL1, MEST, and PEG3) based on previous findings in the literature in any of the group comparisons.

Discussion

This study investigated the potential association between maternal PND with or without co-morbid anxiety, history of depression and SSRIs treatment with cord blood DNA methylation at birth. While no DM CpG site was identified in relation to PND alone, sub-group analyses showed differential DNA methylation associated with comorbidity with anxiety and SSRI treatment.

Prenatal depressive symptoms have been associated with less optimal development and associated alterations in child behavior, suggesting that it might influence fetal epigenetic programming [3]. In this study, no signs of aberrant fetal programming through epigenetic changes in cord blood were observed, since there were no DM found to be associated with PND alone. These negative findings are in agreement with most previous studies [16, 27, 32]. This could imply that PND, together with biological changes associated with this state, such as lower serum allopregnenolone [33], lower levels of peripheral anti-inflammatory markers [34] and lower levels of metabolites such as phosphate and arginine [35], do not affect DNA methylation in cord blood. However, it does not rule out other mechanisms through which PND could be involved in fetal programming effects, such as histone modifications and microRNAs, affecting gene expression. It could also be hypothesized that women with PND constitute a heterogeneous group, and that specific depressive phenotypes [8] need to be clearly outlined in order to find differences in epigenetic programming.

It has been proposed that earlier depression might have a persistent effect on biological systems even years after the depressive episode [9]. We therefore performed sub-group analyses, comparing PND alone with PND together with a history of depression. This sub-group analysis did not however show significant differences in DM based on history of depression. In contrast, the co-existence of PND with anxiety symptoms did show DM in one gene, ABCF1. Interestingly, ABCF1 has been suggested to be one of the putative therapeutic targets for escitalopram, an SSRI treatment used for both depressive symptoms and anxiety [36]. Altered methylation of ABCF1 has also been seen to be associated with autism in an in silico analysis multi-omics data analysis including frontal cortex samples [37]. Previous studies investigating the association between PND with anxiety symptoms and cord blood DNA methylation have not compared with PND alone compared with healthy controls [16, 38].

Of special interest are also the DM CpGs associated with SSRI treatment, especially located in or associated with the promoter region of specific genes. Variation in DNA methylation level in promoter regions is often associated with transcription factor binding where DNA methylation can affect the binding potential of transcriptional factors to promoter regions or transcription factor can protect the DNA from getting methylated. In both scenarios, a lower DNA methylation level is generally associated with more transcription factor binding and higher expression of the gene, higher DNA methylation in gene bodies is associated with increased gene expression. Many of the CpGs identified in the analysis in the SSRI-treated PND group are found in the promoter region of Cereblon (CRBN), MyoD Family Inhibitor Domain Containing (MDFIC), Cyclin-Dependent Kinase 11B (CDK11B), and Mitogen-Activated Protein Kinase Kinase Kinase Kinase 5 (MAP4K5) genes. The CRBN gene has been shown to be involved in intellectual disabilities [39] and has immunomodulatory functions [40]. MDFIC has been shown to be involved in the cross-talk between the glucocorticoid receptor and other proteins important for the glucocorticoid responses [41]. The expression of MDIFC is also associated with the aggressiveness and clinical outcome in neuroblastoma [42], but the mechanisms are not well understood. CDK11B has been shown to be involved in neuroblastoma [43] while MAP4K5 has been suggested to have a role in stress response [44]. Together, the DM of these genes, many related to stress and brain functions, in cord blood, might reflect altered methylation even in the fetal brain, with possible biological consequences in the child.

C10orf11 was one of the genes found to be DM in SSRI-treated PND group compared with both HC and PND without SSRIs. It was also identified as a DMR in girls in the SSRI-treated PND group. Interestingly, in rare occasions, disruption of C10ord11 has been associated with cognitive defects in children [45].

A previous study, by Cardenas and colleagues [16], also identified DM CpG sites in cord blood in relation to SSRI treatment during pregnancy; they did however not consider depressive symptoms with or without SSRI treatment. Amongst the CpG sites they identified, only one, ZNF575, was identified in their validation cohort. ZNF575 was one of the genes of interest separately analyzed even in present study, but no DM was observed in this gene in our cohort. Further, Cardenas and colleagues did not perform any stratified analysis by the sex of the child, while our analysis showed substantial differences between the sexes. Female offspring had considerably more DM sites between groups in comparison with the whole cohort and male offspring. It is however important to note that both the study by Cardenas and colleagues [16] and our study are very limited by the small sample size, and all results must be further validated in larger study cohorts.

Maternal SSRI use during pregnancy has been shown to affect development in the child. However, the effects of SSRI treatment during pregnancy are difficult to investigate without taking into account the underlying effects of depressive symptoms as reviewed by Olivier et al. 2013 [46]. We have previously shown that there are differences in gene expression in the placenta of women with PND being treated with SSRI compared with untreated women with PND [12]. Also, women with PND treated with SSRI have been shown to have higher corticotropin-releasing hormone levels associated with increased risk of pre-term birth [11] compared with untreated women with PND. Our findings add on to the knowledge that there are different biological effects between women experiencing PND with and without SSRI treatment, but it is impossible to distinguish if these are due to the treatment itself, a more severe PND demanding SSRI treatment, or a more prolonged depressive episode. Larger studies with a large number of women with severe depression without medication should be further encouraged. Unfortunately, this is not easy as it is not ethical not to treat depression, and women choosing themselves not to take medication might differ from the rest of the participants.

Strengths of the present study are the population-based prospective assessment of different sub-types of PND with EWAS analysis in cord blood. Moreover, weaknesses of previous studies were addressed in the current study through adjustment for mode of delivery [47] and cell composition known to affect DNA methylation [28, 29], as well as the exclusion of participants with pregnancy complications, also known to affect DNA methylation [48, 49]. Further, we have implemented strict control for multiple testing and stratified our analyses by sex known to affect DNA methylation [30] but also the effect by PND on the child [3]. The results should also be considered in relation to the following limitations: the relatively small sample size, albeit not uncommon in EWAS analyses. Future studies may consider joining the available cohorts investigating the effect of SSRI treatment on cord blood DNA methylation in a meta-analysis in order to gather more statistical strength and perform permutation analysis in a greater cohort. The BASIC study’s relatively low participation rate [50] makes the results not readily generalizable to the background population also depression was assessed through self-reports on a screening instrument and not a clinical diagnosis. In addition, the magnitude of differences in methylation levels was small between groups. The use of surrogate tissue is also a limitation. The cell/tissue-specificity of epigenetic patterns is a major challenge in psychological epigenetics. However, cell-specific analysis would limit research to certain brain areas of suicide victims or cadavers and limit the possibility of investigating early biomarkers affected by maternal depressive symptoms. DNA methylation in blood has been found to be in good agreement with DNA methylation in brain, especially when comparing individual CpGs [22]. Furthermore, we have highlighted in the discussion genes with relevance in stress and brain function with possible biological consequences in the child. In addition, we did not perform any adjustment for cross-reactive probes before group comparisons, which can cause false signals, possibly resulting in invalid conclusions. One of our DM probes were included on the list of genes previously identified cross-reactive probes [51], and this has now been removed from the presented list of DM genes, which should account for this limitation. Finally, this study did not investigate longitudinal follow-up of child behavior in relation to cord blood DNA methylation.

In conclusion, this study investigated potential associations between maternal PND with child epigenetic markers at birth, taking into account history of depression, comorbid anxiety or SSRI treatment, and also sex of the child. While no DM CpG site was identified in relation to PND alone, differences in DNA methylation were observed in the cord blood of children whose mothers had PND together with anxiety and also those who were treated with SSRIs. Several genes, for instance CRBN and MDFIC, previously shown to be associated with brain development and function where identified in those who were exposed to SSRIs. This study can however not distinguish between the effects of SSRI treatment alone or a possibly more severe PND demanding SSRI treatment. Furthermore, the findings must be validated in larger study cohorts including more SSRI exposed samples.

Methods

Study population

This study includes mothers and offspring participating in the BASIC Study (“Biology, Affect, Stress, Imaging and Cognition in pregnancy and puerperium”), a longitudinal study that recruited patients between 2009 and 2019 with the aim to enhance the understanding of perinatal depression [50].

All pregnant women in Uppsala County, over 18 years of age, who speaks Swedish and were scheduled for a routine ultrasound at Uppsala University Hospital, were invited to participate in the BASIC study. Women diagnosed with a pathological pregnancy (malformations leading to termination of pregnancy or miscarriage), with blood-borne disease, or with protected personal data were excluded. The participation rate has been estimated to be 22% of all women giving birth in Uppsala County [50]. Participants contributed via web-based questionnaires during gestational weeks 17 and 32 and postpartum at six weeks, six months, and twelve months. Umbilical cord blood samples were collected from participants during childbirth.

At gestational week 17, participants reported their age, education, smoking habits, employment, and country of birth (a proxy for ethnicity) by use of web-based questionnaire. Participants were also asked questions about their history of depression. Medical information regarding the pregnancy and birth was extracted from medical records.

The Edinburgh Postnatal Depression Scale (EPDS), validated in Swedish [52, 53], was included in the web-based questionnaires at pregnancy week 17 and 32. Anxiety symptoms were measured with the State-Trait Anxiety Inventory for Adults (STAI-AD, state scale) or the Beck Anxiety Inventory (BAI), depending on the time point of inclusion to the BASIC study [50].

For this sub-study, mothers were included if they had been followed throughout pregnancy and the postpartum period until six months postpartum, and a cord blood sample was collected at birth (n = 2,595), as described in Fig. 1. Exclusion criteria were women who were smoking as smoking affects epigenetic modifications. In addition, as samples were collected for additional analyses of child behavior at later in life, women presenting depressive symptoms in the postpartum period were excluded as postpartum depression could affect offspring psychoemotional development.

Only one pregnancy was included if women had participated multiple times in BASIC, and if the mother experienced depressive symptoms in one pregnancy but not the other, the one with depressive symptoms was chosen, otherwise the first pregnancy was included. Mothers defined as controls had EPDS ≤ 11 at both pregnancy week 17 and 32, had to have answered all questionnaires throughout the study and were excluded if they had anxiety defined by use of STAI, BAI or the EPDS anxiety subscale, any history of depression, or had missing information from medical records. Women with depressive symptoms during pregnancy were identified by the use of EPDS > 12 at pregnancy week 17, pregnancy week 32 or both. Final analyses also excluded women with pregnancy complications such as pregnancy hypertension, diabetes, itching, placenta previa, anemia, hepatosis, cystitis, preeclampsia, or HELLP syndrome. This resulted in 373 pregnancies included in the final analysis (Fig. 1).

Mothers were divided into two main study groups, mothers with PND (n = 178) or healthy controls (HC) (n = 195). Sub-groups analyses were also performed between women with PND with (n = 132) or without history of depression (n = 46), in comparison with HC, as well as PND comorbid with anxiety (n = 101) or not (n = 77), in comparison with HC. Finally, also based on PND with (n = 13) or without SSRI treatment (n = 168) during pregnancy, in comparison with HC.

Cord blood collection and DNA extraction

In all steps, from DNA extraction to DNA methylation analysis, PND and HC samples were randomized over batches and arrays. Whole blood from the cord vein was collected at childbirth and stored at − 70 °C until DNA extraction by use of QIAamp DNA blood/mini kit (Qiagen) according to manufacturer’s instructions. DNA concentrations and purity were measured by use of NanoDrop Spectrophotometer (Thermo Scientific).

DNA methylation data

Bisulfite conversion of DNA was performed by use of the Zymo Research DNA Methylation-Lightning Kit (Zymo Research). Genome-wide methylation at 850 000 CpG sites was subsequently analyzed with the Illumina MethylationEPIC kit (Illumina). All reactions were performed according to the manufacturers' instructions at Life and Brain Life & Brain Research Centre University Hospital of Bonn, Germany. Data processing was carried out by National Bioinformatics Infrastructure Sweden (NBIS) in R 3.5.2 using R package minfi 1.28.4 [54]. The workflow followed was as per recommended by the tool guide. Detection p-value (DPV) with a cut-off at 0.01 was used to remove CpGs with poor signal (4.5% of the CpGs had DPV > 0.01 in 10% or more of the samples). No samples, only individual CpGs were removed based on DPV. An additional 3.2% of the CpGs were discarded due to presence of known SNPs in the CpG site or single-base extension. SNP annotation data were obtained through minfi:getSnpInfo, which retrieves SNP data from IlluminaHumanMethylationEPICanno.ilm10b2.hg19 version 0.6.0. All sites with SNPs were removed as recommended by the supplier of the Illumina MethylationEPIC kit (Illumina).

Functions in the R package minfi were used for quantile normalization suitable for datasets where global differences are not expected between samples [55]. Sex information was provided to the function rgf_quantile  ← preprocessQuantile(rgf,sex = sx), beta_quantile ← getBeta(rgf_quantile), m_quantile ← get_mv(beta_quantile). This implements a stratified quantile normalization which is applied to the methylated and unmethylated signal intensities separately and takes into account the different probe types [56]. This to correct for non-biological variation between CpGs and the data was thereafter transformed to beta, representing the proportion of cells methylated at the each CpG site (from 0 to 1 methylation) and M values for statistical analysis. B-values were chosen for visualization as recommended in the literature [57]. CpGs from X and Y chromosomes were discarded and only autosomal CpGs were used for downstream analyses, according to the tool guide instructions. Cell composition in each sample was estimated using R package FlowSorted.CordBlood.450 k [29] as the reference data.

Statistical analyses

PCA was used for exploratory analyses to investigate general structure and patterns in the data. Beta values were used for PCA. The 100,000 CpGs with the highest variance were selected for the PCA analyses. Up to 17% of the total variation in the dataset was explained by the first three principal components (PCs). Benjamini–Hochberg corrected Kruskal–Wallis p-values were used a measure to assess association of top 6 PCs with categorical variables of interest. Benjamini–Hochberg corrected Spearman correlation p-values were used a measure to assess association of top 6 PCs with continuous variables of interest.

Differential methylation (DM) CpGs between the study groups and sub-groups were identified using the R package Limma. The data were also stratified by sex of the child and the analysis was also run excluding participants using SSRI. For each dataset, DM tests were run to assess how many CpGs differed significantly in methylation between various groups. M values were used for DM analysis. For each test, a linear model was built including the variable of interest along with 12 covariates, which were controlled for in the model. The following covariates were used in all the models: sample plate batch, sentrix position on the chip, granulocyte cell estimate, red blood cell estimate, CD8T cell estimate. Also clinical characteristics such as maternal age at delivery, BMI before pregnancy, maternal place of birth (Scandinavia or other), educational level (University degree yes or no), mode of delivery (vaginal, vacuum extraction, emergency cesarean section, and elective cesarean section), gestational age, and parity, based on previous literature [16, 47, 58, 59]. Group differences in DNA methylation in specific CpGs are presented as mean and standard deviation and also log2 fold change (log2 (M-value group1/M-value group 2)) between the groups, meaning that a negative log2 fold changes means lower DNA methylation in group 1. Pathway enrichment analysis on DM genes was performed by use of the FUMA tool by using the mode GENE2FUNC [60].

Differential methylation over regions was run using R package DMRcate 2.4.1 in R 4.0.2 [61]. The parameters (lambda = 500, C = 5) were set as recommended in Mallik et al., 2019 [62]. An FDR cut-off of 0.05 was used to determine significance.

In addition to the overall comparison of all DM between groups, CpGs in genes of interest identified in the literature [16, 18, 25, 26] and were included in a separate analysis, in which the adjustment of p-value for multiple comparisons was based on the number of CpGs in the genes of interest, instead of all included in the EPIC array. These genes of interest were: Nuclear Receptor Subfamily 3 Group C Member 1 (NR3C1), Collagen Type VII Alpha 1 Chain (COL7A1), Zinc Finger Protein 575 (ZNF575), Maternally Expressed Gene 3 (MEG3), Solute Carrier Family 6 Member 4 (SLC6A4), Hydroxysteroid 11-Beta Dehydrogenase 2 (HSD11B2), FKBP Prolyl Isomerase 5 (FKBP5), Brain Derived Neurotrophic Factor (BDNF), Insulin Like Growth Factor 2 (IGF2), PLAG1 Like Zinc Finger 1 (PLAGL1), Mesoderm Specific Transcript (MEST) and Paternally Expressed 3 (PEG3) [16, 18, 25, 26].

Availability of data and materials

The data that support the findings of this study are available from the authors upon reasonable request and with permission of BASIC study cohort.

References

  1. 1.

    Szegda K, Markenson G, Bertone-Johnson ER, Chasan-Taber L. Depression during pregnancy: a risk factor for adverse neonatal outcomes? A critical review of the literature. J Matern Fetal Neonatal Med. 2014;27(9):960–7.

    PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Stein A, Pearson RM, Goodman SH, Rapa E, Rahman A, McCallum M, et al. Effects of perinatal mental disorders on the fetus and child. The Lancet. 2014;384(9956):1800–19.

    Article  Google Scholar 

  3. 3.

    Fransson E, Sorensen F, Kunovac Kallak T, Ramklint M, Eckerdal P, Heimgartner M, et al. Maternal perinatal depressive symptoms trajectories and impact on toddler behaviour—the importance of symptom duration and maternal bonding. J Affect Disord. 2020;273:542–51.

    PubMed  Article  PubMed Central  Google Scholar 

  4. 4.

    Dean DC III, Planalp EM, Wooten W, Kecskemeti SR, Adluru N, Schmidt CK, et al. Association of prenatal maternal depression and anxiety symptoms with infant white matter microstructure. JAMA Pediatrics. 2018;172(10):973–81.

    PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Hay DF, Pawlby S, Waters CS, Sharp D. Antepartum and postpartum exposure to maternal depression: different effects on different adolescent outcomes. J Child Psychol Psychiatry. 2008;49(10):1079–88.

    PubMed  Article  PubMed Central  Google Scholar 

  6. 6.

    Wolford E, Lahti M, Tuovinen S, Lahti J, Lipsanen J, Savolainen K, et al. Maternal depressive symptoms during and after pregnancy are associated with attention-deficit/hyperactivity disorder symptoms in their 3- to 6-year-old children. PLoS ONE. 2017;12(12):e0190248.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  7. 7.

    Goodman JH, Chenausky KL, Freeman MP. Anxiety disorders during pregnancy: a systematic review. J Clin Psychiatry. 2014;75(10):e1153–84.

    PubMed  Article  PubMed Central  Google Scholar 

  8. 8.

    Wikman A, Axfors C, Iliadis SI, Cox J, Fransson E, Skalkidou A. Characteristics of women with different perinatal depression trajectories. J Neurosci Res. 2019.

  9. 9.

    Meltzer-Brody S. New insights into perinatal depression: pathogenesis and treatment during pregnancy and postpartum. Dialogues Clin Neurosci. 2011;13(1):89–100.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Putnam KT, Wilcox M, Robertson-Blackmore E, Sharkey K, Bergink V, Munk-Olsen T, et al. Clinical phenotypes of perinatal depression and time of symptom onset: analysis of data from an international consortium. Lancet Psychiatry. 2017;4(6):477–85.

    PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Hannerfors AK, Hellgren C, Schijven D, Iliadis SI, Comasco E, Skalkidou A, et al. Treatment with serotonin reuptake inhibitors during pregnancy is associated with elevated corticotropin-releasing hormone levels. Psychoneuroendocrinology. 2015;58(Supplement C):104–13.

  12. 12.

    Edvinsson Å, Hellgren C, Kunovac Kallak T, Åkerud H, Skalkidou A, Stener-Victorin E, et al. The effect of antenatal depression and antidepressant treatment on placental tissue: a protein-validated gene expression study. BMC Pregnancy Childbirth. 2019;19(1):479.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Bekkhus M, Lee Y, Nordhagen R, Magnus P, Samuelsen SO, Borge AIH. Re-examining the link between prenatal maternal anxiety and child emotional difficulties, using a sibling design. Int J Epidemiol. 2018;47(1):156–65.

    PubMed  Article  PubMed Central  Google Scholar 

  14. 14.

    Fransson E, Dubicke A, Byström B, Ekman-Ordeberg G, Hjelmstedt A, Lekander M. Negative Emotions and Cytokines in Maternal and Cord Serum at Preterm Birth. Am J Reprod Immunol. 2011;67(6):506–14.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  15. 15.

    Weaver ICG, Cervoni N, Champagne FA, D’Alessio AC, Sharma S, Seckl JR, et al. Epigenetic programming by maternal behavior. Nat Neurosci. 2004;7(8):847–54.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  16. 16.

    Cardenas A, Faleschini S, Cortes Hidalgo A, Rifas-Shiman SL, Baccarelli AA, DeMeo DL, et al. Prenatal maternal antidepressants, anxiety, and depression and offspring DNA methylation: epigenome-wide associations at birth and persistence into early childhood. Clinical epigenetics. 2019;11(1):56.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  17. 17.

    Devlin AM, Brain U, Austin J, Oberlander TF. Prenatal exposure to maternal depressed mood and the MTHFR C677T variant affect SLC6A4 methylation in infants at birth. PLoS ONE. 2010;5(8):e12201.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  18. 18.

    Oberlander TF, Weinberg J, Papsdorf M, Grunau R, Misri S, Devlin AM. Prenatal exposure to maternal depression, neonatal methylation of human glucocorticoid receptor gene (NR3C1) and infant cortisol stress responses. Epigenetics. 2008;3(2):97–106.

    PubMed  Article  PubMed Central  Google Scholar 

  19. 19.

    Meaney MJ, Szyf M. Environmental programming of stress responses through DNA methylation: life at the interface between a dynamic environment and a fixed genome. Dialogues Clin Neurosci. 2005;7(2):103–23.

    PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Hannon E, Spiers H, Viana J, Pidsley R, Burrage J, Murphy TM, et al. Methylation QTLs in the developing brain and their enrichment in schizophrenia risk loci. Nat Neurosci. 2016;19(1):48–54.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Davies MN, Volta M, Pidsley R, Lunnon K, Dixit A, Lovestone S, et al. Functional annotation of the human brain methylome identifies tissue-specific epigenetic variation across brain and blood. Genome Biol. 2012;13(6):R43.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    Braun PR, Han S, Hing B, Nagahama Y, Gaul LN, Heinzman JT, et al. Genome-wide DNA methylation comparison between live human brain and peripheral tissues within individuals. Translational Psychiatry. 2019;9(1):47.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  23. 23.

    Bauer AZ, Kriebel D, Herbert MR, Bornehag C-G, Swan SH. Prenatal paracetamol exposure and child neurodevelopment: A review. Horm Behav. 2018;101:125–47.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Gervin K, Nordeng H, Ystrom E, Reichborn-Kjennerud T, Lyle R. Long-term prenatal exposure to paracetamol is associated with DNA methylation differences in children diagnosed with ADHD. Clin Epigenet. 2017;9(1):77.

    Article  CAS  Google Scholar 

  25. 25.

    Cao-Lei L, de Rooij SR, King S, Matthews SG, Metz GAS, Roseboom TJ, et al. Prenatal stress and epigenetics. Neuroscience & Biobehavioral Reviews. 2017.

  26. 26.

    Mendonça MS, Mangiavacchi PM, De Sousa PF, Crippa JAS, Mendes AV, Loureiro SR, et al. Epigenetic variation at the SLC6A4 gene promoter in mother–child pairs with major depressive disorder. J Affect Disord. 2019;245:716–23.

    PubMed  Article  CAS  Google Scholar 

  27. 27.

    Viuff AC, Sharp GC, Rai D, Henriksen TB, Pedersen LH, Kyng KJ, et al. Maternal depression during pregnancy and cord blood DNA methylation: findings from the Avon Longitudinal Study of Parents and Children. Transl Psychiatry. 2018;8(1):244.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Gervin K, Page CM, Aass HCD, Jansen MA, Fjeldstad HE, Andreassen BK, et al. Cell type specific DNA methylation in cord blood: A 450K-reference data set and cell count-based validation of estimated cell type composition. Epigenetics. 2016;11(9):690–8.

    PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Bakulski KM, Feinberg JI, Andrews SV, Yang J, Brown S, L. McKenney S, et al. DNA methylation of cord blood cell types: Applications for mixed cell birth studies. Epigenetics. 2016;11(5):354–62.

  30. 30.

    Yousefi P, Huen K, Davé V, Barcellos L, Eskenazi B, Holland N. Sex differences in DNA methylation assessed by 450 K BeadChip in newborns. BMC Genomics. 2015;16(1):911.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  31. 31.

    Viuff AC, Pedersen LH, Kyng K, Staunstrup NH, Borglum A, Henriksen TB. Antidepressant medication during pregnancy and epigenetic changes in umbilical cord blood: a systematic review. Clin Epigenetics. 2016;8(1):94.

    PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Rijlaarsdam J, Pappa I, Walton E, Bakermans-Kranenburg MJ, Mileva-Seitz VR, Rippe RCA, et al. An epigenome-wide association meta-analysis of prenatal maternal stress in neonates: A model approach for replication. Epigenetics. 2016;11(2):140–9.

    PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Hellgren C, Åkerud H, Skalkidou A, Sundström-Poromaa I. Cortisol awakening response in late pregnancy in women with previous or ongoing depression. Psychoneuroendocrinology. 2013;38(12):3150–4.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  34. 34.

    Edvinsson Å, Bränn E, Hellgren C, Freyhult E, White R, Kamali-Moghaddam M, et al. Lower inflammatory markers in women with antenatal depression brings the M1/M2 balance into focus from a new direction. Psychoneuroendocrinology. 2017;80:15–25.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  35. 35.

    Henriksson HE, Malavaki C, Brann E, Drainas V, Lager S, Iliadis SI, et al. Blood plasma metabolic profiling of pregnant women with antenatal depressive symptoms. Transl Psychiatry. 2019;9(1):204.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  36. 36.

    Powell TR, Tansey KE, Breen G, Farmer AE, Craig IW, Uher R, et al. ATP-binding cassette sub-family F member 1 (ABCF1) is identified as a putative therapeutic target of escitalopram in the inflammatory cytokine pathway. J Psychopharmacol. 2013;27(7):609–15.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  37. 37.

    Sun Y, Yao X, March ME, Meng X, Li J, Wei Z, et al. Target genes of autism risk loci in brain frontal cortex. Front Genet. 2019;10:707.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Non AL, Binder AM, Kubzansky LD, Michels KB. Genome-wide DNA methylation in neonates exposed to maternal depression, anxiety, or SSRI medication during pregnancy. Epigenetics. 2014;9(7):964–72.

    PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    McDonough PG, Menkveld R, Kruger TF. Sperm morphology—predictive value? Fertil Steril. 1992;57(4):942–3.

    Article  Google Scholar 

  40. 40.

    Shi Q, Chen L. Cereblon: a protein crucial to the multiple functions of immunomodulatory drugs as well as cell metabolism and disease generation. Journal of immunology research. 2017;2017.

  41. 41.

    Oakley RH, Busillo JM, Cidlowski JA. Cross-talk between the glucocorticoid receptor and MyoD family inhibitor domain-containing protein provides a new mechanism for generating tissue-specific responses to glucocorticoids. J Biol Chem. 2017;292(14):5825–44.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Khan FH, Pandian V, Ramraj S, Natarajan M, Aravindan S, Herman TS, et al. Acquired genetic alterations in tumor cells dictate the development of high-risk neuroblastoma and clinical outcomes. BMC Cancer. 2015;15(1):514.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  43. 43.

    White PS, Maris JM, Beltinger C, Sulman E, Marshall HN, Fujimori M, et al. A region of consistent deletion in neuroblastoma maps within human chromosome 1p36.2–36.3. Proc Natl Acad Sci U S A. 1995;92(12):5520–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Shi CS, Tuscano JM, Witte ON, Kehrl JH. GCKR links the Bcr-Abl oncogene and Ras to the stress-activated protein kinase pathway. Blood. 1999;93(4):1338–45.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  45. 45.

    Tzschach A, Bisgaard AM, Kirchhoff M, Graul-Neumann LM, Neitzel H, Page S, et al. Chromosome aberrations involving 10q22: report of three overlapping interstitial deletions and a balanced translocation disrupting C10orf11. Eur J Hum Genet. 2010;18(3):291–5.

    PubMed  Article  PubMed Central  Google Scholar 

  46. 46.

    Olivier J, Åkerud H, Kaihola H, Pawluski J, Skalkidou A, Högberg U, et al. The effects of maternal depression and maternal selective serotonin reuptake inhibitor exposure on offspring. Frontiers in Cellular Neuroscience. 2013;7(73).

  47. 47.

    Almgren M, Schlinzig T, Gomez-Cabrero D, Gunnar A, Sundin M, Johansson S, et al. Cesarean delivery and hematopoietic stem cell epigenetics in the newborn infant: implications for future health? American Journal of Obstetrics and Gynecology. 2014;211(5):502.e1–.e8.

  48. 48.

    Kazmi N, Sharp GC, Reese SE, Vehmeijer FO, Lahti J, Page CM, et al. Hypertensive disorders of pregnancy and DNA methylation in newborns: findings from the pregnancy and childhood epigenetics consortium. Hypertension. 2019;74(2):375–83.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Herzog EM, Eggink AJ, Willemsen SP, Slieker RC, Wijnands KPJ, Felix JF, et al. Early- and late-onset preeclampsia and the tissue-specific epigenome of the placenta and newborn. Placenta. 2017;58:122–32.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  50. 50.

    Axfors C, Bränn E, Henriksson HE, Hellgren C, Kunovac Kallak T, Fransson E, et al. Cohort profile: the Biology, Affect, Stress, Imaging and Cognition (BASIC) study on perinatal depression in a population-based Swedish cohort. BMJ open. 2019;9(10):e031514.

    PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Pidsley R, Zotenko E, Peters TJ, Lawrence MG, Risbridger GP, Molloy P, et al. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 2016;17(1):208.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  52. 52.

    Cox JL, Holden JM, Sagovsky R. Detection of postnatal depression. Development of the 10-item Edinburgh Postnatal Depression Scale. Br J Psychiatry. 1987;150:782–6.

  53. 53.

    Wickberg B, Hwang CP. Screening for postnatal depression in a population-based Swedish sample. Acta Psychiatr Scand. 1997;95(1):62–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  54. 54.

    Fortin JP, Triche TJ Jr, Hansen KD. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics. 2017;33(4):558–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Fortin JP, Labbe A, Lemire M, Zanke BW, Hudson TJ, Fertig EJ, et al. Functional normalization of 450k methylation array data improves replication in large cancer studies. Genome Biol. 2014;15(12):503.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  56. 56.

    Touleimat N, Tost J. Complete pipeline for Infinium((R)) Human Methylation 450K BeadChip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics. 2012;4(3):325–41.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  57. 57.

    Du P, Zhang X, Huang CC, Jafari N, Kibbe WA, Hou L, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010;11:587.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  58. 58.

    Barker ED, Walton E, Cecil CAM. Annual Research Review: DNA methylation as a mediator in the association between risk exposure and child and adolescent psychopathology. J Child Psychol Psychiatry. 2018;59(4):303–22.

    PubMed  Article  PubMed Central  Google Scholar 

  59. 59.

    Merid SK, Novoloaca A, Sharp GC, Kupers LK, Kho AT, Roy R, et al. Epigenome-wide meta-analysis of blood DNA methylation in newborns and children identifies numerous loci related to gestational age. Genome Med. 2020;12(1):25.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  60. 60.

    Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. 2017;8(1):1826.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  61. 61.

    Peters TJ, Buckley MJ, Statham AL, Pidsley R, Samaras K, R VL, et al. De novo identification of differentially methylated regions in the human genome. Epigenetics Chromatin. 2015;8:6.

  62. 62.

    Mallik S, Odom GJ, Gao Z, Gomez L, Chen X, Wang L. An evaluation of supervised methods for identifying differentially methylated regions in Illumina methylation arrays. Brief Bioinform. 2019;20(6):2224–35.

    CAS  PubMed  Article  Google Scholar 

Download references

Acknowledgements

This project would not have been possible without the participating families and the fantastic staff at Uppsala University Hospital maternity ward. All data processing was performed at NBIS and the authors would like to acknowledge Roy Francis and Eva Freyhult for their contribution to the project. The authors would also like to acknowledge support from Science for Life Laboratory and the National Genomics Infrastructure. The computations were performed on resources provided by SNIC through Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX) under Project sens2019515, partially funded by the Swedish Research Council through grant agreement no. 2018-05973.

Funding

This project was funded by The Swedish Brain Foundation, the Göran Gustafsson Foundation, the Swedish Research Council (Grant Nos. 523-2014-2342, 523-2014-07605, and 521-2013-2339) to AS, the Gillbergska Foundation, and the Fredrik and Ingrid Thurings foundation to TKK, E.C. receives funds from SciLifeLab and the Märta Lundqvist foundation.

Author information

Affiliations

Authors

Contributions

TKK, EF, and AS conceived and designed the study. TKK and SL performed parts of the experiments. TKK, EB, EF, ÅJ, SL, EC, RL, and AS analyzed the data and performed statistical analysis in collaboration with NBIS. TKK, EB, EF, and AS wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Theodora Kunovac Kallak.

Ethics declarations

Ethics approval and consent to participate

The study was conducted at the Department of Women’s and Children’s Health, Uppsala University and the Women’s clinic at the Uppsala University Hospital. Ethical approval was granted from the Uppsala Ethical Review Board reference numbers, BASIC 2009/171 with amendments. Written informed consent was obtained from all participants.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1. Figure S1:

Quantile-quantile plot (QQ-plot) displaying the observed and expected Benjamini-Hochberg adjusted –log10 p-values for all group comparisons. a) Healthy control versus Prenatal depressive symptoms. b) History and prenatal depressive symptoms versus Prenatal depressive symptoms only. c) History and prenatal depressive symptoms versus Healthy control. d) Prenatal depressive symptoms only versus Healthy control. e) Anxiety and prenatal depressive symptoms versus Prenatal depressive symptoms only. f) Anxiety and prenatal depressive symptoms versus Healthy control. g) Prenatal depressive symptoms only versus Healthy control. h) Healthy control versus Prenatal depressive symptoms without selective serotonin reuptake inhibitors (SSRI). i) Healthy control versus Prenatal depressive symptoms with SSRI. j) Prenatal depressive symptoms without SSRI versus Prenatal depressive symptoms with SSRI.

Additional file 2. Figure S2:

Manhattan plots displaying differentially DNA methylated genes in cord blood from infants born by mothers who a) suffer from anxiety and prenatal depressive symptoms compared with healthy controls, b) healthy controls compared with prenatal depressive symptoms treated with selective serotonin reuptake inhibitors (SSRIs), and c) women with untreated prenatal depressive symptoms compared with prenatal depressive symptoms treated with SSRIs. Y-axes display the Benjamini-Hochberg adjusted –log10 p-values for specific group comparisons and x-axes shows the chromosomal location.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Kallak, T.K., Bränn, E., Fransson, E. et al. DNA methylation in cord blood in association with prenatal depressive symptoms. Clin Epigenet 13, 78 (2021). https://doi.org/10.1186/s13148-021-01054-0

Download citation

Keywords

  • CBCL
  • Cord blood
  • DNA methylation
  • EPDS
  • EWAS
  • Prenatal depressive symptoms
  • Women