DNA methylation associated with postpartum depressive symptoms overlaps findings from a genome-wide association meta-analysis of depression
Clinical Epigenetics volume 11, Article number: 169 (2019)
Perinatal depressive symptoms have been linked to adverse maternal and infant health outcomes. The etiology associated with perinatal depressive psychopathology is poorly understood, but accumulating evidence suggests that understanding inter-individual differences in DNA methylation (DNAm) patterning may provide insight regarding the genomic regions salient to the risk liability of perinatal depressive psychopathology.
Genome-wide DNAm was measured in maternal peripheral blood using the Infinium MethylationEPIC microarray. Ninety-two participants (46% African-American) had DNAm samples that passed all quality control metrics, and all participants were within 7 months of delivery. Linear models were constructed to identify differentially methylated sites and regions, and permutation testing was utilized to assess significance. Differentially methylated regions (DMRs) were defined as genomic regions of consistent DNAm change with at least two probes within 1 kb of each other. Maternal age, current smoking status, estimated cell-type proportions, ancestry-relevant principal components, days since delivery, and chip position served as covariates to adjust for technical and biological factors. Current postpartum depressive symptoms were measured using the Edinburgh Postnatal Depression Scale. Ninety-eight DMRs were significant (false discovery rate < 5%) and overlapped 92 genes. Three of the regions overlap loci from the latest Psychiatric Genomics Consortium meta-analysis of depression.
Many of the genes identified in this analysis corroborate previous allelic, transcriptomic, and DNAm association results related to depressive phenotypes. Future work should integrate data from multi-omic platforms to understand the functional relevance of these DMRs and refine DNAm association results by limiting phenotypic heterogeneity and clarifying if DNAm differences relate to the timing of onset, severity, duration of perinatal mental health outcomes of the current pregnancy or to previous history of depressive psychopathology.
Perinatal depressive symptoms can occur any time during pregnancy or shortly following birth and even subclinical levels of depressive symptoms have been associated with an increased risk for episodes of major depression with onset in the peripartum (MDP), pregnancy complications, maternal suicide, and adverse infant health outcomes and development [1,2,3,4,5,6,7]. The Diagnostic and Statistical Manual (fifth edition; DSM-5) classifies MDP as a major depressive episode that occurs during pregnancy or within 4 weeks of delivery ; however, in practice, researchers and clinicians may extend this period to up to one year postpartum. Attempts to understand the impact of depressive psychopathology on biological mechanisms have primarily centered on elucidating the relationship between maternal mental health and negative infant outcomes. As a result, relatively little is known regarding how perinatal depressive symptoms affect maternal biological processes. Epidemiological studies suggest that episodes of major depression (MD) may increase the risk for other adverse health outcomes, such as cardiovascular disease, and perturb immune system activities. The biological pathways associated with these persistent changes in immune activity and elevated risk have not yet been identified. More work is needed to uncover the pathways underlying these comorbidities and to determine if the biological impact from depressive symptoms differs by clinical subtype (e.g., perinatal and early-onset). One potential avenue for understanding biological changes associated with depressive psychopathology is through DNA methylation (DNAm) studies.
DNAm is a chemical modification typically found on cytosines bordering guanines (i.e., cytosine-phosphate-guanine [CpG] sites). DNAm can influence gene expression, genomic stability, and chromatin conformation . Inter-individual differences in DNAm have been associated with early mortality , cancer , imprinting disorders [12, 13], childhood trauma exposure , biological age [15, 16], and schizophrenia [17, 18]. Associations between DNAm and clinical MD and/or depressive symptoms proximal to [19,20,21,22] or absent pregnancy [23,24,25,26] have been reported; however, most of the perinatal depression studies have focused on identifying DNAm patterns in fetal tissues associated with maternal mental health, leaving a significant knowledge gap. Characterizing the relationship between DNAm and perinatal depressive symptoms may provide insight regarding the pathoetiology of MDP as well as potentially identify biological markers .
This study sought to identify DNAm patterns associated with perinatal depressive symptoms during the first 7 months postpartum in maternal blood focusing on regional DNAm changes. Genome-wide DNAm and repeated measures of maternal mental health were collected as part of a longitudinal study of preterm birth . The rationale for focusing on differentially methylated regions (DMRs) rather than single CpG site associations was twofold. One, regional changes are thought to represent differences more likely to be biologically meaningful and statistically credible . A single CpG site associated with a trait could be spurious; however, multiple CpG sites within one region each associating with a trait in the same direction is likely to represent a more robust finding. Two, regional analyses reduce the burden of multiple tests and allow one to test for smaller probe effect sizes . The sample size for this study is acceptable to test regional differences, but it is not well-powered to identify individual probe associations.
The data for this analysis come from the Pregnancy, Race, Environment, Genes (PREG) study and its postpartum extension . Both PREG and the extension received IRB approval, and all participants provided written informed consent for both parts. Most participants had the first postpartum visit within 3 months of delivery and the second visit within 9 months. The PREG study recruited an epidemiological sample of women in early pregnancy primarily from two health clinics in Richmond, Virginia. The purpose of the PREG study was to identify factors that influence racial health disparities in preterm birth. The postpartum study extension permitted additional perinatal outcomes, like perinatal depression, to be investigated.
Study eligibility criteria
PREG study enrollment criteria required participants to (1) be < 24 weeks gestation, (2) have a singleton pregnancy, (3) have not used artificial reproductive technology for the current pregnancy, (4) be between 18 and 40 years old, and (5) be absent of major health conditions (e.g., diabetes) . Additionally, the participant and the biological father had to self-identify as either both African-American or both European-American and without Middle Eastern or Hispanic ancestry. Birth exclusion criteria included chromosomal or amniotic abnormalities (e.g., polyhydramnios/oligohydramnios) and congenital birth defects.
Current perinatal depressive symptoms were measured at both postpartum study visits using the Edinburgh Postnatal Depression Scale (EPDS) . The EPDS is a 10-item self-report instrument that is used in clinical practice to screen for probable perinatal depression and in research to assess for perinatal depressive symptoms; it differs from other depressive symptom measures because the items focus on symptoms specifically related to depression that would not be part of a typical pregnancy . For example, difficulty sleeping is common in pregnancies regardless of maternal MD status. EPDS total score was analyzed as a continuous variable. Scores of 13 or more out of 30 indicate probable perinatal depression .
Genome-wide DNAm measurement and processing
Maternal genome-wide DNAm was assayed according to the manufacturer’s protocol (Illumina, San Diego, CA, USA) from peripheral blood using the EPIC beadchip, which includes more than 850,000 probes and interrogates regulatory, genic, and intergenic regions . Blood samples were collected in EDTA tubes at each study visit along with health questionnaires. An aliquot of 1 μg DNA per participant was sent to HudsonAlpha Laboratories for DNAm measurement. Samples were randomized to arrays to minimize potential batch effects related to processing influencing DNAm patterns. Peripheral blood was selected given its accessibility, the availability of cell-type deconvolution methods, and the strong evidence of immune system involvement in MDP pathophysiology [34,35,36,37].
Raw microarray data was processed using Bioconductor packages in the R environment in line with best practices [38,39,40] (version 3.5). Signal intensity and probe failure rate were evaluated using the minfi  package to identify poor quality samples and probes. Samples were removed if either the median unmethylated or methylated signal intensities were less than 10.5. Probes were removed if they failed in > 1% of samples (n = 12,557), overlapped single-nucleotide polymorphisms (SNPs; n = 30,435), or had been identified as cross-hybridizing in the Illumina Infinium HumanMethylation450 beadchip (predecessor technology) . Probes on the sex chromosomes were retained given that the entire sample was female, leaving a total of 782,884 probes. All samples were quantile normalized, and blood cell type proportions were estimated using the Houseman method . Sample identity was confirmed using the 59 control probes on the EPIC beadchip. These probes overlap polymorphic sites and in aggregate can estimate sample relatedness and detect sample duplication. All pairwise sample correlations were calculated. Any sample correlated too poorly with its sister samples (r < 0.8) or too highly with samples from another person (r > 0.6) were removed or relabeled if the correct identity could be ascertained. Only one blood sample per person was used for this analysis. In general, DNAm samples from the first postpartum visit were used; however, if a participant’s first postpartum visit failed quality control and she had a second postpartum DNAm sample within 7 months of delivery, then that sample and the EPDS questionnaire from the same visit was used. For all participants, the phenotypic and DNAm data were from the same study visit. The purpose for using only one sample was to capture the postpartum time period with the highest estimates of depressive symptom prevalence while also limiting potential phenotypic heterogeneity [44,45,46].
Maternal age, number of days postpartum at blood and questionnaire collection, smoking status, and microarray row were selected as covariates a priori based on known or putative associations with DNAm. Principal component analysis (PCA) was applied to the normalized methyl values  of the filtered probe set, and correlations between the top ten principal components (PC) and technical and biological variables were plotted to identify additional potential confounders. Four variables (estimated granulocyte proportion, chip size, self-reported census-based race category, and slide ID number) correlated greater than the absolute value of 0.5 with at least one PC. Slide effects were addressed using ComBat . Genetic ancestry relevant PCs were calculated with the Barfield method . Two PCs strongly independently correlated with self-identified race were included instead of self-identified race because the PCs have been shown to adjust for genomic inflation better than a categorical variable . These two components also correlated strongly with estimated blood cell proportions (r > than the absolute value of 0.5). Estimated granulocyte proportion also was included as a covariate to provide additional adjustment for cell-type heterogeneity.
Identifying single-site and regional DNAm changes
The methyl values for individual probes were regressed onto postpartum EPDS total score using the limma package . Covariates controlling for row, granulocyte proportion, median unmethylated signal intensity maternal age, gestational age at blood collection, smoking status, and allelic ancestry were included. The strength of the association between individual probes and EPDS total score was evaluated using empirical p values derived from a null distribution derived from the dataset itself (k = 20,000 permutations). This analysis strategy provides a more appropriate assessment and adjustment for test statistic inflation than a metric like lambda . The median effect size for probes used in the DMR analysis was assessed using the difference in adjusted R2 values between the full model and a reduced model without EPDS total score.
The probes with the largest observed t statistics (top and bottom 2.5% of tested probes) were used for regional analysis. Differentially methylated regions (DMRs) were defined as contiguous regions of consistent DNAm change (i.e., all hypermethylated or all hypomethylated) that contained at least two probes within 1 kb using a method similar to that described by Ong and Holbrook . This strategy is similar to the DMRcate algorithm in that only the subset of the probes with the best evidence for association is used to construct DMRs .
DMR significance was assessed using a rank-based permutation strategy and an empirical null distribution derived from k = 1000 permutations. DMRs were constructed from both the observed data and 1000 of the 20,000 DMP permutation sets. For each DMR, the area under the curve (AUC) was calculated using the trapezoidal rule, where each probe’s t statistic served as height and the distance between the probes as width. Thus, the magnitude of the AUC reflects both the strength of each probe’s association (height) and the size of the region (width). The significance of analysis microarray (SAM) method was implemented to assign test statistics to each observed DMR . This method ranks all DMRs generated within a permutation by AUC and performs row-wise comparisons between the ranked observed DMRs and the ranked permutation DMRs to calculate a false discovery rate (FDR).
Gene set enrichment and comparison to other genetic findings related to depression
DMPs and DMRs were annotated using AnnotationHub . Gene set enrichment testing for functional and regulatory roles was performed on the combined dataset of DMPs and DMRs using Entrez IDs and clusterProfiler . The rationale for combining DMRs and DMPs into a single group for gene set analysis was to address the issue that not all probes are capable of forming DMRs. In order to give those regions of the genome an opportunity to contribute to gene set enrichment analysis, DMPs and DMRs were analyzed together (see Additional file 1 for full enrichment analysis of the combined DMR-DMP analysis as well as the DMP-only and DMR-only enrichment analyses performed separately with clusterProfiler  and methylGSA ).
The results from this analysis were compared directly to two studies of depression: the latest Psychiatric Genomics Consortium (PGC) meta-analysis of genome-wide association studies of depression and an epigenome-wide association study (EWAS) of early-onset MD. For the PGC study, the 44 significant loci were obtained to determine the extent of overlap with significant DNAm regions . Bootstrap and permutation methods (k = 1000) were used to test if DNAm regions were enriched for PGC loci. For the early-onset MD EWAS, site, regional, and gene enrichment results from the Adolescent and Young Adult Twin Study (AYATS) were compared to findings from this study to determine the extent of overlap and similarity [26, 85].
Sample demographics can be found in Table 1 and are representative of Richmond, Virginia. Approximately half of the participants (46%) self-identified as African-American, and the mean gestational age was 274.5 days (sd = 13.0 days). Most of the women (65%) were primigravida and had full-term pregnancies (94%). Very few participants were current smokers (8%), and 18% of the total sample self-reported a positive lifetime history of MD. Lifetime MD history did not significantly predict postpartum EPDS score, and the distribution of EPDS scores by self-reported MD history are in Additional file 5: Figure S3. The average EPDS score was 5.4, and 13 participants (14%) scored 13 or greater on the EPDS. The average time between delivery and postpartum study visit was 57 days (see Additional file 5).
Differentially methylated probes and regions
After microarray quality control and processing, 782,884 probes remained for analysis. From this filtered probe set, it was possible to create up to 109,340 background DMRs. N = 206,804 probes were ineligible to participate in regional analyses because they did not have a neighboring probe within 1 kb. Individual site analysis of the entire probe set identified 50 DMPs significantly associated with EPDS total score (empirical p = 0 after 20,000 permutations; see Additional file 5 for DMP quantile-quantile plot). Approximately 39,150 probes were taken forward for regional analysis. The median difference in adjusted R2 values for the full versus reduced models for these probes was 0.053 (interquartile range = 0.02). Ninety-eight genomic regions spanning a collective 116.2 kb were significantly differentially methylated by postpartum depressive symptom load (FDR < 5%). The significant regions overlapped 92 genes on 20 chromosomes (none on chromosomes 7 or 18), and on average spanned ~ 1.2 kb in length (see Additional file 2 for full DMR gene list). The number of CpG probes in significant DMRs ranged from 2 to 10 (mean = 3.48 probes).
Gene set enrichment and comparison of DNAm patterns associated with postpartum depressive symptoms and other genetic findings related to depression
Detailed results for the combined gene set enrichment of DMPs and DMRs can be found in Table 2. In short, the combined analysis identified one biological process (BP; cognition) and four cellular components (CC; DNA repair complex, neuron to neuron synapse, axon part, and synapse part) significant at an FDR < 5%. The DMP and DMR only analyses did not identify any categories with an FDR < 5% and produced dissimilar results (see Additional file 1). The DMP only analysis identified multiple BP and CC categories associated with neural processes (e.g., central nervous system neuron development, transmission of nerve impulses, somatodendritic compartment, and neuronal cell body), and the significance of these categories was markedly attenuated using a gene enrichment algorithm that adjusts for number of probes per gene . The DMR only analysis returned Gene Ontology (GO) categories from a variety of biological systems, including platelet formation and morphogenesis, cardiac muscle tissue development, chemical synaptic transmission, inflammatory cell apoptotic process, and tissue morphogenesis.
Three DMRs overlapped PGC GWAS findings on chromosomes 5, 6, and 16 (p = 0.034; see Additional file 3 for more detail). The overlapping region on chromosome 6 occurred in the major histocompatibility complex (MHC) region and neighbored genomic areas previously associated with early-onset major depression (DNAm)  and depression as defined in the PGC genome-wide allelic meta-analysis (see Fig. 1) . None of the sites or regions identified in the AYATS EWAS of early-onset major depression directly overlapped the DMRs associated with postpartum depressive symptom load in this analysis . The DMR with the largest AUC was located on chromosome 15 and spanned five CpG sites (Fig. 2).
Evaluating the credibility and replicability of the significant findings from this analysis is paramount but complicated. Few directly comparable studies exist, and standardized practices for conducting and reporting results from epigenome-wide association studies (EWAS) have not been established. That said, domains to assess the credibility of EWAS results have been proposed, including the level of statistical significance, genomic location, biological relevance, functional relevance, validation of significant associations, and the potential for study design bias or confounding variables to influence the analysis . An indirect effect of evaluating the results with these criteria is that it also highlights a study’s strengths and weaknesses.
This study identified 98 DMRs and 50 DMPs using DNAm measures from the EPIC beadchip (850 k), one of the most comprehensive microarray technologies available to assay DNAm. The combination strategy of site- and region-based analysis identified significant results that overlap a gene set enriched for biological pathways highly relevant to putative depression etiology (e.g., synaptic signaling and dendrite development). Moreover, the significant genome regions include a number of genes that have been previously associated with MD (e.g., RNF145 ) or with estrogen and progesterone signaling (e.g., FOXA1, ARRB2, and ITGB3BP), which may be particularly relevant for MDP and perinatal depressive symptom risk liability . None of the sites or regions from this directly overlapped regions significant in the AYATS EWAS of early-onset MD ; however, both this study and the AYATS EWAS identified significant DNAm regions in the major histocompatibility complex (MHC) region on chromosome 6, which was also identified in the latest PGC GWAS of depression . The two DNAm studies also shared a large proportion of Gene Ontology (GO) terms in their respective gene set enrichment analyses (e.g., synaptic transmission and central nervous system development), suggesting that even if the exact sites differed by study, the biological pathways with associated genes did not.
Comparing these results within the DNAm-MDP literature is difficult because a majority of the studies used candidate gene approaches, which perform best as methods to refine results from genome-wide analyses. The only other genome-wide DNAm-MDP study that used maternal tissue measured prenatal DNAm with the HumanMethylation450 beadchip (450 k). No significant results were found, possibly due to a modest sample size (n = 38 antenatal maternal blood samples)  and using only single probe approach. Modest sample sizes are common in EWAS because of technology costs, but regional analyses can mitigate power issues from small sample sizes by reducing the multiple test burden. Ong and Holbrook estimated that using their regional approach, a two-group 450 k study with 38 people (n = 19 per group) would have 61% power to identify results with an effect size of 2 using a regional approach compared to 39% with a single probe analysis . This power calculation does not map directly onto the analysis described here because Ong and Holbrook’s estimate is for a case-control study design. This study used a continuous measure of depression to avoid losing statistical power from dichotomizing a naturally quantitative trait. As a result, this study likely had at least 61% power but may not have had 80%. Another benefit of regional analyses is that the results are more likely to replicate  in part because significant regional results require multiple nearby probes each to have test statistics greater than a chosen threshold and to exhibit the same direction of effect. Together, the use of genome-wide DNAm, a permutation rank-based approach for assessing significance, and probes with test statistics either in the upper or lower 2.5 percentile to test for DMRs each positively influence the credibility of the results. Further, the genomic locations of the significant DMRs increase the credibility of the results as they overlap genes and genomic regions that either directly corroborate previous findings in the literature or that participate in biological pathways hypothesized to be important for depression risk or onset.
Another strength of this study concerns its design. Significant care was taken to minimize potential biases and confounders from biological, behavioral, and technical factors. The women in this study completed multiple comprehensive questionnaires about perinatal health and behaviors, accessed prenatal care relatively early in gestation, and were generally healthy (e.g., no diabetes) , had healthy singleton pregnancies, and completed a postpartum study visit within 7 months of delivery. These study design aspects allow variation from behavioral, biological, and technical factors to be measured and accounted for (e.g., cell type heterogeneity, tobacco use, slide and positional effects, and signal quality). Additionally, ancestrally-relevant principal components were incorporated to reduce the likelihood of detecting artifacts due to potential population stratification. Visualization of methyl values by self-reported race in significant DMRs suggested that DNAm values did not differ markedly between groups (see Additional file 4; Fig. 2).
Branching out to other depression phenotypes provides both more literature for comparison, but also more uncertainty. For example, it is not immediately clear if the DNAm patterns associated with postpartum depressive symptoms should resemble those associated with other depression phenotypes, including clinical major depression. Though both are depression phenotypes, depressive symptoms and clinical depression are not equivalent [60,61,62,63]. The issue of nonequivalence emerged when comparing the findings in this analysis to published results. For example, Numata et al. reported a significant relationship between the DNAm at cg14472315 and MD case status . No significant relationship existed between that probe and self-reported postpartum depressive symptoms in this study; however, a nominally significant relationship between that probe and lifetime history of MD was present in this study (tested post hoc), reaffirming the difficulty of comparing studies with related but non-identical depressive phenotypes.
For all of its strengths, this study was not without limitations. First, the sample size was modest, which limits the statistical power to conduct robust single-site analysis. The sample also included relatively few women with severe levels of postpartum depressive symptoms, which may have unique DNAm signatures compared to clinical MDP. Moreover, this study could not determine if the DNAm patterns identified in this study were a consequence of previous episodes of MD (which is a risk factor for perinatal depression ) or if they are related to other genetic or environmental factors. Second, detailed information about medication history was not available. Third, the newness of the EPIC microarray means that the probe set has been less well-vetted for cross-hybridization compared to its predecessor [42, 64]. Four, regional analyses are inherently limited because not all probes have the potential to form regions . This analysis attempted to address that weakness by also performing a single site analysis. Furthermore, the regional analysis algorithm used in this study selected an equal number of probes with positive and negative t values, corresponding to hypermethylation and hypomethylation in cases versus controls. This strategy assumes that equal representation of positive and negative t values will yield the most fair results; however, it is possible that this assumption limited the number of DMRs identified, especially if cases had much more hyper- or hypomethylation. Finally, this study was limited in its assessment of functional relevance. Multiple techniques were used to assess gene set enrichment, but none were specifically intended to correct for probe count per region, leading to differences in enrichment results by method. No gene expression, chromatin conformation, or transcription factor binding assays were run concurrently with postpartum DNAm analyses. No algorithms currently exist that can determine the precise change in DNAm necessary to translate into biologically meaningful differences in chromatin shape or gene regulation. As a result, fully understanding the etiology of MDP and depressive symptoms will likely involve integrating repeated measures from multiple biological layers (e.g., genetic sequence, epigenetic mechanisms, transcription, and protein) [65, 66], but no single study could measure every biological layer that might be informative about depression etiology, especially not longitudinally.
Another important consideration for interpreting DNAm results is the tissue source [29, 67]. It remains unclear how detrimental the use of peripheral blood DNAm is for identifying genomic regions associated with a psychiatric phenotype like depression. On one hand, specific brain regions hold intuitive appeal for MD-DNAm studies, and the cross-tissue similarity between brain and blood appears to be modest and tied to allelic variation ; however, given the well-established link between MD pathophysiology and aberrant immune system functioning [34,35,36, 69], peripheral blood may be the best and most feasible option for large or longitudinal studies of stress-related psychiatric traits. Many of the biomarkers associated with MD are transported in the blood (e.g., IL-1, IL-2, IL-6, TNFa, and haptoglobin), and some of these immune-related differences appear to persist after depressive episode remission. Not only does that observation fit with epidemiological studies that find individuals with a positive lifetime history of MD remain sensitive to stress and at higher risk for adverse auto-immune and cardiovascular outcomes , but also it suggests that DNAm patterns detectable in the blood may retain MD-associated differences even after depressive episode remission. Last but not least, peripheral blood may be useful for indexing changes in the relationship between the CNS and immune system. As a sentinel tissue, peripheral blood travels throughout the entire body and can deliver immune cells through the blood-brain barrier. Key players in the CNS like the neurotransmitter serotonin also serve roles in immune-related biological pathways (e.g., leukocyte activation and proliferation, cytokine secretion, chemotaxis, and apoptosis) [36, 70]. The ability of the CNS to modulate and respond to signals from the immune system underscores the intimate relationship between brain function and immune system regulation.
The final category Michels et al. (2013) listed as an important factor in establishing the credibility of DNAm associations with phenotypes is validation . Validating results typically implies either replicating the finding in an independent study (human or animal) or confirming the presence of differentially methylated probes and regions using another technology (e.g., pyrosequencing). While no pyrosequencing was completed, the literature was searched extensively for associations between depressive phenotypes and biological signatures (e.g., allelic, epigenetic, and transcription). As previously mentioned, the genomic regions implicated in this study overlap results from the 2018 Psychiatric Genomics Consortium (PGC) genome-wide allelic association study (GWAS) of MD  and are proximal to significant regions and probes from an EWAS of early-onset MD. Network analysis of the 44 significant loci in the PGC meta-analysis implicated biological pathways associated with neural differentiation, synaptic regulation, risk for schizophrenia, immune response, ion-gated channels, and retinoid X receptors . Similarly, this analysis identified significant DMRs overlapping genes coding for or related to retinoid X receptors (RXRB, ITGB3BP), genes integral to adult neurogenesis and synaptic development and positioning (e.g., AKT2, SYNGAP1, FOXG1, CTNND2, MEF2C, and AIMP1), risk for schizophrenia (e.g., CLSTN3, FARSB, MYOM2, SOX2-OT, SLC39A7, SDCCAG8, and LRRC36), immune response (LRR1, FAM19A2, CMKLR1, and MMP2), and ion-gated channels and binding (e.g., SLC6A4, SLC6A12, SLC39A7, KCNJ5, CLSTN3, and PLS3). While these similarities do not serve as direct replication, they do increase the credibility of the DNAm findings and underscore the potential for DNAm studies to complement GWA studies.
Future work should take note of the apparent differences in depressive symptoms and clinical MD and seek to refine association studies by limiting phenotypic heterogeneity. For MDP, that means not only accounting for whether the depressive symptoms or episodes onset prenatally versus postnatally, but also disentangling which DNAm patterns are associated with perinatal depressive phenotypes and which reflect pre-pregnancy events of depressive psychopathology. It is possible that women who experience their first instance of MD in the peripartum have a unique DNAm profile compared to those who have recurrent MD and happen to onset during the peripartum. Furthermore, it is unknown whether an episode of MD relatively early in life evokes a persistent perturbation in DNAm patterning that then goes on to affect the risk for additional depressive symptoms and episodes as well as other adverse health outcomes frequently comorbid with depression (e.g., cardiovascular disease and diabetes mellitus). Clarifying the phenotypes associated with DNAm patterns will enable wet lab researchers to characterize the functional and biological relevance of implicated genomic sites and regions, which is an essential step not only for understanding the biological mechanisms associated with risk and resilience to depressive psychopathology but also for developing screening tests and identifying novel pharmacotherapeutic targets.
Availability of data and materials
The PREG study research team is committed to participating in open science practices as much as possible. For this specific dataset, the IRB agreements and consent forms restrict the availability of the data due to privacy concerns. These precautions were taken as part of an effort to respect the autonomy and privacy and to encourage individuals from minority populations to feel comfortable participating in racial health disparities research. Anyone interested in access to data from participants who consented to permit data sharing is encouraged to contact Dr. Timothy P. York (corresponding author).
Summary-level data and methods can be found on the Open Science Framework (OSF; landing page https://osf.io/qsc6n). Summary statistics for individual probes and regions can be found at https://osf.io/q2s8d. Detailed information about specimen quality control assessments and covariate selection also are available on the project page in the Supplementary Information component (https://osf.io/8p9q6).
Adolescent and Young Adult Twin Study
Differentially methylated probe
Differentially methylated region
Edinburgh Postnatal Depression Scale
Epigenome-wide association study
False discovery rate
Genome-wide association study
Major depression with onset in the peripartum
Major histocompatibility complex
Principal components analysis/principal component
Psychiatric Genomics Consortium
Pregnancy, Race, Environment, Genes study
Grigoriadis S, VonderPorten EH, Mamisashvili L, et al. The impact of maternal depression during pregnancy on perinatal outcomes: a systematic review and meta-analysis. J Clin Psychiatry. 2013;74(4):e321–41. https://doi.org/10.4088/JCP.12r07968.
Liu CH, Tronick E. Rates and predictors of postpartum depression by race and ethnicity: results from the 2004 to 2007 New York city PRAMS survey (Pregnancy Risk Assessment Monitoring System). Matern Child Health J. 2013;17(9):1599–610. https://doi.org/10.1007/s10995-012-1171-z.
Nolvi S, Karlsson L, Bridgett DJ, et al. Maternal prenatal stress and infant emotional reactivity six months postpartum. J Affect Disord. 2016;199:163–70. https://doi.org/10.1016/j.jad.2016.04.020.
Field T. Prenatal depression effects on early development: a review. Infant Behav Dev. 2011;34(1):1–14. https://doi.org/10.1016/j.infbeh.2010.09.008.
Gentile S. Untreated depression during pregnancy: short- and long-term effects in offspring. A systematic review. Neuroscience. 2017;342:154–66. https://doi.org/10.1016/j.neuroscience.2015.09.001.
Coburn SS, Luecken LJ, Rystad IA, Lin B, Crnic KA, Gonzales NA. Prenatal maternal depressive symptoms predict early infant health concerns. Matern Child Health J. 2018;22(6):786–93. https://doi.org/10.1007/s10995-018-2448-7.
Meaney MJ. Perinatal maternal depressive symptoms as an issue for population health. Am J Psychiatry. 2018;175(11):1084–93. https://doi.org/10.1176/appi.ajp.2018.17091031.
American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders. Fifth. Arlington, VA, US: American Psychiatric Publishing, Inc.; 2013. doi:https://doi.org/10.1176/appi.books.9780890425596.744053
Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92. https://doi.org/10.1038/nrg3230.
Christiansen L, Lenart A, Tan Q, et al. DNA methylation age is associated with mortality in a longitudinal Danish twin study. Aging Cell. 2016;15(1):149–54. https://doi.org/10.1111/acel.12421.
Pan Y, Liu G, Zhou F, Su B, Li Y. DNA methylation profiles in cancer diagnosis and therapeutics. Clin Exp Med. 2018;18(1):1–14. https://doi.org/10.1007/s10238-017-0467-0.
Weksberg R, Shuman C, Caluseriu O, et al. Discordant KCNQ1OT1 imprinting in sets of monozygotic twins discordant for Beckwith-Wiedemann syndrome. Hum Mol Genet. 2002;11(11):1317–25.
Schenkel LC, Rodenhiser DI, Ainsworth PJ, Pare G, Sadikovic B. DNA methylation analysis in constitutional disorders: clinical implications of the epigenome. Crit Rev Clin Lab Sci. 2016;53(3):147–65. https://doi.org/10.3109/10408363.2015.1113496.
Klengel T, Mehta D, Anacker C, et al. Allele-specific FKBP5 DNA demethylation mediates gene-childhood trauma interactions. Nat Neurosci. 2013;16(1):33–41. https://doi.org/10.1038/nn.3275.
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):R115. https://doi.org/10.1186/gb-2013-14-10-r115.
Hannum G, Guinney J, Zhao L, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49(2):359–67. https://doi.org/10.1016/j.molcel.2012.10.016.
Hannon E, Dempster E, Viana J, et al. An integrated genetic-epigenetic analysis of schizophrenia: evidence for co-localization of genetic associations and differential DNA methylation. Genome Biol. 2016;17(1):176. https://doi.org/10.1186/s13059-016-1041-x.
Pries L-K, Guloksuz S, Kenis G. DNA methylation in schizophrenia. Adv Exp Med Biol. 2017;978:211–36. https://doi.org/10.1007/978-3-319-53889-1_12.
Mansell T, Novakovic B, Meyer B, et al. The effects of maternal anxiety during pregnancy on IGF2/H19 methylation in cord blood. Transl Psychiatry. 2016;6(3):e765. https://doi.org/10.1038/tp.2016.32.
Teh AL, Pan H, Chen L, et al. The effect of genotype and in utero environment on interindividual variation in neonate DNA methylomes. Genome Res. 2014;24(7):1064–74. https://doi.org/10.1101/gr.171439.113.
Viuff AC, Sharp GC, Rai D, 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. https://doi.org/10.1038/s41398-018-0286-4.
Nemoda Z, Massart R, Suderman M, et al. Maternal depression is associated with DNA methylation changes in cord blood T lymphocytes and adult hippocampi. Transl Psychiatry. 2015;5(4):e545. https://doi.org/10.1038/tp.2015.32.
Malki K, Koritskaya E, Harris F, Bryson K, Herbster M, Tosto MG. Epigenetic differences in monozygotic twins discordant for major depressive disorder. Transl Psychiatry. 2016;6(6):e839. https://doi.org/10.1038/tp.2016.101.
Byrne EM, Carrillo-Roa T, Henders AK, et al. Monozygotic twins affected with major depressive disorder have greater variance in methylation than their unaffected co-twin. Transl Psychiatry. 2013;3:e269. https://doi.org/10.1038/tp.2013.45.
Numata S, Ishii K, Tajima A, et al. Blood diagnostic biomarkers for major depressive disorder using multiplex DNA methylation profiles: discovery and validation. Epigenetics. 2015;10(2):135–41. https://doi.org/10.1080/15592294.2014.1003743.
Roberson-Nay R, Wolen AR, Lapato DM, et al. Twin study of early-onset major depression finds DNA methylation enrichment for neurodevelopmental genes. bioRxiv. 2018; http://biorxiv.org/content/early/2018/11/01/422345.
Berdasco M, Esteller M. Clinical epigenetics: seizing opportunities for translation. Nat Rev Genet. 2019;20(2):109–27. https://doi.org/10.1038/s41576-018-0074-2.
Lapato DM, Moyer S, Olivares E, et al. Prospective longitudinal study of the pregnancy DNA methylome: the US Pregnancy, Race, Environment, Genes (PREG) study. BMJ Open. 2018;8(5):e019721. https://doi.org/10.1136/bmjopen-2017-019721.
Michels KB, Binder AM, Dedeurwaerder S, et al. Recommendations for the design and analysis of epigenome-wide association studies. Nat Methods. 2013;10(10):949–55. https://doi.org/10.1038/nmeth.2632.
Ong M-L, Holbrook JD. Novel region discovery method for Infinium 450 K DNA methylation data reveals changes associated with aging in muscle and neuronal pathways. Aging Cell. 2014;13(1):142–55. https://doi.org/10.1111/acel.12159.
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.
Milgrom J, Gemmill AW. Screening for perinatal depression. Best Pract Res Clin Obstet Gynaecol. 2014;28(1):13–23. https://doi.org/10.1016/j.bpobgyn.2013.08.014.
Pidsley R, Zotenko E, Peters TJ, et al. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 2016;17(1):208. https://doi.org/10.1186/s13059-016-1066-1.
Maes M, Carvalho AF. The Compensatory Immune-Regulatory Reflex System (CIRS) in depression and bipolar disorder. Mol Neurobiol. 2018;55(12):8885–903. https://doi.org/10.1007/s12035-018-1016-x.
Medina-Rodriguez EM, Lowell JA, Worthen RJ, Syed SA, Beurel E. Involvement of innate and adaptive immune systems alterations in the pathophysiology and treatment of depression. Front Neurosci. 2018;12:547. https://doi.org/10.3389/fnins.2018.00547.
Robson MJ, Quinlan MA, Blakely RD. Immune system activation and depression: roles of serotonin in the central nervous system and periphery. ACS Chem Neurosci. 2017;8(5):932–42. https://doi.org/10.1021/acschemneuro.6b00412.
Guintivano J, Kaminsky ZA. Role of epigenetic factors in the development of mental illness throughout life. Neurosci Res. 2016;102:56–66. https://doi.org/10.1016/j.neures.2014.08.003.
Wright ML, Dozmorov MG, Wolen AR, et al. Establishing an analytic pipeline for genome-wide DNA methylation. Clin Epigenetics. 2016;8:45. https://doi.org/10.1186/s13148-016-0212-7.
R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2016.
Gentleman RC, Carey VJ, Bates DM, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80. https://doi.org/10.1186/gb-2004-5-10-r80.
Aryee MJ, Jaffe AE, Corrada-Bravo H, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30(10):1363–9 https://doi.org/10.1093/bioinformatics/btu049.
Chen Y, Lemire M, Choufani S, et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics. 2013;8(2):203–9. https://doi.org/10.4161/epi.23470.
Houseman EA, Accomando WP, Koestler DC, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13:86. https://doi.org/10.1186/1471-2105-13-86.
Gavin NI, Gaynes BN, Lohr KN, Meltzer-Brody S, Gartlehner G, Swinson T. Perinatal depression: a systematic review of prevalence and incidence. Obstet Gynecol. 2005;106(Feb):1071–83. https://doi.org/10.1097/01.AOG.0000183597.31630.db.
Gaynes BN, Gavin N, Meltzer-Brody S, et al. Perinatal depression: prevalence, screening accuracy, and screening outcomes. Evid Rep Technol Assess (Summ). 2005;(119):1–8 http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4780910/.
Altemus M, Neeb CC, Davis A, Occhiogrosso M, Nguyen T, Bleiberg KL. Phenotypic differences between pregnancy-onset and postpartum-onset major depressive disorder. J Clin Psychiatry. 2012;73(12):e1485–91. https://doi.org/10.4088/JCP.12m07693.
Du P, Zhang X, Huang C-C, et al. Comparison of beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010;11:587. https://doi.org/10.1186/1471-2105-11-587.
Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882–3. https://doi.org/10.1093/bioinformatics/bts034.
Barfield RT, Almli LM, Kilaru V, et al. Accounting for population stratification in DNA methylation studies. Genet Epidemiol. 2014;38(3):231–41. https://doi.org/10.1002/gepi.21789.
Ritchie ME, Phipson B, Wu D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
van Iterson M, van Zwet EW, Heijmans BT, Consortium the B. Controlling bias and inflation in epigenome- and transcriptome-wide association studies using the empirical null distribution. Genome Biol. 2017;18(1):19. https://doi.org/10.1186/s13059-016-1131-9.
Peters TJ, Buckley MJ, Statham AL, et al. De novo identification of differentially methylated regions in the human genome. Epigenetics Chromatin. 2015;8:6. https://doi.org/10.1186/1756-8935-8-6.
Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98(9):5116–21. https://doi.org/10.1073/pnas.091062498.
Morgan M. AnnotationHub: Client to access AnnotationHub resources. 2018.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.
Ren X, Kuan PF. methylGSA: a Bioconductor package and Shiny app for DNA methylation data length bias adjustment in gene set testing. Bioinformatics. 2019;35(11):1958–9. https://doi.org/10.1093/bioinformatics/bty892.
Wray NR, Ripke S, Mattheisen M, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genet. 2018;50(5):668–81. https://doi.org/10.1038/s41588-018-0090-3.
Cecilione JL, Rappaport LM, Hahn SE, et al. Genetic and environmental contributions of negative valence systems to internalizing pathways. Twin Res Hum Genet. 2018;21(1):12–23. https://doi.org/10.1017/thg.2017.72.
Jansen R, Penninx BWJH, Madar V, et al. Gene expression in major depressive disorder. Mol Psychiatry. 2016;21(3):339–47. https://doi.org/10.1038/mp.2015.57.
Kendler KS, Aggen SH, Li Y, et al. The similarity of the structure of DSM-IV criteria for major depression in depressed women from China, the United States and Europe. Psychol Med. 2015;45(9):1945–54. https://doi.org/10.1017/S0033291714003067.
Kendler KS, Aggen SH, Neale MC. Evidence for multiple genetic factors underlying DSM-IV criteria for major depression. JAMA psychiatry. 2013;70(6):599–607. https://doi.org/10.1001/jamapsychiatry.2013.751.
Kendler KS, Gatz M, Gardner CO, Pedersen NL. Personality and major depression: a Swedish longitudinal, population-based twin study. Arch Gen Psychiatry. 2006;63(10):1113–20. https://doi.org/10.1001/archpsyc.63.10.1113.
Kendler KS, Gardner CO, Neale MC, et al. Shared and specific genetic risk factors for lifetime major depression, depressive symptoms and neuroticism in three population-based twin samples. Psychol Med. December 2018:1–9. https://doi.org/10.1017/S003329171800377X.
Price ME, Cotton AM, Lam LL, et al. Additional annotation enhances potential for biologically-relevant analysis of the Illumina Infinium HumanMethylation450 BeadChip array. Epigenetics Chromatin. 2013;6(1):4. https://doi.org/10.1186/1756-8935-6-4.
Stricker SH, Köferle A, Beck S. From profiles to function in epigenomics. Nat Rev Genet. 2016;18(1):51–66. https://doi.org/10.1038/nrg.2016.138.
Lin E, Tsai S-J. Genome-wide microarray analysis of gene expression profiling in major depression and antidepressant therapy. Prog Neuropsychopharmacol Biol Psychiatry. 2016;64:334–40. https://doi.org/10.1016/j.pnpbp.2015.02.008.
Lappalainen T, Greally JM. Associating cellular epigenetic models with human phenotypes. Nat Rev Genet. 2017;18(7):441–51. https://doi.org/10.1038/nrg.2017.32.
Smith AK, Kilaru V, Kocak M, et al. Methylation quantitative trait loci (meQTLs) are consistently detected across ancestry, developmental stage, and tissue type. BMC Genomics. 2014;15:145. https://doi.org/10.1186/1471-2164-15-145.
Jeng J-S, Li C-T, Chen M-H, et al. Repeated low-grade infections predict antidepressant-resistant depression: a nationwide population-based cohort study. J Clin Psychiatry. 2017;79(1). https://doi.org/10.4088/JCP.17m11540.
Liu J, Buisman-Pijlman F, Hutchinson MR. Toll-like receptor 4: innate immune regulator of neuroimmune and neuroendocrine interactions in stress and major depressive disorder. Front Neurosci. 2014;8:309. https://doi.org/10.3389/fnins.2014.00309.
We are extremely grateful to the families, clinicians, and research staff who participated in the Pregnancy, Race, Environment, Genes (PREG) study and its postpartum extension. Additionally, we would like to thank the Virginia Commonwealth University Epigenetics Journal Club members who provided feedback on the manuscript.
The PREG study and its postpartum extension were funded by the National Institute on Minority Health and Health Disparities (P60MD002256, PI: TPY, Strauss), The John and Polly Sparks Foundation and Brain and Behavior Research Foundation (24712, PI: TPY), American Nurses Foundation Research Grant (5232, PI: PAK), and VCU Center for Clinical and Translational Research Endowment Fund (6-40595, PI: PAK). DL is supported by National Institute of Mental Health T32MH020030. The use of REDCap was supported by Clinical and Translational Science Award (CTSA) award No. UL1TR000058 from the National Center for Advancing Translational Sciences.
Ethics approval and consent to participate
Both the PREG study and its postpartum extension were approved by the Virginia Commonwealth University Institutional Review Board (#14000), and participants provided written consent for both.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
DMP-only and DMR-only gene enrichment analysis. This file contains gene ontology enrichment results in tabular format for the combined set of differentially methylated regions and probes (DMR, DMP) results, enrichment analysis of DMRs only, and enrichment analysis of DMPs only. CSV files for these results are on the Open Science Framework (OSF; landing page https://osf.io/qsc6n).
List of genes overlapping differentially methylated regions. List of the 92 genes overlapped by significant differentially methylated regions (DMRs).
Enrichment testing of Psychiatric Genomic Consortium (PGC) supplemental methods. This file contains additional details about how the 95% confidence intervals were calculated for the PGC enrichment analysis.
Supplement to Figure 2. This file contains two additional versions of the differentially methylated region (DMR) highlighted in Fig. 2. The top figure shows the ComBat-adjusted methyl values for each participant colored by self-identified Census-based race category. The points have been jittered left/right to ease visualization by reducing over-plotting. No vertical adjustment was made. The bottom figure shows the mean methyl values for each probe contained in the DMR by self-identified Census-based race category.
Additional Figures. This file contains four figures: quantile-quantile plots for the DMP and DMR analyses and histograms of the distribution of when PREG postpartum study visits occurred (i.e., time since birth) and the distribution of EPDS total scores by self-reported lifetime history of MD (assessed using an extended self-report version of the CIDI-SF).
About this article
Cite this article
Lapato, D.M., Roberson-Nay, R., Kirkpatrick, R.M. et al. DNA methylation associated with postpartum depressive symptoms overlaps findings from a genome-wide association meta-analysis of depression. Clin Epigenet 11, 169 (2019). https://doi.org/10.1186/s13148-019-0769-z