# Association between perinatal factors, genetic susceptibility to obesity and age at adiposity rebound in children of the EDEN mother-child cohort

Aminata Hallimat Cissé, Sandrine Lioret, Blandine de Lauzon-Guillain, Anne Forhan, Ken K Ong, Marie Aline Charles, Barbara Heude

► To cite this version:

Aminata Hallimat Cissé, Sandrine Lioret, Blandine de Lauzon-Guillain, Anne Forhan, Ken K Ong, et al. Association between perinatal factors, genetic susceptibility to obesity and age at adiposity rebound in children of the EDEN mother-child cohort. International Journal of Obesity, 2021, 10.1038/s41366-021-00847-w. hal-03228602v2

HAL Id: hal-03228602

https://hal.science/hal-03228602v2

Submitted on 26 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Association between perinatal factors, genetic susceptibility to obesity and age at adiposity rebound in children of the EDEN mother–child cohort

Aminata Hallimat Cissé<fnref n="1" /> · Sandrine Lioret<fnref n="1" /> · Blandine de Lauzon-Guillain<fnref n="1" /> · Anne Forhan<fnref n="1" /> · Ken K. Ong<fnref n="2" /> · Marie Aline Charles<fnref n="1" /> · Barbara Heude<fnref n="1" />

Received: 15 June 2020 / Revised: 21 April 2021 / Accepted: 27 April 2021
© The Author(s), under exclusive licence to Springer Nature Limited 2021. This article is published with open access

## Abstract

**Background** Early adiposity rebound (AR) has been associated with increased risk of overweight or obesity in adulthood. However, little is known about early predictors of age at AR. We aimed to study the role of perinatal factors and genetic susceptibility to obesity in the kinetics of AR.

**Methods** Body mass index (BMI) curves were modelled by using mixed-effects cubic models, and age at AR was estimated for 1415 children of the EDEN mother–child cohort study. A combined obesity risk-allele score was calculated from genotypes for 27 variants identified by genome-wide association studies of adult BMI. Perinatal factors of interest were maternal age at delivery, parental education, parental BMI, gestational weight gain, maternal smoking during pregnancy, and newborn characteristics (sex, prematurity, and birth weight). We used a hierarchical level approach with multivariable linear regression model to investigate the association between these factors, obesity risk-allele score, and age at AR.

**Results** A higher genetic susceptibility to obesity score was associated with an earlier age at AR. At the most distal level of the hierarchical model, maternal and paternal educational levels were positively associated with age at AR. Children born to parents with higher BMI were more likely to exhibit earlier age at AR. In addition, higher gestational weight gain was related to earlier age at AR. For children born small for gestational age, the average age at AR was 88 [±39] days lower than for children born appropriate for gestational age and 91 [±56] days lower than for children born large for gestational age.

**Conclusion** The timing of AR seems to be an early childhood manifestation of the genetic susceptibility to adult obesity. We further identified low birth weight and gestational weight gain as novel predictors of early AR, highlighting the role of the intrauterine environment in the kinetics of adiposity.

## Introduction

After birth, body mass index (BMI) increases to a peak during the first year of life and declines between age 4 and 6 years, then increases again. This BMI increase following the

minimum is called adiposity rebound (AR) [1–3]. This period has been characterized by cessation in the decline of fat mass index (defined as fat mass divided by height square) [4]. The timing of AR, independent of BMI itself, is considered a predictor of obesity in later childhood and adulthood [5–7]. More specifically, early AR is associated with increased risk of overweight or obesity and cardio-metabolic diseases in adulthood [8–10].

To provide a track to prevent the development of obesity, one must follow children’s BMI growth curve and better understand the determinants of age at AR. Furthermore, the DoHAD theory [11] suggests that the intrauterine and early postnatal environment are crucial to the child’s and adult’s healthy development [12], along with the existence of a nutritional and/or metabolic programming of the development of obesity [13]. As an example, excessive and insufficient gestational weight

<footnote>**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41366-021-00847-w.</footnote>

<footnote>Aminata Hallimat Cissé
aminata-hallimat.cisse@inserm.fr</footnote>

<footnote n="1">¹ Université de Paris, CRESS, INSERM, INRAE, F-75004, Paris, France</footnote>

<footnote n="2">² MRC Epidemiology Unit and Department of Paediatrics, Institute of Metabolic Science, University of Cambridge, Cambridge, UK</footnote>
gain was associated with increased risks of adverse childhood outcomes for the offspring such as low and high birth weight [14] and obesity in later life [15, 16]. Similarly, several studies have shown being small or large for gestational age (LGA) associated with increased risk of overweight and obesity [17, 18]. Although the link between intrauterine and early postnatal environment factors and age at AR has not been thoroughly investigated in the literature, a set of early prenatal and postnatal determinants of age at AR have been identified. A longitudinal study interested in pre-, perinatal, and parental determinants of the kinetics of the BMI curve among 1681 children showed that girls had earlier age at AR than boys [19]. Parental educational level, especially low maternal level of education (<7 years), is an additional important risk factor for early AR [2]. In their retrospective study, Péneau et al. examined the determinants of age at AR and reported that a larger maternal silhouette was associated with early age at AR in both boys and girls, while the paternal silhouette was related with early age at AR in girls only [20]. However, further evidence is needed to confirm these findings with measured parental weight data instead of reported silhouettes.

Genome-wide associations studies have highlighted single-nucleotide polymorphisms (SNPs) implicated in increased risk of obesity [21, 22] and the use of genetic susceptibility scores has been established. In a meta-analysis of four European cohorts [23], Elks et al. showed that association with this genetic susceptibility score emerges progressively associated in childhood with a child's weight, height, and BMI. However, to our knowledge, only one study has investigated the role of this genetic susceptibility to obesity in the kinetics of BMI at the time of AR [24].

In this context, the purpose of this study was to investigate the role of parental education, perinatal factors, and genetic susceptibility to obesity in the kinetics of AR within the French EDEN mother–child cohort. We hypothesize that the association between age at AR and perinatal factors is independent of genetic susceptibility to obesity and that age at AR could be programmed by a set of factors involved as early as in pregnancy.

# Material and methods

## Study population

The EDEN mother–child study is a prospective French cohort that aims to assess prenatal and postnatal determinants of child growth, health, and development. This cohort included 2002 pregnant women (before 24 weeks of amenorrhea) recruited from 2003 to 2006 in two university

hospitals in Nancy and Poitiers. Exclusion criteria were <18 years of age, multiple pregnancies, known diabetes before pregnancy, illiteracy, and intention to give birth elsewhere than these two French hospitals or to move outside the region within the next 3 years. Details of the study protocol have been published elsewhere [25].

The study was approved by the relevant ethical research Committee of Kremlin-Bicêtre Hospital and by the Data Protection Authority. Both parents gave their written consent.

## Data collection and variables generation

### Data collection

Data were collected from medical and obstetrical records or by interviews, face-to-face or self-administered questionnaires and clinical examinations. Biological samples, anthropometric, and clinical parameters measurements were collected during clinical examinations at birth, 1, 3, and 5 years. In all children, weight and height from birth to early adolescence were obtained during these visits or from health records. Children had a mean of 10 weight and height measurements (interquartile range 6–14 and 5–13, respectively) from birth to 13 years. Maternal smoking status during pregnancy (yes/no); parental characteristics such as maternal age (years), maternal and paternal educational level and parental anthropometric measurements were collected during pregnancy. Newborn characteristics as well as sex, gestational age, and birth weight were taken at birth [26]. Preterm birth (yes/no) was defined as babies born before 37 weeks of gestation. Birth weight customized z-scores were calculated according to Gardosi references taking into account physiological fetal (sex and gestational age) and maternal factors (weight, height, parity, origin) [27]. Newborns were classified, according to this customized approach, into three classes: small for gestational age (SGA, ≤10th percentile), appropriate for gestational age (AGA, >10th percentile to ≤90th percentile), and LGA (>90th percentile).

Maternal and paternal education level was considered as the number of years of study. Parental weight and height were collected by interview at inclusion. Parental anthropometric measurements were collected at inclusion were used to calculate parental BMI by dividing the weight (kilograms) by the square of height (meters). Gestational weight gain (kilograms) was collected from obstetric records.

### Genotyping and BMI obesity risk-allele score

DNA samples were extracted from cord-blood at birth. Genotyping for 27 single-nucleotide polymorphisms (SNPs) was performed by the Epidemiology Unit of the Medical
Research Council of Cambridge (iPLEX platform; Sequenom) [28]. Genotyping quality-control criteria (call rate, >95%; Hardy Weinberg balance, p > 0.01) were satisfactory for all variants. In the present study, we considered 27 of 32 SNPs associated with BMI in adults [21], of which also the 16 associated with childhood BMI in a meta-analysis of Elks et al. [23].

For each child, combined obesity risk-allele scores across the 27 SNP loci were calculated as the sum of risk alleles (0, 1, or 2 at each locus) associated with increased BMI. The score indicating genetic predisposition to obesity ranged from 11 to 32. Combined obesity risk-allele scores were available for 1322 children of the EDEN cohort and their mothers (N = 1678) and fathers (N = 1241) [29].

A weighted obesity risk-allele score was also computed by multiplying each risk-allele by the effect estimate on adult BMI as assessed in the genome-wide association study [21], for sensitivity analyses.

## Estimation of AR

Age at AR is defined as the last minimum (nadir) BMI before the continuous increase in BMI over time [1].

We calculated BMI from weight and height measurements. We excluded children with fewer than three BMI measurements from age 18 months to 13 years. Included children had 10 BMI values on average (interquartile range 6–14) from age 18 months to 13 years. Individual growth curves were obtained by using mixed-effects cubic models separately for girls and boys. The methods for growth modeling of age at AR were inspired by the model of Sovio et al. [30]. To this previously described model, we added random effects for the quadratic and cubic terms for age with unstructured covariance matrix, which provided the best fit of individual predicted curves. The model was then fitted for log-transformed BMI, separately for boys and girls, as follows:

$$ \log(\text{BMI}) = \beta_0 + \mu_0 + (\beta_1 + \mu_1)\text{Age} + (\beta_2 + \mu_2)\text{Age}^2 + (\beta_3 + \mu_3)\text{Age}^3 + \varepsilon $$

where BMI is expressed in kilograms per square meter and age in days; $\beta_0$, $\beta_1$, $\beta_2$, $\beta_3$, and $\beta_4$ are the fixed-effects terms; $\mu_0$, $\mu_1$, $\mu_2$, and $\mu_3$ are the individual-level random effects; and $\varepsilon$ is the residual.

For each participant, we estimated age at AR by the first and second derivatives of individual BMI functions. Results of this modeling process are illustrated for a random sample of ten participants in Supplementary Fig. S1.

## Sample selection

Among the 1907 newborns included in the EDEN cohort, we excluded 326 children without BMI measurements after

<figure><img src="image_3.png" /><figcaption>Fig. 1 Flow chart of the population included in the study.</figcaption></figure>

18 months and a further 126 who had fewer than three BMI measurements between age 18 months and 13 years. We also excluded ten children for whom AR could not be assessed because their BMI curve did not display any AR (see Supplementary Fig. S2). Therefore, the final sample included 1415 children; 69.8% (988) had available obesity risk-allele scores (Fig. 1).

## Statistical analysis

Main characteristics of study population were described with mean ± SD or percentage (N) before any imputation. We tested differences between included and excluded populations with chi-square tests for categorical variables and Student t for continuous variables. We also described and compared covariates according to each quintile of age at AR by chi-square tests and ANOVA. The associations between perinatal factors, obesity risk-allele score, and age at AR were investigated by linear multiple regression models.

By using the predictors identified in the literature, we developed a hierarchical level approach to answer our research question (Supplementary Fig. S3), which assumes that proximal variables have more direct effect than distal variables: the relevance of a factor is determined in the model in which the variable was first entered, without considering its performance in other models. This approach allows us to prevent that intermediate variables affect the association between distal variables and the dependent variable.

Therefore, we tested four successive nested models. Model A included the most distal-level variables
corresponding to demographic and socioeconomic factors (maternal age at delivery, maternal and paternal educational level), along with recruitment center; model B further included factors related to obesity risk such as child's obesity risk-allele score and maternal and paternal BMI; model C included also variables corresponding to the intrauterine and prenatal environment (gestational weight gain and maternal smoking during pregnancy); model D corresponded to the most proximal level, with the inclusion of newborn characteristics (sex, prematurity, and birth weight z-score). Birth weight z-score was introduced with a squared term to test for a potential quadratic shape of the association.

To assign a value to missing data for covariates (Supplementary Table S1), we first performed multiple imputation analyses, assuming that data were missing at random. Using the fully conditional specification method, we generated 40 independent datasets. Then, we computed pooled effect estimates. Parental obesity risk-allele scores were used to impute children's scores when these were missing. With multiple imputation, the main analysis involved a sample of 1415 children.

We conducted a first sensitivity analysis running our multivariable models with data for the 866 children with complete cases. Most of the decline in sample size was due to the lack of obesity risk-allele score for many children. As a second sensitivity analysis, we tested the associations between weighted obesity risk-allele score, perinatal determinants, and age at AR in the fully adjusted model with multiple imputation. Finally, we conducted a sensitivity analysis excluding children born prematurely.

As complementary analyses, we first performed association analyses between perinatal factors, obesity risk-allele score, and risk of early AR, using logistic multiple regression models, estimating odds ratios and 95% confidence intervals (OR [95% CI]). We defined early AR as a rebound occurring before the first quintile of the distribution of age at AR. Then, we investigated each polymorphism used to establish the obesity risk-allele score and tested its individual association with age at AR. For this analysis, we assumed that the alleles of a given SNP have co-dominant effects, as is the case for the score computation. Therefore, we considered the SNPs (categorized as 0, 1, and 2) as continuous variables. The analysis of each SNP involved linear regression analysis of complete cases with $\beta$ corresponding to the change in days for one additional mutated allele carried.

We considered nominal $p < 0.05$ as statistically significant. We used R v3.5.1 for modeling BMI trajectories and estimating age at AR and SAS v9.4 (SAS, Cary, NC, on an AIX 7.1 platform) for all other analyses.

<table><caption>**Table 1** Characteristics of the included and excluded EDEN cohort in the study of adiposity rebound.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Mean (SD) or % (N)</th><th rowspan="2">p > |t| or χ<sup>2</sup></th></tr><tr><th>Included<br/>N = 1415</th><th>Excluded<br/>N = 492</th></tr></thead><tbody><tr><td>Sex (boys)</td><td>52.7 (746)</td><td>52.1 (254)</td><td>0.79</td></tr><tr><td>Preterm birth</td><td>5.5 (78)</td><td>6.5 (32)</td><td>0.41</td></tr><tr><td>Birth weight (kg)</td><td>3.29 (0.5)</td><td>3.23 (0.5)</td><td>0.03</td></tr><tr><td>Birth weight z-score</td><td></td><td></td><td>0.04</td></tr><tr><td>SGA</td><td>13.3 (184)</td><td>13.2 (62)</td><td></td></tr><tr><td>AGA</td><td>77.4 (1070)</td><td>81.3 (383)</td><td></td></tr><tr><td>LGA</td><td>9.3 (128)</td><td>5.5 (26)</td><td></td></tr><tr><td>Gestational weight gain (kg)</td><td>13.3 (4.7)</td><td>13.8 (5.2)</td><td>0.03</td></tr><tr><td>Maternal BMI (kg/m<sup>2</sup>)</td><td>23.1 (4.4)</td><td>23.6 (5.1)</td><td>0.08</td></tr><tr><td>Paternal BMI (kg/m<sup>2</sup>)</td><td>25.1 (3.6)</td><td>24.7 (3.7)</td><td>0.03</td></tr><tr><td>Never smoking during pregnancy</td><td>78.6 (1088)</td><td>59.4 (275)</td><td>4.0 × 10<sup>-16</sup></td></tr><tr><td>Maternal age at delivery (years)</td><td>29.9 (4.7)</td><td>28.2 (5.2)</td><td>3.0 × 10<sup>-10</sup></td></tr><tr><td>Maternal educational level (years)</td><td>14.0 (2.6)</td><td>12.5 (2.5)</td><td>4.5 × 10<sup>-30</sup></td></tr><tr><td>Paternal educational level (years)</td><td>13.3 (2.6)</td><td>12.2 (2.4)</td><td>3.7 × 10<sup>-18</sup></td></tr><tr><td>Obesity risk-allele score</td><td>21.6 (3.2)</td><td>21.3 (3.4)</td><td>0.15</td></tr><tr><td>Age at adiposity rebound (years)</td><td></td><td></td><td>0.08<sup>a</sup></td></tr><tr><td>Boys</td><td>5.5 (1.3)</td><td>—</td><td>—</td></tr><tr><td>Girls</td><td>5.4 (1.4)</td><td>—</td><td>—</td></tr></tbody></table>

BMI body mass index.

<sup>a</sup> *p* value of mean age at adiposity rebound comparing boys and girls.

## Results

### Population characteristics

Compared to mothers of excluded children, those of included children (n = 1415) had on average a higher number of years of study and were older at delivery (p = 3.0 × 10<sup>-10</sup>) (Table 1). Included infants were born at higher gestational age (p = 0.03) and had higher birth weight (p = 0.03) than excluded children but were similar in sex (p = 0.79) and maternal pregestational BMI (p = 0.08).

Among the 1415 children included, the mean (±SD) age of AR was estimated at 5.5 (±1.3) and 5.4 (±1.4) years for boys and girls, respectively (p = 0.08). The first quintile of the distribution of age at AR was estimated at 4.4 years.

The description of predictors according to the quintiles of age at AR is provided in Supplementary Table S2. Except for prematurity, gestational weight gain, maternal age at delivery, and obesity risk-allele score, all considered factors were significantly differed by quintile of age at AR.

### Determinants associated with the timing of AR

Multivariable associations with age at AR are displayed in Table 2. Model A showed that high maternal and paternal educational level were associated with later AR (β [±SE] = 15.7 [±6.2] and 14.4 [±6.0] days per year of study, p = 0.01 and 0.02, respectively), with no association with the other
<table><caption>**Table 2** Multivariable linear regression of factors associated with age at adiposity rebound (N = 1415<sup>a</sup>).</caption><thead><tr><th rowspan="2">Characteristics</th><th colspan="2">Model A</th><th colspan="2">Model B</th><th colspan="2">Model C</th><th colspan="2">Model D</th></tr><tr><th>β (SE)</th><th>p value</th><th>β (SE)</th><th>p value</th><th>β (SE)</th><th>p value</th><th>β (SE)</th><th>p value</th></tr></thead><tbody><tr><td colspan="9"><strong>Level 1<sup>b</sup></strong></td></tr><tr><td>Center (ref = Poitiers)</td><td>-0.8 (26.6)</td><td>0.98</td><td>-14.1 (26.2)</td><td>0.59</td><td>-22.3 (26.2)</td><td>0.40</td><td>-14.9 (26.3)</td><td>0.57</td></tr><tr><td>Maternal age at delivery (years)</td><td>2.8 (2.9)</td><td>0.33</td><td>5.1 (2.8)</td><td>0.070</td><td>3.7 (2.9)</td><td>0.19</td><td>3.9 (2.8)</td><td>0.17</td></tr><tr><td>Maternal educational level (years)</td><td>15.7 (6.2)</td><td>0.011</td><td>13.4 (6.1)</td><td>0.027</td><td>11.3 (6.1)</td><td>0.064</td><td>10.0 (6.1)</td><td>0.10</td></tr><tr><td>Paternal educational level (years)</td><td>14.4 (6.0)</td><td>0.017</td><td>8.5 (6.0)</td><td>0.16</td><td>7.4 (6.0)</td><td>0.22</td><td>7.0 (6.0)</td><td>0.24</td></tr><tr><td colspan="9"><strong>Level 2</strong></td></tr><tr><td>Obesity risk-allele score</td><td></td><td></td><td>-9.2 (4.6)</td><td>0.045</td><td>-9.1 (4.6)</td><td>0.045</td><td>-9.0 (4.5)</td><td>0.047</td></tr><tr><td>Maternal BMI (kg/m<sup>2</sup>)</td><td></td><td></td><td>-15.9 (3.1)</td><td>1.9 × 10<sup>-7</sup></td><td>-18.2 (3.1)</td><td>6.3 × 10<sup>-9</sup></td><td>-18.0 (3.1)</td><td>8.9 × 10<sup>-9</sup></td></tr><tr><td>Paternal BMI (kg/m<sup>2</sup>)</td><td></td><td></td><td>-15.6 (3.9)</td><td>5.6 × 10<sup>-5</sup></td><td>-15.3 (3.9)</td><td>6.9 × 10<sup>-5</sup></td><td>-15.3 (3.8)</td><td>6.7 × 10<sup>-5</sup></td></tr><tr><td colspan="9"><strong>Level 3</strong></td></tr><tr><td>Gestational weight gain (kg)</td><td></td><td></td><td></td><td></td><td>-7.5 (2.9)</td><td>0.009</td><td>-8.6 (3.0)</td><td>0.004</td></tr><tr><td>Smoking during pregnancy (ref = no)</td><td></td><td></td><td></td><td></td><td>-72.1 (32.6)</td><td>0.027</td><td>-62.3 (32.7)</td><td>0.057</td></tr><tr><td colspan="9"><strong>Level 4</strong></td></tr><tr><td>Preterm birth (yes vs no)</td><td></td><td></td><td></td><td></td><td></td><td></td><td>39.4 (57.7)</td><td>0.50</td></tr><tr><td>Birth weight z-score</td><td></td><td></td><td></td><td></td><td></td><td></td><td>24.1 (12.4)</td><td>0.053</td></tr><tr><td>Birth weight z-score [2]</td><td></td><td></td><td></td><td></td><td></td><td></td><td>-14.4 (7.3)</td><td>0.047</td></tr><tr><td>Sex (girls vs boys)</td><td></td><td></td><td></td><td></td><td></td><td></td><td>-47.2 (25.8)</td><td>0.067</td></tr></tbody></table>

socio-demographic factors considered in this model. Furthermore, all variables related to parental obesity (Model B) and prenatal and intrauterine factors (Model C) were significantly associated with age at AR: age at AR decreased with mean obesity risk-allele score (−9.2 [±4.6] days per allele, *p* = 0.05), maternal and paternal BMI (−15.8 [±3.1] days per kg/m², *p* = 1.9 × 10⁻⁷, and −15.6 [±3.9] days per kg/m², *p* = 5.6 × 10⁻⁵), gestational weight gain (−7.5 [±2.9] days per kg, *p* = 0.001) and smoking during pregnancy (−72.1 [±32.6] days, *p* = 0.03). Finally, birth weight z-score was the only child characteristic at birth associated with age at AR with a significant quadratic effect (*p* = 0.05).

To illustrate this last result, the mean age of AR and their confidence interval by category of birth weight z-score are presented in Fig. 2. On average, children born SGA had an age at AR 88 [±39] days lower than children born AGA, and 91 [±56] days earlier than children born LGA.

## Sensitivity and complementary analyses

In the sensitivity analysis, results were on the whole consistent with analyses limited to children with no missing data (complete-case analysis) (Supplementary Table S3). However, the associations observed between maternal educational level (Model A), smoking during

pregnancy (Model C) and age at AR were no longer significant. We also observed that the strength of association was higher with obesity risk-allele score and birth weight z-score, but lower with parental BMI. Results with the

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Adjusted means [SE] of age at adiposity rebound according to birth weight z-score categories (N = 1415<sup>a</sup>, model D<sup>b</sup>). <sup>a</sup>Missing data for covariates were imputed with the fully conditional specification method. <sup>b</sup>Adjusted for center, maternal age at delivery, maternal and paternal educational level, maternal and paternal BMI, gestational weight gain, smoking during pregnancy, prematurity, birth weight z-score and sex.</figcaption></figure>

<footnote><sup>a</sup>Missing data for covariates were imputed with the fully conditional specification method.</footnote>

<footnote><sup>b</sup>Relevance of a variable during the hierarchical regression analyses was determined with the corresponding model in which the variable of interest is first entered, regardless of its performance in the subsequent model(s).</footnote>
weighted obesity risk-allele score were consistent with, albeit slightly weaker, than those with the unweighted obesity risk-allele score (Supplementary Table S4). When preterm children were removed from the analyses, the association between birth weight z-score and age at AR became stronger and the risk of early AR for SGA was greater (2.23 [1.52–3.28]) (Supplementary Table S5).

The results of logistic regression analyses testing the association between obesity risk-allele score, perinatal factors, and odds of early AR, are summarized in Supplementary Fig. 4. Higher maternal educational level was associated with reduced odds of early AR. Odds of early AR were increased with each additional obesity risk-allele, higher maternal and paternal BMI, and higher gestational weight gain. Risk of early AR was higher for children born SGA than AGA children. We found no association between LGA children and risk of early AR. Risk of early AR was higher for girls than boys.

In analyses considering each polymorphism separately (Supplementary Fig. S5), we found no significant association. The associations between rs10913469 at the locus harboring SEC16B, rs713586 at the locus harboring RBJ/POMC and age at AR were at the limit of significance.

## Discussion

We found that genetic susceptibility to obesity was negatively associated with age at AR in the EDEN mother–child study. Additionally, we confirmed that some previously described perinatal factors, such as parental BMI, were associated with age at AR [2, 20] but also further identified novel predictors such as birth weight and mother’s gestational weight gain. Probability of early AR was increased with high gestational weight gain and low birth weight.

This is the first study to simultaneously analyze the association between perinatal factors, genetic predisposition to obesity and age at AR, taking the assumed hierarchical structure between the potential predictors into account. At the most distal level, maternal and paternal education were negatively associated with risk of early AR, in accordance with previous studies [2]. This association was attenuated with the inclusion of the more proximal groups of variables. A longitudinal study of 889 representative children in the United Kingdom showed no association between parental educational level or socioeconomic factors and time of AR [31]. This result was obtained after adjustment for parental BMI, which make our two studies concordant. Future studies could use proper mediation analyses to quantify the mediating contribution of each of these factors and help prioritize effective interventions to reduce social inequalities in obesity.

The genetic susceptibility to obesity was associated with an increased probability of early AR. The role of genetic factors in AR timing was investigated recently in a genome-wide association study including up to 7215 children in five European cohorts. In this study, Couto Alves et al. identified an association between two SNPs and age at AR. Although the rs1421085 (at the locus harboring *FTO*) was not associated with age at AR, we noted a statistically significant association for rs10913469 (at the locus harboring SEC16B) and an association at the limit of significance for rs713586 (at the locus harboring RBJ/POMC) but with reduced statistical power. In accordance with our results, the authors reported a negative association between an obesity risk-allele score, based on 97 SNPs, and child age at AR, but did not model the effects of other life factors [24]. Altogether, these results illustrate that the observed association between a score of genetic susceptibility to obesity and timing of AR is the result of a cumulative effect of individual SNPs. The authors also highlighted the relevance of using cumulative scores to address the genetic predisposition to early AR.

Consistent with previous research, we found that high parental BMI predicted early AR in children. A longitudinal study of 248 families enrolled in North Carolina concluded high maternal BMI as the robust predictor of early AR [2]. Another longitudinal study of a New Zealand cohort found that mothers and fathers of children with early BMI rebound had significantly heavier as compared with late rebounders [32]. Two other investigations also demonstrated that children with early AR had at least one parent with obesity [31, 33]. All these analyses recognized maternal BMI as a proxy variable reflecting the child’s susceptibility to obesity. However, we found the association between parental BMI and age at AR was independent of the obesity risk-allele score, which suggests that parental BMI might have an effect independent of genetic predisposition to obesity. Indeed, parental BMI also reflects an obesogenic environment (familial environment, food choices, eating patterns, and physical activity) that may affect the kinetics of adiposity during childhood [34, 35]. Moreover, we observed that the association of mother’s prepregnancy BMI seemed to be slightly higher, even if not significantly (results not shown), than that of paternal BMI, which also suggests a role for nutritional environment during prenatal life.

Our results regarding gestational weight gain support this latter hypothesis. To our knowledge, only one study has considered gestational weight gain as a potential predictor of age at AR, but no significant association was shown [36]. However, previous reports focused on the link between gestational weight gain and childhood overweight at various ages. In the current EDEN cohort, a positive association was found between gestational weight gain and child’s BMI at age 5–6 years [37]. Likewise, a meta-analysis of 37
cohort studies from Europe, North America, and Australia showed that high gestational weight gain was related with increased risk of overweight/obesity in childhood, with stronger associations effects observed at later ages [15]. Our results suggest that the age of AR could in part mediate the association between gestational weight gain and overweight/obesity in late childhood and adolescence described in this meta-analysis, but further research would be necessary to confirm this assumption.

We found a negative association between maternal smoking status during pregnancy and the timing of AR but not risk of early AR. A longitudinal analysis, based on mothers from the longitudinal Growth and Obesity Chilean Cohort Study, also showed no link between smoking status during pregnancy and risk of early AR [38]. In previous analyses of the EDEN cohort, Carles et al. applied a joint Bayesian model to explore the association between maternal smoking status during pregnancy and growth [39]: children from mothers who smoked throughout pregnancy were born smaller but had higher BMI in their first months of life until age 5 as compared with children born to non-smokers. In the Generation R study, maternal smoking during pregnancy was associated with increased BMI at age 4 years for children with normal birth weight and SGA at birth [40]. Conversely, in the ALSPAC cohort, smoking during pregnancy was not associated to child weight status at age 5 years [41]. Thus, the association between maternal smoking during pregnancy and rapid postnatal growth might be due to its impact on fetal development but also to confounding by socioeconomic factors. Further studies are needed to confirm or not any effect of smoking during pregnancy on age at AR.

For the first time, we identified birth weight for gestational age as a predictor of age at AR. Some previous studies did include birth weight in their multiple models but did not find any significant association with age at AR. For the “Children in Focus” group from the ALSPAC cohort, no difference in birth weight was observed between the groups of very early, early, and later AR [42]. In contrast, our results indicate that children born SGA were at increased risk of early AR, whereas children born LGA were not. A cross-sectional study in Lausanne showed an association between low birth weight and anthropometric obesity markers in adulthood [43]. Similarly, the review article of Hong and Chung highlighted the association between SGA children and greater adiposity and obesity in later life [17]. However, Kapral et al. found this association only in preterm children [18]. These findings highlight that SGA newborns tend to rapidly “catch up”, and that this “catch-up” is characterized by an increase in fat mass greater than lean mass [44]. Furthermore, several systematic reviews have shown a strong association between rapid infancy weight gain and greater risk of obesity, overweight or increased adiposity in later life

[45, 46]. The catch-up phenomenon would be due to different mechanisms of adaptation of a newborn exposed to undernutrition in utero, which could have remote impact on age at AR [47]. The respective role of intrauterine growth restriction vs postnatal growth acceleration in this association still needs to be disentangled.

Even if some of our results need confirmation and replication, they highlight a number of prenatal factors, including smoking during pregnancy, maternal prepregnancy BMI, gestational weight gain and birth weight for gestational age, possibly involved in the timing of AR. These findings strongly suggest that some intrauterine programming may also be at work in the etiology of this specific adiposity phenotype. Among the possible mechanisms, altered fetal nutrition could lead to persistent changes in the structure and function of adipose tissue, regulation of appetite and energy metabolism, which could affect the kinetics of adiposity and its rebound [15]. Epigenetic mechanisms could also be involved, but further studies are needed to better understand these [48].

The main strength of this study is that anthropometric data were assessed repetitively and prospectively from birth to age 13, which allowed us to estimate age at AR for a large number of children. As is often the case in cohort studies, the presence of selection bias at inclusion and missing data of covariates but also attrition bias limits the generalization of the results. However, multiple imputations and growth modeling of individual growth curves helped reduce the risk of these biases. Another limitation is that we could not consider ethnicity in our analyses, because the information was not collected in the EDEN cohort. Nonetheless, for 85% of children in the EDEN cohort, both parents were born in France and came from French regions where most of the population has a Caucasian origin. Moreover, the construction of the genetic score relies on the assumption of a co-dominant effect of each allele of a given SNP, which is not always true and could induce measurement errors. Our results, relying on observational data, need to be replicated in other studies with larger samples and with distinct designs to provide additional arguments for a causal interpretation of the observed associations. Further studies are needed to examine postnatal factors involved in the kinetics of BMI growth curve, to better characterize children at risk of early rebound. Indeed, the detection of early age at AR might be a useful tool for tailored and personalized child's growth monitoring aiming at preventing overweight/obesity and chronic disease in later life.

## Conclusion

We confirmed that genetic susceptibility to adult obesity plays a predictive role in the timing of age at AR. Timing of
AR in early childhood may contribute to the path from genetic susceptibility to obesity onset. Independent of this genetic background, we identified several prenatal predictors, including two novel factors, gestational weight gain and birth weight. These results strongly support that intrauterine programming may also be at work in the etiology of this specific adiposity phenotype.

**Acknowledgements** The authors thank the EDEN mother-child cohort study group, whose members are I. Annesi-Maesano, J.Y. Bernard, M.A. Charles, P. Dargent-Molina, B. de Lauzon-Guillain, P. Ducimetière, M. de Agostini, B. Foliguet, A. Forhan, X. Fritel, A. Germa, V. Goua, R. Hankard, B. Heude, M. Kaminski, B. Larroquey, N. Lelong, J. Lepeule, G. Magnin, L. Marchand, C. Nabet, F. Pierre, R. Slama, M.J. Saurel-Cubizolles, M. Schweitzer, and O. Thiebaugeorges.

**Author contributions** AHC and BH designed the research, wrote the manuscript, and analyzed the data. BH, KKO, MAC, AF, and BLG were responsible for data collection and genotyping in EDEN. All authors reviewed drafts, provided critical feedback, approved the final manuscript, and were responsible for the final content of the paper. BH had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

**Funding** The EDEN study was supported by Foundation for medical research (FRM), National Agency for Research (ANR), National Institute for Research in Public health (IRESP: TGIR cohorte santé 2008 program), French Ministry of Health (DGS), French Ministry of Research, INSERM Bone and Joint Diseases National Research (PRO-A), and Human Nutrition National Research Programs, Paris-Sud University, Nestlé, French National Institute for Population Health Surveillance (InVS), French National Institute for Health Education (INPES), the European Union FP7 programs (FP7/2007–2013, HELIX, ESCAPE, ENRIECO, Medall projects), Diabetes National Research Program (through a collaboration with the French Association of Diabetic Patients [AFD]), French Agency for Environmental Health Safety (now ANSES), Mutuelle Générale de l'Education Nationale a complementary health insurance (MGEN), French national agency for food security, French-speaking association for the study of diabetes and metabolism (ALFEDIAM). KKO is supported by the Medical Research Council (Unit program: MC_UU_12015/2). This research benefited from the assistance of the funding partners of the IReSP within the framework of the 2016 General call for projects–Prevention topic (HEUDE-AAP16-PREV-24).

## Compliance with ethical standards

**Conflict of interest** The authors declare no competing interests.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

## References

1. Rolland-Cachera MF, Deheeger M, Bellisle F, Sempé M, Guilloud-Bataille M, Patois E. Adiposity rebound in children: a simple indicator for predicting obesity. Am J Clin Nutr. 1984;39:129–35.

2. Ip EH, Marshall SA, Saldana S, Skelton JA, Suerken CK, Arcury TA, et al. Determinants of adiposity rebound timing in children. J Pediatr. 2017;184:151–6.e2.

3. Koyama S, Ichikawa G, Kojima M, Shimura N, Sairenchi T, Arisaka O. Adiposity rebound and the development of metabolic syndrome. Pediatrics. 2014;133:e114–9.

4. Campbell MW-C, Williams J, Carlin JB, Wake M. Is the adiposity rebound a rebound in adiposity? Int J Pediatr Obes. 2011;6: e207–15.

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

6. Whitaker RC, Pepe MS, Wright JA, Seidel KD, Dietz WH. Early adiposity rebound and the risk of adult obesity. Pediatrics. 1998;101:E5.

7. Rolland-Cachera MF, Cole TJ. Does the age at adiposity rebound reflect a critical period? Pediatr Obes. 2019;14:e12467.

8. Péneau S, Giudici KV, Gusto G, Goxe D, Lantieri O, Hercberg S, et al. Growth trajectories of body mass index during childhood: associated factors and health outcome at adulthood. J Pediatr. 2017;186:64–71.e1.

9. Aris IM, Rifas-Shiman SL, Li L-J, Kleinman KP, Coull BA, Gold DR, et al. Patterns of body mass index milestones in early life and cardiometabolic risk in early adolescence. Int J Epidemiol. 2019;0:11.

10. González L, Corvalán C, Pereira A, Kain J, Garmendia ML, Uauy R. Early adiposity rebound is associated with metabolic risk in 7-year-old children. Int J Obes. 2014;38:1299–304.

11. Barker DJP. The origins of the developmental origins theory. J Intern Med. 2007;261:412–7.

12. Tomar AS, Tallapragada DSP, Nongmaithem SS, Shrestha S, Yajnik CS, Chandak GR. Intrauterine programming of diabetes and adiposity. Curr Obes Rep. 2015;4:418–28.

13. Kwon EJ, Kim YJ. What is fetal programming?: a lifetime health is under the control of in utero health. Obstet Gynecol Sci. 2017;60:506–19.

14. Goldstein RF, Abell SK, Ranasinha S, Misso M, Boyle JA, Black MH, et al. Association of gestational weight gain with maternal and infant outcomes: a systematic review and meta-analysis. JAMA. 2017;317:2207.

15. Voerman E, Santos S, Patro Golab B, Amiano P, Ballester F, Barros H, et al. Maternal body mass index, gestational weight gain, and the risk of overweight and obesity across childhood: an individual participant data meta-analysis. PLoS Med. 2019;16: e1002744.

16. Gaillard R, Santos S, Duijts L, Felix JF. Childhood health consequences of maternal obesity during pregnancy: a narrative review. Ann Nutr Metab. 2016;69:171–80.

17. Hong YH, Chung S. Small for gestational age and obesity related comorbidities. Ann Pediatr Endocrinol Metab. 2018;23:4–8.

18. Kapral N, Miller SE, Scharf RJ, Gurka MJ, DeBoer MD. Associations between birthweight and overweight and obesity in school-age children: birthweight and childhood obesity risk. Pediatr Obes. 2018;13:333–41.
19. Aris IM, Rifas-Shiman SL, Li L-J, Kleinman K, Coull BA, Gold DR, et al. Pre-, perinatal, and parental predictors of body mass index trajectory milestones. J Pediatr. 2018;201:69–77.e8.

20. Péneau S, González-Carrascosa R, Gusto G, Goxe D, Lantieri O, Fezeu L, et al. Age at adiposity rebound: determinants and association with nutritional status and the metabolic syndrome at adulthood. Int J Obes. 2016;40:1150–6.

21. Speliotes EK, Willer CJ, Berndt SI, Monda KL, Thorleifsson G, Jackson AU, et al. Association analyses of 249,796 individuals reveal eighteen new loci associated with body mass index. Nat Genet. 2010;42:937–48.

22. Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518:197–206.

23. Elks CE, Heude B, de Zegher F, Barton SJ, Clément K, Inskip HM, et al. Associations between genetic obesity susceptibility and early postnatal fat and lean mass: an individual participant meta-analysis. JAMA Pediatr. 2014;168:1122.

24. Couto Alves A, De Silva NMG, Karhunen V, Sovio U, Das S, Taal HR, et al. GWAS on longitudinal growth traits reveals different genetic factors influencing infant, child, and adult BMI. Sci Adv. 2019;5:eaaaw3095.

25. Heude B, Forhan A, Slama R, Douhaud L, Bedel S, Saurel-Cubizolles M-J, et al. Cohort profile: the EDEN mother-child cohort on the prenatal and early postnatal determinants of child health and development. Int J Epidemiol. 2016;45:353–63.

26. Regnault N, Botton J, Forhan A, Hankard R, Thiebaugeorges O, Hillier TA, et al. Determinants of early ponderal and statural growth in full-term infants in the EDEN mother-child cohort study. Am J Clin Nutr. 2010;92:594–602.

27. Ego A, Prunet C, Blondel B, Kaminski M, Goffinet F, Zeitlin J. Courbes de croissance in utero ajustées et non ajustées adaptées à la population française. II—Comparaison à des courbes existantes et apport de l'ajustement. J Gynécologie Obstétrique Biol Reprod. 2016;45:165–76.

28. Li S, Zhao JH, Luan J, Luben RN, Rodwell SA, Khaw K-T, et al. Cumulative effects and predictive value of common obesity-susceptibility variants identified by genome-wide association studies. Am J Clin Nutr. 2010;91:184–90.

29. de Lauzon-Guillain B, Clifton EA, Day FR, Clément K, Brage S, Forouhi NG, et al. Mediation and modification of genetic susceptibility to obesity by eating behaviors. Am J Clin Nutr. 2017;106:996–1004.

30. Sovio U, Mook-Kanamori DO, Warrington NM, Lawrence R, Briollais L, Palmer CNA, et al. Association between common variation at the FTO locus and changes in body mass index from infancy to late childhood: the complex nature of genetic association through growth and development. PLoS Genet. 2011;7. https://doi.org/10.1371/journal.pgen.1001307.

31. Dorosty AR, Emmett PM, Cowin SD, Reilly JJ. Factors associated with early adiposity rebound. ALSPAC Study Team. Pediatrics. 2000;105:1115–18.

32. Williams SM, Goulding A. Patterns of growth associated with the timing of adiposity rebound. Obesity. 2009;17:335–41.

33. Magee CA, Caputi P, Iverson DC. Identification of distinct body mass index trajectories in Australian children. Pediatr Obes. 2013;8:189–98.

34. Reilly JJ. Physical activity, sedentary behaviour and energy balance in the preschool child: opportunities for early obesity prevention: symposium on 'behavioural nutrition and energy balance in the young'. Proc Nutr Soc. 2008;67:317–25.

35. Vilchis-Gil J, Galván-Portillo M, Klünder-Klünder M, Cruz M, Flores-Huerta S. Food habits, physical activities and sedentary lifestyles of eutrophic and obese school children: a case-control study. BMC Public Health. 2015;15:124.

36. Linares J, Corvalán C, Galleguillos B, Kain J, González L, Uauy R, et al. The effects of pre-pregnancy BMI and maternal factors on the timing of adiposity rebound in offspring. Obesity. 2016;24:1313–19.

37. Jacota M, Forhan A, Saldanha-Gomes C, Charles MA, Heude B, for the EDEN Mother-Child Cohort Study Group. Maternal weight prior and during pregnancy and offspring's BMI and adiposity at 5-6 years in the EDEN mother-child cohort: Maternal weight history and children's adiposity. Pediatr Obes. 2017;12:320–9.

38. Linares J, Corvalán C, Galleguillos B, Kain J, González L, Uauy R, et al. The effects of pre-pregnancy BMI and maternal factors on the timing of adiposity rebound in offspring: maternal factors and adiposity rebound. Obesity. 2016;24:1313–9.

39. Carles S, Charles M-A, Heude B, Ahmed I, Botton J, the EDEN Mother-Child Study Group. Joint Bayesian weight and height postnatal growth model to study the effects of maternal smoking during pregnancy. Stat Med. 2017;36:3990–4006.

40. Durmuş B, Ay L, Duijts L, Moll HA, Hokken-Koelega ACS, Raat H, et al. Infant diet and subcutaneous fat mass in early childhood: the Generation R Study. Eur J Clin Nutr. 2012;66:253–60.

41. Ong KKL, Preece MA, Emmett PM, Ahmed ML, Dunger DB. Size at birth and early childhood growth in relation to maternal smoking, parity and infant breast-feeding: longitudinal birth cohort study and analysis. Pediatr Res. 2002;52:863–7.

42. Hughes AR, Sherriff A, Ness AR, Reilly JJ. Timing of adiposity rebound and adiposity in adolescence. Pediatrics. 2014;134: e1354–61.

43. Jornayvaz FR, Vollenweider P, Bochud M, Mooser V, Waeber G, Marques-Vidal P. Low birth weight leads to obesity, diabetes and increased leptin levels in adults: the CoLaus study. Cardiovasc Diabetol. 2016;15:73.

44. Wells JCK, Chomtho S, Fewtrell MS. Programming of body composition by early growth and nutrition. Proc Nutr Soc. 2007;66:423–34.

45. Ong K, Loos R. Rapid infancy weight gain and subsequent obesity: Systematic reviews and hopeful suggestions. Acta Paediatr. 2006;95:904–8.

46. Monteiro POA, Victora CG. Rapid growth in infancy and childhood and obesity in later life—a systematic review. Obes Rev. 2005;6:143–54.

47. Ong KKL. Association between postnatal catch-up growth and obesity in childhood: prospective cohort study. BMJ. 2000;320:967–71.

48. Poston L. Maternal obesity, gestational weight gain and diet as determinants of offspring long term health. Best Pract Res Clin Endocrinol Metab. 2012;26:627–39.