- Open Access
Prenatal adverse environment is associated with epigenetic age deceleration at birth and hypomethylation at the hypoxia-responsive EP300 gene
Clinical Epigeneticsvolume 11, Article number: 73 (2019)
Obstetric complications have long been retrospectively associated with a wide range of short- and long-term health consequences, including neurodevelopmental alterations such as those observed in schizophrenia and other psychiatric disorders. However, prospective studies assessing fetal well-being during pregnancy tend to focus on perinatal complications as the final outcome of interest, while there is a scarcity of postnatal follow-up studies. In this study, the cerebroplacental ratio (CPR), a hemodynamic parameter reflecting fetal adaptation to hypoxic conditions, was analyzed in a sample of monozygotic monochorionic twins (60 subjects), part of them with prenatal complications, with regard to (i) epigenetic age acceleration, and (ii) DNA methylation at genes included in the polygenic risk score (PRS) for schizophrenia, and highly expressed in placental tissue.
Decreased CPR measured during the third trimester was associated with epigenetic age deceleration (β = 0.21, t = 3.362, p = 0.002). Exploration of DNA methylation at placentally expressed genes of the PRS for schizophrenia revealed methylation at cg06793497 (EP300 gene) to be associated with CPR (β = 0.021, t = 4.385; p = 0.00008, FDR-adjusted p = 0.11). This association was reinforced by means of an intrapair analysis in monozygotic twins discordant for prenatal suffering (β = 0.027, t = 3.924, p = 0.001).
Prenatal adverse environment during the third trimester of pregnancy is associated with both (i) developmental immaturity in terms of epigenetic age, and (ii) decreased CpG-specific methylation in a gene involved in hypoxia response and schizophrenia genetic liability.
Prenatal environment constitutes the first modulating agent the developing fetus encounters as it progresses through gestation. The tremendous impact of any environmental threat occurring during this period for both short- and long-term consequences is now widely accepted and well-known as the Developmental Origins of Health and Disease (DOHaD) hypothesis . Also known as the theory of fetal programming, the embedding of early life and its ability to exert long-term effects in late-life is thought to rely on epigenetic mechanisms [2, 3].
Recently, several DNA methylation-based epigenetic clocks have been developed in order to predict chronological age with high accuracy [4, 5]; afterward, Knight and colleagues developed a new predictor specifically aimed to predict gestational age (GA) in perinatal samples . Although epigenetic and chronological age robustly show high correlation across studies, the difference between both variables allows the estimation of the so-called age acceleration (i.e., when epigenetic age is higher than chronological age).
On the one hand, epigenetic age acceleration in adult subjects has been associated with cumulative lifetime stress, lifestyle, and all-cause mortality, among others, suggesting its utility as a better predictor for life expectancy than chronological age itself [7,8,9]. On the other hand, epigenetic GA deceleration (i.e., when chronological age is higher than epigenetic age), as measured in cord blood, has been described in newborns born to women with low socioeconomic status, Sjögren syndrome, insulin-treated gestational diabetes mellitus, and experiencing antenatal depressive symptoms [6, 10, 11]. Such findings suggest that newborns exposed to prenatal stressors are born in an immature state independently of their chronological GA. In this regard, boys—but not girls—who exhibited lower epigenetic GA at birth exhibited more internalizing problems, such as anxious-depressive symptoms or somatic complaints, at follow-up (mean age 3.7 years), suggesting they are born with a developmental disadvantage .
Nevertheless, there is a dearth of studies examining the putative relationship between ultrasound parameters acquired during pregnancy and epigenetic GA acceleration. In this regard, the cerebroplacental ratio (CPR) has been reported to be associated with adverse perinatal outcomes not only in growth-restricted fetuses, but also in low-risk population [12, 13]. Briefly, CPR is calculated by dividing the middle cerebral artery (MCA) pulsatility index (PI) by the umbilical artery (UA) PI . The PI is a parameter reflecting vascular impedance or resistance, i.e., decreased blood flow. Specifically, fetal brain blood supply is known to increase in front of hypoxic stimuli thus decreasing PI in the MCA ; while placental insufficiency decreases umbilical blood flow hence increasing UA PI, and has been associated with both short- and long-term detrimental outcomes, including increased cardiovascular risk and deficits in cognition [16, 17]. Consequently, a decreased CPR reflects the combination of both alterations and is an indicator of fetal adaptation to adverse conditions .
Obstetric complications (OCs) constitute one of the risk factors more reliably associated with psychopathology, particularly with neurodevelopmental disorders; specifically, the putative association between OCs and schizophrenia has been debated since the 1970s [18,19,20]. In this regard, a recent umbrella review evaluating all published meta-analysis regarding risk factors and biomarkers for schizophrenia spectrum disorders revealed a history of OCs to significantly increase the risk for developing the disorder with an odds ratio of 2 . Furthermore, exposure to severe OCs together with increased genetic vulnerability, as measured with the polygenic risk score (PRS) for schizophrenia, interact to increase the risk to suffer the disorder up to an odds ratio of 8.36 . In the same study, authors further explored the putative relevance of placental expression of genes included in the PRS; following this approach, they reported (i) an enrichment of PRS genes expressed in placental tissue and (ii) differential expression of PRS genes in placentae from complicated pregnancies (specifically in pre-eclampsia and intrauterine growth restriction). Specifically, the described gene-environment interaction between exposure to OCs and the PRS for schizophrenia was driven by those genes highly expressed in placenta and/or dynamically regulated in complicated pregnancies . Since CPR is a robust indicator of prenatal stress and a predictor of perinatal and long-term morbidity, DNA methylation analysis of genes included in the placental PRS for schizophrenia could shed light on the epigenetic mechanisms mediating the interaction between OCs and neurodevelopmental disorders.
Monozygotic twins have been instrumental for the elucidation of environmental and genetic risks in the etiology of complex traits and disorders. Actually, the differential role of the prenatal environment in shaping psychopathological proneness was first described thanks to monozygotic twin designs [23,24,25]; these pioneering studies focused on dermatoglyphic measures assessed at birth, which can be used as surrogate measures of altered neurodevelopment during the second trimester of pregnancy . Furthermore, monozygotic twin pregnancies and, more specifically, monochorionic twin pregnancies—i.e., those in which both fetuses share the placenta—are at a higher risk of obstetric complications, the more prevalent being twin-to-twin transfusion syndrome (TTTS) and selective intrauterine growth restriction (sIUGR) [27,28,29]. Thus, the thorough and prospective ultrasound assessment of prenatal development through monochorionic twin pregnancies offers a quasi-experimental study design in which the genetic and environmental components of epigenetic variability can be dissected.
The objective of the current study was to investigate whether prenatal adverse environment (i) alters human development in terms of epigenetic age, and if (ii) it can get embedded through epigenetic mechanisms in genes previously identified as risk factors for schizophrenia acting during prenatal stages. We hypothesized that a higher exposure to prenatal adverse environment would be associated with (i) delayed development and (ii) DNA methylation at genes involved in the pathogenesis of schizophrenia. While CPR can have diverse effects on genome-wide DNA methylation, with potential relevance for a multitude of phenotypes, the present study a priori examined how CPR epigenetically regulates risk genes for schizophrenia, previously described to interact with the presence of OCs .
GA estimation using Knight’s epigenetic clock
After exclusion of two twin pairs (see Methods section), the final sample size was 30 twin pairs. The mean GA at birth of our twin cohort (n = 30 twin pairs) was 35.3 weeks (range = 31.7–37.1) and the mean DNA methylation GA at birth was 35 weeks (range = 31.4–37.7). To validate the epigenetic clock predictor in our sample, DNA methylation-based GA was tested for correlation with chronological GA (r = 0.76, p = 1.68 × 10−12; Fig. 1). The average absolute difference between epigenetic GA and chronological GA—hereinafter referred as ΔGA—was 0.9 weeks (range = 0.03–4.02), i.e., 6.3 days. In agreement with previous studies, there was a significant negative correlation between ΔGA and chronological GA (r = − 0.47; p < 0.001).
Association between ΔGA and CPR
ΔGA was tested for associations with CPR measured during the third trimester (mean = 33.8 weeks, range = 28.3–36.4), a few days before childbirth (median = 6.5 days). CPR was significantly associated with ΔGA (β = 0.21, t = 3.362, p = 0.002) when adjusting for sex, birthweight, diagnostic of either TTTS or sIUGR, surgery time interval (when laser fetoscopy had been applied), and gestational age at ultrasound as covariates. The positive association between CPR and ΔGA remained significant after correction for cell type proportion (β = 0.21, t = 2.616, p = 0.01). Figure 2 shows the positive association between third trimester CPR and ΔGA.
Epigenetic exploration of placental PRS for schizophrenia with regard to CPR
Following the approach developed by Ursini and collaborators (2018), association between CPR and DNA methylation was tested in all CpG sites included in the DNA methylation array located within genes of the PRS for schizophrenia expressed in placental tissue (placental PRS) . There were 1400 CpG sites annotated to placental PRS genes out of 866,091 CpG sites included in the array. After FDR correction for multiple testing, methylation at one single CpG site, cg06793497, was significantly associated with CPR (β = 0.021, p = 0.00008, t = 4.385; qFDR adjusted = 0.11; Fig. 3a), such that increased cg06793497 methylation was associated with increased CPR. The top 10 CpG sites yielded by this approach are summarized in Table 1 (all q values > 0.75).
To further explore the association between cg06793497 methylation and CPR, it was analyzed in a monozygotic twin intrapair design. The intrapair twin design further allows controlling for chronological GA, sex, and timing of the Doppler ultrasound, since these variables are shared by co-twins of a pair. Four observations were removed from the analysis due to missingness for any of the variables in one of the co-twins of a pair. Thus, intrapair differences for these measures were calculated for all twin pairs of the sample were both measures were available for both twins of a pair (n = 27 twin pairs). Intrapair differences in cg06793497 methylation and CPR, measured during the third trimester, were significantly correlated (r = 0.64, p < 0.001; Fig. 3b). The association between both variables remained significant after adjusting for cell type count intrapair differences (β = 0.027, t = 3.924, p = 0.001). Intrapair exploration of the top 10 CpG sites (Table 1) revealed significant associations between CPR and DNA methylation at CpG probes cg00262246 (β = 0.012, p = 0.029), cg01024069 (β = − 0.01, p = 0.033), and cg12955069 (β = − 0.021, p = 0.026).
DNA methylation exploration of EP300 gene
To further explore the putative relevance of DNA methylation at other CpG sites located within the EP300 gene and its surrounding regions, DNA methylation at 27 CpG sites included in the array and annotated to this region was also explored with regard to CPR (see Table 2). All analyses were adjusted for cell sex, birthweight, gestational age at ultrasound, and cell type count. In addition to cg06793497, two additional CpG sites—cg12968540 and cg19011939—were significantly associated with CPR (p < 0.05); moreover, methylation at four additional CpG sites—cg04452260, cg24349919, cg11931284 and cg25888227—showed trend associations with CPR (p < 0.01). The intrapair approach was then applied for these newly identified six CpG sites revealing cg11931284 (β = 0.028, t = 2.985, p = 0.008) and cg19011939 (β = − 0.021, t = − 2.343, p = 0.03) to be significantly associated with CPR, when adjusting for cell types intrapair differences.
To the best of our knowledge, this is the first study analyzing the epigenetic age in association with adverse prenatal environment as measured by a hemodynamic ultrasound parameter. Firstly, we describe the significant association between CPR measured during the third trimester of pregnancy with epigenetic age acceleration. Specifically, subjects exhibiting decreased CPR—exposed to prenatal adverse conditions—were born with decelerated epigenetic age, i.e., prenatally stressed subjects were born immature adjusting for their gestational age at birth. Additionally, methylomic exploration of schizophrenia PRS genes known to be expressed in placenta revealed the association between CPR and EP300 gene CpG-specific methylation, at the cg06793497 probe, in our monochorionic twin sample.
Developmental deficits and developmental delays have been previously described in children who would later develop schizophrenia ; although such prodromal symptoms were in accordance with the neurodevelopmental hypothesis for schizophrenia, biological mechanisms mediating these effects remain largely unknown. Epigenetic immaturity in response to prenatal stress could be contributing to this developmental delay. Interestingly, epigenetic age deceleration has been previously described in association with maternal pathologies during pregnancy, such as maternal depression or Sjögren’s syndrome, suggesting it can be a robust biomarker of prenatal suffering [10, 11]. It is worth noting that CPR was measured a few days prior to childbirth; thus, it can be used as a surrogate marker of prenatal adaptation to adverse conditions experienced at the end of the pregnancy, i.e., as a marker of perinatal risk.
Integration of the schizophrenia PRS  with obstetric and placental information , allowed the identification of E1A binding protein p300 (EP300) gene CpG-specific methylation as a putative marker of exposure to prenatal stress. Interestingly, the EP300 gene encodes a histone acetyltransferase (HAT) involved in several cell pathways such as cell proliferation and differentiation. Mutations at EP300 gene have been described to cause Rubinstein-Taybi syndrome, a rare autosomal dominant neurodevelopmental disorder characterized by intellectual disability, psychomotor and language delay, and facial dysmorphisms . Likewise, these symptoms, including developmental delay, learning problems, and cleft palate, characterize the 22q11.2 deletion syndrome, a well-defined congenital condition caused by the deletion of the 22q11.2 segment . Notably, this syndrome is associated with a higher risk to develop schizophrenia, among other psychiatric conditions ; interestingly, EP300 gene is located on chromosome 22 at position 22q13.2.
Further exploration of differential DNA methylation in and around the EP300 gene revealed cg19011939 to be differentially methylated in association with prenatal adversity. While higher exposure to a prenatal adverse environment, as reflected by lower CPR during the third trimester, is associated with decreased methylation at cg06793497 in the hypoxia-responsive EP300 gene, there appears to be increased methylation at cg19011039 at EP300-AS1 gene. Thus, we speculate that higher exposure to prenatal stress might be associated with reciprocal patterns of EP300 and EP300-AS1 epigenetic regulation that could act synergistically, a hypothesis that may be explored in future studies .
Remarkably, EP300 has been identified as a co-activator of the hypoxia-inducible factor 1 alpha (HIF1A). In this regard, hypoxic conditions stimulate EP300 expression, which has a neuroprotective role . Accordingly, genetic variability at EP300 gene has been associated with human adaptations to high altitude regions, e.g., the Tibet . Likewise, pre- and peri-natal hypoxia have been associated with schizophrenia spectrum disorders, particularly by decreasing hippocampal volume [38, 39]; complementarily, a decreased or impaired response to hypoxia via neurotrophic factors has also been implicated in the etiology of schizophrenia . Furthermore, DNA methylation at the IGF2BP1 gene, also involved in prenatal development , has been associated with both adult working memory and birthweight ; further highlighting the advantage of twin study designs to identify environmentally-driven epigenetic consequences of prenatal stress. Overall, these findings point to the existence of a GxE interaction between genetic vulnerability and exposure to prenatal hypoxia, as already highlighted by Ursini and collaborators . In this framework, EP300 methylation could be one of the mediators of such interaction.
A number of limitations of the present study should be noted. First, the moderate sample size (n = 60 subjects, 30 twin pairs) limits the statistical power of the analysis; however, smaller sample sizes (n = 22 MZ twin pairs) have been described to be sufficient to identify methylation differences of 6% with > 80% power . Moreover, a lenient significance threshold after correction multiple testing was used; however, previous epigenetic studies have described FDR values between 5 and 20% as markers of medium-confidence sites . Another limitation regards the moderately small reported effect sizes (around 2%) questioning the biological relevance of our findings ; however, these findings are in agreement with a larger body of evidence regarding cord blood methylation after exposure to a number of prenatal stressors. Such small DNA methylation changes may act in conjunction with a myriad of other epigenetic signatures and biological processes in order to maintain homeostasis in the face of threats. Additionally, although epigenomic information was available from a methylomic array including more than 800,000 CpG sites distributed throughout the whole human genome, only 1400 CpG sites were analyzed; alternative approaches including the total of CpG sites included in the array would have yield different findings, probably pointing to genes involved in other neurodevelopmental disorders besides schizophrenia. Furthermore, while the set of genes analyzed in the current approach were described to be highly expressed in placental tissue , placentae were not available for this sample and cord blood was thus analyzed as the proxy tissue of choice with regard to exposure to prenatal adversity. Finally, MZ twin pregnancies are characterized by lower gestational ages at birth than singleton pregnancies; besides, obstetric scales commonly used in psychiatric studies include twin pregnancies as an obstetric complication. Thus, findings derived from the present design might not be generalizable to the general population.
Further studies are needed to test the time stability of the hereby identified methylation signature. It will be equally relevant to explore neurobehavioral correlates of EP300 methylation during early childhood along with its putative association with neurodevelopmental outcomes, including psychosis liability. Additionally, a longitudinal follow-up is required to test the role of postnatal environment in these phenotypes since both epigenetic age deceleration and CpG-specific differential methylation in association with CPR could return to basal levels after birth. Finally, genetic exploration of these subjects regarding schizophrenia PRS will be instrumental for the study of GxE interactions and genetic liability for an impaired hypoxia response during human development.
This was a prospective study including fetal pairs from monochorionic diamniotic twin pregnancies attended at Hospital Clínic de Barcelona (Spain) during a 2-year recruitment period. Monochorionic monoamniotic twin pregnancies were excluded from the present study to avoid putative confounding with regard to differential exposure to stress in both types of twin pregnancies. The study protocol was approved by the hospital ethics committee (HCB/2016/0046), and all patients provided written informed consent.
We included 32 monochorionic pregnancies (n = 64 samples). The sample was enriched for two monochorionic-specific severe obstetric complications: twin-to-twin transfusion syndrome (TTTS, n = 8) and selective intrauterine growth restriction (sIUGR, n = 9). All TTTS cases were treated upon detection by means of laser fetoscopy .
Maternal age and pre-pregnancy BMI were retrieved from hospital records. Gestational age was dated using first-trimester crown-rump length measurement of the larger twin .
Fetal ultrasound assessment
Ultrasound assessment was performed on a Voluson Expert 8 (General Electrical Medical Systems, Milwaukee, WI, USA) or a Siemens Sonoline Antares (Siemens Medical Systems, Erlangen, Germany) with 8- to 4-MHz or 6- to 4- MHz curved array probes, respectively. All fetuses underwent detailed ultrasound evaluation including fetal anatomy and Doppler measurements such as UAPI, MCAPI and ductus venosus PI. All Doppler evaluations were acquired at a normal fetal heart rate (FHR) in the absence of fetal body or respiratory movements and at an angle of insonation as close to 0° as possible (but always < 15°), and the mechanical and thermal indices were maintained below 1. CPR was calculated as the ratio between MCAPI and UAPI, according to previous studies .
Umbilical vein cord blood samples were obtained from the clamped umbilical cord immediately after delivery of the fetus. All blood samples were collected in EDTA-treated tubes and processed within 1 h. Plasma was separated by centrifugation at 3000 rpm for 10 min at 4 °C, and stored at − 80 °C until further use. Genomic DNA was extracted from fetal cord blood using QIAamp DNA Mini Kit (Qiagen). DNA quality and quantity were assessed by NanoDrop One (Thermo Scientific). Genomic DNA was bisulfite converted using the Zymo EZ-96 DNA Methylation Kit (Zymo Research). Genome-wide DNA methylation levels were assessed over 850,000 CpG sites by means of the Infinium MethylationEPIC BeadChip Kit (Illumina Inc., CA, USA) according to the manufacturer’s protocol. Pre-processing and normalization were performed using the Bioconductor minfi package . CpG probes containing common SNPs were discarded. All probes mapping to the X and Y chromosomes were also removed. Finally, cross-hybridizing probes as previously identified were excluded from further analysis . All samples (n = 64) were run on the same plate.
Absence of maternal contamination was confirmed after retrieving DNA methylation values at 10 CpG sites previously described to identify sample contamination by maternal blood during sample collection . None of the samples assayed exhibited DNA methylation values above the threshold at 5 or more of those CpG sites (see Additional file 1 for specific methylation values). Two samples (from the same twin-pair) were excluded from further analyses due to lack of monozygosity as assessed by 59 SNPs included in the array. One of the samples was removed from analysis due to insufficient DNA concentration, the co-twin sample was also excluded from further analysis.
All statistical analyses were conducted in R version 3.5.0 . DNA methylation-based GA prediction was performed using the R code and statistical pipeline developed by Knight, based on the methylation profile of 148 CpG sites ; this predictor was developed using 15 Illumina DNA methylation datasets (n = 1434 neonates). Following Simpkin et al. recommendations, the Knight clock was preferred for our analysis as it was developed and tested in preterm infants datasets such as our monozygotic twin population, characterized by a mean gestational age at birth of 35.3 weeks . The EPIC array lacks 6 of the CpG sites originally included in the Knight clock, these values were imputed manually as non-available. Interestingly, DNA methylation-based age estimation relying on EPIC array data has already been described to accurately predict age despite the lack of several CpG sites originally included in Horvath’s and Hannum’s clocks .
Gestational age acceleration (ΔGA) was calculated as the absolute difference between epigenetic GA and chronological GA. Since ΔGA was associated with chronological GA (r = − 0.47; p < 0.001), the latter was included as a covariate in all statistical models; this association has been already reported in prior studies exploring epigenetic-based GA estimations at birth [6, 10, 11].
Cell counts of CD4+ T cells, CD8+ T cells, B cells, NK cells, granulocytes, monocytes, and nucleated red blood cells (nRBCs) were estimated using the R code and statistical pipeline developed by Houseman .
A multiple linear regression model was built to analyze the correlation between ΔGA and CPR. Fetal sex, birthweight, diagnostic of either TTTS or sIUGR (binary variable), post-surgery interval (in TTTS cases where laser fetoscopy had been applied), and gestational age at ultrasound were included as independent variables in the model as they are known to influence either DNA methylation (from which ΔGA is calculated) or CPR. This analysis was conducted in the total MZ twin sample (n = 60).
DNA methylation at CpG sites annotated to the 43 genes of the Placental PRS1 as described by Ursini et al.  was retrieved to test their association with CPR. A second multiple linear regression model was then designed to explore putative effects of CPR upon methylation of PRS genes, testing 1,400 associations. The aforementioned confounding variables along with cell types proportions (CD4+ T cells, CD8+ T cells, B cells, NK cells, granulocytes, monocytes, and nRBCs) were included as covariates, as they are known to affect methylation values. False discovery rate (FDR) correction for multiple testing was applied, considering q values under 20% to be indicative of medium-confidence probes following prior studies .
A twin-based approach previously developed in our group  was also applied to refine the association between cg06793497 methylation and CPR. Briefly, intrapair differences for both variables of interest were computed for each twin pair; afterward, a regression model was fitted with an estimated intrapair cg06793497 methylation (Δmethylation) and intrapair CPR (ΔCPR). This last model was not adjusted for either sex or chronological gestational age since both variables are identical for both twins of a pair.
Graignic-Philippe R, Dayan J, Chokron S, Jacquet AY, Tordjman S. Effects of prenatal stress on fetal and child development: A critical literature review. Neurosci Biobehav Rev. 2014;43:137–62 Elsevier Ltd.
Tobi EW, Goeman JJ, Monajemi R, Gu H, Putter H, Zhang Y, et al. DNA methylation signatures link prenatal famine exposure to growth and metabolism. Nat Commun. 2014;5:5592.
Palma-Gudiel H, Cirera F, Crispi F, Eixarch E, Fañanás L. The impact of prenatal insults on the human placental epigenome: A systematic review. Neurotoxicol Teratol Pergamon. 2018;66:80–93.
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14:R115.
Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda SV, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49:359–67 Elsevier Inc.
Knight AK, Craig JM, Theda C, Bækvad-Hansen M, Bybjerg-Grauholm J, Hansen CS, et al. An epigenetic clock for gestational age at birth based on blood methylation data. Genome Biol. 2016;17:1–11.
Zannas AS, Arloth J, Carrillo-Roa T, Iurato S, Röh S, Ressler KJ, et al. Lifetime stress accelerates epigenetic aging in an urban, African American cohort: Relevance of glucocorticoid signaling. Genome Biol. 2015;16:1–12.
Chen B, Marioni RE, Bressler J, Fornage M, Studenski S, Vandiver AR, et al. DNA methylation-based measures of biological age : meta-analysis predicting time to death. Aging (Albany NY). 2016;8:1–22.
Declerck K, Vanden Berghe W. Back to the future: epigenetic clock plasticity towards healthy aging. Mech Ageing Dev. 2018;174:18–29 Elsevier.
Girchenko P, Lahti J, Czamara D, Knight AK, Jones MJ, Suarez A, et al. Associations between maternal risk factors of adverse pregnancy and birth outcomes and the offspring epigenetic clock of gestational age at birth. Clinical Epigenetics. 2017;9:1–14.
Suarez A, Lahti J, Czamara D, Lahti-Pulkkinen M, Knight AK, Girchenko P, et al. The epigenetic clock at birth: associations with maternal antenatal depression and child psychiatric problems. J Am Acad Child Adolesc Psychiatry. 2018;57:321–328.e2.
DeVore GR. The importance of the cerebroplacental ratio in the evaluation of fetal well-being in SGA and AGA fetuses. Am J Obstet Gynecol. 2015;213:5–15 Elsevier Inc.
Morales-Roselló J, Khalil A, Morlando M, Papageorghiou A, Bhide A, Thilaganathan B. Changes in fetal Doppler indices as a marker of failure to reach growth potential at term. Ultrasound Obstet Gynecol. 2014;43:303–10.
Bahado-Singh RO, Kovanci E, Jeffres A, Oz U, Deren O, Copel J, et al. The Doppler cerebroplacental ratio and perinatal outcome in intrauterine growth restriction. Am J Obstet Gynecol. 1999;180:750–6.
Peeters LLH, Sheldon RE, Jones MD Jr, Makowski EL, Meschia G. Blood flow to fetal organs as a function of arterial oxygen content. Am J Obstet Gynecol. 1979;135:637–46 Elsevier.
Mone F, Thompson A, Stewart MC, Ong S, Shields MD. Fetal umbilical artery Doppler pulsatility index as a predictor of cardiovascular risk factors in children—a long-term follow up study. J Matern Neonatal Med. 2014;27:1633–6 Taylor & Francis.
Mone F, McConnell B, Thompson A, Segurado R, Hepper P, Stewart MC, et al. Fetal umbilical artery Doppler pulsatility index and childhood neurocognitive outcome at 12 years. BMJ Open. 2016;6:e008916.
Gottesman II, Shields J. A critical review of recent adoption, twin, and family studies of schizophrenia: behavioral genetics perspectives. Schizophr Bull. 1976;2:360–401.
Cannon M, Ph D, Jones PB, Ph D, Murray RM, Sc D, et al. Obstetric complications and schizophrenia : historical and meta-analytic review. Am J Psychiatry. 2002;159:1080–92.
Bale TL, Baram TZ, Brown AS, Goldstein JM, Insel TR, McCarthy MM, et al. Early life programming and neurodevelopmental disorders. Biol Psychiatry. 2010;68:314–9 Elsevier Inc.
Belbasis L, Köhler CA, Stefanis N, Stubbs B, Os J, Vieta E, et al. Risk factors and peripheral biomarkers for schizophrenia spectrum disorders: an umbrella review of meta-analyses. Acta Psychiatr Scand. John Wiley & Sons, Ltd (10.1111). 2017;137:88–97.
Ursini G, Punzi G, Chen Q, Marenco S, Robinson JF, Porcelli A, et al. Convergence of placenta biology and genetic risk for schizophrenia. Nat Med. 2018;24:1–10.
Bracha HS, Torrey EF, Gottesman II, Bigelow LB, Cunniff C. Second-trimester markers of fetal size in schizophrenia: a study of monozygotic twins. Am J Psychiatry. 1992;149:1355–61 American Psychiatric Publishing.
Torrey EF, Taylor EH, Bracha HS, Bowler AE, McNeil TF, Rawlings RR, et al. Prenatal origin of schizophrenia in a subgroup of discordant monozygotic twins. Schizophr. Bull. US: National Institute of Mental Health; 1994. p. 423–432.
Rosa A, Fañanas L, Bracha HS, Torrey EF, van Os J. Congenital dermatoglyphic malformations and psychosis: a twin study. Am J Psychiatry. American Psychiatric Publishing. 2000;157:1511–3.
Fañanas L, van Os J, Hoyos C, McGrath J, Mellor CS, Murray R. Dermatoglyphic a-b ridge count as a possible marker for developmental disturbance in schizophrenia: replication in two samples. Schizophr Res. 1996;20:307–14.
Valsky DV, Eixarch E, Martinez JM, Crispi F, Gratacós E. Selective intrauterine growth restriction in monochorionic twins: pathophysiology, diagnostic approach and management dilemmas. Semin Fetal Neonatal Med. Elsevier Ltd. 2010;15:342–8.
Sebire NJ, Snijders RJM, Hughes K, Sepulveda W, Nicolaides KH. The hidden mortality of monochorionic twin pregnancies. BJOG An Int J Obstet Gynaecol. John Wiley & Sons, Ltd (10.1111). 1997;104:1203–7.
Chalouhi GE, Stirnemann JJ, Salomon LJ, Essaoui M, Quibel T, Ville Y. Specific complications of monochorionic twin pregnancies: twin-twin transfusion syndrome and twin reversed arterial perfusion sequence. Semin Fetal Neonatal Med. Elsevier Ltd. 2010;15:349–56.
Reichenberg A, Caspi A, Harrington H, Houts R, Keefe RSE, Murray RM, et al. Static and dynamic cognitive deficits in childhood preceding adult schizophrenia: a 30-year study. Am J Psychiatry. American Psychiatric Publishing. 2010;167:160–9.
Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature. Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved. 2014;511:421–7.
López M, García-Oguiza A, Armstrong J, García-Cobaleda I, García-Miñaur S, Santos-Simarro F, et al. Rubinstein-Taybi 2 associated to novel EP300 mutations: deepening the clinical and genetic spectrum. BMC Med Genet. 2018;19:36.
McDonald-McGinn DM, Sullivan KE, Marino B, Philip N, Swillen A, Vorstman JAS, et al. 22q11.2 deletion syndrome. Nat Rev Dis Prim. Macmillan Publishers Limited. 2015;1:15071.
Van L, Boot E, Bassett AS. Update on the 22q11.2 deletion syndrome and its relevance to schizophrenia. Curr Opin Psychiatry. 2017;30:191-6.
Cui I, Cui H. Antisense RNAs and epigenetic regulation. Epigenomics [Internet]. Future Medicine. 2010;2:139–50 Available from: https://doi.org/10.2217/epi.09.46.
Tan XL, Zhai Y, Gao WX, Fan YM, Liu FY, Huang QY, et al. p300 expression is induced by oxygen deficiency and protects neuron cells from damage. Brain Res. 2009;1254:1–9.
Peng Y, Yang Z, Zhang H, Cui C, Qi X, Luo X, et al. Genetic Variations in Tibetan Populations and High-Altitude Adaptation at the Himalayas. Mol Biol Evol. 2011;28:1075–81.
van Erp TGM, Saleh PA, Rosso IM, Huttunen M, Lönnqvist J, Pirkola T, et al. Contributions of genetic risk and fetal hypoxia to hippocampal volume in patients with schizophrenia or schizoaffective disorder, their unaffected siblings, and healthy unrelated volunteers. Am J Psychiatry. American Psychiatric Publishing. 2002;159:1514–20.
Cannon TD, van Erp TGM, Rosso IM, Huttunen M, Lönnqvist J, Pirkola T, et al. Fetal hypoxia and structural brain abnormalities in schizophrenic patients, their siblings, and controls. Arch Gen Psychiatry. 2002;59:35–41.
Cannon TD, Yolken R, Buka S, Torrey EF. Decreased neurotrophic response to birth hypoxia in the etiology of schizophrenia. Biol Psychiatry. 2008;64:797–802.
Parets SE, Conneely KN, Kilaru V, Fortunato SJ, Syed TA, Saade G, et al. Fetal DNA methylation associates with early spontaneous preterm birth and gestational age. PLoS One. Public Library of Science. 2013;8:e67489.
Córdova-Palomera A, Alemany S, Fatjó-Vilas M, Goldberg X, Leza JC, González-Pinto A, et al. Birth weight, working memory and epigenetic signatures in IGF2 and related genes: a MZ twin study. PLoS One. Public Library of Science. 2014;9:e103639.
Dempster EL, Pidsley R, Schalkwyk LC, Owens S, Georgiades A, Kane F, et al. Disease-associated epigenetic changes in monozygotic twins discordant for schizophrenia and bipolar disorder. Hum Mol Genet. 2011;20:4786–96.
Essex MJ, Thomas Boyce W, Hertzman C, Lam LL, Armstrong JM, Neumann SMA, et al. Epigenetic vestiges of early developmental adversity: childhood stress exposure and DNA methylation in adolescence. Child Dev. 2013;84:58–75.
Leenen FAD, Muller CP, Turner JD. DNA methylation: conducting the orchestra from exposure to phenotype? Clin Epigenetics [Internet]. 2016;8:92 Available from: https://doi.org/10.1186/s13148-016-0256-8.
Sago H, Ishii K, Sugibayashi R, Ozawa K, Sumie M, Wada S. Fetoscopic laser photocoagulation for twin–twin transfusion syndrome. J Obstet Gynaecol Res. John Wiley & Sons, Ltd (10.1111). 2018;44:831–9.
Robinson HP. Sonar measurement of fetal crown-rump length as means of assessing maturity in first trimester of pregnancy. Br Med J. 1973;4:28 LP–31.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30:1363–9.
McCartney DL, Walker RM, Morris SW, McIntosh AM, Porteous DJ, Evans KL. Identification of polymorphic and off-target probe binding sites on the Illumina Infinium MethylationEPIC BeadChip. Genomics Data. 2016;9:22–4 Elsevier; [cited 2019 Apr 16];. Available from: https://www.sciencedirect.com/science/article/pii/S221359601630071X?via%3Dihub.
Morin AM, Gatev E, McEwen LM, MacIsaac JL, Lin DTS, Koen N, et al. Maternal blood contamination of collected cord blood can be identified using DNA methylation at three CpGs. Clin Epigenetics. 2017;9:75.
R Development Core Team. R: A Language and environment for statistical computing [Internet]. Team RDC, editor. R Found. Stat. Comput. R Foundation for Statistical Computing; 2011. p. 409. Available from: http://www.r-project.org
Simpkin AJ, Suderman M, Howe LD. Epigenetic clocks for gestational age: Statistical and study design considerations. Clin Epigenetics. Clin Epigenetics. 2017;9:1–2.
Mcewen LM, Jones MJ, Tse D, Lin S, Edgar RD, Husquin LT, et al. Systematic evaluation of DNA methylation age estimation with common preprocessing methods and the Infinium MethylationEPIC BeadChip array. Clin Epigenetics. Clin Epigenetics. 2018;10:123.
Houseman E, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13:86.
Córdova-Palomera A, Fatjó-Vilas M, Palma-Gudiel H, Blasco-Fontecilla H, Kebir O, Fañanás L. Further evidence of DEPDC7 DNA hypomethylation in depression: a study in adult twins. Eur Psychiatry. 2015;30:715–8.
We would like to thank Anna Valldeperas, Miriam Osorio, and Talita Micheletti for their technical assistance with the biological sample processing. We are indebted to the IDIBAPS Biobank, integrated in the Spanish National Biobank Network, for the sample and data procurement.
This work was supported by a NARSAD Distinguished Investigator Grant, awarded to Professor Lourdes Fañanás, the Biomedical Research Networking Center for Mental Health (CIBERSAM) and the Comissionat per a Universitats i Recerca del DIUE, Generalitat de Catalunya (grant number 2017SGR1577). This research was also partially funded by CERCA Programme/Generalitat de Catalunya “la Caixa” Foundation (LCF/PR/GN14/10270005), and AGAUR 2017 SGR grant number 1531. The organizations listed above had no further role in study design, in the collection, analysis and interpretation of data, in the writing of the report, and in the decision to submit the paper for publication.
Availability of data and materials
The datasets used and/or analyzed during the current study together with the associated relevant meta-data are available from the corresponding author upon reasonable request.
Ethics approval and consent to participate
Written informed consent was obtained from all participants after a detailed description of the study aims and design, as approved by the Bioethics Committee of the University of Barcelona. All procedures contributing to this work were performed in accordance with the Helsinki Declaration of 1975, as revised in 2008.
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.
DNA methylation values for CpG probes used to discard the presence of maternal contamination. (DOCX 29.6 kb)