Skip to main content

The temporal relationship between parental concern of overeating and childhood obesity considering genetic susceptibility: longitudinal results from the IDEFICS/I.Family study

Abstract

Background

Many genes and molecular pathways are associated with obesity, but the mechanisms from genes to obesity are less well known. Eating behaviors represent a plausible pathway, but because the relationships of eating behaviors and obesity may be bi-directional, it remains challenging to resolve the underlying pathways. A longitudinal approach is needed to assess the contribution of genetic risk during the development of obesity in childhood. In this study we aim to examine the relationships between the polygenic risk score for body mass index (PRS-BMI), parental concern of overeating and obesity indices during childhood.

Methods

The IDEFICS/I.Family study is a school-based multicenter pan-European cohort of children observed for 6 years (mean ± SD follow-up 5.8 ± 0.4). Children examined in 2007/2008 (wave 1) (mean ± SD age: 4.4 ± 1.1, range: 2–9 years), in 2009/2010 (wave 2) and in 2013/2014 (wave 3) were included. A total of 5112 children (49% girls) participated at waves 1, 2 and 3. For 2656 children with genome-wide data we constructed a PRS based on 2.1 million single nucleotide polymorphisms. Z-score BMI and z-score waist circumference (WC) were assessed and eating behaviors and relevant confounders were reported by parents via questionnaires. Parental concern of overeating was derived from principal component analyses from an eating behavior questionnaire.

Results

In cross-lagged models, the prospective associations between z-score obesity indices and parental concern of overeating were bi-directional. In mediation models, the association between the PRS-BMI and parental concern of overeating at wave 3 was mediated by baseline z-BMI (β = 0.16, 95% CI: 0.10, 0.21) and baseline z-WC (β = 0.17, 95% CI: 0.11, 0.23). To a lesser extent, baseline parental concern of overeating also mediated the association between the PRS-BMI and z-BMI at wave 3 (β = 0.10, 95% CI: 0.07, 0.13) and z-WC at wave 3 (β = 0.09, 95% CI: 0.07, 0.12).

Conclusions

The findings suggest that the prospective associations between obesity indices and parental concern of overeating are likely bi-directional, but obesity indices have a stronger association with future parental concern of overeating than vice versa. The findings suggest parental concern of overeating as a possible mediator in the genetic susceptibility to obesity and further highlight that other pathways are also involved. A better understanding of the genetic pathways that lead to childhood obesity can help to prevent weight gain.

Trial registration

Registry number: ISRCTN62310987 Retrospectively registered 17 September 2018.

Introduction

Obesity is a heritable and highly polygenic chronic disease [1,2,3,4]. Genetic factors explain between 40 to 60% of the proportion of the variability in body mass index (BMI) during childhood [5]. With the advance of genomic research, there has been a shift from twin and family studies towards using measured genes in genome-wide association studies. To date, the largest meta-analysis of genome-wide association studies of BMI among adults has led to the discovery of 751 single nucleotide polymorphisms (SNPs) associated with BMI; these SNPs collectively explain 6% of the phenotypic variance in BMI [4]. Khera et al. (2019) estimated that 23% of the variation in BMI is accounted for by common variants. When rare variants are also assessed by whole genome sequencing, the proportion of variance in BMI explained by measured genetic variants rises to 40% [6], which is close to the heritability estimates from twin and family studies. Using information from 2.1 million measured and imputed common variants irrespective of genome-wide significance [7], the newest polygenic risk score (PRS) is a single measure that quantifies the inherited susceptibility to a disease [8].

Both genes and environment can influence behaviors and physiology and they are involved in the regulation of energy intake and energy expenditure [5]. A low basal metabolic rate is partly heritable (about 40%) [9, 10] and is a risk factor for weight gain [11]. Appetite traits are also partly heritable (50–85%) in infancy and childhood [12,13,14] and are related to weight gain [15]. However, more research is needed to identify what pathway is stronger, as it is not yet established whether obesity susceptibility genes are more likely to influence body weight through metabolism or appetite traits. Since obesity susceptibility genes are highly expressed in the central nervous system [3, 16, 17], most previous research has focused on the role of appetite traits as behavioral pathways to obesity [5]. Results from a British cross-sectional study of 10-year old twins [18] and from a prospective French study of 1 to 5-year old children [19] showed that appetite traits mediate the genetic susceptibility to childhood obesity.

These earlier studies have used PRS-BMI based on a small number of variants (32 or fewer) and thus had a lower predictive power for capturing the obesity genetic risk [7]. Three previous PRS-BMI mediation studies have included all available common genetic variants [20,21,22], but they did not examine the temporal direction of appetite traits and obesity. Longitudinal studies using genetically sensitive designs to investigate the direction of pathways between eating behaviors and weight gain are needed to advance our understanding in the field [23].

Thus, we examined the temporal relationship between parental concern of overeating and obesity indices to examine the pathways by which the genetic susceptibility to obesity expresses itself by using children’s longitudinal data from the IDEFICS/I.Family European multicenter cohort and the latest available PRS for BMI.

Methods

Participants

The IDEFICS/I.Family study is a pan-European multicenter cohort that aimed at investigating eating habits and lifestyle factors in children and adolescents from eight countries (Belgium, Cyprus, Estonia, Germany, Hungary, Italy, Spain, and Sweden) [24]. The baseline examination (wave 1) included 16,229 children aged 2 to 9.9 years old. Follow-up examinations took place 2 years (wave 2, n = 13,586) and 6 years later (wave 3, n = 9639) (Fig. 1). A total of 2656 children were genotyped. Children from Cyprus were not selected for genotyping to minimize population stratification. In waves 1 (baseline) and 2, parents (or legal guardians) completed a self-administered questionnaire about their children’s health and lifestyle. In wave 3, questionnaires were completed by parents or their legal guardians for children up to 11 years old and for themselves and also by children for themselves if aged ≥12 years. Children participating in the study were asked to donate venous or capillary blood and to provide a saliva sample for DNA [25].

Fig. 1
figure 1

Flowchart diagram of participants

Research ethics committees in each country approved the study in accordance with the ethical standards of the Declaration of Helsinki. All children and parents gave informed consent to participate in the study.

Obesity indices

Height was measured with a portable stadiometer (Seca GmbH & Co. KG., Hamburg, Germany) to the nearest 0.1 cm and weight with a TANITA digital scale (TANITA Europe GmbH, Sindelfingen, Germany) to the nearest 0.1 kg. Both measurements were performed in the morning in fasting conditions and light clothing [26]. BMI was calculated from height and weight (kg/m2). Waist circumference was measured with an inelastic tape (Seca 200, Birmingham, UK) in the upright position with a relaxed abdomen at the midway point between the lowest rib and the iliac crest to the nearest 0.1 cm. Age- and sex- specific BMI z-scores (z-BMI) were calculated according to Cole [27] and age- and sex- specific WC z-scores (z-WC) were calculated according to the IDEFICS reference values [28].

Family food environment

Information on the family food environment was collected by questionnaire at all study waves. Parents or legal guardians answered questions about eating- and weight-related behaviors as well as parental worries toward their child’s food style for children aged under 12 years. Children aged 12 or older were not administered this questionnaire.

Questions were selected from the questionnaire from Baughcum et al. [29], with 11-items at waves 1 and 2 and 10-items at wave 3. The questionnaire covered five constructs (difficulty on child feeding; concern about children overeating and becoming overweight; pushing the child to eat more; the situation and structure during feeding; age-inappropriate feeding), with five possible response categories from “never” to “always”. Parents answered two more items: “How concerned are you about your child eating too much when you are not around him/her?” with four possible responses from “unconcerned” to “concerned”, and “How often does your child eat doing something else (e.g. watching TV, playing, sitting at a computer, looking at a book)?” with four possible response categories from “never” to “on several occasions per day”.

Derivation of eating behavior patterns

We derived eating behavior patterns by principal component analyses using 13-items (waves 1 and 2) and 12-items (wave 3) from the family food environment questionnaire described above and shown in Table S1 (Additional files). Based on eigenvalues > 0.80 we retained 4 components which explained 67% (wave 1) and 68% (waves 2 and 3) of the variance of the family food environment questions. Factor loadings ≥0.30 were considered to contribute to the eating behavior pattern and were used to label the 4 components. Out of them, parental concern of overeating was the component with the largest variance (about 30% in all three study waves) and included the following items: had to stop child from eating too much”; “thought about putting child on a diet to keep him/her from becoming overweight”; “worried about child eating too much” and “child eating too much when you are not around” (see Additional files, Table S1 for more details of the different eating behavior patterns).

Genotyping and Quality Control & Polygenic scoring

DNA was extracted from either saliva or blood samples using standard methods. The samples of 3515 children were genotyped using the UK Biobank Axiom array (Thermo-Fisher Scientific, Santa Clara, USA) in two batches (2015 and 2017). Sample and genotype quality control measures were applied following the recommendations of Weale [30], genome-wide imputation was done with Minimac3 v5 (https://genome.sph.umich.edu/wiki/Minimac3) resulting in 3098 children and 3,424,677 genotypes with an estimated posterior genotype probability > 0.8 and a minor allele frequency ≥ 0.05. EMMAX (https://genome.sph.umich.edu/wiki/EMMAX) was applied to calculate a genetic relatedness matrix in order to estimate the degree of relatedness within the study sample. A PRS-BMI was calculated as proposed and validated by Khera et al. [7] using the same reference population from Locke et al. [3]. It consists of 2,100,302 SNPs and is based on genome-wide summary statistics for BMI from European ancestry populations (~ 300,000 samples) [3].

Covariates

Six covariates (parental education, parental income, well-being, screen time duration, playing outside and fruit and vegetable consumption) were assessed by questionnaires. Glycated hemoglobin was analyzed by high-performance liquid chromatography.

Parents were asked to specify their education and income. Parental education was reported by one of the parents according to the International Standard Classification of Education (ISCED) [31], the maximum ISCED level of both parents was calculated. ISCED categories were as follows: low (levels 1–2); medium (levels 3–4) and high (levels 5–6). Parental income was asked according to the monthly net income of the household after taxes and deductions and categorized as follows: low, low-medium, medium, medium-high, and high. Parents filled out the KINDL® questionnaire on their child’s well-being and a score that ranges from 0 to 48 points was calculated [32].

Screen time in hours per week was calculated by summing up the total screen time spent on audiovisual media (TV, video, DVD, computer, game console) on weekdays and weekend days (5 x weekday + 2 x weekend).

Playing outdoors as a proxy for physical activity was calculated in hours per day and derived from hours playing outdoors on weekdays and hours playing outdoors on weekend days (5 x weekday + 2 x weekend)/7).

Information on fruits and vegetables consumption was calculated as the sum of the reported intake of raw and cooked vegetables and fresh fruit, which was obtained by a validated food frequency questionnaire [33, 34].

Glycated hemoglobin (HbA1c) was analyzed by high-performance liquid chromatography (AUTO-GA variant). Details on laboratory analyses can be found in Peplies et al. [35].

Statistical analyses

Descriptive characteristics were presented as means and standard deviations (SDs) for continuous variables and as numbers and percentages for categorical variables. We used chained random forest imputation to replace missing data on covariates. Missing data were randomly imputed to increase the number of observations and to reduce the probability of bias that might result from excluding missing cases. Missing values were assumed to be missing at random and missing values were imputed using chaining random forests [36] as implemented in missRanger (https://CRAN.R-project.org/package=missRanger). The random forest imputation is based on 200 trees and was controlled for a variation of additional predictors. A huge advantage of using random forests is that they produce a single imputed dataset, they are adaptive to interactions and nonlinear relationships not needing to specify an imputation model, and they can handle mixed types of missing data.

We conducted cross-lagged path models (CLPM) in a structural equation modeling framework to examine the temporal associations between parental concern of overeating and obesity indices (z-BMI and z-WC) (Fig. 2). Before CLPM, we regressed parental concern of overeating and obesity indices on age, sex, country, and the seven covariates listed above at each examination wave. We standardized the obtained residuals and used the standardized residuals in CLPM analyses. The effect of related siblings was taken into account by using survey methods with robust standard errors yielded by cluster variance estimators [37]. Goodness-of-fit criteria were a Comparative Fit Index (CFI), a Tucker-Lewis Index (TLI) ≥ 0.90, and a Root Mean Square Error of Approximation (RMSEA) close to 0.06 [38, 39].

To study whether obesity indices at baseline mediated the association between the PRS-BMI and parental concern of overeating at wave 3 (Fig. 3, Model 1) and whether parental concern of overeating at baseline mediated the association between the PRS-BMI and obesity indices at wave 3 (Fig. 3, Model 2), we conducted causal mediation analyses based on the counterfactual approach proposed by VanderWeele [40]. We used the med4way Stata package which allows for decomposing the total effect of an exposure to an outcome into four components: controlled direct effect, pure indirect effect (or mediation effect), reference interaction, and mediated interaction [41]. Sample bias correction estimates and 95% confidence intervals (CIs) were calculated using the bootstrapping approach of 1000 draws (Mediation Model 1) or 2000 draws (Mediation Model 2). We adjusted all models by baseline age, sex, country, the 7 covariates listed above and the first 12 genetic principal components. We a priori selected this set of confounders for which to control in the mediation models. We attempted to control for covariates that may be the cause of the exposure, or of the outcome, or of both [42]. The causal mediation analyses were further used to overcome the limitations of standard approaches, which do not account for mediator-outcome confounding. In case that there is mediator-outcome confounding, standard approaches are biased. For example, screen time, could be causally related to parental concern of overeating and to obesity indices, hence it is reasonable to include it as a confounder. We also conducted sensitivity analyses to examine the presence of unmeasured variables that may confound the relationship between the mediator and the outcome and to assess how strong are the assumptions about confounding that are needed to identify the direct and indirect effects [40].

Description of the analysis dataset

For the cross-lagged analyses, we included all children with data on eating behaviors and obesity indices at waves 1, 2 and 3 (n = 1848 from 1689 families, see Fig. 2A). We also performed the same analyses including children only from waves 1 and 3 (n = 2355 from 2143 families, see Fig. 2B). In mediation analyses, we included all children with genome-wide data, obesity indices at wave 1 and eating behaviors at wave 3 (Mediation Model 1, n = 1246 children) and with genome-wide data, eating behaviors at wave 1 and obesity indices at wave 3 (Mediation Model 2, n = 2386 children). Figure 1 is a flowchart showing the derivation of the sample sizes.

Results

General characteristics

Table 1 summarizes the general characteristics of the children at waves 1 and 3. About half were girls and half of the children came from families with a high parental education level. Children’s z-BMI, z-WC, and screen time increased with age. Table S2 (Additional files) further shows the baseline characteristics of the children with available data on eating behaviors who participated in all three study waves.

Table 1 General characteristics of the study sample (n = 2355)

Prospective associations between obesity indices and parental concern for overeating

Figure 2A shows the CLPMs of children’s parental concern of overeating and obesity indices considering all three available time points (waves 1, 2 and 3). The associations between z-WC to later parental concern of overeating were bi-directional, but the pathways from z-WC to later parental concern of overeating were stronger (β = 0.22, p < 0.001 from wave 1 to 2; β = 0.25, p < 0.001 from wave 2 to 3) than the reverse pathways (β = 0.11, p < 0.001 from waves 1 to 2; β = 0.10, p = 0.04 from waves 2 to 3). Similar associations from z-BMI to later parental concern of overeating were observed, but there were no significant associations in the opposite direction. The associations between z-BMI and z-WC to later parental concern of overeating were stronger than the other way around. The models had a reasonable fit to the data (RMSEA = 0.055, CFI = 0.996, TLI = 0.987 for z-BMI; and RMSEA = 0.078, CFI = 0.991, TLI = 0.967 for z-WC).

Figure 2B presents CLPMs of children’s parental concern of overeating and obesity indices at waves 1 and 3. We observed bi-directional associations between parental concern of overeating and both z-BMI and z-WC. The pathways from z-BMI and z-WC at wave 1 to later parental concern of overeating were stronger (β = 0.31, p < 0.001 and β = 0.28, p < 0.001, respectively) than the reverse pathways (β = 0.06, p = 0.006 for z-BMI; β = 0.09, p < 0.001 for z-WC). The models had a fit to the data (RMSEA = 0.000, CFI = 1.000, TLI = 1.000, for both models).

Fig. 2
figure 2

Cross-lagged path model for parental concern of overeating and obesity indices. Legend: Prospective associations between parental concern of overeating and z-BMI and z-WC. All models were adjusted for age, sex, country, parental education, parental income, well-being, screen time, playing outside, fruits and vegetables consumption and glycated hemoglobin. Model A n = 1848 for z-BMI; n = 1576 for z-WC; Model B n = 2355 for z-BMI; n = 2022 for z-WC. Abbreviations: z-score body mass index (z-BMI), z-score waist circumference (z-WC)

Results from the causal mediation analyses

The results obtained from the causal mediation analyses (Fig. 3) showed that z-BMI and z-WC at wave 1 mediated the association between the PRS-BMI and parental concern of overeating at wave 3 (Fig. 3, Model 1). Parental concern of overeating at wave 1 also mediated the association between the PRS-BMI and both obesity indices at wave 3, but to a lesser extent (Fig. 3, Model 2). The association between the PRS-BMI and parental concern of overeating at wave 3 was also partly due to the interaction with z-WC at wave 1; however, the proportion of the mediated interaction was weaker than the proportion of the mediation effect (9%, p = 0.01 and 42%, p < 0.001, respectively).

Fig. 3
figure 3

Results from the causal mediation models. Legend: Model 1: PRS-BMI (exposure), z-BMI and z-WC at wave 1 (mediators), parental concern of overeating at wave 3 (outcome). Model 2: PRS-BMI (exposure), parental concern of overeating at wave 1 (mediator), z-BMI and z-WC at wave 3 (outcomes). The indirect effect (or mediation effect) is represented by ab, and c represents the direct effect; total effect = c + ab. Results were presented as β-coefficients and their 95% CIs. All models were adjusted for baseline age, sex, country, parental education, parental income, well-being, screen time, playing outside, fruit and vegetable consumption, glycated hemoglobin and genetic principal components. aSample size smaller due to missing values on eating behaviors at wave 3. bSample size smaller due to missing values on z-WC. Abbreviations: polygenic risk score for body mass index (PRS-BMI), z-score body mass index (z-BMI), z-score waist circumference (z-WC)

The overall proportion mediated of the association between the PRS-BMI and parental concern of overeating at wave 3 (Mediation Model 1) was 44% (p < 0.001) and 51% (p < 0.001) through z-BMI and z-WC at wave 1, respectively. Lower mediation percentages were observed in Mediation Model 2, when parental concern of overeating at wave 1 mediated the association between PRS-BMI and both obesity indices at wave 3 (19%, p < 0.001 for z-BMI; 20%, p < 0.001 for z-WC). Sensitivity plots of the models are shown in Fig. S1 (Additional files).

Discussion

In this longitudinal analysis involving European children, we showed that BMI and WC predicted parental concern of overeating to a greater extent than vice versa. Further, we demonstrated that BMI and WC partly mediated the prospective association between the latest PRS-BMI and parental concern of overeating. We also showed that parental concern of overeating partly mediated the prospective association between the PRS-BMI and obesity indices, extending previous cross-sectional findings of eating behaviors as potential behavioral mediators of the genetic susceptibility to obesity.

Our study suggests a bi-directional association between parental concern of overeating and obesity. In particular, obesity indices predicted subsequent parental concern of overeating three to four times stronger than parental concern of overeating predicted subsequent obesity indices. Our results are largely consistent with previous studies of children using CLPMs to assess the directionality between eating behaviors and BMI [43, 44]. However, they used emotional overeating, which refers to eating more than needed when experiencing negative emotions [45], whereas our overeating construct was not related to emotional overeating but rather inquired about parental concern about their children eating too much and putting their child on a diet to prevent them from becoming overweight. In these previous studies, emotional overeating was bi-directionally associated with BMI during childhood, and the pathway from BMI to later emotional overeating was stronger. Similarly, in a study of Finnish adults, BMI and WC predicted greater changes in restrained eating than the other way around [46]. Additionally, another study of children aged 3 to 15 months showed bi-directional associations between appetite traits and weight in early infancy, but the pathway from appetite to subsequent weight was stronger than from weight to appetite [15], suggesting that the temporal associations between eating behaviors and obesity measures may differ during early infancy from those found in later childhood and adulthood.

The molecular and behavioral mechanisms [47] underlying most of the genetic variants associated with childhood obesity are still unknown [48]. The obesity-susceptibility gene FTO, which was discovered by a genome-wide association study [49], has been suggested to be involved in both metabolic and appetite pathways [50,51,52]. To our knowledge, previous research has not yet examined the temporal association between parental concern of overeating and obesity indices considering genetic susceptibility using PRSs. Previous mediation studies on the genetic susceptibility to obesity and eating behaviors have shown that high appetite and lower satiety responsiveness (behaviors similar to overeating) mediate the association between a PRS-BMI and obesity indices in early infancy and childhood [18, 19], suggesting that genes involved in appetite traits [53] may lead to weight gain. Nevertheless, negative reports also exist; for example, a prospective study of Norwegian children did not observe any mediation through appetite traits [54]. A previous study of monozygotic twin pairs discordant for BMI observed a reduced gray matter volume in the heavier twin [55]. These structural differences were involved in valuation and reward processes. Further, a recent study showed that obesity susceptibility genes may influence eating addiction and reward behaviors due to their strong expression in the insula and substantia nigra brain regions [56]. Hence, it seems that genes could affect appetite and satiety through the central nervous system leading to obesity, which further may lead to changes in the central nervous system [57].

Obesity susceptibility genes are not only expressed in the brain [53], but also in the adipose tissue [58], suggesting that there may be pathways other than eating behaviors by which obesity susceptibility genes might exert their influence. A recent American study of adults suggested physical activity, conscientiousness, education and depression as potential pathways by which genetic susceptibility to obesity may lead to weight gain [20]. As a novel aspect, we found that obesity indices mediate the prospective association between the PRS-BMI and parental concern of overeating. In line with this, a recent study of British adolescents showed that the association between different PRS-BMI and disordered eating behaviors was partly mediated through BMI [21]. Therefore, future research should examine pathways other than eating behaviors in the susceptibility to obesity, both on a behavioral but also on a physiological level.

We have to acknowledge some limitations of our study. First, causal mediation analysis relies on the assumption of no unmeasured confounding [40]. Our sensitivity analyses showed that unmeasured confounding might still exist, although we adjusted for a wide range of different confounders. Second, we used the PRS derived from Khera et al. [7] based on BMI data of adults and we applied it to children. The PRS from Khera is age-dependent and shows minimal associations with birthweight, but it was strongly associated with a gradient in weight that emerged in early childhood when the children were 3.5 years old [7]. This is similar to our children’s baseline age, which was 4.4 years old. Further, the use of a PRS with 2.1 million SNPs provided a more comprehensive assessment of genetic susceptibility than other PRSs used in previous studies.

Third, children did not participate in all study waves, and dropout at follow up in our sample is associated with overweight, lower well-being scores or lower parental education which may have led to an underestimation of the associations between overeating and obesity indicators [59]. In addition, parents of overweight children tend to underestimate the weight status of their children [60]. Fourth, parental report of children’s eating behavior may be biased due to concerns about child’s weight or social desirability [61, 62]. It is not surprising that parents who are concerned about their child’s weight are more likely to adopt changes in their child’s feeding practices [62], which may explain why obesity indices predicted subsequent parental concern of overeating more strongly. Future research is needed to carefully check the temporal relationships between obesity indices and objective measures of overeating. Further, the eating behavior questionnaire used in this study aimed to assess parental feeding practices related to overweight during early childhood [29]. In a previous study using the same study population at waves 1 and 2 [63] and the same eating behavior questionnaire [29], the Cronbach’s alpha for the construct “parental concern for overweight and overeating” showed good internal consistency (Cronbach’s alpha = 0.82) [62]. Furthermore, in this previous study the construct “parental concern for overweight or overeating” was related to overweight and obesity longitudinally [63]. Fifth, this study followed a modular approach and some data was only available for subgroups [24], hence some examination modules had fewer study participants resulting in smaller sample sizes. Finally, our sample consisted of European children and replication of our findings in non-European populations is warranted.

The strengths of the present study include its prospective design, which allowed for the examination of the temporal relationship between parental concern of overeating and obesity indices. The inclusion of three time points in our CLPM analyses provided information on eating behaviors and obesity indices when the children were 4, 6 and 10 years old, respectively. Moreover, our analyses were not limited to a single country as we included children from eight European countries.

In conclusion, the findings of this study suggest that the associations between parental concern of overeating and obesity indices are bi-directional, but the pathway from obesity to parental concern of overeating is about three- to four-fold stronger than vice versa. By bridging nutritional and genetic research, this study increases our understanding of the temporal relationships in children’s genetic susceptibility to overeating and obesity. Future prospective studies are needed to replicate and explain these findings in larger samples of culturally and genetically diverse populations of children.

Availability of data and materials

Data described in this study will be made available upon request from the corresponding author.

Abbreviations

BMI:

Body Mass Index

CFI:

Comparative Fit Index

CIs:

Confidence Intervals

CLPM:

Cross-lagged Path Model

ISCED:

International Standard Classification of Education

PRS:

Polygenic Risk Score

RMSEA:

Root Mean Square Error of Approximation

SD:

Standard Deviation

SNP:

Single Nucleotide Polymorphism

TLI:

Tucker-Lewis Index

WC:

Waist Circumference

z-BMI:

Body Mass Index z-score

z-WC:

Waist Circumference z-score

References

  1. Silventoinen K, Jelenkovic A, Sund R, Hur Y, Yokoyama Y, Honda C, et al. Genetic and environmental effects on body mass index from infancy to the onset of adulthood : an individual-based pooled analysis of 45 twin cohorts participating in the COllaborative project of development of anthropometrical measures in twins (CODATwins). Am J Clin Nutr. 2016;104:371–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Silventoinen K, Jelenkovic A, Sund R, Yokoyama Y, Hur Y, Cozen W, et al. Differences in genetic and environmental variation in adult BMI by sex , age , time period , and region : an individual-based pooled analysis of 40 twin cohorts. Am J Clin Nutr. 2017;106:457–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. 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. https://doi.org/10.1038/nature14177.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Yengo L, Sidorenko J, Kemper KE, Zheng Z, Wood AR, Weedon MN, et al. Meta-analysis of genome-wide association studies for height and body mass index in 700 000 individuals of European ancestry. Hum Mol Genet. 2018;27(20):3641–9. https://doi.org/10.1093/hmg/ddy271.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Silventoinen K, Konttinen H. Obesity and eating behavior from the perspective of twin and genetic research. Neurosci Biobehav Rev. 2020;109:150–65. https://doi.org/10.1016/j.neubiorev.2019.12.012.

    Article  CAS  PubMed  Google Scholar 

  6. Wainschtein P, Jain DP, Yengo L, Zheng Z, Group TopmAW, Consortium T-O for PM, et al. Recovery of trait heritability from whole genome sequence data Visscher 2019.pdf. Preprint. Available from: bioRxiv 588020.

  7. Khera AV, Chaffin M, Wade KH, Zahid S, Brancale J, Xia R, et al. Polygenic prediction of weight and obesity trajectories from birth to adulthood. Cell. 2019;177(3):587–96. https://doi.org/10.1016/j.cell.2019.03.028.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Lewis CM, Vassos E. Prospects for using risk scores in polygenic medicine. Genome Med. 2017;9(1):9–11. https://doi.org/10.1186/s13073-017-0489-y.

    Article  Google Scholar 

  9. Bouchard C, Tremblay A, Nadeau A, Després JP, Thériault G, Boulay MR, et al. Genetic effect in resting and exercise metabolic rates. Metabolism. 1989;38(4):364–70. https://doi.org/10.1016/0026-0495(89)90126-1.

    Article  CAS  PubMed  Google Scholar 

  10. Cai G, Cole SA, Butte NF, Voruganti VS, Comuzzie AG. Genome-wide scan revealed genetic loci for energy metabolism in Hispanic children and adolescents. Int J Obes. 2008;32(4):579–85. https://doi.org/10.1038/ijo.2008.20.

    Article  CAS  Google Scholar 

  11. Ravussin E, Lillioja S, Knowler W, Christin L, Freymond D, Abbott W, et al. Reduced rate of energy expenditure as a risk factor for body-weight gain. N Engl J Med. 1988;318(8):467–72.

    Article  CAS  PubMed  Google Scholar 

  12. Llewellyn CH, Van Jaarsveld CHM, Boniface D, Carnell S, Wardle J. Eating rate is a heritable phenotype related to weight in children. Am J Clin Nutr. 2008;88(6):1560–6. https://doi.org/10.3945/ajcn.2008.26175.

    Article  CAS  PubMed  Google Scholar 

  13. Carnell S, Haworth CMA, Plomin R, Wardle J. Genetic influence on appetite in children. Int J Obes. 2008;32(10):1468–73. https://doi.org/10.1038/ijo.2008.127.

    Article  CAS  Google Scholar 

  14. Llewellyn CH, Van Jaarsveld CHM, Johnson L, Carnell S, Wardle J. Nature and nurture in infant appetite: analysis of the Gemini twin birth cohort. Am J Clin Nutr. 2010;91(5):1172–9. https://doi.org/10.3945/ajcn.2009.28868.

    Article  CAS  PubMed  Google Scholar 

  15. Van Jaarsveld CHM, Llewellyn CH, Johnson L, Wardle J. Prospective associations between appetitive traits and weight gain in infancy. Am J Clin Nutr. 2011;94(6):1562–7. https://doi.org/10.3945/ajcn.111.015818.

    Article  CAS  PubMed  Google Scholar 

  16. Turcot V, Lu Y, Highland HM, Schurmann C, Justice AE, Fine RS, et al. Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity. Nat Genet. 2018;50(1):26–35. https://doi.org/10.1038/s41588-017-0011-x.

    Article  CAS  PubMed  Google Scholar 

  17. Calton MA, Vaisse C. Narrowing down the role of common variants in the genetic predisposition to obesity. Genome Med. 2009;1(3):1–4. https://doi.org/10.1186/gm31.

    Article  CAS  Google Scholar 

  18. Llewellyn CH, Trzaskowski M, Van Jaarsveld CHM, Plomin R, Wardle J. Satiety mechanisms in genetic risk of obesity. JAMA Pediatr. 2014;168(4):338–44. https://doi.org/10.1001/jamapediatrics.2013.4944.

    Article  PubMed  PubMed Central  Google Scholar 

  19. de Lauzon-Guillain B, Koudou YA, Botton J, Forhan A, Carles S, Pelloux V, et al. Association between genetic obesity susceptibility and mother-reported eating behaviour in children up to 5 years. Pediatr Obes. 2019;14(5):1–8. https://doi.org/10.1111/ijpo.12496.

    Article  Google Scholar 

  20. Stephan Y, Sutin AR, Luchetti M, Caille P, Terracciano A. An examination of potential mediators of the relationship between polygenic scores of BMI and waist circumference and phenotypic adiposity. Psychol Heal. 2020;35(9):1151–61. https://doi.org/10.1080/08870446.2020.1743839.

    Article  Google Scholar 

  21. Abdulkadir M, Herle M, De Stavola BL, Hübel C, Santos Ferreira DL, Loos RJF, et al. Polygenic score for body mass index is associated with disordered eating in a general population cohort. J Clin Med. 2020;9(4):1187. https://doi.org/10.3390/jcm9041187.

    Article  PubMed Central  Google Scholar 

  22. Masip G, Silventoinen K, Keski-Rahkonen A, Palviainen T, Sipilä PN, Kaprio J, et al. The genetic architecture of the association between eating behaviors and obesity: combining genetic twin modeling and polygenic risk scores. Am J Clin Nutr. 2020;112(4). https://doi.org/10.1093/ajcn/nqaa181.

  23. Herle M, Smith AD, Kininmonth A, Llewellyn C. The role of eating Behaviours in genetic susceptibility to obesity. Curr Obes Rep. 2020. https://doi.org/10.1007/s13679-020-00402-0.

  24. Ahrens W, Siani A, Adan R, De Henauw S, Eiben G, Gwozdz W, et al. Cohort Profile: The transition from childhood to adolescence in European children-how I.Family extends the IDEFICS cohort. Int J Epidemiol. 2017;46(5):1394–1395j. https://doi.org/10.1093/ije/dyw317.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Peplies J, Günther K, Gottlieb A, Lübke A, Bammann K, Ahrens W. Biological samples---standard operating procedures for collection, shipment, storage and documentation. In: Bammann K, Lissner L, Pigeot I, Ahrens W, editors. Instruments for health surveys in children and adolescents. Cham: Springer International Publishing; 2019. p. 57–76. https://doi.org/10.1007/978-3-319-98857-3_4.

    Chapter  Google Scholar 

  26. Bammann K, Peplies J, Marild S, Molnar D, Suling M, Siani A. Physical examinations. In: Bammann K, Lissner L, Pigeot I, Ahrens W, editors. Instruments for health surveys in children and adolescents. Cham: Springer International Publishing; 2019. p. 47–55. https://doi.org/10.1007/978-3-319-98857-3_6.

    Chapter  Google Scholar 

  27. Cole TJ, Lobstein T. Extended international (IOTF) body mass index cut-offs for thinness, overweight and obesity. Pediatr Obes. 2012;7(4):284–94. https://doi.org/10.1111/j.2047-6310.2012.00064.x.

    Article  CAS  PubMed  Google Scholar 

  28. Nagy P, Kovacs E, Moreno LA, Veidebaum T, Tornaritis M, Kourides Y, et al. Percentile reference values for anthropometric body composition indices in European children from the IDEFICS study. Int J Obes. 2014;38:S15–25. https://doi.org/10.1038/ijo.2014.131.

    Article  Google Scholar 

  29. Baughcum AE, Powers SW, Johnson SB, Chamberlin LA, Deeks CM, Jain A, et al. Maternal feeding practices and beliefs and their relationships to overweight in early childhood. J Dev Behav Pediatr. 2001;22(6):391–408. https://doi.org/10.1097/00004703-200112000-00007.

    Article  CAS  PubMed  Google Scholar 

  30. Weale ME. Quality Control for Genome-Wide Association Studies. In: Barnes M, Breen G, editors. Genetic Variation. Methods in Molecular Biology (Methods and Protocols). Totowa: Humana Press; 2010. https://doi.org/10.1007/978-1-60327-367-1_19.

    Chapter  Google Scholar 

  31. UNESCO. The international standard classification of education (ISCED). ISCED 2011. Montreal: UNESCO Institute for Statistics; 2012. https://doi.org/10.1007/BF02207511.

    Book  Google Scholar 

  32. Iguacel I, Michels N, Fernández-Alvira JM, Bammann K, De Henauw S, Felső R, et al. Associations between social vulnerabilities and psychosocial problems in European children. Results from the IDEFICS study. Eur Child Adolesc Psychiatry. 2017;26(9):1105–17. https://doi.org/10.1007/s00787-017-0998-7.

    Article  PubMed  Google Scholar 

  33. Huybrechts I, Börnhorst C, Pala V, Moreno LA, Barba G, Lissner L, et al. Evaluation of the children’s eating habits questionnaire used in the IDEFICS study by relating urinary calcium and potassium to milk consumption frequencies among european children. Int J Obes. 2011;35:S69–78. https://doi.org/10.1038/ijo.2011.37.

    Article  CAS  Google Scholar 

  34. Pala V, Reisch LA, Lissner L. Dietary behaviour in children, adolescents and families: the eating habits questionnaire (EHQ). In: Bammann K, Lissner L, Pigeot I, Ahrens W, editors. Instruments for health surveys in children and adolescents. Cham: Springer International Publishing; 2019. p. 103–33. https://doi.org/10.1007/978-3-319-98857-3_6.

    Chapter  Google Scholar 

  35. Peplies J, Günther K, Bammann K, Fraterman A, Russo P, Veidebaum T, et al. Influence of sample collection and preanalytical sample processing on the analyses of biological markers in the European multicentre study IDEFICS. Int J Obes. 2011;35:S104–12. https://doi.org/10.1038/ijo.2011.41.

    Article  CAS  Google Scholar 

  36. Stekhoven DJ, Bühlmann P. Missforest-non-parametric missing value imputation for mixed-type data. Bioinformatics. 2012;28(1):112–8. https://doi.org/10.1093/bioinformatics/btr597.

    Article  CAS  PubMed  Google Scholar 

  37. Williams RL. A note on robust variance estimation for cluster-correlated data. Biometrics. 2000;56:645–6. https://doi.org/10.1111/j.0006-341X.2000.00645.x.

    Article  CAS  PubMed  Google Scholar 

  38. Hu LT, Bentler PM. Cutoff criteria for fit indexes in covariance structure analysis: conventional criteria versus new alternatives. Struct Equ Model. 1999;6(1):1–55. https://doi.org/10.1080/10705519909540118.

    Article  Google Scholar 

  39. Bentler PM. Comparative fit indexes in structural models. Psychol Bull. 1990 Mar;107(2):238–46. https://doi.org/10.1037/0033-2909.107.2.238.

    Article  CAS  PubMed  Google Scholar 

  40. VanderWeele T. Explanation in causal inference: methods for mediation and interaction. New York: Oxford University Press; 2015.

    Google Scholar 

  41. Discacciati A, Bellavia A, Lee JJ, Mazumdar M, Valeri L. Med4way: a Stata command to investigate mediating and interactive mechanisms using the four-way effect decomposition. Int J Epidemiol. 2019;48(1):15–20. https://doi.org/10.1093/ije/dyy236.

    Article  Google Scholar 

  42. VanderWeele TJ. Principles of confounder selection. Eur J Epidemiol. 2019;34(3):211–9. https://doi.org/10.1007/s10654-019-00494-6.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Derks IPM, Sijbrands EJG, Wake M, Qureshi F, van der Ende J, Hillegers MHJ, et al. Eating behavior and body composition across childhood: a prospective cohort study. Int J Behav Nutr Phys Act. 2018;15(1):1–9. https://doi.org/10.1186/s12966-018-0725-x.

    Article  Google Scholar 

  44. Costa A, Severo M, Vilela S, Fildes A, Oliveira A. Bidirectional relationships between appetitive behaviours and body mass index in childhood: a cross-lagged analysis in the generation XXI birth cohort. Eur J Nutr. 2021;60(1):239–47. https://doi.org/10.1007/s00394-020-02238-9.

    Article  CAS  PubMed  Google Scholar 

  45. Wardle J, Guthrie CA, Sanderson S, Rapoport L. Development of the children’s eating behaviour questionnaire. J Child Psychol Psychiatry Allied Discip. 2001;42(7):963–70. https://doi.org/10.1111/1469-7610.00792.

    Article  CAS  Google Scholar 

  46. Konttinen H, Llewellyn C, Silventoinen K, Joensuu A, Männistö S, Salomaa V, et al. Genetic predisposition to obesity, restrained eating and changes in body weight : a population-based prospective study. Int J Obes. 2018;42(4):858–65. https://doi.org/10.1038/ijo.2017.278.

    Article  CAS  Google Scholar 

  47. Van Der Klaauw AA, Farooqi IS. The hunger genes: pathways to obesity. Cell. 2015;161(1):119–32. https://doi.org/10.1016/j.cell.2015.03.008.

    Article  CAS  PubMed  Google Scholar 

  48. Littleton SH, Berkowitz RI, Grant SFA. Genetic Determinants of Childhood Obesity. Mol Diagnosis Ther. 2020. https://doi.org/10.1007/s40291-020-00496-1.

  49. Loos RJF, Yeo GSH. The bigger picture of FTO - the first GWAS-identified obesity gene. Nat Rev Endocrinol. 2014;10(1):51–61. https://doi.org/10.1038/nrendo.2013.227.The.

    Article  CAS  PubMed  Google Scholar 

  50. Claussnitzer M, Dankel SN, Kim K-H, Quon G, Meuleman W, Haugen C, et al. FTO obesity variant circuitry and adipocyte Browning in humans. N Engl J Med. 2015;373(10):895–907. https://doi.org/10.1056/nejmoa1502214.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Hubáček JA, Pikhart H, Peasey A, Kubínová R, Bobák M. FTO variant, energy intake, physical activity and basal metabolic rate in caucasians. the HAPIEE study. Physiol Res. 2011;60(1):175–83. https://doi.org/10.33549/physiolres.932066.

    Article  PubMed  Google Scholar 

  52. Gerken T, Girard CA, Tung YCL, Webby CJ, Saudek V, Hewitson KS, et al. The obesity-associated FTO gene encodes a 2-oxoglutarate-dependent nucleic acid demethylase. Science. 2007;318(5855):1469–72. https://doi.org/10.1126/science.1151710.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Huvenne H, Dubern B, Clément K, Poitou C. Rare genetic forms of Obesity : clinical approach and current treatments in 2016. Obes Facts. 2016;9:158–73. https://doi.org/10.1159/000445061.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Steinsbekk S, Belsky D, Guzey IC, Wardle J, Wichstrom L. Polygenic risk, appetite traits, andweight gain in middle childhood a longitudinal study. JAMA Pediatr. 2016;170(2). https://doi.org/10.1001/jamapediatrics.2015.4472.

  55. Weise CM, Bachmann T, Pleger B. Brain structural differences in monozygotic twins discordant for body mass index. Neuroimage. 2019;201:116006. https://doi.org/10.1016/j.neuroimage.2019.07.019.

    Article  PubMed  Google Scholar 

  56. Ndiaye FK, Huyvaert M, Ortalli A, Canouil M, Lecoeur C, Verbanck M, et al. The expression of genes in top obesity-associated loci is enriched in insula and substantia nigra brain regions involved in addiction and reward. Int J Obes. 2020;44(2):539–43. https://doi.org/10.1038/s41366-019-0428-7.

    Article  CAS  Google Scholar 

  57. O’Brien PD, Hinder LM, Callaghan BC, Feldman EL. Neurological consequences of obesity. Lancet Neurol. 2017;16(6):465–77. https://doi.org/10.1016/S1474-4422(17)30084-4.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Ahn J, Wu H, Lee K. Integrative analysis revealing human adipose-specific genes and consolidating obesity loci. Sci Rep. 2019;9(1):1–13. https://doi.org/10.1038/s41598-019-39582-8.

    Article  CAS  Google Scholar 

  59. Hense S, Pohlabeln H, Michels N, Mårild S, Lissner L, Kovacs E, et al. Determinants of attrition to follow-up in a multicentre cohort study in children-results from the IDEFICS study. Epidemiol Res Int. 2013;2013:1–9. https://doi.org/10.1155/2013/936365.

    Article  Google Scholar 

  60. Nemecek D, Sebelefsky C, Woditschka A, Voitl P. Overweight in children and its perception by parents: cross-sectional observation in a general pediatric outpatient clinic. BMC Pediatr. 2017;17(1):1–10. https://doi.org/10.1186/s12887-017-0964-z.

    Article  Google Scholar 

  61. Shloim N, Edelson LR, Martin N, Hetherington MM. Parenting styles, feeding styles, feeding practices, and weight status in 4-12 year-old children: a systematic review of the literature. Front Psychol. 2015;6. https://doi.org/10.3389/fpsyg.2015.01849.

  62. Loth KA, Mohamed N, Trofholz A, Tate A, Berge JM. Associations between parental perception of- and concern about-child weight and use of specific food-related parenting practices. Appetite. 2021;160:105068. https://doi.org/10.1016/j.appet.2020.105068.

    Article  PubMed  Google Scholar 

  63. Bammann K, Gwozdz W, Pischke C, Eiben G, Fernandez-Alvira JM, De Henauw S, et al. The impact of familial, behavioural and psychosocial factors on the SES gradient for childhood overweight in Europe. A longitudinal study. Int J Obes. 2017;41(1):54–60. https://doi.org/10.1038/ijo.2016.137.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

For the valuable help with the language editing, we want to thank MPH Alyce Whipp from Language Services, University of Helsinki. We thank the children and parents who participated in the examinations of the IDEFICS/I.Family cohort study.

Funding

This work was done as part of the IDEFICS (http://www.idefics.eu) and I.Family studies (http://www.ifamilystudy.eu/). We gratefully acknowledge the financial support of the European Commission within the Sixth RTD Framework Programme Contract No. 016181 (FOOD), and the Seventh RTD Framework Programme Contract No. 266044. GM has received funding from the University of Helsinki, Doctoral School in Health Sciences of the Doctoral Programme in Population Health. LHB acknowledges funding from the Austrian Science Fund (FWF): M 2449. JK has been supported by the Academy of Finland (grant 312073). The funding sources had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Author information

Authors and Affiliations

Authors

Consortia

Contributions

GM 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. Drafting of the manuscript: GM; Statistical analysis: GM and RF; Study concept and design: LHB, JK, RA and GM; Study supervision: JK, LHB, AK-R, KS; Obtained funding: WA, LL, LM, DM, IP, TV, SDH and JK; Acquisition, analysis, or interpretation of data: all authors; Critical revision of the manuscript for important intellectual content: all authors. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Guiomar Masip.

Ethics declarations

Ethics approval and consent to participate

We certify that all applicable institutional and governmental regulations concerning the ethical use of human volunteers were followed during this research. Ethics approval was obtained from responsible committees in each country in accordance with the ethical standards of the 1964 Declaration of Helsinki and its later amendments (Belgium: Ethics Committee of the Gent University Hospital, 15/10/2007, ref.: No. EC UZG 2007/243 and 19/02/2013, No. B670201316342. Cyprus: Cyprus National Bioethics Committee, 12/07/2007, ref.: No. EEBK/EM/2007/16 and 21/Feb/2013,No.EEBK/ETI/2012/33. Estonia: Tallinn Medical Research Ethics Committee (TMREC), 14/06/2007, ref.: No. 1093 and 17/January 2013, No. 128. Germany: Ethic Commission of the University of Bremen, 16/01/2007 and 11/12/2012. Hungary: Medical Research Council, 21/Jun/2007, ref.: 22–156/2007-1018EKU and 18/12/2012, 4536/2013/EKU. Italy: Ethics Committee of the Local Health Authority (ASL) in Avellino, 19/06/2007, ref.: No. 2/CE and 18/Sep/2012, No. 12/12. Spain: Ethics Committee for Clinical Research of Aragon (CEICA), 20/06/2007, ref.: No. PI07/13 and 13/Feb/2013, No. PI13/0012. Sweden: Regional Ethics Research Board in Gothenburg, 30/07/2007, ref.: No. 264–07 and 10/Jan/2013, No. 927–12). All parents or legal guardians of the participating minor children gave written informed consent. Children older than 16 years provided written informed consent, while children aged 12 and over gave a simplified written consent. Younger children gave oral consent for examinations and sample collection.

Consent for publication

Not applicable.

Competing interests

The authors have no conflicts of interest to declare.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table S1.

Eating behavior patterns and factor loadings in varimax-rotated principal components. Table S2. Baseline characteristics of the study sample of the three study waves (n = 1848). Figure S1. Sensitivity plots from causal mediation analyses.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Masip, G., Foraita, R., Silventoinen, K. et al. The temporal relationship between parental concern of overeating and childhood obesity considering genetic susceptibility: longitudinal results from the IDEFICS/I.Family study. Int J Behav Nutr Phys Act 18, 139 (2021). https://doi.org/10.1186/s12966-021-01205-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12966-021-01205-9

Keywords