INTRODUCTION
Insulin resistance is one of the most important mechanisms in the development of diabetes [1], and it is closely related to the presence and severity of cardiovascular disease [2-4]. Insulin resistance is thought to promote atherogenesis in the absence of hyperglycemia [5], and contributes to atherosclerosis by way of dyslipidemia, hypertension, and other proinflammatory conditions, as well as through the effect of impaired insulin signaling at the cellular level [6,7].
Hyperinsulinemic-euglycemic clamp (HEC) is considered the gold standard for evaluating insulin resistance; however, its practical clinical use is limited by cost, labor, and ethical issues. Various indirect markers of insulin sensitivity, such as the homeostasis model assessment of insulin resistance (HOMA-IR), have been proposed that have great clinical utility [8]. In recent years, the product of triglyceride (TG) fasting plasma glucose (FPG) has been proved an alternative to HOMA-IR for estimating insulin resistance. Especially, triglyceride glucose (TyG) index has a higher sensitivity and specificity compared with HEC for estimating insulin sensitivity [9,10]. Thus, the superiority of TyG in identifying the insulin resistance can be demonstrated [11].
Body mass index (BMI) and waist circumference (WC) are simple, inexpensive, and noninvasive anthropometric parameters that are useful indicators of obesity and other metabolic disease. Evidence increasingly suggests that obesity is intimately associated with insulin resistance; therefore, a combination of obesity and TyG may identify insulin resistance stronger than other surrogate markers. Multiple cohort studies showed that increased insulin resistance was also associated with cardiovascular disease [12-14]. However, there are few studies focusing on the relationship between TyG and subclinical coronary atherosclerosis.
Our study aimed to investigate the association of TyG-related indices and risk of subclinical atherosclerosis in a Korean population. We also evaluated which of the TyG-related most accurately reflected the presence of coronary artery calcification (CAC) when compared with HOMA-IR.
METHODS
Study participants
This cross-sectional study was conducted using a medical health assessment program database obtained from the Health Promotion Center of Kangbuk Samsung Hospital, Sungkyunkwan University School of Medicine, Seoul, Korea. The purpose of the medical health checkup program is to promote the health of employees through a regular health checkup and to enhance early detection of existing disease.
The dataset comprised a total of 33,486 participants who participated in the program between January 2011 and December 2012. Of these participants, 2,692 were excluded from the study owing to self-reported history of ischemic heart disease (n=419), ischemic stroke (n=270), history of aspirin (n=38) or statin (n=1,249) use, and missing data (n=869). Finally, we enrolled 30,776 participants into our study, of which 24,755 were men and 6,021 were women.
This study was reviewed and approved by the Institutional Review Board of Kangbuk Samsung Hospital (KBC13051), and is an accordance with the Helsinki Declaration of 1975.
Anthropometric and laboratory measurements
Height and weight measurements were performed by well-trained nurses according to standardized methods. BMI was calculated by dividing participant’s weight (kg) by the square of participant’s height (m). A single examiner measured WC in participants in the standing position by locating the halfway point between the anterior iliac crest and the lower border of the bottom rib, and then measuring circumferentially. Blood pressure was measured three times consecutively at 1-minute intervals with participants in the seated position after 5 minutes of rest, with the second and third measurements averaged.
Lifestyle habits were assessed through self-questionnaire. A smoker was defined as someone who had smoked at least five packs of cigarettes in his or her lifetime. If a participant drank more than 20 g of alcohol daily, he or she was classified as an alcohol drinker. Having a regular exercise habit was defined as participating in exercise of moderate intensity at least three times a week.
FPG was measured using the hexokinase method (Hitachi Modular D2400, Roche, Tokyo, Japan). The electrochemiluminescence immunoassay (Hitachi Modular E170, Roche) was used to measure fasting insulin levels. The enzymatic colorimetric test was used to assess total cholesterol (TC) and TG concentrations. Selective inhibition method and homogeneous enzymatic colorimetric assay were used to determine high density lipoprotein cholesterol (HDL-C) and low density lipoprotein cholesterol (LDL-C) levels, respectively.
Insulin resistance was measured using the HOMA-IR according to the following equation: fasting insulin (IU/mL)×FPG (mmol/L)/22.5 [8]. Other TyG-related markers were calculated using the following formulae: TyG index=Ln [TG (mg/dL)×FPG (mg/dL)/2]; TyG-BMI=TyG index×BMI; and TyG-WC=TyG index×WC [12-14].
Measurement of coronary artery calcium score
Measurement of coronary artery calcium score (CACS) was performed by 64-slice, spiral multidetector computed tomography (MDCT) scan (GE Health Care, Tokyo, Japan) using the software HEARTBEAT-CS (Philips, Cleveland, OH, USA). The 64-slice MDCT scan was performed using the following protocol: 0.625 mm slice thickness, 120 kVP, 800 effective mAs, and 400 msec rotational speed.
The severity of CAC was assessed by the Agatston score [15]. The total CACS was determined by the sum of individual scores for the four major epicardial coronary arteries: left main, left anterior descending, left circumflex, and right coronary arteries. Technicians who performed the MDCT scans were blinded to all participant information.
The presence of CAC was defined by CACS >0. CACS was analyzed in natural logarithmic form plus 1: Ln(CACS+1). In addition, the participants were divided into four groups according to the severity of CAC: CACS=0; 0<CACS≤100; 100<CACS≤400; and CACS >400.
Statistical analysis
Continuous variables were presented as the mean±standard deviation, and categorical variables were presented as percentages. Student t-test and chi-square test were used to compare the values of the parameters for participants with and without CACS >0. Analysis of variance (ANOVA) test was used to compare the mean values of parameters according to four groups divided by the severity of CAC.
Considering 85.6% of the participants had a CACS=0, natural logarithmic transformation of CACS+1 was used for the bivariate correlation analyses, the latter performed using Pearson correlation and partial correlation analysis [16-18]. Study population was stratified into tertiles based on the following variables: HOMA-IR, TyG, TyG-BMI, and TyG-WC.
Multivariate logistic regression analysis reporting odds ratio (OR) with 95% confidence intervals was performed to estimate the presence of CAC according to tertile group of TyG index and other metabolic parameters. Four regression models were assessed: crude (unadjusted) model, model 1 adjusted for age and sex, model 2 adjusted for age, sex, systolic blood pressure (SBP), and LDL-C, and model 3 adjusted for exercise and smoking status in addition to the parameters of model 2. Statistical significance was defined by a P<0.05. All statistical analyses were performed with SPSS version 24.0 (IBM Co., Armonk, NY, USA).
RESULTS
Baseline characteristics
The baseline biochemical and clinical characteristics of the study subjects according to presence of CAC are shown in Table 1. Among 30,776 participants, the prevalence of CAC (CACS >0) was 14.4% (n=4,420), of which 2.4% (n=737) had a CACS >100. The mean participant age was 41.2±7.3 years, and participants with CAC were more likely to be older than those without CAC.
Compared with individuals without CAC, those with CAC were more likely to be male smokers with higher SBP, FPG, TC, TG, and LDL-C, and have higher metabolic parameters, including HOMA-IR, TyG, TyG-BMI, and TyG-WC (Table 1). These participants were also more likely to have lower HDL-C and participate in physical activity less frequently compared with those without CAC.
Correlation of CACS with variable parameters
Bivariate correlation analysis was performed between CACS and multiple variables (Table 2). Age showed the highest correlation with CACS. After adjustment for age, sex, SBP, and LDL-C, partial correlation analyses revealed a significant correlation between all parameters and CACS. This correlation was unchanged after adjustment for smoking and physical activity. Among metabolic parameters, TyG-BMI showed the strongest correlation with CACS, followed by TyG-WC, TyG, and HOMA-IR, respectively.
Comparison of subclinical atherosclerosis risk according to tertiles of metabolic parameters
When participants were divided into tertiles based on baseline metabolic parameters, the prevalence of CAC significantly increased from the 1st to the 3rd tertile (8.9%, 14.3%, and 19.9% for TyG index, P<0.001) (Fig. 1). In addition, when the participants were divided into four groups according to the severity of CAC, 26,356 participants (82.6%) had CACS=0, 3,683 participants (12.0%) had 0<CACS≤100, 576 participants (1.9%) had 100<CACS≤400, and 161 participants (0.5%) had CACS >400 (Fig. 2). Mean values of all four parameters increased as the severity of CAC increased from CACS=0 to higher than 400, when analyzed with ANOVA test (Fig. 2).
When calculating OR according to baseline metabolic parameter tertiles, the participants with TyG in the 3rd tertile were 1.488 times more likely to have CAC >0 and 1.644 times more likely to have CAC >100 compared with those whose TyG was in the lowest tertile, even after adjustment for confounding variables (e.g., age, sex, SBP, physical activity, and smoking status) (Table 3, Supplementary Table 1). After adjusting for traditional cardiovascular risk factors (e.g., age, sex, HTN, and LDL-C), the OR for CAC >0 and CAC >100 in participants with TyG-BMI in the 3rd tertile were 1.638 and 2.056, respectively, compared with participants whose TyG-BMI was in the lowest tertile (Table 3, Supplementary Table 1). Our results showed that TyG-related parameters were more independently associated with the presence of CAC than HOMA-IR after adjustment for other metabolic parameters.
The TyG-related metabolic parameters achieved moderate prognostic performance for CAC. The highest area under the curve (AUC) was demonstrated by TyG-WC (AUC=0.626), followed by TyG (AUC=0.617), TyG-BMI (AUC=0.616), and HOMA-IR (AUC=0.562) (Table 4).
DISCUSSION
In this large cohort study, TyG, TyG-BMI, and TyG-WC are more independently associated with the presence of subclinical atherosclerosis than HOMA-IR after adjustment for age, sex, and other clinical risk factors. Among these indices, TyG-WC and TyG-BMI, which are associated with adiposity, showed the strongest association with subclinical atherosclerosis as measured by CACS, a marker that predicts coronary heart disease events [19]. Participants with TyG-BMI and TyG-WC in the highest tertile were approximately 1.6 times more likely to have subclinical atherosclerosis after adjustment for confounding variables compared with those who had TyG-related parameters in the lowest tertile. According to receiver operating characteristic analysis, TyG-BMI and TyG-WC were the most desirable surrogate markers of CAC after adjustment of the cofactor, among the IR indices. Therefore, our data provide evidence supporting the use of TyG-related indices for predicting early cardiovascular disease.
Previous invasive studies have reported an association between TyG index and carotid atherosclerosis. Irace et al. [20] showed that, unlike the HOMA index, the TyG index remained significantly associated with carotid atherosclerosis, even after adjustment for all traditional cardiovascular disease risk factors. A study by Tripolino et al. [21] reported that increasing insulin resistance, as assessed by TyG index, is associated with atherosclerosis-prone shear stress reduction in the common carotid artery. These findings suggested that reduced sensitivity to insulin is associated with an increase in shear stress on the arterial wall, predisposing the individual to the development of atheromatous plaques. However, few studies have evaluated the relationship between TyG index and CAC.
In a small study (n=888) performed by Lee et al. [22], individuals with TyG index in the highest tertile were 3.2 times more likely to have CAC than those with TyG index in the lowest tertile. In our study, TyG-WC and TyG-BMI were more efficient than TyG to confirm its association with CAC. Our study results are consistent with previous study findings that high TyG indices are associated with an increased risk of CAC, and that this association is stronger when factoring in patient adiposity. This suggests that particular attention should be paid toward cardiovascular disease among patients with high TyG-related indices.
HOMA-IR reflects insulin resistance in our body, especially the insulin resistance in the liver as it indicates the ability of basal insulin to suppress hepatic glucose production in the fasted state [23]. Contrarily, an increase in plasma TG interferes with normal glucose metabolism in the muscle, resulting in reduced insulin sensitivity. Therefore, TyG index, unlike HOMA, appears to primarily reflect insulin resistance in muscle tissues [24-26]. Adipokine imbalance in the body not only promotes systemic dysfunction among metabolic processes, but also contributes to the development of cardiovascular disease [27]. Therefore, we hypothesized that TyG-related indices reflecting adiposity would be better associated with presence of CAC, which was subsequently supported by our study findings. In our study, we proved the robustness of our finding by analyzing data from a large number of participants (more than 30,000) and evaluating the prevalence of CAC by dividing each metabolic parameter level.
Our study has several limitations. Given that our study was cross-sectional in design, a cause-and-effect relationship could not be determined. A longitudinal study is warranted for further exploration of this relationship. Another limitation of our study is the lack of cultural and ethnic diversity of our study population. As there were interethnic differences in insulin sensitivity and in lipid profile, caution is advised when using a lipid surrogate for insulin resistance [28-30]. Additionally, the relationship among BMI, WC, and adiposity is influenced by ethnicity, supporting the use of ethnic-specific BMI and WC [31]. Another potential study limitation includes the retrospective collection of data, which could introduce various confounding factors. It is possible that our method of participant recruitment introduced selection bias into our study, as these participants were voluntarily recruited during general health assessments. Lastly, our study compared TyG-related indices with HOMA-IR, not the HEC; therefore, we were unable to provide supporting evidence favoring use of TyG-related indices, rather than HEC, as the gold standard measurement of insulin sensitivity. However, studies have shown that HOMA-IR is an appropriate alternative to the HEC [32], and that TyG index has superior sensitivity and specificity compared with the HEC for measuring insulin sensitivity [9]. Despite these limitations, our study has important clinical implications as it is the first to compare the relationship between CAC and various indirect measurements of insulin resistance.
In conclusion, TyG index may be a useful predictor of subclinical coronary artery disease in healthy Korean individuals, and is a more accurate predictor of CAC compared with HOMA-IR. It is important that physicians are aware of the increased risk of cardiovascular disease among patients with an elevated TyG-related index. Upon identification of at-risk patients, it is essential to emphasize lifestyle modification, including weight loss and excellent glycemic control, to prevent progression of cardiovascular disease. We recommend that future studies focusing on the relationship between TyG index and subclinical atherosclerosis be longitudinal in design and include diverse racial and ethnic populations, as this can affect study outcomes.
PDF Links
PubReader
ePub Link
XML Download
Full text via DOI
Download Citation
Supplement
Print
