- Open Access
Development and validation of dietary and lifestyle insulinemic indices among Iranian adult population
Nutrition & Metabolism volume 19, Article number: 5 (2022)
There is no study regarding developing a valid index to predict insulin-related disorders in the Iranian population based on their dietary habits and lifestyle. In the current study, we aimed to develop and validate insulinemic potential indices of diet and lifestyle in Iranian adults.
In this cross-sectional study, we analysed data of 1063 men and women aged ≥ 25 years among participants of the examination three of Tehran lipid and glucose study (TLGS) (2006–2008). Dietary intakes were assessed using a valid semi-quantitative food frequency questionnaire. Dietary and lifestyle indices were developed using stepwise linear regression analysis based on dietary intakes, body mass index, and physical activity data. Fasting serum insulin concentration and homeostatic model assessment for insulin resistance (HOMA-IR) were used as biomarkers of hyperinsulinemia (HI) and insulin resistance (IR). Validation analyses were performed in examination four of TLGS.
We developed four indices related to insulin homeostasis, including the dietary index for HI (DIH), the dietary index for IR (DIR), the lifestyle index for HI (LIH), and the lifestyle index for IR (LIR). Based on multivariable-adjusted models, the relative values of the biomarker in subjects in the highest quartile of indices were 45% for LIH (95% CI 1.36–1.55, Ptrend < 0.001), 28% for DIR (95% CI 1.13–1.42, Ptrend = 0.019), and 51% for LIR (95% CI 1.41–1.61, Ptrend < 0.001), higher than those in the reference quartile, respectively.
We designed and validated indices to determine the insulin potential of diet and lifestyle for the Iranian population, according to Iran’s demographic and dietary intake characteristics.
Insulin is one of the most important hormones maintaining metabolic homeostasis. Insulin affects different tissues in the body, including adipose tissue, liver, muscle, brain, bone , kidneys , and vascular system [3, 4]. Insulin-related disorders, including hyperinsulinemia (HI) and insulin resistance (IR), has been known as an early indicator and central place in metabolic dysfunction [5, 6]. If these initial disorders are left uncontrolled, they will contribute to the majority of chronic diseases such as type 2 diabetes (T2D), metabolic syndrome (MetS) .
Information on IR and HI statistics worldwide and in Iran is limited; however, considering the prevalence of T2D and MetS, as the two diseases most closely related to IR and HI, clarifies the current situation. According to the latest estimates, the prevalence of MetS is about 14–32% worldwide  and 30.8% in Iran . The prevalence of diabetes in Iran in 2021 was 9.1% among adults, and it is projected to reach 10.6% by 2045. It is also estimated that the proportion of undiagnosed T2D in the Iranian adult population is about 36% . This high rate of chronic diseases related to IR and HI globally, especially in Iran, requires special attention to improve the prevention and control of these disorders in the early stages. Lifestyle is the most important modifiable risk factor affecting insulin-related disorders development .
Diet is well-established to be an essential part of lifestyle . Although specific dietary factors have been shown to affect insulin secretion and action [11,12,13], dietary patterns or dietary indices, including several dietary factors and cover the complex interactions between nutrients and foods, can provide a better estimate about the relationship between diet and risk of chronic disease [14, 15]. Other lifestyle factors associated with HI and IR include body mass index (BMI) and physical activity (PA) [16,17,18]. Given that each of these lifestyle factors is a great predictor of metabolic disorders per se, combining diet with them is likely to be more predictive of HI and IR risk than these factors alone.
Recently, Tabung et al. has developed four indices to assess the insulin potential of diet and lifestyle, including empirical dietary and lifestyle indices for HI and IR . To create these indices, the concentration of connecting peptide (C-Peptide) and triglycerides (TGs) / high-density lipoprotein cholesterol (HDL-C) ratio were used as biomarkers for predicting HI and IR respectively. However, this study had some limitations. Questionnaires for medical history information were self-administered, reducing reporting accuracy [20, 21]. The study mentioned above also used blood data of participants in the Nurses’ Health Study (NHS), all of whom are women and highly educated; so, the generalizability of its results to the general population may be limited. Latterly two insulinemic indices were proposed by Mazidi et al. . Similar to Tabung et al., Mazidi et al. have used C-Peptide and TG/HDL-C as biomarkers to develop the indices; however, they derived nutrient patterns, which are limited in covering interactions between different food components, such as those observed in the dietary pattern. Also, in Mazidi et al. study, the dietary intakes have assessed using a single 24-h recall, which may not be appropriate for determining the usual diet or irregularly consumed foods.
Based on the available literature on this topic, there is no study on developing a valid index to predict insulin-related disorders in the Iranian population based on their dietary habits, patterns, and lifestyle. Therefore, this study aims to develop and validate insulin potential indices of diet and lifestyle in the Iranian population.
We performed the current study within the Tehran Lipid and Glucose Study (TLGS), an ongoing population-based prospective study started in 1999, which was implemented to determine the risk factors for non-communicable diseases (NCDs) among a representative urban population in Tehran, including 15,005 participants aged ≥ 3 years. The TLGS data have been collected at three-year intervals; details about TLGS have been reported elsewhere .
In the third examination of the TLGS (2006–2008), of 12,523 participants, 3462 were randomly selected for dietary assessment. For this study, 1286 men and women aged ≥ 25 years were selected with complete data on fasting blood sugar (FBS) and fasting serum insulin concentration. We excluded individuals based on at least one of the following criteria: history of myocardial infarction, cerebrovascular accident or cancer (n = 18), diabetic patients (n = 105), reported daily energy intakes less than 800 or more than 4200 kcal/d in men and out of the range 500–3500 kcal/d in women (n = 80), and pregnancy or lactation (n = 23); it is noticeable that some individuals fell into more than one exclusion category. Ultimately, 1063 participants remained for the analysis. For validation analysis in the fourth examination of the TLGS (2009–11), with applying the exclusion criteria mentioned above, 758 participants remained for the final analysis.
Dietary intake assessment
We use a valid and reliable semi-quantitative 168-item food frequency questionnaire (FFQ) to assess dietary intakes at baseline. The reproducibility and validity of the FFQ have been assessed previously [24, 25]. The consumption frequency of each food item was collected during the previous year, as daily, weekly, or monthly, in a face-to-face interview by trained and skilled dieticians. Portion sizes of consumed foods were reported in household measures and then transformed into grams scale. The energy contents of consumed foods were calculated using the United States Department of Agriculture (USDA) food composition table (FCT). Furthermore, the Iranian FCT was used to analysis of the local food items that were not available in USDA FCT.
Physical activity assessment
The modifiable activity questionnaire (MAQ) was used to assess individuals’ physical activity levels. This questionnaire has previously been modified and validated for the Iranian adult population . Participants were requested to report the frequency and time spent on activities in four intensity levels, light, moderate, heavy, and very heavy, during the previous year, according to a list of daily life common activities. Ultimately, physical activity levels were expressed as metabolic equivalent hours per week (MET-h/wk.).
Demographic and lifestyle-related measurements
Using a pretested questionnaire, trained interviewers collected information on age, sex, medical history, medication use, occupation status, and smoking. The participant’s body weight was measured to the nearest 0.1 kg using digital scales (model 707, Seca, Hamburg, Germany) while lightly clothed and barefoot. Subjects’ height was measured to the nearest 0.5 cm using a stadiometer (model 208 Portable Body Meter Measuring Device; Seca) while subjects were standing without shoes. BMI was computed as weight (kilograms) divided by height (meters squared). Waist circumference (WC) was measured to the nearest 0.1 cm using a non-elastic tape meter, between the lowest chest ribs and the iliac crest at the umbilicus level, over light clothing, without any pressure on the body skin.
Blood samples were taken and transferred into vacutainer tubes between 7:00 and 9:00 a.m, after a 12–14-h overnight fast, while subjects were in a sitting position. Blood samples were centrifuged within 30 to 45 min of collection. All biochemical analyses were performed using a Selectra 2 auto-analyzer at the TLGS research laboratory on the day of blood collection. FBS was measured using an enzymatic colorimetric method with glucose oxidase. Inter- and intra-assay CVs were both 2.2% for FBS . TGs levels were measured using the enzymatic colorimetric method with glycerol phosphate oxidase. Inter- and intra-assay CVs for TGs were 0.6 and 1.6%, respectively. Serum HDL-C was measured after precipitation of the apolipoprotein B-containing lipoproteins with phosphotungstic acid. Fasting Insulin was measured via electrochemiluminescence immunoassay (ECLIA), using Roche Diagnostics kits and Roche/Hitachi Cobas e-411 analyzer (Gmbh, manhim, Germany). Inter- and intra-assay coefficient variations for insulin were 1.2 and 3.5, respectively.
HOMA-IR = FBS (mmol/L) × serum insulin (μU/mL) / 22.5
Visceral adiposity index (VAI): this index was calculated for men and women as below:
Males: VAI = (WC (cm) / (39.68 + (1.88 × BMI (kg/m2))) × (TG (mmol/l) / 1.03) × (1.31 / HDL (mmol/l))
Females: VAI = (WC (cm) / (39.58 + (1.89 × BMI (kg/m2))) × (TG (mmol/l) / 0.81) × (1.52 / HDL (mmol/l))
Waist residual BMI: we regressed WC on BMI to obtain BMI-independent waist circumference values, calculated as differences between each individual’s WC and the WC predicted by BMI.
Development of indices
In the first stage, all food groups and subgroups that their association with the risk of IR, HI, T2D, and MetS have been considered in previous studies (especially based on studies conducted on the Iranian population) were extracted. Finally, 38 food groups were considered for development analysis, including red meat, processed meat, poultry, fish, egg, organ meats, pizza, broth, low-fat dairy, high-fat dairy, doogh, sweetened dairy, refined grain, whole grain, legumes, nuts, dried fruits, fruit juices, other fruits, lemon juices, starchy vegetables, tomato, pickles, leafy vegetables, garlic and onion, other vegetables, fried potato, animal fat, liquid oils, butter, mayonnaise, olives, sweet snacks, artificial beverages, tea, salt, coffee, herbs. The intakes of these food groups or food items were calculated as daily servings per 1000 kcal energy intake.
Stepwise linear regression analysis was conducted to derive dietary patterns associated with HI and IR, using fasting insulin concentration and HOMA-IR index (log-transformed), respectively, as the response variables and the 38 food items mentioned earlier as the predictor. The significance level was considered at P = 0.1 for entry into and retention in the final model. To develop lifestyle indices, BMI and PA were also added to the final model. Ultimately the linear β-regression coefficients (as the weights) and the coefficient of determination (R2), obtained in the final model of stepwise regression analysis, were reported for each food group.
Components derived for each index
The dietary index for HI (DIH), the dietary index for IR (DIR), were developed based on different food items. The DIH had ten components, including refined grains, pickles, doogh (a beverage made with yogurt), sweetened beverages, fish, and lemon juice which directly related to the serum insulin concentration, and broth, red meat, high-fat dairy, starchy vegetables, which inversely related to this biomarker. The DIR was composed based on twelve food items, including six items, which had a direct relationship with the HOMA-IR index: pickles, refined grains, doogh, lemon juice, sweetened beverages, and fish, as well as six items, which had an inverse relationship with this biomarker: starchy vegetables, snacks, low-fat dairy, broth, red meat, high-fat dairy. The items of each food group are defined in detail in Additional file 1: Table S1.
The BMI, PA, and food items were considered for developing the lifestyle indices, including the lifestyle index for HI (LIH), and the lifestyle index for IR (LIR). The LIH had eight components, including BMI, refined grains, doogh, and fish (as positive components associated directly with serum insulin concentration), PA, low-fat dairy, high-fat dairy, and starchy vegetables (with inverse relationship with the serum insulin concentration). The LIR was created according to seven factors, including BMI, refined grains, and doogh, directly related to the HOMA-IR index, low-fat dairy, PA, starchy vegetables, and high-fat dairy that inversely related with HOMA-IR index.
Calculations of indices
Weighted indices were calculated as follows: Each dietary item was converted into serving size per 1000 kcal of energy intake. Then each component, including dietary items, BMI, and physical activity (MET.h/wk.) multiplied by their weights (β coefficient) reported in Table 1, and values for all components were summed as the final score. Food groups with their items used in developed indices were defined in Additional file 1: Table S1. We also calculated unweighted forms of indices (supplementary file).
To evaluate the validity of the indices for predicting the fasting serum insulin and the HOMA-IR levels based on data from the third examination of the TLGS, the developed indices were categorized into quartiles. Concentrations of the biomarkers were back-transformed to their original units. The mean and 95% confidence intervals (95% CI) of serum insulin concentration and HOMA-IR value in the first quartile of each index were assigned as the reference value. The relative values (RVs) of other quartiles were computed by dividing each quartile value by the reference value. Then the adjusted RV and 95% CI of serum insulin and HOMA-IR across quartiles of the weighted insulin-related indices were calculated using univariate linear regression analysis in two adjusted models for potential confounding variables: (1) adjusted for age and sex, and (2) adjusted for all confounding variables, including age, sex, energy intake, smoking, occupation status, medication use, as well as PA and VAI, TGs to HDL-C ratio, and waist residual BMI.
The multivariable linear regression analysis was used to compute the trend of changes in the RVs of fasting serum insulin concentration and the HOMA-IR across quartiles of each weighted index. These analyses were repeated in the next examination of the TLGS (examination 4) to confirm the validity of the weighted DIH, DIR, LIH, and LIR in an independent sample size.
To evaluate the sensitivity of the developed indices, we derived three subgroups with (1) excluding individuals with a family history of T2D, (2) excluding smokers, and (3) including T2D patients in the baseline analysis. The RVs and 95% CI of serum insulin concentration and HOMA-IR were determined for each quartile of related indices, based on two above-mentioned adjusted models using univariate linear regression analysis.
Additionally, we performed a stratified analysis to assess the distribution of the mean fasting serum insulin concentrations across quartiles of DIH and the mean of HOMA-IR across quartiles of DIR by creating two joint categories of BMI and PA as follows: active obese individuals (BMI ≥ 25 kg/m2 and PA ≥ median PA), and non-active obese individuals (BMI ≥ 25 kg/m2 and PA < median PA).
Based on rules of thumb for determining sample size firstly proposed by Roscoe in 1975 that says in multivariate research (including multiple regression analyses), the sample size should be several times (preferably ten times or more) larger than the number of variables in the study [27, 28]. The minimum sample size for development of DIR and DIH (38 variables) and LIR and LIH (40 variables) were 380 and 400 participants, respectively. However, we developed these indices among 1063 participants in the third examination of TLGS. Validation study and other regression analysis were also conducted among sample sizes larger than estimated above.
All analyzes were performed using SPSS software version 20. The normality of the data was assessed using the Kolmogorov–Smirnov test and histogram charts. Stepwise-linear regression was used to develop insulin-related indices, including DIH, DIR, LIH, and LIR. Data are presented as the mean ± SD or the median (interquartile range) for continuous variables and percentages or numbers (percentages) for categorical variables across quartiles of DIH and DIR. Chi-square and linear regression were used to test the trend of qualitative and quantitative variables across these indices.
We computed the Pearson and Spearman correlation coefficients to evaluate the correlation between calculated indices including DIH, DIR, LIH, LIR, and also an unweighted form of these indices with and serum insulin levels and HOMA-IR in the third and fourth examination of TLGS. Furthermore, the correlation of the indices mentioned above with similar indices developed in previous studies, including insulin index, insulin load, the empirical dietary index for HI, the empirical dietary index for IR, the empirical lifestyle indices for HI, and the empirical lifestyle indices for IR, was determined.
Validation, sensitivity, and subgroup analyses were conducted using univariate linear regression analysis adjusted for potential confounders and computed the RVs and 95% CI of serum insulin concentration and HOMA-IR values across quartiles of weighted insulin indices. Then we calculated the P for trend between these RVs across quartiles of each insulin indices using multivariable linear regression. All regression analyses were adjusted for age, sex, energy intake, smoking, occupation status, medication use, PA and VAI (only for DIH and DIR), TGs to HDL-C ratio, and waist residual BMI (only for LIH and LIR). P values < 0.05 were considered as statistically significant.
We also calculated the power of the study for multiple linear regression using G-power software version 220.127.116.11 by considering a two-tailed test at a 5% level of significance, sample size = 1063, 38 predictors for DIR and DIH, 40 predictors for LIR and LIH, and R2 of final model equal to 0.049, 0.049, 0.261, 0.265 for DIR, DIH, LIR, and LIH, respectively. We observed maximum power for all indices.
Components of the four insulinemic dietary and lifestyle indices, including DIH, DIR, LIH, and LIR, with their specific weights and R2 values, were reported in Table 1. Based on the relationship of components with serum insulin and HOMA-IR levels, as response biomarkers for each index, the components were divided into two categories, including (1) food groups that directly related to mentioned biomarkers and (2) food groups that showed an inverse relationship to these biomarkers. The first group components directly related to serum insulin and HOMA-IR levels were given positive weights, while components with inverse relationship with the biomarkers were assigned negative weights. Detailed information about food groups with their items used in developed indices were provided in Additional file 1: Table S1.
Table 2 presents the participants’ baseline characteristics across quartiles of the DIH and DIR. Serum insulin, HOMA-IR, WC, TGs, TG to HDL-C ratio, VAI, and percentage of men were significantly increased across quartiles of DIH and DIR (P for trend = < 0.001). In contrast, participants’ age and HDL-C concentration were decreased across quartiles of DIH and DIR (P for trend = < 0.05).
Pearson and Spearman’s correlations coefficients of the insulinemic dietary and lifestyle indices with serum insulin and HOMA-IR based on examinations 3 and 4 of TLGS were shown in Table 3. In examination 3 (baseline), a significant moderate correlation was observed between all developed indices (DIH, DIR, LIH, and LIR) and serum insulin and HOMA-IR levels. In contrast, in examination 4, this significant correlation was only observed for LIH with serum insulin and LIR with HOMA-IR. In Table 3, we also reported the results on the correlation of previously designed similar insulin-related indices, including insulin index, insulin load, empirical dietary indices for HI and IR, empirical lifestyle indices for HI and IR with serum insulin levels, and HOMA-IR. Although the Spearman correlation coefficients showed that three of four above mentioned indices, including empirical lifestyle index for HI, insulin index, and insulin load, were significantly correlated with serum insulin in both examinations, Pearson correlation was significant only for the including empirical lifestyle index for HI score with serum insulin. Furthermore, a weaker correlation was found between HOMA-IR and empirical dietary index for insulin resistance in examination 3 and HOMA-IR and empirical lifestyle index for insulin resistance in examinations 3 and 4.
The Relative values and 95% CIs of serum insulin concentration and HOMA-IR across quartiles of the DIH, DIR, LIH, and LIR based on the two examinations of TLGS (3 and 4) were shown in Table 4. RVs of serum insulin concentration and HOMA-IR level followed statistically significant upward trends across quartiles of most indices in both examinations. However, a significant trend for biomarkers RVs across DIH quartiles was observed only in examination 3.
In examination 3, participants in the highest quartile of DIH and LIH had 22% and 59% higher serum insulin concentrations than those in the first quartile (95% CI 1.15–1.29, and 1.52–1.67, P for trend < 0.001), respectively. Also, based on multivariable-adjusted models, the HOMA-IR index was 28% and 67% higher among individuals in the highest vs. lowest quartile of DIR and LIR (95% CI 1.20–1.36, and 1.59–1.75, P for trend < 0.001), respectively. In examination 4, the trend of the RV of serum insulin level across quartiles of DIH score was not significant; however, an upward trend was observed in mean serum insulin level and HOMA-IR index across quartiles of the DIR, LIH, and LIR. Based on multivariable-adjusted models, the RV of mean serum insulin has significantly increased across quartile of LIH (RV = 1.45; 95% CI 1.36–1.55, Ptrend < 0.001). Also, in fully adjusted model, the RV of HOMA-IR index level have significantly increased across quartiles of DIR (RV = 1.28; 95% CI 1.13–1.42, Ptrend = 0.019) and LIR (RV = 1.51; 95% CI 1.41–1.61, Ptrend < 0.001).
Using the sensitivity analysis, the trend of the RV and 95% CIs of serum insulin across quartiles of DIH and LIH and the RV and 95% CIs of HOMA-IR across quartiles of DIR and LIR has been shown based on three sub-groups, including (1) excluding the smoked subjects, (2) excluding participants with a family history of T2D, and (3) including type 2 diabetic patients at baseline in Table 5. Our results showed that the RV of serum insulin concentration across quartiles has significantly increased across quartiles DIH and LIH based on age and sex adjusted and as well as multivariable-adjusted models. Also, the RV of HOMA-IR index levels has significantly increased across quartiles DIR and LIR based on both models. Furthermore, the stratified analysis showed that the mean serum insulin concentration and HOMA-IR of active-obese participants with the highest scores of DIH and DIR were 23% and 30% higher than those with the lowest scores. Similarly, non-active obese subjects in the highest quartile of DIH and DIR had 28% and 36% greater HOMA-IR values than those in the lowest quartile (Fig. 1).
In the present study, we developed and validated four indices, including DIH and DIR, consisting of a dietary pattern predictive of HI and IR, as well as LIH and LIR encompassing both diet and lifestyle factors that are significantly related to insulin disorders. Considering the importance of PA and obesity in HI and IR [16, 18], we performed a subgroup analysis based on the physically active obese group and physically inactive obese individuals, separately. Our analysis showed that in the physically active obese subgroup, subjects in the highest quartile of DIH and DIR compared to those in the lowest one had 22 and 31% higher serum insulin and HOMA-IR, respectively. Also, in the physically inactive obese subgroup, serum insulin concentration and HOMA-IR levels of individuals in the highest quartile of DIH and DIR were 26% and 36% higher than those in the lowest quartile, respectively. These analyses showed that individuals with higher DIH and DIR scores had higher serum insulin concentration and HOMA-IR levels, independent of BMI and PA.
In recent years, Tabung et al.  and Mazidi et al. , in two separate studies, introduced several indices related to insulin homeostasis status, using serum C-Peptide concentration and TGs to HDL-C ratio as biomarkers of HI and IR, respectively. Similar to our study, Tabung et al. derived four dietary and lifestyle indices with different food groups, BMI, and PA, which indicate the insulinemic potential of diet and lifestyle; however, Mazidi’s study used a pattern of nutrients to develop insulinemic dietary indices. Using nutrients instead of food groups in Mazidi’s study may not accurately estimate the diet-disease relationship due to ignoring the complex interactions of different substances in whole foods . However, it seems that the most important concern about these indices introduced in Mazidi et al. and Tabung et al.’s studies is that they are all developed in the US population, which may limit the generalizability of their results for other populations. In this study, of all food groups that were entered in the stepwise regression analysis, only those with statistically significant (positive or negative) relationships (P < 0.1) with serum insulin level and HOMA-IR were included in the final patterns. In our study, intakes of red meat, broth, high- and low-fat dairy products, starchy vegetables, snacks, and higher PA level showed an inverse relationship with serum insulin level and HOMA-IR, while consumption of refined grains, sweetened beverages, pickle, doogh, lemon juice, fish and BMI directly associated with the biomarkers as mentioned earlier. Our findings are in line with the results of most of the studies that have suggested the reduced risk of HI and IR in individuals with higher intakes of starchy vegetables, low-fat dairy products, and higher PA levels, as well as lower intake of refined grains and sweetened beverages and obesity [30,31,32,33,34,35].
In Iranian food culture, fish is mainly eaten fried form, which this cooking method can increase the intake of trans fats and weaken the beneficial properties of this food item. Such an effect is also noticed in other studies . In contrast, red meat is mainly eaten unprocessed form along with legumes, vegetables, which may mitigate its potentially harmful effects based on results of previous studies . Based on the extracted patterns in this study, higher consumption of doogh (a beverage made from yogurt, which has relatively high salt content), lemon juice, and pickles directly correlated with serum insulin and HOMA-IR levels. Although doogh is a low-calorie beverage, its high salt content can worsen IR through several mechanisms . Moreover, dietary constituents like doogh, lemon juice, and pickles are often eaten with high-carbohydrate, high-calorie meals, meaning that consuming more indirectly indicates higher energy and carbohydrate intake, which raises IR risk.
In our study, lifestyle indices (LIH and LIR) were 2 to 2.5 times more correlated with biomarker levels than dietary indices (DIH and DIR). This higher correlation indicates the higher potential of lifestyle indices to predict the risk of HI and IR. Dietary indices showed a relatively low but significant correlation with serum insulin and HOMA-IR level. The low correlation of these dietary indices with the biomarkers does not mean that the developed indices are not valid; because the serum insulin and glucose level’ specialty in healthy adults, more influenced by genetic, endogenous factors (like regulatory systems), and other lifestyle factors (such as PA and BMI) compared with a usual diet , so the low correlation of some nutritional factors with serum insulin and HOMA-IR level was expectable. Similar findings have been observed in other studies [19, 39], e.g., in the Fung et al. study, a dietary pattern predicting serum C-peptide showed a relatively low correlation (r = 0.23) with the levels of this biomarker; however, the association of this pattern with the risk of colon cancer was still significant . Also, in the Tabung et al. study, the reported correlation coefficients between developed indices and biomarkers were close to ours, and their indices had good predictive power for the risk of HI and IR .
Various analyses confirmed the robustness and validity of all four developed indices in our study. In the third examination of the TLGS study (baseline), the relative serum insulin concentration and HOMA-IR levels increased across the quartiles of DIH, DIR, LIH, and LIR with a significant trend. Similarly, these significant findings were observed in the fourth examination for all indices mentioned above except DIH. In general, dietary indicators are less capable than lifestyle indicators for predicting the response biomarkers. In the validation study, the adjusted RV of fasting insulin in the highest vs. lowest quartile of DIH was 4% higher; however, it showed no significant trend. In contrast, HOMA-IR showed a significant positive trend across quartiles of DIR. Some possible reasons may explain why DIH has shown a weaker predictive ability than DIR. First, DIH covers a smaller range of dietary components than DIR, as DIR, additional to similar items with DIH, includes snacks and low-fat dairy. Second, fasting insulin concentration (as the response variable for DIH) is less affected by diet compared with HOMA-IR (as the response variable for DIR) because FBS as a component of HOMA-IR is more affected by food intakes . Finally, in our study, about 97% of individuals had a positive DIH score (compared to DIR with 75% positive and 25% negative). This may have undermined the DIH’s ability to discriminate between individuals with different serum insulin levels and, consequently, its predictive ability in our study.
Furthermore, we tested the sensitivity of the indices in three different sub-groups by excluding subjects with a family history of T2D, excluding smoked participants, and including diabetic patients in the analysis. Finally, based on each sub-group analysis, the relative values of serum insulin concentration and HOMA-IR significantly increased across quartiles of all developed indices.
We examined the correlation of the developed indices in the unweighted form and the similar indices developed in previous studies with the serum insulin and HOMA-IR levels. Our results showed that even unweighted forms of DIH and DIR correlate with serum insulin concentration and HOMA-IR more than the empirical dietary index for HI, insulin index, and insulin load. Weighted DIH and DIR showed a lower correlation with biomarkers than unweighted forms. One possible reason could be that in unweighted form, all components were equally involved in the correlation between indices and biomarkers regardless of their dietary intakes, e.g., in the weighted form, the effect of some components such as pickles that were high in weight and could increase correlation was attenuated by low dietary intake, The opposite is true for sweetened beverages. Among lifestyle indices, the large difference in weight and R2 of the BMI compared to other items covered their overall impact and significantly affected increasing correlation. When BMI was considered with the same weight as other factors, its effect is attenuated, which reduces the correlation of index with serum insulin and HOMA-IR. Further analysis showed that the DIH, DIR, LIH, and LIR were significantly correlated (r < 0.40) with their unweighted forms (Additional file 1: Table S2, S3).
The indices introduced in the present study can be used in other populations, especially those with similar habits to the Iranian population, such as people in the Caucasus and the Middle East. The main benefit of our results is that researchers who do not have access to insulin biomarkers data can easily estimate the insulinemic potential of diet and lifestyle using only dietary and lifestyle data. More research and studies are still needed on the derived patterns and the food groups (to determine the optimal intakes). If the results obtained in this study are repeated in future studies, clinical applications for these findings can be considered. Physicians and dietitians can assess the diet and lifestyle of patients at risk of developing insulin-related disorders using the indices proposed in this study.
Our study had several strengths and limitations. This work is the first study conducted to develop and validate dietary and lifestyle indices to predict the risk of HI and IR among the Iranian population. We used serum insulin concentration as a response variable to develop DIH and LIH indices, which is a direct indicator for the diagnosis of HI. HOMA-IR index was also used to detect insulin resistance for DIR and LIR development. This index has been widely used in studies among various populations and on several outcomes due to the high correlation with the euglycemic hyperinsulinemic clamp techniques as the gold standard for diagnosing IR . Our study has some limitations. First, due to the lack of another study on measuring the serum insulin concentration with an acceptable sample size in the Iranian adult population, we had to validate the developed indices using participants’ data from examination 4 of the TLGS study. However, this limitation is somewhat negligible because the validation of the indices was done in the fourth examination of the TLGS study with new entries and different sample sizes, in which their data were collected independently of the third examination (baseline). Second, despite using a valid and reliable FFQ for dietary assessment, measurement errors were inevitable in the estimates of dietary intake using FFQ. However, we entered dietary intakes as energy-adjusted into analysis to reduce the impact of measurement errors to the extent possible. Third, we have no complete data of other lifestyle factors such as income, culture, customs, etc. data to determine better the lifestyle index related to HI and IR. Fourth, although we attempted to adjust for major confounders in the present study, the confounding effect of some unknown and unmeasured residual confounders cannot be excluded. In addition, although direct measurement of PA is not common in large population studies, the use of indirect measurement for PA is another limitation that cannot be ignored completely; however, we collected the PA data using a valid questionnaire in a face-to-face interview to enhance the assessment precision.
In the present study, we have designed and validated four insulinemic indices, including DIH, DIR, LIH, and LIR, to determine the insulinemic potential of diet and lifestyle in the Iranian population, according to their demographic and dietary intakes. The findings of this study can be used for research and clinical purposes. These indicators may help identify individuals, groups, and even populations at risk for diabetes, metabolic syndrome, cardiovascular disease, etc. Researchers who do not have access to insulin biomarker data and researchers in other countries, especially in Southwest Asia, which has a food culture and lifestyle close to Iran, can use these indices in their studies. The results of this study need to be confirmed by large prospective studies and randomized clinical trials.
Availability of data and materials
The data underlying this article will be shared at reasonable request to the corresponding author.
Body mass index
Dietary index for HI
Dietary index for IR
Fasting blood sugar
Food composition table
Food Frequency Questionnaire
High-density lipoprotein cholesterol
Homeostatic model assessment for insulin resistance
Lifestyle index for HI
Lifestyle index for IR
Modified Activity Questionnaire
Nurses’ Health Study
- R2 :
Coefficient of determination
Type 2 diabetes
Tehran lipid and glucose study
United States Department of Agriculture
Visceral adiposity index
Shanbhogue VV, Finkelstein JS, Bouxsein ML, Yu EW. Metabolism: association between insulin resistance and bone structure in nondiabetic postmenopausal women. J Clin Endocrinol Metab. 2016;101:3114–22.
Brands MW, Manhiani MM. Sodium-retaining effect of insulin in diabetes. Am J Physiol-Regul Integr Compar Physiol. 2012;303:R1101–9.
Guo J, Breen DM, Pereira TJ, Dalvi PS, Zhang H, Mori Y, Ghanim H, Tumiati L, Fantus IG, Bendeck MP. The effect of insulin to decrease neointimal growth after arterial injury is endothelial nitric oxide synthase-dependent. Atherosclerosis. 2015;241:111–20.
Fornes R, Ormazabal P, Rosas C, Gabler F, Vantman D, Romero C, Vega MJ. Changes in the expression of insulin signaling pathway molecules in endometria from polycystic ovary syndrome women with or without hyperinsulinemia. Mol Med. 2010;16:129–36.
Shanik MH, Xu Y, Skrha J, Dankner R, Zick Y, Roth J. Insulin resistance and hyperinsulinemia: is hyperinsulinemia the cart or the horse? Diabetes Care. 2008;31(Suppl 2):S262-268.
Thomas DD, Corkey BE, Istfan NW, Apovian CM. Hyperinsulinemia: an early indicator of metabolic dysfunction. J Endocr Soc. 2019;3:1727–47.
Reaven GM. Role of insulin resistance in human disease. Diabetes. 1988;37:1595–607.
Kalan Farmanfarma K, Kaykhaei MA, Adineh HA, Mohammadi M, Dabiri S, Ansari-moghaddam A. Prevalence of metabolic syndrome in Iran: a meta-analysis of 69 studies. Diabetes Metab Syndr. 2019;13:792–9.
IDF Diabetes Atlas. https://www.diabetesatlas.org
Kelly GS. Insulin resistance: lifestyle and nutritional interventions. Altern Med Rev. 2000;5:109–32.
Jensen MK, Koh-Banerjee P, Franz M, Sampson L, Grønbæk M, Rimm EB. Whole grains, bran, and germ in relation to homocysteine and markers of glycemic control, lipids, and inflammation. Am J Clin Nutr. 2006;83:275–83.
Wu T, Willett WC, Hankinson SE, Giovannucci E. Caffeinated coffee, decaffeinated coffee, and caffeine in relation to plasma C-peptide levels, a marker of insulin secretion, in US women. Diabetes Care. 2005;28:1390–6.
Gower BA, Pearson K, Bush N, Shikany JM, Howard VJ, Cohen CW, Tison SE, Howard G, Judd S. Diet pattern may affect fasting insulin in a large sample of black and white adults. Eur J Clin Nutr. 2021;75:628–35.
Hu FB. Dietary pattern analysis: a new direction in nutritional epidemiology. Curr Opin Lipidol. 2002;13:3–9.
Jacques PF, Tucker KL. Are dietary patterns useful for understanding the role of diet in chronic disease? Oxford: Oxford University Press; 2001.
Dubé JJ, Fleishman K, Rousson V, Goodpaster BH, Amati F. Exercise dose and insulin sensitivity: relevance for diabetes prevention. Med Sci Sports Exerc. 2012;44:793.
García-Estévez DA, Araújo-Vilar D, Saavedra-González Á, Fiestras-Janeiro G, Cabezas-Cerrato J. Analysis of the relationship between body mass index, insulin resistance, and beta-cell function: a cross-sectional study using the minimal model. Metabolism. 2004;53:1462–6.
Tatsumi Y, Morimoto A, Miyamatsu N, Noda M, Ohno Y, Deura K. Effect of body mass index on insulin secretion or sensitivity and diabetes. Am J Prev Med. 2015;48:128–35.
Tabung FK, Wang W, Fung TT, Hu FB, Smith-Warner SA, Chavarro JE, Fuchs CS, Willett WC, Giovannucci EL. Development and validation of empirical indices to assess the insulinaemic potential of diet and lifestyle. Br J Nutr. 2016;116:1787–98.
Perez Rodrigo C, Aranceta J, Salvador G, Varela-Moreiras G. Food frequency questionnaires. Nutr Hosp. 2015;31(Suppl 3):49–56.
Brassard D, Lemieux S, Charest A, Lapointe A, Couture P, Labonté M-È, Lamarche B. Comparing interviewer-administered and web-based food frequency questionnaires to predict energy requirements in adults. Nutrients. 2018;10:1292.
Mazidi M, Katsiki N, Mikhailidis DP, Banach M. Effect of dietary insulinemia on all-cause and cause-specific mortality: results from a cohort study. J Am Coll Nutrit. 2020;39:407–13.
Azizi F, Ghanbarian A, Momenan AA, Hadaegh F, Mirmiran P, Hedayati M, Mehrabi Y, Zahedi-Asl S. Prevention of non-communicable disease in a population in nutrition transition: Tehran Lipid and Glucose Study phase II. Trials. 2009;10:5.
Esfahani FH, Asghari G, Mirmiran P, Azizi F. Reproducibility and relative validity of food group intake in a food frequency questionnaire developed for the Tehran Lipid and Glucose Study. J Epidemiol. 2010;20:150–8.
Asghari G, Rezazadeh A, Hosseini-Esfahani F, Mehrabi Y, Mirmiran P, Azizi F. Reliability, comparative validity and stability of dietary patterns derived from an FFQ in the Tehran Lipid and Glucose Study. Br J Nutr. 2012;108:1109–17.
Momenan AA, Delshad M, Sarbazi N, Rezaei GN, Ghanbarian A, Azizi F. Reliability and validity of the Modifiable Activity Questionnaire (MAQ) in an Iranian urban adult population. Arch Iran Med. 2012;15(5):279–82.
Sekaran U, Bougie R. Research methods for business, a skill building approach. Hoboken: Willey; 2003.
VanVoorhis CW, Morgan BL. Understanding power and rules of thumb for determining sample sizes. Tutor Quant Methods Psychol. 2007;3:43–50.
Tapsell LC, Neale EP, Satija A, Hu FB. Foods, nutrients, and dietary patterns: interconnections and implications for dietary guidelines. Adv Nutr (Bethesda, MD). 2016;7:445–54.
Babio N, Becerra-Tomás N, Martínez-González M, Corella D, Estruch R, Ros E, Sayón-Orea C, Fitó M, Serra-Majem L, Arós F, et al. Consumption of yogurt, low-fat milk, and other low-fat dairy products is associated with lower risk of metabolic syndrome incidence in an elderly mediterranean population. J Nutr. 2015;145:2308–16.
Rideout TC, Marinangeli CP, Martin H, Browne RW, Rempel CB. Consumption of low-fat dairy foods for 6 months improves insulin resistance without adversely affecting lipids or bodyweight in healthy adults: a randomized free-living cross-over study. Nutr J. 2013;12:56.
Yannakoulia M, Yiannakouris N, Melistas L, Kontogianni MD, Malagaris I, Mantzoros CS. A dietary pattern characterized by high consumption of whole-grain cereals and low-fat dairy products and low consumption of refined cereals is positively associated with plasma adiponectin levels in healthy women. Metabolism. 2008;57:824–30.
Battista F, Ermolao A, van Baak MA, Beaulieu K, Blundell JE, Busetto L, Carraça EV, Encantado J, Dicker D, Farpour-Lambert N, et al. Effect of exercise on cardiometabolic health of adults with overweight or obesity: Focus on blood pressure, insulin resistance, and intrahepatic fat—a systematic review and meta-analysis. Obes Rev. 2021;22(Suppl 4):e13269.
Hamjane N, Benyahya F, Nourouti NG, Mechita MB, Barakat A. Cardiovascular diseases and metabolic abnormalities associated with obesity: What is the role of inflammatory responses? A systematic review. Microvasc Res. 2020;131:104023.
Matson RI, Perry R, Hunt LP, Chong AH, Beynon R, Hamilton-Shield J, Birch L. Change in obesity-related metabolic abnormalities associated with body mass index improvement through lifestyle intervention: a meta-regression. Pediatr Diabetes. 2020;21:173–93.
Maximova K, Khodayari Moez E, Dabravolskaj J, Ferdinands AR, Dinu I, Lo Siou G, Al Rajabi A, Veugelers PJ. Co-consumption of vegetables and fruit, whole grains, and fiber reduces the cancer risk of red and processed meat in a large prospective cohort of adults from Alberta’s tomorrow project. Nutrients. 2020;12:2265.
Oh H, Lee HY, Jun DW, Lee SM. Low salt diet and insulin resistance. Clin Nutr Res. 2016;5:1–6.
Fu Z, Gilbert ER, Liu D. Regulation of insulin synthesis and secretion and pancreatic Beta-cell dysfunction in diabetes. Curr Diabetes Rev. 2013;9:25–53.
Fung TT, Hu FB, Schulze M, Pollak M, Wu T, Fuchs CS, Giovannucci E. A dietary pattern that is associated with C-peptide and risk of colorectal cancer in women. Cancer Causes Control CCC. 2012;23:959–65.
Wolever TM, Mehling C. Long-term effect of varying the source or amount of dietary carbohydrate on postprandial plasma glucose, insulin, triacylglycerol, and free fatty acid concentrations in subjects with impaired glucose tolerance. Am J Clin Nutr. 2003;77:612–21.
Sarafidis PA, Lasaridis AN, Nilsson PM, Pikilidou MI, Stafilas PC, Kanaki A, Kazakos K, Yovos J, Bakris GL. Validity and reproducibility of HOMA-IR, 1/HOMA-IR, QUICKI and McAuley’s indices in patients with hypertension and type II diabetes. J Hum Hypertens. 2007;21:709–16.
We claim our appreciation to the Tehran Lipid and Glucose Study participants for their generous support and the Research Institute for Endocrine Sciences, Tehran Lipid, and Glucose Study Unit staff for their precious help.
Faculty of Nutrition Sciences and Food Technology, National Nutrition and Food Technology Research Institute, Shahid Beheshti University of Medical Sciences, Tehran, Iran.
Ethics approval and consent to participate
All participants filled a written informed consent. This study was conducted in conformance with good clinical practice standards and was performed according to the Declaration of Helsinki 1975, as subsequent amendments. Also, this study was approved by the ethics research committee of the Research Institute for Endocrine Sciences, Shahid Beheshti University of Medical Sciences, Tehran, Iran (Approval ID: IR.SBMU.NNFTRI.REC.1399.055).
Consent for publication
The authors declared there is no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Table S1.
Food groups with their items used in developed indices. Table S2. Pearson and Spearman’s correlations coefficients between the insulinemic dietary and lifestyle patterns in the two examinations of TLGS. Table S3. Pearson and Spearman’s correlations coefficients between the insulinemic dietary and lifestyle patterns in the two examinations of TLGS. An alternative method for the calculation of indices. Calculation of unweighted forms of indices.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Mokhtari, E., Teymoori, F., Farhadnejad, H. et al. Development and validation of dietary and lifestyle insulinemic indices among Iranian adult population. Nutr Metab (Lond) 19, 5 (2022). https://doi.org/10.1186/s12986-021-00640-6
- Insulin resistance
- Dietary pattern