Research Article - Diabetes Management (2022) Volume 12, Issue 5

Association between triglyceride glucose-body mass index and incident diabetes mellitus in Chinese population: a secondary retrospective analysis

Corresponding Author:
Haofei Hu
Department of Nephrology, Shenzhen Second People’s Hospital, Shenzhen 518000, Guangdong Province, China
E-mail: huhaofei0319@126.com
Corresponding Author:
Xinyu Wang
Shenzhen University Health Science Center, Shenzhen 518000, Guangdong Province, China
E-mail: wxyhorse@126.com

Received: 09-Oct-2022, Manuscript No. FMDM-22-76782; Editor assigned: 11-Oct-2022, PreQC No. FMDM-22-76782 (PQ); Reviewed: 25-Oct-2022, QC No. FMDM-22-76782; Revised: 31-Oct-2022, Manuscript No. FMDM-22-76782 (R); Published: 07-Nov-2022, DOI: 10.37532/1758-1907.2022.12(5).405-416

Abstract

Background: TyG-BMI is one of the insulin resistance markers. However, few studies have investigated the relationship between TyG-BMI and diabetes. The aim of this study is to investigate the relationship between TyG-BMI and diabetes occurrence in Chinese individuals.

Methods: This retrospective study is a population-based cohort of 204978 non-diabetic individuals using healthy screening program data in China between 2010 and 2016. Cox proportional hazards regression analysis evaluated the Hazard Ratio (HR) and 95% Confidence Interval (95% CI) of the relationship between TyG-BMI and diabetes events. Then the association between TyG-BMI and incident diabetes was further explored by generalized additive models and subgroup analyses.

Results: Our study indicated that a increase in diabetes incidence occurred with the rise of TyG-BMI (HR=1.023, 95% CI (1.022, 1.024) after the adjustment of potential confounders. TyG- BMI has a non-linear relationship with the incidence of diabetes, and the inflection point is 232.416. Subgroup analysis showed that the correlation was stronger in the population aged 20-30 (P for interaction<0.0001, HR 1.029, 95%CI:1.024 to 1.035), and the same trend was found in the following populations: age 30-40(HR=1.032), age 40-50(HR=1.029), HDL (high group) (HR=1.024 ), SBP<140(HR=1.025), DBP<90(HR=1.024), current drinker(HR=1.031), and ever drinker(HR=1.032).

Conclusion: TyG-BMI and incident diabetes had a non-linear positive relationship. Before and after TyG-BMI equals 232.416, the risk of diabetes increased by 2.9% and 1.6%, respectively, when TyG-BMI increased one unit. The findings emphasise the importance of baseline TyG-BMI as a measure for identification of patients at risk of diabetes early and improve their outcomes.

Keywords

■ TyG-BMI

■ incident diabetes

■ nonlinearity

Abbreviations

BMI: Body Mass Index; TyG: Triglycerideglucose Index; TyG-BMI: Triglyceride Glucose-body Mass Index; IR: Insulin Resistance; HRs: Hazard Ratios; 95%CIs: 95% Confidence Intervals; IDF: International Diabetes Federation; DM: Diabetes Mellitus; T2DM: Type 2 Diabetes; TG: Triglyceride; FPG: Fasting Plasma Glucose; LDL-C: Low-Density Lipoprotein Cholesterol; HDL-C: High-density Lipoprotein Cholesterol; TC: Total Cholesterol; BUN: Serum Urea Nitrogen; Scr: Serum Creatinine; AST: Aspartate Aminotransferase; ALT: Alanine Aminotransferase; SBP: Systolic Blood Pressure; DBP: Diastolic Blood Pressure; GAM: Generalized Additive Models; SD: Standard Deviation; IQR: Interquartile Range.

Introduction

Diabetes has become a critical worldwide healthy problem with high prevalence. China has the largest number of patients with diabetes in the world [1]. Diabetes Mellitus (DM) and its complications have resulted in a severe economic burden of mortality and disability. Insulin Resistance (IR) and the consequences of compensatory hyperinsulinemia are vital pathological mechanisms of diabetes mellitus and obesity [2]. Therefore, it is important to recognize IR before clinical manifestations of diabetes.

Recently, as Triglyceride Glucose Body Mass Index (TyG-BMI) combines Triglyceride (TG), Fasting Plasma Glucose (FPG), and obesity status, TyG-BMI has been considered to identify IR more reliably than TyG [3,4]. Triglycerideglucose (TyG) index, which is estimated using the formula Ln (fasting triglycerides (mg/dl) × fasting blood glucose (mg/dl)/2), is an alternative to identify insulin resistance in apparently healthy subjects [5-9]. TyG-BMI=TyG × BMI, which considered BMI on the basis of TyG [3]. Body mass index (BMI) is an easily detectable, inexpensive, and non- invasive measurement parameter closely related to IR. Patients with higher levels of TG, FBG, BMI may progress to diabetes over time [10-12]. Several cross-sectional studies have established the relationship between TyG-BMI and DM. At present, there are few cohort study articles about the relationship between TyG-BMI and diabetes.

In this study, we postulated that the TyG-BMI level might serve as an early predictor of incident DM in China. This study was performed based on a previously published data which identified the association of body mass index and age with incident diabetes [12]. To test our hypothesis, Cox proportional hazards regression analysis was used to evaluate the Hazard Ratio (HR) and 95% Confidence Interval (95% CI) of the relationship between TyG-BMI and diabetes events. Non-linear relationships were assessed using generalized additive models. Subgroup analysis helped comprehend other factors that may affect the relationship between TyG-BMI and diabetes in order to discover potential special populations. Our findings demonstrated that tracking the TyG-BMI level aids in the prediction of a diabetes risk. TyG-BMI should be taken into account for clinicians to plan and initiate the management strategies early to improve outcomes for participants with prediabetes.

Methodology

Methods including study population and design, data collection and measurements, statistical analysis are the same as the preprint of this manuscript [13].

■ Study population and design

In this population-based cohort analysis of a medical program established by the Rich Healthcare Group in China between 2010 and 2016, we investigated the effect of TyG- BMI on incident diabetes. We used data from the DATADRYAD website (www.datadryad. org), allowing others to use the database for free. In order to be consistent with the terms of service, this study quoted the data package shared by Chen. The database materials included the following variables: BMI, age, sex, smoking, drinking, family history of diabetes, High-Density Lipoprotein Cholesterol (HDL-C), Total Cholesterol (TC), Low-Density Lipoprotein Cholesterol (LDL-C), TG, Serum urea nitrogen(BUN), FPG, Serum Creatinine (Scr), Alanine Aminotransferase (ALT), Aspartate Aminotransferase (AST), Diastolic Blood Pressure (DBP), Systolic Blood Pressure (SBP), incident diabetes at follow up and follow- up time. The authors claimed that they had given up the copyright of the original paper and related ownership of the database. As for ethics approval, the study was a retrospective analysis approved by the rich healthcare group review committee.

The initial study recruited 685,277 participants ≥ 20 years of age. They conducted at least two visits in 32 locations and 11 cities in China (Shanghai, Nanjing, Beijing, Suzhou, Changzhou, Shenzhen, Guangzhou, Chengdu, Hefei, Nantong, Wuhan) between 2010 and 2016. The data we got has been filtered. Subjects were excluded with the following conditions: (1) missing information about height, weight, fasting plasma glucose, and gender at baseline; (2) extreme BMI values (<15 kg/m2 or >55 kg/ m2); (3) persons were diagnosed with diabetes at baseline and were diagnosed with unclear diabetes status at follow-up; (4) visit intervals less than two years. Finally, in the original study, Ying kept 211,833 participants [12]. Missing values of baseline TG (n=4,887) and zero values of baseline TG(n=860) were excluded from the analysis cohort for further research. And then, TyG-BMI calculation formula was BMI × Ln(fasting triglycerides(mg/dl) × fasting blood glucose(mg/dl)/2). We excluded outliers of TyG- BMI (out of means minus or plus three standard deviation (SD) (n=1,108). Finally, a sum of 204978 participants was selected in our study. FIGURE 1 depicted the participants’ selection process.

Figure 1: Flowchart of study population.

■ Data collection and measurements

The original database contained participants’ clinical history and lifestyle factors based on a standardized questionnaire about lifestyle factors, demographic characteristics, family history of chronic disease, and personal medical history. Trained staff measured height, weight, and blood pressure. The body weight was measured in light clothing without shoes and was accurate to 0.1 kg. The measurement height was accurate to 0.1 cm. The body mass index was calculated by dividing the weight in kilograms by the height in meters. Standard mercury sphygmomanometers measured blood pressure. After fasting for at least 10 hours at each visit, fasting venous blood samples were collected. Triglyceride, LDL-C, HDL-C, total cholesterol, and plasma glucose levels were measured on an automatic analyzer (Beckman 5800). The glucose oxidase method was used to determine the plasma glucose level. The TyG index formula was Ln (fasting triglycerides (mg/dl) × fasting blood glucose (mg/dl)/2)). The target-independent variable is TyG-BMI, which equals the BMI × TyG index. The dependent variable is incident diabetes, defined as fasting plasma glucose ≥ 7 mmol/L, and/or self-reported during the follow-up. This is a retrospective cohort study, so it cuts down the possibility of selection bias and observation bias.

■ Statistical analysis

Normally distributed continuous variables were represented by the mean, Standard Deviation (SD), and non-normally distributed continuous variables were replaced by the median of the Interquartile Range (IQR). Categorical variables were expressed as frequencies and percentages. For handling missing values, missing continuous variables were replaced with a mean or median depending on the distribution. Missing categorical variables could be a new categorical group. Stratified by TyG-BMI index quartiles, statistical differences of the groups were described with Kruskal Wallis H (skewed distribution) test, one-way ANOVA (normal distribution), and chi-square test (categorical variables). In order to explore the relationship between TyG-BMI and the incidence of diabetes, univariate and multivariate Cox proportional hazards regression analysis was used to estimate Hazard Ratios (HRs) and 95% Confidence Intervals (CIs). According to the recommendation of the STROBE statement, We used three models: crude model; model I adjust for: Gender, Age, Smoking Status, DBP, SBP, Family History, Drinking Status; model II adjust for: Gender, Age, Smoking Status, DBP, SBP, Family History, Drinking Status, TC, HDL, LDL, AST, Scr, ALT. To ensure the robustness of the data analysis, we conducted a sensitivity analysis. The TyG-BMI was converted into a categorical variable and calculated the P for trend to perform the linear trend tests. The Cox proportional hazards regression model with cubic spline function was used to identify non-linear relationships. In addition, if there was an apparent smooth curve, the recursive method would automatically calculate the inflection point. The associations of TyG-BMI with incident diabetes in subgroups were also studied using a stratified linear regression model and likelihood ratio test to discover corrections and interactions. The subgroups were classified by age (20 to<30, 30 to<40, 40 to<50, 50 to<60, 60 to<70, ≥ 70), gender (male vs. female), HDL( low, middle, high, not recorded), LDL( low, middle, high, not recorded), SBP(<140, ≥ 140 ), DBP(<90, ≥ 90), Smoking status(current smoker, ever smoker, never smoker, not recorded), Drinking status(current drinker, ever drinker, never drinker, not recorded), Family history of diabetes(no, yes). By using the time to the first event for each endpoint, the Kaplan- Meier method was used to compare survival estimates and cumulative event rates.

All analyses were performed with the statistical software package R (http://www.R-project.org, The R Foundation) and Empower-Stats (http:// www.empowerstats.com, X and Y Solutions, Inc., Boston, MA). Significance was accepted at a two-tailed P<0.05.

Results

■ Baseline characteristics of the study participants

Of the 211,833 subjects recruited in the former study, 204,978 participants were included in the current analysis. The mean age of the population was 42.2 ± 12.7 years old, and 45.17% of participants were women. The average follow- up time was 3.1 ± 0.9 years, and 4093 subjects developed diabetes during the follow-up period.

The mean TyG-BMI was 158.7 ± 32.4, and the mean FPG, BMI, TG were 88.6 ± 10.9 mg/dl, 23.3 ± 3.3 kg/m2 and 24.2 ± 18.6 mg/ dl respectively. Individuals in the highest TyG- BMI group (Q4) were generally older than those in the lowest TyG-BMI group (Q1) and had higher BMI, FPG, TG, SBP, DBP, TC, ALT, SCR values. What’s more, as the value of TyG-BMI increased, the incidence of diabetes gradually increased. (Q1: 0.23% vs. Q2: 0.62% vs. Q3: 1.74% vs. Q4: 5.35%). Compared with the Q1 group, the Q4 group had lower HDL levels, higher AST and LDL levels, higher rates of drinking, smoking, family history (TABLE 1).

TyG-BMI Q1(≤ 167.01) Q2(167.01 to ≤ 191.53) Q3(191.53 to ≤ 219.26) Q4(>219.26) P-value
Participants 51245 51244 51244 51245  
Age(years) 36.65 ± 10.06 41.26 ± 12.13 44.79 ± 13.08 46.16 ± 13.11 <0.001
BMI(kg/m2 19.45 ± 1.36 21.98 ± 1.20 24.11 ± 1.32 27.21 ± 2.14 <0.001
FPG(mmol/L) 4.68 ± 0.53 4.84 ± 0.55 4.98 ± 0.58 5.17 ± 0.63 <0.001
TG(mmol/L) 0.67 (0.52-0.86) 0.91 (0.70-1.18) 1.25 (0.96-1.65)  1.90 (1.40-2.62) <0.001
TyG-BMI 151.98 ± 10.59 179.21 ± 7.07 204.83 ± 7.94 244.14 ± 19.76 <0.001
SBP(mmHg) 110.87 ± 13.45 116.16 ± 14.96 121.63 ± 15.68 127.36 ± 16.40 <0.001
DBP(mmHg) 69.34 ± 9.12 72.00 ± 9.74 75.54 ± 10.30 79.63 ± 11.02 <0.001
TC(mmol/L) 4.36 ± 0.79 4.60 ± 0.85 4.82 ± 0.88 5.05 ± 0.92 <0.001
ALT(U/L) 13.00 (10.30,17.20) 15.90 (12.00-21.90) 20.00 (15.00-28.60) 28.00 (19.30-41.70) <0.001
Scr(umol/L) 63.93 ± 13.20 68.29 ± 15.49 72.58 ± 14.84 75.36 ± 15.64 <0.001
Gender         <0.001
 Male 14983 (29.24%) 24120 (47.07%) 33452 (65.28%) 39551 (77.18%)  
 Female 36262 (70.76%) 27124 (52.93%) 17792 (34.72%) 11694 (22.82%)  
Family history         <0.001
 NO 50327 (98.21%) 50095 (97.76%) 50187 (97.94%) 50179 (97.92%)  
 YES 918 (1.79%) 1149 (2.24%) 1057 (2.06%) 1066 (2.08%)  
HDL(mmol/L)         <0.001
 Low 4321 (8.43%) 7255 (14.16%) 11440 (22.32%) 15213 (29.69%)  
 Middle 8289 (16.18%) 10160 (19.83%) 10449 (20.39%) 9178 (17.91%)  
 High 14117 (27.55%) 11538 (22.52%) 8407 (16.41%) 5957 (11.62%)  
 Not record 24518 (47.84%) 22291 (43.50%) 20948 (40.88%) 20897 (40.78%)  
LDL(mmol/L)         <0.001
 Low 13027 (25.42%) 10258 (20.02%) 8164 (15.93%) 7193 (14.04%)  
 Middle 8693 (16.96%) 10309 (20.12%) 10299 (20.10%) 9935 (19.39%)  
 High 5080 (9.91%) 8594 (16.77%) 12194 (23.80%) 13712 (26.76%)  
 Not record 24445 (47.70%) 22083 (43.09%) 20587 (40.17%) 20405 (39.82%)  
AST(U/L)         <0.001
 Low 10327 (20.15%) 8482 (16.55%) 6057 (11.82%) 3722 (7.26%)  
 Middle 7088 (13.83%) 7635 (14.90%) 7815 (15.25%) 6243 (12.18%)  
 High 3699 (7.22%) 5490 (10.71%) 7956 (15.53%) 11777 (22.98%)  
 Not record 30131 (58.80%) 29637 (57.84%) 29416 (57.40%) 29503 (57.57%)  
Smoking status         <0.001
 Current smoker 1282 (2.50%) 2161 (4.22%) 3248 (6.34%) 4841 (9.45%)  
 Ever smoker 298 (0.58%) 507 (0.99%) 768 (1.50%) 893 (1.74%)  
 Never smoker 11551 (22.54%) 11271 (21.99%) 11017 (21.50%) 10307 (20.11%)  
 Not record 38114 (74.38%) 37305 (72.80%) 36211 (70.66%) 35204 (68.70%)  
Drinking status         <0.001
 Current drinker 108 (0.21%) 209 (0.41%) 386 (0.75%) 615 (1.20%)  
 Ever drinker 1059 (2.07%) 1850 (3.61%) 2572 (5.02%) 3234 (6.31%)  
 Never drinker 11964 (23.35%) 11880 (23.18%) 12075 (23.56%) 12192 (23.79%)  
 Not record 38114 (74.38%) 37305 (72.80%) 36211 (70.66%) 35204 (68.70%)  

TABLE 1: The baseline characteristics of participants.

■ Univariate analysis

The results of the univariate analysis were shown in TABLE 2. TABLE 2 showed that men were more likely to develop diabetes than women, and BMI, SBP, TG, FPG, TyG-BMI, DBP, TC, AST, ALT, LDL, Scr, age, drinking, smoking, and family history were all positively correlated with the incidence of diabetes.

TyG-BMI Q1(≤ 167.01) Q2(167.01 to ≤ 191.53) Q3(191.53 to ≤ 219.26) Q4(>219.26) P-value
Participants 51245 51244 51244 51245  
Age(years) 36.65 ± 10.06 41.26 ± 12.13 44.79 ± 13.08 46.16 ± 13.11 <0.001
BMI(kg/m2 19.45 ± 1.36 21.98 ± 1.20 24.11 ± 1.32 27.21 ± 2.14 <0.001
FPG(mmol/L) 4.68 ± 0.53 4.84 ± 0.55 4.98 ± 0.58 5.17 ± 0.63 <0.001
TG(mmol/L) 0.67 (0.52-0.86) 0.91 (0.70-1.18) 1.25 (0.96-1.65)  1.90 (1.40-2.62) <0.001
TyG-BMI 151.98 ± 10.59 179.21 ± 7.07 204.83 ± 7.94 244.14 ± 19.76 <0.001
SBP(mmHg) 110.87 ± 13.45 116.16 ± 14.96 121.63 ± 15.68 127.36 ± 16.40 <0.001
DBP(mmHg) 69.34 ± 9.12 72.00 ± 9.74 75.54 ± 10.30 79.63 ± 11.02 <0.001
TC(mmol/L) 4.36 ± 0.79 4.60 ± 0.85 4.82 ± 0.88 5.05 ± 0.92 <0.001
ALT(U/L) 13.00 (10.30,17.20) 15.90 (12.00-21.90) 20.00 (15.00-28.60) 28.00 (19.30-41.70) <0.001
Scr(umol/L) 63.93 ± 13.20 68.29 ± 15.49 72.58 ± 14.84 75.36 ± 15.64 <0.001
Gender         <0.001
 Male 14983 (29.24%) 24120 (47.07%) 33452 (65.28%) 39551 (77.18%)  
 Female 36262 (70.76%) 27124 (52.93%) 17792 (34.72%) 11694 (22.82%)  
Family history         <0.001
 NO 50327 (98.21%) 50095 (97.76%) 50187 (97.94%) 50179 (97.92%)  
 YES 918 (1.79%) 1149 (2.24%) 1057 (2.06%) 1066 (2.08%)  
HDL(mmol/L)         <0.001
 Low 4321 (8.43%) 7255 (14.16%) 11440 (22.32%) 15213 (29.69%)  
 Middle 8289 (16.18%) 10160 (19.83%) 10449 (20.39%) 9178 (17.91%)  
 High 14117 (27.55%) 11538 (22.52%) 8407 (16.41%) 5957 (11.62%)  
 Not record 24518 (47.84%) 22291 (43.50%) 20948 (40.88%) 20897 (40.78%)  
LDL(mmol/L)         <0.001
 Low 13027 (25.42%) 10258 (20.02%) 8164 (15.93%) 7193 (14.04%)  
 Middle 8693 (16.96%) 10309 (20.12%) 10299 (20.10%) 9935 (19.39%)  
 High 5080 (9.91%) 8594 (16.77%) 12194 (23.80%) 13712 (26.76%)  
 Not record 24445 (47.70%) 22083 (43.09%) 20587 (40.17%) 20405 (39.82%)  
AST(U/L)         <0.001
 Low 10327 (20.15%) 8482 (16.55%) 6057 (11.82%) 3722 (7.26%)  
 Middle 7088 (13.83%) 7635 (14.90%) 7815 (15.25%) 6243 (12.18%)  
 High 3699 (7.22%) 5490 (10.71%) 7956 (15.53%) 11777 (22.98%)  
 Not record 30131 (58.80%) 29637 (57.84%) 29416 (57.40%) 29503 (57.57%)  
Smoking status         <0.001
 Current smoker 1282 (2.50%) 2161 (4.22%) 3248 (6.34%) 4841 (9.45%)  
 Ever smoker 298 (0.58%) 507 (0.99%) 768 (1.50%) 893 (1.74%)  
 Never smoker 11551 (22.54%) 11271 (21.99%) 11017 (21.50%) 10307 (20.11%)  
 Not record 38114 (74.38%) 37305 (72.80%) 36211 (70.66%) 35204 (68.70%)  
Drinking status         <0.001
 Current drinker 108 (0.21%) 209 (0.41%) 386 (0.75%) 615 (1.20%)  
 Ever drinker 1059 (2.07%) 1850 (3.61%) 2572 (5.02%) 3234 (6.31%)  
 Never drinker 11964 (23.35%) 11880 (23.18%) 12075 (23.56%) 12192 (23.79%)  
 Not record 38114 (74.38%) 37305 (72.80%) 36211 (70.66%) 35204 (68.70%)  

TABLE 2: The results of univariate analysis.

In FIGURE 2, the Kaplan Meier curves of the cumulative risk of diabetes events stratified by TyG-BMI showed that the cumulative risk of diabetes gradually increases with the increase of TyG-BMI. There was a significant difference in the risk of diabetes between the TyG-BMI quartile groups (log-rank test P<0.0001).

Figure

Figure 2: Kaplan-Meier event-free survival curve. Kaplan-Meier analysis of incident of diabetes based on TyG-BMI quartiles (logrank, P<0.0001). Note: Equation Q1; Equation Q2; Equation Q3; Equation Q4.

■ The multivariate analysis of TyG-BMI with DM risk

We applied Cox proportional hazards models to evaluate group differences in the association between TyG-BMI and incident diabetes, and TABLE 3 showed the unadjusted and adjusted models. In the crude model, TyG-BMI was positively correlated with the incidence of diabetes (HR=1.026, 95% Confidence Interval (CI):1.026 to 1.027, P<0.00001). We could draw the same conclusion in the model I(minimally adjusted model, adjusted DBP, SBP, gender, age, drinking status, smoking status, family history) and model II(fully adjusted model, adjusted DBP, SBP, gender, age, drinking status, smoking status, family history, TC, HDL, LDL, ALT, AST, SCR). Model I (HR=1.022, 95% CI: 1.021 to 1.023, P<0.00001), model II (HR=1.023, 95% CI: 1.022 to 1.024, P<0.00001), respectively. We also performed a sensitivity analysis taking TyG-BMI as a categorical variable (quartile) at the same time and calculating P for trend. The result was consistent with that of TyG-BMI as a continuous variable (trend P<0.00001). In the fully adjusted model (model II), the risk of diabetes in the Q4 group increased by 10.261 times compared to the Q1 group, and the trend in the quartile was significant (trend P<0.00001).

Variable Crude model (HR,95%CI,P) Model I (HR,95%CI,P) Model II (HR,95%CI,P)
TyG-BMI 1.026 (1.026, 1.027)<0.00001 1.022 (1.021, 1.023)<0.00001 1.023 (1.022, 1.024)<0.00001
TyG-BMI(quartile)
Q1 Ref Ref Ref
Q2 2.967 (2.405, 3.660)<0.00001 2.020 (1.636, 2.495)<0.00001 2.049 (1.659, 2.532)<0.00001
Q3 8.089 (6.675, 9.801)<0.00001 4.195 (3.450, 5.102)<0.00001 4.437 (3.644, 5.403)<0.00001
Q4 23.876 (19.854, 28.714)<0.00001 10.562 (8.730, 12.779)<0.00001 11.261 (9.277, 13.668)<0.00001
P for trend <0.00001 <0.00001 <0.00001

TABLE 3: Relationship between TyG-BMI and the incident of diabetes in different models.

■ The analyses of the non-linear relationship

Since TyG-BMI was a continuous variable, we used the Generalized Additive Model (GAM) to identify the nonlinear relationship between TyG- BMI and the incidence of diabetes (adjusted gender, age, DBP, SBP, drinking status, smoking status, family history, TC, HDL, LDL, ALT, AST, Scr). In addition, by using two-piecewise linear regression and recursive algorithm, the inflection point of TyG-BMI was calculated as 232.416. The association between TyG-BMI and incident diabetes was positive on either side of the inflection point. The positive potency was slightly weaker on the right side (HR=1.016, 95%CI: 1.014 to 1.018, P<0.0001) of the inflection point than on the left (HR=1.029, 95%CI: 1.027 to 1.031, P<0.0001) (TABLE 4 and FIGURE 3).

  Incident of diabetes (HR,95%CI, P)
Fitting model by standard linear regression 1.023 (1.022, 1.024)<0.0001
Fitting model by two-piecewise linear regression  -
   Inflection point of TyG-BMI 232.416
   <232.416 1.029 (1.027, 1.031)<0.0001
   >232.416 1.016 (1.014, 1.018)<0.0001
P for log likelihood ratio test <0.001

TABLE 4: The result of two-piecewise linear regression model.

Figure 3: The non-linear relationship between TyG-BMI and incident of diabetes after adjusting for Age, Gender, SBP, DBP, Smoking Status, Drinking Status, Family History, TC, HDL, LDL, ALT, AST, SCR.

■ The results of subgroup analyses

TABLE 5 was the subgroup analysis for the association of TyG-BMI with diabetes incidence to explore other risks. The participants were divided into subgroups according to gender, age, HDL, LDL, SBP, DBP, family history of diabetes, drinking status, and smoking status. The association between TyG-BMI and diabetes risk was stable in family history, smoking status, and gender of patients (all P values for interaction>0.05). In contrast, we observed a number of interactions, including age, HDL, LDL, SBP, DBP, Drinking status (all P values of interaction<0.05). The relationship between TyG- BMI and diabetes was stronger in people with age 20-30(HR 1.029, 95%CI :1.024 to 1.035), age 30-40(HR 1.032, 95%CI :1.029 to 1.034), age 40-50(HR 1.029, 95%CI :1.027 to 1.031), HDL(high group) (HR=1.024, 95%CI: 1.022 to 1.026), SBP<140(HR=1.025, 95%CI:1.024 to 1.027), DBP<90(HR=1.024, 95%CI:1.023 to 1.025), current drinker(HR=1.031, 95%CI: 1.022 to 1.041) and ever drinker(HR=1.032, 95%CI: 1.027 to 1.037). In addition, the relationship between TyG-BMI and diabetes risk was weaker in the people with age 60-70 (HR 1.015, 95%CI :1.013 to 1.017), age ≥ 70(HR 1.013, 95%CI :1.011 to 1.016), HDL(low group) (HR=1.020, 95%CI: 1.018 to 1.021), LDL(middle group) (HR=1.022, 95%CI: 1.020 to 1.024), LDL(high group) (HR=1.021, 95%CI: 1.019 to 1.023), SBP ≥ 140(HR=1.017, 95%CI:1.015 to 1.019), DBP ≥ 90(HR=1.018, 95%CI:1.016 to 1.021) and never drinker (HR=1.022, 95%CI: 1.020 to 1.024).

Characteristic No. of participants Effect size(HR,95%CI,P)      P for interaction
Age(years)                                    <0.0001
20 to<30                      27301 1.029 (1.024, 1.035)<0.0001
30 to<40 80043 1.032 (1.029, 1.034)<0.0001  
40 to<50 43888 1.029 (1.027, 1.031)<0.0001  
50 to<60 29252 1.023 (1.021, 1.025)<0.0001  
60 to<70 17278 1.015 (1.013, 1.017)<0.0001  
≥ 70                                                                         7216 1.013 (1.011, 1.016)<0.0001  
Gender   0.9252
Male 112106 1.023 (1.022, 1.024) 0.9252  
Female 92872 1.023 (1.021, 1.024) 0.9252  
HDL(mmol/L)                                                                        0.0012
 Low 38229 1.020 (1.018, 1.021)<0.0001  
 Middle 38076 1.023 (1.021, 1.026)<0.0001  
 High 40019 1.024 (1.022, 1.026)<0.0001  
 Not recorded 88654 1.024 (1.023, 1.026)<0.0001  
LDL(mmol/L)                                                                  0.0222
 Low 38642 1.023 (1.021, 1.025) <0.0001  
 Middle 39236 1.022 (1.020, 1.024) <0.0001  
 High 39580 1.021 (1.019, 1.023) <0.0001  
 Not recorded 87520 1.024 (1.023, 1.026) <0.0001  
SBP(mmHg)     <0.0001
  <140 185128 1.025 (1.024, 1.027)<0.0001  
  ≥140 20958 1.017 (1.015, 1.019)<0.0001  
DBP(mmHg)     <0.0001
  <90 189459 1.024 (1.023, 1.025)<0.0001  
  ≥ 90 16627 1.018 (1.016, 1.021)<0.0001  
Smoking status   0.1151
Current smoker 11532 1.025 (1.022, 1.028)<0.0001  
Ever smoker 2466 1.029 (1.021, 1.036)<0.0001  
Never smoker 44146 1.023 (1.021, 1.025)<0.0001  
Not recorded 146834 1.022 (1.021, 1.023)<0.0001  
Drinking status   0.0002
Current drinker 1318 1.031 (1.022, 1.041)<0.0001
Ever drinker 8715 1.032 (1.027, 1.037)<0.0001  
Never drinker 48111 1.022 (1.020, 1.024)<0.0001  
Not recorded 146834 1.022 (1.021, 1.023)<0.0001  
Family history of diabetes   0.1543 <0.0001
   No 200788 1.023 (1.022, 1.024) <0.0001
   Yes 4190 1.020 (1.016, 1.024) <0.0001

TABLE 5: Effect size of TyG-BMI on diabetes in prespecified and exploratory subgroups.

Discussion

Our results suggest that TyG-BMI is associated with diabates. In view of the fact that IR acts as an core agent of T2DM and that diabetes is believed to be an metabolic disorder disease, these findings are not surprising. During 204978 person years of follow-up, we found a positive non-linear relationship between TyG-BMI and diabetes. Results remained essentially the same when we adjusted for gender, age, DBP, TC, LDL, HDL, AST, ALT, Scr, smoking, SBP. The inflection point value was 232.416, which was consistent in the direction before and after the inflection point, but the effect value was not completely consistent (Right (HR: 1.016, 95% CI: 1.014-1.018, P<0.0001), (left (HR: 1.029, 95% CI: 1.027-1.031, P<0.0001)). To be classified as an metabolic disease, the progress of the disorder must be slowed down or prevented by metabolic therapy. Our data indicate that TyG-BMI, as a new IR predictor index, can predict the occurrence of DM and provide important prognostic information.

The subgroup analysis showed that the association between TyG-BMI and diabetes risk were stronger in the following groups: age 20- 30(HR=1.029, 95%CI:1.024 to 1.035), age 30-40(HR=1.032, 95%CI:1.029 to 1.034), age 40-50(HR 1.029, 95%CI:1.027 to 1.031), high HDL(HR=1.024, 95%CI: 1.022 to 1.026 ), SBP<140(HR=1.025, 95%CI:1.024 to 1.027), DBP<90 (HR=1.024, 95%CI:1.023 to 1.025), current drinker(HR=1.031, 95%CI: 1.022 to 1.041), and ever drinker(HR=1.032, 95%CI: 1.027 to 1.037).

The main events in developing type 2 diabetes are decreased β cell function and insulin resistance. Insulin resistance can cause hyperinsulinemia, which in turn can aggravate insulin resistance [2]. It is generally accepted that insulin resistance is closely related to type 2 diabetes risk. Clinically, the gold standard of insulin resistance is the glucose clamp test, which is inconvenient and expensive. Although Homeostasis Model Assessment of Insulin Resistance (HOMA-IR) has a wide range of clinical applications, its application is limited due to its relatively high cost and low repeatability. Simental-Mendía Proposed the concept of the TyG index, which indicated that the TyG index could be used as an alternative index of insulin resistance in healthy subjects. Many studies have shown that the TyG index is a good alternative marker of insulin resistance [4,14-21]. Compared with the HOMA-IR index, the TyG index has higher sensitivity in recognizing insulin resistance [19]. TyG index is regular and easy to get, including FBG and TG, which have been associated with diabetes risk [8,9,18,22]. In addition to the TyG index, the relationship between obesity and diabetes is also well documented. The incidence rate of type 2 diabetes is lower in non-obese patients [23]. BMI is a simple, economical, and helpful indicator of general obesity. A cross- sectional study of the Taiwan population shows that TyG-BMI is an effective marker for early recognition of insulin resistance [3]. A recent study involving 511 participants showed that TyG-BMI is a stronger IR predictor than TyG- WC [3]. In a Nigerian cross-sectional study, in all 473 participants, TyG-BMI shows larger AUC for metabolic syndrome detection (0.838, 95% CI: 0.802–0.870) than TyG index (0.796, 95% CI: 0.757-0.831). After adjusting for age, gender, smoking, DBP and SBP, only the TyG index and TyG-BMI significantly predicted metabolic syndrome in men [24]. Previous studies have examined the association of the IR indicators such as BMI, TyG with DM risk and their performance in DM risk prediction [8,9,18,22,24]. These conclusions are similar to ours. Firstly, consistent with previous studies, TyG-BMI was positively correlated with incident diabetes. Our findings accord with those of three other studies [14-16].Then we further analyzed and discovered the curvilinear relationship. After adjusting for confounding factors in our study, the association between TyG-BMI and incident diabetes was nonlinear. The inflection point of GAM was 232.416 after adjusting for potential confounding factors (gender, age, DBP, LDL, TC, HDL, AST, ALT, Scr, smoking, SBP). We found that the inflection point had a stronger relationship on the left side of the inflection point. Therefore, controlling TyG-BMI is more valuable for reducing diabetes risk under the inflection point.

Our study has some strengths: (1) As far as we know, this is the first study to evaluate the association between TyG-BMI and diabetes events in the Chinese people based on our preprint outcome; (2) Compared with other researches, our sample size is relatively larger, which can better represent the Chinese people; (3) This study was a retrospective cohort study, which reduced selection and observation biases (4) We found the non-linear relationship and made a deeper discussion, and there are also more confounding factors for adjustment; (5) To make the results more robust, TyG-BMI was treated both as a continuous and categorical variable.

The study also has some potential limitations: (1) The data was from the Rich Healthcare Group in China, representing the Chinese population, and couldn’t be extended to other races and particular groups like pregnant women and children. (2) This research was a secondary analysis based on the published data, and variables were limited to the data of the original study. Other important variables such as medication history, hip circumference, hemoglobin A1C, physical activity, dietary factors were not included. (3) The incidence of diabetes may underestimate because of the study’s diabetes definition, which did not conduct a 2-hour oral glucose tolerance test. But for such a large cohort, improving participants’ oral glucose tolerance test is a considerable project. (4) The study did not differentiate diabetes types. But these conclusions may be more appropriate to type 2 diabetes, which accounts for about 90% of diabetic patients. (5) According to TyG-BMI, we only measured it at baseline, not measured over time. We can consider more variables and longer follow-up times in the future and adopt more sophisticated strategies.

Conclusion

The association between TyG-BMI and incident diabetes is positive and nonlinear after correcting the related confounding factors. The inflection point was 232.416. The relationship between TyG-BMI and diabetes is the most significant on the left side of the inflection point. TyG- BMI and incident diabetes had a non-linear positive relationship. Before and after TyG-BMI equals 232.416, the risk of diabetes increased by 2.9% and 1.6%, respectively, when TyG-BMI increased one unit. The findings emphasise the importance of baseline TyG-BMI as a measure for identification of patients at risk of diabetes early and improve their outcomes.

Declarations

■ Ethics approval and consent to participate

In the previously published article, Ying Chen has clearly stated that: the study was conducted in accordance with the Declaration of Helsinki. Informed consent was obtained from all Participants. Original study was approved by the Rich Healthcare Group Review Board, and the information was retrieved retrospectively.

■ Consent for publication

Not applicable.

■ Availability of data and materials

Data can be downloaded from ‘DATADRYAD’ database (www.Datadryad.org).

■ Competing interests

The authors declare that they have no competing interests.

■ Funding

This study was also supported in part by the International Cooperative Research Project of Shenzhen Municipal Science and Technology Innovation Council (accounts GJHZ20180416164814621).

■ Authors’ contributions

Fan Yang and Xiaohan Ding contributed to the study design, analyzed the data and drafted the manuscript. Zhuangsen CHEN, Yan Liao, Miaoling CHEN contributed to review the manuscript and took part in the interpretation and discussion of results. Weili YAO and Qian LIANG approved the final version and oversaw the progress of the project. Xinyu WANG and Haofei HU had full access to all the data in the study and had primary responsibility for the final content. All authors read and approved the final the manuscript.

Acknowledgement

Not applicable.

References