OBJECTIVE

Fixed genomic variation explains only a small proportion of the risk of adiposity. In animal models, maternal diet alters offspring body composition, accompanied by epigenetic changes in metabolic control genes. Little is known about whether such processes operate in humans.

RESEARCH DESIGN AND METHODS

Using Sequenom MassARRAY we measured the methylation status of 68 CpGs 5′ from five candidate genes in umbilical cord tissue DNA from healthy neonates. Methylation varied greatly at particular CpGs: for 31 CpGs with median methylation ≥5% and a 5–95% range ≥10%, we related methylation status to maternal pregnancy diet and to child’s adiposity at age 9 years. Replication was sought in a second independent cohort.

RESULTS

In cohort 1, retinoid X receptor-α (RXRA) chr9:136355885+ and endothelial nitric oxide synthase (eNOS) chr7:150315553+ methylation had independent associations with sex-adjusted childhood fat mass (exponentiated regression coefficient [β] 17% per SD change in methylation [95% CI 4–31], P = 0.009, n = 64, and β = 20% [9–32], P < 0.001, n = 66, respectively) and %fat mass (β = 10% [1–19], P = 0.023, n = 64 and β =12% [4–20], P = 0.002, n = 66, respectively). Regression analyses including sex and neonatal epigenetic marks explained >25% of the variance in childhood adiposity. Higher methylation of RXRA chr9:136355885+, but not of eNOS chr7:150315553+, was associated with lower maternal carbohydrate intake in early pregnancy, previously linked with higher neonatal adiposity in this population. In cohort 2, cord eNOS chr7:150315553+ methylation showed no association with adiposity, but RXRA chr9:136355885+ methylation showed similar associations with fat mass and %fat mass (β = 6% [2–10] and β = 4% [1–7], respectively, both P = 0.002, n = 239).

CONCLUSIONS

Our findings suggest a substantial component of metabolic disease risk has a prenatal developmental basis. Perinatal epigenetic analysis may have utility in identifying individual vulnerability to later obesity and metabolic disease.

Fixed genomic variations explain only a fraction of the risk of human obesity and metabolic disease (1,2). However, there is increasing epidemiological evidence linking perinatal factors to later adiposity and metabolic disease risk (3,4). For example, famine during pregnancy is associated with obesity in the adult offspring (5) and normal variations in maternal body composition relate to child’s later adiposity (6). Although understanding of the underlying mechanisms is limited, data from animal models suggest that epigenetic processes are an important link between the early life environment, for example maternal diet, and altered metabolism and body composition in the adult offspring (79).

Epigenetic processes such as DNA methylation and histone modifications allow the developmental environment to modulate gene transcription; many of these changes are then stable throughout the lifecourse (10,11). Regulated DNA methylation mostly occurs at a cytosine immediately 5′ to a guanine (CpG sites). Such processes are involved not only in cell differentiation and parental genomic imprinting but also in developmental plasticity through which the environment in early life can affect the developmental trajectory, with long-term effects on gene expression and phenotypic outcome (1214). For example, in the rat unbalanced maternal diet during pregnancy induces changes in DNA methylation and covalent histone modifications in the 5′ regulatory regions of specific nonimprinted genes (1517) and affects the offspring’s later body composition and metabolic phenotype. Induced changes in phenotype can be prevented by nutritional interventions during pregnancy (18), or altered by nutritional interventions during the juvenile-pubertal period (19) or by hormonal interventions during suckling (20).

Although epigenetic processes operating in early development have been implicated in the origins of obesity (3,11), there is as yet no direct evidence for this proposition in humans. Furthermore, there has been considerable debate as to the relative importance of such developmental pathways in determining phenotypic outcomes. Any association between epigenetic state at birth and later phenotype would allow an estimate of the possible developmental contribution, irrespective of whether such an association was causal or simply reflected the developmental state.

We reasoned that the targeted measurement of DNA methylation in human fetal tissues at birth might not only provide evidence that environmental influences have affected prenatal development but, if they then correlated with later phenotype, would provide an approach to demonstrating the role of the prenatal environment in predisposition to adiposity.

Here we first measured the methylation status of CpGs in the promoters of candidate genes in DNA extracted from umbilical cord tissue obtained at birth in children who were later assessed for adiposity at age 9 years (the Princess Anne Hospital [PAH] study) (21). Measurements of perinatal DNA methylation were related to adiposity in later childhood and to information on mother’s diet during pregnancy. Because of the strong associations found, we then sought to replicate associations between umbilical cord CpG methylation and child’s adiposity in a second independent group of children from the Southampton Women’s Survey (SWS) (22).

We studied two prospective cohorts recruited antenatally in Southampton, U.K. In the PAH study, Caucasian women ≥16 years old with singleton pregnancies <17 weeks’ gestation were recruited and a validated food frequency questionnaire (23) was administered at 15 weeks’ gestation; patients with diabetes and hormonally induced conceptions were excluded. When the children approached 9 years, we wrote to 461 still living locally. Two hundred and sixteen (47%) attended a clinic, and adiposity was measured using dual energy X-ray absorptiometry (DEXA); 78 had DNA available from an umbilical cord sample collected at birth and stored at −80°C. None of the mothers were administered glucocorticoids during pregnancy.

In the SWS, women aged 20–34 years were recruited when not pregnant; those who subsequently conceived were followed through pregnancy and the offspring followed up (22). To replicate associations with adiposity found in the PAH cohort, we studied 239 children selected as having both umbilical cord DNA available and childhood adiposity measurements.

Both groups underwent adiposity measurement by DEXA (PAH study, age 9 years: Lunar DPX-L, pediatric software, version 4.7c; GE Corporation, Madison, WI [6]; SWS, age 6 years: Hologic Discovery, pediatric scan mode; Hologic, Bedford, MA). Instruments were calibrated daily; coefficients of variation were <1%. Follow-up of the children and sample collection/analysis was carried out under Institutional Review Board approval (Southampton and SW Hampshire Research Ethics Committee) with written informed consent. Investigations were conducted according to the principles expressed in the Declaration of Helsinki.

Selection of candidate genes.

DNA from 15 randomly selected PAH study umbilical cords underwent methylation-specific chromatin precipitation followed by hybridization to a commercial tiled oligomer microarray (NimbleGen Systems HG17_min_promoter array [24]), which uses 50-mer oligonucleotides positioned on average every 100 bp from −3,750 bp to +500 bp relative to the transcription start site of 24,134 human genes. We focused our analysis initially on a panel of 78 candidate genes, 28 selected on the basis of animal data from our and other laboratories (25,26) and 50 selected from the array as having the highest SD between subjects relative to the SD of within-subject replicates. For this panel we identified those with correlations between overall gene methylation status and DEXA measurements of adiposity at age 9 years and from these genes chose five genes for further study based on individual oligomers showing evidence of correlation with DEXA measurements at age 9 years, biological plausibility and feasibility of designing amplicons suitable for Sequenom analysis; genes selected were retinoid X receptor-α (RXRA), endothelial nitric oxide synthase (eNOS), superoxide dismutase-1 (SOD1), interleukin-8 (IL8), and phosphoinositide-3-kinase, catalytic, δ-polypeptide (PI3KCD). These genes all had associations of comparable strength between childhood body composition and both overall gene promoter methylation and the methylation of individual oligomers. All chromosomal coordinates are based on University of California Santa Cruz, human genome March 2006 assembly (hg18).

Sequenom and single nucleotide polymorphism analysis.

Genomic DNA was purified from frozen cords by proteinase K digestion and phenol:chloroform extraction. DNA methylation was measured with the Sequenom MassARRAY Compact System (www.sequenom.com/). Briefly, gene-specific amplification of bisulfite-treated DNA was followed by in vitro transcription and analysis by matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry (27,28). Sequenom assay design and methods were as per the manufacturer’s protocol. DNA (1 μg) was bisulfite converted using EZ DNA Methylation kit (Zymo Research). PCR primers specific for bisulfite converted DNA were designed (Methprimer). Each reverse primer contained a T7-promoter tag for in vitro transcription (5′-cagtaatacgactcactatagggagaaggct-3′), and the forward primer was tagged with a 10mer to balance Tm (5′-aggaagagag-3′). Supplementary Table 1 lists amplicons, primer sequences, PCR annealing temperature, and genomic coordinates for the extended promoter regions measured. Bisulfite-treated DNA was PCR amplified (Qiagen HotStar Taq Polymerase) in 5 μL reactions and treated with Shrimp Alkaline Phosphatase (Sequenom) (20 min at 37°C), heat inactivated (85°C for 5 min), and simultaneous in vitro transcription/uracil-cleavage reaction was carried out (7 μL reactions) using Sequenom T-cleavage reagent mix. Transcription cleavage products were desalted with 6 mg of Clean-Resin and 20 nL spotted on a 384-pad SpectroCHIP (Sequenom) using a MassARRAY nanodispenser (Samsung). Mass spectra were acquired using a MassARRAY MALDI-TOF MS (Bruker-Sequenom) and peak detection, signal-to-noise calculations, and quantitative CpG site methylation performed using proprietary EpiTyper software v1.0 (Sequenom). Consistent with the manufacturer’s specification (www.sequenom.com) agreement between technical replicates was good, duplicate measures differing by <10% methylation in almost all cases. We excluded from analysis samples that failed to give a reliable PCR product or produced spectra with low confidence scores (<2.9 in EpiTyper). For fragments containing a single CpG site, DNA methylation was calculated by the ratio of methylated to unmethylated fragments. Limitations imposed by Sequenom analysis treat cleavage products containing multiple CpG sites as single units, and the methylation values reported are weighted averages across the unit (referred to as a CpG group). DNA quality and no-template controls, 0%, and 100% methylated DNA were included in all assays. Supplementary Table 2 shows the distributions of methylation values for CpGs and CpG groups and numbers of subjects with measurements available.

Pyrosequencing analysis (PyroMark Q96MD; Qiagen) was undertaken to exclude single nucleotide polymorphisms (SNPs) involving CpG dinucleotides, which would result in loss of potential for methylation. Briefly, DNA was PCR amplified, and the product annealed to streptavidin-coated sepharose beads and denatured (0.2M NaOH) to a single-stranded product. Nucleotides were incorporated to the open 3′ DNA strand in which pyrophosphate is released and used in a sulfurylase reaction emitting ATP, measured in a luciferase reaction, and analyzed using PyroMark MD v1.0 software (Qiagen) (29).

Statistical analysis.

Using Stata 11 (StataCorp) study variables were transformed using Fisher-Yates (30) (for methylation measurements) or logarithmic transformations where necessary to satisfy statistical assumptions of normality. A priori we restricted the analysis to CpGs with median methylation ≥5% and a 5–95% range ≥10%. In PAH subjects we first used Pearson correlation (rp) and linear regression to examine CpG methylation in relation to child’s adiposity. Mutually adjusted regression models were built including CpG sites significantly associated with each outcome in univariate analyses (P < 0.05). Because the measurements of childhood adiposity are on the log scale and the CpG measurements are z scores, the exponentiated results are interpreted as percent change in childhood adiposity measurements per SD change in methylation. To allow assessment of the influence of taking account of the child’s sex, exponentiated regression coefficients are first presented excluding and including adjustment for child's sex. We then preadjusted all adiposity measurements for sex and took account of the mother’s age, adiposity (as continuous variables), and smoking during pregnancy (as a binary variable, yes/no) by including them as covariates in the regression model; we next related mother’s diet to CpG methylation at birth. Finally, we sought replication of the associations between CpG methylation and child’s adiposity in the SWS cohort. SWS measurements, taken over a wider age range than those in PAH subjects, were preadjusted for age and sex.

Median (interquartile range) birth weight, age, fat mass, and %body fat at follow-up values for the PAH subjects were 3,330 g (3,010–3,790), 8.66 years (8.52–8.73), 5.41 kg (3.97–8.41), and 18.3% (14.4–27.0), respectively. Comparable values for the 239 SWS children were 3,485 g (3,205–3,773), 6.59 years (6.41–6.78), 4.60 kg (3.66–5.69), and 23.4% (19.7–27.6), respectively. Similar percentages of mothers smoked in the PAH and SWS cohorts (21 vs. 24%, respectively); median maternal age and prepregnancy BMI were lower in the PAH mothers (28 vs. 31 years and 22.3 vs. 24.3 kg/m2, respectively). The above characteristics are similar to the overall values for participants in the two cohorts.

PAH cohort methylation studies.

Supplementary Table 2 shows the distribution of percent methylation of the 68 CpGs and CpG groups among the 78 PAH subjects; 31 CpGs/CpG groups met a priori criteria for further analysis. Particularly marked interindividual variation in percent methylation was seen at some sites (for example, RXRA chr9:136355885+ had a median of 59% with 5–95th centiles 4–99%; eNOS chr7:150315553+ had a median of 93%, with 5–95th centiles 64–100%). Correlations between methylation of different CpGs were generally low (for example rp= −0.03, P = 0.81, n = 55 between the above two CpGs).

Of the 31 CpGs with variable methylation above the a priori threshold, seven had significant associations with the child’s adiposity at age 9 years (Supplementary Table 3). Table 1 shows the results of uni- and multivariate analyses of CpG sites with significant independent associations with the child’s adiposity. Comparison of the regression coefficients for the univariate and multivariate associations of RXRA chr9:136355885+ and eNOS chr7:150315553+ methylation with child’s fat mass, %fat, and ratio of trunk to limb fat shows that mutual adjustment for each other and sex has little effect on the coefficients; adiposity measurements were therefore preadjusted for sex in further analyses. RXRA chr9:136355885+ methylation had positive associations with childhood fat mass, %fat mass (scatterplots shown in Supplementary Fig. 1), and ratio of trunk to limb fat, with similar positive associations for eNOS chr7:150315553+: exponentiated regression coefficients (β) (95% CI) per SD change in methylation were 17% (4–31), P = 0.009; 10% (1–19), P = 0.023; and 6% (0–12), P = 0.039, respectively (n = 64). Comparable values for eNOS chr7:150315553+ methylation were 20% (9–32), P < 0.001; 12% (4–20), P = 0.002; and 7% (2–13), P = 0.007, respectively (n = 66). Methylation of RXRA chr9:136355885+ and eNOS chr7:150315553+ had similar positive associations with child’s BMI (β = 3% [0–6], P = 0.037; and 4% [2–7], P = 0.001, respectively). Independently of RXRA chr9:136355885+ and eNOS chr7:150315553+ methylation, SOD1 chr21:31853660/63+ methylation was also inversely related to child’s trunk-to-limb fat ratio (P = 0.037). The associations reflect clinically important shifts in adiposity such that the mean fat mass rose from 4.8 kg (17.3% body fat) in the lowest quarter of RXRA chr9:136355885+ methylation to 6.6 kg (21.3% body fat) in the highest quarter of the distribution. Taking account of RXRA chr9:136355885+ methylation and sex explained 26% of the variance in fat mass in PAH children, and simultaneous analyses showed that the association between RXRA chr9:136355885+ methylation and child’s fat mass was not attenuated after adjustment for birth weight and mother’s age, smoking, and BMI.

TABLE 1

Univariate and multivariate analyses of umbilical cord CpG methylation in relation to child’s adiposity measured by DEXA at age 9 years

eNOS chr7:150315553+ (%)PRXRA chr9:136355885+ (%)PSex (%)PVariance explained (%)
Total fat mass (ln kg)        
 Univariate analysis 21 0.001 21 0.003 39 0.001  
 Multivariate analysis 19 0.002 23 0.001 47 0.001 44 
Percentage fat (ln %)        
 Univariate analysis 13 0.006 14 0.006 42 <0.001  
 Multivariate analysis 11 0.009 14 0.005 48 <0.001 47 
Ratio trunk to limb fat (ln)        
 Univariate analysis 0.008 0.016 12 0.009  
 Multivariate analysis 0.004 0.021 11 0.059 28 
eNOS chr7:150315553+ (%)PRXRA chr9:136355885+ (%)PSex (%)PVariance explained (%)
Total fat mass (ln kg)        
 Univariate analysis 21 0.001 21 0.003 39 0.001  
 Multivariate analysis 19 0.002 23 0.001 47 0.001 44 
Percentage fat (ln %)        
 Univariate analysis 13 0.006 14 0.006 42 <0.001  
 Multivariate analysis 11 0.009 14 0.005 48 <0.001 47 
Ratio trunk to limb fat (ln)        
 Univariate analysis 0.008 0.016 12 0.009  
 Multivariate analysis 0.004 0.021 11 0.059 28 

Values are the percent change in outcome variable per SD change in CpG methylation and P value. Univariate analysis rows show the regression coefficients and P values relating each of eNOS chr7:150315553+, RXRA chr9:136355885+, and sex to outcomes individually, whereas multivariate analysis rows show regression coefficients and P values for a combined analysis of eNOS chr7:150315553+, RXRA chr9:136355885+, and sex in relation to outcomes, together with the variance explained by the multivariate model. For variances explained, n = 55.

In this population we have previously linked lower maternal carbohydrate intake in early pregnancy with higher neonatal adiposity (21). Relating umbilical cord CpG methylation to maternal carbohydrate intake in early pregnancy, higher methylation of RXRA chr9:136355885+, but not of eNOS chr7:150315553+, was associated with a lower maternal carbohydrate intake (Fig. 1). The mother’s early pregnancy fat and protein intakes in early pregnancy were not associated with cord RXRA chr9:136355885+ methylation (P = 0.7 and P = 0.6, respectively), and no additional variance in RXRA chr9:136355885+ methylation was explained by simultaneously taking account of maternal fat and protein intakes in early pregnancy or macronutrient intakes in late pregnancy. Relating umbilical cord CpG methylation to the infant’s size at birth, only three CpGs showed a significant correlation with birth weight adjusted for sex and gestational age (PIK3CD chr1:9609980+ [rp= −0.32, P = 0.013, n = 58], PIK3CD chr1:9635676/79+ [rp= −0.36, P = 0.028, n = 37], and SOD1 chr 21:31853827+ [rp= −0.27, P = 0.029, n = 65]). RXRA chr9:136355885+ methylation showed a weak association with the subjects’ sex and age-adjusted height (rp= 0.24, P = 0.054) but was not associated with the infant’s weight or ponderal index at birth.

FIG. 1.

Lower maternal carbohydrate in early pregnancy is associated with higher umbilical cord RXRA chr9:136355885+ methylation in the PAH cohort. Values are means + SEM.

FIG. 1.

Lower maternal carbohydrate in early pregnancy is associated with higher umbilical cord RXRA chr9:136355885+ methylation in the PAH cohort. Values are means + SEM.

Close modal

Methylation analysis in the SWS cohort.

We next sought to replicate the stronger associations found in the first cohort (eNOS chr7:150315553+ and RXRA chr9:136355885+ CpG methylation and child’s adiposity) in the SWS cohort (20). In these 239 SWS children, eNOS chr7:150315553+ showed no association with adiposity, but RXRA chr9:136355885+ showed remarkably similar and statistically significant associations (fat mass, β = 6% [2–10], P = 0.002; %fat mass, β = 4% [1–7], P = 0.002, both n = 239) to those in PAH children. As in the PAH cohort, these associations were little changed after taking account of birth weight and maternal age, smoking, and BMI. Figures 2 and 3 show the graded associations between cord RXRA chr9:136355885+ methylation and child’s fat mass and %fat mass in the two cohorts.

FIG. 2.

Child’s %fat mass and fat mass at age 9 years increase with higher umbilical cord RXRA chr9:136355885+ methylation in the PAH cohort. Values are means + SEM. *Fat mass and percentage fat mass are preadjusted for sex.

FIG. 2.

Child’s %fat mass and fat mass at age 9 years increase with higher umbilical cord RXRA chr9:136355885+ methylation in the PAH cohort. Values are means + SEM. *Fat mass and percentage fat mass are preadjusted for sex.

Close modal
FIG. 3.

In a second independent cohort, child’s %fat mass and fat mass at age 6 years increase with higher umbilical cord RXRA chr9:136355885+ methylation in SWS subjects. Values are means + SEM. *Fat mass and percentage fat mass are preadjusted for age and sex.

FIG. 3.

In a second independent cohort, child’s %fat mass and fat mass at age 6 years increase with higher umbilical cord RXRA chr9:136355885+ methylation in SWS subjects. Values are means + SEM. *Fat mass and percentage fat mass are preadjusted for age and sex.

Close modal

Pyrosequencing across RXRA chr9:136355885+ in SWS subjects showed no C allele SNPs that could account for differences in methylation observed between individuals.

This study provides novel evidence for the importance of the developmental contribution to later adiposity. We found that greater methylation of RXRA chr9:136355885+ measured at birth was strongly correlated with greater adiposity in later childhood in two independent cohorts. Although we studied a subset of children in both cohorts these were selected on the basis of subject and specimen availability, so it is unlikely that selection bias could explain the relationships observed unless the association between RXRA methylation and adiposity was different in the remainder of the cohort. The data build on animal experiments suggesting that the developmental environment acts through epigenetic processes to exert a strong influence on postnatal body composition and metabolic function (711).

Our study shows that specific components of the epigenetic state at birth predict later childhood adiposity. The associations with adiposity were linked to specific CpGs 5′ to the start site of the selected candidate genes. Although some of the CpGs studied were either within the proximal promoter or close to it, others were more distal and may be exerting effects through the regulation of other genes. Nevertheless the data indicate possible mechanistic pathways, suggesting avenues for future study. Our observation that adjacent or nearby CpGs within the same promoter showed differences in the strength of association with child’s adiposity suggests highly specific changes in the transcriptional regulation of these genes induced by the developmental environment, rather than generalized changes in promoter methylation. Both CpG hyper- (eNOS chr7:150315553+ and RXRA chr9:136355885+) and hypomethylation (SOD1 chr21:31853660/63+) at different sites were associated with body fat distribution, again indicating complexity in transcriptional control. The specificity of the associations between methylation of an individual CpG and both maternal diet and child’s phenotype endorses the concept of a fine control of development by environmental factors via epigenetic processes.

Our observation indicates one potential mechanistic pathway involved, because induction of transcription by RXRA is dependent on its binding to ligands including the peroxisome proliferator–activated receptors, involved in insulin sensitivity, adipogenesis, and fat metabolism (31,32). Moreover, RXRA chr9:136355885+ is located in a region considered to contain positive regulatory elements of transcription (33). Figure 4 shows the proximity of RXRA chr9:136355885+ to proposed binding sites for RXR, MAF, NF-κB, and AP1. Retinoid receptor biology is complex, and increased RXRA methylation might be acting through a variety of pathways (34); however, an association between increased RXRA methylation and adiposity is consistent with the observation of strongly diminished RXRA expression in visceral white adipose tissue from obese mice (35). Moreover, a role for retinoid receptor methylation in developmental influences on later metabolic risk is supported by recent experimental data showing an influence of maternal diet during pregnancy on methylation of LXRA, a heterodimeric partner of RXRA (36).

FIG. 4.

Schematic diagram of the RXRA promoter region, showing the position of the CpG group at RXRA chr9:136355885+ (underlined) and of neighboring transcription factor binding sites.

FIG. 4.

Schematic diagram of the RXRA promoter region, showing the position of the CpG group at RXRA chr9:136355885+ (underlined) and of neighboring transcription factor binding sites.

Close modal

Genome-wide association studies suggest that fixed genetic variation makes a relatively small contribution to risk of obesity, heart disease, and diabetes (1,2); our findings raise the possibility that the developmental environment component may be equally or more important. We excluded the presence of a SNP at RXRA chr9:136355885+ by sequencing, but without genome-wide analysis it is not possible to exclude a genetic effect of distant SNPs, which could influence both DNA methylation of a particular sequence and child’s phenotype. However, even if this were the case, our data clearly indicate that epigenetic measures at birth may have prognostic value. Although epigenetic changes can be dynamic, experimental studies have shown that environmental factors acting on the genotype during development relate to epigenetic profile in adulthood (7,15), and there are longitudinal human studies showing that DNA methylation is often stable over time (37). Such changes can be tissue specific, and in this respect the umbilical cord may be advantageous because it contains a high proportion of fetal vascular tissue and mesenchymal cells, which may be relevant to later adiposity. Furthermore, unlike the placenta it is a tissue in which consistency of sampling between individuals is more likely. Although experimental work in the rat suggests that methylation changes induced by maternal diet can be similar in the umbilical cord and liver (38), further work is needed to determine the relevance of epigenetic changes in human umbilical cord tissue. Recent data show that for some genomic regions methylation appears largely independent of tissue of origin, whereas for others there is a clear tissue-specific dependence (39).

Many epidemiological studies have shown associations between fetal development, through the proxy measure of birth size, and later adiposity and metabolic function (3,4), but the developmental contribution to such phenotypic characteristics has remained uncertain and controversial. This study provides the first estimate of the developmental contribution to phenotype associated with human disease risk based on measures of the underpinning biology: our data for RXRA chr9:136355885+ suggest that a substantial proportion of the variation in adiposity in prepubertal children can be explained by epigenetic measurements made at birth. Although our data are correlative and thus can only imply an association between DNA methylation at birth and later phenotype, the importance of the observation stands irrespective of whether the RXRA methylation is causally related to the development of adiposity. Even if it is simply a noncausal association, the changed epigenetic status provides an objective marker of altered developmental trajectory by the time of birth. Despite the limitations of dietary intake assessment tools, the instrument we used is both validated and provides information that can be used to rank the nutrient intakes of individuals (23).

Methylation of eNOS chr7:150315553+ was associated with later adiposity in the initial cohort only. This may reflect a chance finding, or different maternal characteristics in the two cohorts, such as the greater maternal adiposity and maternal folate supplementation in the SWS cohort (22,40,41). Alternatively, adipocyte proliferation is high during the first year of life but then remains low until a second proliferative phase from age 9–14 years (42). Thus the processes that determine adiposity at age 9 years might differ from those at age 6 years. Nitric oxide synthesis by eNOS promotes preadipocyte differentiation (43). Therefore, a further possible explanation is that eNOS methylation in umbilical cord marks capacity for adipogenesis which has a greater net contribution to adiposity at age 9 years than age 6 years.

Our findings show strong associations between epigenetic markers and childhood total and central body fat. Beyond these simple associations, multivariate analysis indicates that the associations explain substantial proportions of the variances in outcomes, emphasizing an important developmental contribution to phenotypes associated with metabolic dysfunction and disease risk. It is noteworthy that the genes for which we report effects are not imprinted. In vitro fertilization increases risk of imprinting disorders (44), and methylation effects on imprinted genes have been reported in offspring of mothers exposed to famine during various periods of pregnancy (45,46), but with no associations with phenotype reported. The current study implicates the human prenatal environment with epigenetic changes in nonimprinted genes and is the first to link epigenetic status at birth with clinically relevant later phenotypic variation.

Variation in the degree of methylation of nonimprinted genes and in later cardiovascular and metabolic physiology can be induced experimentally by manipulation of the developmental environment, for example by altering maternal nutrition or administering glucocorticoids during pregnancy (8,9), often without necessarily affecting the birth size of the offspring. Although we found associations between some epigenetic markers and birth weight, these were weaker associations than those with later phenotype and were for different markers; moreover, previous reports on the same cohort showed only modest associations between birth weight and later body composition (6). Because birth weight in humans is influenced by multiple factors including the mother’s own birth weight and height and by gestational length (40), epigenetic changes may provide a more sensitive index than birth weight of environmentally induced effects on fetal development.

There are potentially important implications of the strong and replicated association between RXRA chr9:136355885+ methylation and later adiposity. First, the effect is considerably greater than that of factors such as birth weight or maternal body composition, suggesting that epigenetic measurements made in the neonate may be useful predictors of later obesity and other phenotypic outcomes. Second, the association between CpG methylation and child’s adiposity operates within the normal ranges of maternal nutritional state and birth size; this supports the argument that developmental programming is the consequence of an evolved and potentially adaptive process involving the mechanisms of developmental plasticity (10). Indeed the data provide strong evidence supporting a role for developmental plasticity in determining individual risk of metabolic disease. Third, the data suggest that developmental factors may be more significant in contributing to phenotypic variation and disease risk than generally considered. Fourth, the association between RXRA chr9:136355885+ methylation and mother’s carbohydrate intake raises the possibility that conditions in early pregnancy could affect child’s adiposity through this pathway. This provides additional support for the argument that all women of reproductive age should have appropriate nutritional, education, and lifestyle support to improve the health of the next generation. Finally, our data suggest that epigenetic measures at birth may have prognostic value and potential utility for monitoring programs to optimize maternal health and nutrition for long-term benefits to the offspring; however, evaluation of this possibility will require further research correlating methylation measurements in early life with those in later life.

This work was supported by WellChild (previously Children Nationwide); the University of Southampton; the U.K. Medical Research Council and National Institute for Health Research; the British Heart Foundation (to M.A.H. and G.C.B.); Arthritis Research U.K.; the New Zealand National Research Centre for Growth and Development (to P.D.G. and A.S.); and the Singapore Agency for Science, Technology and Research (to B.S.E.). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

No potential conflicts of interest relevant to this article were reported.

K.M.G. conceived of the study, undertook the statistical analysis, wrote the first draft of the manuscript, and reviewed and edited the manuscript. A.S. undertook the molecular biology and reviewed and edited the manuscript. P.D.G. conceived of the study, wrote the first draft of the manuscript, and reviewed and edited the manuscript. K.A.L. conceived of the study, undertook the molecular biology, and reviewed and edited the manuscript. G.C.B. conceived of the study and reviewed and edited the manuscript. C.M., J.R., J.L.S.-J., and E.G. undertook the molecular biology and reviewed and edited the manuscript. S.R.C. undertook the statistical analysis and reviewed and edited the manuscript. B.S.E. undertook the molecular biology and reviewed and edited the manuscript. C.R.G. coordinated the phenotyping of the children and reviewed and edited the manuscript. H.M.I. coordinated the phenotyping of the children, assisted in the statistical analysis, and reviewed and edited the manuscript. C.C. assisted in the statistical analysis and reviewed and edited the manuscript. M.A.H. conceived of the study, wrote the first draft of the manuscript, and reviewed and edited the manuscript.

Sue Copland and Farhad Shaffei, University of Auckland, provided technical assistance.

1.
Manolio
TA
,
Collins
FS
,
Cox
NJ
, et al
.
Finding the missing heritability of complex diseases
.
Nature
2009
;
461
:
747
753
[PubMed]
2.
Craddock
N
,
Hurles
ME
,
Cardin
N
, et al
;
Wellcome Trust Case Control Consortium
.
Genome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls
.
Nature
2010
;
464
:
713
720
[PubMed]
3.
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
:
61
73
[PubMed]
4.
Godfrey KM, Barker DJP. Fetal programming and adult health. Public Health Nutrition 2001;2:611–624
5.
Ravelli
GP
,
Stein
ZA
,
Susser
MW
.
Obesity in young men after famine exposure in utero and early infancy
.
N Engl J Med
1976
;
295
:
349
353
[PubMed]
6.
Gale
CR
,
Javaid
MK
,
Robinson
SM
,
Law
CM
,
Godfrey
KM
,
Cooper
C
.
Maternal size in pregnancy and body composition in children
.
J Clin Endocrinol Metab
2007
;
92
:
3904
3911
[PubMed]
7.
Burdge
GC
,
Slater-Jefferies
J
,
Torrens
C
,
Phillips
ES
,
Hanson
MA
,
Lillycrop
KA
.
Dietary protein restriction of pregnant rats in the F0 generation induces altered methylation of hepatic gene promoters in the adult male offspring in the F1 and F2 generations
.
Br J Nutr
2007
;
97
:
435
439
[PubMed]
8.
Drake
AJ
,
Walker
BR
,
Seckl
JR
.
Intergenerational consequences of fetal programming by in utero exposure to glucocorticoids in rats
.
Am J Physiol Regul Integr Comp Physiol
2005
;
288
:
R34
R38
[PubMed]
9.
Champagne
FA
,
Meaney
MJ
.
Transgenerational effects of social environment on variations in maternal care and behavioral response to novelty
.
Behav Neurosci
2007
;
121
:
1353
1363
[PubMed]
10.
Gluckman
PD
,
Hanson
MA
,
Beedle
AS
.
Non-genomic transgenerational inheritance of disease risk
.
Bioessays
2007
;
29
:
145
154
[PubMed]
11.
Godfrey
KM
,
Lillycrop
KA
,
Burdge
GC
,
Gluckman
PD
,
Hanson
MA
.
Epigenetic mechanisms and the mismatch concept of the developmental origins of health and disease
.
Pediatr Res
2007
;
61
:
5R
10R
[PubMed]
12.
Gluckman
PD
,
Hanson
MA
.
Developmental plasticity and human disease: research directions
.
J Intern Med
2007
;
261
:
461
471
[PubMed]
13.
Bateson
P
,
Barker
D
,
Clutton-Brock
T
, et al
.
Developmental plasticity and human health
.
Nature
2004
;
430
:
419
421
[PubMed]
14.
Maleszka
R
.
Epigenetic integration of environmental and genomic signals in honey bees: the critical interplay of nutritional, brain and reproductive networks
.
Epigenetics
2008
;
3
:
188
192
[PubMed]
15.
Lillycrop
KA
,
Phillips
ES
,
Jackson
AA
,
Hanson
MA
,
Burdge
GC
.
Dietary protein restriction of pregnant rats induces and folic acid supplementation prevents epigenetic modification of hepatic gene expression in the offspring
.
J Nutr
2005
;
135
:
1382
1386
[PubMed]
16.
Bogdarina
I
,
Welham
S
,
King
PJ
,
Burns
SP
,
Clark
AJL
.
Epigenetic modification of the renin-angiotensin system in the fetal programming of hypertension
.
Circ Res
2007
;
100
:
520
526
[PubMed]
17.
Lillycrop
KA
,
Slater-Jefferies
JL
,
Hanson
MA
,
Godfrey
KM
,
Jackson
AA
,
Burdge
GC
.
Induction of altered epigenetic regulation of the hepatic glucocorticoid receptor in the offspring of rats fed a protein-restricted diet during pregnancy suggests that reduced DNA methyltransferase-1 expression is involved in impaired DNA methylation and changes in histone modifications
.
Br J Nutr
2007
;
97
:
1064
1073
[PubMed]
18.
Lillycrop
KA
,
Phillips
ES
,
Torrens
C
,
Hanson
MA
,
Jackson
AA
,
Burdge
GC
.
Feeding pregnant rats a protein-restricted diet persistently alters the methylation of specific cytosines in the hepatic PPAR alpha promoter of the offspring
.
Br J Nutr
2008
;
100
:
278
282
[PubMed]
19.
Burdge
GC
,
Lillycrop
KA
,
Phillips
ES
,
Slater-Jefferies
JL
,
Jackson
AA
,
Hanson
MA
.
Folic acid supplementation during the juvenile-pubertal period in rats modifies the phenotype and epigenotype induced by prenatal nutrition
.
J Nutr
2009
;
139
:
1054
1060
[PubMed]
20.
Gluckman
PD
,
Lillycrop
KA
,
Vickers
MH
, et al
.
Metabolic plasticity during mammalian development is directionally dependent on early nutritional status
.
Proc Natl Acad Sci USA
2007
;
104
:
12796
12800
[PubMed]
21.
Godfrey
KM
,
Barker
DJ
,
Robinson
S
,
Osmond
C
.
Maternal birthweight and diet in pregnancy in relation to the infant’s thinness at birth
.
Br J Obstet Gynaecol
1997
;
104
:
663
667
[PubMed]
22.
Inskip
HM
,
Godfrey
KM
,
Robinson
SM
,
Law
CM
,
Barker
DJ
,
Cooper
C
;
SWS Study Group
.
Cohort profile: The Southampton Women’s Survey
.
Int J Epidemiol
2006
;
35
:
42
48
[PubMed]
23.
Robinson
S
,
Godfrey
K
,
Osmond
C
,
Cox
V
,
Barker
DJP
.
Evaluation of a food frequency questionnaire used to assess nutrient intakes in pregnant women
.
Eur J Clin Nutr
1996
;
50
:
302
308
[PubMed]
24.
Weber
M
,
Hellmann
I
,
Stadler
MB
, et al
.
Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome
.
Nat Genet
2007
;
39
:
457
466
[PubMed]
25.
Burdge GC, Lillycrop KA. Nutrition, epigenetics, and developmental plasticity: implications for understanding human disease. Annu Rev Nutr 2010;30:315–339
26.
Lillycrop
KA
,
Rodford
J
,
Garratt
ES
, et al
.
Maternal protein restriction with or without folic acid supplementation during pregnancy alters the hepatic transcriptome in adult male rats
.
Br J Nutr
2010
;
103
:
1711
1719
[PubMed]
27.
Li
LC
,
Dahiya
R
.
MethPrimer: designing primers for methylation PCRs
.
Bioinformatics
2002
;
18
:
1427
1431
[PubMed]
28.
Ehrich
M
,
Nelson
MR
,
Stanssens
P
, et al
.
Quantitative high-throughput analysis of DNA methylation patterns by base-specific cleavage and mass spectrometry
.
Proc Natl Acad Sci USA
2005
;
102
:
15785
15790
[PubMed]
29.
Alderborn
A
,
Kristofferson
A
,
Hammerling
U
.
Determination of single-nucleotide polymorphisms by real-time pyrophosphate DNA sequencing
.
Genome Res
2000
;
10
:
1249
1258
[PubMed]
30.
Armitage
P
,
Berry
G
.
Statistical Methods in Medical Research
.
Oxford
,
Blackwell Science Ltd
,
2002
31.
Alvarez
R
,
Checa
M
,
Brun
S
, et al
.
Both retinoic-acid-receptor- and retinoid-X-receptor-dependent signalling pathways mediate the induction of the brown-adipose-tissue-uncoupling-protein-1 gene by retinoids
.
Biochem J
2000
;
345
:
91
97
[PubMed]
32.
Sugden
MC
,
Holness
MJ
.
Role of nuclear receptors in the modulation of insulin secretion in lipid-induced insulin resistance
.
Biochem Soc Trans
2008
;
36
:
891
900
[PubMed]
33.
Li
G
,
Yin
W
,
Chamberlain
R
, et al
.
Identification and characterization of the human retinoid X receptor alpha gene promoter
.
Gene
2006
;
372
:
118
127
[PubMed]
34.
Mascrez
B
,
Ghyselinck
NB
,
Chambon
P
,
Mark
M
.
A transcriptionally silent RXRalpha supports early embryonic morphogenesis and heart development
.
Proc Natl Acad Sci USA
2009
;
106
:
4272
4277
[PubMed]
35.
Lefebvre
B
,
Benomar
Y
,
Guédin
A
, et al
.
Proteasomal degradation of retinoid X receptor alpha reprograms transcriptional activity of PPARgamma in obese mice and humans
.
J Clin Invest
2010
;
120
:
1454
1468
[PubMed]
36.
van Straten
EM
,
Bloks
VW
,
Huijkman
NC
, et al
.
The liver X-receptor gene promoter is hypermethylated in a mouse model of prenatal protein restriction
.
Am J Physiol Regul Integr Comp Physiol
2010
;
298
:
R275
R282
[PubMed]
37.
Talens
RP
,
Boomsma
DI
,
Tobi
EW
, et al
.
Variation, patterns, and temporal stability of DNA methylation: considerations for epigenetic epidemiology
.
FASEB J
2010
;
24
:
3135
3144
[PubMed]
38.
Burdge
GC
,
Hanson
MA
,
Slater-Jefferies
JL
,
Lillycrop
KA
.
Epigenetic regulation of transcription: a mechanism for inducing variations in phenotype (fetal programming) by differences in nutrition during early life?
Br J Nutr
2007
;
97
:
1036
1046
[PubMed]
39.
Ollikainen
M
,
Smith
KR
,
Joo
EJ
, et al
.
DNA methylation analysis of multiple tissues from newborn twins reveals both genetic and intrauterine components to variation in the human neonatal epigenome
.
Hum Mol Genet
2010
;
19
:
4176
4188
[PubMed]
40.
Godfrey
KM
,
Robinson
S
,
Barker
DJ
,
Osmond
C
,
Cox
V
.
Maternal nutrition in early and late pregnancy in relation to placental and fetal growth
.
BMJ
1996
;
312
:
410
414
[PubMed]
41.
Duggleby
SL
,
Jackson
AA
,
Godfrey
KM
,
Robinson
SM
,
Inskip
HM
;
Southampton Women’s Survey Study Group
.
Cut-off points for anthropometric indices of adiposity: differential classification in a large population of young women
.
Br J Nutr
2009
;
101
:
424
430
[PubMed]
42.
Baum
D
,
Beck
RQ
,
Hammer
LD
,
Brasel
JA
,
Greenwood
MR
.
Adipose tissue thymidine kinase activity in man
.
Pediatr Res
1986
;
20
:
118
121
[PubMed]
43.
Yan
H
,
Aziz
E
,
Shillabeer
G
, et al
.
Nitric oxide promotes differentiation of rat white preadipocytes in culture
.
J Lipid Res
2002
;
43
:
2123
2129
[PubMed]
44.
Amor
DJ
,
Halliday
J
.
A review of known imprinting syndromes and their association with assisted reproduction technologies
.
Hum Reprod
2008
;
23
:
2826
2834
[PubMed]
45.
Heijmans
BT
,
Tobi
EW
,
Stein
AD
, et al
.
Persistent epigenetic differences associated with prenatal exposure to famine in humans
.
Proc Natl Acad Sci USA
2008
;
105
:
17046
17049
[PubMed]
46.
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
:
4046
4053
[PubMed]
Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered. See http://creativecommons.org/licenses/by-nc-nd/3.0/ for details.

Supplementary data