DNA methylation of imprinted genes at birth is associated with child weight status at birth, 1 year, and 3 years
Clinical Epigenetics volume 10, Article number: 90 (2018)
This study assessed the associations between nine differentially methylated regions (DMRs) of imprinted genes in DNA derived from umbilical cord blood leukocytes in males and females and (1) birth weight for gestational age z score, (2) weight-for-length (WFL) z score at 1 year, and (3) body mass index (BMI) z score at 3 years.
We conducted multiple linear regression in n = 567 infants at birth, n = 288 children at 1 year, and n = 294 children at 3 years from the Newborn Epigenetics Study (NEST). We stratified by sex and adjusted for race/ethnicity, maternal education, maternal pre-pregnancy BMI, prenatal smoking, maternal age, gestational age, and paternal race. We also conducted analysis restricting to infants not born small for gestational age.
We found an association between higher methylation of the sequences regulating paternally expressed gene 10 (PEG10) and anthropometric z scores at 1 year (β = 0.84; 95% CI = 0.34, 1.33; p = 0.001) and 3 years (β = 1.03; 95% CI = 0.37, 1.69; p value = 0.003) in males only. Higher methylation of the DMR regulating mesoderm-specific transcript (MEST) was associated with lower anthropometric z scores in females at 1 year (β = − 1.03; 95% CI − 1.60, − 0.45; p value = 0.001) and 3 years (β = − 1.11; 95% CI − 1.98, − 0.24; p value = 0.01). These associations persisted when we restricted to infants not born small for gestational age.
Our data support a sex-specific association between altered methylation and weight status in early life. These methylation marks can contribute to the compendium of epigenetically regulated regions detectable at birth, influencing obesity in childhood. Larger studies are required to confirm these findings.
Understanding factors that influence the risk of obesity in children is crucial to the development of new strategies for obesity prevention. Obesity in early childhood is a risk factor for obesity later in life [1,2,3] and for a number of chronic diseases in both childhood  and adulthood . Birth weight has been associated with weight outcomes later in life, particularly for those who are on the extremes of the birth weight distribution [6,7,8,9]. Early identification of obesity or its risk factors will inform interventions to prevent the progression of obesity and its consequences later in life . Consistent with the developmental origins of disease hypothesis, the intrauterine environment is hypothesized to influence an individual’s later susceptibility for chronic diseases , including obesity [12, 13].
Epigenetic modifications have been proposed as a mechanism for the in utero origin of later obesity, and a growing literature has found supporting evidence [14,15,16]. DNA methylation is the most studied epigenetic mechanism in humans, due in part, to its stability. DNA methylation that controls the monoallelic expression of imprinted genes is established during gametogenesis and is stably maintained throughout somatic division [17,18,19,20,21] and therefore provides a stable “register” of early in utero exposures. A study of famine survivors found that adults who experienced famine in utero had hypo-methylation of the imprinted IGF2 gene compared to their same sex siblings who had not experienced famine in utero . The significance of this locus was reported to not have been replicated in this cohort using alternate techniques, including RRBS. However, this technique is generally biased toward CG-rich areas and may not have covered the specific and limited number of CpGs that comprise the IGF2 DMR. Additional genes, such as INSR and CPT1A, have also been identified in association with exposure to the Dutch famine . Another study found that maternal nutrition, affected by striking seasonal variations in food intake in the Gambia, influenced methylation at RBM46 . A colorectal cancer study found that methylation status of the IGF2/H19 imprinted locus of adult controls was maintained 3 years later . Moreover, a study of NEST children between birth and age 1 year found similar results at the IGF2/H19 locus .
Select imprinted genes have been identified as playing a role in the development of fetal over and undergrowth caused by imprinting defects. The IGF2 locus is used in clinical diagnostic settings to identify Beckwith-Wiedemann syndrome, which is characterized by overgrowth , and the H19 locus has been used in the diagnosis of Silver-Russell syndrome (SRS), which is characterized by undergrowth .
Although epigenetic data linking DNA methylation and childhood obesity has increased exponentially in the last 5 years [29, 30], few regions agnostically identified have been replicated. This could be in part due to differences in the ethnic composition; however, differences could also be due to the sex composition. At imprinted loci, weight has been associated with the IGF2 locus. Studies have found a relationship between the IGF2 domain and fetal growth [31,32,33,34] and children’s body composition or weight [32, 35, 36]. Data with directional consistency in associations between additional differentially methylated regions (DMRs) and weight gain are required.
This study aims to assess the association between methylation at nine DMRs of imprinted genes and birth weight for gestational age (BW/GA) z score, weight-for-length (WFL) z score at 1 year, and BMI z score at 3 years. In this analysis, we include the following DMRs: MEG3 and MEG3-IG, which are involved in regulating the delta-like 1 homolog/maternally expressed gene 3 imprinted domain on chromosome 14q32.2; IGF2 and H19, which are involved in the imprinting of the insulin growth factor 2/H19 domain on chromosome 11p15, which are located upstream of the imprinted promoters of IGF2 and at the imprinting control region for the IGF2/H19 imprinted domain near the H19 promoter, respectively; PLAGL1 at the pleiomorphic adenoma gene-like 1 locus at 6q24.2; MEST at the mesoderm-specific transcript promoter at 7q32.2; NNAT at the neuronatin locus at 20q11.23; PEG3 at the paternally expressed gene 3 promoter region at 19q13.43; and PEG10 at the epsilon sarcoglycan and paternally expressed gene 10 promoter region at 7q21.3. We selected these regions for their association with infant and child growth [32, 34, 37], chronic disease [22, 38], and parental obesity .
Study sample and data collection
We included data from mothers and children in the Newborn Epigenetic Study (NEST). We have described recruitment and enrollment strategies in detail elsewhere . Briefly, between 2009 and 2011, we recruited women from five prenatal clinics and obstetric facilities in Durham, North Carolina. Eligibility criteria included being at least 18 years of age and intention to use one of the qualifying obstetric facilities for delivery. We excluded women if they planned to relinquish custody of the child or planned to move away from the area in the following 3 years. We obtained written informed consent from all participating women. Upon enrollment, mothers completed questionnaires providing information on sociodemographic factors, lifestyle characteristics, and anthropometrics. At delivery, study personnel abstracted birth outcomes from medical records and infant cord blood specimens were obtained to assess offspring methylation. At 1 year, we collected data on child anthropometrics, feeding, and lifestyle. This study was approved by the Institutional Review Board at Duke University Medical Center.
Of the 1700 enrolled, we excluded 396 women for reasons including miscarriage, refusing further participation, moving away from the area, or delivering at a hospital not included in the study. We analyzed DNA methylation data for the first 600 infants in the study. Infants with analyzed DNA methylation were not significantly different than infants whose DNA methylation had been analyzed with respect to race, maternal education, maternal smoking status, maternal pre-pregnancy BMI, maternal age, or weight at age 1 (data not shown).
Among infants with DNA methylation data, birth weight and length measurements were available for 594. At age 1, we used available weight and length measurements for 306 infants, and at age 3, we used available weight and height measurements for 314 children. We calculated BW/GA z scores using an international standard . We classified infants with BW/GA below the 10th percentile as small for gestational age (SGA). We calculated WFL z scores at age 1 year using WHO standards for children’s exact age . We then calculated BMI z scores at age 3 years using CDC standards . We excluded 7 children with a WFL z score greater than 5 or less than − 5 at age 1 year, and 11 children with a BMI z score greater than 5 or less than − 5. In addition, we excluded infants with possible growth disorders or imprinting defects; therefore, we excluded from analysis infants with DNA methylation values ± 4 standard deviations from the mean (n = 6). The current study includes children with available DNA methylation data at birth on at least one of the nine DMRs of interest, and who had plausible length and weight measurements at birth (n = 576), age 1 (n = 288), or age 3 (n = 294). Plausible weight and length was defined as a measurement that fell within the SD limits set for the combined WFL or BMI and that clearly was not a transcript error (e.g., birth weight being copied onto 1-year weight). In addition, we conducted supplemental analysis on 166 children who had non-missing anthropometric values at birth and age 1 and 3 years to assess directional consistency over time.
Specimen collection and DNA methylation methods have been described in detail elsewhere . Briefly, we collected infant cord blood specimens at birth. We collected samples in EDTA-containing vacutainer tubes and centrifuged using standard protocols to allow for collection of plasma and buffy coat, with buffy coat used for DNA extraction (Qiagen; Valencia, CA). We stored specimens at − 80 °C until the time of analysis. We extracted DNA using Puregene reagents according to the manufacturer’s protocol (Qiagen) and assessed quantity and quality using a Nanodrop 1000 Spectrophotometer (Thermo Scientific; Wilmington, DE).
We modified infant genomic DNA (800 ng) by treatment with sodium bisulfite using the EZ DNA Methylation kit (Zymo Research; Irvine, CA). Bisulfite treatment of denatured DNA converts all unmethylated cytosines to uracils, leaving methylated cytosines unchanged, allowing for quantitative measurement of cytosine methylation status. We performed pyrosequencing using a PyroMark Q96 MD pyrosequencer (Qiagen). Pyrosequencing assay design, genomic coordinates, assay conditions, and assay validation are described in detail elsewhere . Briefly, we designed assays to query established imprinted gene DMRs using the PyroMark Assay Design Software (Qiagen). We optimized PCR conditions to produce a single, robust amplification product. We used defined mixtures of fully methylated and unmethylated control DNAs to show a linear increase in detection of methylation values as the level of input DNA methylation increased (Pearson r is 0.99 for all DMRs). Once we defined optimal conditions, we analyzed each DMR using the same amount of input DNA from each specimen (40 ng, assuming complete recovery following bisulfite modification of 800 ng DNA). We determined percentage of methylation for each CpG cytosine using Pyro Q-CpG software (Qiagen). We performed pyrosequencing assays in duplicate for all specimens whose values fell more than two SD above or below the means, in which case we used the average of the two runs. The values obtained represent the mean methylation for the CpG sites contained within the sequence being analyzed (Additional file 1: Figure S1).
We calculated frequencies and means of sociodemographic variables and conducted multiple linear regression to test the association between DNA methylation and early anthropometric outcomes. We determined covariates a priori based on directed acyclic graphs (DAG). We chose sex as a potential effect measure modifier (EMM), as DNA methylation has been previously shown to vary by sex [44, 45]. We tested the following covariates as potential confounders: maternal education (less than a college degree/college degree or greater), maternal gestational diabetes (yes/no), maternal pre-pregnancy BMI, maternal smoking at any time during pregnancy (yes/no), gestational weight gain, parity (primiparous, multiparous), maternal age at delivery, gestational age, paternal race, maternal race, and date of length and weight measurements relative to child’s birthday. We tested potential confounders in the model one at a time and kept variables if they changed the estimate by more than 10%. Final models included maternal race, maternal education, maternal pre-pregnancy BMI, maternal smoking, maternal age, gestational age, and paternal race. As infants who are SGA may have different growth patterns compared to infants who are average for gestational age, we conducted supplemental analysis to determine the effect of excluding infants who were SGA. We also conducted supplemental analysis including maternal alcohol consumption during pregnancy as a covariate, and an additional supplemental analysis, in which we stratified by race/ethnicity to determine possible effect measure modification.
Previously reported Cronbach’s alpha for correlations among methylation values from all CpGs measured at each DMR was > 0.89 ; therefore, we used mean DNA methylation values for each DMR. DNA methylation was assessed in tertiles (low, moderate, high), as both higher and lower levels of methylation have been associated with health outcomes, depending on the DMR [39, 40]. Given the expected 50% methylation of imprinted genes, we used the mid tertile of methylation as the referent category. Thus, results represent the child z scores associated with high or low methylation compared to “moderate” methylation. We conducted all statistical analysis using SAS 9.4 (SAS Institute, Inc., Cary, NC).
Among infants, 37.0% of mothers were African American, 28.3% were White/Caucasian, and 34.7% were of other races/ethnicities including Hispanic and Asian/Pacific Islander (Table 1). For the 1-year sample, 37.9% of mothers were African American, 30.2% were White, and 31.9% were “other” race. For the 3-year sample, 39.4% of mothers were African American, 28.3% were White, and 32.3% were of other races and ethnicities. In all samples, the majority of women in the study completed less than a college degree (70.8% for newborns, 66.4% for age 1 year, and 68.4% for age 3 years) and reported not smoking at any point during pregnancy (83, 85.5, and 85.3% for newborns, age 1, and age 3, respectively). Approximately half of the newborn sample reported some sort of alcohol consumption in early pregnancy (50.9%). The mean (SD) maternal age for women in the birth sample was 28 years (± 5.7). Mothers in the 1-year sample were on average 28.0 (± 5.8) years, and those in the 3-year sample were on average 28.1 (± 5.8) years. The mean maternal pre-pregnancy BMI for women in the birth weight sample was 27.4 (± 7.2), BMI for mothers in the 1-year sample was 28.0, and BMI for mothers in the 3-year sample was 28.1 (± 5.8) years. The mean gestational age for the sample at birth was 38.7 (1.7) weeks. The mean birthweight of infants in the sample was 3304 g (± 540). There were no significant differences in the study sample demographic makeup between children or mothers in the newborn, age 1 year, or age 3 year samples (data not shown). However, women were more likely to be college educated in the sample of 166 complete cases, in which children had anthropometric data for all 3 time points.
Birth weight and DNA methylation by sex
In girls, we observed a statistically significant association between high methylation at MEST and greater birth weight for gestational age (β = 0.45; 95% CI 0.12, 0.78; p value 0.007; data not shown). However, this association did not persist after adjustment. We observed no statistically significant associations between methylation and birth weight for gestational age z scores in boys.
Weight-for-length z scores at 1 year and DNA methylation by sex
After adjustment (Table 2), we observed an association between high PEG10 DMR methylation and greater WFL z scores at 1 year in boys (β = 0.84; 95% CI 0.34, 1.33; p value = 0.001). Alternatively, low methylation at IGF2 DMR was associated with a lower WFL z score at 1 year in boys (β = − 0.63; 95% CI − 1.16, − 0.10; p value = 0.02). In girls, both low and high methylation at the PLAGL1 DMR (low: β = 0.72; 95% CI − 1.19, − 0.25; p value = 0.003; high: β = − 0.81; 95% CI − 1.29, − 0.33; p value = 0.0001) and the MEST DMR (low: β = − 0.99; 95% CI − 1.59, − 0.39; p value = 0.002; high: β = − 1.03; 95% CI − 1.60, − 0.45; p value = 0.001) were associated with lower WFL z scores at age 1 year after adjustment (Fig. 1).
BMI z scores at 3 years and DNA methylation by sex
At age 3 years, the association between high methylation at the PEG10 DMR and greater anthropometric z score persisted after adjustment in boys (β = 1.03; 95% CI 0.37, 1.69; p value 0.003). In girls, the association between high methylation at the MEST DMR and lower anthropometric z score also persisted after adjustment (β = − 1.11; 95%CI − 1.98, − 0.24; p value 0.01). In addition, we observed an association between low methylation at the PEG3 DMR and lower BMI z score at 3 years after adjustment in girls (β = − 0.95; 95% CI − 1.80, − 0.10; p value = 0.03).
Analysis excluding SGA infants
When excluding SGA infants (n = 69 at birth, n = 31 at 1 year, n = 27 at 3 years), additional associations emerged at the MEST DMR at birth and the MEG3, H19, and NNAT DMRs at 1 year of age (data not shown). High methylation at the MEST DMR was associated with a greater BW/GA z score in girls (β = 0.32; 95% CI 0.009, 0.63; p value = 0.04). At 1 year, high and low methylation at the MEG3 DMR were associated with a lower WFL z score in boys (low: β = − 0.81; 95% CI − 1.44, − 0.18; p value = 0.01; high: β = − 0.91; 95% CI − 1.57, − 0.25; p value = 0.008). In girls, high methylation at the H19 DMR was also associated with a lower WFL z score (β = − 0.58; 95% CI − 1.15, − 0.008; p value = 0.047), and lower methylation at the NNAT DMR was associated with a greater WFL z score (β = 0.63; 95% CI 0.13, 1.13; p value = 0.02). All other associations, with the exception of the association between low methylation at the IGF2 DMR and a lower WFL z score at 1 year, persisted after exclusion of SGA infants.
Analysis including maternal alcohol consumption as a covariate
We conducted additional analyses, in which we included maternal alcohol consumption during early pregnancy as a covariate. We found that all associations between methylation and WFL z scores at 1 year remained statistically significant compared to our main analysis (data not shown). Additionally, we observed an association between low IGF2 methylation and greater birth weight for gestational age z score in girls (β = 0.53; 95% CI 0.13, 0.93; p value = 0.01) and an association between both low and high MEST DMR methylation and lower BMI z scores at 3 years in girls (low: β = − 1.31; 95% CI − 2.46, − 0.15; p value = 0.03; high: β = − 1.40; 95% CI − 2.51, − 0.29; p value = 0.01).
Analysis stratified by race/ethnicity
Additionally, we conducted supplemental stratified analyses to see if there were differences by race/ethnicity, as previous studies have found differential methylation in association with race/ethnicity. The only statistically significant association we observed was between low methylation at the MEST DMR and greater birth weight for gestational z scores among Blacks (β = 0.54; 95% CI 0.09, 0.99; p value = 0.02—data not shown). At year 1, we found an association between high PEG10 DMR methylation and greater WFL z scores among Blacks and Whites (Blacks: β = 0.69; 95% CI 0.12, 1.27; p value = 0.02; Whites: β = 0.66; 95% CI 0.11, 1.21; p value = 0.02). We also found an association between high MEG3-IG DMR methylation and lower WFL z scores among Whites (β = − 0.64; 95% CI − 1.25, − 0.03; p value = 0.04). At 3 years, we found an association between low MEG3-IG DMR methylation and lower BMI z scores among Blacks (β = − 0.76; 95% CI − 1.48, − 0.05; p value = 0.04). An association was also observed between low NNAT DMR methylation and greater BMI z scores among Blacks (β = 0.99; 95% CI 0.14, 1.84; p value = 0.02).
Complete case analysis
Supplemental analysis on the 166 infants (n = 15 SGA children) who had non-missing anthropometric data for all 3 time points showed directional consistency in all associations. However, not all associations remained statistically significant in the smaller sample (Table 3). In girls, the association between high methylation at the PLAGL1 DMR and lower WFL z scores remained statistically significant (β = − 0.82; 95% CI − 1.51, − 0.12; p value = 0.02), as did the associations between low and high methylation at the MEST DMR and WFL z scores (low: β = − 1.18; 95% CI − 2.04, − 32; p value = 0.009; high: β = − 1.57; 95% CI − 2.4, − 0.73; p value = 0.0004). At 3 years in boys, the association between high PEG10 DMR methylation and greater BMI z scores remained statistically significant (β = 1.25; 95% CI 0.38, 2.12; p value = 0.006). In girls, the association between lower PEG3 DMR methylation and lower BMI z scores at 3 years also remained statistically significant (β = − 1.30; 95% CI − 2.34, 0.26; p value = 0.02). Notably, additional statistically significant associations were observed between high methylation at the MEG3-IG DMR and lower birth weight for gestational age z scores in boys (β = − 0.77; 95% CI − 1.30, − 0.24; p value = 0.005) and between high methylation at the IGF2 and PEG10 DMRs and greater birth weight for gestational age z scores in girls (IGF2: β = 0.58; 95% CI 0.03, 1.12; p value = 0.04; PEG10: β = 0.63; 95% CI 0.07, 1.20; p value = 0.03).
In these analyses, we examined DNA methylation of nine regulatory regions at birth and anthropometric measures at birth and age 1 and 3 years. No DMR showed a consistent association between methylation and anthropometric z scores at all 3 time points explored. Our key findings were that high methylation of the sequences regulating the PEG10 DMR at birth was associated with a higher age 1-year WFL z score and 3-year BMI z score in boys, low methylation at the NNAT DMR was associated with higher BMI z scores at age 3 in girls, and high methylation was associated with lower WFL z scores at 1 year and BMI z scores at 3 years. These associations persisted after excluding SGA infants. Additional findings included associations between methylation at the PLAGL1 DMR and WFL z scores in girls at age 1, as well as an association at the IGF2 DMR at age 1 among boys. At age 3, we also observed an association between methylation at the PEG3 DMR and BMI z scores in girls. These results suggest that methylation of imprinted genes at birth is associated with anthropometric measures at ages 1 and 3 years, with the PEG10 and NNAT DMRs potentially indicating an early risk for obesity, and MEST potentially indicating a lower risk for obesity.
This study adds to a growing body of epidemiologic evidence on early postnatal growth associated with DNA methylation at birth and suggests that DNA methylation at multiple DMRs may be associated with WFL z score at age 1 year and BMI z score at age 3 years. PLAGL1, MEST, NNAT, and PEG10 have been associated with obesity or weight in previous literature. Paternal obesity has been previously associated with reduced PEG10 transcription in mouse placentas . Our results indicating a higher level of PEG10 methylation is associated with a greater BMI z score show a similar pattern. Previous literature has also shown a possible association between increased MEST expression and inhibition of adipogenesis . The results of this study echo these findings, as boys with higher than average methylation levels had lower BMI z scores. Methylation at the MEST, NNAT, and PEG10 DMRs has also been previously associated with paternal obesity , and small and large for gestational age [48, 49]. In addition, the NNAT gene has been associated with severe obesity in childhood and adulthood . Higher methylation at PLAGL1 has also been associated with maternal obesity , and fetal and postnatal growth . A previous study found a positive correlation between PLAGL1 methylation and BMI z scores at age 1 year. PLAGL1 is thought to be an imprint control region ; however, the implications of this in relation to a potential role in the risk of early obesity are not yet clear. More research is needed to gain a better understanding of the relationship between child BMI and methylation at these DMRs. Previous literature has supported the role of the IGF2 DMR in fetal growth [31,32,33,34] and birth weight [32, 52], as well as infant weight gain [32, 36] and child adiposity . We observed an association between IGF2 methylation and lower BMI z scores; however, this association was not seen at birth or age 3 years.
We observed sex-specific differences in the associations between DNA methylation at birth and anthropometric z scores at ages 1 and 3 years. Sex-specific methylation has been previously observed in relation to nutrition and other environmental exposures [44, 45, 53], as well as in relation to outcomes, such as small for gestational age . However, these studies did not find sex-specific differences in methylation in PEG10, NNAT, or MEST. This study also adds to the growing literature on sex-specific DNA methylation.
We conducted additional analysis to explore the influence of infants who are SGA, as their growth patterns may differ from those of infants who are not SGA, and found that SGA may be associated with DNA methylation. No associations between methylation and BW/GA z scores remained significant after SGA exclusion, suggesting that SGA may have been driving these associations. We also found that associations between PEG3 and BMI z scores at both ages 1 and 3 years became statistically significant after exclusion of SGA infants. This suggests that perhaps the association between PEG3 methylation and SGA is in the opposite direction of the association between PEG3 methylation and BMI z scores for non-SGA infants, thus attenuating the original association. However, a cautious interpretation is warranted, as exclusion of SGA infants also decreased the statistical power, which may have led to unstable estimates. Additional supplemental analysis including only the 166 infants who had non-missing anthropometric values at all 3 time points showed that many of our associations remained statistically significant, including our key findings at the PEG10 and NNAT DMRs. This suggests that these associations were not related to differences in the samples at each time point. Additional associations emerged as statistically significant at MEST, IGF2, PEG10, and NNAT; however, these results must be interpreted with caution, as the analysis was underpowered. Similarly, the results of the supplemental analysis including maternal alcohol consumption and the analysis stratified by race/ethnicity should be interpreted with caution, as our sample size was greatly reduced, and estimates may be underpowered.
The direction of the association between higher methylation at MEST and z scores changed from positive to negative from birth to ages 1 and 3. The reasons for this are unclear; however, it is possible that the modest increase in BW/GA z scores at birth, which we found to be associated with a higher level of methylation at MEST, is also associated with a greater likelihood of becoming lean as the child grows and becomes more mobile. MEST expression has been previously associated with obesity in mice ; however, a study in humans found MEST to possibly inhibit adipogenesis .
This study benefits from an ethnically diverse cohort, and prospectively collected data at multiple time points. This facilitates a better understanding of the timing of methylation with regard to our outcome of interest, weight gain. In addition, the use of BMI z scores provides widely accepted estimate of adiposity that accounts for a child’s age. However, it is not without limitations. Our study only included nine DMRs of imprinted genes. Although these genes were chosen because they have been linked to growth or chronic disease, it is possible that important genes have been left out of this study. In addition, our study’s small sample size may have limited our ability to see statistically significant differences associated with DNA methylation among our population. Although our analyses were hypothesis driven, multiple testing is a limitation in this study, as it may increase the possibility that our results are observed by chance. However, many associations remained (PEG10, MEST, PLAGL1, NNAT) even after the stringent Bonferroni correction. A final limitation was the use of weight status instead of weight gain. Much of the literature has pointed to weight gain in the first year of life as being associated with later obesity ; however, there is some literature indicating high weight status in early childhood as a risk factor for later obesity . It is unclear whether or not these findings are related to later obesity.
In addition, the results of this study show differences in anthropometric z scores in association with DNA methylation that is either lower or higher than the “average” methylation of our study sample. Notably, the actual change in continuous methylation associated with our results is likely small (approximately 1%). However, even a 1% change in methylation has been previously shown to result in a doubling or halving of gene expression at these imprint regulatory regions . Not assessing DNA methylation continuously may be a limitation of our study, as it creates a challenge in comparing our results to those of other studies. However, we believe the results presented in this study may be meaningful for public health, as it provides a range of methylation values that may be associated with anthropometric, and possibly even adiposity.
In summary, our study findings suggest that DNA methylation of the PLAGL1, MEST, PEG10, and NNAT DMRs at birth is associated with BMI z scores in early childhood and varies by sex. Longitudinal assessment of DNA methylation in these DMRs at older age time points is needed to determine whether or not methylation at these DMRs is associated with obesity later in life. Determining the associations between DNA methylation and early obesity risk is important, as DNA methylation of regulatory regions may serve as markers for the assessment of early obesity risk. However, gaining a better understanding of the exposures that affect methylation at these regions is also important, as exposures that modify methylation of regions that are associated with obesity risk may be a good target for early obesity prevention efforts.
Body mass index
Birth weight for gestational age
Differentially methylated region
Insulin-like growth factor 2
Maternally expressed gene 3
Maternally expressed gene 3—intergenic
Paternally expressed gene 10
Paternally expressed gene 3
Pleiomorphic adenoma gene-like 1
Small for gestational age
Monasta L, Batty GD, Cattaneo A, et al. Early-life determinants of overweight and obesity: a review of systematic reviews. Obes Rev. 2010;11(10):695–708.
Rooney BL, Mathiason MA, Schauberger CW. Predictors of obesity in childhood, adolescence, and adulthood in a birth cohort. Matern Child Health J. 2011;15(8):1166–75.
Skilton MR, Marks GB, Ayer JG, et al. Weight gain in infancy and vascular risk factors in later childhood. Pediatrics. 2013;131(6):e1821–8.
Freedman DS, Dietz WH, Srinivasan SR, Berenson GS. The relation of overweight to cardiovascular risk factors among children and adolescents: the Bogalusa Heart Study. Pediatrics. 1999;103:1172-85.
Dietz WH. Health consequences of obesity in youth: childhood predictors of adult disease. Pediatrics. 1998;101:518–25.
Yu ZB, Han SP, Zhu GZ, et al. Birth weight and subsequent risk of obesity: a systematic review and meta-analysis. Obes Rev. 2011;12(7):525–42.
Gallo P, Cioffi L, Limauro R, et al. SGA children in pediatric primary care : what is the best choice, large or small ? A 10-year prospective longitudinal study. Global Pediatric Health. 2016;3(194):1–7.
Baidal JAW, Locks LM, Cheng ER, Blake-lamb TL, Perkins ME, Taveras EM. Risk factors for childhood obesity in the first 1,000 days. Am J Prev Med. 2016;50(6):761–79.
Vidal AC, Overcash F, Murphy SK, et al. Associations between birth and one year anthropometric measurements and IGF2 and IGF2R genetic variants in African American and Caucasian American infants. Journal of Pediatric Genetics. 2013;2:119–27.
Dietz WH. Critical periods in childhood for the development of obesity. Am J Clin Nutr. 1994;1:1–5.
Barker DJP. The fetal and infant origins of adult disease the womb may be more important than the home. BMJ. 1990;301(156):1990–0.
Gluckman PD, Hanson MA, Cooper C, Thornburg KL. Effect of in utero and early-life conditions on adult health and disease. N Engl J Med. 2008;359(1):61–73.
Langley-Evans SC. Developmental programming of health and disease. Proc Nutr Soc. 2006;65(1):97–105.
Dolinoy DC, Weidman JR, Ra W, Jirtle RL. Maternal genistein alters coat color and protects Avy mouse offspring from obesity by modifying the fetal epigenome. Environ Health Perspect. 2006;114(4):567–72.
Godfrey KM, Sheppard A, Gluckman PD, et al. Epigenetic gene promoter methylation at birth is associated with child’s later adiposity. Diabetes. 2011;60(5):1528–34.
Plagemann A, Roepke K, Harder T, et al. Epigenetic malprogramming of the insulin receptor promoter due to developmental overfeeding. J Perinat Med. 2010;38(4):393–400.
Heijmans BT, Tobi EW, Lumey LH, Slagboom PE. The epigenome: archive of the prenatal environment. Epigenetics. 2009;4(8):526–31.
Ja H, Surani MA. DNA methylation dynamics during the mammalian life cycle. Philosophical transactions of the Royal Society of London Series B, Biological sciences. 2013;368(1609):20110328.
Reik W, Walter J. Genomic imprinting: parental influence on the genome. Nature Reviews. 2001;2(January):21–32.
Gluckman PD, Hanson MA, Low FM. The role of developmental plasticity and epigenetics in human health. Birth Defects Research Part C. 2011;93(1):12–8.
Dolinoy DC, Weidman JR, Jirtle RL. Epigenetic gene regulation: linking early developmental environment to adult disease. Reprod Toxicol. 2007;23:297–307.
Heijmans BT, Tobi EW, Stein AD, et al. Persistent epigenetic differences associated with prenatal exposure to famine in humans. Proc Natl Acad Sci U S A. 2008;105:17046–9.
Tobi EW, Goeman JJ, Monajemi R, et al. DNA methylation signatures link prenatal famine exposure to growth and metabolism. Nat Commun. 2014;5:5592.
Dominguez-Salas P, Moore SE, Baker MS, et al. Maternal nutrition at conception modulates DNA methylation of human metastable epialleles. Nat Commun. 2014;5:1–7.
Cruz-Correa M, Zhao R, Oviedo M, et al. Temporal stability and age-related prevalence of loss of imprinting of the insulin-like growth factor-2 gene. Epigenetics. 2009;4(2):114–8.
Murphy SK, Huang Z, Hoyo C. Differentially methylated regions of imprinted genes in prenatal, perinatal and postnatal human tissues. PLoS One. 2012;7(7):e40924.
Hoovers JM, Kalikin LM, Johnson LA, et al. Multiple genetic loci within 11p15 defined by Beckwith-Wiedmann syndrome rearrangement breakpoints and subchromosomal. Proc Natl Acad Sci USA. 1995;92:12456–60.
Bliek J, Terhal P, van den Bogaard M-J, et al. Hypomethylation of the H19 gene causes not only Silver-Russell syndrome (SRS) but also isolated asymmetry or an SRS-like phenotype. Am J Hum Genet. 2006;78(4):604–14.
Pan H, Wu Y, Chen L, et al. HIF3A association with adiposity: the story begins before birth. Epigenomics. 2015;7(6):937–50.
Huang RC, Garratt ES, Pan H, et al. Genome-wide methylation analysis identifies differentially methylated CpG loci associated with severe obesity in childhood. Epigenetics. 2015;10(11):995–1005.
Kadakia R, Josefson J. The relationship of insulin-like growth factor 2 to fetal growth and adiposity. Hormone Res Paediatr. 2016;85:75–82.
Bouwland-Both MI, van Mil NH, Stolk L, et al. DNA methylation of IGF2DMR and H19 is associated with fetal and infant growth: the generation R study. PLoS One. 2013;8(12):e81731.
Su R, Wang C, Feng H, et al. Alteration in expression and methylation of IGF2/H19 in placenta and umbilical cord blood are associated with macrosomia exposed to intrauterine hyperglycemia. PLoS One. 2016;11(2):e0148399.
St-pierre J, Hivert M-f, Perron P, et al. IGF2 DNA methylation is a modulator of newborn ’ s fetal growth and development. Epigenetics. 2012;7(10):1125–32.
Huang R-C, Galati JC, Burrows S, et al. DNA methylation of the IGF2/H19 imprinting control region and adiposity distribution in young adults. Clin Epigenetics. 2012;4(1):21.
Perkins E, Murphy SK, Murtha AP, et al. Insulin-like growth factor 2/H19 methylation at birth and risk of overweight and obesity in children. J Pediatr. 2012;161(1):31–9.
Azzi S, Sas TCJ, Koudou Y, et al. Degree of methylation of ZAC1 (PLAGL1) is associated with prenatal and post-natal growth in healthy infants of the EDEN mother child cohort. Epigenetics. 2014;9(3):338–45.
Kameswaran V, Bramswig NC, McKenna LB, et al. Epigenetic regulation of the DLK1-MEG3 microRNA cluster in human type 2 diabetic islets. Cell Metab. 2014;19(1):135–45.
Soubry A, Murphy SK, Wang F, et al. Newborns of obese parents have altered DNA methylation patterns at imprinted genes. Int J Obes. 2015;39(4):650–7.
Liu Y, Murphy SK, Murtha AP, et al. Depression in pregnancy, infant birth weight and DNA methylation of imprint regulatory elements. Epigenetics. 2012;7(7):735–46.
Villar J, Ismail LC, Victora CG, et al. International standards for newborn weight , length , and head circumference by gestational age and sex : the Newborn Cross-Sectional Study of the INTERGROWTH-21 st Project. Lancet. 2014;384:857–68.
World Health Organization. The WHO Growth Standards. http://www.who.int/childgrowth/standards/en/. Accessed Sept 10, 2016.
Centers for Disease Control and Prevention. A SAS program for the 2000 CDC Growth Charts (ages 0 to <20 years). 2016; http://www.cdc.gov/nccdphp/dnpao/growthcharts/resources/sas.htm. Accessed Sept 10, 2016.
Tobi EW, Lumey LH, Talens RP, et al. DNA methylation differences after exposure to prenatal famine are common and timing- and sex-specific. Hum Mol Genet. 2009;18(21):4046–53.
Murphy SK, Adigun A, Huang Z, et al. Gender-specific methylation differences in relation to prenatal exposure to cigarette smoke. Gene. 2012;494(1):36–43.
Mitchell M, Strick R, Strissel PL, et al. Gene expression and epigenetic aberrations in F1-placentas fathered by obese males. Mol Reprod Dev. 2017;84(November 2016):316–28.
Karbiener M, Glantschnig C, Pisani DF, et al. Mesoderm-specific transcript (MEST) is a negative regulator of human adipocyte differentiation. Int J Obes. 2015;39(12):1733–41.
Kappil MA, Green BB, Armstrong DA, et al. Placental expression profile of imprinted genes impacts birth weight. Epigenetics. 2015;10(9):842–9.
Lim AL, Ng S, Ching S, et al. Epigenetic state and expression of imprinted genes in umbilical cord correlates with growth parameters in human pregnancy. J Med Genet. 2012;49:689–97.
Vrang N, Meyre D, Froguel P, et al. The imprinted gene neuronatin is regulated by metabolic status and associated with obesity. Obesity (Silver Spring, Md). 2010;18(7):1289–96.
Arima T, Ra D, Arney KL, et al. A conserved imprinting control region at the HYMAI/ZAC domain is implicated in transient neonatal diabetes mellitus. Hum Mol Genet. 2001;10(14):1475–83.
Ong KK, Loos RJF. Rapid infancy weight gain and subsequent obesity: systematic reviews and hopeful suggestions. Acta Paediatr. 2006;95(8):904–8.
Gonzalez-Nahm S, Mendez MA, Robinson WR, et al. Low maternal adherence to a Mediterranean diet is associated with increase in methylation at the MEG3-IG differentially methylated region in female infants. Environmental Epigenetics. 2017;3(March):1–10.
Qian Y, Huang X, Liang H, et al. Effects of maternal folic acid supplementation on gene methylation and being small for gestational age. J Hum Nutr Diet. 2016;29:643–51.
Voigt A, Ribot J, Sabater G, Palou A, Bonet ML, Klaus S. Identification of Mest/Peg1 gene expression as a predictive biomarker of adipose tissue expansion sensitive to dietary anti-obesity interventions. Genes Nutr. 2015;10(5):27.
Evensen E, Emaus N, Kokkvoll A, Wilsgaard T, Furberg AS, Skeie G. The relation between birthweight, childhood body mass index, and overweight and obesity in late adolescence: a longitudinal cohort study from Norway, the Tromso study, Fit Futures. BMJ Open. 2017;7(6):e015576.
We thank the NEST participants for their generous contributions to the study and Carole Grenier for her excellent technical assistance.
This research was supported in part by the National Institutes of Health, grant numbers R01ES016772, P30ES025128, R01DK094841, and P01ES022831, and USEPA grant RD-83543701.
Availability of data and materials
Please contact the authors for data requests.
Ethics approval and consent to participate
This study was approved by the Institutional Review Board at Duke University Medical Center. Informed consent was obtained from all participants.
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.
Figure S1. Pyrograms from bisulfite pyrosequencing of representative samples. Example results from three individual cord blood specimens are shown for each DMR. The bisulfite-modified version of the sequence to analyze is shown at the top of each pyrogram, and the actual sequencing output, by base, is shown at the bottom. Pyrosequencing is a sequence-by-synthesis method, and the light generated by the incorporation of each nucleotide (y axis) is proportional to the amount of the template present in the reaction. For mononucleotides in the template sequence, a single peak height is generated, but for runs of two or more of the same nucleotide in the sequence, the peak height is proportional to that number. Blue diamonds at the top of peaks are the expected peak height based on the sequence. Orange diamonds in yellow-brown vertical bars are bisulfite controls, where the original sequence contains a non-CpG cytosine that should be fully converted by the bisulfite, and thus, no cytosine signal should be detected at these positions. There are no diamonds above the potentially methylated cytosine positions (CGs), which are represented by a “T” position (for unmethylated cytosine) followed by a “C” position (for methylated cytosines). These positions are marked by the grey vertical bars, above which shows the percent methylation for that position of the sequence. (PDF 350 kb)
About this article
Cite this article
Gonzalez-Nahm, S., Mendez, M.A., Benjamin-Neelon, S.E. et al. DNA methylation of imprinted genes at birth is associated with child weight status at birth, 1 year, and 3 years. Clin Epigenet 10, 90 (2018). https://doi.org/10.1186/s13148-018-0521-0