DNA methylation variation at HP1BP3 and TTC9B is modified by estrogen exposure in the rodent hippocampus and was previously shown to be prospectively predictive of postpartum depression (PPD) when modeled in antenatal blood. The objective of this study was to replicate the predictive efficacy of the previously established model in women with and without a previous psychiatric diagnosis and to understand the effects of changing hormone levels on PPD biomarker loci. Using a statistical model trained on DNA methylation data from N = 51 high-risk women, we prospectively predicted PPD status in an independent N = 51 women using first trimester antenatal gene expression levels of HP1BP3 and TTC9B, with an area under the receiver operator characteristic curve (AUC) of 0.81 (95% CI: 0.69-0.92, po5 × 10 − 4 ). Modeling DNA methylation of these genes in N = 240 women without a previous psychiatric diagnosis resulted in a cross-sectional prediction of PPD status with an AUC of 0.81 (95% CI: 0.68-0.93, p = 0.01). TTC9B and HP1BP3 DNA methylation at early antenatal time points showed moderate evidence for association to the change in estradiol and allopregnanolone over the course of pregnancy, suggesting that epigenetic variation at these loci may be important for mediating hormonal sensitivity. In addition both loci showed PPD-specific trajectories with age, possibly mediated by age-associated hormonal changes. The data add to the growing body of evidence suggesting that PPD is mediated by differential gene expression and epigenetic sensitivity to pregnancy hormones and that modeling proxies of this sensitivity enable accurate prediction of PPD.

The oxytocin receptor (OXTR) is a key regulator of stress and anxiety and may be regulated by both psychosocial risk factors and gonadal hormones, making it an attractive candidate for study in postpartum depression (PPD). The objective of this study was to investigate both serum hormone and PPD specific DNA methylation variation in the OXTR. Illumina HM450 microarray data generated in a prospective PPD cohort identified significant associations (P=0.014) with PPD in an intronic region in the OXTR located 4bp proximal to an estrogen receptor (ER) binding region. Pyrosequencing confirmed moderate evidence for an interaction of CpGs in the region with childhood abuse status to mediate PPD. These CpGs located on chr3 at positions 8810078 and 8810069 exhibited significant associations with postpartum depression scores from an independent cohort of 240 women with no prior psychiatric history. Hormone analysis suggested a PPD specific negative correlation of DNA methylation in the region with serum estradiol levels. Estradiol levels and OXTR DNA methylation exhibited a significant interaction to associate with the ratio of allopregnanolone to progesterone. Cumulatively, the data corroborate our previous hypotheses of a PPD specific increased sensitivity of epigenetic reprogramming at estrogen target genes and suggests that OXTR epigenetic variation may be an important mediator of mood relevant neuroactive steroid production.

BackgroundSuicide is the second leading cause of death among adolescents in the USA, and rates are rising. Methods to identify individuals at risk are essential for implementing prevention strategies, and the development of a biomarker can potentially improve prediction of suicidal behaviors. Prediction of our previously reported SKA2 biomarker for suicide and PTSD is substantially improved by questionnaires assessing perceived stress or anxiety and is therefore reliant on psychological assessment. However, such stress-related states may also leave a biosignature that could equally improve suicide prediction. In genome-wide DNA methylation data, we observed significant overlap between waking cortisol-associated and suicide-associated DNA methylation in blood and the brain, respectively.ResultsUsing a custom bioinformatic brain to blood discovery algorithm, we derived a DNA methylation biosignature that interacts with SKA2 methylation to improve the prediction of suicidal ideation in our existing suicide prediction model across both blood and saliva data sets. This biosignature was independently validated in the Grady Trauma Project cohort and interacted with HPA axis metrics in the same cohort. The biosignature showed a relationship with immune status by its correlation with myeloid-derived cell proportions in all data sets and with IL-6 measures in a prospective postpartum depression cohort. Three probes showed significant correlations with the biosignature: cg08469255 (DDR1), cg22029879 (ARHGEF10), and cg24437859 (SHP1), of which SHP1 methylation correlated with immune measures.ConclusionsWe conclude that this biosignature interacts with SKA2 methylation to improve suicide prediction and may represent a biological state of immune and HPA axis modulation that mediates suicidal behavior.Electronic supplementary materialThe online version of this article (doi:10.1186/s13148-016-0279-1) contains supplementary material, which is available to authorized users.