Skip to main content

Duration of breastfeeding is associated with leptin (LEP) DNA methylation profiles and BMI in 10-year-old children



Breastfeeding is protective against many long-term diseases, yet the mechanisms involved are unknown. Leptin gene (LEP) is reported to be associated with body mass index (BMI). On the other hand, breastfeeding duration has been found to be associated with DNA methylation (DNAm) of the LEP gene. Therefore, epigenetic regulation of LEP may represent the mechanism underlying the protective effect of breastfeeding duration against obesity.


In the Isle of Wight Birth Cohort, peripheral blood DNAm at 23 cytosine-phosphate-guanine sites (CpGs) in the LEP locus in 10-year-old (n = 297) samples and 16 CpGs in 18-year-old (n = 305) samples, were generated using the Illumina Infinium MethylationEPIC and HumanMethylation450 Beadchips respectively and tested for association with breastfeeding duration (total and exclusive) using linear regression. To explore the association between breastfeeding durations and genome-wide DNAm, epigenome-wide association studies (EWASs) and differential methylation region (DMR) analyses were performed. BMI trajectories spanning the first 18 years of life were used as the outcome to test the association with breastfeeding duration (exposure) using multi-nominal logistic regression. Mediation analysis was performed for significant CpG sites.


Both total and exclusive breastfeeding duration were associated with DNAm at four LEP CpG sites at 10 years (P value < 0.05), and not at 18 years. Though no association was observed between breastfeeding duration and genome-wide DNAm, DMR analyses identified five significant differentially methylated regions (Sidak adjusted P value < 0.05). Breastfeeding duration was also associated with the early transient overweight trajectory. Furthermore, DNAm of LEP was associated with this trajectory at one CpG site and early persistent obesity at another, though mediation analysis was not significant.


Breastfeeding duration is associated with LEP methylation at age 10 years and BMI trajectory. LEP DNAm is also significantly associated with BMI trajectories throughout childhood, though sample sizes were small. However, mediation analysis did not demonstrate that DNAm of LEP explained the protective effect of breastfeeding against childhood obesity.


The multiple health benefits of breastfeeding are well recognised, with the World Health Organization (WHO) recommending that infants should be exclusively breastfed up until 6 months and then breastfed alongside complementary food until 2 years of age [1], though globally, a high proportion of infants are not breastfed according to international guidelines [2, 3].

Infants who are not breastfed are at increased risk of all-cause and infection-related mortality [4]. For infants with higher intensity of breastfeeding, there is a reduction in morbidities such as respiratory infection and diarrhea in their first five years of life [5]. Long-term, adequate breastfeeding has been associated with a range of phenotypes including a lowered risk of type 2 diabetes [6], improved performance in intelligence tests [7], and reduced risk of childhood leukemia [8]. Conversely, associations between breastfeeding and later disease risk are not always as clear cut, such as in the context of allergic disease [9,10,11]. Methodical disparities including defining exposure and adjustment for confounding factors make the comparison between studies difficult [12]. Similarly, conflicting results on the relationship between breastfeeding and obesity have been attributed to heterogeneity between studies, residual confounding, and insufficient statistical power [13]. However, several studies and meta-analyses have recently reported a significant association between breastfeeding duration and childhood obesity [6, 14, 15]. Based on current evidence, the WHO reports that breastfeeding appears to reduce the risk of being overweight and obese in childhood [16].

Early life nutrition, and breastfeeding in particular, up to 2 years of life, are considered to be critical for programming of long-term health. The composition of human breastmilk, as a source of nutrition and bioactive factors, is crucial in understanding the mechanism behind its beneficial effect [17]. For example, breastmilk is high in long-chain polyunsaturated fatty acids (LCPUFAs), which in adults are associated with a lower risk of metabolic syndrome and cardiovascular disease [18], though evidence of their benefit in children is less clear [19, 20]. Moreover, maternal antibodies in human milk are taken up by an infant’s intestinal dendritic cells to support their immature acquired immunity [17]. Exactly how these metabolites and antibodies convey the long-term protective effects of breastmilk needs further investigation [21], though recent evidence suggests that some may be directly related to epigenetic changes [22] and infants’ gut microbiota [23].

A variety of nutritional sources are thought to induce epigenetic changes, including changes in DNA methylation (DNAm) [24]. In particular, early life nutritional exposures, such as maternal protein restriction in pregnancy and methyl donor nutrients including folate, methionine, and some B-vitamins, may leave long-lasting changes to DNAm in humans [25, 26]. Furthermore, components of human breastmilk are thought to induce epigenetic changes such as to LCPUFAs [27], and lactoferrin, which by suppressing the NF-κB signaling pathway may explain why breastmilk is protective against neonatal necrotizing enterocolitis [22].

However, to date, few studies have investigated the association between breastfeeding duration and DNAm levels [28, 29]. In one of these studies, Obermann-Borst and colleagues demonstrated that breastfeeding duration was significantly associated with a reduction in methylation in whole blood DNA of the leptin gene (LEP) in young children (mean age 17 months) [28]. LEP is expressed in adipocytes and encodes the hormone leptin, important in the regulation of energy intake by inducing early satiety. However, defective production of LEP is associated with obesity [30], and LEPR/LEPROT DNAm has been demonstrated to interact with genotype to influence both leptin levels and BMI [31]. This suggests that breastfeeding duration, associated with postnatal DNAm, may contribute to childhood obesity. The aim of this study was to replicate and extend the study by Obermann-Borst et al., in the Isle of Wight Birth Cohort (IOWBC) to test the association between breastfeeding and DNAm of cytosine-phosphate-guanines (CpGs) in LEP gene in peripheral blood at later ages (10 and 18 years) and further explore the association between breastfeeding duration and body mass index (BMI), where DNAm may act as a mediator.


Isle of Wight Birth Cohort

The IOWBC (second generation also known as IoW F1) was recruited between 1989 and 1990 (n = 1536) [32, 33]. The parents (first generation, IoW F0) of all infants born over this period were contacted at birth, and subsequently, 95% of infants (n = 1456) were enrolled following informed consent and exclusion. Follow-ups were conducted at 1, 2, 4, 10, and 18 years, where information was gathered about the participants, such as infant nutrition and breastfeeding practice [34]. Peripheral blood samples were collected from participants at age 10 and 18 years and DNAm measured for 297 and 305 samples, at ages 10 and 18 respectively, of whom 162 participants were matched between the two groups. Descriptive statistics of the sub-cohorts for exposure and covariates are presented in Table 1.

Table 1 Characteristics of Isle of Wight Birth Cohort participants at 10 and 18 years.

Categorization of breastfeeding duration

Two categories, exclusive and total breastfeeding durations, were used to define the duration of breastfeeding for each participant. The exclusive breastfeeding duration was defined as the length of time (in weeks) a child was breastfed, until formula feed and/or solid foods were introduced. The total breastfeeding duration, equates to the total number of weeks a mother breastfed her child irrespective of the introduction of formula feed and/or solid foods.

DNA extraction and arraying

DNA was extracted from peripheral blood samples using a standard salting out procedure. Approximately 1 μg of DNA was bisulfite-treated using the EZ 96-DNA methylation kit (Zymo Research, Irvine, CA, USA) using the manufacturer’s standard protocol. The Infinium MethylationEPIC BeadChips (age 10-year samples) and Infinium HumanMethylation450 (age 18-year samples) from Illumina (Illumina, San Diego, CA, USA) were used to obtain methylation levels following the manufacturer’s standard protocol. Methylation 450k and EPIC data (β values) were pre-processed for quality control, following the CPACOR pipeline [35], and batch effect was corrected using ComBat [36].

Candidate Leptin CpG selection

In Obermann-Borst et al. study, seven CpG sites in the LEP gene were analyzed. The UCSC Batch Coordinate Conversion tool (LiftOver, was used to convert genome coordinates from the hg18 assembly, used by Obermann-Borst et al., to the human genome hg19 assembly [37]. Four of these seven CpG sites were present in our analysis. For all CpGs annotated near LEP gene, for the 18-year samples, a total of 16 CpG sites were analyzed in LEP, whereas 23 CpG sites in LEP were analyzed in the 10-year samples, where the more comprehensive EPIC beadchips were used (see Additional file 1: Table S1).

Confounding factors

The same covariates used by Obermann-Borst et al. were used in the analysis, namely, child sex, birth weight, gestational age, serum leptin (ng/mL at either 10 or 18 years), BMI (at either 10 or 18 years), maternal socioeconomic status (defined using maternal socioeconomic cluster information: high, low, low-low, low-mid, and mid, using household income, number of rooms, and maternal education), and maternal smoking at birth. However, unlike the original study, birth weight, and gestational age were kept as separate variables and z-scores were not considered to be appropriate [38]. Parity and maternal age were added as covariates for these analyses, having been previously found to be associated with breastfeeding duration [39]. Additionally, breastfeeding duration was used as a continuous variable rather than categorical. Cell proportions (CD8T, CD4T, NK, B Cells, monocytes, granulocytes) were estimated using the minfi package [40] and cell composition coefficients were derived using the Houseman method [41].

BMI trajectories

Four BMI trajectories (normal, early persistent obesity, delayed overweight, and early transient overweight) covering the period from 1 to 18 years of age, as proposed by Ziyab et al. [42], were utilized to assess the association between breastfeeding duration and BMI trajectories. These four trajectories were categorized using age-specific BMI thresholds for ‘obesity’ and ‘overweight,’ as defined by the International Obesity Taskforce for ages 4, 10 and 18 years [43], and WHO standards for age 1 year [44]. They were defined in 1240 participants of the IOWBC and the trajectories for 10 and 18-year samples are described in Table 1.

Statistical analyses

Intra-class correlations (ICC) were estimated among 14 common LEP CpGs in 162 matched samples (common between 10 and 18 years) to assess the DNAm stability over time. Four models were run exploring the association between exclusive breastfeeding duration or total breastfeeding duration and LEP DNAm at either 10 or 18 years using linear regression adjusted for covariates and cell types.

Further, we performed an epigenome-wide association study (EWAS) to identify differentially methylated positions (DMPs) that are associated with breastfeeding durations (both exclusive and total), using linear regression adjusted for confounding factors and cell types, in both 10- and 18-year samples as described for the LEP analyses. Multiple hypothesis testing was accounted for by controlling the false discovery rate (FDR), using Benjamini and Hochberg’s algorithm [45]. CpGs with FDR corrected P value < 0.05 were considered statistically significant. To identify differentially methylated regions (DMRs), composed of multiple signals across individual CpG positions, comb-p [46] (Python v2.7) was used. Comb-p identifies regions enriched with low unadjusted P values (here, P value < 0.05) from the EWAS analysis, based on the probe location. For each region, the comb-p algorithm adjusts the CpG P values for auto-correction between probes by using the Stouffer-Liptak-Kechris (slk) correction with multiple testing adjustment using a one-step Sidak correction method. Regions with least two CpG probes within 200 base pairs, having a Sidak-corrected P value < 0.05, were considered statistically significant.

The association between breastfeeding duration (exclusive and total) and BMI trajectories was tested in the subset of children, matched with 10-year old samples, using multi-nominal logistic regression with ‘normal’ BMI trajectory as reference. Furthermore, the association between BMI trajectory and LEP DNAm was investigated with only those CpGs found to be significant with breastfeeding duration. P value < 0.05 was used as the significance level. All analyses were performed using R (version 3.3.2).

Mediation analysis

Model-based causal mediation analysis, using the R-package “mediation” [47], was performed to explore whether the association between breastfeeding duration and BMI trajectory is mediated by DNAm (Fig. 1). Only the LEP CpGs significantly associated with BMI trajectory were utilized for the mediation analysis as the mediator. Firstly, we estimated the effect of breastfeeding duration (exposure) on LEP DNAm (mediator) and then, we estimated the combined effect of breastfeeding duration (exposure) and LEP DNAm (mediator) on the BMI trajectory (outcome). Both of these effects were adjusted for the covariates mentioned above. Bootstrapping was used to test significance level (P value < 0.05).

Fig. 1
figure 1

A mediation model for the association between breastfeeding duration (total and exclusive) and childhood obesity. βM represents the effect estimate for breastfeeding duration on DNA methylation (CpG ~ breastfeeding duration + covariates); βY represents the effect estimate for CpG on childhood obesity (BMI trajectory ~ CpG + covariates; βc represents the direct effect (no mediation) estimate for breastfeeding duration on childhood obesity (BMI trajectory ~ breastfeeding duration + covariates; βc′ represents the total effect estimate on childhood obesity (BMI trajectory = breastfeeding duration + covariates + CpG)


Stability of DNA methylation

Intra-class correlation shows that 64.3% of the common CpGs (9 out of 14) were fairly moderately stable (ICC between 0.05 to 0.75) over 8 years (Additional file 1: Table S2). A moderate-good degree of stability was found in three CpGs between the two measurements. The highest ICC observed (in cg19594666) was 0.69 with a 95% confidence interval from 0.58 to 0.77 (F(161,161) = 3.22, P value = 2.79 × 10−13).

Breastfeeding duration and LEP DNA methylation

At 10 years, the total breastfeeding duration was significantly associated with DNAm at four CpG sites in LEP (P value < 0.05, Table 2). Two of the CpGs (cg03084214 and cg23753947) were available in 18-year samples, which showed significant differences with 10-year samples (Additional file 1: Figure S1). However, for the 18-year samples, there were no significant associations seen between either total breastfeeding or exclusive breastfeeding duration and DNAm at any of the 16 LEP CpG sites interrogated (Additional file 1: Tables S3 and S4).

Table 2 Association between total breastfeeding duration and LEP DNAm at 10 years

Similarly, the duration that a child was exclusively breastfed was associated with significant changes (P value < 0.05) in DNAm at four CpG sites, in samples taken at 10 years (Table 3). Lower methylation of one site (cg23381058) was significantly associated with both exclusive and total breastfeeding duration. Another CpG (cg18603538) was available in 18-year samples, which showed no significant differences with 10-year samples (Additional file 1: Figure S1).

Table 3 Association between exclusive breastfeeding duration and LEP DNAm at 10 years

Genome-wide DMPs and DMRs identifications

The EWAS analyses for both exclusive and total breastfeeding did not identify any significant DMPs (FDR < 0.05) at either age 10 or 18 years (Additional file 1: Figures 2–5). Analysis of differentially methylated regions using comb-p identified five significant DMRs (Sidak-corrected P values 9.6 × 10−9−0.005) (Tables 4 and 5). We identified two DMRs for 10-year samples and one for 18-year samples in the case of total breastfeeding duration. For exclusive breastfeeding duration, we identified two significant DMRs for 18-year samples and none for the 10-year samples. There was no overlap between the DMRs identified for exclusive and total breastfeeding.

Table 4 Statistically significant differentially methylated regions (DMRs) (Sidak P value < 0.05) for total breastfeeding duration
Table 5 Statistically significant differentially methylated regions (DMRs) (Sidak P value < 0.05) for exclusive breastfeeding duration

Breastfeeding duration and BMI trajectories

Total duration of breastfeeding was significantly associated with the early transient overweight BMI trajectory (estimate = − 0.02, P value = 0.002) (Table 6). However, there was no association with the early persistent obesity or delayed overweight trajectories. A similar pattern was seen when testing exclusive breastfeeding duration and obesity (Table 6), though the association between exclusive breastfeeding duration and early transient overweight narrowly missed our threshold for significance (estimate = − 0.03, P value = 0.05).

Table 6 Association between duration of total and exclusive breastfeeding and BMI trajectory

DNA methylation and BMI Trajectories

Of the four CpG sites where total breastfeeding duration was associated with DNAm, one site (cg23381058) also had a significant association between DNAm and the early transient overweight trajectory (P value = 0.034) (Table 7). Note that this CpG site was associated with both total and exclusive breastfeeding duration. Of the four sites where exclusive breastfeeding was associated with DNAm, methylation at another site (cg05091920) was significantly associated with the early persistent obesity trajectory (P value = 0.028). These results have been placed in Additional file 1: Table S7 along with other non-significant CpGs for all of the trajectories.

Table 7 Mediation analysis between cg23381058, total breastfeeding duration and BMI Trajectory 3 (Early transient overweight)

Table 7 Association between early transient overweight (BMI trajectory 3) and LEP DNAm at 10 years previously associated with breastfeeding duration (associated P values are shown)

Columns: Probe name = LEP CpG sites from Illumina Infinium MethylationEPIC Beadchip, which are significantly associated with either total or exclusive breastfeeding; Estimate = coefficient denoting one-unit increase in the methylation level for specific CpG is associated with the decrease in the log odds of being in the BMI trajectory 3; standard error = standard error of coefficient. DNA methylation analyzed against BMI trajectories with BMI trajectory 1 (normal) as a control. P value < 0.05 was considered significant and indicated by *

Mediation analysis

There were no statistically significant effects of cg23381058 methylation as a mediator between total breastfeeding duration and the early transient overweight trajectory (Table 8). Similarly, there was no association when looking at exclusive breastfeeding duration (Table 8). In the case of total breastfeeding duration, no significant direct effect (P value = 0.06) or total effect (P value = 0.07) on early transient overweight trajectory was observed.

Table 8 Mediation analysis between cg23381058, total breastfeeding duration, and BMI trajectory 3 (early transient overweight)


We examined the association between breastfeeding duration and DNAm on the LEP gene, finding a significant association at 10 years but not at 18 years. Exclusive breastfeeding duration was significantly associated with DNAm at four CpGs and total breastfeeding was also associated with DNAm at four sites, and among them, one site, cg23381058, was in common. Previous studies investigating the association between breastfeeding and DNAm, in animal and human models, have used total breastfeeding duration as their exposure [48]. However, there are concerns about the variability of breastfeeding within these parameters. For example, infants could be largely formula-fed, yet still classified as breastfeeding. Using an exclusive breastfeeding category, in comparison with total breastfeeding, gives more confidence over the exposure variable. The importance of this distinction is highlighted by ongoing research into impact of breastfeeding on infant outcomes, shaping health policy [1].

In addition, epigenome-wide association analysis was undertaken to identify other DMPs and DMRs associated with the breastfeeding duration (exclusive and total) and for both 10- and 18-year samples. To our knowledge, this is the first time such an analysis, investigating the association between breastfeeding and DNAm, has been conducted using an EWAS approach. In the case of DMPs, no genome-wide significant (FDR < 0.05) hits were found. This is likely due to lack of power and follow up studies in larger cohorts and/or with DNAm measurements earlier in childhood are warranted. We did however identify five significant DMRs (Sidak-corrected P values < 0.05) associated with breastfeeding duration. One of these genes, GAS5, was reported in a recent paper exploring whether long non-coding RNAs (lncRNAs), important in the development and gene expression, were encapsulated in extracellular vesicles within breastmilk [49]. Among others, GAS5 was detected and indeed highly correlated. Involved in metabolic functions and immune system programming, and as a significant epigenetic regulator, GAS5 is a potential candidate breast-milk compound, that could be important in the effect of breastfeeding. We are not aware of any previous research linking breastfeeding with the remaining DMRs; WBP1L, CYTH4, VSTM1, and FAM228A. Interestingly though, a recent study by Rzehak et al. [50], found an association between DNAm at cg14518658 near CYTH4 gene and body composition in pre-school children, adjacent to one of the significant DMRs (chr22:37678728-37678791) identified in this study.

The study by Obermann–Borst et al. demonstrated that breastfeeding is associated with DNAm in LEP gene at 17 months of age. Obermann-Borst et al. hypothesized that induced epigenetic regulation of leptin expression may be mechanistic in the protective effect of breastfeeding against obesity (Fig. 1). Our results support Obermann-Borst et al. findings, at a later stage in childhood (10 years) but not into young adulthood (18 years). However, though significant, these effect sizes were small and we saw varying direction of effect between probes. Stability of DNAm at cg23381058 over time could not be assessed as it was not present on both DNAm Bead Chip platforms. However, the majority of CpGs exhibited moderate stability over time. This change in DNAm may be a response to other environmental stimuli by 18 years, and indeed it is recognized that DNAm levels in an individual alter throughout childhood [51].

The next stage was to examine whether breastfeeding duration is associated with obesity in the IOWBC. The BMI trajectories, devised by Ziyab et al. [42], incorporate the different age-specific (BMI) thresholds of obesity proposed by the International Obesity Task Force and WHO, allowing analysis of the maintenance of this association throughout childhood. A significant association between total breastfeeding duration and the ‘early transient overweight’ BMI trajectory was seen in the subset of children, matched to 10-year-old samples. It should be noted that the sample size in this trajectory was considerably smaller (n = 33 at age 10) than in the normal reference category (n = 211). No association was seen between total breastfeeding duration and early persistent obesity or delayed overweight trajectory. Again the sample sizes in these trajectories were far smaller (n = 12 at age 10 in early persistent obesity group) than the reference trajectory; therefore, a significant effect of breastfeeding may not have been seen due to low power. Nonetheless, these results suggest breastfeeding duration may be involved in a short-term inverse association with childhood obesity. Another independent study examining the BMI trajectories of 276 infants, found that formula-fed infants gain weight more rapidly and disproportionately than breastfed infants [52]. This was attributed to increased lean mass, suggesting that delayed obesity associated with breastfeeding is not related to increased adiposity in infancy. This evidence contradicts our results, which indicates that breastfeeding is independently associated with both DNAm and obesity in childhood, but not long-term.

Multiple studies have demonstrated that DNAm at a number of loci is significantly associated with childhood obesity, for example in HIF3A gene [53], as well as leptin and adiponectin, where DNAm is also associated with low-density lipoprotein cholesterol levels [54, 55]. However, in this case, LEP DNAm levels have been reported to be negatively associated with BMI, contrasting with the hypothesis by Obermann-Borst et al., who suggested that hypo-methylation on the LEP gene may protect against obesity, i.e., reduce BMI. Indeed, we found a mixture of increased and decreased methylation indicating the relationship may not be simple. We also sought to test the association between DNAm and BMI trajectory in the IOWBC. DNAm at the CpG site cg23381058, which had shown a significant negative association with total and exclusive breastfeeding duration, was found to be negatively associated with early transient overweight, though mediation analysis of this CpG site was not significant (P < 0.05). This CpG site is located within an intron of LEP, in a transcription factor binding site. Conversely, an effect was seen at cg05091920, though here a decrease in the methylation level was associated with the decrease in the log odds of being in early persistent obesity trajectory.

We recognize that there are some limitations to this study. Firstly, we had no measurements of DNAm in infancy and thus our mediation analysis is based on the assumption that DNAm at age 10 years associated with breastfeeding duration is stable from infancy. Most importantly, it is difficult to establish a causal relationship between breastfeeding, DNAm, and obesity. It has been previously reported that changes in DNAm are a consequence of obesity [56,57,58]; however, these studies were largely conducted in adult populations, and inter-study comparison is difficult. Furthermore, breastfeeding may alter DNAm in a tissue-specific manner, and what effect this DNAm has on LEP expression as well as tissue specificity must be explored. Whilst this study and that of Obermann-Borst et al. used peripheral blood samples, and leptin hormone circulates through the blood, the LEP gene is known to be expressed in adipocytes, questioning the biological relevance. However, it has been demonstrated that epigenetic markers in blood accurately reflect those in adipose tissue [54, 55]. Peripheral blood remains an effective tissue from which to measure DNAm in childhood, though the investigation of alternative tissues is of interest.

The findings of this study require replication in other cohorts with childhood data and larger sample sizes to improve our understanding of how the relationship between breastfeeding, DNAm, and obesity varies throughout childhood. Cohorts with methylation data from cord blood and matched participants later in childhood could also help to establish causality. Epigenome-wide association studies have helped identify thousands of differentially methylated genes, associated with disease processes and/or environmental processes. This will be an important step to identify the importance of DNAm as a mechanism propagating the long-term protective effects of breastfeeding.


In line with previous studies, a significant association of breastfeeding duration with altered methylation at seven CpG sites on LEP gene was observed, and an association of breastfeeding duration with childhood obesity was confirmed. Future studies in larger population samples are warranted to further explore the potential of DNAm at LEP and other sites in the genome in mediating the effect of breastfeeding on childhood obesity.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. For access to the full Isle of Wight Cohort data please see:





DNA methylation


Isle of Wight Birth Cohort


Long-chain polyunsaturated fatty acids




World Health Organization


  1. WHO. Essential Nutrition Actions: improving maternal, newborn, infant and young child health and nutrition: World Health Organization; 2013.

  2. McAndrew F, Thompson J, Fellows L, Large A, Speed M, Renfrew MJ. Infant Feeding Survey 2010: Health and Social Care Information Centre; 2012.

  3. Victora CG, Bahl R, Barros AJD, França GVA, Horton S, Krasevec J, et al. Breastfeeding in the 21st century: epidemiology, mechanisms, and lifelong effect. Lancet. 2016;387(10017):475–90.

    Article  PubMed  Google Scholar 

  4. Sankar MJ, Sinha B, Chowdhury R, Bhandari N, Taneja S, Martines J, et al. Optimal breastfeeding practices and infant and child mortality: a systematic review and meta-analysis. Acta Paediatr. 2015;104:3–13.

    Article  PubMed  Google Scholar 

  5. Horta B, Victora CG. Short-term effects of breastfeeding: a systematic review of the benefits of breastfeeding on diarrhoea and pneumonia mortality: World Health Organization; 2013.

  6. Horta BL. Mola CL de, Victora CG. Long-term consequences of breastfeeding on cholesterol, obesity, systolic blood pressure and type 2 diabetes: a systematic review and meta-analysis. Acta Paediatr. 2015;104:30–7.

    Article  CAS  PubMed  Google Scholar 

  7. Horta BL, de MCL, Victora CG. Breastfeeding and intelligence: a systematic review and meta-analysis. Acta Paediatr. 2015 Nov;104:14–9.

    Article  PubMed  Google Scholar 

  8. Amitay EL, Keinan-Boker L. Breastfeeding and Childhood Leukemia Incidence. JAMA Pediatr. 2015 Jun;169(6):e151025.

    Article  PubMed  Google Scholar 

  9. Lodge CJ, Tan DJ, Lau MXZ, Dai X, Tham R, Lowe AJ, et al. Breastfeeding and asthma and allergies: a systematic review and meta-analysis. Acta Paediatr. 2015 Nov;104:38–53.

    Article  CAS  PubMed  Google Scholar 

  10. Yang YW, Tsai CL, Lu CY. Exclusive breastfeeding and incident atopic dermatitis in childhood: a systematic review and meta-analysis of prospective cohort studies. Br J Dermatol. 2009 Aug;161(2):373–83.

    Article  CAS  PubMed  Google Scholar 

  11. Bion V, Lockett GA, Soto-Ramírez N, Zhang H, Venter C, Karmaus W, et al. Evaluating the efficacy of breastfeeding guidelines on long-term outcomes for allergic disease. Allergy. 2016 May;71(5):661–70.

    Article  CAS  PubMed  Google Scholar 

  12. Matheson MC, Allen KJ, Tang MLK. Understanding the evidence for and against the role of breastfeeding in allergy prevention. Clin &amp; Exp Allergy. 2012 Jan;42(6):827–851.

  13. Beyerlein A, von Kries R. Breastfeeding and body composition in children: will there ever be conclusive empirical evidence for a protective effect against overweight? Am J Clin Nutr. 2011 Apr;94(suppl_6):1772S–5S.

    Article  CAS  PubMed  Google Scholar 

  14. Uwaezuoke SN, Eneh CI, Ndu IK. Relationship between exclusive breastfeeding and lower risk of childhood obesity: a narrative review of published evidence. Clin Med Insights Pediatr. 2017 Jan;11:117955651769019.

    Article  Google Scholar 

  15. Yan J, Liu L, Zhu Y, Huang G, Wang PP. The association between breastfeeding and childhood obesity: a meta-analysis. BMC Public Health. 2014 Dec;14(1).

  16. Robinson SM. Infant nutrition and lifelong health: current perspectives and future challenges. J Dev Orig Health Dis. 2015 Jun;6(05):384–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Ballard O, Morrow AL. Human milk composition. nutrients and bioactive factors. Pediatr Clin North Am. 2013 Feb;60(1):49–74.

    Article  PubMed  Google Scholar 

  18. Micha R, Mozaffarian D. Saturated fat and cardiometabolic risk factors, coronary heart disease, stroke, and diabetes: a fresh look at the evidence. Lipids. 2010 Mar;45(10):893–905.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Stroobant W, Braun KV, Jong JCK, Moll HA, Jaddoe VW, Brouwer IA, et al. Intake of different types of fatty acids in infancy is not associated with growth, adiposity, or cardiometabolic health up to 6 years of age. J Nutr. 2017 Jan:jn241018.

  20. Voortman T, van den HEH, Braun KVE, van den BM, Bramer WM, Chowdhurry R, et al. Effects of polyunsaturated fatty acid intake and status during pregnancy, lactation, and early childhood on cardiometabolic health: a systematic review. Prog Lipid Res. 2015 Jul;59:67–87.

    Article  CAS  PubMed  Google Scholar 

  21. Marincola FC, Dessì A, Corbu S, Reali A, Fanos V. Clinical impact of human breast milk metabolomics. Clin Chim Acta. 2015 Dec;451:103–6.

    Article  Google Scholar 

  22. Verduci E, Banderali G, Barberi S, Radaelli G, Lops A, Betti F, et al. Epigenetic Effects of Human Breast Milk. Nutrients. 2014 Apr;6(4):1711–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Pannaraj PS, Li F, Cerini C, Bender JM, Yang S, Rollie A, et al. Association between breast milk bacterial communities and establishment and development of the infant gut microbiome. JAMA Pediatr. 2017 Jul;171(7):647.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Canani RB, Di CM, Leone L, Bedogni G, Brambilla P, Cianfarani S, et al. Epigenetic mechanisms elicited by nutrition in early life. Nutr Res Rev. 2011 Oct;24(02):198–205.

    Article  CAS  PubMed  Google Scholar 

  25. Peter CJ, Fischer LK, Kundakovic M, Garg P, Jakovcevski M, Dincer A, et al. DNA methylation signatures of early childhood malnutrition associated with impairments in attention and cognition. Biol Psychiatry. 2016 Nov;80(10):765–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Choi S-W, Friso S. Epigenetics: a new bridge between nutrition and health. Adv Nutr. 2010 Oct;1(1):8–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Hoile SP, Clarke-Harris R, Huang R-C, Calder PC, Mori TA, Beilin LJ, et al. Supplementation with N-3 long-chain polyunsaturated fatty acids or olive oil in men and women with renal disease induces differential changes in the DNA methylation of FADS2 and ELOVL5 in peripheral blood mononuclear cells. Schunck W-H. PLoS One. 2014;9(10):e109896.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Obermann-Borst SA, Eilers PHC, Tobi EW, de JFH, Slagboom PE, Heijmans BT, et al. Duration of breastfeeding and gender are associated with methylation of the LEPTIN gene in very young children. Pediatr Res. 2013 Jun;74(3):344–9.

    Article  CAS  PubMed  Google Scholar 

  29. Rossnerova A, Tulupova E, Tabashidze N, Schmuczerova J, Dostal M, Rossner P, et al. Factors affecting the 27K DNA methylation pattern in asthmatic and healthy children from locations with various environments. Mutat Res. 2013 Feb;741–742:18–26.

    Article  PubMed  Google Scholar 

  30. Al Maskari MY, Alnaqdy AA. Correlation between serum leptin levels, body mass index and obesity in Omanis. Sultan Qaboos Univ Med J. 2006 Dec;6(2):27–31.

    PubMed  PubMed Central  Google Scholar 

  31. Yousefi M, Karmaus W, Zhang H, Ewart S, Arshad H, Holloway JW. The methylation of the LEPR/LEPROT genotype at the promoter and body regions influence concentrations of leptin in girls and BMI at age 18 years if their mother smoked during pregnancy. Int J Mol Epidemiol Genet. 2013;4(2):86–100.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Arshad SH, Holloway JW, Karmaus W, Zhang H, Ewart S, Mansfield L, et al. Cohort profile: the Isle Of Wight Whole Population Birth Cohort (IOWBC). Int J Epidemiol. 2018.

  33. Arshad SH, Karmaus W, Zhang H, Holloway JW. Multigenerational cohorts in patients with asthma and allergy. J Allergy Clin Immunol. 2017 Feb;139(2):415–21.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Arshad SH, Hide DW. Effect of environmental factors on the development of allergic disorders in infancy. J Allergy Clin Immunol. 1992;90(2):235–41.

    Article  CAS  PubMed  Google Scholar 

  35. Lehne B, Drong AW, Loh M, Zhang W, Scott WR, Tan S-T, et al. A coherent approach for analysis of the Illumina HumanMethylation450 BeadChip improves data quality and performance in epigenome-wide association studies. Genome Biol. 2015;16(1).

  36. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2006;8(1):118–27.

    Article  PubMed  Google Scholar 

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

  38. Delbaere I, Vansteelandt S, De BD, Verstraelen H, Gerris J, De SP, et al. Should we adjust for gestational age when analysing birth weights? The use of z-scores revisited. Hum Reprod. 2007 Jun;22(8):2080–3.

    Article  PubMed  Google Scholar 

  39. Kitano N, Nomura K, Kido M, Murakami K, Ohkubo T, Ueno M, et al. Combined effects of maternal age and parity on successful initiation of exclusive breastfeeding. Prev Med Reports. 2016;3:121–6.

    Article  Google Scholar 

  40. Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: a flexible and comprehensive bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30(10):1363–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Houseman E, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13(1):86.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Ziyab AH, Karmaus W, Kurukulaaratchy RJ, Zhang H, Arshad SH. Developmental trajectories of Body Mass Index from infancy to 18 years of age: prenatal determinants and health consequences. J Epidemiol Community Health. 2014;68(10):934–41.

    Article  PubMed  Google Scholar 

  43. Cole TJ. Establishing a standard definition for child overweight and obesity worldwide: international survey. BMJ. 2000;320(7244):1240–1240.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. WHO Child Growth Standards based on length/height, weight and age. Acta Paediatr (Oslo, Norw 1992) Suppl. 2006;450:76–85.

  45. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B [Internet]. 1995;57(1):289–300. Available from:

    Google Scholar 

  46. Pedersen BS, Schwartz DA, Yang IV, Kechris KJ. Comb-p: software for combining, analyzing, grouping and correcting spatially correlated P-values. Bioinformatics. 2012 Nov;28(22):2986–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Tingley D, Yamamoto T, Hirose K, Keele L, Imai K. mediation: RPackage for causal mediation analysis. J Stat Softw. 2014;59(5).

  48. Hartwig FP, Loret de Mola C, Davies NM, Victora CG, Relton CL. Breastfeeding effects on DNA methylation in the offspring: a systematic literature review. PLoS One [Internet]. 2017;12(3):e0173070 Available from:

    Article  PubMed  Google Scholar 

  49. Karlsson O, Rodosthenous RS, Jara C, Brennan KJ, Wright RO, Baccarelli AA, et al. Detection of long non-coding RNAs in human breastmilk extracellular vesicles: Implications for early child development. Epigenetics [Internet]. 2016;11(10):721–729. Available from:

    Article  PubMed  PubMed Central  Google Scholar 

  50. 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  Google Scholar 

  51. Gervin K, Andreassen BK, Hjorthaug HS, Carlsen KCL, akon CK-H, Undlien DE, et al. Intra-individual changes in DNA methylation not mediated by cell-type composition are correlated with aging during childhood. Clin Epigenetics. 2016;8(1).

  52. Bell KA, Wagner CL, Feldman HA, Shypailo RJ, Belfort MB. Associations of infant feeding with trajectories of body composition and growth. Am J Clin Nutr. 2017;106(2):491–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Wang S, Song J, Yang Y, Zhang Y, Wang H, Ma J. HIF3A DNA methylation is associated with childhood obesity and ALT. Costa M, editor. PLoS One. 2015;10(12):e0145944.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Houde A-A, Légaré C, Biron S, Lescelleur O, Biertho L, Marceau S, et al. Leptin and adiponectin DNA methylation levels in adipose tissues and blood cells are associated with BMI, waist girth and LDL-cholesterol levels in severely obese men and women. BMC Med Genet. 2015;16(1).

  55. Li R, Fein SB, Grummer-Strawn LM. Do infants fed from bottles lack self-regulation of milk intake compared with directly breastfed infants? Pediatrics. 2010;125(6):e1386–93.

    Article  PubMed  Google Scholar 

  56. Dick KJ, Nelson CP, Tsaprouni L, Sandling JK, Aïssi D, Wahl S, et al. DNA methylation and body-mass index: a genome-wide analysis. Lancet. 2014 Jun;383(9933):1990–8.

    Article  CAS  PubMed  Google Scholar 

  57. Ronn T, Volkov P, Gillberg L, Kokosar M, Perfilyev A, Jacobsen AL, et al. Impact of age, BMI and HbA1c levels on the genome-wide DNA methylation and mRNA expression patterns in human adipose tissue and identification of epigenetic biomarkers in blood. Hum Mol Genet. 2015;24(13):3792–813.

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

Download references


We would like to acknowledge the help of all the staff at The David Hide Asthma and Allergy Research Centre in undertaking the assessments of the Isle of Wight Birth Cohort and Nikki Graham for technical support. We would also like to acknowledge the help of the participants and their families who have helped us over the last two and a half decades.


SHA, WK, HZ, and JWH and ongoing studies in the IOWBC cohort are currently supported by the National Institutes of Health (R01AI121226, and R01AI091905, and 1R01HL132321). FIR and JWH acknowledge support from the Ageing Lungs in European Cohorts (ALEC) Study (, which has been funded by the European Union’s Horizon 2020 Research and Innovation programme under grant agreement no. 633212.

Author information

Authors and Affiliations



JWH and FIR conceived and designed the study. WBS analyzed the data and wrote the initial draft. VB, GAL, AHZ, NS, and NM analyzed phenotypes in the cohort. RJK and SHA collected the clinical data. JWH, SE, HZ, and WK designed and oversaw the epigenetic studies in the cohort. SHA is the Chief Investigator for the IOWBC. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Faisal I. Rezwan.

Ethics declarations

Ethics approval and consent to participate

Ethical approvals for both IOWBC was obtained from the relevant authorities as follows:

IOWBC–10-year follow-up: Isle of Wight Local Research Ethics Committee, Isle of Wight Health Authority, Protocol No. 18/98. Date: 20 July 1998. Eighteen-year follow-up: Isle of Wight, Portsmouth, and Southeast Hampshire Research Ethics Committee. REC reference number: 06/Q1701/34, Date: 21 June 2006.

Consent for publication

Not applicable.

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.

Additional file

Additional file 1:

Table S1. Annotation of LEP CpGs used in the analyses. Table S2. Stability of LEP CpGs (between 10 and 18 years).Table S3. Association between total breastfeeding duration and LEP DNAm (16 CpGs) at 18 years. Table S4. Association between exclusive breastfeeding duration and LEP DNAm (16 CpGs) at 18 years. Table S5. Association between total breastfeeding duration and LEP DNAm (23 CpGs) at 10 years. Table S6. Association between exclusive breastfeeding duration and LEP DNAm (23 CpGs) at 10 years. Table S7. Association between BMI trajectories and LEP DNAm (at 10 years) at CpGs previously associated with total or exclusive breastfeeding duration (associated P-values have been shown). Figure S1. Comparisons of methylation level (beta) of significant common CpGs in matched 162 samples (10 and 8 years). (DOCX 3070 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sherwood, W.B., Bion, V., Lockett, G.A. et al. Duration of breastfeeding is associated with leptin (LEP) DNA methylation profiles and BMI in 10-year-old children. Clin Epigenet 11, 128 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: