Rapid Weight Gain, Infant Feeding Practices, and Subsequent Body Mass Index Trajectories: The CALINA Study

We aimed to study growth patterns according to rapid weight gain (RWG) and infant feeding practices during the first 120 days and whether infant feeding practices mediated the association between RWG in the first semester of life and subsequent body mass index (BMI) z-score in children from age 1 to 6. (1) Methods: 862 children from the Growth and Feeding during Lactation and Early Childhood in Children of Aragon study (CALINA in Spanish) were examined. Repeated-measures ANOVA analyses were conducted to assess growth trajectories according to RWG and type of feeding practice. The product of coefficients mediation method was used to assess the potential contribution of infant feeding practices to the association between RWG and BMI z-score. Mediation models were conducted using IBM SPSS-PROCESS Statistics for Windows, Version 26.0. Armonk, NY: IBM Corp. (2) Results: BMI and weight z-score trajectories were significantly higher in the RWG group and the formula-fed group. No significant differences were found regarding height. Infant feeding practices did not mediate the association between RWG and BMI z-score but were associated with BMI at 6 years. (3) Conclusions: Infant feeding practices and RWG determine different growth trajectories of BMI and weight during childhood. Although infant feeding practices did not mediate the association between early RWG and BMI later in life, formula feeding is independently related to higher BMI growth patterns later in childhood.


Introduction
Obesity prevalence has increased by 100% in the last 40 years worldwide and it does not only affect adults, but also children; in fact, over 40 million children under 5 are already overweight and are therefore predisposed to later overweight or obesity during adulthood [1]. According to the World Health Organization (WHO) European Childhood Obesity Surveillance Initiative (COSI) study, overweight and obesity prevalence in European children from 6 to 9 years old oscillates among countries between 19.3-49.0% of boys and 18.4-42.5% of girls when following the 2007 WHO growth references, with Belgium, the Czech Republic, and Norway being the countries with the lowest rates of overweight and obesity and Portugal, Slovenia, and Italy being the countries with the highest rates [2].
Obesity, defined by the WHO as an abnormal or excessive fat accumulation presenting a risk to health [3], has a multifactorial origin, involving mainly genetic and environmental factors [4]. Early-life environmental factors during the perinatal period are recognized as major factors in shaping obesity risk, both in childhood and even later in adulthood [5,6]. In particular, parental exposures such as higher pre-pregnancy body mass index (BMI), prenatal tobacco exposure, excessive gestational weight gain, high birth weight, and infant rapid weight gain (RWG) are factors strongly associated with future childhood obesity [5,7].
From the early stages of life, nutrition plays an important role in growth and development. Breastfeeding is considered by pediatric organizations, such as the Italian [8] and the American [9] Associations of Pediatrics and the WHO [10], as the best method of nourishment for babies and infants during their first year of life. However, there is still controversy about whether breastfeeding can be considered as a protective factor against later obesity [11], as breastfeeding is probably confounded by the effect of socioeconomic status of families [12]. Furthermore, rapid or excessive weight gain during the first months of life seems to be an important risk factor for developing childhood overweight or obesity [5]. In this sense, type of feeding is of relevance, given that formula feeding has been identified as a risk factor for RWG during the first 6 months of life [13], and, on the other hand, breastfeeding has been associated with lower odds of rapid increase in weight when evaluated in different time periods during the first 24 months [14,15]. In addition to isolated differences, there have been differences in growth trajectories in terms of weight and BMI z-scores between breastfed and formula-fed infants, with breastfed infants having lower weight-for-age z-score (WAz) and BMI-for-age z-score (BAz) trajectories [16].
These two factors, exclusive breastfeeding or formula feeding and RWG, occur in parallel and therefore might be inter-related. An example of evidence of this is the fact that formula-fed infants, as opposed to breastfed infants, are more likely to experience RWG [17]. However, there is little information about the impact of the type of infant feeding on the association between RWG and further risk of obesity development during childhood.
In current literature, most studies have shown associations between feeding practices and RWG. Meanwhile, other studies have focused on the effects that RWG has on growth trajectories and later overweight risk. Therefore, in order to try and explain whether type of feeding mediates the established association between RWG and body weight and BMI, the aim of this study is to describe growth patterns according to RWG (whether they were classified as RWgainers or not from 0 to 6 months of life) and type of infant feeding during the first 120 days (exclusive breastfeeding, formula feeding, or mixed feeding) and to assess the mediation effect of the type of feeding during the first 4 months of life on the association between RWG and BMI z-score from 1 to 6 years of age.

Design and Study Population
The study Growth and Feeding during Lactation and Early Childhood in Children of Aragon (Crecimiento y ALimentación en la Infancia en Niños Aragoneses, CALINA in Spanish) is an ongoing birth cohort study whose sampling design is described elsewhere in detail [18]. Briefly, CALINA is a representative birth cohort study of children born between March 2009 and February 2010 in the Autonomous Region of Aragon in Spain. Participants were recruited from primary care centers, where trained pediatricians oversaw measurements. The main objective of the CALINA study was to assess growth patterns, body composition, and feeding aspects in this population and to examine prenatal, postnatal, and socioeconomic factors that could have an influence on them.
In total, 1630 families were contacted during the first 2 weeks of life of the newborns in their primary care centers and invited to participate in the CALINA study, of which 1602 accepted to participate (acceptance rate: 98%). Exclusion criteria included presenting any malformations, diseases, Nutrients 2020, 12, 3178 3 of 14 or physical disabilities. After families accepted to participate in the study, data regarding prenatal factors and birth characteristics were obtained from mothers and their newborn children clinical histories, and direct interviews with the families. Perinatal information from children was obtained after enrollment, and children were periodically re-examined in primary care centers at 2 weeks, monthly (at 1, 2, 4, 6, and 9 months), and yearly (at 1, 2, 3, 4, 5, and 6 years of age). In every assessment, besides clinical evaluation of health indicators, growth indicators (weight and length or height) were measured by pediatricians and nurses that had previously been trained for consistency between measurements. Furthermore, information regarding feeding practices was obtained in order to know what type of feeding children were receiving each month.
For this study, preterm infants (with 36 weeks of gestational age or less at birth) and those without complete data regarding infant feeding practice, weight or height at 6 years (72 months) were excluded from the analysis (n = 740). In the end, the present study included 862 children. Differences between socioeconomic characteristics of the included and the excluded sample were analyzed, and significant differences were observed in maternal age (32.3 ± 5.01 years in the included vs. 31.21 ± 5.3 years in the excluded sample, p = 0.000), origin (14.6% of immigrants in the included sample vs. 33.2% in the excluded sample, p = 0.000), and education (25.2% of low-educated mothers in the included sample vs. 32.3% in the excluded sample and 39.8% of high-educated mothers vs. 32.8% in the excluded sample, p = 0.002).
The study was performed following the ethical guidelines of the Declaration of Helsinki 1964 [19]. Parents or legal guardians provided signed consent for participation. Ethical approval was obtained in June 2008 (P108/0021) from Aragon's Committee of Ethics in Clinical Research (Comité de Ética de la Investigación de la Comunidad de Aragón, CEICA).

Anthropometry (Outcome)
The primary outcomes were z-scores for BMI, weight and length/height assessed at 1, 2, 3, 4, 5, and 6 years old. Body weight was measured to the nearest 10 g using a weight scale/bioelectric impedance instrument (Tanita ® BC-418, Corporation of America, Inc., IL, USA). Child length up to 2 years was measured using a recumbent board. For children older than 2 years, barefoot body height was measured to the nearest 0.1 cm standing up in a fasting state and wearing light clothes using a stadiometer (SECA 225, SECA, Hamburg, Germany). Both height and weight were obtained by trained staff. BMI-, weight-and height-for-age-and-sex z-scores (BAz, WAz, and HAz, respectively) were calculated with the Anthro Software v 3.2.2, Geneva, Switzerland and Anthro Plus Software v 1.0.4, Geneva, Switzerland, according to the WHO growth references [20].

Rapid Weight Gain (Exposure)
The primary exposure was RWG, a dichotomous variable defined as "no" or "yes". RWG is defined as a positive change in the WAz greater than 0.67 between two different ages in childhood [21]. In this study, we determined RWG in the first semester of life, from birth to 6 months of age. Determination of z-score values of weight for age at birth and at 6 months for girls and for boys was performed using the WHO Anthro Software ® , according to the WHO growth standards of 2006-2007 [20].

Type of Infant Feeding (Mediator)
Children were classified according to the type of infant feeding they received during the first 120 days of life into 3 categories: 1 = Formula, 2 = Mixed feeding, and 3 = Breast milk. Formula feeding included children that received only formula during the first 120 days of life, and the third category included children that had exclusively received breast milk during this time. The second category, mixed feeding, included children that received both formula and breast milk. This group included those that received 1 or 2 months of exclusive breastfeeding and then a combination of both or those that received breast milk and formula since they were born.

Covariates
Covariates were identified a priori as potential confounders based on known associations between RWG and further weight gain. These covariates included birth weight, gestational age, maternal and paternal education, maternal and paternal BMI before pregnancy, parental origin, and maternal smoking during pregnancy. (I) Birth weight in grams was a continuous variable. (II) Gestational age in weeks was a continuous variable. (III) Maternal and paternal education level: parents reported their highest level of studies achieved (none, basic, graduate, or postgraduate), which were re-categorized into 3 categories: low (none or basic), medium (graduate), or high (postgraduate). (IV) BMI from parents before pregnancy was obtained by a face-to-face interview in which they reported weight and height and was thereafter calculated as weight in kilograms divided by the square of height in meters (kg/m 2 ). (V) Parental ethnicity: families were classified into 2 categories according to the origin of their parents, which could be Spanish or immigrant, and (VI) Maternal smoking during pregnancy: this variable was a dichotomous variable that indicated whether the mother had smoked during pregnancy, regardless of the number of cigarettes.

Statistical Analyses
Statistical analyses were carried out using IBM SPSS Statistics (IBM SPSS Statistics for Windows, Version 26.0. Armonk, NY: IBM Corp.). A descriptive analysis was conducted using mean and standard deviation for continuous variables and frequencies and percentages for categorical variables according to having developed RWG in the first semester of life.
To test whether the trajectories of BAz, WAz, and HAz differed based on the presence of RWG and the type of infant feeding, we used repeated-measures ANOVA in two separate models [22]. The first model included RWG as the main effect and BAz, WAz, and HAz at birth, 6 months, 1, 2, 3, 4, 5, and 6 years as the repeated-measures variable and the second model included infant feeding practice as the main effect and birth, 6-month, and yearly BAz, WAz, and HAz as the repeated-measures variable. These models were also conducted adjusting by birth weight, gestational age at birth, maternal education, maternal and paternal BMI before pregnancy, parental origin, and maternal smoking during pregnancy.

Mediation Analyses
To assess whether the associations between RWG and BMI z-score were mediated by infant feeding, the following steps were performed.
To qualify as a mediator, the presumed mediator (M) has to be associated with the predictor variable (a-pathway) and also with the outcome variable (b-pathway) [23]. This was assessed by performing bivariate correlations. Then, based on the hypothesis that infant feeding practice could mediate the associations between RWG and BAz during early childhood, the mediation models were examined using the PROCESS macro 3.1.4 software for SPSS [24]. In PROCESS, Model 4 software was applied for simple mediations. The independent variable was RWG, the dependent variable was BAz at 6 years, and the potential mediator variable was infant feeding practice. The covariates included were birth weight, gestational age at birth, maternal education, maternal and paternal BMI at child's birth, parental origin, and maternal smoking during pregnancy. Indirect effects were computed using a bias-corrected bootstrapping procedure (5000); if the 95% confidence interval did not include 0, it meant that the mediation effect was significant.

Results
Characteristics of rapid weight gainers and non-rapid weight gainers included in the analysis are shown in Table 1. Mean of birth weight was significantly higher in the non-rapid weight gainers group (p < 0.001), but from 6 months onward, children in the RWG group presented significantly higher BAz at 6, 12, 24, 36, 60, and 72 months (p < 0.01 in all cases). Differences in BAz at 48 months were not significantly different across groups. Trajectories of BAz, WAz, and HAz in rapid weight gainers and non-rapid weight gainers and according to infant feeding practice are shown in Figures 1 and 2, respectively. Figure 1 shows the non-adjusted repeated-measures ANOVA models of z-scores of BMI, weight, and length-height. p-values from Box's tests are presented on the z-score trajectory of each figure. This value represents the significance of rejecting the assumption of equality variance-covariance which is significant and therefore indicates that there are differences across trajectories of groups. Models were also performed adjusting for covariates that we identified a priori as potential confounders based on known associations with feeding type and BAz.     Figure 1 shows the non-adjusted repeated-measures ANOVA models of z-scores of BMI, weight, and length-height. p-values from Box's tests are presented on the z-score trajectory of each figure. This value represents the significance of rejecting the assumption of equality variance-covariance which is significant and therefore indicates that there are differences across trajectories of groups. Models were also performed adjusting for covariates that we identified a priori as potential confounders based on known associations with feeding type and BAz.
Participants in the RWG group presented significantly higher BAz and WAz from birth to 72 months (p = 0.006 and p = 0.043, respectively). HAz trajectories did not differ across groups (p = 0.828). After adjusting for covariates (birth weight, gestational age at birth, maternal and paternal education, maternal and paternal BMI at child's birth, ethnicity, and maternal smoking during pregnancy), BAz trajectories still differed between groups (p = 0.001). WAz and HAz trajectories were not statistically different between groups (p = 0.16 and p = 0.89, respectively).
Trajectories of growth indicators according to infant feeding practice in the first 120 days of life ( Figure 2) showed that formula-fed infants had higher BAz and WAz compared to infants fed with mixed-feeding or breast milk (p < 0.001 and p = 0.005, respectively). HAz trajectories did not differ across categories (p = 0.560). After adjusting for birth weight, gestational age, maternal and paternal education, maternal and paternal BMI, parental origin, and maternal smoking, the trajectories of BAz and WAz remained statistically different among infant feeding practice groups (p < 0.01 in both cases). HAz trajectory was not statistically different between groups (p > 0.05) in the adjusted model.

Mediation Results
Given that the presumed mediator (M) has to be associated with the predictor variable (a-pathway) and also with the outcome variable (b-pathway) to qualify as a mediator; Pearson correlations were performed to confirm this by assessing if there are associations between infant feeding practices and RWG development in the first semester of life, infant feeding practices, and BMI z-score at 1, 2, 3, 4, 5, and 6 years ( Table 2). In this step, we obtained that RWG was significantly correlated with BAz at 1, 2, 3, 5, and 6 years and that infant feeding practices (1 = Breast milk, 2 = Mixed feeding, and 3 = Formula milk) were positively associated with RWG development, which indicated that the mediation analyses could only be performed considering BAz at 6 years as the dependent variable.  ALimentación en la Infancia en Niños Aragoneses in Spanish), b BMI z-score (BAz), body mass index z-score for age and sex according to World Health Organization (WHO) growth standards; RWG, rapid weight gain; BMI, body mass index; y, years. c RWG, Rapid weight gain was defined as a positive change of >0.67 in the BMI z-score a . d Maternal and paternal education were classified as none or basic, medium, or high . e Maternal and paternal BMI (body mass index) before pregnancy was calculated from self-reported weight and height with the formula kg/m 2 . f Infant feeding practice was a categorical variable classified as 1 = Breast Milk, 2 = Mixed feeding (formula milk and breast milk), and 3 = Formula milk, indicating the type of practice received during the first 120 days. Note: Bivariate Pearson correlations are presented. * Correlation is significant at the 0.05 level, ** Correlation is significant at the 0.01 level, *** Correlation is significant at the 0.001 level.
Results of mediation can be found in Table 3 and a graphical explanation of the model is presented in Figure 3. Non-adjusted mediation analyses showed that Pathway A' (association between RWG and infant feeding practice), indicating that formula intake was positively associated with RWG development. Pathway B' (association between type of feeding and BAz at 6 years) was also positive and significant. Pathway C, which represents the association between RWG and BAz was positive and significant (p < 0.02), as was Pathway C', which represents the previous associations but considering the mediator as a covariate. After adjustments, the significance in the pathways was the same as in the unadjusted model, except for Pathway A', which was not significant.
Nutrients 2020, 12, x FOR PEER REVIEW 10 of 15 Figure 3. Graphical illustration of the significant interactions between rapid weight gain, type of feeding, and BAz at 6 years.  Simple mediation analysis adjusted by birth weight, gestational age at birth, maternal education, maternal and paternal BMI at child's birth, ethnicity, and maternal smoking during pregnancy. Pathway A': Association between RWG and infant feeding practice. Pathway B': Association between infant feeding practice and BAz at 6 years. Pathway C': Direct association of RWG and BAz after adjustment for mediator (infant feeding practice). Pathway C: The total effect (c) shows the association between RWG in the 1st semester of life and BAz at 6 years. A*B: Indirect effect of the infant feeding practice on the association between RWG and BAz. Rapid weight gain was defined as a positive change of >0.67 in the z-score of BMI for age and sex according to World Health Organization (WHO) growth standards. Infant feeding practice was a categorical variable classified as 1 = Breast milk, 2 = Mixed feeding (formula milk and breast milk), and 3 = Formula milk. BMI z-score was calculated according to the WHO growth standards. Abbreviations: RWG, rapid weight gain; BAz, body mass index z-score.; CI: confidence interval. Statistically significant results are shown in bold font. Note. Unstandardized coefficients are presented, SE; standard error, CI; confidence interval, * p-value is significant at the 0.05 level and *** p-value is significant at the 0.001 level, -= Not significant indirect effects.
The indirect effect of the potential mediator, type of feeding, indicated that the mediating effect of infant feeding practice on the relationship between RWG and BMI z-score (AB coefficient) was not statistically significant in either the unadjusted or the adjusted model, which means that there is no significant difference in the association between RWG and BAz at 6 years, depending on the type of feeding children received in the first 120 days of life. In general, there is a statistically significant direct effect between RWG and BAz from 0 to 6 years old, regardless of covariates, but infant feeding practice during the first 120 days did not mediate this association.

Discussion
The present study aimed to describe growth patterns according to RWG (whether they were  Simple mediation analysis adjusted by birth weight, gestational age at birth, maternal education, maternal and paternal BMI at child's birth, ethnicity, and maternal smoking during pregnancy. Pathway A': Association between RWG and infant feeding practice. Pathway B': Association between infant feeding practice and BAz at 6 years. Pathway C': Direct association of RWG and BAz after adjustment for mediator (infant feeding practice). Pathway C: The total effect (c) shows the association between RWG in the 1st semester of life and BAz at 6 years. A*B: Indirect effect of the infant feeding practice on the association between RWG and BAz. Rapid weight gain was defined as a positive change of >0.67 in the z-score of BMI for age and sex according to World Health Organization (WHO) growth standards. Infant feeding practice was a categorical variable classified as 1 = Breast milk, 2 = Mixed feeding (formula milk and breast milk), and 3 = Formula milk. BMI z-score was calculated according to the WHO growth standards. Abbreviations: RWG, rapid weight gain; BAz, body mass index z-score.; CI: confidence interval. Statistically significant results are shown in bold font. Note. Unstandardized coefficients are presented, SE; standard error, CI; confidence interval, * p-value is significant at the 0.05 level and *** p-value is significant at the 0.001 level, -= Not significant indirect effects.
The indirect effect of the potential mediator, type of feeding, indicated that the mediating effect of infant feeding practice on the relationship between RWG and BMI z-score (AB coefficient) was not statistically significant in either the unadjusted or the adjusted model, which means that there is no significant difference in the association between RWG and BAz at 6 years, depending on the type of feeding children received in the first 120 days of life. In general, there is a statistically significant direct effect between RWG and BAz from 0 to 6 years old, regardless of covariates, but infant feeding practice during the first 120 days did not mediate this association.

Discussion
The present study aimed to describe growth patterns according to RWG (whether they were classified as RWgainers or not from 0 to 6 months of life) and type of infant feeding during the first 120 days (exclusive breastfeeding, formula feeding, or mixed feeding) and to study whether these feeding practices mediate the association between RWG in the first semester of life and later BAz in children from 1 to 6 years. We found that BAz and WAz trajectories were significantly higher in rapid weight gainers compared to non-rapid weight gainers. Regarding infant feeding practices, BAz and WAz trajectories were significantly higher in the formula-fed group than breastfed or mixed-fed children. However, no significant differences were found for height in either of the analyses. Infant feeding practices for the 120 first days of life did not mediate the association between RWG and BAz at 6 years, although it appeared to be independently associated with BAz at this age.
An important aspect of the present study is that z-scores have been used to assess growth in children, which is widely recognized as the best system for analysis and presentation of anthropometric data because of its advantages compared to other methods. Furthermore, given that they are sex-independent, their use permits the evaluation of children's growth status by combining sex and age groups [25] and has been widely used in recent studies with similar aims [26][27][28].
In the present investigation, and after adjusting for birth weight, gestational age at birth, parental education, parental BMI before pregnancy, ethnicity, and maternal smoking during pregnancy, the type of infant feeding (breastfeeding, mixed feeding, or formula feeding) during the first 120 days did not have a direct effect on RWG. However, an important direct effect of breastfeeding on long term BAz (BMI at 6 years old) was visible. Even though this study did not find the causes behind this association, results showed a small but significant protective effect of breastfeeding on the development of excess adiposity (BAz). Additionally, this effect appeared to be independent of the effect that breastfeeding had on RWG.
According to some studies, formula-fed infants are more likely to have RWG than breastfed infants. However, in a recent meta-analysis, only limited evidence (if any) was found in the association between formula feeding practices and RWG in infancy [29]. In our analysis, the association between type of infant feeding and RWG from 0 to 6 months of age showed a statistically significant effect only in the non-adjusted models. In fact, this association was no longer significant after adjusting for the abovementioned covariates. The increased risk of RWG in formula-fed infants may be due to a number of different factors, such as feeding to schedule versus feeding on demand [17]. It may also be due to the mode of milk delivery [30], since breastfeeding may be associated with better appetite regulation in the long term in comparison to bottle-feeding, regardless of the type of milk in the bottle [31]. It can be suggested that formula feeding is not directly associated with RWG, but some behaviors related to it, such as the use of high protein formula or adding cereals into the bottle, may be related [29].
Regarding the relationship between infant feeding practices and BAz, there was only a negative significant association between breastfeeding and BAz at 6 years old in the fully adjusted model. In other studies, breastfeeding was only significant in raw models, but when adjusting for important confounders, this association did not remain significant [32]. In our study, although this protective effect of breastfeeding could only be seen after several years, the effect was small but significant. Similarly, delaying the introduction of bottle feeding had a protective effect against obesity at 6 years of age in a sample of Spanish children, highlighting the need for greater support of breastfeeding to avoid future childhood obesity [11,33].
Similar to other studies [34][35][36], we found an important direct effect of RWG on children's BMI and, moreover, not mediated by breastfeeding during the first 120 days. This effect persisted over the years and remained significant at 6 years old. Even though the mechanism through which infant RWG programs subsequent adiposity remains unclear, it has been speculated that this association may be influenced by maternal factors and birth weight.
In this study, the average birth weight was significantly lower among rapid weight gainers compared to non-rapid weight gainers. Because of the potential influence of low birth weight on accelerated postnatal catch-up growth [37], which occurs typically in the first 24 months of postnatal life in infants with low birth weight, birth weight was included as a covariate in all models.
Nevertheless, from 6 months onwards, those rapid weight gainers presented a higher BAz at all ages (which remained significant after adjustments) and WAz (which did not remain significant after adjustments). Particularly, it has been shown that infants of low birth weight are more likely to have higher adrenal androgen levels, insulin resistance, and central fat deposition. They are therefore prone to weight gain [38], which may explain at least part of our results. In this line, it is important to note that our sample also differed in terms of birth weight, and even though it was included in the models as a covariate, it may be another starting point that determines different growth trajectories between groups.
We also adjusted for different maternal factors and birth weight and found a statistically significant result between RWG and BMI, suggesting that there must be other factors that are playing a role independently. Beyond birth weight, gestational age at birth, parental education, parental BMI before pregnancy, ethnicity, and maternal smoking during pregnancy, factors that may play an important role include dietary intake [39], eating behaviors [40], parental feeding practices [41], and physical activity levels [42].
Finally, even though there is a significant difference in growth trajectories, especially in the BAz, according to RWG in the first semester of life and infant feeding practices, no mediation effect of infant feeding practice was observed in this association.

Strengths and Limitations
To the best of our knowledge, this is the first study that has analyzed the mediation effect that infant feeding practices might have in the association with BAz during childhood in a large cohort of Spanish children followed from birth to 6 years old, taking into account perinatal factors from their parents. The prospective collection of data on a wide range of risk factors extending from pregnancy through infancy and the ability to adjust for several important confounding socioeconomic factors are among the strengths of this study.
However, there are some limitations to this research that should be acknowledged. Firstly, the included children are not representative of the Spanish population since the region in which the study took place (Aragon) covers a limited geographic area within the country (northeast) and results might not be extrapolated to the whole Spanish population. Secondly, an analysis of further confounders highly associated with obesity levels, such as dietary patterns, sedentary behaviors, physical activity, sleep duration, or family income, were not included in this study. Thirdly, there are some measures that were reported by parents that could imply a reliance problem (such as parental weight and height and their education). Furthermore, information regarding specific feeding practices, such as feeding to schedule or feeding on demand and the addition of cereals to bottles, was not included in the analyses, and it may be important to consider their effect on growth and weight gain. Moreover, a selection bias cannot be precluded as there were participants (mainly children whose parents were originally from Eastern European countries, Africa, and Latin America and had lower parental education) who did not complete all information required or did not continue the study at follow-up. Furthermore, feeding practice report may be prone to bias, given that self-administered questionnaires and health staff registrations on a monthly or less frequent basis have shown to be measurements that may overestimate the duration of exclusive breastfeeding [43]. Future studies should consider addressing other potential factors, such as dietary intake, physical activity, and sleep, which are lifestyle factors that may also affect growth trajectories in addition to type of feeding during the first months of life.

Conclusions
Both infant feeding practices and RWG determine different growth trajectories in terms of BMI and weight during childhood. Even though infant feeding practices do not seem to mediate the association between early RWG and BMI later in life, formula feeding is independently related to higher BMI growth patterns later in childhood. Assessment and Promotion and the European Regional Development Fund (ERDF). P. F-B received financial support by a grant from the Aragón's Regional Government (Diputación General de Aragón, DGA). Publication expenses for this article have been supported by Cátedra ORDESA-University of Granada, Spain, as part of the Special Issue Early Nutrition and Re-programming of Health and Disease.