Homeostasis Model Assessment of Insulin Resistance in Relation to the Incidence of Cardiovascular Disease

RESEARCH DESIGN AND METHODS

The SAHS (Phase II) is a population-based study of diabetes and CVD in Mexican-Americans and non-Hispanic whites. A detailed description of the study has been published previously (22). In brief, between January 1984 and December 1988, we randomly selected households from low income (barrio), middle income (transitional), and high income (suburban) census tracts in San Antonio, Texas. All men and nonpregnant women aged 25–64 years who resided in the randomly sampled households were eligible to participate. Mexican-Americans were defined as those whose ancestry derived from a Mexican national origin. All subjects gave informed consent, and the study was approved by the Institutional Review Board of the University of Texas Health Sciences Center at San Antonio.

Between October 1991 and December 1996, we conducted follow-up examinations of the baseline cohort (median follow-up time 7.5 years). The present report includes individuals who were free of documented CVD and diabetes at baseline. Vital status was obtained on 98.1% of the cohort, whereas 92% received a telephone interview, and 65% received a full clinic examination.

Clinical and laboratory measurements

At both baseline and follow-up examinations, blood specimens were obtained after a 12- to 14-h fast for determination of serum lipid and lipoprotein levels and insulin and plasma glucose concentrations. Insulin was measured with a solid-phase radioimmunoassay (Diagnostic Products, Los Angeles) (22), which shows a high degree of cross-reactivity (70–100%) with proinsulin (23). Lipid, lipoprotein, and glucose levels were measured as previously described (24). We estimated IR using the HOMA-IR index (25), which is defined as fasting insulin (μU/ml) times fasting glucose (mmol/l) divided by 22.5. HOMA-IR has been validated in normoglycemic subjects against insulin sensitivity measured directly from the euglycemic-hyperinsulinemic clamp technique (n = 12, r = 0.83, P < 0.01) (25) and has been widely used in epidemiological studies. A 75-g oral glucose tolerance test (Orangedex; Custom Laboratories) was administered, and blood specimens were obtained 30 min, 1 h, and 2 h later for determination of plasma glucose and serum insulin concentrations. Diabetes was diagnosed according to World Health Organization (WHO) criteria (26). Subjects who did not meet WHO plasma glucose criteria but who were undergoing treatment with oral antidiabetic agents or insulin were considered to have diabetes.

Anthropometric measurements (height, weight, and waist circumference) were made after participants had removed their shoes and upper garments and donned an examination gown (23). Systolic (first phase) and diastolic (fifth phase) blood pressure were measured to the nearest even digit by use of a random-zero sphygmomanometer (Hawksley-Gelman). Three readings were made for each individual, and the average of the second and third readings was used in the analysis. Subjects were considered to have hypertension if they were taking antihypertensive medication or if they had a systolic blood pressure ≥140 or a diastolic blood pressure ≥90. Lifestyle risk factors, including smoking, alcohol consumption, and participation in planned leisure-time exercise, were measured using standardized interviewer-administered questionnaires. Smoking status was categorized for the present analysis as never/former or current smoker. Alcohol consumption was parameterized as grams per month and exercise as number of sessions per week.

Determination of CVD outcomes and vital status

In the present study, we defined incident CVD as death from CVD, incident heart attack, stroke, or heart surgery, or incident angina. Cardiovascular mortality was identified and confirmed by follow-up interviews with next of kin, after which death certificates were obtained (27). Death certificates were coded according to ICD-9 by a certified nosologist; ICD-9 codes for CVD include 390–459. Cause of death was defined as the underlying cause of death. Incident cases of heart attack, stroke, and heart surgery were determined by self-report during either the telephone interview or the clinic examination. Heart attack and stroke were defined as self-reports of a physician-confirmed event. Heart surgery was defined as coronary revascularization or angioplasty. Incident angina was determined during the clinic examination (65% of the cohort) using the WHO Rose chest pain questionnaire. We confirmed 31 of 32 nonfatal CVD outcomes using chart reviews in local area hospitals. However, this validation was not complete, so our CVD outcomes should be considered as self-reports of physician-diagnosed events.

Statistical analysis

Means and SDs or proportions were presented for subjects stratified by ethnicity and sex. Differences between groups were assessed using t tests or χ² tests, with natural log transformations of skewed variables used as appropriate (Table 1). Associations between baseline metabolic and anthropometric variables were determined using Spearman correlation analysis, adjusted for age, sex, and ethnicity. Mean levels of cardiovascular risk factors, adjusted for age, sex, and ethnicity, were estimated for quintiles of HOMA-IR and fasting insulin using ANCOVA (Table 2).

Eight-year cumulative incidence rates of CVD outcomes were calculated as the proportion of disease-free subjects at baseline who developed the given CVD end point over the follow-up period. These outcomes included CVD death (0.61%), heart attack (2.01%), heart surgery (1.49%), stroke (1.19%), or angina (5.36%). Given the small number of events for individual end points, we considered combined CVD events as the outcome variable of interest. In this combined outcome variable, all subjects known to have experienced any of the CVD incidents listed above were defined as having had an event, whereas those known to have experienced none of these CVD events were defined as disease free. The incidence rate for combined CVD death, heart attack, heart surgery, and stroke was similar to the incidence rate for angina (4.47 vs. 5.36%). The inclusion of angina in the combined outcome was supported by a similar magnitude and direction of association for high versus low HOMA-IR (<2 vs. ≥2; 2 being close to the median value in the population) and risk of both angina (odds ratio [OR] 2.2, 95% CI 1.4–3.5) and nonangina CVD end points (1.6, 1.1–2.4). ORs for associations of HOMA-IR and insulin with the 8-year incidence of combined CVD outcomes were estimated using unconditional logistic regression analysis. In these models, the independent variables (HOMA-IR and fasting insulin) were analyzed as quintiles, with the lowest quintile treated as the reference category. Logistic models were adjusted for 1) age, sex, and ethnicity; and 2) age, sex, ethnicity, LDL, triglyceride, HDL, systolic blood pressure, smoking, alcohol consumption, physical activity, and waist circumference (median split). Trends across quintiles of the independent variables were assessed in separate models using a 5-level categorical variable (Figs. 1A and B; Figs. 2A and B). To examine the association of other variables with risk of CVD, we present models adjusted for 1) age, sex, and ethnicity; 2) model 1 plus smoking, alcohol, LDL, and systolic blood pressure; and 3) model 2 plus waist, HDL, triglyceride, and physical activity (Table 3). To determine whether ethnicity (Hispanic versus non-Hispanic white), waist circumference (sex-specific median split), glucose tolerance status (baseline normal glucose tolerance [NGT] versus impaired glucose tolerance [IGT]), dyslipidemia (triglyceride <200 or HDL >35 vs. triglyceride ≥200 or HDL ≤35), hypertension (no versus yes), or sex modified the association between HOMA-IR and risk of CVD, we included interaction terms in separate demographically adjusted models. In addition, we examined this issue by plotting the ORs and 95% CIs for each strata of the interaction variable under consideration (Fig. 3).

RESULTS

Baseline characteristics of nondiabetic participants in the SAHS (Phase II), stratified by sex and ethnicity, are presented in Table 1. Mexican-Americans were younger but had higher BMI and systolic blood pressure and higher concentrations of 2-h glucose, insulin, HOMA-IR, and triglyceride compared with non-Hispanic whites. Male subjects had higher waist circumference and systolic blood pressure and higher concentrations of fasting glucose and triglyceride but lower concentrations of 2-h glucose and HDL.

Anthropometric variables, fasting and 2-h glucose and fasting insulin, increased in a linear fashion across quintiles of HOMA-IR after adjustment for age, sex, and ethnicity (Table 2). A similar pattern was observed across quintiles of fasting insulin concentration (data not shown). HOMA-IR showed almost perfect correlation with fasting insulin concentration (r = 0.99). In addition, HOMA-IR was significantly associated with baseline measures of anthropometry (BMI, r = 0.51; waist circumference, r = 0.49; and waist-to-hip ratio, r = 0.32), glucose (fasting, r = 0.38 and 2-h, r = 0.30; both P < 0.0001), and cardiovascular risk factors (HDL, r = −0.32; cholesterol, r = 0.09; triglyceride, r = 0.36; systolic blood pressure, r = 0.21; and diastolic blood pressure, r = 0.24; all P < 0.0001). As would be expected, associations of fasting insulin with anthropometric and metabolic variables were very similar (data not shown), although the magnitude of each of these correlation coefficients was slightly lower compared with those for HOMA-IR.

Figure 1 presents results of multiple logistic regression analyses, with HOMA-IR modeled in quintiles. After adjustment for age, sex, and ethnicity, there was a significant trend in increasing risk for CVD across quintiles of HOMA-IR (P for trend < 0.0001) (Fig. 1A), and this association remained significant with the addition of a sex-specific median split waist circumference variable to the model (data not shown). Subjects in the highest quintile of HOMA-IR experienced a greater than twofold increased cardiovascular risk compared with those in the lowest quintile (OR 2.52, 95% CI 1.46–4.36). Although the magnitude of these associations was reduced slightly after additional adjustment for LDL, triglyceride, HDL, systolic blood pressure, smoking, alcohol consumption, exercise, and waist circumference, the trends and quintile comparisons remained statistically significant (P for trend 0.02; quintile 5 vs. quintile 1, OR 1.94, 95% CI 1.05–3.59) (Fig. 1B). The magnitude and direction of these associations were similar when stroke and coronary artery disease (CAD) were analyzed separately (HOMA-IR ≥2 vs. <2; stroke, OR 2.07, 95% CI 0.89–4.74; CAD, 1.54, 0.91–2.62). In addition, the results were not materially different when Cox regression was used to analyze the association of HOMA-IR with nonfatal CVD events (data not shown). The patterns of association for fasting insulin were similar to those for HOMA-IR (Figs. 2A and B). Although the magnitude of the insulin associations was slightly weaker in each case compared with HOMA-IR, the areas under the receiver-operating characteristic curves for HOMA-IR and fasting insulin were not significantly different (χ² = 0.20, P = 0.65).

In demographically adjusted analyses, increased risk of CVD was also predicted by 1 SD differences in triglyceride (OR 1.21, 95% CI 1.06–1.37), HDL (0.80, 0.67–0.95), systolic blood pressure (1.18, 1.01–1.38), and smoking (no versus yes: 1.69, 1.20–2.38). In a full multivariate model, HOMA-IR, age, and smoking were significant independent predictors of incident CVD in this population (Table 3). When HOMA-IR was excluded from the model, older age and smoking were significantly associated with risk of CVD (both P < 0.01), and elevated triglyceride concentration, elevated systolic blood pressure, and high waist circumference showed associations of borderline significance (P = 0.12, 0.07, and 0.12, respectively). A similar analysis using fasting insulin and glucose in place of HOMA-IR yielded very similar results (data not shown).

In Fig. 3, results of demographically adjusted logistic regression analyses are presented after stratification for potential effect-modifying variables. In this figure, HOMA-IR is modeled as a dichotomous variable (<2 vs. ≥2). Tests of heterogeneity were nonsignificant at the 0.05 level between strata of ethnicity, sex, hypertension, dyslipidemia, glucose tolerance (NGT versus IGT), and waist circumference (all P > 0.10). However, there were some differences in the point estimates of risk between categories of certain variables. Risk of progression appeared to be lower among non-Hispanic whites (compared with Mexican-Americans) and individuals with NGT, hypertension, dyslipidemia, and low waist circumference.

CONCLUSIONS

In the present study, we have demonstrated significant associations between baseline levels of both HOMA-IR and insulin and subsequent risk of CVD outcomes in a large population-based study. These associations remained significant after adjustment for multiple potential confounding variables, and there was no strong evidence of effect modification by adiposity, glucose tolerance status, dyslipidemia, hypertension, ethnicity, or sex. In separate analysis of stroke and CAD, the magnitude of association with HOMA-IR was similar. Furthermore, waist circumference was not a significant variable in multivariate analysis, suggesting that IR increases CVD risk independent of central adiposity. Finally, we were able to investigate these associations in both men and women, and this was the first published study to have included Hispanic Americans. In addition to observing significant associations between HOMA-IR and risk of CVD, we also found that fasting insulin demonstrated similar (and only marginally weaker) associations with CVD risk. These results are not surprising given the almost perfect correlation between the two variables (0.99), and they indicate that the HOMA index is probably reflecting hyperinsulinemia in this population. We chose HOMA for this analysis because it is the most widely used surrogate measure of IR. Other surrogate indexes have recently been proposed that show superior correlation with criterion measures in small groups of subjects (28–30), although HOMA has also recently been shown to perform well in validation studies (31,32).

Previous studies have reported both the presence and the absence of significant associations between insulin levels and risk of CVD (11,13,15–21). A number of review papers have highlighted possible reasons for these inconsistencies in the literature (11,14, 15). A recent meta-analysis by Ruige et al. (15) reported that the relationship between insulin and CVD was modified by ethnic background and type of insulin assay, and that there was inconsistency across studies in the adjustment for potentially confounding variables (particularly measures of adiposity). In our analysis (Figs. 1 and 2), we presented both minimally and fully adjusted models, given the present uncertainty regarding confounders in this relationship. It has also been suggested that the association is confined to middle-aged subjects (11).

The present report is one of only a few prospective studies of IR or insulin concentrations and risk of CVD to have included female participants and the second to demonstrate a significant association in women independent of other cardiovascular risk factors. Folsom et al. (16) reported a significant trend for risk of coronary heart disease across categories of fasting insulin level in women (P = 0.02) but not in men. In the four other previous studies that have included women (20,33–35), one did not present analyses separately by sex (20), and the remaining three did not find a significant effect among women (33–35), although these studies may have lacked statistical power to detect associations in stratified analyses. Additional studies of insulin concentrations and IR among women are required.

The biological mechanism to explain the association between HOMA-IR or insulin levels and risk of CVD is not well understood. It is possible that elevated circulating insulin levels have a direct detrimental effect by promoting the proliferation of vascular smooth muscle cells, a suggestion supported by investigations in animal models (36,37). Alternatively, elevated insulin levels may be a marker of the underlying cluster of abnormalities that characterizes the IR syndrome. In a recent population-based study from Finland, the IR syndrome (parameterized as a score generated by factor analysis of a collection of metabolic variables) was significantly associated with coronary heart disease and stroke (38). In contrast, Yarnell et al. (19) used cluster analysis to identify subjects with metabolic syndrome and found that there was no excess risk associated with this construct beyond that explained by the individual defining variables themselves. The most convincing evidence to date comes from the Insulin Resistance Atherosclerosis Study, in which insulin sensitivity index (S_I) was directly measured using the frequently sampled intravenous glucose tolerance test with minimal model analysis (39). In this study, S_I showed a significant inverse association with ultrasonographically measured carotid artery atherosclerosis that was independent of insulin concentrations and other CVD risk factors. Finally, it is possible that the detrimental effect of IR on CVD might operate through a mechanism of chronic subclinical inflammation (40,41).

The lack of association of male sex, LDL, and blood pressure with CVD was unexpected. This somewhat unique pattern of association may be related to the fact that a large proportion of the SAHS cohort was Mexican-American, an ethnic group that may have a CVD risk profile that differs slightly from that in European Americans. In addition, male sex was not significant because women are at greater risk of angina than men, who in turn are at greater risk of other events included in our outcome variable, and thus they offset each other.

In conclusion, we have demonstrated that HOMA-IR was significantly and independently associated with risk of CVD outcomes in Mexican-American and non-Hispanic white men and women in the SAHS. In addition, we found that there were only minor differences in the prediction of CVD when fasting insulin was used in place of HOMA-IR. The literature remains highly inconsistent on the topic of IR and CVD, however, and additional studies are required, particularly among women and minority populations.