Skip to main content

Triglyceride-glucose index in the development of peripheral artery disease: findings from the Atherosclerosis Risk in Communities (ARIC) Study

Abstract

Background

It remains unclear whether triglyceride-glucose (TyG) index, a surrogate marker of insulin resistance, is prospectively associated with incident peripheral arterial disease (PAD).

Methods

We included 12,320 Atherosclerosis Risk in Communities Study participants (aged 54.3 ± 5.7 years) free of a history of PAD at baseline (visit 1: 1987–1989). The TyG index was determined using ln (fasting triglycerides [mg/dL] × fasting glucose [mg/dL]/2), and measured at 5 visits between 1987 and 2013. Incident PAD was defined as the first hospitalization with PAD diagnosis or a new onset of measured ABI < 0.90 during follow-up visits. We quantified the association of both baseline and trajectories of TyG index with incident PAD using Cox regression and logistic regression analysis, respectively.

Results

Over a median follow-up of 23 years, 1300 participants developed PAD. After adjustment for traditional PAD risk factors, each 1-SD (0.58) increase in TyG index was associated with an 11.9% higher risk of incident PAD [hazard ratio, 1.119 (95% CI, 1.049–1.195)]. Results were similar when individuals were categorized by TyG index quartiles [hazard ratio, 1.239 (95% CI, 1.028–1.492); comparing extreme quartiles]. Four distinct trajectories of stable TyG indexes at various levels along the follow-up duration were identified [low (22.2%), moderate (43.2%), high (27.5%), and very high (7.1%) trajectory groups]. Compared with those with a TyG index trajectory at a low level, those participants with TyG index trajectories at high and very high levels had an even greater risk of future incident PAD [odds ratio (95%CI): 1.404 (1.132–1.740) and 1.742 (1.294–2.344), respectively] after multivariate adjustments for traditional PAD risk factors.

Conclusions

Higher TyG index is independently associated with an increased risk of incident PAD. Long-term trajectories of TyG index help identify individuals at a higher risk of PAD who deserve specific preventive and therapeutic approaches.

Trial registration: Clinical trial registration number: The ARIC trial was registered at clinicaltrials.gov as NCT00005131.

Background

Peripheral artery disease (PAD) is an important manifestation of systemic atherosclerosis affecting an estimated over 200 million people worldwide [1, 2]. It is associated with significant cardiovascular morbidity and mortality, with a variable spectrum of symptoms from none to severe when patients present with claudication or critical limb ischemia [3]. Despite its high prevalence and well-described adverse outcomes, the pathobiology of PAD is incompletely understood. There is, therefore, the need to identify potential biomarkers that could predict the risk of PAD to facilitate diagnosis and timely intervention at early stages of atherosclerosis.

Insulin resistance (IR), a pathophysiological state characterized by the attenuated insulin sensitivity of peripheral tissues, is the key feature of metabolic syndrome and type 2 diabetes [4]; and it contributes significantly to the development of atherosclerotic cardiovascular disease [5]. However, the role of IR in PAD has been inadequately explored, compared with that of other atherogenic mechanisms such as inflammation [6,7,8]. Moreover, risk factors for PAD have not been as thoroughly investigated as those for coronary heart disease (CHD) [9]. The triglyceride-glucose (TyG) index, which is calculated using fasting triglycerides (TG) and fasting glucose, is a reliable measure of IR [10, 11]. Growing evidence has demonstrated that the TyG index is related to morbidity and mortality of cardiovascular disease in the general population, patients with and those without diabetes [12, 13]. This is possibly because elevated TyG index by itself contributes to systemic arterial atherosclerosis, including carotid atherosclerosis and coronary artery calcification, an established marker of subclinical atherosclerosis [14, 15]. Previous studies evaluating association of the TyG index with markers of PAD, incident PAD or their complexity are inherently limited by small sample sizes, a retrospective study or a cross-sectional analysis, and by the use of measured TyG index at a single time point [16,17,18]. So far long-term specific prospective studies on the relationship between PAD and the TyG index, and its trajectory derived from the multiple measurements over time have not been performed.

We hypothesized that dynamic changes in IR over decades might modify the development of PAD. Accordingly, we used the data from the Atherosclerosis Risk in Communities (ARIC) study to evaluate the association of the TyG index with PAD and to determine the influence of baseline TyG index and different trajectories of its change over 20 years on the development of PAD.

Methods

Study population

The ARIC Study is a prospective cohort study that enrolled 15,792 participants aged 45 to 64 years, recruited between 1987 to 1989 from 4 US communities (Forsyth County, North Carolina; Jackson, Mississippi; eight northern suburbs of Minneapolis, Minnesota; and Washington County, Maryland), aimed at investigating the natural history, etiology, and clinical manifestations of atherosclerotic disease in black and white men and women. Cohort exams were conducted at visit 1 (1987–1989), visit 2 (1990–1992), visit 3 (1993–1995), visit 4 (1996–1998), visit 5 (2011–2013) and visit 6 (2016–2017); other exams are ongoing. Details about the study design have been previously described [19]. Written informed consent was obtained from all ARIC participants, and the ARIC study was approved by the institutional review boards at each site.

We excluded participants who had PAD diagnosis at baseline (n = 613); those who had missing data regarding PAD (n = 555); and those who had missing data regarding other covariates of interest (n = 2038). We also excluded participants who had no follow-up information on PAD (n = 266). This resulted in a final sample of 12,320 participants for the analysis of association between baseline TyG index and incident PAD. We further excluded those participants with fewer than three valid TyG index during follow-up visits; the remaining 9097 participants were included in the analysis of association between TyG index group-based trajectory and incident PAD (Fig. 1).

Fig. 1
figure1

Flowchart for selecting the Atherosclerosis Risk in Communities Study participants for analysis

Data collection at baseline

Trained interviewers collected information using standardized questionnaires on demographic, lifestyle, and detailed medical information at visit 1. Age, sex, race, educational level, physical activity, drinking and smoking status were self-reported. Specifically, physical activity was measured using the modified Baecke Physical Activity Questionnaire [20]. For each activity type, information related to the duration (hours per week) and frequency (number of weeks per month) were collected. Based on the Baecke Questionnaire scoring of summary estimates, index scores for sport, and leisure-time ranging from 1 to 5 (reflecting the highest activity level) were calculated. Leisure-time index score was based on the frequency of watching television, walking and bicycling (including to work), or shopping during leisure time whereas sport index score was based on the frequency of any sports and associated with sweat during leisure time [20]. Educational attainment was categorized as basic (less than high school), intermediate (high school graduate or vocational school), and advanced (college, graduate school, or professional school). Smoking and drinking status were classified as current, former, or never. Body mass index (BMI) was calculated as weight in kilograms divided by height in meters squared. Seated blood pressure represented the mean of the last two of three measurements using a random-zero sphygmomanometer after a 5-min rest. Hypertension was defined as systolic blood pressure readings ≥ 140 mmHg or diastolic blood pressure readings ≥ 90 mmHg, or use of antihypertensive drugs in the previous two weeks. Diabetes was defined as a fasting glucose level ≥ 126 mg/dL (≥ 7.0 mmol/L), a non-fasting glucose level ≥ 200 mg/dL (≥ 11.1 mmol/L), self-reported physician diagnosis of diabetes, or any use of antidiabetic drugs. History of CHD and stroke were both determined according to participants’ self-report or relevant measurements [8]. Medications were determined through self-reported usage in the previous two weeks and inspection of medication containers that participants brought to the visit.

Four biomarkers of inflammation and coagulation were measured at visit 1: fibrinogen, von Willebrand factor, factor VIII activity, and white blood cell count, from blood stored at − 70 °C using standardized protocols that have been described previously [21]. Plasma total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), and TG were also measured using automated enzymatic procedures, and low-density lipoprotein cholesterol (LDL-C) was calculated using the Friedewald equation when the concentration of TG is < 400 mg/dL [22]. The estimated glomerular filtration rate (eGFR) was calculated using the Chronic Kidney Disease Epidemiology Collaboration creatinine equation [23]. The TyG index was calculated as ln(fasting TG [mg/dL] × fasting glucose [mg/dL]/2). The ankle‐brachial index (ABI), as a surrogate marker of PAD, was calculated as the average of 2 resting ankle systolic pressure readings divided by the average of 2 resting brachial systolic pressure readings.

Ascertainment of incident PAD

Based on previous literature [8, 24, 25], PAD-related hospitalizations were identified by the following International Classification of Diseases Ninth Revision (ICD-9) discharge codes: peripheral vascular disease, unspecified (443.9); atherosclerosis of native arteries of the extremities, unspecified (440.20); atherosclerosis of native arteries of the extremities with intermittent claudication (440.21); atherosclerosis of native arteries of the extremities with rest pain (440.22); atherosclerosis of native arteries of the extremities with ulceration (440.23); atherosclerosis of native arteries of the extremities with gangrene (440.24); other atherosclerosis of native arteries of the extremities (440.29); atherosclerosis of bypass graft of the extremities (440.3); chronic total occlusion artery extremities (440.4); atherosclerosis of other specified arteries (440.8); coexisting leg amputation (84.11, 84.12, 84.15, 84.17); leg artery revascularization (38.18, 39.25, 39.29, 39.50); lower extremity ulcer and gangrene (707.1x). Incident PAD was defined as the first hospitalization with diagnosis of PAD or a new onset of measured ABI < 0.90 during follow-up visits [3]. Critical limb ischemia (CLI), the severe form of PAD, was based on the discharge codes (84.11, 84.12, 84.15, 84.17, 707.1x). The follow-up period was set as the time from visit 1 (baseline) to the incidence of PAD, or loss to follow-up, or September 30, 2015, whichever came first. September 2015 was chosen as the last point of follow-up to maintain diagnostic consistency since ICD codes switched from ICD-9 to ICD-10 in October 2015.

Statistical analysis

Normally distributed continuous data were expressed as mean ± SD, and the non-normally distributed continuous data, otherwise, were expressed as the median (interquartile range). Categorical data were expressed as numbers (percentage). Differences among groups were evaluated using analysis of variance (ANOVA) or Kruskal–Wallis h-test when appropriate for the continuous variables, and the χ2 test for the categorical variables. Kaplan–Meier estimates were used to compute cumulative incidence of incident PAD by TyG index quartiles and the differences in estimates were compared using the log-rank procedure. Cox proportional hazards regression model was used to calculate hazard ratios and 95% CIs between TyG index and time to incident PAD. Three multivariate models with progressive degrees of adjustment were used to adjust for potential confounders of PAD. Model 1 was adjusted for age, sex, and race. Model 2 was further adjusted for other clinical variables, including antihypertensive drugs, BMI, diabetes, drinking status, education level, sport and leisure time index scores, lipid-lowering drugs, SBP, and smoking status. Model 3 was additionally further adjusted for serum parameters, including eGFR, factor VIII activity, fibrinogen, LDL-C, von Willebrand factor, and white blood cell count. To determine whether our results were independent of the effects of potential drugs on TyG index, we also conducted a sensitivity analysis and repeated the aforementioned analyses by excluding those participants that were on any lipid- or glucose-lowing drugs. We further used a restricted cubic spline regression model with 3 knots to assess the nonlinear dose–response association between baseline TyG index and incident PAD. Subgroup analyses were performed stratifying by age, sex, race, smoking status, BMI, hypertension, and diabetes at baseline, respectively. Group-based trajectory analysis was designed to identify clusters of individuals with similar patterns of change over time [26]. We used latent class models to identify different patterns of longitudinal TyG index levels within the ARIC Study participants who had at least 3 TyG index measurements during follow-up visits. Models were fit using R 3.6.1 based on R package tidyLPA. We tested models with groups ranging from 2 to 5 and examined different criteria including Bayesian information criteria to assess the optimal number of trajectories so that no group included less than 5% of participants (Additional file: Table S1). Participants were assigned to the trajectory group for which they had the greatest posterior predictive probability. All the final models classified participants into trajectory groups with good discrimination: the mean probability of final group membership was 0.89. To estimate the association of TyG index trajectory groups with incident PAD, trajectory group was included as an independent variable in a logistic regression model examining predictors incident PAD at follow-up.

All analyses were conducted in SPSS version 23 (SPSS, Inc, Chicago, Illinois). A two-sided P value of < 0.05 was considered statistically significant.

Results

Baseline characteristics according to quartiles of TyG index

The average age of all the participants was 54.3 ± 5.7 years, 5693 (46.2%) were men, 3131 (25.4%) were current smokers, 1378 (11.2%) were diabetic, and 5192 (42.1%) had hypertension. The mean TyG index was 8.7 ± 0.6 (Table 1). We categorized the included population into 4 groups based on the quartiles of baseline TyG index (Table 1). Participants with a higher TyG index were older, and more often male and white; they had higher levels of BMI, SBP, DBP, TC, LDL-C, TG, fasting glucose, fibrinogen, Von Willebrand factor, factor VIII activity and white blood cell count; they had lower levels of education, physical activity scores (leisure time and sport index scores), HDL-C, and eGFR; they were also less frequently current drinkers (all P < 0.001). Likewise, participants in a higher TyG index quartile had a higher prevalence of hypertension, diabetes, CHD, and stroke; and were more prone to take antihypertensive drugs, and lipid-lowering drugs (all P < 0.001).

Table 1 Baseline characteristics of study participants by quartiles of TyG index

Association between baseline TyG index and incident PAD

During a follow-up of 23.0 (14.7, 25.5) years, 1300 incident cases (10.6%) of PAD were observed. As Table 1 shows, the risk of incident PAD increased with increasing quartiles of TyG index [quartiles 1–4: 222 (7.2%) vs. 278 (9.0%) vs. 359 (11.7%) vs. 441 (14.3%); P < 0.001]. A similar trend in incident CLI was observed. In the multivariate model that measured TyG index as a continuous variable, a 1-SD increase (corresponding to 0.58) in TyG index was associated with an 11.9% higher risk of incident PAD after full adjustment for the potential confounders [hazard ratio, 1.119 (95% CI, 1.049–1.195); P = 0.001; Table 2]. Results were similar when we categorized individuals by TyG index quartiles: the highest risk of incident PAD was observed in the participants with the highest TyG index quartile, in 3 different adjusted models (all  P < 0.05, Table 2). In the final model, the hazard ratios (95%CIs) for incident PAD comparing the second, third, fourth quartiles of TyG index with the first quartile were 1.040 (95% CI, 0.869–1.245), 1.208 (95% CI, 1.013–1.441), and 1.239 (95% CI, 1.028–1.492), respectively (model 3 in Table 2; Fig. 2). These findings were not modified when the analysis included only participants without any lipid- or glucose-lowering medication (Additional file 1: Tables S2). Figure 3 shows the restricted cubic splines of the risk of incident PAD across levels of TyG index. Consistent with the analysis using quartiles of sample distribution, the risk of incident PAD increased in participants with a higher TyG index. However, there was no significant difference for the risk of incident PAD in participants with TyG index < 8.6 (Fig. 3).

Table 2 Risk of incident PAD for baseline TyG index
Fig. 2
figure2

Cumulative incidence of incident PAD by quartiles of baseline TyG index. Cumulative incidence curves are statistically different (log-rank P < 0.001)

Fig. 3
figure3

Adjusted hazard ratios of incident peripheral artery disease (PAD) by baseline triglyceride-glucose (TyG) index. Each hazard ratio was computed with a TyG index level of 8.6 as the reference. The hazard ratio was adjusted for age, antihypertensive medication, body mass index, diabetes, drinking status, education level, estimated glomerular filtration rate, factor VIII activity, fibrinogen, leisure time and sport index scores, lipid-lowering medication, low-density lipoprotein cholesterol, race, sex, smoking status, systolic blood pressure, von Willebrand factor, and white blood cell count. Red solid line represents the hazard ratio of TyG index across the whole range. Red dotted lines represent the 95% CI. Blue dotted line is the reference line as hazard ratio = 1. Histograms represent the frequency distribution of baseline TyG index

When participants were stratified by age (≤ 54 or > 54 years), sex (male or female), race (white or black), smoking status (current, former or never), body mass index (< 30 or ≥ 30 kg/m2), hypertension (yes or no), and diabetes (yes or no), the association between TyG index and incident PAD remained consistent (all P for interactions > 0.05; Fig. 4).

Fig. 4
figure4

Subgroup analysis of the association between baseline TyG index and incident PAD. Cox regression after adjustment for antihypertensive medication, drinking status, education level, estimated glomerular filtration, factor VIII activity, fibrinogen, leisure time and sport index scores, lipid-lowering medication, low-density lipoprotein cholesterol, systolic blood pressure, von Willebrand factor, and white blood cell count was performed in subgroups according to age (≤ 54 or > 54 years), gender (male or female), race (White or black), smoking status (current or former or never), body mass index (BMI; < 30 or ≥ 30 kg/m2), hypertension (yes or no), and diabetes (yes or no)

Association between TyG index trajectories and incident PAD

A total of 9097 participants were included for further trajectory analysis (Fig. 1). Four discrete trajectories with stable TyG indexes at various levels from visit 1 to visit 5 were identified (Fig. 5): low (n = 2019, 22.2%), moderate (n = 3934, 43.2%), high (n = 2499, 27.5%), and very high (n = 645, 7.1%) TyG index trajectory groups. The median (interquartile range) changes in TyG index level during the visits for these trajectory groups were 0.06 (− 0.02 to 0.14) in the low trajectory group, 0.05 (− 0.04 to 0.14) in the moderate trajectory group, 0.03 (− 0.07 to 0.14) in the high trajectory group, 0.05 (− 0.11 to 0.18) in the very high trajectory group (Additional file 1: Tables S3). As shown in Fig. 6, the rates of incident PAD were 8.1%, 9.9%, 14.1%, 21.2%, in the low, moderate, high and very high. TyG index trajectory groups, respectively (P < 0.001). Multivariate logistic regression analyses identified those with TyG index trajectory at high and very high levels as having an even greater risk of incident PAD in 3 different adjusted models (all P < 0.05, Table 3). In the fully adjusted model, compared with those with a low trajectory at a low level, the odds ratios (95% CIs) for associations of those participants with TyG index trajectories at the moderate, high, and very high levels with the risk of incident PAD were 1.063 (95% CI, 0.871–1.298; P = 0.547), 1.404 (95% CI, 1.132–1.740; P = 0.002), and 1.742 (95% CI, 1.294–2.344; P < 0.001), respectively (model 3 in Table 3).

Fig. 5
figure5

Trajectories by TyG index in the Atherosclerosis Risk in Communities Study

Fig. 6
figure6

Prevalence of incident PAD across the TyG index trajectory groups

Table 3 Risk of incident PAD for various levels of TyG index trajectory groups

Discussion

In this large-scale, community-based prospective cohort of middle-aged adults, we show for the first time that higher levels of TyG index are significantly associated with an increased risk of PAD over a median follow-up of 23 years. Furthermore, we identify that the 4 distinct trajectories of TyG index confer different risk of PAD, and a two-decade trajectory with elevated TyG index carries a greater risk of future incident PAD. These findings suggest a potential role for long-lasting high level of IR in the pathogenesis of PAD.

IR has been considered as an important risk factor for cardiovascular disease [27, 28]. Methods to directly measure IR are invasive, complex, and costly [10, 11]. Therefore, a number of surrogate markers of IR have been proposed and compared with the gold standard of the hyperinsulinemic-euglycemic clamp [29]. Homeostatic model assessment of IR (HOMA-IR), which is calculated by fasting insulin and glucose, is commonly used for assessing IR. However, the insulin concentrations are not routinely measured in clinical settings. As for the ARIC cohort, there were many missing insulin values due to low detection rate and the absence of insulin measurement at visit 2 and visit 3. Therefore, in this large-scale, community-based prospective cohort study, we utilized the TyG index as a biomarker of IR. The TyG index has been proved to be highly correlated with the euglycemic-hyperinsulinemic clamp test [10], and thus has a validity similar to HOMA-IR [11]. The immense advantage of using such a simple method of IR identification is obviously that it is easily accessible in any clinical settings, making our findings immediately usable by clinicians.

Among the multiple pathological consequences of atherosclerosis, PAD has generally been paid far less attention than CHD or stroke. Based on the updated estimates of PAD prevalence at global regional levels in 2015, 236.62 million (5.56%) people aged 25 years and older had PAD, among whom 73% were in low-income and middle-income countries [30]. However, only about 10% of patients with PAD demonstrate the typical symptomatology of intermittent claudication; and the majority of patients with PAD are thus asymptomatic and underdiagnosed [31]. Therefore, it is of great importance to regularly measure markers of risk for PAD and take preventive measures at an earlier clinical stage. In addition to age, significant atherosclerotic risk factors for PAD include cigarette smoking, dyslipidemia, and diabetes [32]. It has been well-established that IR and coexisting hyperinsulinemia are implicated in the development of dyslipidemia, hypertension, hypercoagulability, and atherosclerosis [33, 34]. However, there is a paucity of prospective data regarding the association between IR assessed by HOMA-IR and PAD [35, 36]. A cross-sectional study of 3242 adults from data in the National Health and Nutrition Examination Survey has identified a positive association between IR and PAD [35]. Only one community-based longitudinal study, which enrolled 4208 participants over the age of 65 years in the Cardiovascular Health Study, showed that IR was associated with a higher risk of clinical PAD [36]. In line with previous studies, our study of a larger sample size showed that the metabolic risk factors such as hypertension, diabetes, and hyperlipidemia, were more obvious among participants of higher quartiles of TyG index. Meanwhile, individuals with the highest quartile of baseline TyG index had a 2.23-fold higher risk for developing PAD than those with the lowest quartile. The association remained statistically significant after adjusting for all the aforementioned PAD risk factors. These findings suggested that the clinical management of TyG index may bring additional effect on PAD development even under vigorous control of traditional risk factors. Further studies are needed to unravel this aspect. Most previous studies based on the TyG index measured at a single time point [12, 13, 16,17,18], which may not reflect long-term exposure, for the TyG index levels may vary over time. Therefore, measurements of long-term trajectories of TyG index provide more reliable and robust results. Our study is the first, to our knowledge, to investigate the impact of long-lasting IR at various levels assessed by TyG index on future PAD incidence. We highlighted the fact that within the ARIC population there were heterogeneous patterns of trends in TyG index, which cannot be fully identified by the baseline TyG index levels and that such risk groups may change during follow-up. Different from the initial TyG index levels, long-term TyG index trajectory actually reflects the chronic impact of TyG index on incident PAD. Our results further suggest that those trajectory groups with long-term high and very high TyG index levels beginning in midlife are at a greater risk of incident PAD over 20 years after adjustment for traditional PAD risk factors. The participant’s electronic medical record allows a rapid integration of data across multiple time points. Thus, longitudinal measurements and recording of TyG indexes to identify TyG index trajectories are feasible and it represents an added value to the baseline levels to plan and monitor participants’ follow-up. In clinical practice, we can graph trends in TyG index to identify high-risk individuals who behave similarly to those with TyG index trajectories at high and very high levels (both mean TyG index > 8.6) that were observed in the present analysis. Such population may benefit from earlier and more frequent screening for PAD and aggressive risk factor management (e.g., blood pressure control, smoking cessation, maintaining metabolic health, etc.).

Our findings at least in part support the important role of impaired systemic glucolipid metabolism in the pathophysiology of PAD [37]; higher TyG index is a surrogate marker of IR, which is responsible for chronic hyperglycemia and also dyslipidemia including high plasma TG [38, 39]. However, many studies have indicated that IR could promote atherosclerosis not only through mechanisms that involve systemic factors, such as dyslipidemia, hypertension, and a proinflammatory state, but also through the effect of perturbed insulin signaling at the level of the intimal cells [34, 40]. Therefore, further studies are warranted to elucidate the precise mechanism for the observed association.

Study limitations

Several limitations of this investigation are worth noting. The study included only whites and blacks aged 45–64 years at baseline, results may differ outside this age range and in other ethnicities. Due to the missing records of insulin levels in the ARIC study, we cannot compare trajectories of TyG index with HOMA-IR for predicting incident PAD. Inflammation-induced endothelial dysfunction was demonstrated to be at least one possible biological pathway between IR and atherosclerosis [41]. High-sensitivity C-reactive protein was not collected at baseline in the ARIC database and thus we unfortunately could not study the possible association with this sensitive marker of systemic inflammation [8, 42, 44]. However, we could analyze other inflammation biomarkers-albeit less sensitive-(fibrinogen and white blood cell count) as well as coagulation factors (von Willebrand factor, factor VIII activity) that were available from the ARIC database [42], as studies also support the role of these hemostatic factors in PAD [37, 43]. Moreover, because of the nature of any observational studies, we cannot exclude the possibility of residual confounders despite our careful adjustment for the well-known and suspected risk factors.

Conclusions

Higher TyG index is associated with incident PAD, which is independent of other traditional atherosclerotic risk factors; suggesting that IR is actually involved in the pathogenesis of PAD. Trajectories denoting long-term exposure to high IR assessed by TyG index provide additional information about the cumulative burden of risk for future PAD. The results support the contribution of higher TyG index to the development of PAD and have implications on its prevention and treatment.

Availability of data and materials

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Code availability

Data described in the manuscript, code book, and analytic code will be made available upon request pending application and approval. But we will make the data used in the manuscript, code book, and analytic code available to editors upon request either before or after publication for checking.

Abbreviations

ABI:

Ankle‐brachial index

ARIC:

Atherosclerosis Risk in Communities Study

BMI:

Body mass index

CHD:

Coronary heart disease

CLI:

Critical limb ischemia

DBP:

Diastolic blood pressure

eGFR:

Estimated glomerular filtration rate

HDL-C:

High-density lipoprotein cholesterol

HOMA-IR:

Homeostatic model assessment of insulin resistance

IR:

Insulin resistance

LDL-C:

Low-density lipoprotein cholesterol

PAD:

Peripheral artery disease

SBP:

Systolic blood pressure

TC:

Total cholesterol

TG:

Triglycerides

TyG:

Triglyceride-glucose

References

  1. 1.

    Fowkes FG, Rudan D, Rudan I, Aboyans V, Denenberg JO, McDermott MM, Norman PE, Sampson UK, Williams LJ, Mensah GA, et al. Comparison of global estimates of prevalence and risk factors for peripheral artery disease in 2000 and 2010: a systematic review and analysis. Lancet. 2013;382(9901):1329–40.

    PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Criqui MH, Aboyans V. Epidemiology of peripheral artery disease. Circ Res. 2015;116(9):1509–26.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  3. 3.

    Gerhard-Herman MD, Gornik HL, Barrett C, Barshes NR, Corriere MA, Drachman DE, Fleisher LA, Fowkes FG, Hamburg NM, Kinlay S, et al. 2016 AHA/ACC Guideline on the management of patients with lower extremity peripheral artery disease: executive summary: a report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines. Circulation. 2017;135(12):e686–725.

    PubMed  PubMed Central  Google Scholar 

  4. 4.

    Faerch K, Vaag A, Holst JJ, Hansen T, Jørgensen T, Borch-Johnsen K. Natural history of insulin sensitivity and insulin secretion in the progression from normal glucose tolerance to impaired fasting glycemia and impaired glucose tolerance: the Inter99 study. Diabetes Care. 2009;32(3):439–44.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    DeFronzo RA. Insulin resistance, lipotoxicity, type 2 diabetes and atherosclerosis: the missing links. The Claude Bernard Lecture 2009. Diabetologia. 2010;53(7):1270–87.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Brevetti G, Giugliano G, Brevetti L, Hiatt WR. Inflammation in peripheral artery disease. Circulation. 2010;122(18):1862–75.

    PubMed  Article  PubMed Central  Google Scholar 

  7. 7.

    Dann R, Hadi T, Montenont E, Boytard L, Alebrahim D, Feinstein J, Allen N, Simon R, Barone K, Uryu K, et al. Platelet-derived MRP-14 induces monocyte activation in patients with symptomatic peripheral artery disease. J Am Coll Cardiol. 2018;71(1):53–65.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    Ding N, Yang C, Ballew SH, Kalbaugh CA, McEvoy JW, Salameh M, Aguilar D, Hoogeveen RC, Nambi V, Selvin E, et al. Fibrosis and inflammatory markers and long-term risk of peripheral artery disease: the ARIC Study. Arterioscler Thromb Vasc Biol. 2020;40(9):2322–31.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  9. 9.

    Kullo IJ, Leeper NJ. The genetic basis of peripheral arterial disease: current knowledge, challenges, and future directions. Circ Res. 2015;116(9):1551–60.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Guerrero-Romero F, Simental-Mendía LE, González-Ortiz M, Martínez-Abundis E, Ramos-Zavala MG, Hernández-González SO, Jacques-Camarena O, Rodríguez-Morán M. The product of triglycerides and glucose, a simple measure of insulin sensitivity. Comparison with the euglycemic-hyperinsulinemic clamp. J Clin Endocrinol Metab. 2010;95(7):3347–51.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  11. 11.

    Simental-Mendía LE, Rodríguez-Morán M, Guerrero-Romero F. The product of fasting glucose and triglycerides as surrogate for identifying insulin resistance in apparently healthy subjects. Metab Syndr Relat Disord. 2008;6(4):299–304.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  12. 12.

    Ma X, Dong L, Shao Q, Cheng Y, Lv S, Sun Y, Shen H, Wang Z, Zhou Y, Liu X. Triglyceride glucose index for predicting cardiovascular outcomes after percutaneous coronary intervention in patients with type 2 diabetes mellitus and acute coronary syndrome. Cardiovasc Diabetol. 2020;19(1):31.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Barzegar N, Tohidi M, Hasheminia M, Azizi F, Hadaegh F. The impact of triglyceride-glucose index on incident cardiovascular events during 16 years of follow-up: Tehran Lipid and Glucose Study. Cardiovasc Diabetol. 2020;29(19):155.

    Article  CAS  Google Scholar 

  14. 14.

    Irace C, Carallo C, Scavelli FB, De Franceschi MS, Esposito T, Tripolino C, Gnasso A. Markers of insulin resistance and carotid atherosclerosis. A comparison of the homeostasis model assessment and triglyceride glucose index. Int J Clin Pract. 2013;67(7):665–72.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    Park K, Ahn CW, Lee SB, Kang S, Nam JS, Lee BK, Kim JH, Park JS. Elevated TyG index predicts progression of coronary artery calcification. Diabetes Care. 2019;42(8):1569–73.

    CAS  Article  Google Scholar 

  16. 16.

    Chiu TH, Tsai HJ, Chiou HC, Wu PY, Huang JC, Chen SC. A high triglyceride-glucose index is associated with left ventricular dysfunction and atherosclerosis. Int J Med Sci. 2021;18(4):1051–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Duran Karaduman B, Ayhan H, Keleş T, Bozkurt E. The triglyceride-glucose index predicts peripheral artery disease complexity. Turk J Med Sci. 2020;50:1217–22.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  18. 18.

    Zhao S, Yu S, Chi C, Fan X, Tang J, Ji H, Teliewubai J, Zhang Y, Xu Y. Association between macro- and microvascular damage and the triglyceride glucose index in community-dwelling elderly individuals: the Northern Shanghai Study. Cardiovasc Diabetol. 2019;18(1):95.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  19. 19.

    The Atherosclerosis Risk in Communities (ARIC) Study: design and objectives. The ARIC investigators. Am J Epidemiol. 1989;129(4):687–702.

  20. 20.

    Baecke JA, Burema J, Frijters JE. A short questionnaire for the measurement of habitual physical activity in epidemiological studies. Am J Clin Nutr. 1982;36(5):936–42.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Papp AC, Hatzakis H, Bracey A, Wu KK. ARIC hemostasis study—I. Development of a blood collection and processing system suitable for multicenter hemostatic studies. Thromb Haemost. 1989;61(1):15–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. 1972;18(6):499–502.

    CAS  Article  Google Scholar 

  23. 23.

    Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF 3rd, Feldman HI, Kusek JW, Eggers P, Van Lente F, Greene T, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–12.

    PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Wattanakit K, Folsom AR, Selvin E, Coresh J, Hirsch AT, Weatherley BD. Kidney function and risk of peripheral arterial disease: results from the Atherosclerosis Risk in Communities (ARIC) Study. J Am Soc Nephrol. 2007;18(2):629–36.

    PubMed  Article  PubMed Central  Google Scholar 

  25. 25.

    Liang J, Zhang H, Sun X, Liao L, Li X, Hu X, Du J, Zhuang X, Liao X. Association between calf girth and peripheral artery disease in the Atherosclerosis Risk in Communities Study. J Cardiol. 2020;76(3):273–9.

    PubMed  Article  PubMed Central  Google Scholar 

  26. 26.

    Nagin DS, Odgers CL. Group-based trajectory modeling in clinical research. Annu Rev Clin Psychol. 2010;6:109–38.

    PubMed  Article  PubMed Central  Google Scholar 

  27. 27.

    Laakso M, Kuusisto J. Insulin resistance and hyperglycaemia in cardiovascular disease development. Nat Rev Endocrinol. 2014;10(5):293–302.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  28. 28.

    Adeva-Andany MM, Martínez-Rodríguez J, González-Lucán M, Fernández-Fernández C, Castro-Quintela E. Insulin resistance is a cardiovascular risk factor in humans. Diabetes Metab Syndr. 2019;13(2):1449–55.

    PubMed  Article  PubMed Central  Google Scholar 

  29. 29.

    Bastard JP, Lavoie ME, Messier V, Prud’homme D, Rabasa-Lhoret R. Evaluation of two new surrogate indices including parameters not using insulin to assess insulin sensitivity/resistance in non-diabetic postmenopausal women: a MONET group study. Diabetes Metab. 2012;38(3):258–63.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  30. 30.

    Song P, Rudan D, Zhu Y, Fowkes FJI, Rahimi K, Fowkes FGR, Rudan I. Global, regional, and national prevalence and risk factors for peripheral artery disease in 2015: an updated systematic review and analysis. Lancet Glob Health. 2019;7(8):e1020–30.

    PubMed  Article  PubMed Central  Google Scholar 

  31. 31.

    Hirsch AT, Criqui MH, Treat-Jacobson D, Regensteiner JG, Creager MA, Olin JW, Krook SH, Hunninghake DB, Comerota AJ, Walsh ME, et al. Peripheral arterial disease detection, awareness, and treatment in primary care. JAMA. 2001;286(11):1317–24.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  32. 32.

    Barnes JA, Eid MA, Creager MA, Goodney PP. Epidemiology and risk of amputation in patients with diabetes mellitus and peripheral artery disease. Arterioscler Thromb Vasc Biol. 2020;40(8):1808–17.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  33. 33.

    Wang T, Zhao Z, Xu Y, Qi L, Xu M, Lu J, Li M, Chen Y, Dai M, Zhao W, et al. Insulin resistance and β-cell dysfunction in relation to cardiometabolic risk patterns. J Clin Endocrinol Metab. 2018;103(6):2207–15.

    PubMed  Article  PubMed Central  Google Scholar 

  34. 34.

    Bornfeldt KE, Tabas I. Insulin resistance, hyperglycemia, and atherosclerosis. Cell Metab. 2011;14(5):575–85.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Pande RL, Perlstein TS, Beckman JA, Creager MA. Association of insulin resistance and inflammation with peripheral arterial disease: the National Health and Nutrition Examination Survey, 1999 to 2004. Circulation. 2008;118(1):33–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Britton KA, Mukamal KJ, Ix JH, Siscovick DS, Newman AB, de Boer IH, Thacker EL, Biggs ML, Gaziano JM, Djoussé L. Insulin resistance and incident peripheral artery disease in the Cardiovascular Health Study. Vasc Med. 2012;17(2):85–93.

    PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Garg PK, Biggs ML, Carnethon M, Ix JH, Criqui MH, Britton KA, Djoussé L, Sutton-Tyrrell K, Newman AB, Cushman M, et al. Metabolic syndrome and risk of incident peripheral artery disease: the cardiovascular health study. Hypertension. 2014;63(2):413–9.

    CAS  PubMed  Article  Google Scholar 

  38. 38.

    Ferrannini E, Barrett EJ, Bevilacqua S, DeFronzo RA. Effect of fatty acids on glucose production and utilization in man. J Clin Invest. 1983;72(5):1737–47.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    Ormazabal V, Nair S, Elfeky O, Aguayo C, Salomon C, Zuñiga FA. Association between insulin resistance and the development of cardiovascular disease. Cardiovasc Diabetol. 2018;17(1):122.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Giacco F, Brownlee M. Oxidative stress and diabetic complications. Circ Res. 2010;107(9):1058–70.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Brevetti G, Schiano V, Chiariello M. Endothelial dysfunction: a key to the pathophysiology and natural history of peripheral arterial disease? Atherosclerosis. 2008;197(1):1–11.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  42. 42.

    Walker KA, Gottesman RF, Wu A, Knopman DS, Gross AL, Mosley TH Jr, Selvin E, Windham BG. Systemic inflammation during midlife and cognitive change over 20 years: The ARIC Study. Neurology. 2019;92(11):e1256–67.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Small AM, Huffman JE, Klarin D, Sabater-Lleal M, Lynch JA, Assimes TL, Sun YV, Miller D, Freiberg MS, Morrison AC, et al. Mendelian randomization analysis of hemostatic factors and their contribution to peripheral artery disease-brief report. Arterioscler Thromb Vasc Biol. 2021;41(1):380–6.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank the investigators, the staffs and participants of the ARIC Study for their highly valued contributions.

Funding

This work was supported by grants from National Natural Science Foundation of China [81900379, 81870315, 81970683]. The ARIC Study has been funded in whole or in part with Federal funds from the National Heart, Lung, and Blood Institute, National Institutes of Health, of the United States. Department of Health and Human Services, under Contract Nos. (HHSN268201700001I, HHSN268201700003I, HHSN268201700005I, HHSN268201700004I, HHSN2682017000021).

Author information

Affiliations

Authors

Contributions

JWG, QYH and PML contributed to the study concept, and drafted the manuscript. JWG and QYH has full access to the data and performed the analyses; JWG and QYH assisted in data interpretation; MG, KZ and XZL helped in the data methods and presentation; JWG, DAV, SLZ and PML critically revised the manuscript. JFW and PML supervised the study analyses. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Shao-Ling Zhang or Pin-Ming Liu.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the institutional review boards at all field centers of ARIC study, and informed consent was obtained from all participants.

Consent for publication

The consent to publish was obtained from all participants in this study.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. Group-based trajectory model fit summary (N=9097). Table S2. Risk of incident PAD for baseline TyG index among participants without any lipid- or glucose-lowering medication. Table S3. TyG index at examination visits by trajectory groups of TyG index.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gao, JW., Hao, QY., Gao, M. et al. Triglyceride-glucose index in the development of peripheral artery disease: findings from the Atherosclerosis Risk in Communities (ARIC) Study. Cardiovasc Diabetol 20, 126 (2021). https://doi.org/10.1186/s12933-021-01319-1

Download citation

Keywords

  • Triglyceride-glucose index
  • Insulin resistance
  • Peripheral artery disease
  • Risk factors
  • Cardiovascular disease