- Open Access
Maternal eating disorders affect offspring cord blood DNA methylation: a prospective study
Clinical Epigeneticsvolume 9, Article number: 120 (2017)
Eating disorders (ED) are chronic psychiatric disorders, common amongst women of reproductive age. ED in pregnancy are associated with poor nutrition and abnormal intrauterine growth. Increasing evidence also shows offspring of women with ED have adverse developmental and birth outcomes. We sought to carry out the first study investigating DNA methylation in offspring of women with ED. We compared cord blood DNA methylation in offspring of women with active ED (n = 21), past ED (n = 43) and age- and social class-matched controls (n = 126) as part of the Avon Longitudinal Study of Parents and Children.
Offspring of women with both active and past ED had lower whole-genome methylation compared to controls (active ED 49.1% (95% confidence intervals 50.5–47.7%), past ED 49.2% (95% CI 50.7–47.7.0%), controls 52.4% (95% CI 53.0%–51.0%)). Amongst offspring of ED women, those born to women with restrictive-type and purging-type ED had lower methylation levels compared to those of controls. Offspring of women with an active restrictive ED in pregnancy had lower whole-genome methylation compared to offspring of women with past restrictive ED. We observed decreased methylation at the DHCR24 locus in offspring of women with active pregnancy ED (effect size (ES) = − 0.124, p = 6.94 × 10−8) and increased methylation at the LGALS2 locus in offspring of women with past ED (ES = 0.07, p = 3.74 × 10−7) compared to controls.
Maternal active and past ED are associated with differences in offspring whole-genome methylation. Our results show altered DNA methylation in loci relevant to metabolism; these might be biomarkers of disrupted metabolic pathways in offspring of ED mothers. Further work is needed to examine potential mechanisms and functional outcomes of the observed methylation patterns.
Eating disorders (ED; anorexia nervosa, bulimia nervosa, binge eating disorder, and other specified feeding and eating disorders) are severe psychiatric disorders that affect about 10% of women of reproductive age [1,2,3]. Maternal ED, both active—during pregnancy—and past, have been shown to affect child eating, growth, psychopathology and development [4,5,6,7,8]. However, little is known about biological risk mechanisms . Women with active and lifetime ED differ from non-ED women in terms of their food intake during pregnancy, with differences highlighted in micronutrient and macronutrient intakes [10, 11]. Nutrient intake in pregnancy is affected by ED type, with differences in diet composition seen in women with anorexia nervosa (AN) compared to bulimia nervosa (BN) and binge eating disorder (BED) [10, 11].
Maternal food intake in utero, in terms of quality and quantity, is likely an important environmental determinant affecting foetal metabolic adaptation to the prevailing environmental conditions, as shown by studies on maternal nutritional deprivation as a result of widespread famine [12, 13]. Increasing evidence has highlighted the effect of maternal nutrition in pregnancy on epigenetic modifications, in particular DNA methylation, both in animals and humans [14, 15]. Low micronutrient (methyl donors) and protein intake, as well as low maternal pre-pregnancy BMI, have been prospectively associated with lower global methylation levels and hypomethylation at specific sites [14, 16].
No studies to date have investigated the effect of maternal ED on offspring cord blood methylation. Given that nutrition in pregnancy is likely to differ in women with ED compared to healthy women, our first aim was to carry out an epigenome-wide association study (EWAS) in offspring of mothers with ED and offspring of non-ED mothers. There is evidence that maternal nutrition might affect offspring methylation not just in utero but also preconception ; we therefore also aimed to disentangle the effect of active ED in pregnancy vs. ED at any point prior to pregnancy, by investigating both global methylation levels and methylation at CpG sites in offspring of women with active ED and past ED and offspring of age- and social class-matched controls. Lastly given that ED differ in terms of manifestation of specific behaviours, we explored the effect of maternal restrictive eating, purging (self-induced vomiting and laxative use for weight control) and binge eating on overall offspring methylation.
The Avon Longitudinal Study of Parents and Children (ALSPAC) is a longitudinal, population-based, extensive prospective study of women and their children, set up in the 1990s to investigate the effects of environment, genetic and other factors on child health and development [17, 18]. All pregnant women living in the geographical area of Avon, UK, who were expected to deliver their baby between 1 April 1991 and 31 December 1992 were invited to take part in the study. Uptake was high and those enrolled represented approximately 85% of the eligible population. ALSPAC recruited 14,541 pregnant women; all women gave informed and written consent. The study website contains details of all the data that is available through a fully searchable data dictionary: http://www.bris.ac.uk/alspac/researchers/data-access/data-dictionary/. Data are accessible to all bona fide researchers via the ALSPAC study. However, ALSPAC is required by its ethical approvals to manage access through an application process (which applied equally to the authors of this manuscript).
We selected women whose offspring had available DNA methylation data from cord blood (n = 1018) as part of the Accessible Resource for Integrated Epigenomic Studies (ARIES; for details, see http://www.ariesepigenomics.org.uk/), a representative sub-set of ALSPAC [18,19,20]. Offspring were not excluded from the study based on gestational age, prematurity, low birth weight and non-singleton pregnancies.
Exposure (maternal ED)
Amongst women with offspring methylation data from cord blood (n = 1018), maternal active ED was determined via three sources: (a) questionnaires at 12 weeks enquiring about having had a recent ED; (b) questionnaires at 18 and 32 weeks in pregnancy on ED behaviours and cognitions [including specific questions about loss of control eating, purging (self-induced vomiting or laxative use for weight loss) and the weight and shape concern subscales of the Eating Disorder Examination Questionnaire (EDE-q), a well-validated instrument , as previously described in [21, 22]; and (c) validated semi-structured interviews recently carried out as part of a two-phase lifetime ED prevalence study .
Women were defined as having had an active ED in pregnancy (n = 21) if they (a) reported having had an active ED in pregnancy (AN or BN) (n = 11) or (b) reported symptoms consistent with an ED in pregnancy (loss of control eating, purging, high weight and shape concern in the context of underweight (n = 10)); or (c) reported symptoms consistent with an ED diagnosis in pregnancy on interviews (n = 7) (see Table 1). The majority (52.4%) of women were categorised as having active ED in pregnancy based on two out of three of the above sources.
Women who self-reported having suffered from an ED prior to pregnancy, had no ED symptoms in pregnancy, and did not report any symptoms consistent with an ED in pregnancy were categorised as having suffered from a past ED (n = 43).
We also categorised women according to broad ED diagnostic type (restrictive, purging, binge eating, binge-purge) according to symptoms reported during and prior to pregnancy. Amongst women with active ED, 14 (66.7%) had a restrictive ED, 3 (14.3%) a binge-purge-type ED, 3 (14.3%) a purge-type ED and 1 (4.8%) a binge eating ED. Amongst women with past ED, 25 had a restrictive-type ED (58.1%), 10 (23.3%) a binge-purge-type ED and 4 (9.3%) a purge-type ED.
Two age- and social class-matched controls were selected for each case (n = 126).
Socio-demographic (maternal social class, maternal age) data, smoking in pregnancy, and weight and height pre-pregnancy were obtained by self-completion questionnaire at 12, 18 and 32 weeks’ gestation; body mass index (BMI) was calculated as pre-pregnant weight/height2. Total energy intake was collected via food frequency questionnaire (FFQ), completed by the women at 32 weeks’ gestation . Offspring gender was obtained from obstetric records.
Methylation measurements and quality control
DNA was extracted and was bisulfite converted using the Zymo EZ DNA MethylationTM kit (Zymo, Irvine, CA). After the conversion, DNA methylation level was measured using the Illumina Infinium® HumanMethylation450k BeadChip in line with standard protocols at the University of Bristol, UK.
The arrays were scanned using an Illumina iScan, and initial quality review was assessed using GenomeStudio (version 2011.1). Further quality checks included CpG detection p value test which removed the samples whose average detection p value is > 0.01. We also used several study specific variables and created multidimensional scaling (MDS) plots to identify any potential outliers, but all samples were within the thresholds. The measures of DNA methylation were computed as betas (β-value) using the minfi package in the R statistical programming language . Methylation beta values range from 0 (no cytosine methylation) to 1 (completely cytosine methylated). Raw beta values were normalised using the functional normalisation method of the minfi package . We excluded SNPs, probes with a detection p value > 0.05 in more than 5% samples and features on sex chromosomes. After excluding these features, 468,782 probes were left for analysis.
We used multiple linear regression models to explore the effect of maternal ED on offspring cord blood methylation using the CpGassoc R package .
Firstly, cord blood in offspring of mothers with active ED, past ED and controls were compared. We profiled overall methylation levels amongst these three groups. We calculated a median beta value for each individual and then computed a median beta for each exposure group for comparisons. Crude and adjusted analyses were performed. Adjusted analyses included bisulfite conversion batch (BCD) to adjust for batch effects, and principal components (computed using principal component analysis (PCA) [26,27,28]. PCA was carried out to adjust for any unknown batch effects. Only those principal components (PC) were included which captured most of the variance and were not associated with the exposure. Maternal smoking status and child gender were included in all adjusted analyses as covariates. These covariates were used as fixed effect to minimise the computational burden.
Cell-type correction was applied using the reference-based Houseman method  in the minfi package in R. This method estimates the relative proportions of six white blood cell subtypes (CD4+ T lymphocytes, CD8+ T lymphocytes, NK (natural killer) cells, B lymphocytes, monocytes and granulocytes), based on a standard reference population . Effect sizes and % methylation are reported.
We analysed four linear regression models including model 1: a crude regression of offspring DNA methylation on maternal ED group (as a predictor); model 2: as model 1 with adjustment for BCD and first six PCs (batches); model 3: as model 2 with additional adjustment for covariates (smoking and child gender); and model 4: as model 3 with additional adjustment for cell counts. We show results of Model 4, the fully adjusted model, throughout the paper. In EWAS results, we highlighted cross-hybridising probes as these probes can interfere with accurate detection of methylation levels .
For EWAS analyses, we computed the genomic inflation factor and generated quantile-quantile (QQ), Manhattan and volcano plots to compare the genome-wide distribution of p values and effect sizes across different analyses.
Methylation of two CpGs sites (uncorrected p value < 10−7) identified in fully adjusted EWAS analysis of active ED versus controls was profiled in the three exposure groups: controls, active ED, and past ED.
In order to determine whether any effects of active ED might be due to maternal underweight, we carried out sensitivity analyses to compare our top hits with those found in cord blood of offspring of underweight mothers in ALSPAC  and in a meta-analysis of sustained maternal smoking during pregnancy and cord blood methylation .
Linear regression analyses were performed to investigate the effects of specific ED type on offspring cord blood methylation after removing cross-hybridising probes (n = 440,643).
Post-hoc analyses investigated whether cord blood methylation differed in offspring of women with active vs. past restrictive-type ED.
Overall, 189 women (n = 21 with active ED, n = 43 with past ED, n = 125 healthy controls) had data on exposure, outcomes and relevant confounders under study. One healthy control was dropped from the study sample due to missing data on smoking in pregnancy.
Socio-demographic characteristics of the women under study are shown in Table 1. Women with current and past ED had lower average daily energy intake compared to controls (F = 2.8, p = 0.06). The majority of women in the active and past ED group had a restrictive-type ED.
We removed cross-hybridising probes and calculated medians in order to plot and compare the overall methylation levels across the three groups (controls, active ED and past ED) (Additional file 1: Figure S1). Cord blood of offspring of women with active ED had the lowest methylation level (49.1%; 95% confidence intervals 50.5–47.7%); offspring of controls had the highest global methylation levels (52.4%; 95% CI 53.0–51.0%). The offspring of women with past ED had a lower global methylation level compared to those of controls (49.2%; 95% CI 50.7–47.7.0%) and slightly higher than the active ED group. An ANOVA showed medians were statistically different across groups (p = 0.008).
Identification of differentially methylated CpGs
We carried out EWAS analyses to identify the association between maternal eating disorder and epigenome-wide DNA methylation in offspring cord blood measured using the Illumina HumanMethylation450 Beadchip. After adjustment for potential confounders and estimated white cell counts, differences were shown in two CpG sites between offspring of women with active ED and controls (all p < 10−7) (see Table 2) (Additional file 1: Figure S2 shows the relevant QQ and Manhattan plot), and in two CpGs between offspring of women with past ED and controls (all p < 10−6) (see Table 2). The top two CpGs identified in analysis of active ED and controls also passed Bonferroni correction (1.06 × 10−7 corrected for 468,782 associations). No significantly different CpGs were identified in offspring of mothers with active vs. past ED.
In particular, offspring of women with active ED had a hypomethylated region in cg10177197, mapping to the DHCR24 gene, which corresponds to the gene encoding 3-β-hydroxysterol delta-24 reductase, an enzyme that catalyses the last step in the biosynthesis of cholesterol (also known as seladin 1). The second CpG site identified was cg21146184 on chromosome 22, which is near the ATXN10 gene (the gene encoding Ataxin 10—a protein that may be involved in neuron survival, neuron differentiation and neuritogenesis) .
Two loci were differentially methylated in offspring of women with past ED vs. controls, and these trended for significance following Bonferroni correction: cg05963962, corresponding to LOC84856 (LINC00839 a long intergenic non-protein coding element) on chromosome 10, and cg11081833, corresponding to a soluble beta-galactoside binding lectin (LGALS2) gene. This gene and its polymorphisms have been implicated in the metabolic syndrome , inflammation and immune response. Additional file 1: Figure S3 and Figure S4 shows methylation levels across groups in these CpGs.
In order to determine whether weight status partially explained our findings, we compared our EWAS top CpG sites (p value < 10−4) with results from a recent ALSPAC study on the effects of maternal over- and underweight and offspring cord blood methylation . We considered the fully adjusted model reported in  in relation to maternal pre-pregnancy BMI (n = 284,972 probes) to compare the CpG sites identified in our analysis of active ED versus controls, and our top two hits were not present in . The direction of effect was not consistent for the majority of the sites, the previously reported effect sizes were smaller, and p values were larger ranging from 0.01 to 0.99 (Additional file 2: Table S1). We also compared the CpG sites identified in our study in offspring of past ED versus controls with CPGs identified in , and the top CpG sites in our study were not found in the previously reported study . In agreement with the previous comparison, the effect sizes in  were smaller and p values were larger ranging from 0.07 to 0.9 (Additional file 2: Table S2).
We compared the CpG sites identified in fully adjusted EWAS of active ED versus controls (p value < 10−4) with FDR significant sites found in a study of sustained maternal smoking in pregnancy (n = 6074 probes), and our top two hits were not present in . There were only three overlapping CpGs (Additional file 2: Table S3). We also compared the CpG sites identified in our study in offspring of past ED versus controls with CPGs identified in  and four CpGs overlapped that include our second top site (Additional file 2: Table S4).
We investigated the association between ED diagnostic type (restrictive ED, purging ED, and binge eating and purging ED)—across active and past ED—and whole-genome methylation (results are shown in Table 3). In adjusted analyses, all ED diagnostic types predicted lower offspring cord blood methylation levels compared to controls: restrictive-type ED (ß = − 0.011, p = 0.15), purging-type ED (ß = − 0.026, p = 0.08), and binge + purge-type ED (− 0.009, p = 0.4) (see Table 3). Results were comparable when cross-hybridising probes were removed.
As the restrictive ED group was the largest, we were able to also explore differences in whole-genome methylation amongst women with active restrictive-type ED and past restrictive-type ED. Offspring of women with active restrictive-type ED had significantly lower average whole-genome methylation (ß = − 0.07 (95% CI: − 0.069, − 0.077), p = 0.01) compared to women with past restrictive-type ED in adjusted analyses.
This is the first study to investigate cord blood methylation in offspring born to mothers with eating disorders. As such, it provides novel and unique preliminary evidence of likely effects of maternal ED on global methylation and specific methylation signatures in domains relevant to metabolism and neuronal development. This study builds not only on existing evidence on effects of maternal ED on child development but also on novel findings highlighting the importance of genes relevant to metabolism in risk for ED [34, 35], and the first locus to be identified as being significantly associated with AN at genome-wide level .
Maternal eating disorders are associated with negative developmental outcomes across a range of behavioural and physical aspects. Pregnant women with active eating disorders also differ from healthy women in terms of their eating, diet and nutritional intakes and stress/anxiety levels [10, 11, 37,38,39]. Hence, we hypothesised that maternal eating disorders, especially if active during pregnancy, would be associated with differential methylation patterns in offspring cord blood. Offspring of mothers with active ED had in fact lower global methylation, by 3.3%, compared to offspring of age- and socio-demographic status-matched non-ED women. Interestingly, offspring of women with past ED had global methylation levels slightly higher than offspring of women with active ED, but nevertheless 3% lower than offspring of non-ED women. Children of women with restrictive-type ED and purge-type ED had the largest differences in global methylation compared to controls, although the number of women with purging behaviours only was small. In line with results indicating active pregnancy ED was associated with lowest global methylation levels, exploratory analyses indicated that restrictive eating in pregnancy was associated with lower offspring cord blood methylation compared to restrictive eating prior to pregnancy. We also identified intriguing significant differentially methylated CpGs in regions/genes relevant to metabolism, and neuronal development were associated with both maternal active and past ED.
Maternal diet and nutritional intake are amongst the major factors influencing the developing foetus . Methylation depends on the availability of methyl group donors, mostly derived from diet. Therefore, hypomethylation has been shown to ensue diets with multiple methyl-related compounds nutrient deficiencies . Epigenetic mechanisms have been suggested as potential mediating factors between maternal exposures (under and over-nutrition) and offspring adiposity. Evidence from the ALSPAC study also suggests that cord blood global hypomethylation is associated with maternal underweight and that hypomethylated sites [associated with maternal underweight] were negatively associated with offspring adiposity in later life. This suggests a possible mechanism for intergenerational foetal programming. Given the overlap between underweight and some ED (particularly restrictive ones), the fact that our top hypomethylated sites were not present amongst those associated with maternal underweight in ALSPAC strengthens our hypotheses that the effects seen in our study might not simply be secondary to maternal underweight. Cord blood hypomethylation has also been shown as a result of maternal pregnancy smoking  in humans. Although women with ED have higher prevalence of smoking in pregnancy [7, 22], smoking was adjusted for in all analyses. Corroboration that the effect of maternal ED on cord blood methylation might be driven by eating behaviours or their epiphenomena comes from our findings on ED-type showing lower methylation levels in offspring of mothers with restrictive- and purge-type ED (and within the restrictive ED subgroup, those women who had restrictive ED in pregnancy).
Cord blood methylation has also been suggested as a potential biomarker of altered child neurodevelopment , and research linking specific gene methylation and child neurodevelopmental outcomes has shown promising results. Our study highlights interesting differences in differentially methylated CpGs between offspring of ED mothers and controls. In particular, we identified differential methylation in genes relevant to lipid synthesis (DHCR24, LGALS2) . DHCR24 is responsive to lipid levels  and oestrogen , both likely to be abnormal in active ED and therefore possibly relevant to foetal programming. In fact, there is evidence that expression of DHCR24 is relevant to foetal programming. DHCR24 hepatic transcription is differentially induced by maternal undernutrition and low-protein diet in rats  and mice ; expression of DHCR24 and related genes in the placenta was associated with foetal growth . Two differentially methylated CpGs identified were in and near genes relevant to general pathways of normal embryonic and placental development (LGALS2) and neurogenesis (ATXN10). These findings provide initial preliminary insights into possible biological effects of maternal ED; mapping onto existing evidence of both (a) developmental effects of maternal ED  and (b) risk pathways likely to be relevant to the intergenerational cycle of risk of ED .
Our findings have to be interpreted in the context of strengths and limitations. This is the first study to investigate methylation in offspring of women with ED. ED were ascertained using a triangulation method relying on self-report and interviews; we were able to test the effect of active ED vs. past ED. Given the nature of the sample, we were able to match controls to cases and adjust for confounders known to affect methylation levels, e.g. smoking. This study, however, is limited by the relatively small sample size, which nevertheless allowed exploration of the effect of ED type on the outcome. Other limitations include that no technical validation (e.g. pyrosequencing) was performed; studies have shown that the Illumina microarrays are less sensitive to extreme methylation levels [49, 50].
This study provides for the first time evidence of the effect of maternal ED on early biological markers of offspring development, i.e. cord blood methylation. We showed important effects on genome-wide methylation levels and methylation of specific genes relevant to metabolism. Replication of our findings will be important. Recent studies have shown persistence of methylation levels into childhood ; hence, understanding whether differences identified in methylation are persistent and whether they are associated with adverse developmental outcomes will allow a more in-depth understanding of developmental risk pathways relevant to the intergenerational transmission of ED and broader psychopathology.
Avon longitudinal study of parents and children
Bisulfite conversion batch
Binge eating disorder
Body mass index
Epigenome-wide association study
Micali N, Hagberg KW, Petersen I, Treasure JL. The incidence of eating disorders in the UK in 2000-2009: findings from the General Practice Research Database. BMJ Open. 2013;3(5)
Easter A, Bye A, Taborelli E, Corfield F, Schmidt U, Treasure J, et al. Recognising the symptoms: how common are eating disorders in pregnancy? Eur Eat Disord Rev. 2013;21(4):340–4.
Micali N, Martini MG, Thomas JJ, Eddy KT, Kothari R, Russell E, et al. Lifetime and 12-month prevalence of eating disorders amongst women in mid-life: a population-based study of diagnoses and risk factors. BMC Med. 2017;15(1):12.
Easter A, Howe LD, Tilling K, Schmidt U, Treasure J, Micali N. Growth trajectories in the children of mothers with eating disorders: a longitudinal study. BMJ Open. 2014;4(3):e004453.
Easter A, Naumann U, Northstone K, Schmidt U, Treasure J, Micali N. A longitudinal investigation of nutrition and dietary patterns in children of mothers with eating disorders. J Pediatr. 2013;163(1):173–8. e1
Micali N, Simonoff E, Stahl D, Treasure J. Maternal eating disorders and infant feeding difficulties: maternal and child mediators in a longitudinal general population study. J Child Psychol Psychiatry. 2011;52(7):800–7.
Micali N, Stemann Larsen P, Strandberg-Larsen K, Nybo Andersen AM. Size at birth and preterm birth in women with lifetime eating disorders: a prospective population-based study. BJOG. 2016;123(8):1301–10.
Micali N, De Stavola B, Ploubidis GB, Simonoff E, Treasure J. The effects of maternal eating disorders on offspring childhood and early adolescent psychiatric disorders. Int J Eat Disord. 2014;47(4):385–93.
Micali N, Treasure J. Biological effects of a maternal ED on pregnancy and foetal development: a review. Eur Eat Disord Rev. 2009;17(6):448–54.
Micali N, Northstone K, Emmett P, Naumann U, Treasure JL. Nutritional intake and dietary patterns in pregnancy: a longitudinal study of women with lifetime eating disorders. Br J Nutr. 2012;108(11):2093–9.
Siega-Riz AM, Haugen M, Meltzer HM, Von Holle A, Hamer R, Torgersen L, et al. Nutrient and food group intakes of women with and without bulimia nervosa and binge eating disorder during pregnancy. Am J Clin Nutr. 2008;87(5):1346–55.
Roseboom TJ, van der Meulen JH, Ravelli AC, Osmond C, Barker DJ, Bleker OP. Effects of prenatal exposure to the Dutch famine on adult disease in later life: an overview. Twin Res. 2001;4(5):293–8.
Schulz LC. The Dutch Hunger Winter and the developmental origins of health and disease. Proc Natl Acad Sci U S A. 2010;107(39):16757–8.
Dominguez-Salas P, Moore SE, Baker MS, Bergen AW, Cox SE, Dyer RA, et al. Maternal nutrition at conception modulates DNA methylation of human metastable epialleles. Nat Commun. 2014;5:3746.
Zhang N. Epigenetic modulation of DNA methylation by nutrition and its mechanisms in animals. Animal Nutrition. 2015;1(3):144–51.
Sharp GC, Lawlor DA, Richmond RC, Fraser A, Simpkin A, Suderman M, et al. Maternal pre-pregnancy BMI and gestational weight gain, offspring DNA methylation and later offspring adiposity: findings from the Avon Longitudinal Study of Parents and Children. Int J Epidemiol. 2015;44(4):1288–304.
Golding J, Pembrey M, Jones R, Team AS. ALSPAC—the Avon Longitudinal Study of Parents and Children. I. Study methodology. Paediatr Perinat Epidemiol. 2001;15(1):74–87.
Boyd A, Golding J, Macleod J, Lawlor DA, Fraser A, Henderson J, et al. Cohort profile: the 'children of the 90s'—the index offspring of the Avon Longitudinal Study of Parents and Children. Int J Epidemiol. 2013;42(1):111–27.
Fraser A, Macdonald-Wallis C, Tilling K, Boyd A, Golding J, Davey Smith G, et al. Cohort profile: the Avon Longitudinal Study of Parents and Children: ALSPAC mothers cohort. Int J Epidemiol. 2013;42(1):97–110.
Relton CL, Gaunt T, McArdle W, Ho K, Duggirala A, Shihab H, et al. Data resource profile: Accessible Resource for Integrated Epigenomic Studies (ARIES). Int J Epidemiol. 2015;44(4):1181–90.
Micali N, Treasure J, Simonoff E. Eating disorders symptoms in pregnancy: a longitudinal study of women with recent and past eating disorders and obesity. J Psychosom Res. 2007;63(3):297–303.
Micali N, Simonoff E, Treasure J. Risk of major adverse perinatal outcomes in women with eating disorders. Br J Psychiatry. 2007;190:255–9.
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 (Oxford, England). 2014;30(10):1363–9.
Fortin JP, Labbe A, Lemire M, Zanke BW, Hudson TJ, Fertig EJ, et al. Functional normalization of 450k methylation array data improves replication in large cancer studies. Genome Biol. 2014;15(12):503.
Barfield R, Conneely K, Kilaru V. CpGassoc: association between methylation and a phenotype of interest. R package version 2.55. 2016.
Becker RA, Champbers JM. The new S language. Wilks AR: Wadsworth & Brooks/Cole; 1988.
Mardia KV, Kent JT, Bibby JM. Multivariate analysis. London: Academic Press; 1979.
Venables WN, Ripley BD. Modern applied statistics with S. 2002.
Houseman EA, 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.
Reinius LE, Acevedo N, Joerink M, Pershagen G, Dahlen SE, Greco D, et al. Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PLoS One. 2012;7(7):e41361.
Chen YA, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics. 2013;8(2):203–9.
Joubert BR, Felix JF, Yousefi P, Bakulski KM, Just AC, Breton C, et al. DNA Methylation in newborns and maternal smoking in pregnancy: genome-wide consortium meta-analysis. Am J Hum Genet. 2016;98(4):680–96.
Christensen MB, Lawlor DA, Gaunt TR, Howell WM, Davey Smith G, Ebrahim S, et al. Genotype of galectin 2 (LGALS2) is associated with insulin-glucose profile in the British Women's Heart and Health Study. Diabetologia. 2006;49(4):673–7.
Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh PR, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015;47(11):1236–41.
Hinney A, Kesselmeier M, Jall S, Volckmar AL, Focker M, Antel J, et al. Evidence for three genetic loci involved in both anorexia nervosa risk and variation of body mass index. Molecular Psychiatry. 2017;22(2):321–2. doi:10.1038/mp.2016.126. Epub 2016 Jul 26.
Duncan L, Yilmaz Z, Gaspar H, Walters R, Goldstein J, Anttila V, et al. Significant Locus and Metabolic Genetic Correlations Revealed in Genome-Wide Association Study of Anorexia Nervosa. Am J Psychiatry. 2017;174(9):850–8. doi:10.1176/appi.ajp.2017.16121402.
Easter A, Solmi F, Bye A, Taborelli E, Corfield F, Schmidt U, et al. Antenatal and postnatal psychopathology among women with current and past eating disorders: longitudinal patterns. Eur Eat Disord Rev. 2015;23(1):19–27.
Easter A, Taborelli E, Bye A, Zunszain PA, Pariante CM, Treasure J, et al. Perinatal hypothalamic-pituitary-adrenal axis regulation among women with eating disorders and their infants. Psychoneuroendocrinology. 2017;76:127–34.
Micali N, Simonoff E, Treasure J. Pregnancy and post-partum depression and anxiety in a longitudinal general population cohort: the effect of eating disorders and past depression. J Affect Disord. 2011;131(1–3):150–7.
Chmurzynska A. Fetal programming: link between early nutrition, DNA methylation, and complex diseases. Nutr Rev. 2010;68(2):87–98.
Jacob RA. Folate, DNA methylation, and gene expression: factors of nature and nurture. Am J Clin Nutr. 2000;72(4):903–4.
Guerrero-Preston R, Goldman LR, Brebi-Mieville P, Ili-Gangas C, LeBron C, Hernandez-Arroyo M, et al. Global DNA hypomethylation is associated with in utero exposure to cotinine and perfluorinated alkyl compounds. Epigenetics. 2010;5(6):539–46.
Hodyl NA, Roberts CT, Bianco-Miotto T. Cord Blood DNA Methylation Biomarkers for Predicting Neurodevelopmental Outcomes. Genes (Basel). 2016;7(12).
Daimiel Lidia A, Fernández-Suárez María E, Rodríguez-Acebes S, Crespo L, Lasunción Miguel A, Gómez-Coronado D, et al. Promoter analysis of the DHCR24 (3β-hydroxysterol Δ(24)-reductase) gene: characterization of SREBP (sterol-regulatoryelement-binding protein)-mediated activation. Biosci Rep. 2013;33(1):e00006.
Dekkers KF, van Iterson M, Slieker RC, Moed MH, Bonder MJ, van Galen M, et al. Blood lipids influence DNA methylation in circulating cells. Genome Biol. 2016;17(1):138.
Ellis PJ, Morris TJ, Skinner BM, Sargent CA, Vickers MH, Gluckman PD, et al. Thrifty metabolic programming in rats is induced by both maternal undernutrition and postnatal leptin treatment, but masked in the presence of both: implications for models of developmental programming. BMC Genomics. 2014;15(1):49.
Zheng J, Xiao X, Zhang Q, Wang T, Yu M, Xu J. Maternal low-protein diet modulates glucose metabolism and hepatic microRNAs expression in the early life of offspring †. Nutrients. 2017;9(3):205.
Kappil MA, Green BB, Armstrong DA, Sharp AJ, Lambertini L, Marsit CJ, et al. Placental expression profile of imprinted genes impacts birth weight. Epigenetics. 2015;10(9):842–9.
Dedeurwaerder S, Defrance M, Calonne E, Denis H, Sotiriou C, Fuks F. Evaluation of the Infinium Methylation 450K technology. Epigenomics. 2011;3(6):771–84.
Pidsley R, Zotenko E, Peters TJ, Lawrence MG, Risbridger GP, Molloy P, et al. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 2016;17(1):208.
Herbstman JB, Wang S, Perera FP, Lederman SA, Vishnevetsky J, Rundle AG, et al. Predictors and consequences of global DNA methylation in cord blood and at three years. PLoS One. 2013;8(9):e72824.
We are extremely grateful to all the families who took part in this study, the midwives for their help in recruiting them, and the whole ALSPAC team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists and nurses.
This work was supported by a National Institute of Health Research (NIHR) clinician scientist award to Dr. N Micali (DHCS/08/08/012), by a Wellchild project grant (N/A), and by the UK Medical Research Council (MC_UU_12013/2 and MC_UU_12013/8). The UK Medical Research Council and the Wellcome Trust (Grant ref.: 092731) and the University of Bristol provide core support for ALSPAC.
Availability of data and materials
Data are available from the ALSPAC study; access is subject to ALSPAC data access policies.
Ethics approval and consent to participate
Ethical approval for the study was obtained from the ALSPAC Ethics and Law Committee and the Local Research Ethics Committees. All women gave informed consent at enrolment in the study.
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.
Whole-genome cord blood methylation in offspring of active ED, past ED, and control women. Figure S2. QQ and Manhattan plots: association between ED and offspring cord blood DNA methylation1. Figure S3. Methylation levels of cg21146184 with p value < 10−7 identified in fully adjusted EWAS analysis of active ED versus controls. Figure S4. Methylation levels of cg10177197 with p value < 10−7 identified in fully adjusted EWAS analysis of active ED versus controls. (DOCX 398 kb)
Comparison of top CpG sites (p-value < 10−4) found in EWAS of active ED versus controls with results from ALSPAC study of maternal pre-pregnancy BMI. Table S2. Comparison of top CpG sites (p-value < 10−4) found in EWAS of past ED versus controls with results from ALSPAC study of maternal pre-pregnancy BMI. Table S3. Comparison of top CpG sites (p-value < 10−4) found in EWAS of active ED versus controls with results of sustained maternal smoking in pregnancy. Table S4. Comparison of top CpG sites (p-value < 10−4) found in EWAS of past ED versus controls with results of sustained maternal smoking in pregnancy. (DOCX 49.2 kb)