Triglyceride-Glucose Index as a Novel Biomarker in the Occurrence of Kidney Stones: A Cross-Sectional Population-Based Study.

BACKGROUND: Triglyceride-glucose (TyG) index has been considered as the reliable marker of insulin resistance (IR), which is one risk factor of kidney stone. This study aimed to evaluate the TyG index in the occurrence of kidney stones among the United States (US) population.
METHODS: Participants who received assessment were retrieved from the National Health and Nutrition Examination Survey (NHANES) between 2007 and 2018. The logistic regression analysis was conducted to assess the relationship between the TyG index and kidney stones occurrence. A 1:1 matched-pair analysis was established to optimize the bias in kidney stones by propensity score matching (PSM). The dose-response curve was performed to verify the association between the TyG index and risk of kidney stones.
RESULTS: Of the 14,158 eligible enrolled participants, 1346 (9.5%) self-reported a history of kidney stones. All participants were divided into two groups (high TyG index group and low TyG index group) based on the median TyG index. The dose-response curve exhibited a positive non-linear correlation between the TyG index and kidney stones risk. High TyG index was related to increased kidney stones occurrence in participants, with adjusted odds ratios (AOR) of 1.14 (95% confidence intervals (CI): 1.01-1.30, P = 0.038) compared with the low TyG index subgroup before PSM. After PSM, the risk of kidney stones was 19% higher in the high TyG group compared with the low TyG group (AOR = 1.19, 95% CI: 1.02-1.38, P = 0.026), and the dose-response curve still showed a positive association between TyG index and kidney stone risk.
CONCLUSION: The TyG index was independently associated with kidney stones and would be a novel biomarker in predicting occurrence for clinical decision.

Introduction

Kidney stones are the current most urinary stones in patients globally,1 with high prevalence and a worldwide increase over the past few decades. The prevalence of kidney stones in the USA had a nearly threefold relative increase over a 30-year period, with National Health and Nutrition Examination Survey (NHANES) data reflecting from 3.2% in 1976–1980 to 8.8% in 2007–2010.2 Moreover, the estimated prevalence is 7–13% in North America, 4% in South America, 5–10% in Europe and 1–19% in Asia.2 The recurrences are also common in a symptomatic kidney stone patient after the initial stone episode within the first 5 years,3 while the recurrence rate is up to 50% during the whole lifetime in both genders.4

Accounting for the increasing prevalence and high recurrence, kidney stone disease brings a tremendous economic burden for health systems, especially in a working-age population.5 As considered the most expensive urological disease, kidney stone disease has cost more than 10 billion dollars annually in health care system.6 Therefore, it is necessary to pay great attention to preventing the occurrence of kidney stones. Numerous systemic factors have been reported to be associated with an increased risk of kidney stones, suggesting that genetic, environmental, and nutritional factors play an important role in stone formation.7 In particular, kidney stones may be affected by modifiable factors, like dietary and lifestyle factors. Obesity, metabolic syndrome, and diabetes are identified as independent factors that contribute to increasing the stone risk.8 One of the underlying pathophysiologies of stone formation for these patients is thought to be insulin resistance (IR).9

The gold standard to assess IR is the hyperinsulinaemic-euglycaemic glucose clamp test.10 However, it is difficult to be popularized in clinical settings due to the time-consuming, complex and expensive characteristics. Recently several previous studies have proposed the triglyceride–glucose (TyG) index, calculated from fasting levels of triglycerides (TG) and glucose, as a simple, reliable and reproducible predictive indicator marker for the assessment of IR.11 Therefore, in this study, we aimed to evaluate the value of the TyG index in the occurrence of kidney stones among the United States (US) population.

Materials and Methods

Study Population

The baseline clinical data evaluated in this study were retrieved from NHANES from 2007 to 2018. As a cross-sectional survey, NHANES is conducted every other year in the US population by the Centers for Disease Control and Prevention (CDC). Study protocols in NHANES were reviewed and approved by the Institutional Review Board of the National Center for Health Statistics (NCHS), and consent from participants was signed during the survey. This study used previously collected deidentified data, which was deemed exempt from review by the Ethics Committee of the Affiliated Zhongda Hospital of Southeast University.

Our work assessed six consecutive 2-year survey cycles (2007–2008, 2009–2010, 2011–2012, 2013–2014, 2015–2016, 2017–2018), including the questionnaires of kidney stones. All the participants received the KIQ026 survey (Have you ever had kidney stones) for the assessment. The exclusion criteria were as follows: 1) unknown stone (n=25,163); 2) unknown triglyceride, high-density lipoprotein cholesterol (HDL-C) and low-density lipoprotein cholesterol (LDL-C) (n=19,950); 3) unknown blood glucose (n=54); 4) unknown body mass index (BMI) (n=181); 5) unknown uric acid (n=4); 6) unknown insulin or urine creatinine (n=332) (Figure 1).

Figure 1

Schematic flowchart of exclusion criteria in our study.

Data Collection and Definitions

The evaluation of kidney stones was conducted by the KIQ026 survey, which asked the participants whether ever had kidney stones or not. Participants who answered “No” to the “Have you ever had kidney stones?” were considered to have no history of kidney stones. Blood samples were collected for further detection. Clinical variables in this study included gender (male and female), age (20–39, 40–59 and 60+), race (non-Hispanic white, non-Hispanic black, Mexican American, other Hispanic and other), education (less than high school, high school or equivalent, college or above and other), marital status (married and unmarried/other), BMI, hypertension, diabetes, recreational activities (vigorous and moderate), kidney stones, urine creatinine, blood urea nitrogen, creatinine, uric acid, fasting glucose, total cholesterol, TG, HDL-C and LDL-C. The TyG index was determined by ln [(TG level mg/L × fasting glucose level mg/dL)/2], and the homeostasis model assessment for IR (HOMA-IR) was calculated as (fasting insulin [μIU/mL] × fasting glucose [mmol/L])/22.5. With a median of quartiles as 8.49, the TyG index was classified into low and high groups.

Statistical Analysis

All the continuous variables were recorded by mean ± standard deviation (SD) and the classified variables were presented by proportion. Chi-square analysis was applied in the assessment of clinical variables in this study. Logistic regression was conducted to calculate the adjusted odds ratios (AOR) with a 95% confidence interval (CI) for TyG-related kidney stones risk. Univariate and multivariate models were established and included different variables. The basic multivariate model consisted of TyG, age and gender. In the core multivariate model, the characteristics of TyG, age, gender, race, education levels, marital status, and BMI. Additionally, clinical variables of diabetes, hypertension, vigorous recreational activities, moderate recreational activities, urine creatinine, blood urea nitrogen, creatinine, and uric acid were adjusted in the extended multivariate model based on the core model. The 1:1 propensity score matching (PSM) was performed to normalize potential baseline confounders due to the imbalance existed in kidney stones occurrence. A restricted cubic spline was subsequently conducted to determine the dose–response relationship between the TyG index and kidney stones occurrence. Statistical analyses were conducted by Statistical Package for the Social Sciences software (version 24.0) with P-value <0.05 as statistically significant. The PSM and restricted cubic spline were established through R (version 3.5.1).

Results

A number of 14,158 eligible participants retrieved from NHANES from 2007 to 2018 were included in our study. For the eligible participants, there were 1346 (9.5%) reported have kidney stones. The baseline characteristics are represented in Table 1. Of the total participants, the majority were female (7357, 52.0%), 60+ years (4806, 33.9%), non-Hispanic white (5833, 41.2%), college or above education (7440, 52.5%), married (7358, 52.0%), without hypertension (9010, 63.6%) and prediabetes (11,974, 84.6%). The average BMI was 29.10 ± 6.91 kg/m2 in our eligible participants. Participants with kidney stones had higher levels of HOMA-IR, blood urea nitrogen, creatinine, uric acid, fasting glucose, TG, LDL-C, TyG, and urine creatinine than those without kidney stones.

Table 1

Baseline Characteristics of 14,158 Participants in NHANES 2007–2018 According to Kidney Stone Status Before PSM

Characteristic	Total	Kidney Stone		P value
		No	Yes
	No. (%)	No. (%)	No. (%)
Total patients	14,158	12,812 (90.5)	1346 (9.5)
Gender				<0.001
Male	6801 (48.0)	6061 (47.3)	740 (55.0)
Female	7357 (52.0)	6751 (52.7)	606 (45.0)
Age				<0.001
20–39 years	4608 (32.5)	4357 (34.0)	251 (18.6)
40–59 years	4744 (33.5)	4278 (33.4)	466 (34.6)
60+ years	4806 (33.9)	4177 (32.6)	629 (46.7)
Race				<0.001
Non-Hispanic white	5833 (41.2)	5113 (39.9)	720 (53.5)
Non-Hispanic black	2847 (20.1)	2685 (21.0)	162 (12.0)
Mexican American	2165 (15.3)	1986 (15.5)	179 (13.3)
Other Hispanic	1562 (11.0)	1407 (11.0)	155 (11.5)
Other	1751 (12.4)	1621 (12.7)	130 (9.7)
Education				0.650
Less than high school	3516 (24.8)	3171 (24.8)	345 (25.6)
High school or equivalent	3192 (22.5)	2887 (22.5)	305 (22.7)
College or above	7440 (52.5)	6744 (52.6)	696 (51.7)
Other	10 (0.1)	10 (0.1)	0 (0.0)
Marital status				<0.001
Married	7358 (52.0)	6471 (51.3)	787 (58.5)
Unmarried/Other	6800 (48.0)	6241 (48.7)	559 (41.5)
BMI (kg/m2)	29.10, 6.91	28.94, 6.87	30.63, 7.16	<0.001
Hypertension				<0.001
Yes	5148 (36.4)	4485 (35.0)	663 (49.3)
No/Unknown	9010 (63.6)	8327 (65.0)	683 (50.7)
Diabetes				<0.001
Yes	1836 (13.0)	1535 (12.0)	301 (22.4)
Prediabetes	11,974 (84.6)	10,977 (85.7)	997 (74.1)
No/Unknown	348 (2.5)	300 (2.3)	48 (3.6)
Vigorous recreational activities				<0.001
Yes	3030 (15.0)	2828 (22.1)	202 (15.0)
No	11,128 (78.6)	9984 (77.9)	1144 (85.0)
Moderate recreational activities				<0.001
Yes	5654 (39.9)	5163 (40.3)	491 (36.5)
No	8504 (60.1)	7649 (59.7)	855 (63.5)
HOMA-IR	4.02, 6.79	3.92, 6.76	4.90, 6.98	<0.001
Blood urea nitrogen	13.64, 5.87	13.52, 5.82	14.85, 6.19	<0.001
Creatinine	0.88, 0.38	0.88, 0.39	0.91, 0.32	0.004
Uric acid	5.49, 1.44	5.47, 1.43	5.64, 1.47	<0.001
Fasting glucose	102.93, 34.07	102.17, 33.42	110.16, 38.96	<0.001
TC	190.89, 40.60	191.14, 40.70	188.46, 39.56	0.021
TG	116.29, 65.46	115.09, 65.06	127.65, 68.14	<0.001
HDL-C	54.26, 15.95	54.59, 16.03	51.12, 14.72	<0.001
LDL-C	113.37, 35.49	113.53, 35.59	111.81, 34.41	0.090
TyG	8.52, 0.64	8.50, 0.64	8.68, 0.64	<0.001
Urine creatinine	125.76, 77.00	125.43, 77.33	128.88, 73.75	0.118

Notes: For categorical variables, P values were analyzed by chi-square tests; For continuous variables, the t-test for slope was used in generalized linear models; All the continuous variables were recorded by mean ± standard deviation (SD).

Abbreviations: PSM, propensity score matching; BMI, body mass index; HOMA-IR, homeostasis model assessment for IR; TC, total cholesterol; TG, triglycerides; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; TyG, triglyceride–glucose.

All participants were divided into two groups (high TyG index group and low TyG index group) based on the median TyG index. The chi-square analysis results indicated that all included variables (except urinary creatinine) were related with the TyG index (All P < 0.001) and a higher proportion (11.6%) of kidney stones in the high TyG index group (Table 2). After included gender, age, race, education, marital status, BMI, hypertension, diabetes, urine creatinine, blood urea nitrogen, creatine, uric acid, vigorous and moderate recreational activities for adjustment, the dose–response curve exhibited the positive non-linear correlation between the TyG index and kidney stones risk (Figure 2).

Table 2

Baseline Characteristics of 14,158 Participants in NHANES 2007–2018 According to TyG Status Before PSM

Characteristic	Total	TyG Index		P value
		Low (<8.49)	High (≥8.49)
	No. (%)	No. (%)	No. (%)
Total patients	14,158	7080 (50.0)	7078 (50.0)
Gender				<0.001
Male	6801 (48.0)	3047 (43.0)	3754 (53.0)
Female	7357 (52.0)	4033 (57.0)	3324 (47.0)
Age				<0.001
20–39 years	4608 (32.5)	2906 (41.0)	1702 (24.0)
40–59 years	4744 (33.5)	2206 (31.2)	2538 (35.9)
60+ years	4806 (33.9)	1968 (27.8)	2838 (40.1)
Race				<0.001
Non-Hispanic white	5833 (41.2)	2790 (39.4)	3043 (43.0)
Non-Hispanic black	2847 (20.1)	1875 (26.5)	972 (13.7)
Mexican American	2165 (15.3)	847 (12.0)	1318 (18.6)
Other Hispanic	1562 (11.0)	672 (9.5)	890 (12.6)
Other	1751 (12.4)	896 (12.7)	855 (12.1)
Education				<0.001
Less than high school	3516 (24.8)	1456 (20.6)	2060 (29.1)
High school or equivalent	3192 (22.5)	1549 (21.9)	1643 (23.2)
College or above	7440 (52.5)	4073 (57.5)	3367 (47.6)
Other	10 (0.1)	2 (0.0)	8 (0.1)
Marital status				<0.001
Married	7358 (52.0)	3417 (48.3)	3941 (55.7)
Unmarried/Other	6800 (48.0)	3663 (51.7)	3137 (44.3)
BMI (kg/m2)	29.10, 6.91	27.50, 6.69	30.70, 6.76	<0.001
Hypertension				<0.001
Yes	5148 (36.4)	2010 (28.4)	3138 (44.3)
No/Unknown	9010 (63.6)	5070 (71.6)	3940 (55.7)
Diabetes				<0.001
Yes	1836 (13.0)	386 (5.5)	1450 (20.5)
Prediabetes	11,974 (84.6)	6563 (92.7)	5411 (76.4)
No/Unknown	348 (2.5)	131 (1.9)	217 (3.1)
Vigorous recreational activities				<0.001
Yes	3030 (15.0)	1914 (27.0)	1116 (15.8)
No	11,128 (78.6)	5166 (73.0)	5962 (84.2)
Moderate recreational activities				<0.001
Yes	5654 (39.9)	3079 (43.5)	2575 (36.4)
No	8504 (60.1)	4001 (56.4)	4503 (63.6)
Kidney stone				<0.001
Yes	1346 (9.5)	523 (7.4)	823 (11.6)
No	12,812 (90.5)	6557 (92.6)	6255 (88.4)
HOMA-IR	4.72, 8.59	2.44, 3.04	5.59, 8.84	<0.001
Blood urea nitrogen	13.64, 5.87	13.11, 5.39	14.17, 6.27	<0.001
Creatinine	0.88, 0.38	0.86, 0.35	0.90, 0.41	<0.001
Uric acid	5.49, 1.44	5.14, 1.34	5.83, 1.45	<0.001
Fasting glucose	102.93, 34.07	91.68, 12.68	114.17, 43.68	<0.001
TC	190.89, 40.60	181.38, 36.93	200.40, 41.86	<0.001
TG	116.29, 65.46	70.03, 21.37	162.56, 61.91	<0.001
HDL-C	54.26, 15.95	60.31, 16.43	48.21, 12.86	<0.001
LDL-C	113.37, 35.49	107.07, 32.04	119.67, 37.59	<0.001
TyG	8.52, 0.64	8.01, 0.35	9.02, 0.42	<0.001
Urine creatinine	125.76, 77.00	126.66, 79.97	124.86, 73.90	0.163

Notes: For categorical variables, P values were analyzed by chi-square tests; For continuous variables, the t-test for slope was used in generalized linear models; All the continuous variables were recorded by mean ± standard deviation (SD).

Abbreviations: PSM, propensity score matching; BMI, body mass index; HOMA-IR, homeostasis model assessment for IR; TC, total cholesterol; TG, triglycerides; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; TyG, triglyceride–glucose.

Figure 2

The dose–response analysis between triglyceride–glucose (TyG) index and kidney stones occurrence before propensity score matching (PSM).

Univariate and multivariate logistic regression were used to analyze the ability of TyG index and HOMA-IR to predict the risk of kidney stones. The results indicated higher TyG index was positively related with kidney stones risk in the univariate model (high vs low: AOR = 1.65, 95% CI = 1.47–1.85, P < 0.001), basic model (high vs low: aOR = 1.43, 95% CI = 1.27–1.61, P < 0.001), core model (high vs low: AOR = 1.21, 95% CI = 1.07–1.36, P = 0.003) and extended model (high vs low: AOR = 1.14, 95% CI = 1.01–1.30, P = 0.038) (Table 3). However, HOMA-IR was not statistically associated with the risk of kidney stones in the extended model.

Table 3

Adjusted Odds Ratios for Associations Between the TyG and the Risk of Kidney Stone in NHANES 2007–2018

Characteristic	Univariate Analysis		Basic Model		Core Model		Extended Model
	aOR (95% CI)	P value	aOR (95% CI)	P value	aOR (95% CI)	P value	aOR (95% CI)	P value
Before PSM
TyG
Low group	1.00 (Reference)	<0.001	1.00 (Reference)	<0.001	1.00 (Reference)	0.003	1.00 (Reference)	0.038
High group	1.65 (1.47–1.85)		1.43 (1.27–1.61)		1.21 (1.07–1.36)		1.14 (1.01–1.30)
HOMA-IR
Low group	1.00 (Reference)	<0.001	1.00 (Reference)	<0.001	1.00 (Reference)	0.015	1.00 (Reference)	0.145
High group	1.51 (1.35–1.70)		1.41 (1.25–1.58)		1.18 (1.03–1.34)		-
After PSM
TyG
Low group	1.00 (Reference)	0.008	1.00 (Reference)	0.008	1.00 (Reference)	0.008	1.00 (Reference)	0.026
High group	1.24 (1.06–1.44)		1.24 (1.06–1.44)		1.24 (1.06–1.44)		1.19 (1.02–1.38)
HOMA-IR
Low group	1.00 (Reference)	0.001	1.00 (Reference)	0.001	1.00 (Reference)	0.252	1.00 (Reference)	0.275
High group	1.30 (1.11–1.52)		1.30 (1.11–1.52)		-		-

Notes: Adjusted covariates: Basic model: age, gender; Core model: basic model plus race, education levels, marital status, and BMI; Extended model: core model plus diabetes, hypertension, vigorous recreational activities, moderate recreational activities, blood urea nitrogen, creatinine, and uric acid.

Abbreviations: CI, confidence interval; aOR, adjusted odds ratio; PSM, propensity score matching.

According to the significant difference in the occurrence of kidney stones (9.5% vs 90.5%), a 1:1 matched cohort analysis was conducted to minimize the potential bias through PSM (Figure 3). We confirmed 2614 participants in the matched cohort through PSM. The chi-square analysis results evidenced that significant differences remained in most variables in the high TyG index group and low TyG index group (Table 4). Similarly, multivariate logistic regression analysis revealed that the risk of kidney stones was 19% higher in the high TyG group compared with the low TyG group (AOR = 1.19, 95% CI: 1.02–1.38, P = 0.026) in the extended model (Table 3), and the dose–response curve still showed a positive association between TyG index and kidney stone risk (Figure 4).

Figure 3

The standardized mean difference (SMD) results of different variables after PSM.

Table 4

Baseline Characteristics of 2614 Participants in NHANES 2007–2018 After PSM

Characteristic	Total	TyG Index		P value
		Low (<8.49)	High (≥8.49)
	No. (%)	No. (%)	No. (%)
Total patients	2614	1089 (41.7)	1525 (58.3)
Gender				<0.001
Male	1450 (55.5)	542 (49.8)	908 (59.5)
Female	1164 (44.5)	547 (50.2)	617 (40.5)
Age				<0.001
20–39 years	533 (20.4)	303 (27.8)	230 (15.1)
40–59 years	852 (32.6)	348 (32.0)	504 (33.0)
60+ years	1229 (47.0)	438 (40.2)	791 (51.9)
Race				<0.001
Non-Hispanic white	1364 (52.2)	551 (50.6)	813 (53.3)
Non-Hispanic black	357 (13.7)	216 (19.8)	141 (9.2)
Mexican American	340 (13.0)	113 (10.4)	227 (14.9)
Other Hispanic	293 (11.2)	107 (9.8)	186 (12.2)
Other	260 (9.9)	102 (9.4)	158 (10.4)
Education				<0.001
Less than high school	659 (25.2)	233 (21.4)	426 (27.9)
High school or equivalent	598 (22.9)	244 (22.4)	354 (23.2)
College or above	1356 (51.9)	611 (56.1)	745 (48.9)
Other	1 (0.0)	1 (0.1)	0 (0.0)
Marital status				0.053
Married	1529 (58.5)	613 (56.3)	916 (60.1)
Unmarried/Other	1085 (41.5)	476 (43.7)	609 (39.9)
BMI (kg/m2)	29.92, 6.95	28.10, 6.72	31.23, 6.82	<0.001
Hypertension				<0.001
Yes	1264 (48.4)	429 (39.4)	835 (54.8)
No/Unknown	1350 (51.6)	660 (60.6)	690 (45.2)
Diabetes				<0.001
Yes	551 (21.1)	98 (9.0)	453 (29.7)
Prediabetes	1993 (76.2)	966 (88.7)	1027 (67.3)
No/Unknown	70 (2.7)	25 (2.3)	45 (3.0)
Vigorous recreational activities				<0.001
Yes	407 (15.6)	234 (21.5)	173 (11.3)
No	2207 (84.4)	855 (78.5)	1352 (88.7)
Moderate recreational activities				0.002
Yes	934 (35.7)	426 (39.1)	508 (33.3)
No	1680 (64.3)	663 (60.9)	1017 (66.7)
Kidney stone				0.008
Yes	1307 (50.0)	729 (47.8)	578 (53.1)
No	1307 (50.0)	796 (52.2)	511 (46.9)
HOMA-IR	4.72, 8.59	2.48, 2.57	6.33, 10.75	<0.001
Blood urea nitrogen	14.64, 6.52	14.15, 6.51	15.00, 6.50	<0.001
Creatinine	0.92, 0.37	0.91, 0.43	0.92, 0.32	0.001
Uric acid	5.64, 1.46	5.31, 1.43	5.88, 1.43	0.234
Fasting glucose	109.15, 38.29	93.45, 13.57	120.36, 45.61	<0.001
TC	189.33, 40.75	181.85, 37.63	194.66, 42.05	<0.001
TG	124.61, 67.13	71.42, 20.73	162.60, 62.88	<0.001
HDL-C	52.12, 15.31	59.30, 16.05	47.00, 12.44	<0.001
LDL-C	112.29, 35.08	108.29, 32.05	115.14, 36.84	<0.001
Urine creatinine	126.50, 72.93	123.18, 75.91	128.86, 70.67	0.050

Notes: For categorical variables, P values were analyzed by chi-square tests; For continuous variables, the t-test for slope was used in generalized linear models; All the continuous variables were recorded by mean ± standard deviation (SD).

Abbreviations: PSM, propensity score matching; BMI, body mass index; HOMA-IR, homeostasis model assessment for IR; TC, total cholesterol; TG, triglycerides; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; TyG, triglyceride–glucose.

Figure 4

The dose–response analysis between triglyceride–glucose (TyG) index and kidney stones occurrence after PSM.

Discussion

Kidney stones as the common occurring disease pose a significant health burden in public.5 This work introduced the relationship between the TyG index and kidney stones occurrence in this cross-sectional population-based study. The high TyG index was examined to be positively associated with the kidney stones occurrence in the eligible participants. Meanwhile, the positive dose–response relationship was also proved in the 1:1 matched cohort after PSM. All the statistical results provided the favorable evidences in confirmation of the positive association existed in the TyG index and kidney stones.

Recent reports had revealed the IR as the risk factor contributing to kidney stones occurrence.12,13 Together with individual habits, genetic factors and environmental conditions, metabolic imbalance was considered as the major risk factor for kidney stones. In practice, IR would reduce urinary pH and further damage the ability of generating renal ammonium in response to an acid load.12 Kim et al14 conducted a retrospective study based on 278,628 adults for the evaluation of glycemic status and IR in nephrolithiasis. The results indicated that IR was positively related with elevated risk of nephrolithiasis, which in consistent with our fundings. Previous research evidences from Weinberg et al15 also confirmed the potential relationship between IR and kidney stones. It was supposed that uric acid excretion functions as the mechanism that may illustrate the stone formation.16 All the results supported that IR was the risk factor for kidney stones and would contribute to the occurrence prediction.

The TyG index, a novel indicator determined by TG and fasting glucose, had been recently accepted as an intuitive and reliable predictor of IR, assisting in clinical decision.17–19 Different from HOMA-IR, the TyG index was universally adaptable in all patients due to the calculation without quantification of insulin.20,21 It was well distinguished that an increased TyG index would contributing to diagnosing T2DM and predicting related adverse outcomes.22,23 Kidney as a metabolic organ is easy to be damaged by glucose metabolism disorders and diabetic nephropathy has been considered as a huge burden in population. Whether the TyG index would contribute to predicting renal diseases remains unclear. In a study of 11,712 participants, Okamura et al24 investigated the predicting performance of the TyG index in chronic kidney disease (CKD) risks. And it was reported that the TyG index predicts significant risks for CKD in population with a hazard ratio (HR) of 1.32 for men, and 1.50 for women. The results promoted us to explore the application of the TyG index in renal diseases, especially the ones related with IR. Furthermore, our study found that the TyG index was a better predictor of kidney stone risk than HOMA-IR.

In this study, we estimated the relationship between the TyG index and kidney stones through various analyses. Meanwhile, we established the dose–response relationships before and after the PSM for the evaluation of the TyG index predicting performance in stone occurrence. To our knowledge, this work was the first attempt to assess the function of the TyG index in kidney stones occurrence. The aORs (95% CI) of high TyG group were 2.32 (2.06–2.66) and 2.37 (2.14–2.66) in extended model before/after PSM among the 14,490 eligible participants. Our results revealed that TyG index could function as the predictor of kidney stones occurrence with an excellent performance. And this work would expand the application of the TyG index in renal diseases prediction and contribute to individualized therapy and clinical decision.

Nevertheless, this study has limitations as follows. Firstly, this survey as a retrospective study, could not avoid the happen of bias, such as recording, inquiring, etc. Secondly, the data retrieved from NHANES only involved the population from America which was unable to represent other regions including Asia, Europe, etc. Thirdly, the detailed clinical variables were not disclosed in the database including the personal drug history and kidney stones types, which remaining further exploration. Therefore, the prospective international multicenter studies desired to establish for the evaluation of our results. Despite these limitations, there are strengths of the study including proposing a novel index for kidney stones occurrences and revealing the relationship between the TyG index and kidney stones, which assisting the accurate clinical previsions and treatments.

Conclusion

This study reported high TyG index was examined to be positively associated with the kidney stones occurrence, and TyG index would be a novel biomarker in predicting occurrence for individualized therapy and clinical decision.