Study design and participants
Participants enrolled in this study were from the CHARLS, which is an ongoing longitudinal survey for accessing the social, economic and health status of community residents aged 45 years or older in China. Details of the study design of CHARLS have been described elsewhere [16]. Briefly, the CHARLS adopted a multistage probability sampling and investigated 17708 individuals in 28 provinces through random selection of 10257 households to cover the overall population in China in the first wave (W1, 2011-2012), making the sample nationally representative. The response rate was up to 81% in the baseline survey. Well-trained interviews collected sociodemographic, physical, biological, and health-related information of study participants complying with the standard study plan. To date, the follow-up surveys have been conducted twice, including the second wave (W2) in 2013 and the third wave (W3) in 2015. For the present study, 10111 individuals with laboratory measurement at W1 were initially enrolled. Participants aged less than 45y, missing TG, FBG, and glycated hemoglobin were excluded, leaving 9663 participants. Of those, 1760 with diabetes at baseline were further excluded. Next, we excluded subjects who died (n=103) or lost follow-up (n=372) in the subsequent waves (W2 and W3) of the study. Finally, 7428 participants were enrolled for the final analysis (Figure 1). The Ethics Review Committee of Peking University approved CHARLS (IRB00001052–11015) and all participants gave informed consent before participation.
Data collection and definitions
After taking a rigorous training program, the researchers visited and interviewed participants in their homes using computer-assisted personal interview technology to collect sociodemographic information [including age, gender, education level (primary school or lower, secondary school, and higher), and marriage status (current married or not)], health behavior (including current smoking and drinking), medical history [including self-reported hypertension, diabetes, and cardiovascular disease (CVD)] and medication usage (including antihypertensive drugs, antidiabetic drugs, and lipid-lowering drugs).
Anthropometric indicators included systolic blood pressure (SBP), diastolic blood pressure (DBP), body mass index (BMI, calculated as weight in kilograms divided by the square of height in meters), and waist circumference (WC). Blood pressure were measured three times in a seated position by trained nurses using the HEM-7200 electronic monitor (Omron, Dalian, Japan). Hypertension was defined as SBP≥140 mmHg or DBP ≥ 90 mmHg or self-reported prior diagnosis of hypertension by a doctor or using antihypertensive drugs in the past two weeks [17]. Height and WC were accurate to 0.1 cm and weight was accurate to 0.1 kg.
For biomarkers assessment, the CHARLS researcher collected fasting blood samples from every participant. These samples were transported from all over the country to Beijing and were stored at minus 80℃ at the Chinese Center for Disease Control and Prevention. The biological determination of FBG, Hemoglobin A1c, TG, total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), and high-density lipoprotein cholesterol (HDL-C) were conducted by professional staff. TyG index was calculated as ln [TG (mg/dl) × FBG (mg/dl)/2] [18]. At baseline and follow-up, diabetes was defined as FBG > 125 mg/dL or Hemoglobin A1c > 6.5% or self-reported prior diagnosis of hypertension by a doctor or using antidiabetic medications. Non-diabetic participants whose FBG was at 100–125 mg/dL or Hemoglobin A1c was at 5·7–6·4% were classified as having prediabetes [19]. Participants without diabetes and prediabetes were normoglycemic.
Statistical Analysis
Data were presented as mean standard deviation (SD) for continuous variables and percentage for categorical variables. Baseline characteristics between groups according to quartiles of TyG index (Q1, Q2, Q3, Q4) or developing diabetes or not were reported and compared using the One-Way ANOVA, Kruskal-Wallis H test and chi-square tests, as appropriate. We initially conducted Cox proportional hazards models to estimate HR with 95% confidence interval (CI) of diabetes for per SD increase as well as quartiles of TyG index. Three models were fitted. Model 1 was univariate. Models adjusted for age and gender. Fully adjusted model (Model 3) incorporated covariates including age, gender, education, marriage, smoking, drinking, BMI, WC, SBP, history of hypertension, history of CVD, and usage of lipid-lowering drugs. Next, the shape of association between TyG index and incident diabetes was examined by multivariate adjusted Cox restricted cubic spline regression model. We choose three knots at quartiles 25th, 50th, and 75th. Finally, we conducted subgroup analyses (multivariate Cox proportional hazards models) including age (< 65 or ≥ 65 years), gender (male or female), BMI (< 25 or ≥ 25 kg/m2), and glycemic status (normoglycemia or prediabetes). P < 0.05 was considered statistically significant. R version 3.3.2 (R Foundation for Statistical Computing, Vienna, Austria) was used for all statistical analyses.