Genetic Polymorphisms, Mediterranean Diet and Microbiota-Associated Urolithin Metabotypes can Predict Obesity in Childhood-Adolescence.

Environmental and genetic factors are associated with pandemic obesity since childhood. However, the association of overweight-obesity with these factors, acting as a consortium, has been scarcely studied in children. We aimed here to assess the probabilities of being overweighed-obese in a randomly recruited cohort of Spanish children and adolescents (n = 415, 5−17 years-old) by estimating the odds ratios for different predictor variables, and their relative importance in the prediction. The predictor variables were ethnicity, age, sex, adherence to the Mediterranean diet (KIDMED), physical activity, urolithin metabotypes (UM-A, UM-B and UM-0) as biomarkers of the gut microbiota, and 53 single-nucleotide polymorphisms (SNPs) from 43 genes mainly related to obesity and cardiometabolic diseases. A proportional-odds logistic ordinal regression, validated through bootstrap, was used to model the data. While every variable was not independently associated with overweight-obesity, however, the ordinal logistic model revealed that overweight-obesity prevalence was related to being a young boy with either UM-B or UM-0, low KIDMED score and high contribution of a consortium of 24 SNPs, being rs1801253-ADRB1, rs4343-ACE, rs8061518-FTO, rs1130864-CRP, rs659366-UCP2, rs6131-SELP, rs12535708-LEP, rs1501299-ADIPOQ, rs708272-CETP and rs2241766-ADIPOQ the top-ten contributing SNPs. Additional research should confirm and complete this model by including dietary interventions and the individuals’ gut microbiota composition.


Methods
Study population. This research ('The PolyMicroBio study') was included in the Spanish National Project AGL2015-64124-R and complied with the ethical guidelines outlined in the Declaration of Helsinki and ethical principles for medical research involving human subjects (Seoul, Korea, 2008). The study was conceived to stratify the participants according to their urolithin metabotypes after three days of walnuts or pomegranate juice consumption 26 and was not intended to modify any variable in the children. The trial was registered at clinicaltrials.gov (NCT03318042), and the Spanish National Research Council's Bioethics Committee (Madrid, Spain) approved the protocol. Inclusion criteria were ages from 5 to 17 years old and good health status. Exclusion criteria were diagnosed pathology, previous gastrointestinal surgery, chronic medication and antibiotic intake one month before participating. A total of 415 children and adolescents were randomly recruited. Children within the 5 to 12 years old group (n = 202) were recruited from the public primary school 'CEIP Jara Carrillo' (Alcantarilla, Murcia, Spain) and adolescents aged from 13 to 17 (n = 213) from the public high school 'IES Alcántara' (Alcantarilla, Murcia, Spain). Parents were fully informed and gave their written informed consent before the participation of all students.
Urolithin metabotypes. Children and adolescents consumed 25 g peeled raw walnuts daily or 250 mL of pomegranate juice daily (in the case of individuals allergic to nuts) for three days. Packs of peeled walnuts were kindly provided by Borges International Group, S.L. (Reus, Tarragona, Spain) and pomegranate juice by the AMC Group (Espinardo, Murcia, Spain). In the morning of the fourth day, a sample of urine was provided for its analysis by high-performance liquid chromatography with diode array detection coupled to electrospray ionisation and ion-trap tandem mass spectrometry (HPLC-DAD-ESI-IT-MS/MS), and ultra-high performance liquid chromatography coupled with electrospray ionization-quadrupole-time-of-flight-mass spectrometry (UPLC-ESI-QTOF-MS) as described elsewhere 27 . This allowed the stratification of the participants according to their different capacity to metabolise ellagic acid derivatives into urolithins, i.e., urolithin metabotypes UM-A, UM-B or UM-0 as previously described 26 . Anthropometric measurements and validated questionnaires. The determinations of height, weight, and waist and hip circumference were performed always by the same research staff, using the same equipment in all cases, and in the presence of teachers from the educational centres. The child growth standards from the World Health Organization (WHO) were used to define the BMI (kg/m 2 ) cut-offs for underweight, normoweight, overweight and obese individuals as a function of sex and age 5 . The students were asked to record possible incidences (medication, protocol compliance, etc.), and also their physical activity level 28 , which took into account the two hours of physical activity a week in their schools (low activity) and the practice of additional extracurricular sports at least three days a week (high activity). Besides, a validated questionnaire to assess the adherence to the Mediterranean diet in children (KIDMED) was used 29 . The score in this questionnaire (ranging www.nature.com/scientificreports www.nature.com/scientificreports/ from 1 to 13, from very poor to optimum adhesion, respectively) was grouped as 'Low' (score from 1 to 4), 'Medium' (from 5 to 8), and 'Good' (from 9 to 13). Candidate genes and polymorphisms were identified after browsing   the Single Nucleotide Polymorphism Database (dbSNP) and examining the published literature regarding each  known gene and variant (favourable and unfavourable) associations [30][31][32] , and analysing their potential regulatory and biological functions with the Regulome DB and HaploReg v4.1 databases 33,34 . On the same day of the anthropometric evaluation, saliva samples were obtained by gently rubbing the inside part of the cheek with a sterile swab, free of human RNA and DNA (Deltalab, Barcelona, Spain). Children were asked to clean their mouths and avoid eating or drinking 60 min before collection of samples to prevent contaminations. Two samples were obtained per student. The swabs were immediately stored in refrigeration and further frozen at −80 °C until their processing. Genomic DNA extraction and genotyping were carried out at the GENYAL Platform (IMDEA-Food, Madrid, Spain) using the OpenArray ™ AccuFill ™ System (Life Technologies Inc. Carlsbad, CA, USA) as described elsewhere 35 . Data analysis was made by TaqMan Genotyper Software v1.3 (autocaller confidence level >90%).

Statistical analysis.
A proportional-odds logistic ordinal regression was used to model the data with the software R version 3.5.1 (www.r-project.org). Nine subjects with >40% missing SNPs were removed, resulting in a final sample size of 406 students ( Supplementary Fig. 2). The missing data was singly imputed using the missForest R package. Redundant predictors (rs9928094-FTO, rs9935401-FTO) were identified and removed using the Hmisc R package. SNPs with either favourable or unfavourable genotype frequencies below 5% were also removed (rs4994-ADRB3, rs7913948-ALOX5, rs7412-APOE, rs328-LPL, rs16139-NPY, rs6008259-PPARγ, rs2066826-PTGS2) (Supplementary Table 1). Using the 'n/15 rule' 36 , only 21 predictors could be used to get reliable estimates. Consequently, data reduction was applied to the SNPs by applying Multiple Correspondence Analysis (MCA) and using only the first 15 MCA dimensions. Ethnic groups representing less than 1% each were merged into the 'Other' category.
Wald tests for all the predictors in the model were generated, and they were further ranked by importance based on the χ 2 -degrees of freedom (-df) score. A simplified model was obtained by applying a 'fast-backwards' variable elimination approach 37 based on the Akaike's Information Criterion (AIC) 38 . The approximate βs and ORs (and their 95% confidence intervals, CI) of the remaining variables were reported. The full model was validated through bootstrap to provide estimates of the performance of the model in new data in comparison with the training data, in the form of Sommers D xy , R 2 , intercept, slope, E max , and Briers B score. The significant contributions (coefficients of determination R 2 with p values < 0.05) of SNPs to the essential MCA dimension were plotted to deconvolute it.
We used the proportional odds assumption in the model. 'Physical activity' , 'KIDMED' and 'Urolithin metabotype' showed some deviation from this assumption. However, alternative extended continuation ratio models with this assumption relaxed for these variables did not result in better models according to the AIC, i.e., the higher complexity of the model was not compensated by the increase in the fit. The inclusion of transformations of some predictors, including a restricted cubic spline for both ' Age' and 'KIDMED' , did not result in improved models, as judged by the AIC. Finally, a genetic score (computed as the sum of risk alleles) was also tested as a possible surrogate for the SNPs variables, but this did not result in a better replacement for the MCA dimensions. All tests were bilateral, with a significance level of 0.05.
Other statistical analyses were carried out using the SPSS software, v23.0 (SPSS Inc., Chicago, IL, USA). When more than two groups were compared, analyses of variance (ANOVA), followed by Bonferroni-corrected t-test (for post-hoc analysis) or the Kruskal-Wallis followed by Dunn's test were used for normally and non-normally distributed data, respectively (KIDMED score vs FTO genotype TT, AT or AA, etc.). Comparison of non-normally distributed quantitative variables between two clusters was approached using the Mann-Whitney U-test (FTO TT genotype vs BMI or waist, etc.). Comparison of categorical variables was assessed using the Pearson's χ 2 test. Spearman's rank or Pearson correlations were applied to explore possible associations between variables (BMI vs hip/height, etc.). Plots of data were performed using Sigma Plot 13.0 (Systat Software, San Jose, CA, USA). Table 1 shows characteristics of the cohort as a function of age and sex, including anthropometric values (hip, waist, weight and BMI), the distribution of urolithin metabotypes (A, B and 0), the KIDMED scores (grouped as low, medium and good adherence), physical activity and the percentage of normoweight, overweight and obesity. The participants were mainly Caucasian-Europeans (93.5%), with a small proportion of Arabs (2.9%) and Amerindians (2.2%), and a marginal presence of Black-Africans (0.96%), Asian-Chinese (0.22%), and Indo-Aryans (0.22%).

characteristics of the cohort and associations with overweight-obesity.
The hip-to-height ratio was the best anthropometric index associated with BMI (r = 0.78, p = 1.3 × 10 -86 ) vs the waist-to-hip ratio (r = 0.13, p = 0.007) and the waist-to-height ratio (r = 0.32, p = 0.001) ( Supplementary  Fig. 3). As expected in growing children, BMI values increased on average from 5 to 17 years (Fig. 1A). The percentage of overweight-obesity decreased from 5 to 17 years (from 50% to 25%, respectively) with the exemption of boys from 9 to 12 years old (n = 64) that reached the highest prevalence of overweight-obesity (∼70%) (Fig. 1B). This categorized representation (Fig. 1B) was coincident with the quantitative evolution of BMI z-scores for girls ( Fig. 1C) and boys (Fig. 1D) as a function of age.
Regarding the SNPs analysed, after correcting for multiple tests, three SNPs did not satisfy the Hardy-Weinberg equilibrium (HWE) (rs1801253, rs5082, rs11868035) (Supplementary Table 1). The rest of the SNPs were in equilibrium and were close to European frequencies. In the case of these three SNPs, we can speculate that the reason could be the association established by the presence of certain consanguinity (several sibling groups in the      www.nature.com/scientificreports www.nature.com/scientificreports/ cohort). Nevertheless, this did not affect our results, since we aimed to compare variables from different domains and not only in estimating the particular effect of a single SNP. Overall, it is not absolutely necessary to have HWE in our approach, i.e., to estimate odds ratios of an MCA dimension and rank the predictors.
The distribution of urolithin metabotypes in this cohort (Fig. 1E) revealed that both a lower prevalence of UM-A and a higher occurrence of UM-0 were associated with an increased percentage of overweight-obesity after bivariate analysis (p = 0.015) (Fig. 1F).
No significant association between physical activity and overweight-obesity distribution was found (results not shown). On the contrary, there were many overweight-obese children with high physical activity, which should be explained as a consequence of their overweight-obesity status (results not shown). Regarding the diet, the KIDMED scores showed mean values of 6.9 ± 2.1 for the entire cohort (Supplementary Table 2) and ranged from the lowest value of 5.7 ± 1.4 in 15-year-old girls to the highest value of 8.8 ± 2.2 in 6-year-old boys ( Fig. 2A). However, no significant differences were found between boys and girls as well as through the range of age ( Fig. 2A). Besides, no significant association was found between overweight-obesity and KIDMED in this group (results not shown).
We next explored the potential association of both adherence to the Mediterranean diet and percentage of overweight-obesity with the well-known obesity risk allele A of the rs9939609 SNP in the FTO gene. Figure 2B shows the KIDMED scores and the percentage of normoweight and overweight-obesity in children depending on their genotype TT, AT or AA. Remarkably, the adherence to the Mediterranean diet was significantly lower in children with the risk-associated genotypes AT and AA vs the TT genotype (Fig. 2B). Although the proportion of overweight-obesity was the highest for the AA genotype (45.5%), however, the difference vs that of the TT genotype (38.7%) did not reach statistical significance (Fig. 2B). Once again, many exemptions prevented the usefulness of the combination of KIDMED scores and rs9939609 genotypes as unique predictors of overweight-obesity in this cohort.
Therefore, although both the urolithin metabotypes and rs9939609 SNP-FTO could partially contribute to the overweight-obesity distribution as independent predictor variables in this cohort of children and adolescents; however, all the possible SNP-SNP interactions together with the rest of variables had not been taken into account. Therefore, we next developed an ordinal logistic model to identify the consortium of variables that could estimate the odds ratios of the overweight-obesity distribution in this cohort. An ordinal logistic model to identify the consortium of variables associated with overweight-obesity. In this model, we used as the ordinal response the normoweight, overweight and obesity classification for children, based on sex and age to estimate the odds ratios (and the corresponding 95% confidence intervals) for different predictor variables, as well as their relative importance in the prediction of the response. In this holistic approach, we used as predictors the 'Ethnicity' , 'Urolithin metabotypes' , 'KIDMED score' , 'Physical activity' , and genetic polymorphisms (44 SNPs were finally included and compressed into 15 MCA dimensions, hereafter termed 'SNP.Dim. '), together with sex and age. Although the distribution of normoweight, overweight and obesity WHO-based categories was apparently adjusted by sex and age, we still observed a trend for decreasing average age when moving from normoweight, overweight and obese children, as well as enrichment in boys in the same order. Therefore, we also included these variables in the model. Figure 3 displays the predictors used in the ordinal logistic model ranked by their apparent importance, as measured by their χ 2 -df score. The model was highly significant (p < 0.0001), yielding the components ' Age' , and 'SNP.Dim.14' as the two most important predictors. Therefore, in the prediction of overweight-obesity in this population, apart from ' Age' and 'Sex' , the most critical contributing variables were SNPs (through the variable SNP.Dim.14), followed by 'Ethnicity' , 'Urolithin metabotype' and 'KIDMED' . The 'Physical activity' seemed to be irrelevant in this sample.
We next applied a fast-backwards approach based on AIC to obtain a reduced model and estimate the corresponding βs and ORs (as well as their 95% CI). The variables kept were 'Sex' , ' Age' , 'Urolithin metabotype' , 'KIDMED' , and the genetic components 'SNP.Dim.3' , 'SNP.Dim.11' , and 'SNP.Dim.14' ( Table 2). The variable 'Ethnicity' was not retained in the simplified model, probably due to its high complexity (4 levels). Our results reveal that being a boy either with UM-B or UM-0 and having a higher SNP.Dimension.14, all increased the chances of overweight-obesity in this study population. On the contrary, ageing, better adherence to the Mediterranean diet (higher KIDMED score), and being UM-A was associated with lower probabilities of being overweighed-obese. Remarkably, all these variables operated additively to build the final probability of overweight-obesity for each subject.  www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 4 shows the R 2 values of the 24 SNPs that significantly contributed to SNP.Dim.14 (p < 0.05). The most contributing SNP to this consortium was rs1801253-ADRB1, as well as other obesity risk-associated SNPs (UCP2, ADIPOQ, LEP, MC4R, etc.). However, there were other SNPs with less known involvement in obesity. Table 3 shows all the SNPs contributing to SNP.Dim.14 with their definitions and main processes in which they are involved.

Model validation.
We used bootstrap to validate the model ( Table 4). The model showed some degree of overfitting, reflected in a decrease of the indexes after correction for optimism. The discriminative capacity of the model was modest, although not negligible according to the Dxy, R 2 indexes (it must be taken into account that the typical values of R 2 in ordinal models are much lower than those observed in linear regression models), and the Brier score of 0.23. The calibration showed some degree of shrinkage, as seen from the deviation of the     www.nature.com/scientificreports www.nature.com/scientificreports/ intercept and slope from 0 and 1, respectively. E max , the maximum calibration error in predicting p (Y > normal), was 0.11, showing some degree of miscalibration. Overall, we obtained a significant model that still showed predictive power in external data.

Discussion
The present study shows that the probabilities of childhood overweight-obesity cannot be explained by specific isolated variables (that is, one or a few SNPs, or just physical activity, or diet, etc.), but by complex, multifactorial associations of environmental and genetic components. Besides, a genetic score, widely used in the literature 39 , and understood as the theoretical equal-additive contribution of every single SNP could not explain the distribution of overweight-obesity in our cohort. In this regard, our analysis adds significance to previous studies that have associated the genetic background of adult individuals with the adherence to the Mediterranean diet and some markers related to obesity and metabolic syndrome [39][40][41] . However, to the best of our knowledge, there are no previous studies that explore the occurrence of overweight-obesity in children and adolescents by estimating the odds ratios for the predictor variables 'urolithin metabotypes' as gut microbiota biomarkers, 'age' , 'sex' , 'adherence to the Mediterranean diet' , and an identified consortium of 24 SNPs from 22 genes, mainly related to obesity and cardiometabolic diseases.
The purpose of the ordinal logistic regression model was to estimate the odds ratios of the different explanatory variables and rank them by importance, rather than a tool to predict the probability of overweight and obesity of new individuals. There are in the literature several previous reports that model the probabilities of overweight and obesity in children through ordinal regression models [42][43][44][45][46] . However, to the best of our knowledge, this is the first time that this type of model is used with such a wide set of predictors comprising different putatively influencing domains (genetic, diet, exercise, microbiome urolithin metabotype, and ethnicity, in addition to sex and age) and ranks them according to their relative importance. The χ 2 -df score returned the following importance ranking: age >SNP.Dim.14 >ethnicity >sex >urolithin metabotype >KIDMED; while physical activity seemed to have a negligible importance. The use of an AIC-based fast-backwards reduction of the model removed 'Ethnicity' from the final model, probably because of the large number of categories in this predictor 42 , that would be highly penalised by the AIC criterion. Besides, our sample was mainly of Caucasian-European origin (93.5%), and it is expected that in a sample with a more balanced distribution of ethnicity, this predictor would have high importance. Furthermore, the modest predictive capability of the model suggested the need for additional predictors, which is not unexpected, given the multifaceted aetiology of obesity. Still, the model, derived for estimating purposes, allowed us to rank predictors of different domains by their importance, as well as to estimate odds ratios for them.
The classification of overweight and obesity in adults is rather simple and independent of age and sex (i.e. 25 ≤ BMI < 30 for overweight, and BMI ≥ 30 for obesity). However, the WHO-based definition for overweight-obesity in children and teens uses age-and sex-based percentiles 5 , and a higher BMI does not necessarily correspond with overweight-obesity in growing children. Therefore, using overweight and obesity as response variables instead of BMI values is especially useful in children and adolescents. This is due to the difficulty in comparing them at different ages and of the two sexes, given the remarkable, sex-dependent change in weight and height, during this period of development (i.e., while our model predicted that overweight-obesity decreases upon ageing, the use of BMI yielded the opposite result, masking the real result). Still, sex and age were essential predictors in our model, with estimated ORs of 1.67 for boys, and 0.88 for a one-year increase, respectively, which means an increased percentage of overweighed-obese boys, and a reduction of overweight-obesity upon ageing.
The gut microbiota contributes to the pathophysiology of obesity 3 , and a recent report shows the role of gut microbial metabolites in the expression of the microRNA-181 family, which regulates white adipose tissue inflammation and obesity in children 14 . However, the gut microbiota, as a predictor of childhood obesity, has been scarcely approached and not usually considered together with SNPs and other variables. In a targeted approach,  50 . Specifically, the genus Collinsella was more abundant in obese adolescents, a microbial group that could be related to the endothelial dysfunction 50,51 . Overall, all these results are remarkable since we recently reported that UM-B was enriched in the Coriobacteriaceae family, which was positively correlated with blood total-cholesterol, LDL-cholesterol, and BMI in adults 25 Table 4. Validation of the model using Bootstrap. *D xy = bias-corrected Somers rank correlation coefficient that goes from −1 to 1; R 2 = Nagelkerke R 2 that goes from 0 to 1; Intercept and Slope of a logistic calibration equation (should be 0 and 1, respectively, for a perfect fit); E max = maximum calibration error for p (Y = 0) based on the linear-logistic recalibration; B = Brier's quadratic probability score, which goes from 0 (the best score) to 1 (the worst score). (2020) 10:7850 | https://doi.org/10.1038/s41598-020-64833-4 www.nature.com/scientificreports www.nature.com/scientificreports/ we had previously observed a trend between UM-B occurrence and overweight-obesity in adult individuals, however, we could not establish a definitive link probably due to the lack of other interacting factors such as those included in the present study (SNPs, diet, etc.) 26,52 . Therefore, the Coriobacteriaceae family, and probably the Proteobacteria phylum, more abundant in obese children as well as in UM-B, could be the link between UM-B occurrence and overweight-obesity, which suggests that the microbiota associated with urolithin metabotypes could contribute in the prediction of the probability of being overweighed or obese. Regarding UM-0, its occurrence in the population is approximately constant (~10%), although there is both a higher occurrence and variability in the childhood 26 . The microbiota associated with UM-0 has been reported to show lower diversity than UM-B and UM-A, which could be indicative of an obesity-prone microbiota 4 . Nevertheless, we cannot exclude a possible shift of metabotype in children from UM-0 to either UM-A or UM-B determined by ageing 26 or after ellagitannin-rich diets as previously described in adults 21 .
Regarding the possible health implications of belonging to one or another metabotype, we have recently reported that urolithin metabotypes determined the different restoration capacity of the gut microbiota and the anthropometric values (weight, waist and hip) of healthy women up to 12 months after delivery 53 . The gut microbiota of pregnant women is in dysbiosis, which persists at least 1 month after delivery. We observed that the gut microbiota associated with UM-B was more resilient than that of UM-A, which would have negative implications in the dysbiotic-prone UM-B. In contrast, the gut microbiota of UM-A women progressively became normal during the year after childbirth. Therefore, we suggested that the determination of urolithin metabotypes in pregnant and lactating women could be a useful tool to predict their predisposition to the recovery of the gut microbiota and anthropometric values, significantly altered during pregnancy and after childbirth 53 .
Many studies describe the association, or lack of association, of specific SNPs with obesity. The rs9939609 SNP-FTO has been reported to confer a predisposition to obesity by regulating the control of food intake and food choice, suggesting a link to a hyperphagic phenotype or a preference for energy-dense foods in Scottish children 54 . This agrees with the connection between the low adherence to the Mediterranean diet and the risk allele A of rs9939609 SNP-FTO in our cohort. However, we did not observe a clear association of this SNP with overweight-obesity, but only a trend with many exemptions, in agreement with other studies where some SNPs, previously reported to be involved in obesity, such as the rs17782313-MC4R and rs9939609-FTO, exerted weak effects and very scarce contribution to obesity in 773 pre-pubertal Portuguese children 55 . In the present study, a consortium of 24 SNPs was identified as the second contributing predictor to overweight-obesity in our cohort. The rs1801253 SNP (also called Arg389Gly), located in the gen ADRB1, was the most contributing SNP within this consortium. It is known the involvement of the ADRB1 gene polymorphisms in cardiovascular diseases 56 but also in obesity 57,58 .
The ADRB1 gene codifies a G-coupled protein (the β1-adrenergic receptor) that binds the catecholamines epinephrine and norepinephrine and controls sympathetic responses in the heart, kidney and adipocytes. Interestingly, Dionne et al. reported that the rs1801253 SNP-ADRB1 was associated with higher body weight and BMI in a cohort of Caucasian women (n = 931) 57 . In the same line, Aradillas-García et al. reported that the rs1801253 SNP-ADRB1, but not the Trp64Arg ADRB3, was associated with obesity in Mexican children 59 . The connection between ADBR1 SNPs and obesity could rely on the catecholamines, which are considered significant lipolysis regulators 60 and affect differentiation and proliferation of adipocytes 61 . In this regard, Lee et al. 62 established the association between impaired urinary epinephrine and norepinephrine excretion and obesity, insulin resistance, and metabolic syndrome in a cohort of 577 Chinese subjects.
Overall, all the above highlights again the need to consider SNPs consortia, interacting with other variables, instead of few SNPs in those studies aimed to associate SNPs with obesity. We are aware that the present study is an exploratory validation for a proof-of-concept, i.e., an ordinal logistic model that associates child overweight-obesity with a consortium of SNPs potentially interacting with the urolithin metabotypes-associated microbiota, adherence to the Mediterranean diet, age, and sex. Although we claim for the rationale of our approach and its potential usefulness, however, our results should be confirmed with additional research. We also acknowledge some limitations that should be considered in further studies, which also could improve its prediction capability. For example, it would be interesting to include other possible variables, such as the detailed composition and functionality of the individuals' gut microbiomes, and dietary interventions to evaluate not only associations but also individuals' responses. The latter would be even better than the use of validated questionnaires. Besides, a higher number of SNPs (or many SNPs associated with a specific gene) should be explored, especially in children from other geographical origins and ethnicities. Also, the inclusion of serobiochemical variables and traits related to obesity and its comorbidities (blood lipid profile, blood pressure, glucose homeostasis, etc.) could yield relevant information. Finally, a validation cohort (i.e., a parallel-group with all the children either normoweight or obese) should also be considered in further studies to confirm our model fully.

conclusions
The present research highlights the need for a holistic approach to unravel the predictors of overweight-obesity in children. Our results confirm, in agreement with the multifaceted aetiology of obesity, the link of childhood overweight-obesity to multifactorial associations of environmental and genetic components. The ordinal logistic model revealed that child overweight-obesity prevalence was related to being a young boy with either UM-B or UM-0, low KIDMED score and high contribution of a consortium of 24 SNPs, being rs1801253-ADRB1, rs4343-ACE, rs8061518-FTO, rs1130864-CRP, rs659366-UCP2, rs6131-SELP, rs12535708-LEP, rs1501299-ADIPOQ, rs708272-CETP and rs2241766-ADIPOQ the top-ten contributing SNPs. Therefore, it is of particular relevance the evaluation of interactive SNPs consortia along with the stratification of the children according to their urolithin metabotypes, which could be early biomarkers, in the case of UM-B and UM-0, of a dysbiotic-prone obesity-associated microbiota.