Multiple Biomarkers Improved Prediction for the Risk of Type 2 Diabetes Mellitus in Singapore Chinese Men and Women

Article information

Diabetes Metab J. 2020;44(2):295-306
Publication date (electronic) : 2019 November 22
doi : https://doi.org/10.4093/dmj.2019.0020
1Health Services and Systems Research, Duke-NUS Medical School, Singapore.
2Saw Swee Hock School of Public Health, National University of Singapore and National University Health System, Singapore.
3UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA, USA.
4Department of Epidemiology, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, PA, USA.
5Department of Epidemiology and Biostatistics, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China.
Corresponding author: An Pan. Department of Epidemiology and Biostatistics, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, 13 Hangkong Road, Qiaokou District, Wuhan 430060, China. panan@hust.edu.cn
Corresponding author: Yeli Wang. Health Services and Systems Research, Duke-NUS Medical School, 8 College Road, 169857, Singapore. yeli.wang@duke-nus.edu.sg
Received 2019 January 27; Accepted 2019 March 19.

Abstract

Background

Multiple biomarkers have performed well in predicting type 2 diabetes mellitus (T2DM) risk in Western populations. However, evidence is scarce among Asian populations.

Methods

Plasma triglyceride-to-high density lipoprotein (TG-to-HDL) ratio, alanine transaminase (ALT), high-sensitivity C-reactive protein (hs-CRP), ferritin, adiponectin, fetuin-A, and retinol-binding protein 4 were measured in 485 T2DM cases and 485 age-and-sex matched controls nested within the prospective Singapore Chinese Health Study cohort. Participants were free of T2DM at blood collection (1999 to 2004), and T2DM cases were identified at the subsequent follow-up interviews (2006 to 2010). A weighted biomarker score was created based on the strengths of associations between these biomarkers and T2DM risks. The predictive utility of the biomarker score was assessed by the area under receiver operating characteristics curve (AUC).

Results

The biomarker score that comprised of four biomarkers (TG-to-HDL ratio, ALT, ferritin, and adiponectin) was positively associated with T2DM risk (P trend <0.001). Compared to the lowest quartile of the score, the odds ratio was 12.0 (95% confidence interval [CI], 5.43 to 26.6) for those in the highest quartile. Adding the biomarker score to a base model that included smoking, history of hypertension, body mass index, and levels of random glucose and insulin improved AUC significantly from 0.81 (95% CI, 0.78 to 0.83) to 0.83 (95% CI, 0.81 to 0.86; P=0.002). When substituting the random glucose levels with glycosylated hemoglobin in the base model, adding the biomarker score improved AUC from 0.85 (95% CI, 0.83 to 0.88) to 0.86 (95% CI, 0.84 to 0.89; P=0.032).

Conclusion

A composite score of blood biomarkers improved T2DM risk prediction among Chinese.

INTRODUCTION

The global prevalence of diabetes has increased explosively in the past three decades, particularly in Asian countries that have experienced rapid transitions in nutritional status and economic growth [1]. According to the International Diabetes Federation, a total of 425 million people had diabetes worldwide in 2017, with China ranking top (114.4 million) and Southeast Asia contributing to around one-fifth of world's population with diabetes [2]. Type 2 diabetes mellitus (T2DM) accounts for around 90% of diabetes cases [3], and several clinical trials have shown that intensive lifestyle modifications focusing on healthy diet and physical activity can prevent T2DM development [45]. Therefore, it is vital to identify individuals at higher T2DM risk so that intensive lifestyle and dietary interventions can be introduced to lower their risks of developing T2DM.

The etiology of T2DM involves complex yet poorly-understood interactions between genetic and environmental factors. Body mass index (BMI), family history of T2DM, fasting blood glucose and glycosylated hemoglobin (HbA1c) are clinical risk factors that are commonly used for identifying people at higher risk of developing T2DM [6]. However, the utility of these factors as predictors of T2DM in the general population is not high; for example, it has been estimated that less than 20% of obese Caucasians will progress to T2DM over the next few years [7], and genotype scores provide minimal improvement in the risk prediction on top of common risk factors [89]. In addition, by the time abnormal glucose regulation has been detected by blood measurements, a significant proportion of such individuals have already had complications due to the chronicity of T2DM progression [10]. In addition, diagnostic markers such as high blood glucose and HbA1c levels are outcome indicators of T2DM and provide limited insights into the underlying etiology [11].

Recent studies have shown that dysregulation of several metabolic pathways such as a failure of intracellular lipid homeostasis [12], abnormal adipocyte signaling [13], dyslipidemia [14], inflammation [15], and iron overload [16] preceded the development of T2DM. In addition, some biomarkers implicated in the abovementioned pathways, such as alanine transaminase (ALT) [17], fetuin-A [18], retinol-binding protein 4 (RBP4) [1920], adiponectin [21], triglyceride (TG) [2223], high density lipoprotein cholesterol (HDL-C) [2223], high-sensitivity C-reactive protein (hs-CRP) [24], and ferritin [25], have been found to be significantly associated with T2DM risk. Furthermore, epidemiological studies in Western populations have shown that combination of these biomarkers could improve T2DM prediction over traditional clinical risk factors [11262728293031], while such evidence is scarce in Asian populations. To our best knowledge, only one cross-sectional study has conducted joint analysis of multiple biomarkers for identifying T2DM in a Chinese population [32]. Compared with the Western populations, Asians develop T2DM at lower BMIs, have lower β-cell function and are more insulin resistant [33]. Thus, it is of scientific interest to assess whether joint analysis of multiple biomarkers could improve T2DM prediction in a prospective cohort study among Asian populations.

Thus, we conducted this nested case-control study within the Singapore Chinese Health Study to create a weighted biomarker score, and prospectively evaluate the association between the biomarker score and T2DM risk with adjustment for established T2DM risk factors. We further explored whether this biomarker score could improve T2DM prediction in this Chinese population.

METHODS

Study population

The Singapore Chinese Health Study is a population-based cohort study. At baseline between 1993 and 1998, a total of 63,257 Chinese living in Singapore (age range, 45 to 74 years) were recruited. At recruitment interviews, trained research staff collected information on medical history, diet intakes and lifestyle habits using structured questionnaires. From 1999 to 2004, we conducted the first follow-up interview by telephone, and re-interviewed 52,322 subjects to update their information on medical history and selected lifestyle habits. In addition, 32,535 participants consented and donated bio-specimens for the purpose of research. Subsequently, we carried out the second follow-up via phone from 2006 to 2010 and re-contacted 39,528 subjects, including 25,477 who donated their bio-specimens at the first follow-up. The Institutional Review Boards at the University of Pittsburgh and the National University of Singapore approved the study protocol (NUS-IRB Ref Code: 06-027), and all subjects signed the consent forms after hearing the full explanation of the purpose and nature of all procedures used. In addition, all procedures complied with the Declaration of Helsinki. The detailed description of the study design has been previously [34].

Ascertainment of T2DM and other covariates

At baseline and two follow-ups, we asked participants if a doctor had told them that they had diabetes. If participants answered yes, we also asked for the age of first diagnosis. We further validated the accuracy of the self-reported T2DM data in a separate study [35].

At recruitment, we collected information on age, body height and weight, education levels, smoking status, alcohol intakes, weekly activity levels, and history of hypertension, cardiovascular disease (CVD), stroke, and cancer. BMI was calculated using body weight (kg) divided by square of height (m).

Establishment of nested case-control study

Among subjects without physician-diagnosed T2DM, CVD and cancer during blood collection (1999 to 2004), 571 subjects reported that they were diagnosed with T2DM at the second follow-up visit (2006 to 2010), and thus were considered as incident T2DM cases for the current study. In the remaining subjects without T2DM at the second follow-up and had HbA1c <6.0%, controls were randomly chosen and matched with cases on a 1:1 ratio for age (±3 years), date of blood collection (±6 months), sex, and dialect group (Cantonese, Hokkien). The matching on dialect group was to reduce the confounding by culture, lifestyle, and genetic factors that may associated with different dialect group. A total of 86 pairs were excluded because either cases or controls in these pairs lacked the information for ferritin measurements. The final sample size for the current study was 485 case-control pairs. The flowchart of the study design is presented in Supplementary Fig. 1.

Laboratory measurements

We collected 20-mL morning blood from each consented participant at random time slots, and transported it to the laboratory with ice immediately. We separated the blood into different components including plasma, serum, red blood cells and buffy coat, and put all components in −80℃ freezers for long-term storage. For the current study, plasma samples from each case-control pair were analyzed in the same batch at the National University Hospital Reference Laboratory, Singapore.

Plasma levels of total cholesterol, HDL-C, TG, ALT, hs-CRP were determined by the colorimetric method (AU5800 Analyzer; Beckman Coulter, Brea, CA, USA). In addition, ferritin levels were measured by the Sandwich enzyme-linked immunosorbent assay (ELISA) method (Dxl 800 Analyzer; Beckman Coulter). Furthermore, levels of fetuin-A, adiponectin, and RBP4 were both determined via Sandwich ELISA (Bio-Rad Laboratories, Hercules, CA, USA). Moreover, HbA1c from the red blood cells was measured by the HPLC method (Bio-Rad Laboratories) in red blood cells. The within- and between-assay coefficients of variation of all measures were <10%.

Statistical analysis

Conditional logistic regression models were used to compute the odds ratio (OR) of each biomarker (quartiles) in association with T2DM risk. Model 1 was adjusted for age (years), education levels (none, primary, secondary, or higher), alcohol intakes (never or monthly, weekly, or daily), weekly physical activity levels (<0.5, 0.5 to 3.9, and ≥4.0 hr/wk), smoking status (never, ever smoker), BMI (kg/m2), history of hypertension and fasting status. Of note, since cases and controls were matched on age with a range of ±3 years, we adjusted for age as a continuous variable to further account for the residual confounding. Model 2 was further adjusted for all the blood biomarkers measured for this study. As a result, four biomarkers (TG-to-HDL ratio, ALT, ferritin, and adiponectin) were significantly associated with T2DM risk and were selected, and the biomarker score was derived using the equation: the biomarker score=k×(β1×biomarker12×biomarker2+…+βk×biomarkerk)/(β12+…+βk), where k=4 in this study [32]. To construct the composite score, we first used each biomarker as an ordinal variable, which has been used in studies to create weighted biomarker score [32] and genetic risk scores [3637]. The ordinal values of each biomarker in quartile categories (coded as 0, 1, 2, and 3 from the lowest to the highest quartile, except for adiponectin which was assigned reversely as 3, 2, 1, and 0 due to its inverse association with T2DM risk) were used as continuous variables in the conditional logistic regression models, and the corresponding weight of each biomarker was assigned by its respective β coefficient from the regression models. We chose quartile categories because they yielded better model fit compared to tertile and quintile categories. Subsequently, the association between the composite biomarker score and T2DM risk was assessed using multivariable conditional logistic regression models with adjustment for the abovementioned covariates in the total population. To test for the influence of potential selection bias, we further assessed the association both among T2DM cases with HbA1c <6.0% and <6.5% at baseline along with their matched controls. In addition, potential interactions between the biomarker score and age (<60 or ≥60 years), alcohol intakes (never or monthly, weekly or daily), weekly activity levels (<0.5 or ≥0.5 hr/wk), BMI (<23 or ≥23 kg/m2), and fasting status (yes, no) were tested by adding an interaction term (interaction variable×biomarker score) to the regression models. For stratified analyses by variables other than matching factors, unconditional logistic regression models with further adjustment for sex and dialect group were used. In addition, since our study did not have a validation cohort and thus may subject to overfitting, we used 10-fold cross validation method in the sensitivity analysis to evaluate the goodness-of-fit of the model, and used the root mean squared error as the assessment metric [38].

To assess the predictive performance of the biomarker score, we established three logistic regression models. Base model 1 was established using a forward selection procedure (P<0.05), and included age, history of hypertension, smoking status, and BMI. Sex and dialect group were not included in the prediction model because cases and controls were perfectly matched on these factors. Base model 2 additionally included the levels of random glucose and random insulin. Since HbA1c has been shown to outperform glucose levels in T2DM prediction [39], we further substituted the random glucose levels with HbA1c levels and established model 3. We have tested the added value of each individual biomarker and the biomarker score to each base model for T2DM prediction, and the improvement in the discrimination between the base models and the model plus the biomarker/biomarker score was compared by area under receiver operating characteristic curve (AUC) [40]. Moreover, net reclassification improvement (NRI) and integrated discrimination improvement (IDI) statistics were further used as complementary methods [4142]. NRI calculates the overall differences between higher predicted risk for T2DM cases and lower predicted risk for controls, where a higher value indicates better predictive capacity. Furthermore, Akaike information criterion (AIC) was used to evaluate the model fitness, where a lower value indicated better model fit. In addition, we repeated the abovementioned analyses only among cases with HbA1c <6.0% or <6.5% at baseline and their matched controls.

In the sensitivity analysis, we created another score using each biomarker as a continuous variable (per log change) which captured more information compared to treating biomarkers as ordinal variables. We compared the performance of the two scores in terms of association with and prediction for T2DM risk. We used STATA version 14.0 (Stata Corp., College Station, TX, USA) for all analyses, and considered two-sided P values <0.05 as statistically significant.

RESULTS

Compared with controls, cases had worse metabolic profiles: they were more likely to be heavier and hypertensive, as well as having higher levels of blood biomarkers that were T2DM risk factors (ALT, TG, TG-to-HDL ratio, hs-CRP, ferritin, fetuin-A, RBP4, HbA1c) and lower levels of protective biomarkers (adiponectin, HDL-C). However, distributions of other baseline characteristics such as education levels, weekly activity levels, smoking status, alcohol consumption, and fasting status of blood samples were similar between cases and controls (Table 1).

Baseline characteristics and biomarker levels of diabetes cases and matched controls, the Singapore Chinese Health Study

The associations between all biomarkers and T2DM risks are presented in Table 2. In a multivariable model that included all the biomarkers and potential confounders, hs-CRP (OR with per quartile increment was 1.16; 95% confidence interval [CI], 0.98 to 1.38), fetuin-A (OR with per quartile increment was 1.10; 95% CI, 0.92 to 1.32) and RBP4 (OR with per quartile increment was 1.00; 95% CI, 0.83 to 1.21) were not significantly associated with T2DM risk and thus were not included in the composite biomarker score. In addition, TG-to-HDL ratio, ALT and ferritin were positively associated with T2DM risk, and the respective OR (95% CI) per quartile increment was 1.48 (95% CI, 1.21 to 1.82), 1.30 (95% CI, 1.08 to 1.57), and 1.24 (95% CI, 1.04 to 1.48), while adiponectin was inversely associated with T2DM risk (OR with per quartile increment was 0.72; 95% CI, 0.60 to 0.86).

Associations between per quartile increment of all the biomarkers and risk of type 2 diabetes mellitusa

We created a biomarker score based on the four biomarkers with significant associations, using their β coefficients associated with per quartile increments. The biomarker score ranged from 0 to 12 with a median value of 8.0. The median of the score was 9.1 (interquartile range [IQR], 7.2 to 10.9) among cases and 6.3 (IQR, 3.9 to 8.9) among controls (P<0.001) (Table 1). The distribution of the biomarker score among cases and controls is shown in Fig. 1, and the distributions of the biomarker score and included biomarkers according to the quartiles of the biomarker score are shown in Supplementary Table 1.

Fig. 1

Odds ratio for type 2 diabetes mellitus by the biomarker score and percentages of participants in each biomarker score category. The solid line represents the point estimates of relative risk for the association between the biomarker score and the risk of incident type 2 diabetes mellitus using conditional logistic regression model after adjustment for age at blood taken (years), smoking (never, ever smoker), alcohol intake (never, weekly, or daily), weekly activity (<0.5, 0.5 to 3, and ≥4 hr/wk), education level (primary school and below, secondary or above), history of hypertension (yes, no), body mass index (kg/m2), fasting status (yes, no), and levels of random glucose and random insulin, and the dotted lines represent the upper and lower bound of 95% confidence interval (CI). The light grey bars represent the percentage of controls within each category for the biomarker score (n=485), and the dark grey bars represent the percentage of cases within each category for the biomarker score (n=485).

The association between the biomarker score established based on the ordinal variables and T2DM risk is shown in Table 3 and Fig. 1. In the total population (485 case-control pairs), the biomarker score was positively associated with T2DM risk in a dose-dependent manner (P trend <0.001). Compared to those in the lowest quartile of the score, the OR was 12.0 (95% CI, 5.43 to 26.6; P trend <0.001) for those in the highest quartile. In addition, among 246 cases with baseline HbA1c <6.5% or 129 cases with HbA1c <6.0% and their respective matched controls, the strong positive association between the biomarker score and T2DM risk remained largely unchanged, and the OR comparing the highest versus lowest quartile of the biomarker score for T2DM risk was 8.62 (95% CI, 3.32 to 22.4; P trend <0.001) and 10.1 (95% CI, 2.47 to 41.4; P trend <0.001), respectively. Furthermore, in the stratified analysis, although the association was slightly stronger among participants with older age (≥60 years), higher BMI (≥23 kg/m2), more physical activity (≥0.5 hr/wk) and non-fasted samples compared to their respective counterparts, no significant interaction has been observed (Supplementary Table 2). In addition, the 10-fold cross validation test has suggested similar model fit from 10 attempts, and the root mean squared error ranged from 0.63 to 0.72.

Odds ratios (95% confidence intervals) of type 2 diabetes mellitus associated with quartile levels of the ordinal biomarker scorea

The predictive performance of the biomarker score is presented in Table 4. The AUC of the base model comprising age, smoking status, BMI and history of hypertension was 0.70 (95% CI, 0.66 to 0.73). When adding the composite biomarker score to the base model, the AUC was improved to 0.76 (95% CI, 0.73 to 0.79; P<0.001). In addition, when the base model was further included random glucose and insulin levels, adding the biomarker score increased AUC from 0.81 (95% CI, 0.78 to 0.83) to 0.83 (95% CI, 0.81 to 0.86; P=0.002). When substituting the random glucose levels and HbA1c in the base model, adding the biomarker score still yield a small yet statistically significant increment in AUC from 0.85 (95% CI, 0.83 to 0.88) to 0.86 (95% CI, 0.84 to 0.89; P=0.032). Moreover, when restricting to cases with baseline HbA1c <6.0% or <6.5% and their controls, significant improvements in AUCs were also observed. Furthermore, NRIs were all positive values, thus suggested that adding the biomarker score significantly improved the prediction of T2DM risk (all P<0.001).

Summary statistics to assess the biomarker score in predicting type 2 diabetes mellitus, the Singapore Chinese Health Study

To further stratify the NRI analyses by T2DM status, adding the biomarker score on top of base model 3 correctly assigned 32.0% of T2DM cases to a higher predicted T2DM risk and 15.5% of controls to a lower predicted T2DM risk in the total population (Supplementary Table 3). In the subgroup of cases with HbA1c <6.5% and matched controls, the biomarker score correctly re-assigned 31.7% cases and 14.6% controls. Moreover, in cases with HbA1c <6.0% and their matched controls, the biomarker score correctly reassigned 24.0% cases to higher T2DM risk and 8.6% controls to lower T2DM risk (Supplementary Table 3).

Since TG-to-HDL ratio and ALT are routine clinical tests, we also compared the performance the biomarker score with TG-to-HDL ratio and ALT alone. Adding to the base model comprising age, smoking status, BMI and history of hypertension, the biomarker score had stronger association than TG-to-HDL ratio and ALT alone (OR comparing the highest vs. lowest quartile: 13.3 [95% CI, 6.79 to 26.0] vs. 8.03 [95% CI, 4.45 to 14.5]), better discrimination (AUC, 0.76 [95% CI, 0.73 to 0.79] vs. 0.75 [95% CI, 0.72 to 0.78]), better model fit (AIC, 468 vs. 489), and additionally reclassified 12.4% participants to the correct T2DM risk category (Supplementary Table 4).

In the sensitivity analysis, we further created a biomarker score using each biomarker as a continuous variable, and calculated the corresponding β coefficient associated with per log increment of each biomarker (Supplementary Table 5). The continuous biomarker score ranged from −1.3 to 5.8 with a median value of 2.5, and the median of the score was 3.0 (IQR, 2.1 to 3.9) among cases and 1.8 (IQR, 0.6 to 2.9) in controls (P<0.001; data not shown). In terms of association, the continuous biomarker score has shown similar strength of association with T2DM risk compared to the ordinal score, and the OR comparing the extreme quartile of the score was 14.2 (95% CI, 6.03 to 33.3) in the total population (Supplementary Table 6). Regarding prediction, the continuous biomarker score had the same performance compared to the ordinal score: adding the continuous biomarker score to the base Model 1, 2, and 3 increased AUCs from 0.70 (95% CI, 0.66 to 0.73) to 0.76 (95% CI, 0.73 to 0.79), 0.81 (95% CI, 0.78 to 0.83) to 0.83 (95% CI, 0.81 to 0.86), and 0.85 (95% CI, 0.83 to 0.88) to 0.86 (95% CI, 0.84 to 0.89), respectively (all P<0.05) (Supplementary Table 7). In addition, the abovementioned association and prediction of the continuous score were essentially the same among cases with baseline HbA1c <6.0% or <6.5% and their matched controls.

DISCUSSION

In this middle-aged and older Chinese adults living in Singapore, we found that a weighted biomarker score, comprising TG-to-HDL ratio, ALT, ferritin, and adiponectin, was strongly associated with T2DM risk. The participants in the highest versus lowest quartile had an approximately 12-fold increased risk of T2DM. In addition, the biomarker score significantly improved the T2DM prediction, and it correctly reassigned 32.0% cases to higher T2DM risk and 15.5% controls to lower T2DM risk. Collectively, our results suggested that joint analysis of multiple biomarkers may be a useful tool to predict T2DM risk.

A few studies that have examined the predictive utility of multiple biomarkers for T2DM prediction have also included biomarkers that were common to the current study [1126272829303132]. A cross-sectional study in a Chinese population showed that a weighted biomarker score that included adiponectin, ferritin and two other inflammatory biomarkers (plasminogen activator inhibitor-1 and interleukin 6) substantially improved discrimination (AUC, NRI, and IDI) of T2DM risk [32]. Similarly, studies conducted across Western countries such as Denmark [11], Germany [2627], Switzerland [28], France [29], Finland [30], and United States [31] have also shown that adding adiponectin and/or ferritin along with other biomarkers have improved T2DM prediction on top of lifestyle characteristics, glucose and/or insulin parameters. Of note, the results from a prospective German cohort had an overlap of four biomarkers (TG, HDL-C, ALT, and adiponectin) with the current study, and the study showed that these biomarkers along with GGT and hs-CRP significantly improved T2DM on top of fasting glucose and HbA1c [26]. In addition, other risk scores developed in Asian populations (China [4344] and Korea [4546]) that mainly comprised of clinical information and lifestyle factors also included biomarkers such as lipids [43444546] or liver enzyme [46] that were similar to the current study.

The significant association between each biomarker included in the score (TG-to-HDL ratio, adiponectin, ferritin, and ALT) were consistent with findings across populations [172122232425]. Experimental studies have provided the biological plausibility for the underlying mechanisms of these biomarkers in the pathogenesis of T2DM. High TG or low HDL, resulting in a high TG-to-HDL ratio as a marker for dyslipidemia, may increase T2DM risk through lipotoxicity, inflammation, and endoplasmic reticulum stress [14]. Adiponectin, a well-known adipokine, has been shown to improve insulin sensitivity and glucose metabolism [47]. Ferritin, a widely-used marker for body iron, has been linked to T2DM risk through insulin resistance and β-cell function impairment [16]. ALT, the most specific marker for liver injury and a surrogate marker for nonalcoholic fatty liver disease [48], has been shown to impact T2DM risk by increasing hepatic insulin resistance, inflammation and dysregulation of lipid and glucose metabolism [49].

Although the clinical usefulness of the biomarker score has been shown in the current study, whether the biomarker score is also cost-effectiveness is unclear yet. Among the four biomarkers, TG-to-HDL ratio and ALT are routinely tested in the clinical setting, while measurements for adiponectin and ferritin are less common and could be expensive. Though evidence is scarce regarding the cost-effectiveness of utilization of adiponectin and ferritin in clinical settings, a study in the United States has shown that a biomarker score (PreDx Diabetes Risk Score) comprised of adiponectin, ferritin, and five other biomarkers (glucose, insulin, HbA1c, C-reactive protein, and interleukin 2 receptor-α) was more cost-effective than using the impaired fasting glucose alone over a 10-year period [50]. However, since no other study has examined this issue, future studies are warranted to evaluate the cost-effectiveness of the biomarker score developed in the current study.

To our best knowledge, this is the first prospective study in a Chinese population to evaluate the value of a weighted biomarker score for T2DM prediction. In addition, we comprehensively used AUC, NRI, and IDI to evaluate the predictive performance of the biomarker score. However, several limitations merit consideration. First of all, biomarkers were measured in one-time collection of blood samples from the participants and we could not study subsequent changes in these biomarkers. Furthermore, residual confounding may exist since body weight and height, and history of hypertension were self-reported. In addition, information on family history of T2DM and waist circumference was not collected; therefore, we could not apply the basic prediction model that were used in other studies which comprised of these variables [2830].

In conclusion, we have shown the utility of a composite score of four biomarkers as a predictor of T2DM in a Chinese population, and our finding suggest that this score is a promising marker as a screening tool to identify at-risk individuals for targeted diet and lifestyle intervention. Future studies are warranted to further validate these findings in other populations to determine cut-off values for more precise risk prediction of T2DM and to investigate the upstream factors that influence the levels of the biomarkers, as well as to evaluate the cost-effectiveness of the biomarker score in the local setting and other Asian countries.

ACKNOWLEDGMENTS

This work was supported by the National Medical Research Council, Singapore (NMRC/CIRG/1354/2013) and National Institutes of Health, USA (RO1 CA144034 and UM1 CA182876). Woon-Puay Koh is supported by the National Medical Research Council, Singapore (NMRC/CSA/0055/2013). An Pan is supported by the National Key Research and Development Program of China (2017YFC0907504).

We thank Siew-Hong Low of the National University of Singapore for supervising the fieldwork of the Singapore Chinese Health Study, and Renwei Wang for the maintenance of the cohort study database. We also thank the founding principal investigator of the Singapore Chinese Health Study, Mimi C. Yu as well as all the cohort participants.

Notes

CONFLICTS OF INTEREST: No potential conflict of interest relevant to this article was reported.

AUTHOR CONTRIBUTIONS:

  • Conception or design: W.P.K., X.S., A.P.

  • Acquisition, analysis, or interpretation of data: Y.W., W.P.K., X.S., J.M.Y., A.P.

  • Drafting the work or revising: Y.W., W.P.K., X.S., J.M.Y., A.P.

  • Final approval of the manuscript: Y.W., W.P.K., X.S., J.M.Y., A.P

References

1. Chan JC, Malik V, Jia W, Kadowaki T, Yajnik CS, Yoon KH, Hu FB. Diabetes in Asia: epidemiology, risk factors, and pathophysiology. JAMA 2009;301:2129–2140. 19470990.
2. International Diabetes Federation. IDF Diabetes Atlas 8th edth ed. Brussels: International Diabetes Federation; 2017.
3. Holman N, Young B, Gadsby R. Current prevalence of type 1 and type 2 diabetes in adults and children in the UK. Diabet Med 2015;32:1119–1120. 25962518.
4. Knowler WC, Barrett-Connor E, Fowler SE, Hamman RF, Lachin JM, Walker EA, Nathan DM. Diabetes Prevention Program Research Group. Reduction in the incidence of type 2 diabetes with lifestyle intervention or metformin. N Engl J Med 2002;346:393–403. 11832527.
5. Pan XR, Li GW, Hu YH, Wang JX, Yang WY, An ZX, Hu ZX, Lin J, Xiao JZ, Cao HB, Liu PA, Jiang XG, Jiang YY, Wang JP, Zheng H, Zhang H, Bennett PH, Howard BV. Effects of diet and exercise in preventing NIDDM in people with impaired glucose tolerance. The Da Qing IGT and Diabetes Study. Diabetes Care 1997;20:537–544. 9096977.
6. Meigs JB. Multiple biomarker prediction of type 2 diabetes. Diabetes Care 2009;32:1346–1348. 19564478.
7. Meigs JB, Wilson PW, Fox CS, Vasan RS, Nathan DM, Sullivan LM, D'Agostino RB. Body mass index, metabolic syndrome, and risk of type 2 diabetes or cardiovascular disease. J Clin Endocrinol Metab 2006;91:2906–2912. 16735483.
8. Meigs JB, Shrader P, Sullivan LM, McAteer JB, Fox CS, Dupuis J, Manning AK, Florez JC, Wilson PW, D'Agostino RB Sr, Cupples LA. Genotype score in addition to common risk factors for prediction of type 2 diabetes. N Engl J Med 2008;359:2208–2219. 19020323.
9. Talmud PJ, Cooper JA, Morris RW, Dudbridge F, Shah T, Engmann J, Dale C, White J, McLachlan S, Zabaneh D, Wong A, Ong KK, Gaunt T, Holmes MV, Lawlor DA, Richards M, Hardy R, Kuh D, Wareham N, Langenberg C, Ben-Shlomo Y, Wannamethee SG, Strachan MW, Kumari M, Whittaker JC, Drenos F, Kivimaki M, Hingorani AD, Price JF, Humphries SE. UCLEB Consortium. Sixty-five common genetic variants and prediction of type 2 diabetes. Diabetes 2015;64:1830–1840. 25475436.
10. Wong TY, Liew G, Tapp RJ, Schmidt MI, Wang JJ, Mitchell P, Klein R, Klein BE, Zimmet P, Shaw J. Relation between fasting glucose and retinopathy for diagnosis of diabetes: three population-based cross-sectional studies. Lancet 2008;371:736–743. 18313502.
11. Kolberg JA, Jorgensen T, Gerwien RW, Hamren S, McKenna MP, Moler E, Rowe MW, Urdea MS, Xu XM, Hansen T, Pedersen O, Borch-Johnsen K. Development of a type 2 diabetes risk model from a panel of serum biomarkers from the Inter99 cohort. Diabetes Care 2009;32:1207–1212. 19564473.
12. Unger RH. Minireview: weapons of lean body mass destruction: the role of ectopic lipids in the metabolic syndrome. Endocrinology 2003;144:5159–5165. 12960011.
13. Guilherme A, Virbasius JV, Puri V, Czech MP. Adipocyte dysfunctions linking obesity to insulin resistance and type 2 diabetes. Nat Rev Mol Cell Biol 2008;9:367–377. 18401346.
14. Li N, Fu J, Koonen DP, Kuivenhoven JA, Snieder H, Hofker MH. Are hypertriglyceridemia and low HDL causal factors in the development of insulin resistance? Atherosclerosis 2014;233:130–138. 24529133.
15. Cefalu WT. Inflammation, insulin resistance, and type 2 diabetes: back to the future? Diabetes 2009;58:307–308. 19171748.
16. Haap M, Machann J, von Friedeburg C, Schick F, Stefan N, Schwenzer NF, Fritsche A, Haring HU, Thamer C. Insulin sensitivity and liver fat: role of iron load. J Clin Endocrinol Metab 2011;96:E958–E961. 21430023.
17. Kunutsor SK, Apekey TA, Walley J. Liver aminotransferases and risk of incident type 2 diabetes: a systematic review and meta-analysis. Am J Epidemiol 2013;178:159–171. 23729682.
18. Guo VY, Cao B, Cai C, Cheng KK, Cheung BMY. Fetuin-A levels and risk of type 2 diabetes mellitus: a systematic review and meta-analysis. Acta Diabetol 2018;55:87–98. 29127490.
19. Luft VC, Pereira M, Pankow JS, Ballantyne C, Couper D, Heiss G, Duncan BB. ARIC Investigators. Retinol binding protein 4 and incident diabetes: the Atherosclerosis Risk in Communities Study (ARIC Study). Rev Bras Epidemiol 2013;16:388–397. 24142010.
20. Sun L, Qi Q, Zong G, Ye X, Li H, Liu X, Zheng H, Hu FB, Liu Y, Lin X. Elevated plasma retinol-binding protein 4 is associated with increased risk of type 2 diabetes in middle-aged and elderly Chinese adults. J Nutr 2014;144:722–728. 24647386.
21. Li S, Shin HJ, Ding EL, van Dam RM. Adiponectin levels and risk of type 2 diabetes: a systematic review and meta-analysis. JAMA 2009;302:179–188. 19584347.
22. Chien K, Cai T, Hsu H, Su T, Chang W, Chen M, Lee Y, Hu FB. A prediction model for type 2 diabetes risk among Chinese people. Diabetologia 2009;52:443–450. 19057891.
23. Wilson PW, Meigs JB, Sullivan L, Fox CS, Nathan DM, D'Agostino RB Sr. Prediction of incident diabetes mellitus in middle-aged adults: the Framingham Offspring Study. Arch Intern Med 2007;167:1068–1074. 17533210.
24. Wang X, Bao W, Liu J, Ouyang YY, Wang D, Rong S, Xiao X, Shan ZL, Zhang Y, Yao P, Liu LG. Inflammatory markers and risk of type 2 diabetes: a systematic review and meta-analysis. Diabetes Care 2013;36:166–175. 23264288.
25. Sun L, Zong G, Pan A, Ye X, Li H, Yu Z, Zhao Y, Zou S, Yu D, Jin Q, Hu FB, Lin X. Elevated plasma ferritin is associated with increased incidence of type 2 diabetes in middle-aged and elderly Chinese adults. J Nutr 2013;143:1459–1465. 23902953.
26. Schulze MB, Weikert C, Pischon T, Bergmann MM, Al-Hasani H, Schleicher E, Fritsche A, Haring HU, Boeing H, Joost HG. Use of multiple metabolic and genetic markers to improve the prediction of type 2 diabetes: the EPIC-Potsdam Study. Diabetes Care 2009;32:2116–2119. 19720844.
27. Herder C, Baumert J, Zierer A, Roden M, Meisinger C, Karakas M, Chambless L, Rathmann W, Peters A, Koenig W, Thorand B. Immunological and cardiometabolic risk factors in the prediction of type 2 diabetes and coronary events: MONICA/KORA Augsburg case-cohort study. PLoS One 2011;6e19852. 21674000.
28. Marques-Vidal P, Schmid R, Bochud M, Bastardot F, von Kanel R, Paccaud F, Glaus J, Preisig M, Waeber G, Vollenweider P. Adipocytokines, hepatic and inflammatory biomarkers and incidence of type 2 diabetes. the CoLaus study. PLoS One 2012;7e51768. 23251619.
29. Julia C, Czernichow S, Charnaux N, Ahluwalia N, Andreeva V, Touvier M, Galan P, Fezeu L. Relationships between adipokines, biomarkers of endothelial function and inflammation and risk of type 2 diabetes. Diabetes Res Clin Pract 2014;105:231–238. 24931702.
30. Salomaa V, Havulinna A, Saarela O, Zeller T, Jousilahti P, Jula A, Muenzel T, Aromaa A, Evans A, Kuulasmaa K, Blankenberg S. Thirty-one novel biomarkers as predictors for clinically incident diabetes. PLoS One 2010;5e10100. 20396381.
31. Raynor LA, Pankow JS, Duncan BB, Schmidt MI, Hoogeveen RC, Pereira MA, Young JH, Ballantyne CM. Novel risk factors and the prediction of type 2 diabetes in the Atherosclerosis Risk in Communities (ARIC) study. Diabetes Care 2013;36:70–76. 22933437.
32. Wu H, Yu Z, Qi Q, Li H, Sun Q, Lin X. Joint analysis of multiple biomarkers for identifying type 2 diabetes in middle-aged and older Chinese: a cross-sectional study. BMJ Open 2011;1e000191.
33. Chan JC, Yeung R, Luk A. The Asian diabetes phenotypes: challenges and opportunities. Diabetes Res Clin Pract 2014;105:135–139. 24947419.
34. Hankin JH, Stram DO, Arakawa K, Park S, Low SH, Lee HP, Yu MC. Singapore Chinese Health Study: development, validation, and calibration of the quantitative food frequency questionnaire. Nutr Cancer 2001;39:187–195. 11759279.
35. Odegaard AO, Pereira MA, Koh WP, Arakawa K, Lee HP, Yu MC. Coffee, tea, and incident type 2 diabetes: the Singapore Chinese Health Study. Am J Clin Nutr 2008;88:979–985. 18842784.
36. Lin X, Song K, Lim N, Yuan X, Johnson T, Abderrahmani A, Vollenweider P, Stirnadel H, Sundseth SS, Lai E, Burns DK, Middleton LT, Roses AD, Matthews PM, Waeber G, Cardon L, Waterworth DM, Mooser V. Risk prediction of prevalent diabetes in a Swiss population using a weighted genetic score: the CoLaus Study. Diabetologia 2009;52:600–608. 19139842.
37. De Jager PL, Chibnik LB, Cui J, Reischl J, Lehr S, Simon KC, Aubin C, Bauer D, Heubach JF, Sandbrink R, Tyblova M, Lelkova P, Havrdova E, Pohl C, Horakova D, Ascherio A, Hafler DA, Karlson EW. Steering committee of the BENEFIT study. Steering committee of the BEYOND study. Steering committee of the LTF study. Steering committee of the CCR1 study. Integration of genetic risk factors into a clinical algorithm for multiple sclerosis susceptibility: a weighted genetic risk score. Lancet Neurol 2009;8:1111–1119. 19879194.
38. Daniels B. CROSSFOLD: Stata module to perform k-fold cross-validation Boston: Boston College Department of Economics; 2012.
39. Selvin E, Steffes MW, Zhu H, Matsushita K, Wagenknecht L, Pankow J, Coresh J, Brancati FL. Glycated hemoglobin, diabetes, and cardiovascular risk in nondiabetic adults. N Engl J Med 2010;362:800–811. 20200384.
40. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics 1988;44:837–845. 3203132.
41. Pencina MJ, D'Agostino RB Sr, D'Agostino RB Jr, Vasan RS. Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med 2008;27:157–172. 17569110.
42. Pencina MJ, D'Agostino RB Sr, Steyerberg EW. Extensions of net reclassification improvement calculations to measure usefulness of new biomarkers. Stat Med 2011;30:11–21. 21204120.
43. Luo S, Han L, Zeng P, Chen F, Pan L, Wang S, Zhang T. A risk assessment model for type 2 diabetes in Chinese. PLoS One 2014;9e104046. 25101994.
44. Wang A, Chen G, Su Z, Liu X, Liu X, Li H, Luo Y, Tao L, Guo J, Liu L, Chen S, Wu S, Guo X. Risk scores for predicting incidence of type 2 diabetes in the Chinese population: the Kailuan prospective study. Sci Rep 2016;6:26548. 27221651.
45. Lim NK, Park SH, Choi SJ, Lee KS, Park HY. A risk score for predicting the incidence of type 2 diabetes in a middle-aged Korean cohort: the Korean genome and epidemiology study. Circ J 2012;76:1904–1910. 22640983.
46. Ha KH, Lee YH, Song SO, Lee JW, Kim DW, Cho KH, Kim DJ. Development and validation of the Korean diabetes risk score: a 10-year national cohort study. Diabetes Metab J 2018;42:402–414. 30113144.
47. Kwon H, Pessin JE. Adipokines mediate inflammation and insulin resistance. Front Endocrinol (Lausanne) 2013;4:71. 23781214.
48. Ballestri S, Zona S, Targher G, Romagnoli D, Baldelli E, Nascimbeni F, Roverato A, Guaraldi G, Lonardo A. Nonalcoholic fatty liver disease is associated with an almost twofold increased risk of incident type 2 diabetes and metabolic syndrome. Evidence from a systematic review and meta-analysis. J Gastroenterol Hepatol 2016;31:936–944. 26667191.
49. Vozarova B, Stefan N, Lindsay RS, Saremi A, Pratley RE, Bogardus C, Tataranni PA. High alanine aminotransferase is associated with decreased hepatic insulin sensitivity and predicts the development of type 2 diabetes. Diabetes 2002;51:1889–1895. 12031978.
50. Sullivan SD, Garrison LP Jr, Rinde H, Kolberg J, Moler EJ. Cost-effectiveness of risk stratification for preventing type 2 diabetes using a multi-marker diabetes risk score. J Med Econ 2011;14:609–616. 21740291.

SUPPLEMENTARY MATERIALS

Supplementary materials related to this article can be found online at https://doi.org/10.4093/dmj.2019.0020.

Supplementary Table 1

Distribution of the biomarker score and included biomarkers according to the quartiles of the biomarker score

dmj-44-295-s001.pdf

Supplementary Table 2

Odds ratios (95% confidence intervals) of type 2 diabetes mellitus by stratified analysis

dmj-44-295-s002.pdf

Supplementary Table 3

Reclassification of diabetes cases and controls with no categories based on their ordinal biomarker scorea

dmj-44-295-s003.pdf

Supplementary Table 4

Summary statistics to assess the biomarker score versus TG-to-HDL ratio and ALT in predicting type 2 diabetes mellitus, the Singapore Chinese Health Study

dmj-44-295-s004.pdf

Supplementary Table 5

Associations between per log increment of all the biomarkers and risk of type 2 diabetes mellitusa

dmj-44-295-s005.pdf

Supplementary Table 6

Odds ratios (95% confidence intervals) of type 2 diabetes mellitus associated with quartile levels of the continuous biomarker scorea

dmj-44-295-s006.pdf

Supplementary Table 7

Summary statistics to assess the continuous biomarker score in predicting type 2 diabetes mellitus, the Singapore Chinese Health Study

dmj-44-295-s007.pdf

Supplementary Fig. 1

Flowchart of the Singapore Chinese Health Study. HbA1c, glycosylated hemoglobin.

dmj-44-295-s008.pdf

Article information Continued

Funded by : National Medical Research Councilhttps://doi.org/10.13039/501100001349
Award ID : NMRC/CIRG/1354/2013
Award ID : NMRC/CSA/0055/2013
Funded by : National Institutes of Healthhttps://doi.org/10.13039/100000002
Award ID : RO1 CA144034
Award ID : UM1 CA182876
Funded by : National Key Research and Development Program of China
Award ID : 2017YFC0907504

Fig. 1

Odds ratio for type 2 diabetes mellitus by the biomarker score and percentages of participants in each biomarker score category. The solid line represents the point estimates of relative risk for the association between the biomarker score and the risk of incident type 2 diabetes mellitus using conditional logistic regression model after adjustment for age at blood taken (years), smoking (never, ever smoker), alcohol intake (never, weekly, or daily), weekly activity (<0.5, 0.5 to 3, and ≥4 hr/wk), education level (primary school and below, secondary or above), history of hypertension (yes, no), body mass index (kg/m2), fasting status (yes, no), and levels of random glucose and random insulin, and the dotted lines represent the upper and lower bound of 95% confidence interval (CI). The light grey bars represent the percentage of controls within each category for the biomarker score (n=485), and the dark grey bars represent the percentage of cases within each category for the biomarker score (n=485).

Table 1

Baseline characteristics and biomarker levels of diabetes cases and matched controls, the Singapore Chinese Health Study

Characteristic Cases Controls P valuea
Age at blood taken, yr 59.4±5.94 59.4±6.05 -
Female sex 273 (56.3) 273 (56.3) -
 Dialect -
 Cantonese 242 (49.9) 242 (49.9)
 Hokkien 243 (50.1) 243 (50.1)
Body mass index, kg/m2 24.9±3.62 22.8±3.26 <0.001
 Level of education 0.423
 No formal education 77 (15.9) 73 (15.1)
 Primary school 220 (45.4) 204 (42.1)
 Secondary and above 188 (38.8) 208 (42.9)
History of hypertension 229 (47.2) 121 (25.0) <0.001
Cigarette smoking 0.107
 Never smokers 344 (70.9) 360 (74.2)
 Former smoker 58 (12.0) 65 (13.4)
 Current smokers 83 (17.1) 60 (12.4)
Weekly moderate-to-vigorous activity, hr/wk 0.058
 <0.5 387 (79.8) 384 (79.2)
 0.5–3.9 72 (14.9) 58 (12.0)
 ≥4.0 26 (5.4) 42 (8.9)
Alcohol intake 0.950
 Abstainers 424 (87.4) 421 (86.8)
 Weekly drinkers 47 (9.7) 50 (10.3)
 Daily drinkers 14 (2.9) 14 (2.9)
Fasting status (yes) 157 (32.4) 145 (29.9) 0.405
TG, mmol/L 2.12 (1.45–2.85) 1.51 (1.05–2.14) <0.001
HDL-C, mmol/L 1.09 ± 0.24 1.24±0.32 <0.001
TG-to-HDL ratio 1.94 (1.27–2.93) 1.23 (0.78–2.01) <0.001
ALT, IU/L 27 (21–37) 20 (15–27) <0.001
hs-CRP, mg/L 1.8 (1.0–3.5) 1.2 (0.6–2.2) <0.001
Ferritin, µg/L 185 (106–283) 131 (77–201) <0.001
Adiponectin, µg/mL 6.8 (5.2–8.4) 8.4 (6.5–10.7) <0.001
Fetuin-A, µg/mL 730 (564–931) 650 (506–861) <0.001
RBP4, µg/mL 28 (23–34) 26 (23–32) <0.001
HbA1c, % 6.4 (5.9–7.2) 5.6 (5.4–5.7) <0.001
HbA1c, mmol/mol 46 (41–55) 38 (36–39) <0.001
Biomarker score 9.1 (7.2–10.9) 6.3 (3.9–8.9) <0.001

Values are presented as mean±standard deviation, number (%), or median (interquartile range). Cases and controls are matched on age at blood taken (±3 years), gender, dialect, and date of blood collection (±6 months).

TG, triglycerides; HDL-C, high density lipoprotein cholesterol; ALT, alanine aminotransferase; hs-CRP, high-sensitivity C-reactive protein; RBP4, retinol-binding protein 4; HbA1c, glycosylated hemoglobin.

aP values were based conditional logistic regression models.

Table 2

Associations between per quartile increment of all the biomarkers and risk of type 2 diabetes mellitusa

Variable OR (95% CI) per quartile increment β Coefficient from model 2c
Model 1b Model 2c
TG-to-HDL ratio 1.90 (1.58–2.28) 1.48 (1.21–1.82) 0.39
ALT 1.68 (1.42–1.98) 1.30 (1.08–1.57) 0.26
hs-CRP 1.37 (1.17–1.59) 1.16 (0.98–1.38) 0.15
Ferritin 1.40 (1.20–1.63) 1.24 (1.04–1.48) 0.21
Adiponectin 0.58 (0.49–0.68) 0.72 (0.60–0.86) −0.33
Fetuin-A 1.27 (1.08–1.49) 1.10 (0.92–1.32) 0.10
RBP4 1.12 (0.96–1.32) 1.00 (0.83–1.21) 0.01

OR, odds ratio; CI, confidence interval; TG-to-HDL ratio, the ratio of triglycerides to high density lipoprotein cholesterol; ALT, alanine aminotransferase; hs-CRP, high-sensitivity C-reactive protein; RBP4, retinol-binding protein 4.

aThe sample size for all the biomarkers was 485 type 2 diabetes mellitus cases and 485 controls. Cases and controls were matched on age at blood taken (±3 years), sex, dialect, and date of blood collection (±6 months), bModel 1 was calculated using conditional logistic regression model with adjustment for age at blood taken (continuous), smoking (never, former, and current smoker), alcohol intake (never, weekly, or daily), weekly activity (<0.5, 0.5 to 3.9, and ≥4.0 hr/wk), education level (primary school and below, secondary or above), history of hypertension (yes, no), body mass index (continuous), and fasting status (yes, no), cModel 2: model 1 plus adjustment for all the other biomarkers (per quartile increment).

Table 3

Odds ratios (95% confidence intervals) of type 2 diabetes mellitus associated with quartile levels of the ordinal biomarker scorea

Variable Biomarker score P trendb
Q1 Q2 Q3 Q4
Total dataset
 Median (range) 2.29 (0–3.87) 5.24 (3.93–6.29) 7.73 (6.32–8.87) 10.5 (8.91–12.0)
 No. of cases/controls 28/123 58/123 132/118 267/121
 Model 1c 1.00 2.63 (1.36–5.09) 7.24 (3.76–13.9) 13.3 (6.79–26.0) <0.001
 Model 2d 1.00 2.71 (1.24–5.92) 6.40 (2.94–13.9) 12.0 (5.43–26.6) <0.001
Limited to cases with baseline HbA1c <6.5% and matched controls
 Median (range) 2.28 (0–3.87) 5.14 (3.93–6.29) 7.60 (6.32–8.87) 10.6 (8.91–12.0)
 No. of cases/controls 20/65 34/65 71/63 121/53
 Model 1c 1.00 2.46 (1.05–5.76) 4.90 (2.15–11.2) 9.68 (4.03–23.3) <0.001
 Model 2d 1.00 2.88 (1.14–7.28) 4.14 (1.60–10.7) 8.62 (3.32–22.4) <0.001
Limited to cases with baseline HbA1c <6.0% and matched controls
 Median (range) 2.12 (0–3.87) 5.21 (3.93–6.29) 8.00 (6.45–8.87) 10.7 (8.91–12.0)
 No. of cases/controls 14/34 25/30 35/35 55/30
 Model 1c 1.00 3.17 (1.07–9.39) 4.89 (1.56–15.4) 8.25 (249–27.4) 0.001
 Model 2d 1.00 3.71 (1.12–12.4) 4.94 (1.18–20.7) 10.1 (2.47–41.4) 0.002

HbA1c, glycosylated hemoglobin.

aThe biomarker score was constructed using each biomarker (triglycerides to high density lipoprotein cholesterol ratio, alanine aminotransferase, ferritin, and adiponectin) as ordinal variables, bLinear trend was tested by using the median level of each quartile of the biomarker score, cModel 1 was calculated using conditional logistic regression models after adjusting for age at blood taken (continuous), smoking (never, former, and current smoker), alcohol intake (never, weekly, or daily), weekly activity (<0.5, 0.5 to 3.9, and ≥4.0 hr/wk), education level (primary school and below, secondary or above), history of hypertension (yes, no), body mass index (continuous), and fasting status (yes, no), dModel 2: model 1 plus levels of random glucose and random insulin (both in quartiles).

Table 4

Summary statistics to assess the biomarker score in predicting type 2 diabetes mellitus, the Singapore Chinese Health Study

Variable Multivariable model
Discrimination AUC (95% CI) Calibration (AIC) NRI IDI
Total dataset
 Base model 1a 0.70 (0.66–0.73) 558
 Base model 1a+biomarker scoreb 0.76 (0.73–0.79)c 501 0.56 0.08
 Base model 2d 0.81 (0.78–0.83) 418
 Base model 2d+biomarker scoreb 0.83 (0.81–0.86)c 369 0.59 0.06
 Base model 3e 0.85 (0.83–0.88) 338
 Base model 3e+biomarker scoreb 0.86 (0.84–0.89)c 303 0.47 0.04
Limited to cases with baseline HbA1c <6.5% and matched controls
 Base model 1a 0.70 (0.66–0.75) 280
 Base model 1a+biomarker scoreb 0.75 (0.70–0.79)c 252 0.54 0.06
 Base model 2d 0.75 (0.71–0.80) 261
 Base model 2d+biomarker scoreb 0.78 (0.74–0.82)c 242 0.50 0.05
 Base model 3e 0.81 (0.78–0.85) 215
 Base model 3e+biomarker scoreb 0.83 (0.79–0.87)f 204 0.46 0.04
Limited to cases with baseline HbA1c <6.0% and matched controls
 Base model 1a 0.65 (0.59–0.72) 177
 Base model 1a+biomarker scoreb 0.68 (0.62–0.75)c 172 0.36 0.03
 Base model 2d 0.71 (0.65–0.77) 194
 Base model 2d+biomarker scoreb 0.73 (0.67–0.79) 189 0.37 0.03
 Base model 3e 0.71 (0.65–0.78) 194
 Base model 3e+biomarker scoreb 0.74 (0.68–0.80) 186 0.33 0.03

AUC, area under the receiver operating characteristic curve; CI, confidence interval; AIC, Akaike information criterion; NRI, net reclassification improvement; IDI, integrated discrimination improvement; HbA1c, glycosylated hemoglobin.

aBase model 1 included age at blood taken (continuous), smoking (never, former, and current smoker), history of hypertension (yes, no), and body mass index (continuous), bThe biomarker score was constructed using each biomarker (triglycerides to high density lipoprotein cholesterol ratio, alanine aminotransferase, ferritin, and adiponectin) as ordinal variables, and was used as categorical variables (in quartiles) in the prediction model, cCompared with the base model, the increment in AUC value was statistically significant (P<0.05), dBase model 2: base model 1 plus random levels of glucose and insulin (both in quartiles), eBase model 3: base model 1 plus levels of HbA1c and random insulin (both in quartiles), fCompared with the base model, the increment in AUC value was marginally significant (P=0.052).