© 2002 by the American Diabetes Association, Inc.
Development of a Prediction Equation for Insulin Sensitivity From Anthropometry and Fasting Insulin in Prepubertal and Early Pubertal ChildrenFrom the Institute for Health Promotion and Disease Prevention Research, Departments of Preventive Medicine and Physiology & Biophysics, Keck School of Medicine, University of Southern California, Los Angeles, California
OBJECTIVETo test the utility of homeostasis model assessment (HOMA) in predicting insulin sensitivity [x10- 4 min-1/(µIU/ml)] in children and to develop and compare two new prediction equations for insulin sensitivity in children using demographic and anthropometric measures in the presence or absence of fasting insulin. RESEARCH DESIGN AND METHODSWe studied 156 white and African-American children with complete data (mean age 9.7 ± 1.8 years, 87.8% Tanner Stage 1 or 2). For development of new equations, two-thirds of the children were randomly assigned to a development group, whereas the remaining children were assigned to a cross-validation group. RESULTSA modified HOMA equation accurately predicted insulin sensitivity, but its utility is similar to fasting insulin alone. Demographic and anthropometric measures alone did not predict insulin sensitivity accurately, even when precise measures of body composition were included in the prediction model. Ethnicity, calf skinfold, and fasting insulin together explained 73% of the variance in insulin sensitivity and accurately predicted insulin sensitivity. The regression of measured versus predicted insulin sensitivity in the cross-validation group was not significantly different from the line of identity (P > 0.05). Mean difference between measured and predicted insulin sensitivity was also not significant (P > 0.05). Some bias was apparent, particularly in white boys. CONCLUSIONSEthnicity, calf skinfold, and fasting insulin can accurately predict insulin sensitivity with greater precision than HOMA or fasting insulin alone (R2 = 0.73). Future studies, however, are needed to examine whether a universal equation is possible. A cross-validated prediction equation may be useful in population-based studies when complex measures of insulin sensitivity are not available.
Abbreviations: CV, coefficient of variation DEXA, dual X-ray absorptiometry FSIGTT, frequently sampled intravenous glucose tolerance test HOMA, homeostasis model assessment
Insulin sensitivity is measured by analyzing the degree of glucose uptake at induced hyperinsulinemia (1). However, in many clinical and research settings, it cannot always be easily and directly obtained due to the high cost and complexity of precise techniques. As a result, the use of available proxy measures, such as measures of fasting insulin and glucose, can be very beneficial. Because no single proxy measure can account for all the variance in insulin sensitivity, a combination of fasting blood measures with anthropometric and/or demographic measures may be appropriate. However, few studies to date, and none in children, have cross-validated any prediction equations for insulin sensitivity (2 4). Most studies have only examined the associations of adiposity, fasting glucose, and fasting insulin with insulin sensitivity in a correlative fashion. Therefore, development of a validated prediction equation may be useful for future research and clinical needs. The homeostasis model assessment (HOMA), which is based on fasting insulin and glucose, has been used frequently to predict insulin resistance (5). Although the HOMA method has been used in children to predict insulin resistance (68), its accuracy and precision have never been examined relative to more sophisticated measures of insulin sensitivity in children, such as the clamp or intravenous glucose tolerance test. Given the fact that direct and precise measures of insulin sensitivity are often difficult to obtain, especially in children, the current study first aimed to examine whether the HOMA method could accurately predict insulin sensitivity in our sample of children. In addition, to estimate insulin sensitivity with more accuracy and precision, we aimed to develop and compare two prediction equations of insulin sensitivity in children. The first equation included demographic and anthropometric variables, whereas the second included the same variables as well as fasting insulin and glucose. Demographic measures included ethnicity, gender, age, and family history of type 2 diabetes. Anthropometric measures included weight, height, skinfold thicknesses, waist circumference, and sexual maturity. The significance of this study is that a cross-validated prediction equation for insulin sensitivity may be useful, especially in population-based studies, when complex and expensive techniques are not available.
Subjects A total of 156 children were recruited by newspaper and radio advertisements and by word of mouth. Subjects were screened by medical history and were ineligible if they were 1) <4 years of age; 2) taking medications known to affect body composition or physical activity (e.g., prednisone, Ritalin, or growth hormone); 3) previously diagnosed with syndromes known to affect body composition or fat distribution (e.g., Cushings syndrome, Downs syndrome, type 1 diabetes, or hypothyroidism); or 4) diagnosed previously with any major illness. Because the intent was to recruit a heterogeneous group of children, there were no criteria for other characteristics such as obesity. This study was approved by the Institutional Review Board at the University of Alabama at Birmingham. Parents provided informed consent before testing began.
Protocol and measurements of anthropometry An evening snack was allowed but only water and energy-free, noncaffeinated beverages were permitted after 2000. On the following morning after an overnight fast, blood was collected and a tolbutamide-modified frequently sampled intravenous glucose tolerance test (FSIGTT) was performed.
Tolbutamide-modified FSIGTT
Assay of glucose and insulin Insulin was assayed in duplicate 200-µl aliquots with Coat-A-Count kits (Diagnostic Products, Los Angeles, CA). According to the supplier, cross-reactivity of this assay with proinsulin is 40% at midcurve; C-peptide was not detected. In our laboratory, this assay had a sensitivity of 11.4 pmol/l (1.9 µIU/ml), a mean intra-assay CV of 5%, and a mean interassay CV of 6%. Commercial quality control sera of low, medium, and high insulin concentration (Lyphochek; Bio-Rad, Anaheim, CA) were included in every assay to monitor variation over time.
Assessment of insulin sensitivity by HOMA
Statistics Two-thirds of the sample (n = 104) was randomly assigned to a development group from which the new equations were derived. The remaining one-third of the group (n = 52) was used for cross validation. Prediction equations were considered to cross-validate if the regression between measured and predicted insulin sensitivity in the cross-validation group was not statistically different from the line of identity (i.e., slope = 1, intercept = 0), and the paired Students t test was not significantly different from 0. A plot of (measured-predicted) versus measured insulin sensitivity was used to examine possible bias in the prediction equations. There was one observation per subject in the current analysis, and all subjects provided complete data. In all analyses, a type I error of 0.05 was used. All procedures were conducted using SAS software (Version 8.01; SAS Institute, Cary, NC).
The descriptive characteristics of the sample with regard to ethnicity, gender, age, Tanner Stage, weight, height, waist circumference, calf skinfold thickness, suprailiac skinfold thickness, fasting glucose, fasting insulin, and insulin sensitivity are shown in Table 2. The development group and the cross-validation group did not differ on any of these variables (P > 0.05). Because the current sample is part of a larger, heterogeneous observational study on childhood obesity, the range of insulin sensitivity reflects that heterogeneity. The range of insulin sensitivity in the current study is consistent with what we have published previously. Although insulin sensitivity obtained from the minimal model cannot be compared directly with that obtained from the clamp technique, the two methods have been shown to correlate as highly as 0.89 (13). Therefore, they both give valid assessments of overall insulin sensitivity (13).
Fitting HOMA in children Regression of FSIGTT-measured insulin sensitivity on HOMA-estimated insulin resistance in the development group yielded an equation where loge insulin sensitivity = 2.393 - (0.306 * HOMA insulin resistance). Insulin sensitivity predicted by HOMA accounted for 63.4% of the variance in observed insulin sensitivity. This equation was then subsequently used to predict insulin sensitivity in the cross-validation group. Regression of measured insulin sensitivity against insulin sensitivity predicted by the equation containing HOMA in the cross-validation group showed that there was not a significant deviation from the line of identity (intercept ± SE = -0.27 ± 0.21, P > 0.05; slope ± SE = 1.19 ± 1.14, P > 0.05; see Fig. 1). The mean difference between measured and predicted insulin sensitivity was also not significantly different from 0 (untransformed mean ± SD = 5.4 ± 3.8 vs. 4.6 ± 1.7 x 10-4 min-1/(µIU/ml), P > 0.05). In addition, bias was detected by regressing the difference between FSIGTT-measured insulin sensitivity and insulin sensitivity predicted by the equation containing HOMA on FSIGTT-measured insulin sensitivity (intercept ± SEM = -0.70 ± 0.09, P < 0.001; slope ± SEM = 0.50 ± 0.06, P < 0.001; r = 0.77, P < 0.001; see Fig. 2).
Development of prediction equation with demographic and anthropometric measures For children randomly assigned to the development group, the potential demographic and anthropometric predictors of insulin sensitivity were entered into a stepwise prediction model, and an equation including calf skinfold thickness (P < 0.001), ethnicity (P < 0.001), weight (P < 0.05), and gender (P < 0.05) was defined (Table 3). In this equation, calf skinfold thickness accounted for 42.5% of the variance in insulin sensitivity and ethnicity accounted for 13.5% of the variance, whereas weight and gender accounted for 2.5 and 1.8% of the variance, respectively (total R2 = 0.60).
Cross-validation of the equation was performed in the randomly assigned validation group. Through regression analysis of the measured versus predicted insulin sensitivity, the equation with only demographic and anthropometric measures was significantly different from the line of identity (intercept ± SEM = 0.64 ± 0.15, P < 0.001; slope ± SEM = 0.83 ± 0.14, P > 0.05; see Fig. 3). By paired Students t tests, there was a significant mean difference between measured and predicted insulin sensitivity, using the equation with only demographic and anthropometric measures (untransformed mean ± SD = 5.4 ± 3.9 vs. 3.0 ± 1.5 x 10-4 min-1/(µ IU/ml), P < 0.001).
Development of prediction equation with demographic, anthropometric, and fasting blood measures A second stepwise regression model was performed using potential demographic, anthropometric, as well as fasting insulin and glucose measures. This analysis yielded an equation with fasting insulin (P < 0.001), ethnicity (P < 0.001), and calf skinfold thickness (P < 0.001) as the significant predictors (Table 4). In this equation, fasting insulin accounted for 63.8% of the variance in insulin sensitivity, ethnicity accounted for 5.1% of the variance, and calf skinfold thickness accounted for 3.7% of the variance (total R2 = 0.73).
Cross-validation of this equation showed that the regression of measured versus predicted insulin sensitivity was not significantly different from the line of identity (intercept ± SEM = 0.07 ± 0.15; slope ± SEM = 0.97 ± 0.09, P > 0.05; see Fig. 4). In addition, paired Students t tests showed no significant difference between measured and predicted insulin sensitivity (untransformed mean ± SD = 5.4 ± 3.9 vs. 4.8 ± 2.4 x 10-4 min-1/ (µIU/ml), P > 0.05). The equation including fasting insulin was therefore successfully cross-validated. Table 5 shows the means of measured and predicted insulin sensitivity by gender and ethnicity.
An examination of the discrepancy between measured and predicted insulin sensitivity as a function of measured insulin sensitivity in the cross-validation group is shown in Fig. 5 (intercept ± SEM = -0.40 ± 0.11, P < 0.001; slope ± SEM = 0.30 ± 0.07, P < 0.001; r = 0.53, P < 0.001). In this case, an intercept of 0 and a slope of 0 would indicate absence of bias. Because some bias was apparent in the group as a whole (i.e., significant slope), we conducted a stratified analysis by gender and ethnicity. We found that the bias occurred largely in white boys (intercept ± SEM = -1.24 ± 0.18, P < 0.001; slope ± SEM = 0.73 ± 0.09, P < 0.001) and, to a lesser extent, in African-American girls (intercept ± SEM = -0.48 ± 0.20, P < 0.05; slope ± SEM = 0.49 ± 0.16, P < 0.01). Bias was not detected in either white girls (intercept ± SEM = -0.26 ± 0.17, P > 0.05; slope ± SEM = 0.13 ± 0.10, P > 0.05) or African-American boys (intercept ± SEM = -0.04 ± 0.24, P > 0.05; slope ± SEM = 0.19 ± 0.18, P > 0.05).
The current study is the first in children to examine the possibility of predicting insulin sensitivity (as based on an FSIGTT and the minimal model) by using simple measures of demography and anthropometry, in the presence or absence of fasting insulin. The main findings of this paper are: 1) the modified HOMA method accurately predicted insulin sensitivity in this sample of children but did not account for more of the variance in observed insulin sensitivity than fasting insulin alone; 2) insulin sensitivity in children could not be accurately predicted by demographic and anthropometric measures alone; 3) insulin sensitivity could be accurately predicted in children by an equation with demographic, anthropometric, and fasting insulin measures, and this equation was better than HOMA or fasting insulin alone; and 4) although the combination of demographic, anthropometric, and fasting insulin measures accurately predicted insulin sensitivity in the group as a whole, the equation lacked individual precision in white boys and, to a lesser extent, in African-American girls.
In the first part of the study, we attempted to cross-validate an equation containing the HOMA index in a cohort of children. We found that, in children, insulin sensitivity could be validly predicted using this method. However, similar to findings from previous studies in adults, the modified HOMA equation may not be very precise and stable in children, in whom predicted values of insulin sensitivity from the equation containing HOMA accounted for only
The current study was not able to cross-validate a prediction equation of insulin sensitivity, using basic demographic and anthropometric measures. We decided to use only conventional measures of adiposity because our goal was to develop a prediction equation of insulin sensitivity that was easily accessible. Use of BMI as a potential predictor variable, rather than simple weight and height measurements, did not yield different results. However, it is possible that more precise and complex measures of adiposity, such as dual energy X-ray absorptiometry (DEXA), would better correlate with insulin sensitivity. Nevertheless, in a separate stepwise regression analysis in a subset of children, we found that the inclusion of total body fat and lean tissue mass measured by DEXA did not improve the prediction equation significantly (
An equation with demographic, anthropometric, and fasting insulin measures validly predicted insulin sensitivity in our sample. Comparing the equation with demographic and anthropometric measures alone and the one with demographic, anthropometric, as well as fasting insulin measures, the latter accounted for 13% more of the variance in insulin sensitivity. Fasting insulin was shown to be a major predictor of insulin sensitivity, accounting for It is also noteworthy that because fasting insulin alone accounted for most of the variance in observed insulin sensitivity, it did yield a univariate equation (Step 1, Table 4) that was successfully cross-validated in the validation group (regression results not shown). However, fasting insulin alone was not as accurate as the combination of fasting insulin with demographic and anthropometric measures (Step 3, Table 4). In addition, to rely on fasting insulin alone, the correlation between (measured-predicted) and measured insulin sensitivity was 0.75, suggesting that there was substantially more bias than the new equation developed in this study (i.e., recall that a correlation of 0 represents no bias in this instance; in other words, at any given level of measured insulin sensitivity, the difference between measured and predicted insulin sensitivity is 0). The additional measures provided more accuracy and slightly more precision in the prediction by accounting for an additional 9% of the variance in observed insulin sensitivity. These measures are simple and inexpensive, therefore we believe they should be included. We also note that fasting glucose was not a significant predictor of insulin sensitivity. This is not too surprising because fasting glucose does not vary greatly in healthy children. This may partly explain why HOMA is only as useful as fasting insulin alone in children, because it estimates insulin resistance by relying solely on the combination of fasting insulin and fasting glucose. Although our equation with demographic, anthropometric, and fasting insulin measures was successfully cross-validated in the current sample, there seemed to be some bias in the prediction. In a separate stratified analysis by gender and ethnicity, where the deviation of measured and predicted insulin sensitivity was regressed against measured insulin sensitivity, we found that the bias occurred largely in white boys and, to a lesser extent, in African-American girls. Bias was not detected in either white girls or African-American boys. It is not clear why this may be the case, particularly given the fact that the current sample was limited by its size for us to draw any conclusions per any gender by racial subgroup. However, compared with white girls and African-American boys combined, white boys in our sample had significantly higher insulin sensitivity and African-American girls had significantly lower insulin sensitivity (mean 5.4, 7.8, 3.5 x 10-4 min-1/(µIU/ ml), respectively, P < 0.001). Therefore, it may be that at higher and lower ranges of insulin sensitivity, the current prediction equation loses some degree of precision. Future studies are needed, however, to examine whether this remains a problem in larger samples. There are some limitations of the study that should be considered. First, because only African-American and white children were included in the current study, generalization to other ethnic groups is not possible. Second, our sample size did not allow us to more carefully examine whether it would be useful to develop a separate equation for each gender by ethnicity subgroup or for different ranges of insulin sensitivity. This warrants further investigation in the future, when larger samples of children are available. Third, the fact that Tanner Stage was not selected as a significant predictor may be due to the fact that only 12% of our sample was far enough into puberty. Therefore, in prepubertal and early pubertal children, fat composition and ethnicity (other than fasting insulin) may be most important in predicting insulin sensitivity. Given that we have previously shown, in a longitudinal study, that even in the most obese children insulin sensitivity falls by 33% from Tanner I to Tanner III (14), we recognize that our equations may not be generalizable to all stages of maturation, and additional analysis will be required to incorporate Tanner stage into any future equations. Finally, readings of insulin concentrations may be different, depending on the type of assay used. In the current study, polyclonal antibodies were used in the insulin assay. However, as long as the discrepancy from different antibodies is consistent across individuals, we believe that our equation would still be valid using other insulin assays. Previously, we had compared the proinsulin levels between African-American and Caucasian children and found that there was no difference (Gower BA, Goran MI, unpublished data). Our equation could make a difference if African-Americans had more proinsulin than Caucasians (and therefore had higher nonspecific insulin levels but not higher specific insulin levels). However, this was not the case; in fact, the data were identical regardless of the use of a specific insulin assay or a nonspecific one. Future studies in different laboratories should determine whether our equations are equally valid using other forms of assay. In conclusion, insulin sensitivity in children cannot be validly predicted by the combination of merely demographic and anthropometric measures. An equation containing HOMA can validly predict insulin sensitivity, but its precision is not better than an equation containing fasting insulin alone. A combination of demographic (i.e., ethnicity, gender), anthropometric (i.e., calf skinfold, suprailiac skinfold, weight), and fasting insulin measures predicted insulin sensitivity in white and African-American children better than the equation with HOMA or fasting insulin alone. Our validated equation of insulin sensitivity may be useful in population-based studies when complex techniques of measuring insulin sensitivity are not available. Because some bias exists in the current equation, its use in diagnosing insulin sensitivity on an individual basis is not encouraged. Future studies need to further examine whether a universal equation is feasible. Nevertheless, the current study is the first of such in children and suggests that a simple and easily accessible prediction equation of insulin sensitivity in children is possible.
This study was supported by National Institutes of Health Grant R01-HD/HL-33064 and the University of Alabama at Birmingham General Clinical Research Center (RR-00032).
Address correspondence and reprint requests to Michael I. Goran, PhD, Institute for Prevention Research, University of Southern California, 1540 Alcazar St, CHP208D, Los Angeles, CA 90089. E-mail: goran{at}usc.edu. Received for publication 18 September 2001 and accepted in revised form 22 May 2002. A table elsewhere in this issue shows conventional and Système International (SI) units and conversion factors for many substances.
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||