Skip to main content

Childhood DNA methylation as a marker of early life rapid weight gain and subsequent overweight

Abstract

Background

High early postnatal weight gain has been associated with childhood adiposity; however, the mechanism remains unknown. DNA methylation is a hypothesised mechanism linking early life exposures and subsequent disease. However, epigenetic changes associated with high early weight gain have not previously been investigated. Our aim was to investigate the associations between early weight gain, peripheral blood DNA methylation, and subsequent overweight/obese. Data from the UK Avon Longitudinal study of Parents and Children (ALSPAC) cohort were used to estimate associations between early postnatal weight gain and epigenome-wide DNA CpG site methylation (Illumina 450 K Methylation Beadchip) in blood in childhood (n = 125) and late adolescence (n = 96). High weight gain in the first year (a change in weight z‐scores > 0.67), both unconditional (rapid weight gain) and conditional on birthweight (rapid thrive), was related to individual CpG site methylation and across regions using the meffil pipeline, with and without adjustment for cell type proportions, and with 5% false discovery rate correction. Variation in methylation at high weight gain-associated CpG sites was then examined with regard to body composition measures in childhood and adolescence. Replication of the differentially methylated CpG sites was sought using whole-blood DNA samples from 104 children from the UK Southampton Women’s Survey.

Results

Rapid infant weight gain was associated with small (+ 1% change) increases in childhood methylation (age 7) for two distinct CpG sites (cg01379158 (NT5M) and cg11531579 (CHFR)). Childhood methylation at one of these CpGs (cg11531579) was also higher in those who experienced rapid weight gain and were subsequently overweight/obese in adolescence (age 17). Rapid weight gain was not associated with differential DNA methylation in adolescence. Childhood methylation at the cg11531579 site was also suggestively associated with rapid weight gain in the replication cohort.

Conclusions

This study identified associations between rapid weight gain in infancy and small increases in childhood methylation at two CpG sites, one of which was replicated and was also associated with subsequent overweight/obese. It will be important to determine whether loci are markers of early rapid weight gain across different, larger populations. The mechanistic relevance of these differentially methylated sites requires further investigation.

Background

Children are becoming obese at younger ages [1], suggesting that factors in early life may play a role in obesity development. The developmental origins of health and disease (DOHaD) hypothesis proposes that early life environmental exposures have the potential to modify the risk of later-life diseases, such as obesity [2]. Rapid weight gain (RWG) is an early life factor that has been consistently associated with childhood adiposity both dependently and independently of birthweight [3,4,5]. Weight gain in the first year specifically (opposed to change in weight over periods greater or less than 1 year) has been found to be most predictive of childhood obesity [6], suggesting this is a critical period.

Given the responsiveness to environmental stimuli, the capacity to alter gene expression and their stability over time, epigenetic changes are a proposed mechanism underlying the DOHaD hypothesis. Through programming effects, epigenetic marks laid down at an early developmental stage could elicit effects at a later stage [7]. DNA methylation (DNAm) is the most stable and widely studied epigenetic modification and is a key mechanism regulating gene expression. DNAm involves the covalent addition of a methyl group to cytosine residues adjacent to guanine in DNA (CpG sites) and is associated with changes in gene transcription [8]. If early life factors lead to stable changes in DNAm, these changes could be used as biomarkers and to identify individuals who may benefit from intervention prior to disease onset.

BMI has been associated with variation in DNAm from birth to adulthood [9,10,11,12]. Epigenome-Wide Association Studies (EWAS) are a comprehensive approach to identify epigenetic variation associated with a biological trait or exposure [13, 14]. In EWAS, other early life risk factors for childhood obesity such as birthweight and maternal BMI have been associated with variation in DNAm [15, 16]. To our knowledge, there have been no EWAS to date on early life rapid growth and DNAm.

Our first aim was to identify DNAm changes associated with early life growth. In this study we hypothesised that early life rapid growth is associated with DNAm changes. Using epigenome-wide DNAm array data from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort, we investigated if early life rapid growth is associated with variation in childhood methylation, and if methylation changes persist into adolescence. As catch up growth is more likely in low birthweight infants, we investigated rapid growth both adjusted (rapid thrive, RT) [17] and unadjusted for birthweight (rapid weight gain, RWG) [4]. We also examined differential methylation in a subset of known BMI-associated CpG loci [12] with the aim of identifying differentially methylated loci more likely to be related to body composition, and by analysing fewer loci, to offset the multiple comparison problem often associated with null findings in EWAS.

An important consideration is that not all children with rapid infancy weight gain will have increased adiposity in childhood [18], and previous studies have highlighted the necessity to distinguish infants at greatest risk of overweight/obesity. Therefore, we also aimed to explore if differential methylation was associated with later life BMI and overweight/obesity in those who experienced early rapid weight gain to determine potential risk markers. Finally, we sought replication in an independent cohort, the UK Southampton Women’s survey.

Methods

Cohort

We performed our initial analysis in the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort, which has detailed early life, anthropometric and epigenome-wide DNAm data at multiple time points. This ongoing longitudinal birth cohort, based in Bristol, England (UK), initially invited pregnant women resident in Avon, UK, with expected dates of delivery 1 April 1991 to 31 December 1992 [19, 20]. There were 14,541 initial pregnancies enrolled (for these at least one questionnaire has been returned or a “Children in Focus” clinic had been attended by 19/07/99), with a total of 14,676 foetuses, resulting in 14,062 live births and 13,988 children who were alive at 1 year of age.

The cohort has had extensive questionnaires as well as clinical assessments (including measures of height & weight) throughout childhood. The Accessible Resource for Integrated Epigenomic Studies (ARIES) is a subset of the ALSPAC cohort [21], for which epigenome-wide DNAm analysis was carried out for 1,018 mother and child pairs. Ethical approval for the study was obtained from the ALSPAC Ethics and Law Committee and local research ethics committees. The study website contains details of all the data that are available through a fully searchable data dictionary accessible at www.bris.ac.uk/alspac/researchers/data-access/data-dictionary/.

Early life data

Birthweight and gestational age were taken from medical records. At 12 months, infants were weighed using the Seca 724 (or Seca 835 for children who could only be weighed with a parent). Birthweight and weight-for-age [12 months] z-scores were calculated using the British 1990 growth reference [22] and were used to determine the rapid growth variables: RWG and RT. Whilst birthweight is an important factor in childhood obesity [23], it has previously been examined in the ALSPAC cohort [15, 24]. However, in order to differentiate the effects of birthweight on early life rapid growth, both RWG (not conditional on birthweight) and RT (conditional on birthweight) in the first year were examined (Fig. 1). Conditional weight gain (also known as thrive index) accounts for normal catch-up growth from low birthweight as a linear measure of weight gain adjusted for regression to the mean [17]. Rapid thrive was determined as z-score12m − r × z-scorebirth [17], where r is the cohort regression coefficient (r = 0.35) of birthweight-z on weight-z (12 months). RWG is most often defined as a > + 0.67 standard deviation change in weight-for-age z-score, equivalent to crossing a growth centile band on a standard child growth chart [4]. Both RWG and RT were analysed as a dichotomised variables of a > + 0.67 standard deviation change (Fig. 1).

Fig. 1
figure 1

Description and measurement of rapid infant weight gain exposures, and epigenome-wide DNA methylation and anthropometric outcomes, and an overview of the analysis. Participants were measured at around age 7 (mean 7.4 (SD 0.10) years) and 17 (mean 17.7 (SD 0.31) years). RWG, rapid weight gain; RT, rapid thrive; BMIz, BMI z-scores; OWOB, overweight/obese; DNAm, DNA methylation; EWAS, Epigenome-wide association study; UK90, the British 1990 growth reference

Anthropometric data

Anthropometric measures at (approximately) age 7 and 17 were also analysed as outcomes in analyses. At age 7, height was measured to the nearest millimetre without shoes or socks using a Holtain stadiometer (Holtain Ltd, Crymych, Pembs, UK), whilst weight was measured using Tanita THF 300GS body fat analyser and weighing scales (Tanita UK Ltd, Yewsley, Middlesex, UK). At age 17, height was measured with a Harpenden stadiometer to the nearest mm, and weight using the Tanita Body Fat Analyser (Model TBF 401A) to the nearest 50 g.

BMIz and overweight/obese were examined as outcomes at both time points. Body mass index (BMI) was calculated as weight (kg) divided by height (m) squared and was transformed to age and sex standardised BMI z-scores using the British 1990 growth reference with the zanthro program in Stata [25]. Clinical cut-offs were used to determine weight categories, whereby healthy weight was between the 2nd and 91st centiles, and overweight/obese as greater than the 91st centile [26].

Epigenetic data

ALSPAC collected (peripheral) blood at ages 7 and 17 and DNA was extracted. Epigenome-wide DNAm at specific CpG sites was measured for ~ 1000 individuals using the Infinium® HumanMethylation450K BeadChip assay (Illumina, Inc., CA, USA). DNAm data were pre-processed, including background correction and subset quantile normalization using the pipeline described by Touleimat and Tost [27] (further details in the ARIES cohort profile [21]). Estimation of white blood cell counts (CD8T cells, CD4T cells, Natural Killer, B cells, Monocytes and Granulocytes) was done using the Houseman algorithm [28]. Cross-reactive and polymorphic probes identified by Chen et al. [29] and probes on sex chromosomes were removed prior to downstream analysis (n = 453,723 probes). Due to few participants with non-Caucasian or missing ethnicity, and few non-singleton births, these were removed.

Statistical analysis

To examine if childhood or adolescent DNAm in peripheral blood (around ages 7 and 17) was associated with early life rapid growth, three different analyses were undertaken, including: analysis of differentially methylation positions, differentially methylated regions, and differentially methylated positions in a subset of candidate loci.

For each analysis, DNAm was the outcome and the independent variable was either RWG or RT, with models adjusted for sex and age at blood collection. Cell type composition is a significant source of variation in DNAm analysis [30]; however chronic, low level inflammation is a component of the obesity phenotype; therefore, to find novel biomarkers associated with this phenotype, DNAm was investigated in models with and without adjustment for cell composition (all 6 cell types). Correction for multiple testing was applied using a false discovery rate (FDR) threshold of p < 0.05 [31].

First, epigenome-wide association studies (EWAS) were conducted for rapid growth (first year) and DNAm outcomes in childhood (age 7.4) and adolescence (age 17.6) in the Meffil R package [32]. Estimation of differentially methylated sites was carried out using the beta values as the outcome and rapid growth (RWG/RT) as the exposure adjusted for age and sex. Surrogate variable analysis (SVA) and independent surrogate variable analysis (ISVA) methods were utilised to control for unmodelled or unknown confounding factors (such as batch) [33, 34]. Meffil simultaneously computes unadjusted, adjusted, SVA and ISVA models, thereby allowing results to be compared [32]. In order to minimise the influence of outliers in methylation data, beta values were winsorised at the level of 5% (95th percentile cut-off).

Second, in order to detect a DNAm signature of rapid growth, differentially methylated regions (DMRs) were analysed. DMRs, which are stretches or clusters of neighbouring CpG probes, may have more of a functional effect on gene expression than individual CpG loci [35]. Additionally, if changes in DNAm are small but persistent across a region, there may be more statistical power to detect them collectively as DMRs [36]. The R package DMRcate was used for the estimation of DMRs [37], using the M-values and default settings. The surrogate variables (that were calculated using Meffil in the EWAS models) were included as covariates in the DMR models.

Finally, in order to focus on loci with anticipated associations with adiposity, a candidate gene approach was taken using a subset of CpG sites robustly associated with BMI. The candidate sites were selected from a large-scale meta-analysis EWAS which utilised data from multiple cohorts of European and Indian-Asian descent [12]. After validation, 187 CpG sites were associated with BMI. We conducted EWAS for both RWG and RT using the 187 identified CpG sites as candidates, at both time points using the Meffil R package [32].

Any significantly differentially methylated sites were analysed further to determine whether DNAm was also associated with body composition (at age 7 and 17). Using linear regression, significant CpG sites were examined with respect to childhood and adolescent BMI (dependent variable). Similarly, as differentially methylated loci were identified in the RWG EWAS, the CpG sites were also examined with regard to RT in adjusted linear models. All models were adjusted for age and sex. Differences in DNAm by phenotype (RWG, overweight/obese) were assessed using ANOVA tests between groups (with Bonferroni correction for multiple testing).

All EWAS and bioinformatic analyses were done in Rstudio version 3.3.2. The human reference genome (GRCh37/hg19 assembly) was used to determine the location and features of the gene region using the UCSC Genome Browser [38]. Recoding of the variables and statistical analysis of differentially methylated sites was done in Stata version 15.1 (StataCorp, College Station, Texas, USA).

Replication analysis

CpG sites with FDR P-values < 0.05 were carried forward for replication using DNA methylation data from the children from the Southampton Women’s Survey (SWS), a similar UK-based cohort [39]. Blood methylation measures (EPIC array) and early life weight data for 104 of the SWS children (age = 11–13 years, all Caucasian) were available. Further details of the replication cohort and methylation data processing are provided in the Additional file 1. The exact same models were estimated for the differentially methylated CpG sites in childhood (age 11–13) and early life rapid growth (RWG and RT, 0–12 months) using the meffil R package with both SVA and ISVA. Models were all adjusted for child’s sex and age at DNAm measurement and were run with and without adjustment for cell type composition.

Results

Descriptive characteristics

Early life growth and 450 K array data were available for 125 ALSPAC children at age 7 and 96 at age 17 (Table 1). Weight at 12 months was measured in a fraction (n = 1,432, 9.3%) of ALSPAC children, thereby limiting the sample size for the methylation analysis. The ARIES sub-sample was mostly representative of the main study population; however, ARIES mothers were slightly older, less likely to have a manual occupation and were less likely to smoke during pregnancy [21]. At age 7, 13% of study members were overweight/obese, and 19% were at age 17 (Additional file 2: Table 1).

Table 1 The proportion of individuals in the study sample with RWG or RT at age 7 and 17

EWAS results

We observed associations (PFDR < 0.05) between RWG and individual CpG loci at age 7. Across the adjusted models, there were 4 associations identified for RWG (PFDR < 0.05) corresponding to 2 unique CpG sites (Table 2). These loci were cg01379158 (NT5M) and cg11531579 (CHFR), and both were associated with a 1% increase in methylation (PFDR = 0.02) in those who had RWG. In the models without cell counts, two of the model p values were also below the Bonferroni p value threshold (1.04 × 10–7). There were no associations (PFDR < 0.05) between RT and individual CpG loci at age 7 in the EWAS. We examined whether methylation at the RWG-associated CpG sites was also associated with RT using regression analysis; however the magnitude of the coefficients was lower and less statistically robust than for RWG (Additional file 2: Table 2). There was no evidence that RWG or RT was associated (PFDR < 0.1) with differential DNAm in adolescence for the EWAS, or for the 2 CpG sites identified as differentially methylated at age 7 using linear regression (Additional file 3). In the DMR analysis, there were no overall DMRs identified; all Stouffer corrected p values were non-significant, suggesting a lack of consistency in the direction of the methylation changes.

Table 2 Associations (FDR p < 0.05) between individual CpG sites (age 7, n = 453,723) and the early life growth in models

The candidate gene analysis

The aim of the candidate gene analyses was to select CpG loci already known to be associated with the outcome phenotype of interest (body composition). Using a smaller subset of loci as candidates has the advantage of reducing the stringent p value threshold when correcting for multiple tests. The candidate gene analysis utilised findings from a consortium, which integrated data from 4 discovery cohorts and replicated findings in 9 cohorts, and found 187 validated methylation markers associated with BMI [12].

The associations between the candidate CpG loci (n = 187) and early life rapid growth were examined using the ALSPAC methylation childhood and adolescent data; however there were no associations identified (Bonferroni p value > 3 × 10–4).

Investigating phenotypic differences in DNAm associated with RWG

Methylation at either CpG site was not directly associated with BMIz at age 7 or 17 in regression analyses. For both CpG sites, highest methylation was in those who had RWG and were subsequently overweight/obese (compared to healthy weight), both in childhood and adolescence (Additional file 2: Table 3). At the cg11531579 site, childhood methylation was higher in those who were subsequently overweight/obese (age 17) and had RWG compared to those who did not have RWG (Fig. 2, ANOVA p < 0.05). However, the sample sizes for these groups were small and therefore results are inconclusive.

Table 3 Associations between early life rapid growth and childhood methylation in the replication cohort (SWS children age 11–13, n = 104)
Fig. 2
figure 2

Childhood methylation and adolescent body composition. Methylation at the cg01379158 locus (left plot), and cg11531579 (right plot). There were no differences between groups for cg01379158 ANOVA (p > 0.05). The cg11531579 plot presents tests for significance (p values) from the Kruskal–Wallis test with Bonferroni correction for multiple testing. RWG, rapid weight gain; OWOB, overweight/obese

Furthermore, those who were healthy weight at age 7 but were overweight/obese at age 17 had higher methylation at age 7 (Fig. 3), whereas those who had RWG but were a healthy weight (at either time point) had consistently lower levels of methylation. On average, methylation was lower in those who did not have RWG regardless of weight status. Although group sizes were small (n = 6), methylation at age 7 could have indicated future risk of overweight/obesity in the ‘high-risk’ group.

Fig. 3
figure 3

Pathways of mean methylation (cg11531579) at age 7 (%) and body composition (at ages 7 and 17). Dotted outline indicates the ‘high-risk’ individuals (i.e. those who had RWG and were subsequently overweight/obese (OWOB). Sample sizes are for those with complete data. Group sizes were small for some phenotypes.

Methylation change over time within individuals

Overall, the two CpG sites which were positively associated with RWG tended to increase in methylation over time from childhood to adolescence within individuals (Fig. 4). However, when stratifying by RWG, in those who had RWG there was a decrease in methylation over time compared to those who did not experience RWG, particularly for cg11531579 (p < 0.001, Fig. 4).

Fig. 4
figure 4

Change in methylation at the loci within individuals from age 7 to 17 by RWG. Test for differences using the Student’s t-test. Those who did not have RWG (n = 60) demonstrated small mean increases (cg01379158, + 0.6%; cg11531579, + 1.2%) in methylation, whereas those who had RWG (n = 34) demonstrated small (cg01379158, − 0.4%; cg11531579, − 0.6%) decreases in methylation between ages 7 and 17. Error bars represent standard deviation.

Childhood methylation in the replication cohort

Replication of the significant CpG sites was carried out using a similar UK-based cohort with data on growth in early life and epigenetic data in childhood (n = 104 at age 12). Compared to ALSPAC children, fewer experienced rapid growth in the first year in the SWS cohort; 29.7% had RWG (30/101) and 22.8% (23/101) had RT. Similar to the findings in ALSPAC, there was evidence of an association between RWG and DNA methylation at cg11531579 in the ISVA (p = 0.02) and SVA (p = 0.04) models, although the coefficients were smaller (0.005 and 0.004 in the ISVA and SVA models, respectively) (Table 3). There was no association for the either the cg01379158 site or the RT models.

Genomic location of the differentially methylated CpG sites

The CpG site: cg01379158 was located upstream of the transcriptional start site in a CpG island (chr17:17,206,527–17,207,306). The nearest gene to cg01379158 is NT5M, also known as 5′,3′-Nucleotidase, Mitochondrial. The second differentially methylated CpG site (cg11531579, pFDR < 0.05, Table 2) is located within a CpG island on chromosome 12: upstream 30+ kilobases is the protein coding gene Checkpoint With Forkhead And Ring Finger Domains (CHFR) (Table 2), whilst 558 base pairs downstream is a small (2 exons) non-coding region (AK055957), for which there is limited information.

Discussion

Summary

In this study, we identified that RWG in the first year of life was associated with small but significant increases in childhood DNAm (age 7) at two CpG sites (cg01379158 and cg11531579). The highest levels of methylation at the cg11531579 locus (age 7) were in those who had RWG and were either currently (age 7) or subsequently overweight/obese (age 17). Furthermore, there was suggestive evidence that this site was differentially methylated in the replication cohort. We did not find evidence of differentially methylated regions, of differential methylation associated with RT, or of differentially methylated BMI-associated candidate sites.

Interpretation and comparison with previous findings

To our knowledge this is the first EWAS to identify differential DNAm associated with early life rapid weight gain. DNAm at the locus near CHFR (cg11531579) was higher in those who had RWG and who were overweight/obese in childhood or adolescence. Early life RWG was associated with small changes in childhood methylation, but not with methylation in adolescence, which could have been partly due to a smaller sample size or a lack of persistence in the differential DNAm seen in early life. Indeed, in those who had RWG there was a decrease in methylation over time (age 7 to 17), which may reflect the ‘recovery’ of hypermethylation in childhood. This phenomenon of attenuation of DNAm over time from signals detected in early life has been reported previously [24]. There is the possibility that methylation changes may be greater earlier in childhood closer to the timing of the exposure.

There are inherent links between birthweight and postnatal growth, and birthweight associated DNAm changes are also often related to growth control [40]. Rapid thrive accounts for catch-up growth from low birthweight, whereas RWG includes some of the effects of low birthweight. Although birthweight influences RWG, associations between RWG and adiposity remain after adjustment for birthweight [41]. As associations were stronger between DNAm (at the identified CpG sites) and RWG (rather than RT), it is plausible that methylation at these sites also encompasses some of the effects of catch-up growth from low birthweight.

Neither early life growth nor birthweight have been previously associated with DNAm at either of the identified sites. We searched the EWAS Catalogue (https://ewascatalog.org/) to assess whether any of the CpG sites had been previously identified in other EWAS, with results suggesting that both of these loci may be linked to bone composition and cholesterol metabolism, factors which could plausibly be linked to growth. As the DNAm changes identified were small, it is perhaps speculative to discuss the impact on gene expression.

The cg11531579 site, which was positively associated with RWG in ALSPAC and the SWS children, has nearby transcripts with cancer-associated roles [42,43,44,45,46,47,48]. The CHFR (cg11531579) gene encodes a E3 ubiquitin-protein ligase which regulates the cell cycle at the antephase checkpoint (prior to cell division) [42]. Differential epigenetic regulation of CHFR has been identified in cancer as a result of promoter hypermethylation [43, 44], or deacetylation of histones in the promoter region [45]; however it is unclear whether changes in expression are a cause or consequence of cancer. The CpG site is located within a DNAse I hypersensitivity cluster, which may suggest a transcription factor binding region, and a H3K27Ac histone mark, which is often found near regulatory elements and is thought to be a transcription enhancer, suggesting regulatory functions. Downstream of cg11531579 is AK055957, a small non-coding RNA regulatory sequence, with an uncharacterised biological role. Recently, this CpG (in combination with others) has been identified as a potential DNAm biomarker for use in detection panels for hepatocellular carcinoma [46] and pancreatic ductal adenocarcinoma [47], and is differentially methylated in children with acute myeloid leukaemia post-chemotherapy (− 0.24 change in beta value, p = 0.004) [48]. These findings suggest this locus may have a role in carcinogenesis, which may suggest a weak link with rapid growth.

The cg01379158 site is located in the transcriptional start site of the NT5M gene, a gene involved in nucleotide metabolism. The gene is located on chromosome 17 in the Smith-Magenis syndrome-critical region, which is a rare condition characterised by inverse circadian rhythm and disturbed sleep, factors which have also been linked to child obesity [49, 50]. At this CpG site there was no evidence of differential methylation in the SWS cohort, and methylation at the cg11531579 site was around half that observed in ALSPAC. There are several possible reasons for this, such as the moderate sample sizes, or smaller proportion of those who had rapid growth in the SWS. From the ALSPAC data it was evident that there is loss of methylation over time (age 7 to 17) at these loci in those with RWG; therefore as the SWS were also slightly older (+ 5 years) this may also explain less differential methylation. This may suggest that a biomarker for early rapid growth could have greater utility in early childhood and warrants further investigation in younger cohorts.

There were no associations between previously reported BMI-related CpG sites and RWG or RT arising from our candidate analysis, which could be for various reasons. First, the candidate loci mapped to genes with specific roles, which could be different to the mechanisms and pathways of RWG. Secondly, although some of the associations have been replicated in pre-school children [9], primarily the candidate loci were relevant to an adult population, whereas this cohort was sampled and analysed at a much younger age. Finally, RWG has been associated with subsequent changes in BMI, whereas Wahl and colleagues conclude that the majority of the identified BMI-related CpGs were a consequence (rather than a cause) of changes in BMI [12]. Thus, if rapid growth was associated with DNA changes in this subset of CpGs, this perhaps would have been more likely to have been as a consequence of changes in BMI. Indeed, current evidence suggests that the direction of the effect is from BMI to DNAm [12, 51]; therefore RWG (i.e. early life increases in BMI) associated DNAm changes may also be consequential of the phenotype rather than causal.

Similar to Reed et al,. [51], we did not identify strong direct associations between childhood methylation at these CpGs and later BMI. They did however identify associations between a DNAm score for BMI and health outcomes where BMI is a risk factor, suggesting methylation may have more utility as a biomarker of BMI-related morbidity than as a predictor of BMI itself [51].

Strengths and limitations

A limitation of our study is the small sample size; however despite this, results were independently replicated at 1 of the CpGs. It will be important to replicate these findings in other cohorts with much larger sample sizes and a range of ages. This study has a number of strengths, principally the rare combination of detailed early life phenotypic and anthropometric data, as well as epigenetic data, as cohorts with longitudinal and epigenetic data of this nature are scarce. Other studies will undoubtedly also be limited by the lack of available early life weight measures, and future birth cohort studies should strive to collect these vital data. Whilst our analysis may have been underpowered, robust associations were still identified, although it is possible that other differentially methylated sites may have been missed.

Here we investigated RWG in the first year; however the entire childhood period could also be a critical period for growth and development of obesity [52]. There is the possibility that associations may be stronger in early life (before age 7), closer to the timing of the exposure, although this is not possible to test without multiple measures of DNAm throughout childhood.

There were differences in the associations identified in the models with and without adjustment for cell types. Whole blood represents a mixed cell population with varying proportions of white blood cells. Phenotypic variation in cell-type composition could confound analyses, or it could also represent an important physiological change in response to an exposure or disease, which may be related to the phenotype of interest. Obesity is an acknowledged chronic, inflammatory condition, and has been associated with inflammatory indicators including C-reactive protein [53] and white blood cell counts [54, 55]. When investigating biomarkers (related to an exposure) that are associated with an inflammatory disease outcome, to remove variation from cell counts could potentially disregard important loci.

We utilised SVA and ISVA to remove unwanted variation from confounders (which cannot always be adequately corrected for), whilst retaining differences due to the exposure of interest [56]. SVA finds sources of variation from the methylation data itself, and models these as linearly uncorrelated singular vectors (surrogate variables) which are then included as covariates in the regression model [33]. ISVA is a modified version of SVA where surrogate variables are deemed independent. In support of ISVA, known confounding factors such as age and batch are linearly uncorrelated statistically independent variables, therefore it would be appropriate to model these as independent variables. ISVA was shown to perform best at capturing a known specific biological signature when compared to other adjustment methods [34], although this may not hold true for all datasets [34, 57]. A thorough simulation study compared each of the common adjustment methods (Houseman’s reference-based method, RefFreeEWAS, SVA, ISVA, EWASher and RUV) and found no method performed perfectly for all parameters, but concluded SVA was most robust (and ‘safest’) [58]. In summary, as high-performing adjustment methods both were utilised in these analyses.

Conclusion

Our findings suggest that differential DNAm at 2 loci could be markers of early weight gain. At one CpG site, the highest levels of childhood methylation were in those who had RWG in early life and were subsequently overweight/obese in childhood or adolescence; therefore this site may have use as a biomarker of subsequent overweight/obesity in those who experience RWG. The EWAS identified 2 potentially important candidate sites, which could be the focus of further investigation. Further work is required to determine whether these CpG sites are consistently, differentially methylated in different populations, time points, and ages.

References

  1. Johnson W, Li L, Kuh D, Hardy R. How has the age-related process of overweight or obesity development changed over time? co-ordinated analyses of individual participant data from five United Kingdom Birth Cohorts. PLOS Med. 2015;12(5):e1001828.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Barker DJP. The developmental origins of adult disease. J Am Coll Nutr. 2004;23(sup6):588S-S595.

    Article  CAS  PubMed  Google Scholar 

  3. Stettler N, Kumanyika SK, Katz SH, Zemel BS, Stallings VA. Rapid weight gain during infancy and obesity in young adulthood in a cohort of African Americans. Am J Clin Nutrit. 2003;77(6):1374–8.

    Article  CAS  PubMed  Google Scholar 

  4. Ong KK, Loos RJF. Rapid infancy weight gain and subsequent obesity: Systematic reviews and hopeful suggestions. Acta Pædiat. 2006;95(8):904–8.

    Article  PubMed  Google Scholar 

  5. Druet C, Stettler N, Sharp S, Simmons RK, Cooper C, Davey Smith G, et al. Prediction of childhood obesity by infancy weight gain: an individual-level meta-analysis. Paediatr Perinat Epidemiol. 2012;26(1):19–26.

    Article  PubMed  Google Scholar 

  6. Zheng M, Lamb KE, Grimes C, Laws R, Bolton K, Ong KK, et al. Rapid weight gain during infancy and subsequent adiposity: a systematic review and meta-analysis of evidence. Obes Rev. 2018a;19(3):321–32.

    Article  CAS  PubMed  Google Scholar 

  7. Mathers JC, McKay JA. Epigenetics - potential contribution to fetal programming. Adv Exp Med Biol. 2009;646:119–23.

    Article  PubMed  Google Scholar 

  8. Suzuki MM, Bird A. DNA methylation landscapes: provocative insights from epigenomics. Nat Rev Genet. 2008;9(6):465–76.

    Article  CAS  PubMed  Google Scholar 

  9. Rzehak P, Covic M, Saffery R, Reischl E, Wahl S, Grote V, et al. DNA-methylation and body composition in preschool children: epigenome-wide-analysis in the european childhood obesity project (CHOP)-study. Sci Rep. 2017;7(1):14349.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. van Dijk SJ, Peters TJ, Buckley M, Zhou J, Jones PA, Gibson RA, et al. DNA methylation in blood from neonatal screening cards and the association with BMI and insulin sensitivity in early childhood. Int J Obes (Lond). 2018;42(1):28–35.

    Article  CAS  Google Scholar 

  11. van Dijk SJ, Tellam RL, Morrison JL, Muhlhausler BS, Molloy PL. Recent developments on the role of epigenetics in obesity and metabolic disease. Clinical epigenetics. 2015;7:66-.

  12. Wahl S, Drong A, Lehne B, Loh M, Scott WR, Kunze S, et al. Epigenome-wide association study of body mass index, and the adverse outcomes of adiposity. Nature. 2017;541(7635):81–6.

    Article  CAS  PubMed  Google Scholar 

  13. Rakyan VK, Down TA, Balding DJ, Beck S. Epigenome-wide association studies for common human diseases. Nat Rev Genet. 2011;12(8):529–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Flanagan JM. Epigenome-wide association studies (EWAS): past, present, and future. Methods Mol Biol (Clifton, NJ). 2015;1238:51–63.

    Article  Google Scholar 

  15. Küpers LK, Monnereau C, Sharp GC, Yousefi P, Salas LA, Ghantous A, et al. Meta-analysis of epigenome-wide association studies in neonates reveals widespread differential DNA methylation associated with birthweight. Nat Commun. 2019;10(1):1893.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Sharp GC, Salas LA, Monnereau C, Allard C, Yousefi P, Everson TM, et al. Maternal BMI at the start of pregnancy and offspring epigenome-wide DNA methylation: findings from the pregnancy and childhood epigenetics (PACE) consortium. Hum Mol Genet. 2017;26(20):4067–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wright C, Matthews J, Waterston A, Aynsley-Green A. What is a normal rate of weight gain in infancy? Acta Paediatr. 1994;83(4):351–6.

    Article  CAS  PubMed  Google Scholar 

  18. Wright CM, Cox KM, Sherriff A, Franco-Villoria M, Pearce MS, Adamson AJ. To what extent do weight gain and eating avidity during infancy predict later adiposity? Public Health Nutrition. 2012;15(4):656–62.

    Article  PubMed  Google Scholar 

  19. 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.

    Article  PubMed  Google Scholar 

  20. 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.

    Article  PubMed  Google Scholar 

  21. 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.

    Article  PubMed  Google Scholar 

  22. Cole TJ, Freeman JV, Preece MA. Body mass index reference curves for the UK, 1990. Arch Dis Child. 1995;73(1):25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Yu ZB, Han SP, Zhu GZ, Zhu C, Wang XJ, Cao XG, et al. Birth weight and subsequent risk of obesity: a systematic review and meta-analysis. Obes Rev. 2011;12(7):525–42.

    Article  CAS  PubMed  Google Scholar 

  24. Simpkin AJ, Suderman M, Gaunt TR, Lyttleton O, McArdle WL, Ring SM, et al. Longitudinal analysis of DNA methylation associated with birth weight and gestational age. Hum Mol Genet. 2015;24(13):3752–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Vidmar S, Carlin J, Hesketh K, Cole T. Standardizing anthropometric measures in children and adolescents with new functions for egen. Stata J. 2004;4(1):50–5.

    Article  Google Scholar 

  26. SACN. RCPCH Expert Group. Consideration of issues around the use of BMI centile thresholds for defining underweight, overweight and obesity in children aged 2–18 years in the UK. 2012.

  27. Touleimat N, Tost J. Complete pipeline for Infinium® Human Methylation 450K BeadChip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics. 2012;4(3):325–41.

    Article  CAS  PubMed  Google Scholar 

  28. 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 Bioinform. 2012;13:86.

    Article  Google Scholar 

  29. Chen Y-A, 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Jaffe AE, Irizarry RA. Accounting for cellular heterogeneity is critical in epigenome-wide association studies. Genome Biol. 2014;15(2):R31.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc: Ser B (Methodol). 1995;57(1):289–300.

    Google Scholar 

  32. Min J, Hemani G, Davey Smith G, Relton CL, Suderman M. Meffil: efficient normalisation and analysis of very large DNA methylation samples. bioRxiv. 2017.

  33. Leek JT, Storey JD. Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genet. 2007;3(9):e161.

    Article  PubMed Central  CAS  Google Scholar 

  34. Teschendorff AE, Zhuang J, Widschwendter M. Independent surrogate variable analysis to deconvolve confounding factors in large-scale microarray profiling studies. Bioinformatics (Oxford, England). 2011;27(11):1496–505.

    CAS  Google Scholar 

  35. Jones PA, Baylin SB. The fundamental role of epigenetic events in cancer. Nat Rev Genet. 2002;3(6):415.

    Article  CAS  PubMed  Google Scholar 

  36. Robinson MD, Kahraman A, Law CW, Lindsay H, Nowicka M, Weber LM, et al. Statistical methods for detecting differentially methylated loci and regions. Front Genet. 2014;5:324.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Peters TJ, Buckley MJ, Statham AL, Pidsley R, Samaras K, V Lord R, et al. De novo identification of differentially methylated regions in the human genome. Epigenetics & Chromatin. 2015;8(1):6.

  38. Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, et al. The human genome browser at UCSC. Genome Res. 2002;12(6):996–1006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Inskip HM, Godfrey KM, Robinson SM, Law CM, Barker DJ, Cooper C. Cohort profile: the Southampton women’s survey. Int J Epidemiol. 2006;35(1):42–8.

    Article  PubMed  Google Scholar 

  40. Turan N, Ghalwash MF, Katari S, Coutifaris C, Obradovic Z, Sapienza C. DNA methylation differences at growth related genes correlate with birth weight: a molecular signature linked to developmental origins of adult disease? Bmc Medical Genomics. 2012;5.

  41. Zheng M, Lamb KE, Grimes C, Laws R, Bolton K, Ong KK, et al. Rapid weight gain during infancy and subsequent adiposity: a systematic review and meta-analysis of evidence. Obes Rev. 2018b;19(3):321–32.

    Article  CAS  PubMed  Google Scholar 

  42. Scolnick DM, Halazonetis TD. Chfr defines a mitotic stress checkpoint that delays entry into metaphase. Nature. 2000;406(6794):430.

    Article  CAS  PubMed  Google Scholar 

  43. Sanbhnani S, Yeong FM. CHFR: a key checkpoint component implicated in a wide range of cancers. Cell Mol Life Sci CMLS. 2012;69(10):1669–87.

    Article  CAS  PubMed  Google Scholar 

  44. Derks S, Cleven AH, Melotte V, Smits KM, Brandes JC, Azad N, et al. Emerging evidence for CHFR as a cancer biomarker: from tumor biology to precision medicine. Cancer Metastasis Rev. 2014;33(1):161–71.

    CAS  PubMed  Google Scholar 

  45. Oh YM, Kwon YE, Kim JM, Bae SJ, Lee BK, Yoo SJ, et al. Chfr is linked to tumour metastasis through the downregulation of HDAC1. Nat Cell Biol. 2009;11:295.

    Article  CAS  PubMed  Google Scholar 

  46. Kisiel JB, Dukek BA, Ghoz HM, Yab TC, Berger CK, et al. Hepatocellular carcinoma detection by plasma methylated DNA: discovery, phase I pilot, and phase II clinical validation. Hepatology (Baltimore, Md). 2019;69(3):1180–92.

    CAS  PubMed Central  Google Scholar 

  47. Majumder S, Taylor WR, Foote PH, Berger CK, Wu CW, Yab TC, et al. Mo1370 - pancreatic cancer detection by plasma assay of novel methylated Dna markers: a case-control study. Gastroenterology. 2019;156(6):754.

    Article  Google Scholar 

  48. Gore L, Triche TJ, Farrar JE, Wai D, Legendre C, Gooden GC, et al. A multicenter, randomized study of decitabine as epigenetic priming with induction chemotherapy in children with AML. Clin Epigenet. 2017;9(1):108.

    Article  CAS  Google Scholar 

  49. Froy O. Metabolism and circadian rhythms—implications for obesity. Endocr Rev. 2010;31(1):1–24.

    Article  CAS  PubMed  Google Scholar 

  50. Woo Baidal JA, Locks LM, Cheng ER, Blake-Lamb TL, Perkins ME, Taveras EM. Risk factors for childhood obesity in the First 1,000 days: a systematic review. Am J Prev Med. 2016;50(6):761–79.

    Article  PubMed  Google Scholar 

  51. Reed ZE, Suderman MJ, Relton CL, Davis OSP, Hemani G. The association of DNA methylation with body mass index: distinguishing between predictors and biomarkers. Clin Epigenet. 2020;12(1):50.

    Article  CAS  Google Scholar 

  52. Cole T. Children grow and horses race: is the adiposity rebound a critical period for later obesity? BMC Pediatr. 2004;4(1):6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Visser M, Bouter LM, McQuillan GM, Wener MH, Harris TB. Elevated C-reactive protein levels in overweight and obese adults. JAMA. 1999;282(22):2131–5.

    Article  CAS  PubMed  Google Scholar 

  54. Farhangi MA, Keshavarz S-A, Eshraghian M, Ostadrahimi A, Saboor-Yaraghi A-A. White blood cell count in women: relation to inflammatory biomarkers, haematological profiles, visceral adiposity, and other cardiovascular risk factors. J Health Popul Nutrit. 2013;31(1):58–64.

    Google Scholar 

  55. Bastard J-P, Maachi M, Lagathu C, Kim MJ, Caron M, Vidal H, et al. Recent advances in the relationship between obesity, inflammation, and insulin resistance. Eur Cytokine Netw. 2006;17(1):4–12.

    CAS  PubMed  Google Scholar 

  56. Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Gayther SA, Apostolidou S, et al. An epigenetic signature in peripheral blood predicts active ovarian cancer. PLoS ONE. 2009;4(12):e8274.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  57. Teschendorff AE, Zheng SC. Cell-type deconvolution in epigenome-wide association studies: a review and recommendations. Epigenomics. 2017;9(5):757–68.

    Article  CAS  PubMed  Google Scholar 

  58. McGregor K, Bernatsky S, Colmegna I, Hudson M, Pastinen T, Labbe A, et al. An evaluation of methods correcting for cell-type heterogeneity in DNA methylation studies. Genome Biol. 2016;17(1):84.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  59. 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):1–17.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We are extremely grateful to all the participants who took part in the study, the midwives for their help in recruiting them, and the whole ALSPAC team, cohort administrators and all of the entire ALSPAC team involved in data collection, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists and nurses. The authors are grateful to the women of Southampton and their offspring who gave their time to take part in the Southampton Women’s Survey and to the research nurses and other staff who collected and processed the data.

Funding

The UK Medical Research Council and Wellcome (Grant ref: 217065/Z/19/Z) and the University of Bristol provide core support for ALSPAC. Methylation data in the ALSPAC cohort were generated as part of the UK BBSRC funded (BB/I025751/1 and BB/I025263/1) Accessible Resource for Integrated Epigenomic Studies (ARIES, https://www.ariesepigenomics.org.uk). CLR is supported by the UK Medical Research Council Integrative Epidemiology Unit at the University of Bristol (Grant ref: MC_UU_00011/5). Full details of ALSPAC grant funding can be found on the ALSPAC website (https://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). The SWS is supported by grants from the Medical Research Council, National Institute for Health Research Southampton Biomedical Research Centre, University of Southampton and University. Hospital Southampton National Health Service Foundation Trust, and the European Union’s Seventh. Framework Programme (FP7/2007–2013), project EarlyNutrition (grant 289346). Study participants were drawn from a cohort study funded by the Medical Research Council and the Dunhill Medical Trust. This publication is the work of the authors, who will serve as guarantors for the contents of this paper. This research was supported by a BBSRC DTP studentship (Grant ref: BB/M011186/1) at Newcastle University. CLR is supported by the UK Medical Research Council Integrative Epidemiology Unit at the University of Bristol (Grant ref: MC_UU_00011/5). KAL, EA and RM are part of an academic consortium that has received research funding from Abbott Nutrition, Nestec, BenevolentAI Bio Ltd. and Danone. EA is supported by Diabetes UK (16/0005454). KAL is supported by NIHR Southampton Biomedical Research Centre (IS-BRC-1215–20004)), British Heart Foundation (RG/15/17/3174), and Diabetes UK (16/0005454). KMG is supported by the UK Medical Research Council (MC_UU_12011/4), the National Institute for Health Research (NIHR Senior Investigator (NF-SI-0515–10042) and NIHR Southampton Biomedical Research Centre (IS-BRC-1215–20004)), the European Union (Erasmus + Programme ImpENSA 598488-EPP-1–2018-1-DE-EPPKA2-CBHE-JP), British Heart Foundation (RG/15/17/3174), the US National Institute On Aging of the National Institutes of Health (Award No. U24AG047867) and the UK ESRC and BBSRC (Award No. ES/M00919X/1). KMG has received reimbursement for speaking at conferences sponsored by companies selling nutritional products, and is part of an academic consortium that has received research funding from Abbott Nutrition, Nestec, BenevolentAI Bio Ltd. and Danone. MAH is supported by the British Heart Foundation.

Author information

Authors and Affiliations

Authors

Contributions

NR was responsible for the data analysis, prepared the tables/figures, drafted the manuscript, and reviewed and edited the manuscript. JAM, HB, VA and MSP were responsible for the concept and design, and critical revision of the manuscript. CLR contributed to the early appraisal of the methods. KMG, KAL, MAH, RM are responsible for the SWS study design, concept and/or data collection, EA performed the statistical replication analysis in the SWS. All authors read and approved the final manuscript.

Corresponding author

Correspondence to N. Robinson.

Ethics declarations

Ethics approval and consent to participate

Informed consent for the use of ALSPAC data collected via questionnaires and clinics was obtained from participants following the recommendations of the ALSPAC Ethics and Law Committee at the time. Follow-up of the children and sample collection/analysis was carried out under Institutional Review Board approval (Southampton and SW Hampshire Local Research Ethics Committee) with written informed consent obtained from parents or guardians.

Consent for publication

Manuscript was approved by the ALSPAC Executive prior to journal submission.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Southampton's Women’s Survey supplementary methods.

Additional file 2.

Supplementary tables and figures.

Additional file 3.

Top CpG hits by model for DNA methylation at ages 7 and 17.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Robinson, N., Brown, H., Antoun, E. et al. Childhood DNA methylation as a marker of early life rapid weight gain and subsequent overweight. Clin Epigenet 13, 8 (2021). https://doi.org/10.1186/s13148-020-00952-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-020-00952-z

Keywords