The prevalence of insulin resistance (IR) in children and adolescents is increasing around the world(Reference Viner, Segal and Lichtarowicz-Krynska1–Reference Valerio, Licenziati and Iannuzzi5). Although there is no universally accepted definition of IR, all studies evaluating the prevalence of IR indicate that more than one out of three obese children or adolescents display IR. For instance, in Greece, the prevalence of IR was found to be 9·2 % (2·9 % in normal-weight, 10·5 % in overweight and 31·0 % in obese children), using the threshold of homeostasis model assessment–insulin resistance index (HOMA-IR) > 2·10 (i.e. 97·5th percentile in normal-weight participants)(Reference Manios, Moschonis and Kourlaba2). In Italy, an IR prevalence of 3·0 % in normal-weight and 40·8 % in obese children (using HOMA-IR > 2·5 as the threshold) was reported(Reference Valerio, Licenziati and Iannuzzi5). In USA, the prevalence of IR was 3·1 % in normal-weight, 15·0 % in overweight and 52·1 % in obese adolescents (using HOMA-IR > 4·39 as a threshold)(Reference Lee, Okumura and Davis4). Finally, a recent study conducted among Bolivian children and adolescents with obesity revealed an IR prevalence of 39·4 % by using a threshold of HOMA-IR > 3·5(Reference Caceres, Teran and Rodriguez3).
IR in children and adolescents has been associated with CVD and metabolic disorders such as hypertension, dyslipidaemia, hepatic steatosis and endothelial dysfunction(Reference Chiarelli and Marcovecchio6, Reference Lee7). All of these risk factors can track into adulthood, increasing the risk of cardiovascular morbidity and mortality(Reference Raitakari, Juonala and Kahonen8, Reference Margolis9). Therefore, it is important to determine factors associated with IR in order to design and implement appropriate preventive programmes.
Although our genetic background has not changed, changes in environmental parameters (such as abundance of food and sedentary lifestyle) have triggered the expression of genes towards obesity and diabetes – the same genes that once helped our ancestors survive periods of food shortage(Reference Barness, Opitz and Gilbert-Barness10–Reference Kahn, Imperatore and Cheng12). Unhealthy dietary patterns and lack of physical activity (PA) seem to be the most important risk factors(Reference Lee7, Reference Mrdjenovic and Levitsky13, Reference Gutin, Johnson and Humphries14). High energy intake coming from increased consumption of simple carbohydrates and dietary fat, especially saturated and trans fatty acids, in conjunction with low consumption of foods rich in dietary fibre, seems to play a role in the early appearance and development of IR in childhood – even in normal-weight children(Reference Artz and Freemark15). The limited data available on food consumption and IR indicate that increased consumption of wholegrain cereals, dairy products (especially low-fat ones), some fish, fruits and all types of vegetables is inversely associated with IR(Reference Hirschler, Ruiz and Romero16–Reference Lara-Castro and Garvey19), while consumption of sugar-sweetened beverages (SSB) and energy-dense foods such as fast foods seems to be positively associated with IR(Reference Bremer, Auinger and Byrd20, Reference Isganaitis and Lustig21).
The aim of the present study was to evaluate the relationship between the consumption of selected food groups and IR, with an emphasis on SSB, in European adolescents.
Methods
Research design
The HELENA-CSS (Healthy Lifestyle in Europe by Nutrition in Adolescence Cross-Sectional Study) is a multicentre investigation carried out in ten European cities: Athens (Greece), Dortmund (Germany), Ghent (Belgium), Heraklion (Greece), Lille (France), Pécs (Hungary), Rome (Italy), Stockholm (Sweden), Vienna (Austria) and Zaragoza (Spain). The main aim of the HELENA-CSS was to obtain reliable and comparable data on a broad battery of relevant nutrition- and health-related parameters such as dietary intake, anthropometry, PA, fitness, haematological and biochemical indices(Reference Moreno, De Henauw and Gonzalez-Gross22). Data collection from the HELENA-CSS took place in 2006–2007.
All participants were recruited at schools and met the general HELENA-CSS inclusion criteria: age range 12·5–17·5 years, not participating simultaneously in another clinical trial, being free of any acute infection lasting less than 1 week before inclusion and having information on weight and height(Reference Mesa, Ortega and Ruiz23). The present analyses were conducted in a subset of the original sample (n 546) with complete data on glucose, insulin and FFQ. This subset was representative of the original sample in terms of gender, age and BMI. Ethics committees from each country approved the HELENA-CSS protocols(Reference Beghin, Castera and Manios24).
Physical examination
The anthropometric methods followed in the HELENA-CSS were described in detail by Nagy et al.(Reference Nagy, Vicente-Rodriguez and Manios25). Weight was measured in underwear and without shoes with an electronic scale (type SECA 861) to the nearest 0·05 kg, and height was measured barefoot in the Frankfort plane with a telescopic height-measuring instrument (type SECA 225) to the nearest 0·1 cm. BMI (kg/m2) was calculated as body weight (in kilograms) divided by the square of height (in metres). Age- and sex-standardized BMI cut-off points according to the International Obesity Task Force were used to define normal weight, overweight and obesity(Reference Cole, Bellizzi and Flegal26). Pubertal stage was recorded by a researcher of the same sex as the child, after brief observation according to Tanner and Whitehouse(Reference Tanner and Whitehouse27).
Physical activity assessment
More details regarding PA assessment can be found in Hagstromer et al.(Reference Hagstromer, Bergman and De Bourdeaudhuij28). In brief, PA was measured by using an accelerometer (Actigraph MTI, model GT1M; Manufacturing Technology Inc., Fort Walton Beach, FL, USA) placed on each individual for several days. The monitor was secured underneath clothing at the lower back using an elastic belt and was worn for seven consecutive days. Adolescents were also instructed to wear the accelerometer during all time awake and only to remove it during water-based activities and sleep time. It was initialized as described by the manufacturer and a 15 s epoch was used. The sum of accelerations was transformed into counts. Low PA was considered when the mean of time spent in activity was from 500 to 1999 counts. Moderate PA was considered when the mean of time spent in activity was from 2000 to 3999 counts. Vigorous PA was considered when the mean of time spent in activity was more than 4000 counts. The moderate-to-vigorous physical activity (MVPA) represents the time spent on at least 2000 counts or more for PA per d.
Dietary assessment
Food consumption
Eating habits were assessed using a mini FFQ from the Health Behaviour in School-Aged Children (HBSC) study. The frequency of consumption of selected food items was recorded by asking the respondent how many times weekly he/she usually eats or drinks the following: fruits, vegetables, sweets (candy or chocolate), coke or other soft drinks that contain sugar (SSB), diet coke or diet soft drinks, low-fat/semi-skimmed milk, whole-fat milk, cheese, other milk products (e.g. yoghurt, chocolate milk, pudding, quark), cereals (e.g. cornflakes, muesli, choc pops), white bread, brown bread, crisps, chips and fish. The response categories were ‘never’, ‘less than once a week’, ‘once a week’, ‘2–4 days a week’, ‘5–6 days a week’, ‘once a day, every day’ and ‘every day, more than once’. The particular food items were selected as indicators of fat, sugar, Ca and dietary fibre intake. In a validity study performed in Belgium, comparison of the FFQ with 7d food diaries showed no overestimation for soft drinks(Reference Vereecken and Maes29). No specific quantities were recorded; therefore, collected data were only used for assessing the frequency of consumption.
Blood samples
Serum concentrations of glucose and insulin were measured after an overnight fast. The HOMA-IR was calculated as [fasting insulin (μIU/ml) × fasting glucose (mmol/l)]/22·5 (to convert fasting insulin values in μIU/ml to pmol/l, multiply by 6·945)(Reference Matthews, Hosker and Rudenski30). A detailed description of the blood analysis has been reported elsewhere(Reference Gonzalez-Gross, Breidenassel and Gomez-Martinez31).
Statistical analysis
Normally distributed continuous variables are expressed as mean values and standard deviations, while skewed variables are reported as median (25th, 75th percentile). Normality of distribution was evaluated through the Kolmogorov–Smirnov test. HOMA-IR was not normally distributed and thus log-transformed values were used. Categorical variables are summarized as relative frequencies and percentages. Associations between categorical variables were tested using the χ 2 test. The associations between the continuous and binary variables (i.e. sex) were evaluated through Student's t test or the Mann–Whitney test when the former were normally or skewed distributed, respectively. Comparisons of log-transformed HOMA-IR values among the categories of food group intake were performed by using one-way ANOVA, after testing for equality of variances. The results are presented as geometric means and 95 % confidence intervals. Bonferroni correction was used to account for increase in type I error due to multiple comparisons.
Multiple linear regression analysis was conducted in order to determine the association of selected food groups with HOMA-IR after adjusting for sex, Tanner stage, total energy intake, PA and BMI percentile. Food groups entered in the model were those found to be significantly associated with HOMA-IR at a univariate level. The results are presented as β coefficients and 95 % confidence intervals. Stratified analysis by sex was also conducted. P values < 0·05 from two-sided hypotheses are considered as statistically significant. The SPSS statistical software package version 18·0 (SPSS Inc., Chicago, IL, USA) was used to conduct all statistical analyses.
Results
Table 1 presents descriptive statistics of the anthropometric parameters, fasting glucose levels, fasting insulin levels, HOMA-IR, PA, total energy intake, carbohydrate intake and fat intake for the total study population and by sex. Median (25th, 75th percentile) HOMA-IR was 0·62 (0·44, 0·87) and this was significantly higher among girls than boys. Moreover, fasting insulin levels were found to be significantly higher among girls than boys, while total energy, fat and carbohydrate intakes were found to be significantly lower in the former compared with the latter (P < 0·001).
P25, 25th percentile; P75, 75th percentile; MVPA, moderate-to-vigorous physical activity; HOMA-IR, homeostasis model assessment–insulin resistance index.
Values were significantly different compared with those for males:*P < 0·05.
†Data are presented as median and (P25, P75).
‡Data are presented as mean and sd.
§Data are presented as n and %.
∥Data are presented as %.
¶MVPA represents time spent on at least 2000 counts or more for physical activity per d.
Table 2 presents the frequency of consumption of several food groups in the study population, as well as the geometric means (95 % CI) of HOMA-IR for each consumption category of the various foods selected in the present study. The most common frequency of consumption for fruits, vegetables, cheese and milk products except for whole-fat and skimmed milk was found to be 2–6 times/week, while the most common frequency of consumption for cereals, white and brown bread, whole-fat milk, skimmed milk, chips, soft drinks, sweets and fish was ≤1 time/week. It was found that HOMA-IR levels increased as the consumption of white bread increased. However, a statistically significant difference was detected only between the means of HOMA-IR in the very low (≤1 time/week) and very high (>1 time/d) consumption categories of white bread (P < 0·05). Moreover, it was detected that the mean HOMA-IR was statistically significantly higher among adolescents consuming brown bread ≤1 time/week compared with those consuming 2–6 times/week (P < 0·01). Finally, the analysis showed that HOMA-IR was also higher in the adolescents consuming SSB >5 times/week compared with those consuming less frequently, although a statistically significant difference was detected between those consuming SSB 5–6 times/week and 2–4 times/week (P = 0·049). No other statistically significant associations were detected between HOMA-IR levels and the other food groups.
HOMA-IR, homeostasis model assessment–insulin resistance index; SSB, sugar-sweetened beverages.
† Log-transformed values.
Association between frequency of consumption of selected food and HOMA-IR was significant, based on one-way ANOVA after Bonferroni correction to account for increase in type I error due to multiple comparisons: *P ≤ 0·05, **P < 0·01.
Table 3 illustrates the results of the multiple linear regression model using the log-transformed values of HOMA-IR as dependent variable and the consumption frequency of SSB, white and brown bread as independent variables, controlling for sex, Tanner stage, total energy intake, PA and BMI percentile. The results of this analysis indicate that among the three food groups included in the model, only the consumption of SSB was significantly associated with HOMA-IR. In particular, it was found that HOMA-IR levels were higher among adolescents consuming SSB 5–6 times/week and ≥1 time/d compared with those consuming ≤1 time/week by 0·281 and 0·191 units, respectively (P = 0·009 and 0·046, respectively).
HOMA-IR, homeostasis model assessment–insulin resistance index; SSB, sugar-sweetened beverages.
† Log-transformed values.
‡ After controlling for sex, Tanner stage, total energy intake, physical activity and BMI percentile.
Similar findings were detected when stratified analysis by sex was conducted (Tables 4 and 5).
HOMA-IR, homeostasis model assessment–insulin resistance index; SSB, sugar-sweetened beverages.
† Log-transformed values.
‡ After controlling for sex, Tanner stage, total energy intake, physical activity and BMI percentile.
HOMA-IR, homeostasis model assessment–insulin resistance index; SSB, sugar-sweetened beverages.
† Log-transformed values.
‡ After controlling for sex, Tanner stage, total energy intake, physical activity and BMI percentile.
Discussion
Although there are plenty of studies examining the association between SSB consumption and obesity, type 2 diabetes and other CVD risk factors(Reference Malik, Popkin and Bray32), limited data are available regarding the relationship between SSB consumption and IR(Reference Bremer, Auinger and Byrd20). To the best of our knowledge, the present work is the first large European study examining this association.
The present findings indicate that SSB consumption is significantly associated with HOMA-IR levels even after controlling for white and brown bread consumption, obesity indices (i.e. BMI percentiles), PA, total energy intake and other potential confounders. In particular, it was found that adolescents with SSB consumption equal to or higher than once daily had higher HOMA-IR than adolescents consuming SSB less than once weekly. Similar association was detected in both genders.
This result could be partly explained by the contribution of SSB consumption to weight gain(Reference Malik, Popkin and Bray32, Reference Malik, Schulze and Hu33). However, the fact that the increased SSB consumption is associated with higher levels of IR-related indices, even after adjusting for BMI percentiles, indicates that an independent effect may also stem from the large quantities of rapidly absorbable carbohydrates used to flavour these beverages. In particular, it has been shown that SSB consumption raises blood glucose and insulin concentrations rapidly and dramatically(Reference Janssens, Shapira and Debeuf34). Therefore, when consumed in large amounts, SSB contribute to a high dietary glycaemic load, which has been shown to induce glucose intolerance and IR(Reference Ludwig35).
Current evidence suggests that a high dietary intake of fructose that comes from either sucrose-sweetened beverages or other foods may lead individuals to develop IR and metabolic syndrome through several mechanisms(Reference Montonen, Jarvinen and Knekt36–Reference Meyer, Kushi and Jacobs38). First, it has been suggested that high intake of fructose is associated with higher concentrations of C-peptide, a marker of insulin secretion and IR(Reference Wu, Giovannucci and Pischon39). It is already known that fructose does not stimulate insulin secretion and also reduces circulating leptin concentrations. The combined effects of lowered leptin and insulin concentrations could induce the likelihood of weight gain and its associated metabolic sequelae(Reference Elliott, Keim and Stern40). Finally, fructose seems to induce weight gain due to inadequate compensation of energy intake from solid foods when SSB are ingested(Reference Bray41).
Although our findings indicate a low consumption of SSB in European adolescents compared with other similar populations (i.e. US), these findings are in accordance with those reported from a similar study carried out in the USA among 3831 students in 6th to 12th grade(Reference Bremer, Auinger and Byrd20). Bremer et al. observed that increased SSB consumption was independently associated with increased HOMA-IR, LDL and TAG concentrations and decreased HDL concentrations. The innovative finding of that study was that the increased SSB consumption was related to increased HOMA-IR and TAG concentrations in boys but to BMI and waist circumference in girls(Reference Bremer, Auinger and Byrd20). Moreover, a recent meta-analysis revealed that adults consuming more than one serving of SSB daily had higher risk of developing type 2 diabetes than those consuming less than one serving monthly, indicating that high SSB consumption is a risk factor for metabolic disorders not only in childhood but also in adulthood(Reference Malik, Popkin and Bray42).
The current study is cross-sectional and hence no causal relationship between SSB consumption and IR can be extracted. Recent study has shown that the amount of fructose and/or glucose additives in beverages is associated with IR development(Reference Basciano, Federico and Adeli37). However, this association was not evaluated in the current study due to lack of related data. In addition, in the current study SSB included only coke or other soft drinks that contain sugar, while in general SSB is a wider group including additionally fruit drinks, sweetened teas, sport drinks, etc. Therefore, the effect of SSB on IR may have been underestimated. Finally, an important limitation of the analysis is that food frequency consumption data were collected by a simple FFQ, even though more detailed 24 h recall data are available in the present study. However, the use of standardized methodologies and tools for the collection of dietary data within HELENA-CSS strengthens the value of our findings. At this point it should be highlighted that this is an exploratory analysis only, and these associations should be further explored using more detailed dietary data.
Conclusions
Frequent consumption of SSB of more than once daily seems to increase fasting glucose and HOMA-IR levels, which could lead to an increased risk for early development of type 2 diabetes in adolescents. If further research findings on the potential effect of SSB are consistent with those of the current study, then it is important that educational programmes aiming to improve adolescents’ dietary habits, including the reduction of SSB consumption, are designed and implemented in European adolescent populations in order to prevent the increment in type 2 diabetes prevalence.
Acknowledgements
The HELENA Study was carried out with the financial support of the European Community Sixth RTD Framework Programme (Contract FOOD-CT-2005-007034). The content of this article reflects the authors’ views only, and the European Community is not liable for any use that may be made of the information contained herein. The authors declare that there are no conflicts of interest. The authors’ contributions are as follows: K.K. is the main author; E.G. is the second author; D.J.-P. helped with the discussion; S.D.H. provided guidelines on cut-off points for statistics; M.G.-G. provided blood data; M.S. provided the analysis of anthropometric statistics; F.G. and D.M. are members of the steering committee; L.A.M. is the study coordinator; A.K. is the data provider for southern Europe; C.G. is the data provider for northern Europe; M.K. is a member of the steering committee; Y.M. is the database provider for the Greece samples. The authors wish to thank the members of the Core Group of the HELENA Study.