Failure of Homeostatic Model Assessment of Insulin Resistance to Detect Marked Diet-Induced Insulin Resistance in Dogs
- Marilyn Ader1,
- Darko Stefanovski1,
- Joyce M. Richey2,
- Stella P. Kim1,
- Cathryn M. Kolka1,
- Viorica Ionut1,
- Morvarid Kabir1 and
- Richard N. Bergman1
- 1Diabetes and Obesity Research Institute, Cedars-Sinai Medical Center, Los Angeles, CA
- 2Keck School of Medicine, University of Southern California, Los Angeles, CA
- Corresponding author: Marilyn Ader, .
Accurate quantification of insulin resistance is essential for determining efficacy of treatments to reduce diabetes risk. Gold-standard methods to assess resistance are available (e.g., hyperinsulinemic clamp or minimal model), but surrogate indices based solely on fasting values have attractive simplicity. One such surrogate, the homeostatic model assessment of insulin resistance (HOMA-IR), is widely applied despite known inaccuracies in characterizing resistance across groups. Of greater significance is whether HOMA-IR can detect changes in insulin sensitivity induced by an intervention. We tested the ability of HOMA-IR to detect high-fat diet–induced insulin resistance in 36 healthy canines using clamp and minimal model analysis of the intravenous glucose tolerance test (IVGTT) to document progression of resistance. The influence of pancreatic function on HOMA-IR accuracy was assessed using the acute insulin response during the IVGTT (AIRG). Diet-induced resistance was confirmed by both clamp and minimal model (P < 0.0001), and measures were correlated with each other (P = 0.001). In striking contrast, HOMA-IR ([fasting insulin (μU/mL) × fasting glucose (mmol)]/22.5) did not detect reduced sensitivity induced by fat feeding (P = 0.22). In fact, 13 of 36 animals showed an artifactual decrease in HOMA-IR (i.e., increased sensitivity). The ability of HOMA-IR to detect diet-induced resistance was particularly limited under conditions when insulin secretory function (AIRG) is less than robust. In conclusion, HOMA-IR is of limited utility for detecting diet-induced deterioration of insulin sensitivity quantified by glucose clamp or minimal model. Caution should be exercised when using HOMA-IR to detect insulin resistance when pancreatic function is compromised. It is necessary to use other accurate indices to detect longitudinal changes in insulin resistance with any confidence.
Insulin sensitivity can be accurately quantified from the hyperinsulinemic glucose clamp (1) or intravenous glucose tolerance test (IVGTT) (2), while indices derived from the oral glucose tolerance test (e.g., Refs. 3,4) may be confounded by the influence of glucose absorption (5). The use of surrogate indices based solely on fasting measurements are often considered cost-effective in both small- and large-scale clinical and genetic studies (e.g., Refs. 6,7). One widely used surrogate is the homeostatic model assessment of insulin resistance (HOMA-IR) (8). HOMA-IR ([fasting insulin × fasting glucose]/22.5) is based on computer estimation of fasting values expected for a given degree of insulin resistance and pancreatic β-cell function and was derived from hyperglycemic clamp data obtained from normal healthy subjects (9). A value of 1.0 is considered normal, with higher indices reflecting a greater degree of insulin resistance. Subjects undergoing severe caloric restriction appear to be extremely insulin-sensitive, as reflected in HOMA-IR as low as 0.29 (10), while morbidly obese subjects with type 2 diabetes may exhibit >20-fold greater values of the HOMA-IR index (11).
It is unknown whether HOMA-IR can accurately quantify longitudinal improvements in, or deteriorations of, insulin sensitivity. HOMA-IR is calculated from fasting glucose and insulin, but resistance often develops in the absence of elevated glycemia (12). If insulin resistance was accompanied by a β-cell defect, fasting hyperinsulinemia may not occur, and hormone levels may not exceed normal values (13). Furthermore, when HOMA-IR is used under conditions when fasting glucose is unchanged, the index is equivalent to fasting insulin, and it is unclear whether the surrogate measure offers any significant advantage over the hormone measurement alone.
In this study, we test the accuracy of HOMA-IR to detect longitudinal development of insulin resistance in dogs fed a fat-supplemented diet. Resistance was quantified by euglycemic clamps and minimal model analysis and compared with diet-associated changes in HOMA-IR.
Research Design and Methods
Procedures were performed on 36 male mongrel dogs (28.2 ± 0.5 kg); data were pooled from previous published studies (14–16). All procedures were approved by the University of Southern California Institutional Animal Care and Use Committee.
Dogs were fed a weight-maintaining diet of 3,885 kcal/day (38% carbohydrates, 26% protein, and 35% fat). After 2 to 3 weeks of baseline metabolic testing, animals were switched to a hypercaloric high-fat diet (5,236 kcal/day, derived from 28% carbohydrates, 19% protein, and 51% fat). Animals were fed at 9 a.m. and given 3 h to consume available ration.
All dogs underwent comprehensive assessment of insulin sensitivity prior to and after a 6-week period of high-fat feeding, with both the hyperinsulinemic-euglycemic clamp (EGC) and minimal model analysis of the IVGTT. Also, HOMA-IR was calculated from fasting glucose and insulin at each of these two study periods. Development of diet-induced obesity was confirmed by magnetic resonance imaging. All procedures were performed after overnight fasting, and all sensitivity tests occurred in conscious, unstressed animals.
Tracer (25 µCi + 0.25 µCi/min; PerkinElmer NEN) was infused into a peripheral vein and continued for a 90-min equilibration period. After basal blood sampling, somatostatin (1 μg/min/kg; Bachem) was infused i.v. along with insulin (regular purified pork; Lilly) from t = 0 to 180 min at either 0.75 mU/min/kg (n = 19) or 1.15 mU/min/kg (n = 17) to induce hyperinsulinemia. Euglycemia was maintained by variable rate 50% dextrose infusion, spiked with [3-3H]glucose (specific activity: 2.2 μCi/g). Blood samples were obtained from a jugular vein catheter (n = 8) or sampling port (n = 17) surgically implanted >1 week prior to testing or via peripheral limb intracatheters (n = 11). Samples were drawn every 10–15 min for 180 min and assayed for glucose, insulin, and tracer.
Glucose (0.3 g/kg) and insulin (0.03 units/kg) were injected intravenously at t = 0 and 20 min, respectively, and 31 blood samples were drawn from t = −30 to 180 min and assayed for glucose and insulin.
Blood Sampling and Assays
Blood was collected in lithium- and heparin-coated tubes containing EDTA, centrifuged, and plasma stored at −80°C. Glucose was assayed by the glucose oxidase technique (YSI Model 2300; YSI Life Sciences), with intra-assay coefficient of variation (CV) <1%. Insulin was measured in duplicate by ELISA (Linco Research, St. Charles, MO), with detection limit of 5 pmol/L and intra- and interassay CVs of 2 ± 1 and 5 ± 1%, respectively.
Calculations and Data Analysis
Pooled data from multiple studies with similar feeding regimens were used to test the ability of HOMA-IR to detect changes in insulin sensitivity induced by fat feeding, as described below.
EGCs were used to quantify whole-body insulin sensitivity (SICLAMP), as well as sensitivity of periphery and liver (17). SICLAMP was calculated as:where ΔGINF and ΔINS are respective increments in glucose infusion rate and insulin during exogenous insulin infusion, and GLUss is the steady-state glucose concentration (final 30 min of the EGC). In a subset of animals (n = 31), rates of glucose uptake (Rd) and hepatic glucose output (HGO) were calculated using Steele’s equations modified for use with labeled glucose infusion (18). Peripheral (SIpCLAMP) and hepatic (SIHGOCLAMP) sensitivity indices were defined as:where ΔRd and ΔHGO are the changes in Rd and HGO from basal to steady state.
The insulin sensitivity index (SImm) was calculated from minimal model analysis of the IVGTT (MINMOD Millennium, ver. 6.02; Ref. 19).
where insulin and glucose are plasma concentrations after an overnight fast (8) (in μU/mL and mmmol, respectively; mean of three samples, assayed in duplicate).
Glucose-Stimulated Insulin Response
We calculated glucose-stimulated insulin response (AIRG) from the IVGTT, as the incremental insulin area under the curve from 0–10 min after glucose injection.
Statistics (t test and ANOVA, with Tukey post hoc analysis when overall significance was detected) were performed using MINITAB (ver. 13.32; State College, PA). Statistical significance was set at P ≤ 0.05.
Dogs had fasting glucose and insulin concentrations in the normal range (89.4 ± 1.2 mg/dL and 11.6 ± 1.0 μU/mL, respectively), with body weight within a narrow range (28.2 ± 0.5 kg; CV = 10%). Body adiposity varied widely (three- to sixfold range for visceral and subcutaneous fat mass, respectively). SICLAMP ranged from 5.9–72.9 dL/min/kg/μU/ml (12.4-fold variation; mean: 34.5 ± 2.6), and similar variability was observed in SImm (9.1-fold range; mean: 4.53 ± 0.39 × 10−4 min−1/μU/mL). HOMA-IR was 2.67 ± 0.16, and the range from most-sensitive (0.86) to most-resistant animal (4.64) was narrow (5.4-fold). Baseline HOMA-IR did not correlate with either clamp- or IVGTT-based indices (Fig. 1).
Effect of Fat Feeding
Consumption of a high-fat diet induced substantial increases in both total (64 ± 7%) and regional adiposity (44 ± 5 and 110 ± 12% for visceral and subcutaneous depots, respectively; P < 0.0001). Diet-induced obesity caused significant insulin resistance (Fig. 2A). SICLAMP declined from 34.5 ± 2.6 to 22.9 ± 1.6 dL/min/kg/μU/mL (P < 0.0001), and resistance developed in 32 of 36 animals tested, both in insulin action on Rd (SIpCLAMP; −51 ± 12%; P < 0.0001) and HGO (SIHGOCLAMP; −96 ± 45%, P < 0.0001). Insulin resistance was confirmed by IVGTT (SImm: 4.5 ± 0.4 to 2.9 ± 0.2 × 10−4 min−1/μU/mL; P < 0.0001; Fig. 2B). Clamp-based whole-body insulin sensitivity and SImm were well-correlated (P = 0.001).
Diet-induced resistance was not detected by HOMA-IR (Fig. 2C). Twenty-three of 36 animals showed an increase in HOMA-IR, denoting increasing insulin resistance (ΔHOMA-IR = 1.1 ± 0.2); 13 (36% of all animals) showed an artifactual decrease (ΔHOMA-IR = −0.8 ± 0.1). The overall change in HOMA-IR was minimal, with a large SE (ΔHOMA-IR = 0.4 ± 0.2); change in HOMA-IR failed to attain statistical significance despite the large number of animals (P = 0.22).
The small increase in HOMA-IR with fat feeding was, by definition, the result of observed effects of diet on fasting glucose and insulin (Eq. 2). Glycemia was unaffected by fat feeding, but animals did develop modest fasting hyperinsulinemia (15.1 ± 1.3 μU/mL; P = 0.0005). HOMA-IR offered no additional advantage over insulin values per se in predicting development of insulin resistance.
Impact of AIRG on HOMA-IR Accuracy
Development of diet-induced insulin resistance led to increased AIRG (prefat: 609 ± 38; postfat: 850 ± 52 μU/mL). HOMA-IR displayed a strong, positive correlation with AIRG (Fig. 3), consistent with the known reciprocal relationship between insulin sensitivity and insulin responsiveness (20).
The ability of HOMA-IR to detect insulin resistance is centered on the development of hyperinsulinemia, which is proportional to insulin responsiveness (i.e., AIRG). Using both baseline and post–fat feeding data, we examined the influence of AIRG on HOMA-IR accuracy by segregating data into groups below or above the mean AIRG (730 μU/mL, mean of 72 assessments; Fig. 4). In animals with a robust AIRG (Fig. 4A), HOMA-IR was correlated with SICLAMP (r = −0.36; P = 0.046), but HOMA-IR failed to accurately estimate insulin sensitivity in animals with lesser β-cell function (P = 0.575; Fig. 4B). Failure to detect changes in sensitivity was also observed with the quantitative insulin sensitivity check index (QUICKI) (21), a related surrogate (P = 0.014 and P = 0.268 for correlations in animals with high and low AIRG, respectively), and detection was not improved with second-generation HOMA2-IR estimates (http://www.dtu.ox.ac.uk/homacalculator). The ability of HOMA-IR (or QUICKI) to detect development of insulin resistance is dependent upon β-cell function.
Accurate quantitation of insulin sensitivity is critical for comprehensive metabolic phenotyping essential in the search for genetic loci for diabetes risk. The euglycemic clamp (1) measures the steady-state response to hyperinsulinemia to quantify whole-body insulin sensitivity, as well as indices of both peripheral and hepatic sensitivity. The IVGTT-based minimal model (2,19) analyzes the dynamic glucose-insulin relationship following glucose injection to yield the insulin sensitivity index, as well as measures of β-cell function and insulin clearance (e.g., Refs. 22,23). Equivalence of the clamp- and IVGTT-based indices of sensitivity has been demonstrated (17). Surrogate measures of insulin sensitivity (8,21) require only measurement of fasting glucose and insulin and represent the simplest approach. But the current study reveals serious shortcomings of surrogate measures HOMA-IR, as well as QUICKI, to detect confirmed development of insulin resistance.
At baseline, animals exhibited a wide range of insulin sensitivity, and clamp-based values were well-correlated with minimal model indices. But there was no significant relationship between either clamp- or minimal model–based measures of baseline sensitivity and corresponding HOMA-IR (Fig. 1). HOMA-IR detected insulin resistance in animals where clamps revealed normal sensitivity. The median value of HOMA-IR was 2.50, but animals considered resistant by virtue of HOMA-IR exceeding the median were not resistant by clamp-based SICLAMP (32.7 ± 3.8 vs. 36.2 ± 3.6 dL/min/kg/μU/mL for animals with HOMA-IR values below the median; P = 0.81).
The interventional methods were both effective in detecting deterioration of whole-body insulin sensitivity induced by high-fat diet. SICLAMP and SImm decreased by 34 and 37%, respectively (Fig. 2A and B), and clamp data indicate resistance developed at both hepatic and peripheral tissues. In stark contrast, HOMA-IR was of surprisingly limited utility in detecting this diet-induced resistance (Fig. 2C). There was a modest increase in HOMA-IR after fat feeding; the mean change was small, with a large SE, due in part to small, diet-induced reduction in fasting glycemia in 18 of 36 animals (ΔHOMA-IR = 0.2 ± 0.4), but this change was not statistically significant (P = 0.51). Moreover, 13 of 36 animals tested showed an artifactual decrease in HOMA-IR (ΔHOMA-IR = −0.8 ± 0.1).
We observed a strong correlation between HOMA-IR and the AIRG (Fig. 3). We segregated HOMA-IR and SICLAMP data into two groups, reflecting animals with high (higher than mean; n = 32) or low AIRG (lower than or equal to mean; n = 40). HOMA-IR accuracy was indeed dependent on the magnitude of the insulin secretory response (Fig. 4). When β-cell function was robust (i.e., AIRG exceeded the mean), HOMA-IR yielded an accurate estimation of insulin sensitivity. The limitations of HOMA-IR are exposed in subjects that display a lesser insulin secretory response. No animals in this study were diabetic; the low AIRG is a likely reflection of high insulin sensitivity in those animals.
HOMA-IR accuracy is weakest for estimating insulin sensitivity under conditions in which insulin secretory function is in the normal range, but less robust. Thus, HOMA-IR should be viewed with skepticism if β-cell function is not known a priori in all individuals. Relationship between HOMA-IR as a quantitative trait and specific variants may not indicate genetic signals for insulin resistance per se, but may more represent signals for β-cell function, as HOMA-IR may reflect islet cell function or metabolic clearance of insulin, rather than insulin resistance itself (24).
Acknowledgments. The authors thank Ed Zuniga and Edgardo Paredes for gentle care and treatment of the experimental animals and Dr. Erlinda Kirkman for expert veterinary care and advice.
Duality of Interest. No potential conflicts of interest relevant to this article were reported.
Author Contributions. M.A. conceived of the study, designed the experiments, analyzed and interpreted the data, and wrote the manuscript. D.S. assisted with data and statistical analysis. J.M.R., S.P.K., C.M.K., V.I., and M.K. performed experiments and MRI analysis. R.N.B. conceived of the study, interpreted the data, and assisted with manuscript preparation. M.A. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Prior Presentation. This study was presented in poster form at the 71st Scientific Sessions of the American Diabetes Association, San Diego, California, 24–28 June 2011.
- Received August 8, 2013.
- Accepted December 13, 2013.
- © 2014 by the American Diabetes Association.
Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered. See http://creativecommons.org/licenses/by-nc-nd/3.0/ for details.