Urinary metabolic biomarkers of diet quality in European children are associated with metabolic health
======================================================================================================

* Nikolaos Stratakis
* Alexandros P. Siskos
* Eleni Papadopoulou
* Anh N. Nguyen
* Yinqi Zhao
* Katerina Margetaki
* Chung-Ho E. Lau
* Muireann Coen
* Lea Maitre
* Silvia Fernández-Barrés
* Lydiane Agier
* Sandra Andrusaityte
* Xavier Basagaña
* Anne Lise Brantsaeter
* Maribel Casas
* Serena Fossati
* Regina Grazuleviciene
* Barbara Heude
* Rosemary R C McEachan
* Helle Margrete Meltzer
* Christopher Millett
* Fernanda Rauber
* Oliver Robinson
* Theano Roumeliotaki
* Eva Borràs
* Eduard Sabidó
* Jose Urquiza
* Marina Vafeiadi
* Paolo Vineis
* Trudy Voortman
* John Wright
* David V. Conti
* Martine Vrijheid
* Hector C. Keun
* Leda Chatzi

## Abstract

Urinary metabolic profiling is a promising powerful tool to reflect dietary intake and can help understand metabolic alterations in response to diet quality. Here, we used 1H-NMR spectroscopy in a multi-country study in European children (1147 children from 6 different cohorts) and identified a common panel of 4 urinary metabolites (hippurate, *N*-methylnicotinic acid, urea and sucrose) that was predictive of Mediterranean diet adherence (KIDMED) and ultra-processed food (UPF) consumption and also had higher capacity in discriminating children’s diet quality than that of established sociodemographic determinants. Further, we showed that the identified metabolite panel also reflected the associations of these diet quality indicators with C-peptide, a stable and accurate marker of insulin resistance and future risk of metabolic disease. This methodology enables objective assessment of dietary patterns in European child populations, complementary to traditional questionary methods, and can be used in future studies to evaluate diet quality. Moreover, this knowledge can provide mechanistic evidence of common biological pathways that characterize healthy and unhealthy dietary patterns, and diet-related molecular alterations that could associate to metabolic disease.

Keywords
*   Metabolomics
*   NMR spectroscopy
*   Mediterranean Diet Adherence
*   Ultra-processed Food Intake
*   European Children

## Introduction

Dietary habits are considered a key element for the prevention of chronic non-communicable diseases.1 In 2017, 22% of all deaths among adults, were attributed to dietary risks, with type 2 diabetes being among the top causes of diet-related deaths.2 However, it is notoriously difficult to measure diet accurately in large population studies and the absence of accurate dietary assessment methods is hampering the evidence linking diet and disease.3, 4 There is a need for novel approaches to better elucidate diet-related metabolic alterations and their association with disease risk.

Metabolomics is the systematic study of small-molecule metabolites in a biological system and has recently emerged as a powerful top-down approach providing a comprehensive phenotype of biological status. Urinary metabolic phenotypes carry rich information on environmental, lifestyle and nutritional exposures, physiological and metabolic status, and disease risks on an individual and population level.5–7 Urine specimens have high concentrations of food-derived metabolites and studies have shown that urinary metabolic profiles could provide an objective measure of dietary intake.8

Previous research identifying diet-related metabolic profiles, has largely focused on selected food groups including fruits, vegetables, meat and seafood,5, 9–12 while dietary patterns and food processing are far less studied.6, 7, 13, 14 Ultra-processed foods (UPFs), which are industrial formulations undergoing a series of physical and chemical processes and typically lack intact healthy food components and include various additives, can result to cumulative intake of salt, added sugars, and fats.15 UPF consumption has been increasing worldwide16–19 and, to our knowledge, no previous report exists on its metabolic signature. In contrast with the study of overall diets, including the UPF diet, exploring intakes of single foods, which is what is traditionally done in nutrition research, might be failing to provide a realistic image of dietary metabolic footprint. Moreover, most previous studies have focused on adults, and little is known about body’s metabolic response to diet during childhood. This is important as age is a major source of variation in metabolite profiling, with large differences being observed between children and adults.20

Alterations in metabolite profiling could provide a mechanistic link between diet and disease development as early as in childhood. Alterations in the levels of several metabolic biomarkers and pathways have been associated with insulin resistance in childhood including, branched-chain and aromatic amino acid metabolism, urea cycle, glucose and carbohydrate metabolism.21–23 In addition, dietary patterns characterized by high consumption of UPFs, such soft drinks, sweet and savoury snacks, have been related to higher risk of insulin resistance in children.24, 25 All previous studies have included measurement of serum insulin levels as marker of insulin resistance. An alternative marker is C-peptide, a protein that is co-secreted with insulin on an equimolar basis from pancreatic beta cells and has been shown to strongly predict metabolic disease progression.26 C-peptide has a longer half-life than insulin and is recognized as a stable and accurate marker of endogenous insulin secretion, even in non-fasting conditions.27–29 Previous studies in children have shown that higher carbohydrate intake is related to higher C-peptide concentrations.30, 31 However, there is little evidence on the metabolic signatures underlying the association of diet with C-peptide levels in children. We conducted a multi-country study in European children within the Human Early-Life Exposome (HELIX) project,32 aiming *1*) to identify urinary metabolites associated with Mediterranean diet adherence and UPF consumption, and *2)* to determine the extent to which these metabolites were associated with C-peptide, used as an early marker of metabolic health.

## Methods

### Study population

This study is embedded within the Human Early-Life Exposome (HELIX) project,32 a collaborative project across six established and ongoing longitudinal population-based birth cohort studies in Europe: Born in Bradford (BiB, United Kingdom),33 Étude des Déterminants pré et postnatals du développement et de la santé de l’Enfant (EDEN, France),34 Kaunas Cohort (KANC, Lithuania),35 INfancia y Medio Ambiente (INMA, Spain),36 Norwegian Mother, Father and Child Cohort Study (MoBa, Norway),37 and RHEA (RHEA, Greece).38 Participating cohorts covered singleton deliveries from 2003 to 2008. As part of HELIX, a sub-cohort of 1301 children (approximately 200 children in each cohort) were followed in 2014 to 2015 for a clinical examination, a computer-assisted interview with the parents, and the collection of biological samples. Data collection was standardized across cohorts and performed by trained staff. A full description of the HELIX follow-up methods and study population are provided by Maitre et al.32

Our study population consisted of 1147 children with available information on dietary intake, plasma C-peptide levels, and metabolomic biomarkers in urine gathered during the HELIX follow-up at a mean age of 7.9 years (range: 5.4-12.0 years) (Supplementary Figure 1). All participating families provided written informed consent. Approval for the HELIX project was obtained from the local ethical committees at each site. Additionally, the current study was approved by the University of Southern California Institutional Review Board.

### Dietary assessment

Information about the children’s habitual diet was collected via a semi quantitative food-frequency questionnaire (FFQ) covering the child’s habitual diet, which was filled in by the parent attending the examination appointment. The FFQ, covering the past year, was developed by the HELIX research group, translated and applied to all cohorts.32 It included 43 questions of intake of food items, which were aggregated in 16 main food groups: meat & meat products, fish & seafood, sweets, beverages, potatoes, vegetables, dairy products, fruits, bread & cereal, sweet bakery products, added fats, eggs, nuts, salty snacks, pulses and dressings. It also included 15 specific questions to examining the degree of adherence to Mediterranean diet. Diet quality was assessed using two different approaches, i) by assessing the degree of adherence to a Mediterranean diet based on the KIDMED index (Mediterranean Diet Quality Index for children and adolescents)39 and ii) by assessing the proportion of UPF in the overall diet.15

For the KIDMED index,39 items positively associated with the Mediterranean diet pattern (11 items) were assigned a value of +1, while those negatively associated with the Mediterranean diet pattern (4 items) were assigned a value of −1 (Supplementary Table 1). The scores for all 15 items were summed, resulting in a total KIDMED score ranging from −4 to 11, with higher scores reflecting greater adherence to a Mediterranean diet. We categorized the score into three groups: low (<1); moderate (1-4); and high (>4).

For ultra-processed food (UPF) intake, we identified foods and drinks as “ultra-processed” by using the NOVA classification, a food classification system based on the nature, extent, and purpose of industrial food processing.15 We identified the following “ultra-processed” foods: cookies, pastries, sugar-sweetened, low-sugar and artificially-sweetened beverages, cold meat cuts; ham, dairy desserts, sugar-sweetened and other breakfast cereals, crispbread and rusks; chocolate, sweets, margarine, dressings and salty snacks. For some food items, our FFQ did not provide enough information on food processing to determine if a specific item belongs to one processing category or another. We discussed the classification of each food item with a team of nutritionists and used a conservative approach, such that the lower level of processing was chosen-for instance, we made the assumption that fries are homemade from fresh potatoes, and therefore, they were not classified as UPF. For each child, we calculated the daily proportion of all UPF in the total diet as the ratio between the sum of daily servings of UPF to the total daily sum of all food and drink servings. More details on the categorization of foods according to the NOVA classification are presented in Supplementary Table 2.

### Urine metabolite profiling

Two urine samples, representing last night-time and first morning voids, were collected on the evening and morning before the clinical examination, kept in a fridge and transported in a temperature-controlled environment, and aliquoted and frozen within 3 h of arrival at the clinics. They were subsequently pooled to generate a more representative sample of the last 24 h for metabolomic analysis 9.

Urinary metabolic profiles were acquired using 1H NMR spectroscopy according to 9. In brief one-dimensional 600 MHz 1H NMR spectra of urine samples from each cohort were acquired on the same Bruker Avance III spectrometer operating at 14.1 Tesla within a period of 1 month. The spectrometer was equipped with a Bruker SampleJet system, and a 5-mm broad-band inverse configuration probe maintained at 300K. Prior to analysis, cohort samples were randomised. Deuterated 3-(trimethylsilyl)-[2,2,3,3-d4]-propionic acid sodium salt (TSP) was used as internal reference. Aliquots of the study pooled quality control (QC) sample were used to monitor analytical performance throughout the run and were analysed at an interval of every 23 samples (i.e. 4 QC samples per well plate). The 1H NMR spectra were acquired using a standard one-dimensional solvent suppression pulse sequence. 44 metabolites were identified and quantified as described in 9. The urinary NMR showed excellent analytical performance, the mean coefficient of variation across the 44 NMR detected urinary metabolites was 11%. For the statistical analysis we have used creatinine-normalized metabolite concentrations (μmol/mmol of creatinine).

### Plasma C-peptide

Blood was collected at the end of the clinical examination during the HELIX follow-up visit. The median postprandial interval (time between last meal and blood collection) was 3.3 h (IQR: 2.8–4.0).

For each cohort, concentration of C-peptide was assessed in child plasma at the CRG/UPF Proteomics Unit (Barcelona, Spain) using the xMAP and Luminex System multiplex platform according to the manufacturer’s protocol. Blood samples were randomized and blocked by cohort prior to measurement to ensure a representation of each cohort in each plate (batch). For protein quantification, an 8-point calibration curve per plate was performed with protein standards provided in the Luminex kit and following the procedures described in the standard procedures described by the vendor. Commercial heat inactivated, sterile-filtered plasma from human male AB plasma (Sigma Cat #. H3667) was used as constant controls to control for intra-and inter-plate variability. Four control samples were added per plate. Raw intensities obtained with the xMAP and Luminex system for each sample were converted to pg/mL using the calculated standard curves of each plate and accounting for the dilutions that were made prior measurement. The coefficient of variation for C-peptide was 16%. The LOD was determined and the lower and upper quantification limits (LOQ1 and LOQ2, respectively) were obtained from the calibration curves. C-peptide concentrations were log2-transformed to achieve normal distribution. Plate batch effect was corrected by subtracting for each individual and each protein the difference between the overall protein average minus the plate specific protein average. Finally, values below LOQ1 and above LOQ2 were imputed using a truncated normal distribution using the truncdist R package.

### Covariates

Adjustment factors were selected *a priori* based on literature40–42 and included: maternal age (in years), maternal education level (low, middle, high), maternal pre-pregnancy BMI (in kg/m2), family affluence score (cohort-specific definition of low, middle, high), child sex, child age (in years), child BMI (in kg/m2), child sedentary behavior (minutes/day of time spent watching TV, playing computer games or other sedentary games), child ethnicity (White European, Asian, other), and postprandial interval (in hours). We also included a cohort indicator as a fixed effect in the models, as this, in the context of an observational study, is expected to control for cohort effects.43 We imputed missing values for covariates (ranging from 0 to 4%) using the method of chained equations with the R package *mice*. Details about the imputation process in HELIX, diagnostics, and comparison between imputed and complete-case values have been reported in detail elsewhere.44

### Statistical analysis

As a first step in our analysis, we conducted a metabolome-wide association study (MWAS) to assess the associations of urinary metabolites with diet quality. Creatinine-normalized metabolite concentrations (μmol/mmol of creatinine) were log10 transformed prior to statistical analyses to improve model fit. We fitted separate multivariable regression models for each metabolite with the KIDMED score or UPF intake (expressed as per 5% change of total daily food intake). To account for multiple hypothesis testing, we applied the Benjamin–Hochberg false discovery rate (FDR) correction; an FDR-corrected *P* value <0.05 denoted statistical significance. For metabolites identified to be associated with the diet quality indicators, we assessed between-cohort heterogeneity with the I2 statistic and χ2 test from Cochran’s Q.

To examine the ability of the identified metabolite panels in discriminating children with low *vs.* high KIDMED scores (<1 *vs* >4) and low *vs.* high UPF intakes (Quartile 1: <18% *vs* Quartile 4: ≥29%), we plotted receiver operating characteristic (ROC) curves and estimated area under the ROC curve (AUC) values, indicative of the discriminative performance of the metabolite models, based on ten-fold cross validations. We also repeated the ROC analysis for a set of established sociodemographic factors (maternal education level, maternal pre-pregnancy BMI, family affluence score, child sedentary behavior, ethnicity, age, and sex) linked to childhood diet quality both previously40–42 and in our study population (all P<0.05) and compared the discriminative performance of this sociodemographic set with that of the metabolites.

Next, we examined the association of the KIDMED score and of UPF intake (as independent variables) with plasma C-peptide concentration (as dependent variable) using multivariable linear regression models. No departures from linearity in the associations of these diet quality indicators with C-peptide concentration were observed both visually and statistically (P for linearity > 0.42) using generalized additive models. We examined the KIDMED score both as continuous (per score unit increase) and in categories of low (score<1, reference), moderate (score=1-4), and high (score>4). Likewise, UPF intake was assessed both as continuous (per 5% change of total daily food intake) and in quartiles (Q1: <18%, reference; Q2: 18%-<23%; Q3: 23%-<29%; and Q4: ≥29% of total daily food intake). We also included a product term between KIDMED and UPF intake in the regression analysis to assess their interaction; to simplify interpretation of this model, we categorized the KIDMED score as low/moderate *vs* high and UPF intake based on the median population intake (<23% *vs.* ≥23%). We conducted two sets of sensitivity analyses to assess the robustness of the results. First, we calculated cohort-specific effect estimates and assessed heterogeneity with the I2 statistic and χ2 test from Cochran’s Q. Second, we examined potential effect modification by child sex and by child weight status (IOTF-defined normal weight vs overweight/obese) on C-peptide by testing the multiplicative interaction term between the potential effect modifier and each diet quality measure. Finally, we fitted regression models with the metabolites found to be associated with diet quality indicators, and C-peptide in order to assess whether diet-related metabolites are related to β-cell function.

We performed analyses with both complete (missingness <4% in each covariate) and imputed data. Results were similar across raw and imputed data analyses, and hence, we present those using the imputed covariate data. For easier interpretation of effect estimates for the log transformed C-peptide and metabolite values, we back-transformed regression coefficients and present results as percent change (% change =(back-transformed[beta] – 1) × 100).

Analyses were conducted using STATA version 14.2 (StataCorp LLC, TX) and R software version 3.5.3. Linear regression analyses were performed in STATA with the command “*regress” and using “mi estimate”* to account for the imputed covariate data.45 ROC analysis was performed in R with the caret package.46 Visualizations of the results were carried out using the ggplot2 package in R.47

## Results

### Study population

Among participating children (n=1147), 626 (54.6%) were boys and 1028 (89.6%) were white (Table 1). The mean (SD) age at assessment was 7.9 (1.6) years. Median (IQR) C-peptide concentration was 1.26 (0.03, 1.95) ng/mL.

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T1)

Table 1. Characteristics of the study population

Seventeen percent of children (n=195) had a high KIDMED score (>4), indicative of high adherence to the Mediterranean diet. Mean (SD) UPF intake in the overall study population was 24.2 (8.7) % of total daily food intake. KIDMED and UPF intake were negatively correlated (Spearman r = −0.44); children with high and low (<1) KIDMED score had a mean (SD) UPF intake of 18.8 (6.6) % and 33.4 (9.9) %, respectively. The two diet quality scores were associated with most recorded food intakes in opposite directions (Supplementary Tables 3 and 4). Daily fruit and vegetable consumption, weekly fish consumption and not skipping breakfast were the major dietary habits differentiating children with low and high KIDMED score (Supplementary Table 1). Intakes of pastries and bakery products, dairy desserts, margarine and dressings were the major determinants of UPF intake (Supplementary Figure 2). Children with the highest KIDMED scores were mostly from Norway and Spain (Supplementary Figure 3), while those with the highest UPF intake were mostly from Lithuania and the UK (Supplementary Figure 4).

### Diet quality and urinary metabolome

Figure 1 and Supplementary Tables 5 and 6 show the association of diet quality indicators with urinary metabolite levels in childhood. KIDMED and UPF intake exhibited an opposite pattern of association for most of the metabolites (30 out of 43). After controlling for false discovery rate (FDR-corrected P value<0.05), we found that a panel of 4 metabolites related to both diet quality indicators. Specifically, a higher KIDMED score was associated with higher levels of hippurate, *N*-methylnicotinic acid, and urea and with lower levels of sucrose; UPF intake exhibited opposite associations with these 4 metabolites. A higher KIDMED score was also associated with higher acetate and pantothenic acid concentrations, while tyrosine and valine concentrations were inversely associated with UPF intake. There was no evidence of significant between-cohort heterogeneity in the associations between the diet quality indicators and these metabolites (I2 <30%; p for heterogeneity >0.2). ROC curve analyses showed that the combination of 4 metabolites associated with both diet quality indicators performed better than individual metabolites in discriminating children with high and low KIDMED scores and UPF intake (Figure 2). The discriminative ability of this metabolite panel was not improved with the addition of urinary metabolites specifically linked to each diet quality indicator and was equal or even greater to that of established sociodemographic factors previously linked to diet quality in childhood. The regression formulas (scores) for predicting children’s diet quality indicators based on the urinary metabolites are given in Supplementary Table 7.

![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/09/14/2021.09.09.21263335/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/F1)

Figure 1. Adjusted associations of diet quality indicators with urinary metabolites in childhood.
Linear regression models were adjusted for maternal age, maternal education level, maternal pre-pregnancy BMI, family affluence status, child sex, child age, child BMI, child sedentary behavior, child ethnicity, and a cohort indicator. The blue line represents a P-value of 0.05. The red line represents an FDR-adjusted P-value of 0.05. 2-HIB, 2-hydroxyisobutyrate; 3-AIB, 3-aminoisobutyrate; 3-HB/3-AB, 3-hydroxybutyrate/3-aminoisobutyrate; 3-HIB, 3-hydroxyisobutyrate; 3-HIS, 3-hydroxyisovalerate; me-NAM, N1-methyl-nicotinamide; *N*-Acet-NA, *N*-acetyl neuraminic acid; *N*-Me-2-pyr-5-Carb, *N*-methyl-2-pyridone-5-carboxamide; *N*-me-NA, *N*-methylnicotinic acid; *N*-me-PA, *N*-methylpicolinic acid; TMAO, Trimethylamine *N*-oxide.

![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/09/14/2021.09.09.21263335/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/F2)

Figure 2. Receiver operating characteristic (ROC) curves reflecting the ability of urinary metabolites of interest in discriminating adherence to diet quality in childhood.
Panel A illustrates the ability of urinary metabolites of interest in discriminating high adherence to the Mediterranean diet (KIDMED > 4) from low adherence (KIDMED <1). Panel B illustrates the ability of urinary metabolites of interest in discriminating high ultra-processed food consumption (UPF ≥29% of total intake) from low consumption (UPF <18% of total intake). The area under the receiver operating characteristic curve (AUC) values is displayed in the legend. The common panel of 4 metabolites includes the metabolites associated with both diet quality indicators (hippurate, sucrose, urea and *N*-methylnicotinid acid). The panel of 6 metabolites includes the metabolites associated with each diet quality indicator (common panel of 4 plus acetate and pantothenic acid for KIDMED, and plus valine and tyrosine for UPF). The panel of sociodemographic factors includes maternal education level, maternal pre-pregnancy BMI, family affluence score, child sedentary behavior, ethnicity, age, and sex.

### Diet quality and C-peptide levels

Table 2 presents the associations of diet quality indicators with C-peptide concentration in childhood. We found that a higher KIDMED score was associated with lower C-peptide. Specifically, compared to a low KIDMED score (<1), children with a moderate score (1-4) had a 28% lower C-peptide concentration (percent change: −27.7, 95% CI: −49.6 to 3.9) and those with a high score (>4) had a 39% lower C-peptide concentration (percent change: −39.0, 95% CI: −60.6 to −5.7) (P-trend =0.03). An opposite association was observed for UPF intake. Compared to children at the lowest quartile of UPF intake (<18% of total daily food intake), children at the second quartile (18-<23% of total daily food intake) had a 24% higher C-peptide concentration (percent change: 24.3, 95% CI: −6.4 to 65.2), those at the third quartile (23-<29% of total daily food intake) had a 39% higher concentration (percent change: 38.5, 95% CI: 3.8 to 84.9), and those at the fourth quartile (≥29% of total daily food intake) had a 46% higher concentration (percent change: 46.0, 95% CI: 8.1 to 97.3) (P-trend=0.01). There was no evidence of interaction between the diet quality indicators (Supplementary Table 8). When we examined the associations separately in each cohort, we found no significant between-cohort heterogeneity (I2<18%, p for heterogeneity >0.29) (Supplementary Figures 5 and 6). We also found no evidence that the observed associations differed by the sex of the children or their weight status (Supplementary Table 9).

View this table:
[Table 2.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T2)

Table 2. Adjusted associations of diet quality with C-peptide levels in childhood1

### Diet-related urinary metabolites and C-peptide levels

When we examined the associations with the metabolite scores for each diet quality indicator, we found that the scores for KIDMED were associated with lower C-peptide, while opposite associations were observed with the metabolite scores for UPF (Supplementary Table 10). In the pairwise associations between the individual urinary metabolites linked to KIDMED or UPF intake and C-peptide, we found that higher levels of sucrose were associated with higher C-peptide levels

## Discussion

This is the first study of European children that identified a panel of urinary metabolites associated with diet quality, as assessed by a Mediterranean diet adherence score (KIDMED) and ultra-processed food (UPF) intake. These metabolites include food constituents (or their metabolic products), vitamin related compounds and metabolites related to amino acid, protein and carbohydrate metabolism. For both diet quality indices, KIDMED score and UPF intake, there was a common panel of 4 metabolites exhibiting an opposite pattern of association: higher levels of hippurate, *N*-methylnicotinic acid and urea, and lower levels of sucrose were concomitant with a higher KIDMED score and lower UPF intake. Moreover, these urinary metabolites that reflected these diet quality indicators, were also associated with C-peptide levels, a well-known marker of-cell function and insulin resistance.26

Previous studies assessing the metabolic response of the body to diet have largely focused on specific foods groups,6 yet humans eat a combination of foods that may have additive or interactive effects on human physiology. In our study, we assessed two diet quality indices to better describe the complexity of overall diet. The Mediterranean diet is characterized by high intake of fruits, vegetables, legumes, nuts, and whole grain products, fish and low intakes of red meat and sweets, and it has long been appraised for its cardiometabolic benefits, even in children.48 Moreover, we assessed UPF intake to characterize the intake of industrial processed and prepared food items, with altered food structure, nutritional content, and taste49. UPF intake assessment diverts from the traditional strategy of studying nutrients, foods or dietary patterns to identify the link between diet with health and disease, and it reflects the cumulative intake of artificial substances (i.e. flavorings, colorings, emulsifiers, and other additives) and processing by-products. Effects of UPF consumption are likely to be the result of synergistic effects of many food ingredient compound and characteristics, and high consumption of UPF has been proposed to lead to metabolic dysregulation.50

For both KIDMED and UPF, there was a common panel of 4 urinary metabolites (hippurate, *N*-methylnicotinic acid, urea, sucrose), exhibiting an opposite pattern of association. In addition, KIDMED was positively associated with pantothenic acid and acetate, whereas UPF was negatively associated with two amino acids, valine and tyrosine. These metabolic signatures demonstrate a core common metabolic biomarker panel reflecting habitual dietary intake, but also show that the two diet quality indices and their metabolic signatures act complementary and can highlight different aspects of human metabolism and physiology. Our finding are consistent with previous metabolomics studies conducted in adults.13, 51, 52 Four dietary interventions were developped13 with similar energy content and within the WHO healthy eating guidelines, but with varying macro-and micronutrient intake. The urine metabolic profile for the diet most concordant with the guidelines, characterized by high intakes of dietary fiber (through fruits, vegetables and whole grain cereal products) and low intakes of fat, sugar and salt showed systematic differences for a total of 28 metabolites, including increased levels of acetate, hippurate, *N*-methylnicotinic acid and urea, compared to the diet that diverted from the healthy guidelines.13 The PREDIMED study also used NMR to define urinary biomarkers associated with a high adherence to a Mediterranean diet pattern in adults and the proposed biomarkers also included higher levels of urea.52 To our knowledge, our study is the first to show the urinary metabolic footprint of total ultra-processed food intake, and also the first to show that urine NMR-derived scores of diet quality are reflective of a key biomarker of metabolic health, c-peptide, in healthy children. The ability to use a rapid, non-invasive biofluid screen to measure objective biomarkers of diet quality in children opens new avenues for exploring the significance of nutritional patterns to healthy development early in life.

Among the diet-related urinary metabolites, sucrose individually reflected the associations observed between diet quality and C-peptide.53 Elevated levels of added sugars could lead to elevated glucose load in the human body, which in turn can lead to an increase in glycemic response and, thus, C-peptide production in healthy populations. To the best of our knowledge, there is only one previous study from Mexico examining the relation of diet to C-peptide in children.54 Similar to our findings, this study showed that adherence to a prudent dietary pattern characterized by food groups commonly consumed in Mediterranean diet (vegetables, fruit, fish, and legumes) was associated with low C-peptide levels in boys. Regarding other markers of glucose regulation, our findings are in line with previous studies in children and adults reporting that low adherence to a healthy dietary pattern or high consumption of specific ultra-processed foods (sugar-sweetened beverages, and ultra-processed meat) were associated with impaired glucose regulation and insulin resistance.55–60 Overall, our findings are consistent with the public health concerns raised by the World Health Organization (WHO), relating poor overall quality with high intake level of sugars and added sugars, and with poor metabolic health and a high risk of several non-communicable diseases.53, 61–63

Regarding the other urinary metabolites reflecting diet quality, hippurate is a normal component of urine, a metabolic product of phenolic compounds which are present in various dietary sources. It is also a biomarker of fruit/vegetable intake, as it has been confirmed in previous studies in healthy children, adolescents, 9, 64 and adults.65 *N*-methylnicotinic acid (trigonelline) is a product of the metabolism of niacin (vitamin B3) and a biomarker of various dietary sources like legumes66 and fruits.9 Urea is the principal end product of amino acids and protein catabolism and a potential marker of protein intake from food.51

In addition, adherence to the Mediterranean diet was also positively associated with urinary levels of pantothenic acid and acetate. Pantothenic acid (vitamin B5, necessary to form coenzyme-A) is present in many foods, and we have previously reported a positive association between consumption of dairy products and urinary pantothenic acid, in the same study population.9 Acetate has also multiple dietary sources, it is produced by acetate producing bacteria in foodstuff and urinary acetate is also modulated by human gut microbial metabolism. We have previously reported a positive association between potato consumption and urinary acetate levels, in the same study population.9 Both compounds have a central role in human biochemistry and the metabolism and synthesis of carbohydrates, proteins, and fats.

Moreover, UPF intake was negatively associated with two amino acids, valine and tyrosine. Tyrosine is regarded as a conditionally essential amino acid in adults and essential in children. Foods high in dietary tyrosine include dairy, meat, eggs, beans, nuts, grains. Tyrosine is a precursor for neurotransmitters and hormones, increases dopamine availability which in turn could enhance cognitive performance.67 Valine is an essential branch chain amino acid critical to energy homeostasis, protein and muscle metabolism.68, 69 Neither of these metabolites have been linked with habitual dietary intake before.

### Strengths and limitations

The main strengths of the study are the multicentric design which included children from 6 countries spanning north to south in Europe, the use of standardized data collection and biomarker measurement protocols across cohorts, and the fairly large sample size with biomarker data. Identified panels of urinary metabolites had similar or higher capacity in discriminating children’s diet quality to that of established sociodemographic determinants. We chose 1H NMR spectroscopy for urinary analysis as this is an inherently high reproducible, high throughput technique suitable for the identification and quantification of urine metabolites, which are typically of high concentrations, without complex sample preparation which could potentially introduce analytical biases. Also, urinary 1H NMR spectroscopy has been applied in many other cohort studies and proposed as an objective method for dietary assessment13 potentially facilitating comparative studies in the future. We have used a pooled urine sample collection design which combined the last sample before bedtime with the first morning void sample of the following day, and we have shown in our preliminary work that this sample collection strategy has the advantage of reducing diurnal variations.9, 70

Our study has also some limitations. As in any observational study, there is the possibility of unmeasured residual confounding. In our analyses, we took into account a number of sociodemographic and lifestyle factors in childhood (e.g., socioeconomic status, ethnicity, sedentary behavior) that are associated with both diet quality and glycemic response. We did not have data available to control for energy intake. However, in all our models, we included BMI of the children, a measure strongly correlated to energy intake,71 and assessed ultra-processed food intake as proportion of total food intake. Moreover, the absence of heterogeneity across cohorts with different correlation structures and confounding patterns in their data32, provide evidence to support that unmeasured confounding is unlikely to have influenced the observed associations. Since the data collected are cross-sectional, and there is no temporality in the observed associations, further longitudinal studies examining metabolic and glycemic alterations in relation to diet quality are needed. Although 1H NMR spectroscopy had the advantage of improving the specificity of the quantitation and provided explicit metabolite identification, it limited the number of metabolites being measured and provided partial coverage of the urine metabolome. Absolute quantification of some metabolites with exchangable protons such as urea could also be negatively impacted by the solvent suppression methods required for 1H NMR spectroscopy of urine. Supplementing the current study with other complementary untargeted and targeted metabolomic approaches in future, such as mass spectrometry, would help enhance identification and robust quantification of urinary metabolites associated with diet quality in children.

In summary, this multi-center European study showed that urinary metabolic profiles related to food constituents (or their metabolic products), to amino acid and carbohydrate metabolism reflect adherence to the Mediterranean diet and ultra-processed food intake in childhood. Higher adherence to Mediterranean diet, lower ultra-processed food intake, and lower levels of the diet-related carbohydrate sucrose were associated with lower C-peptide levels, a marker of β-cell function. These results provide evidence to support efforts by public health authorities to recommend increased adherence to the Mediterranean diet and limiting UPF consumption in childhood. Further prospective studies examining the association of diet quality and related metabolomic profiles with C-peptide and other surrogates of insulin resistance are needed to replicate our findings.

## Data Availability

The HELIX data warehouse has been established as an accessible resource for collaborative research involving researchers external to the project. Access to HELIX data is based on approval by the HELIX Project Executive Committee and by the individual cohorts. Further details on the content of the data warehouse (data catalog) and procedures for external access are described on the project website (http://www.projecthelix.eu/index.php/es/data-inventory).

## Funding

The research leading to these results has received funding from the European Community’s Seventh Framework Programme (FP7/2007–2013) under grant agreement no 308333 – the HELIX project. The STOP project ([http://www.stopchildobesity.eu/](http://www.stopchildobesity.eu/)) received funding from the European Union’s Horizon 2020 research and innovation programme under Grant Agreement No. 774548. The STOP Consortium is coordinated by Imperial College London and includes 24 organisations across Europe, the United States and New Zealand. The content of this publication reflects only the views of the authors, and the European Commission is not liable for any use that may be made of the information it contains.

INMA data collections were supported by grants from the Instituto de Salud Carlos III, CIBERESP, and the Generalitat de Catalunya-CIRIT. KANC was funded by the grant of the Lithuanian Agency for Science Innovation and Technology (6-04-2014_31V-66). For a full list of funding that supported the EDEN cohort, refer to: Heude B et al. Cohort Profile: The EDEN mother-child cohort on the prenatal and early postnatal determinants of child health and development. Int J Epidemiol. 2016 Apr;45(2):353-63. The Norwegian Mother, Father and Child Cohort Study (MoBa) is supported by the Norwegian Ministry of Health and Care Services and the Ministry of Education and Research. The Rhea project was financially supported by European projects, and the Greek Ministry of Health (Program of Prevention of Obesity and Neurodevelopmental Disorders in Preschool Children, in Heraklion district, Crete, Greece: 2011–2014; ‘Rhea Plus’: Primary Prevention Program of Environmental Risk Factors for Reproductive Health, and Child Health: 2012–2015). Born in Bradford received funding from the Wellcome Trust (101597). Professor Wright and McEachan receive funding from the National Institute for Health Research Applied Research Collaboration for Yorkshire and Humber. The views expressed are those of the author(s) and not necessarily those of the NIHR or the Department of Health and Social Care. Dr. Maribel Casas received funding from Instituto de Salud Carlos III (Ministry of Economy and Competitiveness) (MS16/00128). Dr. Leda Chatzi was supported by NIH/NIHES R01 ES029944, R01ES030691, R01ES030364, R21 ES029681, R21 ES028903, and P30 ES007048-23. Dr. David Conti was supported by P01CA196569, R01CA140561, R01 ES016813, R01 ES029944, R01ES030691, R01ES030364. Dr. Nikos Stratakis was supported by NIH/NIHES R21 ES029681 and P30 ES007048-23, and NIH/NIDDK P30 DK048522-24. Dr. Hector Keun & Dr. Alexandros Siskos were also supported by the European Union’s Horizon 2020 research and innovation programme under grant agreement No 874583 (“ATHLETE”). Dr. Eleni Papadopoulou was supported by the Research Council of Norway, under the MILJØFORSK program (project number: 268465). Dr. Fernanda Rauber was supported by the Fundação de Amparo à Pesquisa do Estado de São Paulo 2016/14302-7 and 2018/19820-1. The CRG/UPF Proteomics Unit is part of the Spanish Infrastructure for Omics Technologies (ICTS OmicsTech) and it is a member of the ProteoRed PRB3 consortium which is supported by grant PT17/0019 of the PE I+D+i 2013-2016 from the Instituto de Salud Carlos III (ISCIII) and ERDF. ISGlobal acknowledges support from the Spanish Ministry of Science, Innovation and Universities, “Centro de Excelencia Severo Ochoa 2013-2017”, SEV-2012-0208, and “Secretaria d’Universitats i Recerca del Departament d’Economia i Coneixement de la Generalitat de Catalunya” (2017SGR595).

## Author Contributions

Nikos Stratakis (USC), Alexandros P. Siskos (Imperial), Eleni Papadopoulou (NIPH), Hector C. Keun (Imperial) and Leda Chatzi (USC) designed analytical and statistical methods. Alexandros P. Siskos, Chung-Ho E. Lau (Imperial), and Hector C. Keun generated and reviewed the primary metabolomics data. Nikos Stratakis, Anh N. Nguyen (Erasmus), and Katerina Margetaki (USC) analyzed the data. Nikos Stratakis, Alexandros P. Siskos, Eleni Papadopoulou, Hector C. Keun and Leda Chatzi wrote the first draft of the paper. Martine Vrijheid (ISGlobal) coordinated the HELIX project. Martine Vrijheid, Lea Maitre, and Oliver Robinson (ISGlobal) coordinated the HELIX data collection. All authors contributed to the data collection and/or interpretation of the results. All authors contributed to and approved the manuscript.

## Supplementary material

View this table:
[Supplementary Table 1.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T3)

Supplementary Table 1. KIDMED items, and their scoring, used to assess adherence to the Mediterranean diet in HELIX children

View this table:
[Supplementary Table 2.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T4)

Supplementary Table 2. Items included in the FFQ by ultra-processed food (UPF) inclusion and information regarding the extent and purpose of food processing

View this table:
[Supplementary Table 3.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T5)

Supplementary Table 3. Intakes of food groups (in servings/week) by categories of the KIDMED score1

View this table:
[Supplementary Table 4.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T6)

Supplementary Table 4. Intakes of food groups (in servings/week) by quartiles of ultra-processed food intake1

View this table:
[Supplementary Table 5.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T7)

Supplementary Table 5. Associations of KIDMED score with urinary metabolites in childhood1

View this table:
[Supplementary Table 6.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T8)

Supplementary Table 6. Associations of ultra-processed food consumption with urinary metabolites in childhood1

View this table:
[Supplementary Table 7.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T9)

Supplementary Table 7. Regression formulas (scores) for predicting diet quality in childhood based on panels of urinary metabolites (creatinine-normalized concentrations, μmol/mmol of creatinine)

View this table:
[Supplementary Table 8.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T10)

Supplementary Table 8. Interaction of diet quality indicators in association to Cpeptide concentration in childhood1

View this table:
[Supplementary Table 9.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T11)

Supplementary Table 9. Associations of diet quality with C-peptide concentration in childhood after stratifying by sex and by weight status, respectively1

View this table:
[Supplementary Table 10.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/T12)

Supplementary Table 10. Associations between urinary metabolites (μmol/mmol of creatinine) linked to diet quality and C-peptide in childhood1

![Supplementary Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/09/14/2021.09.09.21263335/F3.medium.gif)

[Supplementary Figure 1.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/F3)

Supplementary Figure 1. Participant flowchart

![Supplementary Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/09/14/2021.09.09.21263335/F4.medium.gif)

[Supplementary Figure 2.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/F4)

Supplementary Figure 2. Median intakes (servings/day) of the foods and drinks identified as ultra-processed (UPF) by quartile of UPF intake.
Quartiles of ultra-processed food intake (expressed as % of total daily food intake) were as follows: Q1, ≤18%; Q2, >18% but <23%; Q3, ≥23% but <29%; and Q4, ≥29%.

![Supplementary Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/09/14/2021.09.09.21263335/F5.medium.gif)

[Supplementary Figure 3.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/F5)

Supplementary Figure 3. Levels of childhood adherence to the Mediterranean in each HELIX cohort.
Adherence levels of Mediterranean diet were defined as follows: low, KIDMED score, <1; moderate, KIDMED score, 1-4; and high, >4.

![Supplementary Figure 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/09/14/2021.09.09.21263335/F6.medium.gif)

[Supplementary Figure 4.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/F6)

Supplementary Figure 4. Levels of ultra-processed food consumption in each HELIX cohort.
Levels of ultra-processed food intake (expressed as % of total daily food intake) are based on quartile cutoffs according to the overall intake distribution.

![Supplemental Figure 5.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/09/14/2021.09.09.21263335/F7.medium.gif)

[Supplemental Figure 5.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/F7)

Supplemental Figure 5. Cohort-specific associations of the KIDMED score with C-peptide in childhood.
Beta coefficients (95% CIs) by cohort were obtained using linear regression models adjusted for maternal age, maternal education level, maternal pre-pregnancy BMI, family affluence status, child sex, child age, child BMI, child sedentary behavior, child ethnicity, and postprandial interval. Combined estimates were obtained by using a fixed-effects meta-analysis. Effect estimates were exponentiated and reported as percent change in C-peptide levels. Squares represent the cohort-specific effect estimates; diamond represents the combined estimate; and horizontal lines denote 95%CIs.

![Supplemental Figure 6.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/09/14/2021.09.09.21263335/F8.medium.gif)

[Supplemental Figure 6.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263335/F8)

Supplemental Figure 6. Cohort-specific associations of ultra-processed food intake with Cpeptide in childhood.
Beta coefficients (95% CIs) by cohort were obtained using linear regression models adjusted for maternal age, maternal education level, maternal pre-pregnancy BMI, family affluence status, child sex, child age, child BMI, child sedentary behavior, child ethnicity, and postprandial interval. Combined estimates were obtained by using a fixed-effects meta-analysis. Effect estimates were exponentiated and reported as percent change in C-peptide levels. Squares represent the cohort-specific effect estimates; diamond represents the combined estimate; and horizontal lines denote 95%CIs.

## Acknowledgments

The authors declare no conflict of interest. We acknowledge the input of the entire HELIX consortium. We are grateful to all the participating families in the five cohorts (BiB, EDEN, INMA, MoBa, KANC and RHEA cohorts), that took part in this study. We are equally grateful to all the fieldworkers for their dedication and efficiency in this study. A full roster of the INMA and RHEA study investigators can be found at [http://www.proyectoinma.org/en/inma-project/inma-project-researchers/](http://www.proyectoinma.org/en/inma-project/inma-project-researchers/) and [http://www.rhea.gr/en/about-rhea/the-rheateam/](http://www.rhea.gr/en/about-rhea/the-rheateam/), respectively. The Born in Bradford study is only possible because of the enthusiasm and commitment of the participating children and parents. We are grateful to all the participants, health professionals and researchers who have made Born in Bradford happen. We are also grateful to all the participating families in Norway who take part in the on-going MoBa cohort study. We thank all the children and families participating in the EDEN-HELIX mother-child cohort. We are grateful to Joane Quentin, Lise Giorgis-Allemand, and Rémy Slama (EDEN study group) for their work on the HELIX project. We thank Sonia Brishoual, Angelique Serre and Michele Grosdenier (Poitiers Biobank, CRB BB-0033-00068, Poitiers, France) for biological sample management and Prof Frederic Millot (principal investigator), Elodie Migault, Manuela Boue and Sandy Bertin (Clinical Investigation Center, Inserm CIC1402, CHU de Poitiers, Poitiers, France) for planification and investigational actions. We are also grateful to Veronique Ferrand-Rigalleau, Celine Leger and Noella Gorry (CHU de Poitiers, Poitiers, France) for administrative assistance. We also acknowledge the commitment of the members of the EDEN Mother-Child Cohort Study Group: I. Annesi-Maesano, J.Y. Bernard, J. Botton, M.A. Charles, P. Dargent-Molina, B. de Lauzon-Guillain, P. Ducimetière, M. de Agostini, B. Foliguet, A. Forhan, X. Fritel, A. Germa, V. Goua, R. Hankard, M. Kaminski, B. Larroque, N. Lelong, J. Lepeule, G. Magnin, L. Marchand, C. Nabet, F Pierre, M.J. Saurel-Cubizolles, M. Schweitzer, O. Thiebaugeorges.

## Footnotes

*   * Co-first authors

*   † Co-senior authors

*   Received September 9, 2021.
*   Revision received September 9, 2021.
*   Accepted September 14, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.Grosso G, Mateo A, Rangelov N, Buzeti T, Birt C. Nutrition in the context of the Sustainable Development Goals. Eur J Public Health. 2020;30(Supplement_1):i19- i23.
    
    

2.  2.Collaborators GBDD. Health effects of dietary risks in 195 countries, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2019;393(10184):1958–1972.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

3.  3.Posma JM, Garcia-Perez I, Frost G, et al. Nutriome-metabolome relationships provide insights into dietary intake and metabolism. Nat Food. 2020;1(7):426–436.
    
    

4.  4.Ioannidis JPA. The Challenge of Reforming Nutritional Epidemiologic Research. JAMA. 2018;320(10):969–970.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2018.11025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30422271&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

5.  5.Gibson R, Lau CE, Loo RL, et al. The association of fish consumption and its urinary metabolites with cardiovascular risk factors: the International Study of Macro-/Micronutrients and Blood Pressure (INTERMAP). Am J Clin Nutr. 2020;111(2):280–290.
    
    

6.  6.Collins C, McNamara AE, Brennan L. Role of metabolomics in identification of biomarkers related to food intake. Proc Nutr Soc. 2019;78(2):189–196.
    
    

7.  7.Rebholz CM, Lichtenstein AH, Zheng Z, Appel LJ, Coresh J. Serum untargeted metabolomic profile of the Dietary Approaches to Stop Hypertension (DASH) dietary pattern. Am J Clin Nutr. 2018;108(2):243–255.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ajcn/nqy099&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29917038&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

8.  8.O’Gorman A, Brennan L. The role of metabolomics in determination of new dietary biomarkers. Proc Nutr Soc. 2017;76(3):295–302.
    
    

9.  9.Lau CE, Siskos AP, Maitre L, et al. Determinants of the urinary and serum metabolome in children from six European populations. BMC Med. 2018;16(1):202.
    
    

10. 10.Guertin KA, Moore SC, Sampson JN, et al. Metabolomics in nutritional epidemiology: identifying metabolites associated with diet and quantifying their potential to uncover diet-disease relations in populations. Am J Clin Nutr. 2014;100(1):208–217.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo5OiIxMDAvMS8yMDgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzM1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

11. 11.Playdon MC, Sampson JN, Cross AJ, et al. Comparing metabolite profiles of habitual diet in serum and urine. Am J Clin Nutr. 2016;104(3):776–789.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo5OiIxMDQvMy83NzYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzM1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

12. 12.Scalbert A, Brennan L, Manach C, et al. The food metabolome: a window over dietary exposure. Am J Clin Nutr. 2014;99(6):1286–1308.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo5OiI5OS82LzEyODYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzM1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

13. 13.Garcia-Perez I, Posma JM, Gibson R, et al. Objective assessment of dietary patterns by use of metabolic phenotyping: a randomised, controlled, crossover trial. Lancet Diabetes Endocrinol. 2017;5(3):184–195.
    
    

14. 14.Martinez JA, Keun HC, Siskos AP. Effects of Mediterranean Diet on the Metabolome. In: Romagnolo DF, Selmin OI, eds. Mediterranean Diet: Dietary Guidelines and Impact on Health and Disease. Cham: Springer International Publishing; 2016:121–137.
    
    

15. 15.Monteiro CA, Cannon G, Levy RB, et al. Ultra-processed foods: what they are and how to identify them. Public Health Nutr. 2019;22(5):936–941.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S1368980018003762&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30744710&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

16. 16.Monteiro CA, Moubarac JC, Cannon G, Ng SW, Popkin B. Ultra-processed products are becoming dominant in the global food system. Obes Rev. 2013;14 Suppl 2:21–28.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/obr.12107&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24102801&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

17. 17.Monteiro CA, Cannon G, Moubarac JC, Levy RB, Louzada MLC, Jaime PC. The UN Decade of Nutrition, the NOVA food classification and the trouble with ultra-processing. Public Health Nutrition. 2018;21(1):5–17.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=doi:10.1017/S1368980017000234&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28322183&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

18. 18.Monteiro CA, Moubarac JC, Levy RB, Canella DS, Louzada M, Cannon G. Household availability of ultra-processed foods and obesity in nineteen European countries. Public Health Nutr. 2018;21(1):18–26.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S1368980017001379&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28714422&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

19. 19.Baker P, Machado P, Santos T, et al. Ultra-processed foods and the nutrition transition: Global, regional and national trends, food systems transformations and political economy drivers. Obes Rev. 2020;21(12):e13126.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

20. 20.Ellul S, Wake M, Clifford SA, et al. Metabolomics: population epidemiology and concordance in Australian children aged 11-12 years and their parents. BMJ Open. 2019;9(Suppl 3):106–117.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMzoiOS9TdXBwbF8zLzEwNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA5LzE0LzIwMjEuMDkuMDkuMjEyNjMzMzUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

21. 21.Rauschert S, Kirchberg FF, Marchioro L, Koletzko B, Hellmuth C, Uhl O. Early Programming of Obesity Throughout the Life Course: A Metabolomics Perspective. Ann Nutr Metab. 2017;70(3):201–209.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000459635&link_type=DOI) 

22. 22.Zhao X, Gang X, Liu Y, Sun C, Han Q, Wang G. Using Metabolomic Profiles as Biomarkers for Insulin Resistance in Childhood Obesity: A Systematic Review. J Diabetes Res. 2016;2016:8160545.
    
    

23. 23.Martos-Moreno GA, Mastrangelo A, Barrios V, et al. Metabolomics allows the discrimination of the pathophysiological relevance of hyperinsulinism in obese prepubertal children. Int J Obes (Lond*).* 2017;41(10):1473–1480.
    
    

24. 24.Karatzi K, Moschonis G, Barouti AA, et al. Dietary patterns and breakfast consumption in relation to insulin resistance in children. The Healthy Growth Study. Public Health Nutr. 2014;17(12):2790–2797.
    
    

25. 25.Romero-Polvo A, Denova-Gutierrez E, Rivera-Paredez B, et al. Association between dietary patterns and insulin resistance in Mexican children and adolescents. Ann Nutr Metab. 2012;61(2):142–150.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000341493&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23037180&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

26. 26.Patel N, Taveira TH, Choudhary G, Whitlatch H, Wu WC. Fasting serum C-peptide levels predict cardiovascular and overall death in nondiabetic adults. J Am Heart Assoc. 2012;1(6):e003152.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiYWhhb2EiO3M6NToicmVzaWQiO3M6MTE6IjEvNi9lMDAzMTUyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDkvMTQvMjAyMS4wOS4wOS4yMTI2MzMzNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

27. 27.Hope SV, Knight BA, Shields BM, Hattersley AT, McDonald TJ, Jones AG. Random non-fasting C-peptide: bringing robust assessment of endogenous insulin secretion to the clinic. Diabet Med. 2016;33(11):1554–1558.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27100275&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

28. 28.Vezzosi D, Bennet A, Fauvel J,  Caron P. Insulin, C-peptide and proinsulin for the biochemical diagnosis of hypoglycaemia related to endogenous hyperinsulinism. Eur J Endocrinol. 2007;157(1):75–83.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiZWplIjtzOjU6InJlc2lkIjtzOjg6IjE1Ny8xLzc1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDkvMTQvMjAyMS4wOS4wOS4yMTI2MzMzNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

29. 29.Polonsky KS, Licinio-Paixao J, Given BD, et al. Use of biosynthetic human C-peptide in the measurement of insulin secretion rates in normal volunteers and type I diabetic patients. J Clin Invest. 1986;77(1):98–105.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI112308&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3511094&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1986AYB1800016&link_type=ISI) 

30. 30.Sunehag AL, Toffolo G, Treuth MS, et al. Effects of dietary macronutrient content on glucose metabolism in children. J Clin Endocrinol Metab. 2002;87(11):5168–5178.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2002-020674&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12414888&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000179121800050&link_type=ISI) 

31. 31.Buyken AE, Kellerhoff Y, Hahn S, Kroke A, Remer T. Urinary C-peptide excretion in free-living healthy children is related to dietary carbohydrate intake but not to the dietary glycemic index. J Nutr. 2006;136(7):1828–1833.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjEwOiIxMzYvNy8xODI4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDkvMTQvMjAyMS4wOS4wOS4yMTI2MzMzNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

32. 32.Maitre L, de Bont J, Casas M, et al. Human Early Life Exposome (HELIX) study: a European population-based exposome cohort. BMJ Open. 2018;8(9):e021311.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMToiOC85L2UwMjEzMTEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzM1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

33. 33.Wright J, Small N, Raynor P, et al. Cohort Profile: the Born in Bradford multi-ethnic family cohort study. Int J Epidemiol. 2013;42(4):978–991.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dys112&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23064411&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325167800010&link_type=ISI) 

34. 34.Heude B, Forhan A, Slama R, et al. Cohort Profile: The EDEN mother-child cohort on the prenatal and early postnatal determinants of child health and development. International Journal of Epidemiology. 2016;45(2):353–363.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyv151&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26283636&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

35. 35.Grazuleviciene R, Danileviciute A, Nadisauskiene R, Vencloviene J. Maternal smoking, GSTM1 and GSTT1 polymorphism and susceptibility to adverse pregnancy outcomes. International journal of environmental research and public health. 2009;6(3):1282–1297.
    
    

36. 36.Guxens M, Ballester F, Espada M, et al. Cohort Profile: The INMA--INfancia y Medio Ambiente--(Environment and Childhood) Project. IntJEpidemiol. 2012;41(4):930–940.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyr054&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21471022&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308232200009&link_type=ISI) 

37. 37.Magnus P, Birke C, Vejrup K, et al. Cohort Profile Update: The Norwegian Mother and Child Cohort Study (MoBa). International Journal of Epidemiology. 2016.
    
    

38. 38.Chatzi L, Leventakou V, Vafeiadi M, et al. Cohort Profile: The Mother-Child Cohort in Crete, Greece (Rhea Study). Int J Epidemiol. 2017;46(5):1392–1393k.
    
    

39. 39.Serra-Majem L, Ribas L, Ngo J, et al. Food, youth and the Mediterranean diet in Spain. Development of KIDMED, Mediterranean Diet Quality Index in children and adolescents. Public Health Nutr. 2004;7(7):931–935.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1079/PHN2004556&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15482620&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224447300017&link_type=ISI) 

40. 40.Aranceta J, Perez-Rodrigo C, Ribas L, Serra-Majem L. Sociodemographic and lifestyle determinants of food patterns in Spanish children and adolescents: the enKid study. Eur J Clin Nutr. 2003;57 Suppl 1:S40–44.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.ejcn.1601813&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12947451&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185372100010&link_type=ISI) 

41. 41.Scaglioni S, De Cosmi V, Ciappolino V, Parazzini F, Brambilla P, Agostoni C. Factors Influencing Children’s Eating Behaviours. Nutrients. 2018;10(6).
    
    

42. 42.Patrick H, Nicklas TA. A review of family and social determinants of children’s eating patterns and diet quality. J Am Coll Nutr. 2005;24(2):83–92.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/07315724.2005.10719448&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15798074&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228529000002&link_type=ISI) 

43. 43.Basagana X, Pedersen M, Barrera-Gomez J, et al. Analysis of multicentre epidemiological studies: contrasting fixed or random effects modelling and meta-analysis. Int J Epidemiol. 2018;47(4):1343–1354.
    
    

44. 44.Agier L, Basagana X, Maitre L, et al. Early-life exposome and lung function in children in Europe: an analysis of data from the longitudinal, population-based HELIX cohort. Lancet Planet Health. 2019;3(2):e81–e92.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30737192&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

45. 45.StataCorp. Stata multiple-imputation reference manual. College Station, TX: StataCorp LLC; 2021.
    
    

46. 46.Kuhn M. Building Predictive Models in R Using the caret Package. Journal of Statistical Software. 2008;28(5).
    
    

47. 47.Wickham H. ggplot2: Elegant Graphics for Data Analysis, Springer-Verlag New York. Available at: [https://ggplot2.tidyverse.org](https://ggplot2.tidyverse.org). 2016.
    
    

48. 48.Chatzi L, Rifas-Shiman SL, Georgiou V, et al. Adherence to the Mediterranean diet during pregnancy and offspring adiposity and cardiometabolic traits in childhood. Pediatr Obes. 2017;12 Suppl 1:47–56.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28160450&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

49. 49.Elizabeth L, Machado P, Zinocker M, Baker P, Lawrence M. Ultra-Processed Foods and Health Outcomes: A Narrative Review. Nutrients. 2020;12(7).
    
    

50. 50.Srour B, Fezeu LK, Kesse-Guyot E, et al. Ultraprocessed Food Consumption and Risk of Type 2 Diabetes Among Participants of the NutriNet-Sante Prospective Cohort. JAMA Intern Med. 2020;180(2):283–291.
    
    

51. 51.Garcia-Perez I, Posma JM, Chambers ES, et al. Dietary metabotype modelling predicts individual responses to dietary interventions. Nature Food. 2020;1(6):355- 364.
    
    

52. 52.Almanza-Aguilera E, Urpi-Sarda M, Llorach R, et al. Microbial metabolites are associated with a high adherence to a Mediterranean dietary pattern using a (1)H-NMR-based untargeted metabolomics approach. J Nutr Biochem. 2017;48:36–43.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jnutbio.2017.06.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28692847&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

53. 53.Guideline: Sugars Intake for Adults and Children. Geneva: © World Health Organization, 2015.; 2015.
    
    

54. 54.Perng W, Fernandez C, Peterson KE, et al. Dietary Patterns Exhibit Sex-Specific Associations with Adiposity and Metabolic Risk in a Cross-Sectional Study in Urban Mexican Adolescents. J Nutr. 2017;147(10):1977–1985.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjExOiIxNDcvMTAvMTk3NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA5LzE0LzIwMjEuMDkuMDkuMjEyNjMzMzUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

55. 55.Manios Y, Kourlaba G, Grammatikaki E, et al. Development of a lifestyle-diet quality index for primary schoolchildren and its relation to insulin resistance: the Healthy Lifestyle-Diet Index. European Journal of Clinical Nutrition. 2010;64(12):1399–1406.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ejcn.2010.172&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20808332&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

56. 56.Asghari G, Yuzbashian E, Mirmiran P, Hooshmand F, Najafi R, Azizi F. Dietary Approaches to Stop Hypertension (DASH) Dietary Pattern Is Associated with Reduced Incidence of Metabolic Syndrome in Children and Adolescents. J Pediatr-Us. 2016;174:178-+.
    
    

57. 57.Ping-Delfos WLS, Beilin LJ, Oddy WH, Burrows S, Mori TA. Use of the Dietary Guideline Index to assess cardiometabolic risk in adolescents. Brit J Nutr. 2015;113(11):1741–1752.
    
    

58. 58.McKeown NM, Dashti HS, Ma J, et al. Sugar-sweetened beverage intake associations with fasting glucose and insulin concentrations are not modified by selected genetic variants in a ChREBP-FGF21 pathway: a meta-analysis. Diabetologia. 2018;61(2):317–330.
    
    

59. 59.Ley SH, Sun Q, Willett WC, et al. Associations between red meat intake and biomarkers of inflammation and glucose metabolism in women. Am J Clin Nutr. 2014;99(2):352–360.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo4OiI5OS8yLzM1MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA5LzE0LzIwMjEuMDkuMDkuMjEyNjMzMzUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

60. 60.Fiorito LM, Marini M, Francis LA, Smiciklas-Wright H, Birch LL. Beverage intake of girls at age 5 y predicts adiposity and weight status in childhood and adolescence. Am J Clin Nutr. 2009;90(4):935–942.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo4OiI5MC80LzkzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA5LzE0LzIwMjEuMDkuMDkuMjEyNjMzMzUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

61. 61.Stanhope KL, Schwarz JM, Havel PJ. Adverse metabolic effects of dietary fructose: results from the recent epidemiological, clinical, and mechanistic studies. Current opinion in lipidology. 2013;24(3):198–206.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/MOL.0b013e3283613bca&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23594708&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

62. 62.Diet, nutrition, and the prevention of chronic diseases. Report of a WHO Study Group. World Health Organ Tech Rep Ser. 1990;797:1–204.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2124402&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

63. 63.Diet, nutrition and the prevention of chronic diseases. World Health Organ Tech Rep Ser. 2003;916:i–viii, 1-149, backcover.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12768890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

64. 64.Krupp D, Doberstein N, Shi L, Remer T. Hippuric acid in 24-hour urine collections is a potential biomarker for fruit and vegetable consumption in healthy children and adolescents. J Nutr. 2012;142(7):1314–1320.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjEwOiIxNDIvNy8xMzE0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDkvMTQvMjAyMS4wOS4wOS4yMTI2MzMzNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

65. 65.Edmands WM, Beckonert OP, Stella C, et al. Identification of human urinary biomarkers of cruciferous vegetable consumption by metabonomic profiling. J Proteome Res. 2011;10(10):4513–4521.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/pr200326k&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21770373&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 

66. 66. Sri Harsha PSC, Wahab RA, Garcia-Aloy M, et al. Biomarkers of legume intake in human intervention and observational studies: a systematic review. Genes Nutr. 2018;13:25–25.
    
    

67. 67.Kühn S, Düzel S, Colzato L, et al. Food for thought: association between dietary tyrosine and cognitive performance in younger and older adults. Psychological Research. 2019;83(6):1097–1106.
    
    

68. 68.Brosnan JT, Brosnan ME. Branched-Chain Amino Acids: Enzyme and Substrate Regulation. The Journal of Nutrition. 2006;136(1):207S–211S.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjEwOiIxMzYvMS8yMDdTIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDkvMTQvMjAyMS4wOS4wOS4yMTI2MzMzNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

69. 69.Nie C, He T, Zhang W, Zhang G, Ma X. Branched Chain Amino Acids: Beyond Nutrition Metabolism. Int J Mol Sci. 2018;19(4).
    
    

70. 70.Maitre L, Lau CE, Vizcaino E, et al. Assessment of metabolic phenotypic variability in children’s urine using (1)H NMR spectroscopy. Sci Rep. 2017;7:46082.
    
    

71. 71.Jakes RW, Day NE, Luben R, et al. Adjusting for energy intake--what measure to use in nutritional epidemiological studies? Int J Epidemiol. 2004;33(6):1382–1386.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyh181&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15333618&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000226198800038&link_type=ISI) 

72. 72.Cole TJ, Lobstein T. Extended international (IOTF) body mass index cut-offs for thinness, overweight and obesity. Pediatr Obes. 2012;7(4):284–294.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.2047-6310.2012.00064.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22715120&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263335.atom)