|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
INNOVATIVE METHODOLOGY
Facultades de Farmacia y Medicina, Universidad CEU (Centro de Estudios Universitarios) San Pablo, Madrid, Spain
Submitted 4 February 2008 ; accepted in final form 15 September 2008
| ABSTRACT |
|---|
|
|
|---|
homeostasis model assessment of insulin resistance; quantitative insulin sensitivity check index; fasting glucose-to-insulin ratio; hyperinsulinemic isoglycemic clamp; calibration model
In the past, surrogate measures of insulin resistance have been developed based on measurements of fasting glucose and insulin concentrations (15, 21, 26). These indexes have been validated in humans by comparison with the hyperinsulinemic-euglycemic clamp, and they were found to correlate reasonably well with whole body insulin sensitivity determined with the clamp technique. The homeostasis model assessment of insulin resistance (HOMA-IR) was first described by Matthews et al. (26) as a measure of basal insulin sensitivity. The fasting glucose-to-insulin ratio (FGIR) has become popular since its first description as an accurate index of insulin sensitivity in women with polycystic ovary syndrome (21). The most recently proposed derivation, using simple fasting measures, is the quantitative insulin sensitivity check index (QUICKI), which is based on a log transform of the insulin glucose product (15). The log HOMA-IR and QUICKI are simply related by inversion and differ only by the normalizing constant used to calculate the HOMA-IR. To date, the best direct validation studies of simple surrogate indexes of insulin sensitivity, including HOMA-IR and QUICKI, were based on examining correlations with the reference glucose clamp method (15, 18, 25, 38). They provide a simple estimate for whole body insulin sensitivity with variability and discriminant power comparable to those of the euglycemic-hyperinsulinemic clamp (25), the minimal model (37), or the OGTT (16). Furthermore, QUICKI and logHOMA have been found to be excellent measures to predict the insulin sensitivity index obtained in the clamp (SIClamp) (10).
In animal studies, including different rat models, one or several of these indexes have been applied to quantify insulin sensitivity (19, 28, 34–36). Although a recent paper has analyzed the correlation between surrogate indexes of insulin resistance in mice (20), no study has been designed so far to validate these indexes in rats. This lack of information, therefore, violates the assumptions of the model (38). Furthermore, to our knowledge, no analysis of the discriminant power of these indexes has been performed in animals. The purpose of the present study was to evaluate whether the HOMA-IR, QUICKI, and FGIR indexes can be used to accurately estimate insulin sensitivity in nonpregnant and pregnant rats. We have performed this study in Wistar and Sprague-Dawley rats since they are the most commonly used strains for insulin-stimulated glucose measurements, and they are known to exhibit variations in whole body or tissular insulin sensitivity (22).
| MATERIALS AND METHODS |
|---|
|
|
|---|
Hyperinsulinemic isoglycemic clamp studies.
To assess insulin sensitivity, nonpregnant rats and rats at day 20 of gestation were subjected to a hyperinsulinemic isoglycemic clamp after fasting for 6 h. In brief, blood samples were obtained from the tail tips for determination of basal glucose and insulin levels. Subsequently, animals were anesthetized with ketamine cocktail anesthesia (50 mg/ml ketamine, 5 mg/ml diazepam, and 1 mg/ml atropine; 5:4:1 vol/vol/vol at a final volume of 1 ml/ kg), and a Silastic brand catheter (0.02 in. ID, 0.037 in. OD; Dow Corning, Midland, MI) was placed in the right jugular vein and another catheter in the right femoral vein. Each catheter was connected to an infusion pump (Precidor Infusion Pump Type 5003; Infors HT, Denkendorf, Germany). Human insulin (Actrapid monocomponent; Novo, Copenhagen, Denmark) was infused by means of one pump at a constant rate of 16 µl/min (0.8 IU·h–1·kg body wt–1) for 60 min. Glucose (20%) was infused at a variable rate through the second pump to maintain the blood glucose concentration constant at basal levels. Blood samples were collected from the tail tip at different time points to monitor the glycemia of the animals. Steady-state glucose infusion was generally achieved within 30 min after starting the clamp experiment. Additional blood samples (200 µl) were collected to determine the glucose and insulin concentration at the steady state in EDTA-plasma samples. The glucose disposal rate (M) was estimated as the rate of glucose infusion at the steady state normalized to body weight. SIClamp, as proposed by Ader and Bergman (1), was calculated as M/(G x
I), where G is the steady-state blood glucose concentration, and
I is the increment of insulin concentration from basal levels to steady state.
Plasma analysis. Blood glucose during the clamp was measured by an immobilized glucose oxidase method (Reflolux IIM; Boehringer-Mannheim) (4). In EDTA-treated plasma samples, glucose was determined by an enzymatic colorimetric test (GOD-PAP; Roche Diagnostics, Barcelona, Spain) and insulin (standard curve range 0.1–10 µg/l; interassay: 8.5–9.4%; intra-assay: 1.4–4.6%; rat C-peptide not detectable) by a specific RIA kit for rats (Linco).
Calculation of insulin sensitivity indexes. From the short-term fasting plasma glucose (FPG) and insulin (FPI) values obtained in each animal before the clamp, the following indexes were calculated as estimates of insulin sensitivity: HOMA-IR, QUICKI, and FGIR.
HOMA-IR was calculated as the product of the FPG and FPI levels, divided by a constant, assuming that control young adult rats have an average HOMA-IR of 1, analogous to the assumptions applied in the development of HOMA-IR in humans (26). The equation was as follows HOMA-IR = (FPG x FPI)/2,430, where FPI was in microunits per milliliter and FPG in milligram per deciliter. QUICKI was calculated according to the original formula (15) as the inverse log sum of fasting insulin in microunit per milliliter and fasting glucose in milligram per deciliter. QUICKI = 1/[log(FPG) + log(FPI)]. Finally, FGIR was calculated as the ratio of FPG divided by FPI levels (21). FGIR = FPG/FPI, where FPG was in milligrams per deciliter and FPI in microunits per milliliter.
Calibration model analysis of surrogate insulin sensitivity indexes.
To evaluate the ability of surrogate indexes to accurately predict insulin sensitivity as determined by the reference glucose clamp method, we used a calibration model to compare the ability of HOMA-IR, QUICKI, and FGIR to predict SIClamp as previously described by others (10, 20). In brief, calibration is the inverse of regression (7), thus using an estimated model y = f (x;
), where x is the independent variable, y is the dependent variable, and
is an unknown parameter, predicting a new y* for a given x* is regression. Conversely, predicting a new x* for a given y* is calibration. Accordingly, in the present study, we fitted a calibration model xi =
+ βyi +
i, where xi is the SIClamp, yi is each surrogate index, and
i is the random error for the ith subject. Even though SIClamp is measured with error, the assumption can be made that the measurement error of SIClamp (determined from a direct and data-intensive protocol) is very small relative to that of the indexes determined from single fasting measurements. Therefore, to simplify the analysis, we neglected the measurement error for SIClamp in our calibration model. For each surrogate index, two types of predicted residuals were considered. The first one is derived from the calibration model with all animals included and represents the difference between the measured SIClamp (xi, for the i subject) and the fitted SIClamp for the same i subject. The second one is the residual obtained from a leave-one-out cross-validation model and represents the difference between the measured SIClamp (xi, for the i subject) and the predicted SIClamp from the calibration model that excludes the i subject. Next, predictive accuracy was evaluated by root mean squared error of prediction (RMSE) and leave-one-out cross-validation type root mean squared error of prediction (CVPE). Smaller values of RMSE and CVPE indicate better prediction. The distribution of the obtained residuals for each index was displayed in box-and-whisker plots that ended at the first and third quartiles (27).
Data analysis and statistical evaluation. Results are expressed as means ± SE of 12–20 animals/group. Regarding the lognormal distribution of the insulin concentration, the statistical analyses were applied to the natural logarithm (log) of this parameter. All variables were evaluated for normality of distribution with the Kolmogorov-Smirnov goodness of fit. Where indicated in Tables 1–5 and the legends for Figs. 1–4, statistical comparisons between two groups were made with the Student's t-test. The level for statistical significance was set at 0.05 (P < 0.05). The relationship between insulin sensitivity, determined during the clamp SIClamp and indexes obtained from fasting glucose and insulin, was based on correlation analysis (Pearson coefficient) between pairs of indexes. Assessment of the performance of the various models was made using the receiver operating characteristic (ROC) curves by plotting the sensitivity against the corresponding false-positive rate (100-specificity) (14). The area under the ROC curve (AUC) was used as a measure of how well a continuous variable predicts the development of insulin resistance. A test with perfect discrimination power yields a ROC curve that passes through the upper left corner with an AUC of one (100% sensitivity and 100% specificity). Thus the closer the ROC area to one, the higher the discriminant power of the method. To construct the ROC curves, the presence of insulin resistance was defined according to the World Health Organization (European Group Insulin Resistance) (5) as a SIClamp value below the 25th percentile of the normal distribution in the nonpregnant animals (normal insulin sensitivity). To establish potential cutoff values for HOMA-IR, QUICKI, and FGIR, we determined the optimal decision point from the ROC curve, assigning equal weights to the sensitivity and specificity of the test. Statistical comparison of the areas under the ROC curves, derived from the same set of animals, was performed as described by Hanley and McNeil (13), taking into account the correlation between the areas that is induced by the paired nature of the data. Pearson correlation coefficients and ROC analysis were calculated using GraphPad programs (version 5.0 for Macintosh). Statistical comparisons of areas under ROC curves and Pearson correlation coefficients from the same sample were made using the SimpleStat software. Calibration model and leave-one-out cross-validation analysis were performed by MATLAB version 7.
|
|
|
| RESULTS |
|---|
|
|
|---|
Validation studies of insulin sensitivity indexes. Because the direct measurement of insulin sensitivity with the glucose clamp is complex, we evaluated whether indirect estimates based on easy-to-measure fasting glucose and insulin levels can be used as valid predictors of insulin sensitivity in rats. For this purpose, first, we analyzed whether the indexes obtained from fasting glucose and insulin levels correlate with each other. Both in nonpregnant and late-pregnant Wistar rats, all indexes correlated significantly with each other. The analysis reveals a higher degree of correlation between QUICKI and HOMA-IR (r = –0.970 and –0.997, P < 0.001 for nonpregnant and pregnant rats, respectively) than between these latter indexes and FGIR (r = –0.653 and 0.781 for correlations of HOMA-IR to FGIR and of QUICKI to FGIR, respectively, in nonpregnant rats; r = –0.558 and 0.568 for correlations of HOMA-IR to FGIR and of QUICKI to FGIR, respectively, in late-pregnant rats). Correlation analysis including all Wistar animals shows the same pattern. Similar results were obtained with Sprague-Dawley rats, obtaining Pearson correlation coefficients between HOMA-IR and QUICKI close to 1 (–0.958 and –0.972 for nonpregnant and late-pregnant rats, respectively), and significantly lower (P < 0.001) between these indexes and FGIR (r = –0.669 and 0.676 for correlations of HOMA-IR to FGIR and of QUICKI to FGIR, respectively, in nonpregnant rats; r = –0.705 and 0.673 for correlations of HOMA-IR to FGIR and of QUICKI to FGIR, respectively, in late-pregnant rats). When the analysis was performed with the whole group of Sprague-Dawley animals (nonpregnant and late-pregnant rats), the results were very similar to those obtained with the Wistar strain.
The utility of these fasting indexes in estimating insulin resistance depends on the underlying correlation of these estimates with directly determined experimental data, being the euglycemic-hyperinsulinemic clamp the gold standard for quantifying insulin resistance. Table 2 shows the relationship between SIClamp and the HOMA-IR, QUICKI, and FGIR indexes, as estimated from plasma glucose and insulin levels, both in Wistar and Sprague-Dawley rats. First, the correlation analysis was performed separately with the nonpregnant and the late-pregnant group. The relationship (Pearson correlation coefficient) between SIClamp and the three indexes was statistically significant (P < 0.05 for all comparisons), independent of whether the rats were pregnant or not (Table 2). Because the slopes for the correlations of nonpregnant or pregnant rats in each strain were not significantly different, all data were pooled. When the entire group of animals, i.e., both nonpregnant and pregnant rats, was included in the analysis, the association of each of the indexes with SIClamp was even higher. QUICKI was the index that showed the best correlation to SIClamp in nonpregnant and late-pregnant rats in both strains of animals (nonpregnant: 0.745 and 0.869, late pregnant: 0.727 and 0.725, for Wistar and Sprague-Dawley rats, respectively; P < 0.05). In general, the weakest relationship was found for FGIR vs. SIClamp (nonpregnant: 0.497 and 0.614, late pregnant: 0.705 and 0.556, for Wistar and Sprague-Dawley rats, respectively; P < 0.05 compared with the correlation obtained between QUICKI and HOMA-IR vs. SIClamp).
|
Correlation may be excellent even when prediction of reference values by the surrogate is poor. Thus it is important to evaluate the ability of surrogate indexes to accurately predict insulin sensitivity as determined by the reference glucose clamp method. In the present study, we used a calibration model to compare the ability of HOMA-IR, QUICKI, and FGIR to predict SIClamp. To this end, we regressed measured SIClamp for each animal on each surrogate index and fitted these data to a calibration model. Next, the predictions of SIClamp obtained from the calibration model (using the leave-one-out cross-validation approach) were plotted as a function of measured SIClamp by each index (Fig. 2). When a surrogate index perfectly predicts SIClamp, the values fall on a straight line with a slope of one and a y-intercept of zero. As shown in Fig. 2, the three indexes generated accurate predictions of SIClamp (with slopes of 1.02 ± 0.13, 1.02 ± 0.12, and 1.01 ± 0.12; intercepts of –0.19 ± 1.06, –0.10 ± 0.98, and –0.08 ± 0.95 for fitting between SIClamp vs. SIClamp predicted by HOMA-IR, QUICKI, and FGIR, respectively). In fact, statistical analysis indicates that the data did not differ significantly from a straight line. In addition, a linear least-squares fit between predicted SIClamp and measured SIClamp derived from the different indexes yielded correlation coefficients (r = 0.692, 0.720, and 0.728, for predictions by HOMA-IR, QUICKI, and FGIR, respectively, P < 0.0001 for all analysis) that did not significantly differ between each other.
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
First, we performed a correlation study of theses indexes with SIClamp. The HOMA-IR index is based on the premise that circulating glucose and insulin levels are determined by a feedback loop between the liver and the pancreas; thus, this index essentially reflects on changes in hepatic insulin sensitivity. This model has been used for many years and has been validated in different physiological and pathological conditions in humans showing a very good correlation with the clamp method (6, 38). This index has also been used in some animal studies (19, 28, 34–36). However, to our knowledge, no validation with the clamp has been made in rats, and only two studies have been performed to validate HOMA-IR in animal models. The first study was conducted in cats, comparing glucose and insulin-based indexes with the minimal model. In this study, the authors found that the most useful predictors were basal plasma insulin and HOMA-IR (2). In a very recent report, performed in the mouse, it has been found that QUICKI and HOMA-IR were modestly correlated with SIClamp. The authors pointed out that this may be due to inherent technical difficulties in performing clamps in mice (20).
A criticism of the HOMA model is its deviation from linearity with increasing insulin resistance in human pathologies such as gestational diabetes (18) or type 2 diabetes (6). However, in our study performed with rats, similar correlations have been obtained for the overall insulin sensitivity derived from HOMA-IR or log HOMA-IR with the SIClamp, and these results are similar to those obtained in human pregnancy and in subjects with different degrees of insulin sensitivity (6). Therefore, despite its possible limitations in assessing peripheral insulin sensitivity, HOMA-IR is a good predictor of total body insulin sensitivity during pregnancy both in humans (18) and in rats (this study). In addition, our findings support the notion that HOMA-IR may be an useful tool to assess maternal insulin status independent of the rat strain used.
In the present study, insulin sensitivity assessed with the QUICKI index showed the strongest correlation with direct measurements of insulin sensitivity using the glucose clamp, being as significant as those reported previously in human studies (3, 15). The strength of the relation was maintained when we examined the data in nonpregnant and late-pregnant animals. Similar results have been observed in human pregnancy, where QUICKI has been found to be a good estimate of insulin sensitivity during both early and late pregnancy (18). Thus QUICKI provides an excellent alternative to the rather laborious and complex glucose clamps for assessing insulin sensitivity in rats. Although a comparison of insulin sensitivity between Sprague-Dawley and Wistar rats was beyond the scope of our study, it should be noted that the fasting indexes obtained here are in agreement with the previously reported lower insulin sensitivity in Sprague-Dawley compared with Wistar rats (31).
FGIR became a popular index (11) since its first application in women with polycystic ovary syndrome (21), but, as yet, this index has not been validated in animal models. In this study, although FGIR also correlated significantly with SIClamp, the correlation was weaker than the correlation found for HOMA-IR or QUICKI, both in nonpregnant and pregnant animals. In fact, it has been shown that FGIR is a conceptually flawed index of insulin sensitivity (29), since it does not appropriately reflect the physiology underlying insulin sensitivity, in particular when fasting glucose levels are not in the normal range. As shown above and according to the normal physiology of pregnancy, glucose is lower in late-pregnant rats than in the nonpregnant animals, which could artificially decrease the FGIR in the former group.
Thus HOMA-IR, QUICKI, and FGIR correlated significantly with each other, and, although each of these indexes correlated significantly with the SIClamp, statistical analysis of the correlation coefficients showed that the QUICKI provided the stronger correlations followed by HOMA-IR and FGIR.
To evaluate the predictive accuracy of these indexes, we used a calibration model, obtaining two criterion functions, the RMSE and a CVPE. CVPE is more robust than RMSE because it uses an estimate that excludes the ith subject when predicting results for the ith subject. CVPE also handles extreme data in a more rigorous way. In our study, both RMSE and CVPE were similar, suggesting that there were no extreme outliers that could have introduced a bias into the obtained results. Furthermore, we did not detect any differences of theses parameters comparing the different surrogate indexes. This suggests that, in rats, HOMA-IR, QUICKI, and FGIR provide similar accuracy in predicting SIClamp. Consistent with this finding, the distribution of residuals was very similar, showing only one outlier for HOMA-IR. In summary, this calibration model corroborates that, both in Wistar and Sprague-Dawley rats, the three indexes HOMA-IR, QUICKI, and FGIR provide comparable accuracy in predicting SIClamp.
Despite their use in animal models, the predictive performance of these indexes to identify insulin-resistant animals has not been examined so far. Therefore, we analyzed the performance of HOMA-IR, QUICKI, and FGIR against the SIClamp by ROC analysis, defining insulin resistance as an SIClamp value <7.1 (10–14·dl·min–1·kg–1)/(µU/ml) according to previously established criteria (5). Our data on the validity of these fasting indexes are robust and are based on data from two different rat strains and from animals with different degrees of insulin sensitivity (nonpregnant and late-pregnant rats). As evidenced by closely similar AUC values, we could not establish the superiority of any of the studied indexes of insulin sensitivity in Wistar or Sprague-Dawley rats. From the ROC analysis, different cutoff values and the corresponding values for sensitivity and specificity were obtained. Sensitivity was comparable for the three studied indexes, yielding the highest sensitivity for FGIR. This index, however, exhibited lower specificity when compared with HOMA-IR or QUICKI. HOMA-IR and QUICKI showed a comparable performance in nonpregnant and pregnant animals independent of whether they were Wistar or Sprague-Dawley rats.
All of these fasting indexes are highly dependent on glucose and fasting insulin levels. Although plasma glucose assays are very reproducible, insulin values have been reported to vary considerably between different laboratories (32). Although insulin assays have been improved during the last years, a proper standardization of insulin assays is still lacking (24). In addition, the variability of insulin measurements is further increased by the high biological variability of insulin levels, a consequence of its short serum half-life and pulsatile secretion. Thus, taking into account the lack of standardization of insulin assays, it is not possible to determine absolute and universal cutoff values that define insulin resistance by using an index that depends on insulin measurements. Despite these limitations, cutoff values for the surrogate indexes can serve as reference points in long-term studies provided that they are determined under identical conditions.
In conclusion, the present study shows that simple mathematical indexes derived from a single blood fasting sample, namely HOMA-IR, QUICKI, and FGIR, can provide an easy but accurate measure of insulin sensitivity in both nonpregnant and late-pregnant Wistar and Sprague-Dawley rats. Although not intended to replace the clamp, these fasting-based indexes, in particular QUICKI and HOMA-IR, offer important advantages in estimating insulin sensitivity. Because they are obtained from single fasting blood samples, they provide a useful tool for assessing insulin sensitivity in experimental settings in which the use of anesthesia is not recommended, such as pregnancy, as well as for long-term studies in which insulin resistance has to be assessed at different time points.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. Muniyappa, H. Chen, R. H. Muzumdar, F. H. Einstein, X. Yan, L. Q. Yue, N. Barzilai, and M. J. Quon Comparison between surrogate indexes of insulin sensitivity/resistance and hyperinsulinemic euglycemic clamp estimates in rats Am J Physiol Endocrinol Metab, November 1, 2009; 297(5): E1023 - E1029. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. A. Aziz, L. S. Kenney, B. Goulet, and E.-S. Abdel-Aal Dietary Starch Type Affects Body Weight and Glycemic Control in Freely Fed but Not Energy-Restricted Obese Rats J. Nutr., October 1, 2009; 139(10): 1881 - 1889. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Mather Surrogate measures of insulin resistance: of rats, mice, and men Am J Physiol Endocrinol Metab, February 1, 2009; 296(2): E398 - E399. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |