Previous studies have reported an association between circulating insulin and incident cardiovascular disease, but limited knowledge is available on the association across subgroups. We examined the associations of fasting insulin with incident coronary heart disease (CHD) and ischemic stroke in multiple subgroups of a biracial, middle-age cohort. A total of 12,323 subjects were included in the analysis. The incidence of CHD (n = 960) and ischemic stroke (n = 445) through 2005 was determined through annual interviews, repeat examinations, and community surveillance. Serum insulin was measured at baseline. Cox regression analysis was used to estimate the hazard ratios by quintile of fasting insulin at baseline and to determine the significance of effect modification. In the minimally adjusted models (age, gender, race, and field center), the baseline fasting insulin quintile was positively associated with both incident CHD (hazard ratio per quintile insulin = 1.12, p-trend <0.0001) and ischemic stroke (hazard ratio per quintile insulin = 1.11, p = 0.0018). The adjustment for high-density lipoprotein completely attenuated the association of insulin with CHD but not with stroke. The associations of insulin with CHD were stronger in nonsmokers (p-interaction = 0.018) and in those without hypertension (p-interaction = 0.0087). The associations of insulin with stroke were stronger in women (p-interaction = 0.037), whites (compared to blacks; p-interaction = 0.036), and those without hypertension (p-interaction = 0.0027).
It is well established that type 2 diabetes mellitus increases the risk of coronary heart disease (CHD) and stroke, but the mechanisms of these associations are incompletely understood. A recent meta-analysis reviewed prospective cohort studies examining the association between circulating insulin and incident CHD and concluded that circulating insulin is a modest risk factor for CHD but that larger studies are required, in particular so the association could be examined in subgroups and effect modification (interaction) evaluated. Fewer studies have examined the prospective association of circulating insulin with stroke as a discrete end point, and all, including an earlier Atherosclerosis Risk In Communities (ARIC) Study analysis had relatively few events (<150 stroke cases). The results from these studies have been inconclusive, and no previous studies have included extensive subgroup analysis. The ARIC study has follow-up through 2005 with large numbers of CHD and stroke events. Therefore, we examined the association of baseline circulating insulin levels with incident CHD and incident ischemic stroke in the ARIC study.
Methods
The ARIC study is a multicenter prospective investigation of cardiovascular disease. White and black men and women aged 45 to 64 years were recruited in 1987 to 1989 from 4 communities: Forsyth County, North Carolina; Jackson, Mississippi; suburban areas of Minneapolis, Minnesota; and Washington County, Maryland. A total of 15,792 subjects participated in the baseline examination. Three triennial follow-up examinations were performed. The institutional review board at each field center approved the study, and all participants gave informed consent.
Of the 15,792 ARIC participants, we excluded anyone self-reporting race other than white or black and nonwhite persons from Minneapolis and Washington County (n = 103). We also excluded anyone with prevalent cardiovascular disease (n = 1,394) or prevalent diabetes (n = 1,721) at baseline. Prevalent cardiovascular disease was defined as a self-reported history of physician-diagnosed myocardial infarction (MI) or stroke; previous MI detected by electrocardiogram; or previous cardiovascular surgery or coronary angioplasty. Prevalent diabetes was defined as the presence of any of the following: fasting serum glucose of ≥126 mg/dl (7.0 mmol/L); nonfasting serum glucose level of ≥200 mg/dl (11.1 mmol/L); self-reported physician diagnosis of type 2 diabetes; or pharmacologic treatment of diabetes in the past 2 weeks. Finally, we excluded anyone who had fasted <8 hours (n = 251). The total sample size for the analysis was 12,323.
Blood samples were drawn from an antecubital vein into tubes containing ethylenediaminetetraacetic acid (lipids) or a serum separator gel (glucose and insulin). Serum insulin was measured with a radioimmunoassay (125Insulin Kit, Cambridge Medical Diagnostics, Baillerica, Massachusetts). The assay was not specific for insulin, and thus, some cross-reactivity occurred with pro-insulin. The details of the validity and reproducibility of this assay have been previously published. The total plasma cholesterol and triglyceride levels were determined using enzymatic methods. High-density lipoprotein (HDL) cholesterol was measured after dextran sulfate-magnesium precipitation of non-HDL lipoproteins. The low-density lipoprotein (LDL) cholesterol and non-HDL cholesterol levels were calculated. Trained technicians measured blood pressure 3 times. The mean of the last 2 measurements was used for analysis. Body mass index (BMI) was calculated from the participants’ height and weight measured in scrub suits. A questionnaire was used to assess the participants’ cigarette smoking history. Study participants were asked to bring all medication, vitamins, and supplements taken in the 2 weeks before the examination. Information on pharmacologic treatment of hypertension and diabetes was determined from the participants’ self-reported use of any medication to treat high blood pressure or elevated glucose and the transcription and coding of all medication names. Hypertension was defined as systolic blood pressure of ≥140 mm Hg, diastolic blood pressure of ≥90 mm Hg, or use of antihypertensive medication. Left ventricular hypertrophy was determined using the Cornell voltage criteria for the resting electrocardiogram.
Ascertainment of stroke and CHD events in the ARIC study has previously been described in detail. Events were identified using several methods. Interviewers contacted all participants annually by telephone to identify hospitalizations and death. Local death certificates and discharge lists, including International Classification of Diseases codes from local hospitals were also surveyed by the ARIC staff. For hospitalized patients, the records were abstracted for diagnostic information by trained staff. At the ARIC examinations, the technicians visually coded up to 3 12-lead electrocardiograms and evaluated the wave form evolution using side-by-side comparisons to identify silent MIs. For out-of-hospital coronary deaths, an interview with the next of kin, physician questionnaires, coroner reports, and autopsy reports were obtained when possible.
The CHD incidence was defined as definite, probable, or silent MI or definite CHD death. A committee of physicians reviewed and adjudicated all potential clinical CHD events using published criteria. For strokes, the relevant abstracted information was classified by computer algorithm and physician review. Differences were adjudicated by a second physician. Details on quality assurance for the ascertainment and classification of ischemic stroke events have been previously published. Incident ischemic stroke events were defined as validated definite or probable embolic or thrombotic brain infarctions.
All analyses were performed using Statistical Analysis Systems, version 9.1 (SAS Institute, Cary, North Carolina). For participants dying from CHD or with definite or probable CHD or stroke, the follow-up period was the interval from the baseline examination to the first event. For those with silent MI, the follow-up period ended at the midpoint between the ARIC examination that had revealed the electrocardiogram changes and the previous examination. For noncases, the follow-up period ended at the date of death, date of last contact, or December 31, 2005. Adjusted incidence rates were estimated using Poisson regression analysis.
The baseline fasting insulin levels were divided into 5 groups according to the quintile cut points of the entire distribution, and a variable coded 1 to 5 was created for the 5 categories. To test the association of fasting insulin with covariates, the covariates were regressed on this quintile variable to perform a test of linear trend. To test the association of the quintile variable with dichotomous outcomes, a chi-square test was used.
We used Cox proportional hazards regression analysis to estimate the hazard ratios (HRs) for ischemic stroke and CHD by baseline level of fasting insulin. To test the association of fasting insulin with incident CHD and ischemic stroke, the quintile variable was included in the model to perform a test of linear trend. In separate models, the fasting insulin quintile was coded as a 4-degree of freedom class variable to determine the HR for CHD or stroke for each quintile of fasting insulin compared to the lowest quintile. Finally, regression models were run using log-transformed insulin as a continuous predictor variable. Minimally adjusted models were adjusted for age, race, gender, and ARIC field centers. In other models, the association was additionally adjusted for significant predictors of incident CHD or incident stroke from analyses of ARIC-specific CHD and ischemic stroke risk scores. These covariates included systolic blood pressure, HDL cholesterol, use of antihypertensive medication, current smoking for CHD and systolic blood pressure, current smoking, use of antihypertensive medication, and prevalent left ventricular hypertrophy for ischemic stroke. Additionally, we adjusted for BMI because of the strong relation between BMI and fasting insulin and for triglycerides and non-HDL cholesterol because of the hypothesized relations among insulin, triglycerides, and very-low-density lipoprotein. We evaluated the effect modification of the fasting insulin and CHD or stroke association by including interaction terms with each covariate of interest in the minimally adjusted linear trend model of fasting insulin and by performing regression analyses stratified by the potential covariates. We considered p <0.05 to be evidence of a significant effect modification to balance the typically low power for interaction analyses with the number of interaction tests performed.
Results
Of the 12,323 subjects included in the analysis, 56.8% were women and 24.5% were black. The mean age was 53.8 years at baseline. In the 16 to 18 years of follow-up, we identified 960 incident CHD events and 445 incident ischemic stroke events. Adjusted for age, gender, race, and ARIC center, the incidence rate for CHD was 4.59 per 1,000 person-years, and the incidence rate for ischemic stroke was 2.04 per 1,000 person-years. All proposed covariates were significantly associated with the quintiles of baseline fasting insulin ( Table 1 ).
Variable | Fasting Insulin (μU/ml) | p Value ⁎ | ||||
---|---|---|---|---|---|---|
<6 | 6 to <8 | 8 to <11 | 11 to <15 | ≥15 | ||
Participants (n) | 2,680 | 2,243 | 2,736 | 2160 | 2,504 | |
Age (years) | 53.5 | 53.7 | 54.0 | 54.1 | 54.0 | 0.0002 |
Blacks | 16.0% | 18.5% | 22.2% | 28.2% | 38.2% | <0.0001 |
Women | 60.8% | 59.2% | 54.2% | 53.1% | 56.2% | <0.0001 |
Total cholesterol (mg/dl) | 207.2 | 213.7 | 214.9 | 215.8 | 217.1 | <0.0001 |
High-density lipoprotein cholesterol (mg/dl) | 61.0 | 56.5 | 52.7 | 48.7 | 45.3 | <0.0001 |
Non–high-density lipoprotein cholesterol (mg/dl) | 146.1 | 157.2 | 162.2 | 176.2 | 171.8 | <0.0001 |
Low-density lipoprotein cholesterol (mg/dl) | 127.5 | 136.1 | 138.5 | 140.2 | 140.9 | <0.0001 |
Triglycerides (mg/dl) | 93.0 | 106.3 | 119.5 | 137.0 | 162.1 | <0.0001 |
Systolic blood pressure (mm Hg) | 115 | 117 | 120 | 123 | 126 | <0.0001 |
Body mass index (m/kg 2 ) | 23.9 | 25.3 | 26.7 | 28.7 | 31.8 | <0.0001 |
Antihypertensive medication | 10.9% | 14.0% | 18.2% | 25.2% | 38.2% | <0.0001 |
Current smoker | 31.9% | 28.0% | 24.5% | 23.4% | 21.1% | <0.0001 |
Prevalent left ventricular hypertrophy (Cornell definition) | 1.2% | 1.1% | 1.7% | 2.1% | 2.9% | <0.0001 |
⁎ Test for linear trend for continuous variables; chi-square test for overall difference for dichotomous variables.
Table 2 lists the HRs of CHD and ischemic stroke by quintile of fasting insulin. We tested the proportional hazards assumption for the association of fasting insulin with CHD and ischemic stroke in minimally adjusted models and found that this assumption was not violated in either case. In a minimally adjusted model (model 1), the baseline fasting insulin quintile was a significant predictor of incident CHD (p <0.0001). The HRs for CHD were increased across every quintile of fasting insulin. When additional covariates were included individually in the regression models, most attenuated the association of insulin and CHD somewhat, although the association remained significant. However, the inclusion of HDL in the regression model completely attenuated the association between insulin and CHD. A similar pattern was seen when log-transformed insulin was used as a continuous predictor variable. In the minimally adjusted model, log-insulin was a significant predictor (p <0.0001) of CHD; however, this association was completely attenuated after adjustment for HDL. To investigate this unexpected attenuation, we ran several additional models, adjusting for ethanol intake and exercise, removing all those receiving statin treatment or other cholesterol-lowering treatment, and stratifying by baseline HDL level (data not shown). Significant attenuation of the association after adjustment for HDL level was seen in all these models. We also ran additional Cox models excluding only patients with treated diabetes (data not shown). For that analysis, in the minimally adjusted model, the HRs were increased; however, nearly complete attenuation of the association in a model adjusted for HDL was still seen.
Variable | Fasting Insulin (μU/ml) | Per Quintile | p Value ⁎ | ||||
---|---|---|---|---|---|---|---|
<6 | 6 to <8 | 8 to <11 | 11 to <15 | ≥15 | |||
Coronary heart disease | |||||||
Events (n) | 164 | 147 | 208 | 173 | 268 | — | |
Person-years | 42,075 | 35,496 | 43,157 | 33,718 | 38,623 | — | |
Hazard ratio, model 1 † (95% CI) | 1.00 (Referent) | 1.02 (0.82–1.28) | 1.13 (0.92–1.39) | 1.17 (0.94–1.44) | 1.59 (1.31–1.94) | 1.12 (1.07–1.17) | <0.0001 |
Hazard ratio, model 2 ‡ (95% CI) | 1.00 (Referent) | 1.04 (0.83–1.30) | 1.09 (0.88–1.34) | 1.05 (0.83–1.34) | 1.28 (1.02–1.62) | 1.05 (1.00–1.11) | 0.058 |
Hazard ratio, model 3 § (95% CI) | 1.00 (Referent) | 0.90 (0.72–1.13) | 0.92 (0.75–1.14) | 0.85 (0.67–1.07) | 1.01 (0.80–1.29) | 1.00 (0.95–1.06) | 0.99 |
Ischemic stroke | |||||||
Events (n) | 75 | 61 | 97 | 87 | 125 | — | |
Person-years | 42,590 | 36,192 | 43,839 | 34,340 | 39,537 | — | |
Hazard ratio, model 1 † (95% CI) | 1.00 (Referent) | 0.90 (0.64–1.26) | 1.13 (0.84–1.53) | 1.21 (0.89–1.65) | 1.45 (1.09–1.94) | 1.11 (1.04–1.19) | 0.0018 |
Hazard ratio, model 4 ¶ (95% CI) | 1.00 (Referent) | 0.89 (0.63–1.25) | 1.11 (0.82–1.51) | 1.10 (0.80–1.53) | 1.28 (0.92–1.78) | 1.07 (0.99–1.16) | 0.079 |