Skip to main content

Changes in the triglyceride glucose-body mass index estimate the risk of stroke in middle-aged and older Chinese adults: a nationwide prospective cohort study

Abstract

Background

Stroke was reported to be highly correlated with the triglyceride glucose-body mass index (TyG-BMI). Nevertheless, literature exploring the association between changes in the TyG-BMI and stroke incidence is scant, with most studies focusing on individual values of the TyG-BMI. We aimed to investigate whether changes in the TyG-BMI were associated with stroke incidence.

Methods

Data were obtained from the China Health and Retirement Longitudinal Study (CHARLS), which is an ongoing nationally representative prospective cohort study. The exposures were changes in the TyG-BMI and cumulative TyG-BMI from 2012 to 2015. Changes in the TyG-BMI were classified using K-means clustering analysis, and the cumulative TyG-BMI was calculated as follows: (TyG-BMI2012 + TyG-BMI2015)/2 × time (2015–2012). Logistic regressions were used to determine the association between different TyG-BMI change classes and stroke incidence. Meanwhile, restricted cubic spline regression was applied to examine the potential nonlinear association of the cumulative TyG-BMI and stroke incidence. Weighted quantile sum regression was used to provide a comprehensive explanation of the TyG-BMI by calculating the weights of FBG, triglyceride-glucose (TG), and BMI.

Results

Of the 4583 participants (mean [SD] age at baseline, 58.68 [9.51] years), 2026 (44.9%) were men. During the 3 years of follow-up, 277 (6.0%) incident stroke cases were identified. After adjusting for potential confounders, compared to the participants with a consistently low TyG-BMI, the OR for a moderate TyG-BMI with a slow rising trend was 1.01 (95% CI 0.65–1.57), the OR for a high TyG-BMI with a slow rising trend was 1.62 (95% CI 1.11–2.32), and the OR for the highest TyG-BMI with a slow declining trend was 1.71 (95% CI 1.01–2.89). The association between the cumulative TyG-BMI and stroke risk was nonlinear (Passociation = 0.017; Pnonlinearity = 0.012). TG emerged as the primary contributor when the weights were assigned to the constituent elements of the TyG-BMI (weight2012 = 0.466; weight2015 = 0.530).

Conclusions

Substantial changes in the TyG-BMI are independently associated with the risk of stroke in middle-aged and older adults. Monitoring long-term changes in the TyG-BMI may assist with the early identification of individuals at high risk of stroke.

Introduction

Stroke, the foremost cause of mortality and disability in numerous nations, poses a progressively escalating burden on the global healthcare system [1]. The Global Burden of Disease 2013 Study reported a substantial number of stroke survivors, with approximately 25.7 million individuals alive, along with 6.5 million stroke-related fatalities and 113 million stroke-related disability-adjusted life-years observed worldwide [2]. While the age-standardized mortality rate of stroke has exhibited a decline on a global scale, the incidence and prevalence of stroke in China continue to exhibit an alarming increase [2,3,4]. Consequently, there exists an immediate imperative to develop cost-effective and reproducible indicators that can enhance the stratification of stroke risk.

Insulin resistance (IR), recognized as a novel risk factor for stroke, is considered an early indication of type-2 diabetes, extending beyond patients with diabetes to include nondiabetic individuals as well [5]. Various methods exist to evaluate IR, with the hyperinsulinaemic glucose clamp (HEC) being the established reference standard [6]. However, HEC necessitates intravenous administration of glucose and insulin, accompanied by multiple blood samples [7]. This intricate and costly procedure is not widely employed in clinical practice. In contrast, the Homeostasis Model Assessment of IR (HOMA-IR) has gained widespread usage and demonstrated effectiveness in predicting cardiovascular and cerebrovascular diseases [7,8,9,10]. Nonetheless, HOMA-IR requires measurement of fasting insulin levels, which possesses limited clinical practicality. As a feasible alternative for assessing IR, the triglyceride-glucose (TyG) index can be effortlessly derived from routine clinical laboratory tests and is reported to be associated with stroke occurrence and recurrence. Furthermore, recent investigations have revealed that the TyG index outperforms HOMA-IR in predicting stroke [11, 12].

Recently, a study introduced a novel index called the triglyceride glucose-body mass index (TyG-BMI), which is computed as ln [TG (mg/dl) × FBG (mg/dl)/2] × BMI (kg/m2), with the objective of using it as a potential and straightforwards marker for IR [13]. TyG-BMI incorporates BMI into its calculation. This modification may allow TyG-BMI to better account for the influence of obesity on IR, and the combination of obesity and TyG can potentially identify IR more strongly than other surrogate markers, since obesity is a well-established risk factor for IR [14]. Several studies have also shown that TyG-BMI has a better predictive performance than TyG in metabolic diseases or cardiovascular disease [15,16,17,18]. However, only a cross-sectional study conducted on two population-based samples in China demonstrated a positive association between a high TyG-BMI and ischaemic stroke [19]. However, the exploration of the association between changes in TyG-BMI and stroke incidence has been infrequently reported in the literature, with most studies focusing on individual TyG-BMI values. Despite the valuable insights provided by earlier studies, there remains a notable gap in our understanding of how changes in TyG-BMI may relate to stroke risk, especially in a nationwide context. In addition, the relative contributions of TG, FBG, and BMI to stroke incidence is not clear. Our study aims to fill this gap by exploring this relationship in greater depth, utilizing data sourced from the "China Health and Retirement Longitudinal Study (CHARLS)", an ongoing nationwide cohort study designed to represent the population.

Methods

Study population

This study utilizes a secondary analysis of data from the CHARLS, an ongoing nationwide cohort study designed to represent the population [20]. The study design has been previously described [20]. In summary, a total of 17,708 participants residing in 10,257 households were selected using a multistage stratified probability-proportional-to-size sampling technique. These participants were recruited from 150 counties or districts and 450 villages within 28 provinces in China, spanning the period from June 2011 to March 2012. A standardized questionnaire was administered to collect information on sociodemographic and lifestyle factors, as well as health-related data. The baseline survey (Wave 1) achieved a response rate of 80.5%, and subsequent follow-up assessments were conducted every 2 years, with Wave 2 in 2013, Wave 3 in 2015 and Wave 4 in 2018. Blood samples were also collected at baseline and Wave 3. For this analysis, participants had to be aged 45 years and older, and complete data on fasting blood glucose (FBG), triglycerides (TG) and body mass index (BMI) were needed. People were excluded if they had a stroke before 2015.

The CHARLS study was approved by the institutional review board of Peking University. Written informed consent was obtained from all participants. This study was conducted following the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guidelines [21].

Assessment of the change in TyG-BMI

The exposure of this study was the change in TyG-BMI values between 2012 and 2015. The TyG-BMI was calculated by the formula ln[TG (mg/dl) × FBG (mg/dl)/2] × BMI (kg/m2) [13]. We calculated the cumulative TyG-BMI with reference to the cumulative TyG change formula [22]: (TyG-BMI2012 + TyG-BMI2015)/2 × time (2015 − 2012). Height and weight were measured by a trained nurse. BMI was calculated as weight in kilograms divided by height in metres squared.

Ascertainment of incident stroke events

Stroke was the main outcome of this study. In accordance with previous studies [23, 24], stroke events were assessed by individuals who self-reported “yes” to the question of “Have you been diagnosed with stroke by a doctor?” or selected specific answers to questions regarding the treatment of stroke “by Chinese Traditional Medicine/Western Modern Medicine/Physical Therapy/Acupuncture and Moxibustion/Occupational Therapy/None of the Above” were regarded as people with stroke. The date of stroke diagnosis was recorded as being between the date of the last interview and that of the interview reporting an incident stroke [23, 24].

Covariates

At baseline (Wave 1), trained interviewers collected information on sociodemographic status and health-related factors using a structured questionnaire, including age, sex, living residence, marital status, and educational level. Educational level was classified as no formal education, primary school, middle or high school, and college or above. Health-related factors included self-reported smoking and drinking status (never, former, or current), self-reported physician-diagnosed medical conditions (diabetes, hypertension, heart problems, kidney disease, and dyslipidaemia), and use of medications for diabetes, hypertension, and dyslipidaemia. Laboratory examination contained total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), estimated glomerular filtration ratio (eGFR), and glycosylated haemoglobin, type A1c (HbA1c) [25]. Marital status was classified into 2 groups: married and other marital status (never married, separated, divorced, and widowed). Diabetes was defined as fasting plasma glucose ≥ 126 mg/dl (to convert to millimoles per litre, multiplied by 0.0555), current use of antidiabetic medication, or self-reported history of diabetes. Hypertension was defined as systolic blood pressure ≥ 140 mmHg, diastolic blood pressure ≥ 90 mmHg, current use of antihypertensive medication, or self-reported history of hypertension. Dyslipidaemia was defined as total cholesterol ≥ 240 mg/dl (to convert to millimoles per litre, multiplied by 0.0259), triglycerides ≥ 150 mg/dl, low-density lipoprotein cholesterol ≥ 160 mg/dl, high-density lipoprotein cholesterol < 40 mg/dl, current use of lipid-lowering medication, or self-reported history of dyslipidaemia.

Statistical analyses

Data were analysed using various statistical methods. For normally distributed continuous variables, means and standard deviations (SDs) were reported. Categorical variables were described in terms of frequency and percentage. The χ2 test and analysis of variance were utilized to compare differences in the baseline characteristics between the different groups, as appropriate. Approximately 4.3% (197/4583) of the total data items were found to be missing. These missing values were assumed to be missing at random and were addressed through the multiple imputation of chained equations method using the baseline characteristics. To account for the missing data, five imputed datasets were created. The results were then pooled using R statistical software along with the mice package.

We utilized an unsupervised machine learning technique called K-means with Euclidean distance to group patients based on their TyG-BMI measurements in 2012 and 2015. We opted for the k-means algorithm due to its computational efficiency and ability to generate easy-to-understand visualizations of the data points [26]. The K-means algorithm, a centroid-based clustering approach, divides a dataset into K clusters by minimizing the sum of squared distances within each cluster [26, 27]. To execute the process, we followed a three-step procedure: first, we specified the desired number of clusters; then, we randomly selected k patients as the initial cluster centres; next, we assigned each patient to the nearest centroid and sequentially updated the cluster centroids [27]. This iterative process continued until the total within-cluster sum of squares was minimized, and each patient was assigned to a specific cluster based on their distance from the centres, as determined by the Euclidean distance [26]. To determine the appropriate number of clusters, we visually analysed the reduction in the sum of squared distances resulting from varying the number of clusters. We presented a visual representation of the resulting clustering. Importantly, throughout the entire clustering process, the algorithm remained unaware of the outcome variables, ensuring unbiased analysis.

During a median follow-up of 36.7 months (range: 29.6–41.8 months) from baseline to Wave 3, the data set of transition of the TyG-BMI was analysed and classified into 4 classes using K-means clustering (Fig. 1A), and the paired t test was employed to assess the changes within each class: for Class 1 (n = 1273), the TyG-BMI ranged from 163.42 ± 15.13 in 2012 to 163.36 ± 15.45 in 2015 (P = 0.914), and the mean (SD) cumulative TyG-BMI was 490.18 ± 38.64, representing a consistently low TyG-BMI; for Class 2 (n = 1664), the TyG-BMI ranged from 197.39 ± 14.40 in 2012 to 200.17 ± 14.45 in 2015 (P < 0.001), and the mean (SD) cumulative TyG-BMI was 596.34 ± 30.55, representing a moderate TyG-BMI with a slow rising trend; for Class 3 (n = 1171), the TyG-BMI ranged from 233.00 ± 17.59 in 2012 to 236.43 ± 16.93 in 2015 (P < 0.001), and the mean (SD) cumulative TyG-BMI was 704.14 ± 34.37, representing a high TyG-BMI with a slow rising trend; for Class 4 (n = 475), the TyG-BMI ranged from 283.42 ± 29.91 in 2012 to 279.56 ± 27.16 in 2015 (P = 0.037), and the mean (SD) cumulative TyG-BMI was 844.48 ± 62.90, representing the highest TyG-BMI with a slow declining trend (Fig. 1B). The distribution of TyG-BMI according to the change in TyG-BMI classes is shown in Fig. 1C, D. The normal distribution of TyG-BMI within each class is observable, and there exists a statistically significant disparity in the mean TyG-BMI values among these classes.

Fig. 1
figure 1

Clustering of the change in the TyG-BMI from 2012 to 2015. A Four clusters were found using the K-means method with Euclidean distance: the x- and y-axes are principal components of the change in the TyG-BMI; B data visualization for the classes of the change in the TyG-BMI: in Class 1 (n = 1273), the TyG-BMI ranged from 163.42 ± 15.13 in 2012 to 163.36 ± 15.45 in 2015 (P = 0.914); in Class 2 (n = 1664), the TyG-BMI ranged from 197.39 ± 14.40 in 2012 to 200.17 ± 14.45 in 2015 (P < 0.001); in Class 3 (n = 1171), the TyG-BMI ranged from 233.00 ± 17.59 in 2012 to 236.43 ± 16.93 in 2015 (P < 0.001); and in Class 4 (n = 475), the TyG-BMI ranged from 283.42 ± 29.91 in 2012 to 279.56 ± 27.16 in 2015 (P = 0.037); C, D Distribution for the TyG-BMI at 2012 or 2015: the normal distribution of the TyG-BMI within each class is observable, and there exists a statistically significant disparity in the mean TyG-BMI values among these classes. BMI: body mass index; PC: principal component; TyG: triglyceride-glucose

Upon establishing classes of the change in TyG-BMI, we examined the association between different classes with changes in TyG-BMI and stroke events, and binary logistic models were used to calculate odds ratios (ORs) with 95% CIs. Four models were estimated: Model 1 adjusted for age and sex; Model 2 adjusted for age, sex, marital status, residence, educational level, smoking status, and drinking status; Model 3 adjusted for the variables in Model 2 and history of hypertension, diabetes, heart disease, dyslipidaemia, kidney disease, medication use for hypertension, medication use for diabetes, medication use for dyslipidaemia, systolic blood pressure, and diastolic blood pressure; and Model 4 adjusted for the variables in Model 3 and total cholesterol, HDL-C, LDL-C, HbA1c, and the eGFR.

To examine the association between the cumulative TyG-BMI and stroke events, the cumulative TyG-BMI was split into quartiles and then included in binary logistic models with the first quartile as the reference group. We searched for a linear trend by modelling the median value of each quantile to test ordered relations across quantiles of cumulative TyG-BMI. In addition, we explored the potential nonlinear association using a restricted cubic spline (RCS) regression model, and the model was conducted with 4 knots at the 5th, 35th, 65th, and 95th percentiles of cumulative TyG-BMI (reference is the 5th percentile). We further applied a two-piecewise linear regression model to examine the threshold effect of the cumulative TyG-BMI on stroke using a smoothing function [28, 29]. Subgroup analyses were conducted to examine whether the potential association of the change in TyG-BMI and cumulative TyG-BMI with stroke was moderated by the following demographic and clinical characteristics: age, sex, marital status, residence, educational level, smoking status, drinking status, history of hypertension, diabetes, heart disease, dyslipidaemia, and kidney disease. P values for interaction were evaluated using interaction terms and likelihood ratio tests.

To evaluate the performance of the TyG-BMI in predicting stroke, a receiver-operating characteristic (ROC) curve analysis was conducted, the area under the ROC curves (AUCs) were calculated, and the curves were generated with a smooth kernel density. To evaluate the extent to which TyG-BMI improved the prediction performance over TyG, Delong’s test was performed. TyG-BMI was derived using a mathematical formula incorporating TG, FBG, and BMI variables. To provide a comprehensive explanation of the formula, we employed the weighted quantile sum (WQS) regression model, employing bootstrap resampling methods for 1000 iterations. The WQS model facilitated the determination of weights assigned to FBG, TG, and BMI, quantifying their respective contributions to the overall effect. These weights were constrained within the range of 0 to 1, with a cumulative sum of 1 [30]. Higher weights indicated greater significance of the corresponding indicator in stroke prediction.

Five sensitivity analyses were conducted as follows: (1) repeating primary analyses using the complete data set (4373 participants) without multiple imputations; (2) repeating primary analyses excluding participants who had heart disease (538 participants) to account for loss to follow-up due to cardiovascular disease; (3) using the competing risk model to account for competing risks due to mortality (188 participant deaths); (4) using the Cox proportional hazards models to account for censored data; and (5) treating cumulative TyG-BMI as a continuous variable to examine the linear relationship between cumulative TyG-BMI and the risk of stroke without imposing predefined categories. Considering that the alteration in effect size per unit of TyG-BMI is small, we undertook the normalization of TyG-BMI to assess the effect in terms of a per-SD change in TyG-BMI. A two-sided P < 0.05 was considered to indicate statistical significance. All analyses were performed using R statistical software version 4.2.2 (R Foundation).

Results

Baseline characteristics of study participants

Of the 17,708 CHARLS participants at study baseline, we excluded 11,770 individuals who lacked FBG data at Waves 1 and 3. Additionally, 982 participants were excluded due to incomplete information on the TyG index and BMI at Waves 1 and 3. Furthermore, we excluded 136 participants younger than 45 years and 237 individuals who reported a history of stroke at Waves 1 and 3. Consequently, a total of 4583 participants met the inclusion criteria and were included in the subsequent analysis (Fig. 2). A comparison of baseline characteristics between participants included and those who were not included in the analysis is shown in Additional file 1: Table S1.

Fig. 2
figure 2

Flowchart of the study population. BMI: body mass index; FBG: fasting blood glucose; TyG: triglyceride-glucose

A total of 4583 adults were included in the analyses. The mean age at baseline was 58.68 ± 9.51 years; 2056 (44.9%) of the participants were men, and 2527 (55.1%) were women. Table 1 shows the characteristics of the participants. The mean TyG-BMI was 205.97 ± 40.60 in 2012 and 207.44 ± 40.06 in 2015, and the mean cumulative TyG-BMI was 620.12 ± 115.68.

Table 1 Baseline characteristics of 4583 participants according to the change in the TyG-BMI

When compared with Class 1, participants in the other classes were more likely to be older, female, and married; to reside in a rural setting; to have a middle or high school education, fewer current smokers and drinkers, and higher systolic and diastolic blood pressure; to have a higher prevalence of hypertension, diabetes, dyslipidaemia and heart disease; to have a history of medication use for hypertension, diabetes and dyslipidaemia; to have higher TC, LDL-C, and HbA1 levels; and to have lower HDL-C levels (Table 1).

Odds ratios for incident stroke

During a median follow-up of 36.5 months between 2015 and 2018, 277 participants experienced incident stroke, and the incidence rate of stroke was 6.0%. Table 2 shows the associations between the change in TyG-BMI and incident stroke events. After adjusting for potential confounders (in Model 4), compared with Class 1, the adjusted ORs (95% CIs) for incident stroke were 1.01 (0.65–1.57) for Class 2, 1.62 (1.11–2.32) for Class 3, and 1.71 (1.01–2.89) for Class 4. The results were found when modelling the cumulative TyG-BMI as quantiles (Table 2), and the baseline characteristics of the participants according to the quantiles of cumulative TyG-BMI are shown in Additional file 1: Table S2. After adjusting for confounders (in Model 4), when compared with Quartile 1, the adjusted ORs (95% CIs) for incident stroke were 1.66 (1.11–2.50) for Quartile 2, 1.41 (0.91–2.17) for Quartile 3, and 1.36 (0.85–2.18) for Quartile 4. Notably, no statistically significant differences were observed between Quartiles 3, 4, and 1 (P for trend = 0.563). Moreover, ROC analyses (Additional file 1: Fig. S1) demonstrated that baseline TyG-BMI had greater accuracy in predicting stroke compared with TyG (AUC, 0.62 vs. 0.57; P < 0.001).

Table 2 Associations of different classes of the TyG-BMI with stroke incidence

In the RCS model, the association between cumulative TyG-BMI and risk of incident stroke was nonlinear (for association, P = 0.017; for nonlinearity, P = 0.012) (Fig. 3). The risk of stroke increased with a cumulative TyG-BMI < 570 (OR per SD increased: 1.32, 95% CI 1.00–1.74). The risk decreased when the cumulative TyG-BMI was between 570 and 720 (OR per SD decreased: 0.78, 95% CI 0.64–0.95). The risk increased with a cumulative TyG-BMI > 720 (OR per SD increased: 1.15, 95% CI 0.90–1.47); however, the difference was not statistically significant (P = 0.251).

Fig. 3
figure 3

Nonlinear association between cumulative TyG-BMI and stroke. A Distribution for cumulative TyG-BMI from 2012 to 2015; B graphs show ORs for stroke. Data were fitted by a restricted cubic spline (RCS) logistics regression model, and the models were conducted with 4 knots at the 5th, 35th, 65th, and 95th percentiles of the cumulative TyG-BMI (reference is the 5th percentile). TyG-BMI ranged from 374.60 to 1032.59, due the RCS model handles outliers or extremes. Solid lines indicate ORs, and shadow shapes indicate 95% CIs. C Forest plot for the association between cumulative TyG-BMI and stroke. OR was evaluated by a 1‑SD increase in TyG-BMI. All models were adjusted for age, sex, marital status, residence, educational level, smoking status, drinking status, history of hypertension, diabetes, heart disease, dyslipidaemia, kidney disease, medication use for hypertension, medication use for diabetes, medication use for dyslipidaemia, systolic blood pressure, diastolic blood pressure, total cholesterol, HDL-C, LDL-C, HbA1c, and the eGFR. BMI: body mass index; CI: confidence interval; SD: standard deviation; OR: odds ratio; TyG: triglyceride-glucose

Subgroup analyses

Tables 3 and 4 show the association of the change in TyG-BMI and cumulative TyG-BMI with incident stroke events stratified by potential risk factors. After adjusting for potential confounders, no interaction was found between changes in TyG-BMI classes and subgroup variables. Marital status and educational level moderated the association of cumulative TyG-BMI with stroke.

Table 3 Associations of different classes of change in TyG-BMI with stroke incidence stratified by different factors
Table 4 Associations of the cumulative TyG-BMI with stroke stratified by different factors

WQS analyses

The WQS regression model was utilized to offer a thorough elucidation of the TyG-BMI. The weights assigned to the constituent elements of TyG-BMI, encompassing the cumulative influence of stroke while controlling for potential confounding variables, are presented in Fig. 4. Importantly, TG emerged as the primary contributor in both 2012 and 2015, with weights of 0.466 and 0.530, respectively.

Fig. 4
figure 4

Estimated weights assigned to TyG-BMI with the WQS model. Weights in a positive direction obtained when the effect parameter of the WQS model was constrained to the positive direction with 1 000 repeated holdout validations for TyG-BMI in 2012 (A) and 2015 (B). WQS models were adjusted for age, sex, marital status, residence, educational level, smoking status, drinking status, history of hypertension, diabetes, heart disease, dyslipidaemia, kidney disease, medication use for hypertension, medication use for diabetes, medication use for dyslipidaemia, systolic blood pressure, diastolic blood pressure, total cholesterol, HDL-C, LDL-C, HbA1c, and the eGFR. BMI: body mass index; FBG: fasting blood glucose; TyG: triglyceride-glucose; WQS: weighted quantile sum

Sensitivity analyses

Sensitivity analyses consistently yielded congruent results when performing complete data analyses (Additional file 1: Table S3), excluding participants who had heart disease (Additional file 1: Table S4), employing competing risk regression to account for competing risks due to mortality (Additional file 1: Table S5), or employing Cox regression to account for censored data (Additional file 1: Table S6). Furthermore, the nonlinear nature of the association between cumulative TyG-BMI and the risk of incident stroke persisted when subjecting the data to complete data analyses (Additional file 1: Fig. S2) or treating cumulative TyG-BMI as a continuous variable (Additional file 1: Table S7).

Discussion

Our study revealed that substantial changes in the TyG-BMI are independently associated with the risk of stroke in individuals aged 45 years and above from the CHARLS national data. Notably, the association between the cumulative TyG-BMI and stroke occurrence exhibited a nonlinear pattern; specifically, the risk of stroke increased with a cumulative TyG-BMI < 570. Considering that the TyG-BMI incorporates FBG, TG, and BMI, our investigation identified TG as the primary contributor to the observed association. Meanwhile, our findings have important clinical implications, suggesting that monitoring long-term changes in the TyG-BMI may assist with the early identification of individuals at high risk of stroke. Furthermore, to mitigate the risk of stroke, prioritizing the management of TG levels may be worth considering.

A previous investigation provided evidence demonstrating the rapid induction of IR in humans through elevated levels of free fatty acids in plasma [31]. In addition, when exposed to high-glucose conditions, glucose molecules form conjugates with proteins resembling the insulin receptor on the cytoplasmic membrane, resulting in the prompt formation of advanced glycation end products [32]. Consequently, the binding of serum insulin to the deposited advanced glycation end products on the insulin receptor becomes imperfect, thereby impeding the mediation of insulin's glucose transport stimulation and triggering IR [33]. Recent studies have suggested that the product of plasma TyG holds promise as an effective measure for detecting IR [6, 34]. Moreover, BMI, a straightforwards anthropometric parameter commonly employed as an indicator of obesity and IR, is noteworthy. In individuals with obesity, adipose tissue lipolysis intensifies, leading to the release of substantial amounts of free fatty acids, which represents a crucial factor influencing insulin sensitivity modulation [35]. Consequently, it is reasonable to hypothesize that the TyG-BMI, derived from the anthropometric BMI and TyG parameters, may also serve as a valuable marker for IR.

Indeed, the association between the TyG-BMI and HOMA-IR has been empirically established [13]. Consequently, the TyG-BMI has been advocated as a reliable metric for assessing IR and IR-related ailments in numerous studies. A cross-sectional investigation carried out in rural Beijing, China affirmed that the TyG-BMI exhibited superior efficacy in detecting IR [36]. Moreover, the 2015 Health, Well-Being, and Ageing Study proposed the utility of TyG-BMI as an assessment tool for prediabetes, albeit not as the optimal index [37]. Findings from a substantial cross-sectional survey involving 11,149 participants, the Korean National Health and Nutrition Examination Survey, demonstrated that TyG-BMI served as a viable alternative marker for evaluating IR when compared to other IR parameters [38]. Additionally, the associations between TyG-BMI and prehypertension, as well as hypertension, have been conclusively validated [39, 40]. Despite the identification of an association between TyG-BMI and ischaemic stroke in a recent cross-sectional study [19], certain limitations, such as the study's cross-sectional design and focus on stroke subtype, warrant caution in drawing definitive conclusions. Consequently, the precise association between TyG-BMI and ischaemic stroke remains equivocal.

Based on our current understanding, this study represents a novel approach in utilizing cluster analysis to categorize the changes in the TyG-BMI values. Each category within the analysis corresponded to distinct subpopulations, wherein individuals with a consistently low TyG-BMI exhibited the lowest risk, while those with a highest TyG-BMI and a slow declining trend displayed the highest risk. Previous studies predominantly relied on a single TyG-BMI value to predict the occurrence of stroke, often yielding different results on different occasions [19]. Furthermore, it is worth noting that our research sample consisted of a representative cross-section of healthy individuals from various regions of China. By focusing on dynamic processes, our investigation has contributed further evidence to elucidate the association between TyG-BMI and stroke. Specifically, our RCS model elucidated a nonlinear association between the cumulative TyG-BMI and stroke. However, the underlying mechanistic explanation for this observed association remains uncertain. Finally, we employed the WQS regression method to augment the interpretability of the TyG-BMI, wherein we observed that TG emerged as the primary contributor to the observed effects.

The precise mechanism by which IR contributes to the development of stroke remains poorly elucidated, with several potential pathways warranting consideration. First, IR is believed to induce endothelial dysfunction, foam cell formation, and the formation of vulnerable plaques, thereby playing a critical role in the pathogenesis of atherosclerosis [41,42,43,44]. Additionally, IR, characterized by a low-grade inflammatory state, facilitates the progression of atherosclerosis and stimulates the production of inflammatory markers [45, 46]. Second, IR exerts an impact on platelet adhesion, activation, and aggregation [47,48,49,50], culminating in stroke occurrence through arterial stenosis or occlusion. Third, IR has been associated with heightened sympathetic nervous system activity [34] and impaired cardiac autonomic function [51], both of which contribute to the pathophysiology of acute cardiovascular and cerebrovascular diseases. Last, individuals with IR often exhibit a larger waist circumference and BMI, hypertension, diabetes, cardiovascular disease, and a history of dyslipidaemia. Moreover, they frequently present with elevated fasting blood glucose levels, triglycerides, and glycosylated haemoglobin, all of which represent established risk factors for stroke [51,52,53,54].

Our study contributes to the existing body of knowledge by providing evidence supporting the utilization of the dynamic change in the TyG-BMI as a clinically valuable marker for identifying individuals at a heightened risk of cardiovascular disease. The TyG-BMI was derived through the calculation of FBG, TG, and BMI. However, apart from measuring TG and FBG, body weight and height are also needed to calculate the TyG-BMI, which results in a more complicated formula. We still assert that the TyG-BMI holds promising prospects for accurately identifying patients at a heightened risk of stroke. First, these biochemical parameters can be conveniently obtained from a single sample at the same time, presenting a cost-effective and convenient alternative to the euglycaemic-hyperinsulinaemic clamp method. Furthermore, the widespread availability and routine performance of height, weight, FBG, and TG measurements in primary healthcare settings enhance the generalizability and practicality of the TyG-BMI in both clinical and epidemiological studies. Second, the TyG-BMI is designed to offer a more comprehensive assessment of insulin resistance, particularly in populations where obesity is prevalent. By incorporating BMI, which reflects body adiposity, the TyG-BMI provides a more holistic evaluation that accounts for the influence of obesity on IR [13]. Third, obesity is a well-established risk factor for IR and metabolic disorders [14, 55]. In clinical practice, addressing weight-related factors is often a key aspect of managing IR. The TyG-BMI aligns with this clinical relevance, potentially guiding interventions and treatment strategies for individuals with obesity-related IR. Last, the ROC analysis demonstrated that baseline TyG-BMI had greater accuracy in predicting stroke compared with TyG.

However, it is essential to address several limitations associated with this study. First, while the TyG-BMI has demonstrated reliability and convenience as a surrogate marker for insulin resistance, establishing a direct association between IR and stroke necessitates a comparison with the gold standard diagnostic method. Unfortunately, this study did not undertake such a comparison, thereby limiting its ability to provide a direct explanation of the insulin resistance-stroke association. Second, the study only incorporated two blood tests, precluding a comprehensive assessment and refinement of the TyG-BMI. Third, it should be noted that, similar to other studies, the diagnosis of stroke in this research relied on self-reporting, which introduces a logistical constraint. Due to the absence of medical records in the CHARLS dataset, validation and confirmation of self-reported incident stroke were not feasible. However, it is worth mentioning that other large-scale studies, such as the English Longitudinal Study of Ageing, have reported a satisfactory level of agreement between self-reported incident stroke and medical records [56]. Fourth, this study adopts an observational design, and it is essential to acknowledge the presence of selection bias arising from the loss to follow-up. This bias is influenced by factors such as nonresponse due to severe cardiovascular disease or competing events resulting from mortality, which may potentially lead to an underestimation of the association between the TyG-BMI and stroke. Nevertheless, it is noteworthy that our findings remained consistent even after we conducted exclusions involving participants with pre-existing heart disease. Furthermore, we employed Cox regression or competing risk regression analyses, which also yielded results consistent with our primary findings. Last, only participants from China were involved in this study; thus, the findings may not be fully generalizable to other countries.

Conclusions

In this study, we revealed that substantial changes in the TyG-BMI are independently associated with the risk of stroke in individuals aged 45 years and above from the CHARLS national data. Consequently, monitoring long-term changes in the TyG-BMI should prioritize stroke prevention strategies. Furthermore, our findings elucidated the underlying mechanisms of the TyG-BMI by highlighting TG as the primary contributor to the observed effects.

Data availability

Online repositories contain the datasets used in this investigation. The names of the repositories and accession numbers can be found at http://charls.pku.edu.cn/en.

References

  1. O’Donnell MJ, Chin SL, Rangarajan S, et al. Global and regional effects of potentially modifiable risk factors associated with acute stroke in 32 countries (INTERSTROKE): a case-control study. Lancet. 2016;388(10046):761–75.

    PubMed  Google Scholar 

  2. Feigin VL, Krishnamurthi RV, Parmar P, et al. Update on the global burden of ischemic and hemorrhagic stroke in 1990–2013: the GBD 2013 study. Neuroepidemiology. 2015;45(3):161–76.

    PubMed  Google Scholar 

  3. Wang W, Jiang B, Sun H, et al. Prevalence, incidence, and mortality of stroke in china: results from a nationwide population-based survey of 480 687 adults. Circulation. 2017;135(8):759–71.

    PubMed  Google Scholar 

  4. Guan T, Ma J, Li M, et al. Rapid transitions in the epidemiology of stroke and its risk factors in China from 2002 to 2013. Neurology. 2017;89(1):53–61.

    PubMed  Google Scholar 

  5. Kernan WN, Inzucchi SE, Viscoli CM, et al. Impaired insulin sensitivity among nondiabetic patients with a recent TIA or ischemic stroke. Neurology. 2003;60(9):1447–51.

    PubMed  CAS  Google Scholar 

  6. Guerrero-Romero F, Simental-Mendia LE, Gonzalez-Ortiz M, et al. 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.

    PubMed  CAS  Google Scholar 

  7. Smith GI, Mittendorfer B, Klein S. Metabolically healthy obesity: facts and fantasies. J Clin Invest. 2019;129(10):3978–89.

    PubMed  PubMed Central  Google Scholar 

  8. Cai X, Zhang Y, Li M, et al. Association between prediabetes and risk of all cause mortality and cardiovascular disease: updated meta-analysis. BMJ. 2020;370: m2297.

    PubMed  PubMed Central  Google Scholar 

  9. Chang Y, Kim CK, Kim MK, Seo WK, Oh K. Insulin resistance is associated with poor functional outcome after acute ischemic stroke in non-diabetic patients. Sci Rep. 2021;11(1):1229.

    PubMed  PubMed Central  CAS  Google Scholar 

  10. Ago T, Matsuo R, Hata J, et al. Insulin resistance and clinical outcomes after acute ischemic stroke. Neurology. 2018;90(17):e1470–7.

    PubMed  CAS  Google Scholar 

  11. Pan Y, Jing J, Chen W, et al. Post-glucose load measures of insulin resistance and prognosis of nondiabetic patients with ischemic stroke. J Am Heart Assoc. 2017;6(1): e004990.

    PubMed  PubMed Central  Google Scholar 

  12. Wang S, Shi J, Peng Y, et al. Stronger association of triglyceride glucose index than the HOMA-IR with arterial stiffness in patients with type 2 diabetes: a real-world single-centre study. Cardiovasc Diabetol. 2021;20(1):82.

    PubMed  PubMed Central  CAS  Google Scholar 

  13. Er LK, Wu S, Chou HH, et al. Triglyceride glucose-body mass index is a simple and clinically useful surrogate marker for insulin resistance in nondiabetic individuals. PLoS ONE. 2016;11(3): e0149731.

    PubMed  PubMed Central  Google Scholar 

  14. Kahn BB, Flier JS. Obesity and insulin resistance. J Clin Invest. 2000;106(4):473–81.

    PubMed  PubMed Central  CAS  Google Scholar 

  15. Wang X, Liu J, Cheng Z, Zhong Y, Chen X, Song W. Triglyceride glucose-body mass index and the risk of diabetes: a general population-based cohort study. Lipids Health Dis. 2021;20(1):99.

    PubMed  PubMed Central  CAS  Google Scholar 

  16. Li Y, Gui J, Liu H, et al. Predicting metabolic syndrome by obesity- and lipid-related indices in mid-aged and elderly Chinese: a population-based cross-sectional study. Front Endocrinol (Lausanne). 2023;14:1201132.

    PubMed  Google Scholar 

  17. Wang R, Dai L, Zhong Y, Xie G. Usefulness of the triglyceride glucose-body mass index in evaluating nonalcoholic fatty liver disease: insights from a general population. Lipids Health Dis. 2021;20(1):77.

    PubMed  PubMed Central  CAS  Google Scholar 

  18. Raimi TH, Dele-Ojo BF, Dada SA, et al. Triglyceride-glucose index and related parameters predicted metabolic syndrome in nigerians. Metab Syndr Relat Disord. 2021;19(2):76–82.

    PubMed  PubMed Central  CAS  Google Scholar 

  19. Du Z, Xing L, Lin M, Sun Y. Estimate of prevalent ischemic stroke from triglyceride glucose-body mass index in the general population. BMC Cardiovasc Disord. 2020;20(1):483.

    PubMed  PubMed Central  CAS  Google Scholar 

  20. Zhao Y, Hu Y, Smith JP, Strauss J, Yang G. Cohort profile: the China Health and Retirement Longitudinal Study (CHARLS). Int J Epidemiol. 2014;43(1):61–8.

    PubMed  Google Scholar 

  21. von Elm E, Altman DG, Egger M, et al. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Epidemiology. 2007;18(6):800–4.

    Google Scholar 

  22. Cui H, Liu Q, Wu Y, Cao L. Cumulative triglyceride-glucose index is a risk for CVD: a prospective cohort study. Cardiovasc Diabetol. 2022;21(1):22.

    PubMed  PubMed Central  CAS  Google Scholar 

  23. Zheng F, Yan L, Zhong B, Yang Z, Xie W. Progression of cognitive decline before and after incident stroke. Neurology. 2019;93(1):e20–8.

    PubMed  Google Scholar 

  24. Li H, Zheng D, Li Z, et al. Association of depressive symptoms with incident cardiovascular diseases in middle-aged and older Chinese adults. JAMA Netw Open. 2019;2(12): e1916591.

    PubMed  PubMed Central  Google Scholar 

  25. Chen X, Crimmins E, Hu PP, et al. Venous blood-based biomarkers in the china health and retirement longitudinal study: rationale, design, and results from the 2015 wave. Am J Epidemiol. 2019;188(11):1871–7.

    PubMed  PubMed Central  Google Scholar 

  26. Singh A, Yadav A, Rana A. K-means with three different distance metrics. Int J Comput Appl. 2013;67(10):13.

    Google Scholar 

  27. Sinaga KP, Yang M-S. Unsupervised K-means clustering algorithm. IEEE Access. 2020;8:80716–27.

    Google Scholar 

  28. Vieth E. Fitting piecewise linear regression functions to biological responses. J Appl Physiol. 1989;67(1):390–6.

    PubMed  CAS  Google Scholar 

  29. Nakamura T. BMDP program for piecewise linear regression. Comput Methods Programs Biomed. 1986;23(1):53–5.

    PubMed  CAS  Google Scholar 

  30. Tanner EM, Bornehag CG, Gennings C. Repeated holdout validation for weighted quantile sum regression. MethodsX. 2019;6:2855–60.

    PubMed  PubMed Central  Google Scholar 

  31. Roden M, Price TB, Perseghin G, et al. Mechanism of free fatty acid-induced insulin resistance in humans. J Clin Invest. 1996;97(12):2859–65.

    PubMed  PubMed Central  CAS  Google Scholar 

  32. Ruiz HH, Ramasamy R, Schmidt AM. Advanced glycation end products: building on the concept of the “Common Soil” in metabolic disease. Endocrinology. 2020;161(1): bqz006.

    PubMed  Google Scholar 

  33. Shulman GI. Ectopic fat in insulin resistance, dyslipidemia, and cardiometabolic disease. N Engl J Med. 2014;371(12):1131–41.

    PubMed  Google Scholar 

  34. Simental-Mendia LE, Rodriguez-Moran 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  CAS  Google Scholar 

  35. Morigny P, Houssier M, Mouisel E, Langin D. Adipocyte lipolysis and insulin resistance. Biochimie. 2016;125:259–66.

    PubMed  CAS  Google Scholar 

  36. Tuo X, Yuan J, Wang XH, Xin Z. Identifying the insulin resistance index in nondiabetic Chinese subjects. Medicine (Baltimore). 2020;99(5): e19023.

    PubMed  Google Scholar 

  37. Ramirez-Velez R, Perez-Sousa MA, Gonzalez-Ruiz K, et al. Obesity- and lipid-related parameters in the identification of older adults with a high risk of prediabetes according to the American Diabetes Association: an analysis of the 2015 Health, Well-Being, and Aging Study. Nutrients. 2019;11(11):2654.

    PubMed  PubMed Central  CAS  Google Scholar 

  38. Lim J, Kim J, Koo SH, Kwon GC. Comparison of triglyceride glucose index, and related parameters to predict insulin resistance in Korean adults: an analysis of the 2007–2010 Korean National Health and Nutrition Examination Survey. PLoS ONE. 2019;14(3): e0212963.

    PubMed  PubMed Central  CAS  Google Scholar 

  39. Zeng ZY, Liu SX, Xu H, Xu X, Liu XZ, Zhao XX. Association of triglyceride glucose index and its combination of obesity indices with prehypertension in lean individuals: a cross-sectional study of Chinese adults. J Clin Hypertens (Greenwich). 2020;22(6):1025–32.

    PubMed  CAS  Google Scholar 

  40. Bala C, Gheorghe-Fronea O, Pop D, et al. The association between six surrogate insulin resistance indexes and hypertension: a population-based study. Metab Syndr Relat Disord. 2019;17(6):328–33.

    PubMed  CAS  Google Scholar 

  41. Guo W, Zhu W, Wu J, et al. Triglyceride glucose index is associated with arterial stiffness and 10-year cardiovascular disease risk in a Chinese population. Front Cardiovasc Med. 2021;8: 585776.

    PubMed  PubMed Central  CAS  Google Scholar 

  42. Miao M, Zhou G, Bao A, et al. Triglyceride-glucose index and common carotid artery intima-media thickness in patients with ischemic stroke. Cardiovasc Diabetol. 2022;21(1):43.

    PubMed  PubMed Central  CAS  Google Scholar 

  43. Wu S, Xu L, Wu M, Chen S, Wang Y, Tian Y. Association between triglyceride-glucose index and risk of arterial stiffness: a cohort study. Cardiovasc Diabetol. 2021;20(1):146.

    PubMed  PubMed Central  CAS  Google Scholar 

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

    PubMed  PubMed Central  CAS  Google Scholar 

  45. Rohm TV, Meier DT, Olefsky JM, Donath MY. Inflammation in obesity, diabetes, and related disorders. Immunity. 2022;55(1):31–55.

    PubMed  PubMed Central  CAS  Google Scholar 

  46. Jin A, Wang S, Li J, et al. Mediation of systemic inflammation on insulin resistance and prognosis of nondiabetic patients with ischemic stroke. Stroke. 2023;54(3):759–69.

    PubMed  CAS  Google Scholar 

  47. Du T, Yuan G, Zhang M, Zhou X, Sun X, Yu X. Clinical usefulness of lipid ratios, visceral adiposity indicators, and the triglycerides and glucose index as risk markers of insulin resistance. Cardiovasc Diabetol. 2014;13:146.

    PubMed  PubMed Central  Google Scholar 

  48. Ferreiro JL, Gomez-Hospital JA, Angiolillo DJ. Platelet abnormalities in diabetes mellitus. Diab Vasc Dis Res. 2010;7(4):251–9.

    PubMed  Google Scholar 

  49. Ueno M, Fujita K, Yamamoto H, et al. Impact of impaired glucose tolerance on clopidogrel response in patients with coronary artery disease. J Thromb Thrombolysis. 2015;40(2):174–81.

    PubMed  CAS  Google Scholar 

  50. Santilli F, Vazzana N, Liani R, Guagnano MT, Davi G. Platelet activation in obesity and metabolic syndrome. Obes Rev. 2012;13(1):27–42.

    PubMed  CAS  Google Scholar 

  51. Shi W, Xing L, Jing L, et al. Value of triglyceride-glucose index for the estimation of ischemic stroke risk: insights from a general population. Nutr Metab Cardiovasc Dis. 2020;30(2):245–53.

    PubMed  CAS  Google Scholar 

  52. Zhou Y, Pan Y, Yan H, et al. Triglyceride glucose index and prognosis of patients with ischemic stroke. Front Neurol. 2020;11:456.

    PubMed  PubMed Central  Google Scholar 

  53. Wang A, Tian X, Zuo Y, et al. Change in triglyceride-glucose index predicts the risk of cardiovascular disease in the general population: a prospective cohort study. Cardiovasc Diabetol. 2021;20(1):113.

    PubMed  PubMed Central  CAS  Google Scholar 

  54. Zhao Y, Sun H, Zhang W, et al. Elevated triglyceride-glucose index predicts risk of incident ischaemic stroke: the Rural Chinese cohort study. Diabetes Metab. 2021;47(4): 101246.

    PubMed  CAS  Google Scholar 

  55. Karra P, Winn M, Pauleck S, et al. Metabolic dysfunction and obesity-related cancer: beyond obesity and metabolic syndrome. Obesity (Silver Spring). 2022;30(7):1323–34.

    PubMed  CAS  Google Scholar 

  56. Steptoe A, Breeze E, Banks J, Nazroo J. Cohort profile: the English longitudinal study of ageing. Int J Epidemiol. 2013;42(6):1640–8.

    PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank all the members of the CHALRS for their contributions and the participants who contributed their data.

Funding

None.

Author information

Authors and Affiliations

Authors

Contributions

XY conceived the study; all authors participated in the acquisition of the data; RH, LZ and QL analysed the data; RH and LZ drafted and revised the manuscript; and all the authors have read and approved the final version of the manuscript.

Corresponding author

Correspondence to Xue-Mei You.

Ethics declarations

Ethics approval and consent to participate

The studies involving human participants were reviewed and approved by the Ethics Review Committee of Peking University. The patients/participants provided their written informed consent to participate in this study.

Competing interests

The authors declare no potential conflicts of interest.

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: Figure S1

. Receiver operating characteristic curves (ROCs) for baseline TyG-BMI predicting stroke. Figure S2. Nonlinear association between cumulative TyG-BMI and stroke in subpopulations of 4373 participants with complete data. Table S1. Baseline characteristics between participants included and not included. Table S2. Baseline characteristics of 4583 participants according to the quartile of cumulative TyG-BMI. Table S3. Associations of different classes of TyG-BMI with stroke incidence in subpopulations of 4373 participants with complete data. Table S4. Associations of different classes of TyG-BMI with stroke incidence in subpopulations of 4045 participants without heart disease. Table S5. Associations of different classes of TyG-BMI with stroke incidence using competing risk regression. Table S6. Associations of different classes of TyG-BMI with stroke incidence using the Cox proportional hazards regression. Table S7. Associations of cumulative TyG-BMI with stroke incidence when treating cumulative TyG-BMI as a continuous variable.

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Huo, RR., Zhai, L., Liao, Q. et al. Changes in the triglyceride glucose-body mass index estimate the risk of stroke in middle-aged and older Chinese adults: a nationwide prospective cohort study. Cardiovasc Diabetol 22, 254 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s12933-023-01983-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12933-023-01983-5

Keywords