Specific triacylglycerol, diacylglycerol, and lyso-phosphatidylcholine species for the prediction of type 2 diabetes: a ~ 16-year prospective study in Chinese
Cardiovascular Diabetology volume 21, Article number: 234 (2022)
Bioactive lipids play an important role in insulin secretion and sensitivity, contributing to the pathophysiology of type 2 diabetes (T2D). This study aimed to identify novel lipid species associated with incident T2D in a nested case–control study within a long-term prospective Chinese community-based cohort with a median follow-up of ~ 16 years.
Plasma samples from 196 incident T2D cases and 196 age- and sex-matched non-T2D controls recruited from the Hong Kong Cardiovascular Risk Factor Prevalence Study (CRISPS) were first analyzed using untargeted lipidomics. Potential predictive lipid species selected by the Boruta analysis were then verified by targeted lipidomics. The associations between these lipid species and incident T2D were assessed. Effects of novel lipid species on insulin secretion in mouse islets were investigated.
Boruta analysis identified 16 potential lipid species. After adjustment for body mass index (BMI), triacylglycerol/high-density lipoprotein (TG/HDL) ratio and the presence of prediabetes, triacylglycerol (TG) 12:0_18:2_22:6, TG 16:0_11:1_18:2, TG 49:0, TG 51:1 and diacylglycerol (DG) 18:2_22:6 were independently associated with increased T2D risk, whereas lyso-phosphatidylcholine (LPC) O-16:0, LPC P-16:0, LPC O-18:0 and LPC 18:1 were independently associated with decreased T2D risk. Addition of the identified lipid species to the clinical prediction model, comprised of BMI, TG/HDL ratio and the presence of prediabetes, achieved a 3.8% improvement in the area under the receiver operating characteristics curve (AUROC) (p = 0.0026). Further functional study revealed that, LPC O-16:0 and LPC O-18:0 significantly potentiated glucose induced insulin secretion (GSIS) in a dose-dependent manner, whereas neither DG 18:2_22:6 nor TG 12:0_18:2_22:6 had any effect on GSIS.
Addition of the lipid species substantially improved the prediction of T2D beyond the model based on clinical risk factors. Decreased levels of LPC O-16:0 and LPC O-18:0 may contribute to the development of T2D via reduced insulin secretion.
The global epidemic of type 2 diabetes (T2D) has been growing faster than expected over the last two decades, rising from over 150 million cases in 2000 to more than 530 million cases in 2021 . This number was projected to further increase to more than 640 million in 2030 and over 780 million in 2045 . To reduce the global burden of diabetes, it is important to discover novel biomarkers both for the early prediction of T2D onset and a comprehensive understanding of the pathophysiological changes preceding the disease.
Bioactive lipids play an important role in the regulation of glucose homeostasis [2, 3]. It is well established that the dysregulation of these lipids leads to the development of T2D. The overflow of lipids leads to lipotoxicity, causing impaired insulin signaling in skeletal muscle, altered secretion of hepatokines in the liver , and pancreatic β-cell dysfunction or cell loss . On the other hand, certain lipids have shown a beneficial effect on insulin secretion . Traditional indicators of lipid metabolism measured in the blood, such as triacylglycerols (TGs), high-density lipoprotein (HDL), low-density lipoprotein (LDL), TG/HDL ratio and triglyceride-glucose (TyG) index, have been identified as predictors of T2D in previous studies . However, the plasma lipidome comprises thousands of lipid species. Clinical assays of blood lipids failed to address the complexity and diversity of lipid species. With the recent advances in lipidomics, the identification resolution and coverage of lipid species have been unprecedentedly enhanced and expanded.
Accumulative studies have investigated the association between human plasma lipidome and the pathogenesis of T2D [8,9,10,11,12,13,14]. Lipid species identified by lipidomics have been shown to be superior over the traditional lipotoxicity indicators in the prediction of T2D development. Prospective lipidomic studies in Asians assessing the link between lipid species and the risk of T2D have recently emerged [8, 9, 15,16,17,18]. However, most of them have relatively short follow-up duration and have adopted the targeted approach which has limited ability to provide a comprehensive investigation of novel predictive lipid species. Here we used a cohort of Chinese residents in Hong Kong (mean age ± standard deviation: 51.13 ± 11.57 years) with a long follow-up period of almost 16 years to identify novel lipid species that are associated with incident T2D.
In the current study, we aimed to: (i) discover novel lipid species associated with incident T2D by performing an untargeted lipidomic analysis in a nested case–control study within a prospective population-based cohort; (ii) evaluate the improvement of prediction for incident T2D provided by the identified lipid species; and (iii) investigate the effect of the identified lipid species on insulin secretion.
A nested case–control study was performed in the Hong Kong Cardiovascular Risk Factor Prevalence Study (CRISPS). CRISPS is a prospective, population-based longitudinal cohort designed to study the development of cardiovascular risk factors, including diabetes, in Hong Kong from 1995 to 2018. Details of the CRISPS cohort were previously described elsewhere [19, 20]. Briefly, in 1995–1996 (CRISPS-1), 2895 Hong Kong Chinese, aged 25–74, were randomly recruited from the general population of Hong Kong to undergo a comprehensive assessment. Subjects were followed up in CRISPS-2 (2000–2004), CRISPS-3 (2005–2008), CRISPS-4 (2010–2012) and CRISPS-5 (2016–2018). All subjects had blood taken at each visit after an overnight fast of at least 10 h. Details of anthropometric measurements and methods for the measurement of biochemical parameters were described previously . Hypertension was defined as blood pressure ≥ 140/90 mmHg or the use of antihypertensive medications. A 75 g oral glucose tolerance test (OGTT) was performed in all except those taking antidiabetic medications. T2D was defined as fasting plasma glucose ≥ 7.0 mmol/l or 2-h plasma glucose ≥ 11.1 mmol/l during OGTT, or on antidiabetic medications, according to the World Health Organization (WHO) 1998 diagnostic criteria . Impaired fasting glucose (IFG) was defined as fasting plasma glucose 6.1–6.9 mmol/L and 2-h plasma glucose < 7.8 mmol/L. Impaired glucose tolerance (IGT) was defined as fasting plasma glucose < 7 mmol/L and 2-h plasma glucose ≥ 7.8 and < 11.1 mmol/L . Prediabetes was defined as the presence of IFG and/or IGT. TG/HDL ratio was calculated by the following formula: TG/HDL ratio = TG (mg/dl)/HDL (mg/dl). Homeostatic model assessment of insulin resistance (HOMA-IR) was calculated using the following formula: HOMA-IR = fasting insulin (µU/l) × fasting glucose (mmol/l)/22.5 . Homeostatic model assessment of β-cell function (HOMA-β) was calculated using the following formula: HOMA-β = [20 × fasting insulin (µU/ml)]/[fasting glucose (mmol/l) − 3.5] . All participants had given written informed consent. Ethics approval was obtained from the Institutional Review Board of the University of Hong Kong/Hospital Authority, Hong Kong West Cluster. The lipidomic analysis was performed on the stored plasma samples collected at CRISPS-2, therefore, this visit was considered as the baseline of the current study.
Subjects who were free from T2D at CRISPS-2 (i.e., baseline of the current study) were followed for their glycaemic status at the subsequent visits. Incident T2D cases included those who had developed T2D when assessed at the CRISPS-3, CRISPS-4 or CRISPS-5 visits. The persistent non-T2D controls were those who remained free from T2D at CRISPS-5. In this prospective nested case–control study, the eligible baseline population was limited to those who did not have T2D at baseline and with plasma samples available for the lipidomic analysis. Subjects who were on lipid-lowering drugs at baseline were further excluded to minimize the possible drug effects on lipid metabolism. Each of the 196 incident T2D cases who fulfilled the inclusion criteria was age- and sex-matched with a non-T2D control who remained free of T2D at the end of the follow-up (CRISPS-5) using the propensity score matching method.
All lipid standards were purchased from Avanti Polar Lipids except for triacylglycerol (TG) 15:0_15:0_15:0, lyso-phosphatidylcholine (LPC) O-16:0, LPC O-18:0 and palmitic acid. TG 15:0_15:0_15:0 and palmitic acid were obtained from Sigma-Aldrich. DG 18:2_22:6 and TG 12:0_18:2_22:6 were synthesized using a method described by Halldorsson et al. , as detailed in Additional file 1: Methods. LPC O-16:0 and LPC O-18:0 were purchased from Cayman Chemical (Michigan, USA) and diluted in phosphate-buffered saline (PBS) at a stock concentration of 5 mmol/l prior to the analysis. Palmitic acid was dissolved in ethanol at a stock concentration of 20 mmol/l. DG 18:2_22:6 and TG 12:0_18:2_22:6 were dissolved in dimethyl sulfoxide at a stock concentration of 20 mmol/l.
Lipids were extracted from 40 μL of plasma based on a modified method of Matyash et al. . Untargeted lipidomic analysis was performed on a Thermo UltiMate 3000 ultra-high performance liquid chromatography tandem mass spectrometry (UHPLC) system coupled to a Thermo Orbitrap Fusion mass spectrometer . Targeted lipidomic analysis of the selected lipids was performed on a Vanquish™ UHPLC Systems coupled to TSQ Altis™ Triple Quadrupole Mass Spectrometer (Thermo Scientific, USA). More details are provided in Additional file 1: Methods.
Mouse islet isolation
Mouse pancreatic islets were isolated, cultured overnight and picked under a microscope for the glucose-stimulated insulin secretion assay (GSIS). The detail of islet isolation is described in Additional file 1: Methods. All animal experimental protocols were approved by the Animal Ethics Committee of The University of Hong Kong.
Glucose-stimulated insulin secretion assay
The isolated mouse islets were fasted for 120 min with glucose-free Krebs buffer containing 10 mmol/l HEPES, 129 mmol/l NaCl, 4.8 mmol/l KCl, 1.2 mmol/l MgSO4∙7H2O, 1.2 mmol/l KH2PO4, 2.5 mmol/l CaCl2∙2H2O, 5 mmol/l NaHCO3 and 0.1% fatty acid-free bovine serum albumin (pH 7.4). Next, the islets were incubated with Krebs buffer containing 2 mmol/l glucose, 20 mmol/l glucose, or 20 mmol/l glucose with 10 μmol/l or 50 μmol/l of LPC O-16:0, LPC O-18:0, DG 18:2_22:6, TG 12:0_18:2_22:6, and palmitic acid for another 30 min (n = 9 in each experimental group with 10 islets per well). The buffer samples were then collected for the measurement of insulin level using an in-house mouse high-sensitive insulin ELISA kit (ImmunoDiagnostics Limited, The University of Hong Kong).
All analyses were conducted in R v4.3.0. Normally distributed data were presented as means ± standard deviation. Non-normally distributed data (determined by Kolmogorov–Smirnov test) were transformed by natural logarithm to near normality before the analysis and presented as median (interquartile). Missing data were rare (Table 1) and median imputation was used for the missing values. Student t-tests were used to compare continuous variables while Pearson χ2 tests and Fisher’s exact tests were applied to compare categorical variables at baseline examination between the incident case and non-T2D control groups.
Lipid species levels were transformed by natural logarithm to near normality before the analysis. To investigate the connection within the lipid species, the R package “weighted gene co-expression network analysis (WGCNA)” was used to determine modules of highly interconnected lipid species. It includes feature co-expression, network construction, module identification, module-phenotype correlation recognition, and key driver gene identification. Highly connected lipid species were defined in a colored module by topological overlap measure (TOM). Finally, a weighted co-expression network generated by WGCNA was visualized by Cytoscape (version 3.7.1).
Boruta analysis was applied to select features that were most important for the prediction of T2D. The algorithm generates randomly shuffled copies of lipid species concentrations (shadow features). Lipid species with higher Z-scores than the maximum Z-score of their shadow features were categorized either as confirmed or tentative features . These lipid species were then selected and quantified using the targeted lipidomic analysis.
Stepwise model selection approach was used to identify the best model that comprised the most important clinical variables for incident T2D. Multiple conditional logistic regression analyses were then performed to assess the independent associations of the lipid species with incident T2D after adjustment for the most important clinical predictors identified. The false discovery rate (FDR) using the Benjamin-Hochberg method was employed for the correction of multiple testing. A two-sided p-value < 0.05 and an FDR-adjusted p-value (q-value) < 0.1 was considered statistically significant. Area under the receiver operating curves (AUROCs) was estimated to evaluate the predictive ability of the lipid species that showed statistically significant and independent associations with incident T2D. The Delong test was used for comparing the AUROCs. The improvement of the predictive ability given by the lipid species was further quantified by continuous net reclassification index (cNRI) and integrated discrimination index (IDI) . A two-tailed p-value < 0.05 was considered statistically significant.
For the GSIS study, data were presented as the mean ± mean standard error (SEM) and were analyzed using student’s t-test or one-way analysis of variance with the post-hoc test. A two-sided p-value < 0.05 was considered statistically significant.
The baseline characteristics of the study participants are presented in Table 1. The median follow-up time of the 196 incident T2D cases and 196 age- and sex-matched non-T2D controls was 15.57 years (interquartile range [IQR] 9.98–16.97). As expected, subjects who had developed T2D showed significantly higher BMI, fasting glucose, 2-h glucose, the presence of prediabetes, HOMA-IR, TG, and TG/HDL ratio at baseline compared to the non-T2D controls. Based on stepwise model selection, the best clinical model predictive of incident T2D comprised of BMI, TG/HDL ratio and the presence of prediabetes (Additional file 2: Table S1).
Lipid profiling and co-expression network
A total of 301 lipid species across 13 classes were detected from the untargeted lipidomic analysis (Additional file 2: Table S2). To understand the complex processes of lipid dysregulation preceding T2D, a lipid co-expression network was constructed using WGCNA and 13 modules were generated (Additional file 1: Fig. S1a–c and Additional file 2: Table S3). The lipid species that were not co-expressed with other lipid species were grouped in the grey module, which was ignored in further analyses. The threshold analyses and the cluster dendrogram of the network construction were shown in Additional file 1: Fig. S1a–c. After adjustment for BMI, TG/HDL ratio and the presence of prediabetes, the eigenvalues of the blue, pink, green-yellow, black, and turquoise modules containing lipid classes of TG, diacylglycerol (DG), phosphatidylcholine (PC) and phosphatidylethanolamine (PE), were positively correlated with incident T2D (p < 0.05), whereas the red, green, yellow and brown modules comprising of the ether-PCs, ether-PEs, sphingomyelins (SMs), acyl-LPCs and ether-LPCs were negatively associated with T2D risk (Additional file 1: Fig. S1d, e). An interaction network generated by Cytoscape 3.8.2 based on these modules was shown in Additional file 1: Fig. S1f. Lipid species with higher connectivity were closer to the center of the network. In the multiple conditional logistic regression analysis based on the quartiles of module eigenvalues, the highest quartile of module blue (mainly TGs and PCs) was positively associated with T2D compared to the lowest quartile (OR [95% CI] 2.18 [1.13–4.22]; p for trend = 0.001) after adjustments for BMI, TG/HDL ratio and the presence of prediabetes. On the other hand, module brown (mainly acyl-LPCs and ether-LPCs) was inversely associated with T2D (lowest OR [95% CI] 0.21 [0.09–0.47]; p for trend = 0.035) (Table 2).
Lipid species selection
Boruta analysis was used to select the most important lipid species for the prediction of T2D. The levels of importance of the lipid species were compared with the shadow-max value and 16 confirmed lipid species were identified as potential predictors of incident T2D (Additional file 2: Table S4 and Additional file 1: Fig. S2). Among these, 11 lipid species were also highlighted in the interaction network analysis (Additional file 1: Fig. S1f). These 16 selected lipid species, including 4 PCs, 5 TGs, 2 DGs, 4 LPCs and one PE, were then carried forward for absolute quantification using the targeted lipidomic analysis.
Lipid species and the risk of T2D
All of the selected lipid species quantified by the targeted lipidomic analysis, except DG 18:0_18:0, were significantly different between the incident T2D cases and non-T2D controls (all unadjusted p < 0.01; q-value < 0.10) (Additional file 1: Fig. S3). Lipid species of the same category were highly correlated and the LPCs were inversely correlated with DGs as well as TGs (Additional file 1: Fig. S4).
Multiple conditional logistic regression analysis was adopted to identify the independent lipid predictors for incident T2D (Table 3). After adjustment for BMI, TG/HDL ratio, and the presence of prediabetes, TG 12:0_18:2_22:6, TG 16:0_11:1_18:2, TG 49:0, TG 51:1 and DG 18:2_22:6 were shown to be positively associated with incident T2D, and LPC O-16:0, LPC P-16:0, LPC O-18:0 and LPC 18:1 were inversely associated with incident T2D (p < 0.05, q-value < 0.10) (Table 3). Recognizing the potential collinearity between TG and TG/HDL ratio, we replaced TG/HDL ratio with TG in the analysis and the result remained similar.
Next, we assessed whether the independent lipid species showed an incremental prediction value on the development of T2D using ROC curves analyses (Fig. 1). Adding TG 12:0_18:2_22:6, TG 16:0_11:1_18:2, TG 16:0_16:0_17:0, TG 16:0_17:0_18:1, DG 18:2_22:6, LPC O-16:0, LPC P-16:0, LPC O-18:0 and LPC 18:1 to the clinical risk model significantly increased the AUROC from 0.785 to 0.823 (improved by 3.8%, p = 0.0026). The improvement was further confirmed by both cNRI (37.8%, [95% CI 18.4–57.1%], p < 0.001) and IDI (3.7% [95% CI 1.9–5.6%], p < 0.001).
LPC O-16:0 and LPC O-18:0, but not DG 18:2_22:6 and TG 12:0_18:2_22:6, potentiated glucose-stimulated insulin secretion in mouse islets
The majority of East Asian patients with T2D, including Chinese, show prominent defects in insulin secretion relative to insulin resistance, and β-cell dysfunction plays a key role in the development of T2D among East Asian populations [28, 29]. We observed an inverse association of several LPCs but positive correlation of some glycerolipids with the risk of incident T2D. To interrogate whether the newly identified LPC O-16:0, LPC O-18:0, DG 18:2_22:6 and TG 12:0_18:2_22:6 modulate insulin secretion, we further investigated their possible effects on GSIS using mouse islets. Isolated mouse islets were incubated with 20 mmol/l glucose only or further supplemented with 10 µmol/l or 50 µmol/l of LPC O-16:0, LPC O-18:0, DG 18:2_22:6, TG 12:0_18:2_22:6, or palmitic acid. Under high (20 mmol/l) glucose condition, LPC O-16:0 and LPC O-18:0 significantly potentiated GSIS in a dose-dependent manner, 50 µmol/l of LPC O-16:0 and LPC O-18:0 increased GSIS by 4- and 16-fold, respectively. On the other hand, DG 18:2_22:6, TG 12:0_18:2_22:6 and palmitic acid showed no effect on GSIS (Fig. 2).
In this study, we identified specific TG, DG, and LPC species that were independently associated with incident T2D in a phenotypically well-characterized Chinese population-based cohort with a long follow-up of ~ 16 years. These lipid predictors are reported in a prospective study for the first time and were able to achieve a substantial increase in the prediction of incident T2D beyond a predictive model based on the most important clinical risk factors identified in this cohort. We demonstrated their pathophysiological involvement in the onset of the disease and identified potential therapeutic targets for the treatment of T2D.
Our untargeted lipidomic analysis covered over 300 lipid species, which is larger than most of the previous studies in Asian populations [8, 9, 15,16,17,18]. The network analysis revealed a collective effect of TGs, DGs and PCs with a common docosahexaenoic acid chain (22:6), which was associated with increased T2D risks. TGs with a common linoleic acid (18:2), as well as LPCs and ether LPCs containing long-chain fatty acids and very-long-chain fatty acids showed a collective effect in the association with the onset of T2D. Our network analysis further showed that most of the lipid species independently associated with T2D were also inter-connected with a large number of other lipids that contain similar fatty acid chains, indicating that they may play a nexus role in the metabolic pathways of those lipids.
In this study, several LPCs were detected to be inversely associated with incident T2D after adjustments for BMI, TG/HDL ratio and the presence of prediabetes, the most important clinical predictors identified with the stepwise model selection approach. Among them, LPC P-16:0 and LPC 18:1 have previously been reported to be associated with T2D [11, 30], while our study is the first to show the inverse associations of LPC O-16:0 and LPC O-18:0 with incident T2D. Our findings are in agreement with previous cross-sectional studies which reported a significantly lower level of LPC O-16:0 in the skeletal muscle of T2D patients  and a negative association between plasma LPC O-18:0 and dysglycemia . LPC O-16:0 and LPC O-18:0 are ether-linked LPCs which are also known as lyso-platelet-activating factors (lyso-PAFs). Three other lyso-PAFs, LPC O-22:0, O-24:1 and O-24:2, have been shown to be associated with a decreased risk of T2D in a prospective study conducted in the Australian Diabetes, Obesity and Lifestyle Study . Lyso-PAF has long been considered as an inactive precursor of PAF, but emerging evidence suggests that lyso-PAF has its unique functions such as inhibiting the PAF-potentiated NADPH oxidase activation in neutrophils . However, the role of lyso-PAFs in the pathogenesis of T2D remains unknown. Our functional study demonstrated that LPC O-16:0 and LPC O-18:0 dose-dependently potentiated GSIS in mouse islets, thereby providing evidence for their potential effect on β-cell function. It has been observed that Asians with T2D tend to have prominent impairment in β-cell function, better insulin sensitivity and lower BMI compared to Caucasians [28, 29, 35]. In this collection, our study raised the possibility that lower LPC O-16:0 and LPC O-18:0 levels may contribute to impaired insulin secretion and T2D in Chinese. This is further supported by our observation that LPC O-16:0 and LPC O-18:0 potently stimulated GSIS at the physiological range whereas the same concentrations of palmitic acid, were unable to demonstrate such an effect. Indeed, consistent with our findings, previous studies reported that palmitic acid stimulated GSIS only at a high concentration of 500 µmol/l, which is out of the physiological range [36, 37]. Such observation suggests that LPC O-16:0 and LPC O-18:0 may stimulate GSIS through a possibly unknown molecular pathway independent of the classical glycerolipid/free fatty acid pathway . Further investigations are warranted to explore the signaling pathways whereby LPC O-16:0 and LPC O-18:0 potentiate GSIS and to explore their therapeutic potential for T2D.
We also detected independent associations of TG 12:0_18:2_22:6, TG 16:0_11:1_18:2, TG 49:0, TG 51:1 and DG 18:2_22:6 with increased risk of incident T2D. Our findings were consistent with previous clinical studies showing positive associations of TGs with the risk of T2D [10, 13]. In contrast, some other studies suggested that certain TGs were inversely associated with the risk of T2D [9, 11, 14]. However, these studies were conducted over relatively shorter follow-up durations, thus the cumulative detrimental effect of TGs was not prominent. Indeed, previous studies which are consistent with our findings also had a longer follow-up, with an average of 20 years [10, 13]. Our present study identified TG 12:0_18:2_22:6 and DG 18:2_22:6, which are structurally similar, as novel lipid species associated with incident T2D. Although the exact role of these two lipid species in T2D is unknown, Szendroedi et al. previously demonstrated a strong positive relationship of several DG lipid species (DG 16:0_18:2, 18:1_18:2, 18:2_18:2, and 18:2_18:0) comprised of a common α-linoleic acid (fatty acid 18:2) chain with protein kinase C theta (PKC-θ) activation in human muscle cells . The activated PKC-θ would lead to a decrease in insulin-stimulated insulin receptor substrate-1 (IRS-1)/IRS-2 tyrosine phosphorylation and subsequently disturbed downstream insulin signaling, inducing insulin resistance in the muscle. In our functional analysis, we did not observe a direct effect of the newly identified TG 12:0_18:2_22:6 and DG 18:2_22:6 on GSIS. It is possible that these lipid species may affect the development of T2D via a pathway independent of insulin secretion. Previous studies have reported that the downstream products of TGs, such as acyl-CoA, DG and ceramide, impair the insulin signaling pathways thereby leading to T2D risk . Further functional analyses to elucidate the potential role of these novel lipid species in T2D are warranted.
With the addition of only 9 lipid species, we were able to achieve a significantly increased predictive power compared to the clinical predictive model comprising of BMI, TG/HDL ratio and the presence of prediabetes (enhanced by 3.8%, p = 0.0026). Our results appeared to show a better improvement in prediction with the addition of lipid species compared to previous studies which commonly increased the prediction value by 1.0–3.2%, even with over 50 lipid species [10, 11].
The major strength of this study was the use of a phenotypically well-characterized prospective population-based CRISPS cohort which has a relatively long follow-up period compared to most of the population-based studies among the Asian population. Another strength was the use of the high-throughput untargeted lipidomic profiling, which is a robust, non-biased strategy for lipidomic depiction and discovery of novel biomarkers. On the other hand, targeted lipidomic analyses have higher sensitivity, specificity, and quantification ability, which is applicable to known compounds. Merging of untargeted and targeted analyses provides an alternative way to combine the advantages of novel biomarker discovery and accurate quantification. Given the enormous complexities of the human plasma lipidome, UHPLC-mass spectrometry based-lipidomics is the preferred method for in-depth studies of lipid-related pathological mechanisms and the identification of predictive biomarkers for diabetes. Moreover, the employment of untargeted lipidomics, which has provided comprehensive coverage of various categories of lipids, has allowed us to analyze the interactions among lipids using WGCNA and select potential factors from a large pool of lipids using a machine-learning method.
We acknowledge that the limitations of our study included the relatively small sample size and lack of an external validation cohort. Furthermore, hemoglobin A1c (HbA1c) was not evaluated in our cohort because the 1998 WHO diagnostic criteria was employed in CRISPS-2 (2000–2004), when the baseline data for the current study were collected. Nonetheless, recent studies showed that HbA1c only added 0.5% to the prevalence of diabetes diagnosed by OGTT in the Chinese population . Finally, all the participants were of Chinese ancestry, and hence our findings may not be generalizable to other populations.
In conclusion, our results provided novel insights into the underlying mechanisms of lipid species in the development of T2D. We discovered several novel lipid species associated with the pathophysiological changes before T2D onset. The identified lipid species substantially increased the predictive value beyond the traditional risk factors. Further validations in independent prospective population-based cohorts are required to confirm our findings. Future functional studies to elucidate the mechanistic pathways of the lipid species in T2D are warranted.
Availability of data and materials
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
Sun H, Saeedi P, Karuranga S, Pinkepank M, Ogurtsova K, Duncan BB, et al. IDF diabetes atlas: global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Res Clin Pract. 2022;183: 109119.
Yazıcı D, Sezer H. Insulin resistance, obesity and lipotoxicity. Adv Exp Med Biol. 2017;960:277–304.
Unger RH, Zhou YT, Orci L. Regulation of fatty acid homeostasis in cells: novel role of leptin. Proc Natl Acad Sci USA. 1999;96(5):2327–32.
Meex RCR, Blaak EE, van Loon LJC. Lipotoxicity plays a key role in the development of both insulin resistance and muscle atrophy in patients with type 2 diabetes. Obes Rev. 2019;20(9):1205–17.
Imai Y, Cousins RS, Liu S, Phelps BM, Promes JA. Connecting pancreatic islet lipid metabolism with insulin secretion and the development of type 2 diabetes. Ann N Y Acad Sci. 2020;1461(1):53–72.
Yea K, Kim J, Yoon JH, Kwon T, Kim JH, Lee BD, et al. Lysophosphatidylcholine activates adipocyte glucose uptake and lowers blood glucose levels in murine models of diabetes. J Biol Chem. 2009;284(49):33833–40.
Chamroonkiadtikun P, Ananchaisarp T, Wanichanon W. The triglyceride-glucose index, a predictor of type 2 diabetes development: a retrospective cohort study. Prim Care Diabetes. 2020;14(2):161–7.
Yun H, Sun L, Wu Q, Zong G, Qi Q, Li H, et al. Associations among circulating sphingolipids, beta-cell function, and risk of developing type 2 diabetes: a population-based cohort study in China. PLoS Med. 2020;17(12): e1003451.
Lu J, Lam SM, Wan Q, Shi L, Huo Y, Chen L, et al. High-coverage targeted lipidomics reveals novel serum lipid predictors and lipid pathway dysregulation antecedent to type 2 diabetes onset in normoglycemic chinese adults. Diabetes Care. 2019;42(11):2117–26.
Fernandez C, Surma MA, Klose C, Gerl MJ, Ottosson F, Ericson U, et al. Plasma lipidome and prediction of type 2 diabetes in the population-based malmo diet and cancer cohort. Diabetes Care. 2020;43(2):366–73.
Razquin C, Toledo E, Clish CB, Ruiz-Canela M, Dennis C, Corella D, et al. Plasma lipidomic profiling and risk of type 2 diabetes in the PREDIMED trial. Diabetes Care. 2018;41(12):2617–24.
Muilwijk M, Goorden SMI, Celis-Morales C, Hof MH, van der Vlugt KG, Beers-Stet FS, et al. Contributions of amino acid, acylcarnitine and sphingolipid profiles to type 2 diabetes risk among South-Asian Surinamese and Dutch adults. BMJ Open Diabetes Res Care. 2020. https://doi.org/10.1136/bmjdrc-2019-001003.
Merino J, Leong A, Liu CT, Porneala B, Walford GA, von Grotthuss M, et al. Metabolomics insights into early type 2 diabetes pathogenesis and detection in individuals with normal fasting glucose. Diabetologia. 2018;61(6):1315–24.
Rhee EP, Cheng S, Larson MG, Walford GA, Lewis GD, McCabe E, et al. Lipid profiling identifies a triacylglycerol signature of insulin resistance and improves diabetes prediction in humans. J Clin Invest. 2011;121(4):1402–11.
Lu Y, Wang Y, Zou L, Liang X, Ong CN, Tavintharan S, et al. Serum lipids in association with type 2 diabetes risk and prevalence in a Chinese population. J Clin Endocrinol Metab. 2018;103(2):671–80.
Chew WS, Torta F, Ji S, Choi H, Begum H, Sim X, et al. Large-scale lipidomics identifies associations between plasma sphingolipids and T2DM incidence. JCI Insight. 2019. https://doi.org/10.1172/jci.insight.126925.
Sun L, Liang L, Gao X, Zhang H, Yao P, Hu Y, et al. Early prediction of developing type 2 diabetes by plasma acylcarnitines: a population-based study. Diabetes Care. 2016;39(9):1563–70.
Lu Y, Wang Y, Ong CN, Subramaniam T, Choi HW, Yuan JM, et al. Metabolic signatures and risk of type 2 diabetes in a Chinese population: an untargeted metabolomics study using both LC-MS and GC-MS. Diabetologia. 2016;59(11):2349–59.
Janus ED. Epidemiology of cardiovascular risk factors in Hong Kong. Clin Exp Pharmacol Physiol. 1997;24(12):987–8.
Lui DTW, Lee CH, Woo YC, Fong CHY, Tso AWK, Cheung BMY, et al. Cohort profile: The Hong Kong Cardiovascular Risk Factor Prevalence Study (CRISPS) and the follow-up studies. Int J Epidemiol. 2021;50(4):1069-h.
Alberti KG, Zimmet PZ. Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet Med. 1998;15(7):539–53.
Singh B, Saxena A. Surrogate markers of insulin resistance: a review. World J Diabetes. 2010;1(2):36–47.
Halldorsson A, Magnusson CD, Haraldsson GG. Chemoenzymatic synthesis of structured triacylglycerols by highly regioselective acylation. Tetrahedron. 2003;59(46):9101–9.
Matyash V, Liebisch G, Kurzchalia TV, Shevchenko A, Schwudke D. Lipid extraction by methyl-tert-butyl ether for high-throughput lipidomics. J Lipid Res. 2008;49(5):1137–46.
Zhang H, Shao X, Zhao H, Li X, Wei J, Yang C, et al. Integration of metabolomics and lipidomics reveals metabolic mechanisms of triclosan-induced toxicity in human hepatocytes. Environ Sci Technol. 2019;53(9):5406–15.
Acharjee A, Larkman J, Xu Y, Cardoso VR, Gkoutos GV. A random forest based biomarker discovery and power analysis framework for diagnostics research. BMC Med Genomics. 2020;13(1):178.
Pencina MJ, D’Agostino RB Sr, Steyerberg EW. Extensions of net reclassification improvement calculations to measure usefulness of new biomarkers. Stat Med. 2011;30(1):11–21.
Yabe D, Seino Y, Fukushima M, Seino S. β cell dysfunction versus insulin resistance in the pathogenesis of type 2 diabetes in East Asians. Curr Diab Rep. 2015;15(6):602.
Qian L, Xu L, Wang X, Fu X, Gu Y, Lin F, et al. Early insulin secretion failure leads to diabetes in Chinese subjects with impaired glucose regulation. Diabetes Metab Res Rev. 2009;25(2):144–9.
Shi L, Brunius C, Lehtonen M, Auriola S, Bergdahl IA, Rolandsson O, et al. Plasma metabolites associated with type 2 diabetes in a Swedish population: a case-control study nested in a prospective cohort. Diabetologia. 2018;61(4):849–61.
Diamanti K, Cavalli M, Pan G, Pereira MJ, Kumar C, Skrtic S, et al. Intra- and inter-individual metabolic profiling highlights carnitine and lysophosphatidylcholine pathways as key molecular defects in type 2 diabetes. Sci Rep. 2019;9(1):9653.
Yin X, Willinger CM, Keefe J, Liu J, Fernandez-Ortiz A, Ibanez B, et al. Lipidomic profiling identifies signatures of metabolic risk. EBioMedicine. 2020;51: 102520.
Meikle PJ, Wong G, Barlow CK, Weir JM, Greeve MA, MacIntosh GL, et al. Plasma lipid profiling shows similar associations with prediabetes and type 2 diabetes. PLoS ONE. 2013;8(9): e74341.
Welch EJ, Naikawadi RP, Li Z, Lin P, Ishii S, Shimizu T, et al. Opposing effects of platelet-activating factor and lyso-platelet-activating factor on neutrophil and platelet activation. Mol Pharmacol. 2009;75(1):227–34.
Iwahashi H, Okauchi Y, Ryo M, Noguchi M, Morita S, Kishida K, et al. Insulin-secretion capacity in normal glucose tolerance, impaired glucose tolerance, and diabetes in obese and non-obese Japanese patients. J Diabetes Investig. 2012;3(3):271–5.
Cen J, Sargsyan E, Bergsten P. Fatty acids stimulate insulin secretion from human pancreatic islets at fasting glucose concentrations via mitochondria-dependent and -independent mechanisms. Nutr Metab. 2016;13(1):59.
Gravena C, Mathias PC, Ashcroft SJ. Acute effects of fatty acids on insulin secretion from rat and human islets of Langerhans. J Endocrinol. 2002;173(1):73–80.
Campbell JE, Newgard CB. Mechanisms controlling pancreatic islet cell function in insulin secretion. Nat Rev Mol Cell Biol. 2021;22(2):142–58.
Szendroedi J, Yoshimura T, Phielix E, Koliaki C, Marcucci M, Zhang D, et al. Role of diacylglycerol activation of PKCθ in lipid-induced muscle insulin resistance in humans. Proc Natl Acad Sci USA. 2014;111(26):9597–602.
Muoio DM, Newgard CB. Mechanisms of disease: Molecular and metabolic mechanisms of insulin resistance and beta-cell failure in type 2 diabetes. Nat Rev Mol Cell Biol. 2008;9(3):193–205.
Wang L, Gao P, Zhang M, Huang Z, Zhang D, Deng Q, et al. Prevalence and ethnic pattern of diabetes and prediabetes in China in 2013. JAMA. 2017;317(24):2515–23.
The authors thank all the study participants, clinical and research staff of the CRISPS cohort study for their contribution to this research study.
This work was supported by Hong Kong Research Grants Council Area of Excellence (AoE/M/707-18) and National Natural Science Foundation of China (82161138026 and 82070860).
Ethics approval and consent to participate
All participants had given written informed consent. Ethics approval was obtained from the Institutional Review Board of the University of Hong Kong/Hospital Authority, Hong Kong West Cluster.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Methods. Figure S1. WGCNA of the lipid species profile. Figure S2. Boruta analysis result. Figure S3. Violin Plot of the Targeted Lipidomics. Figure S4. Correlation among identified lipid species.
: Table S1. Stepwise model selection coefficients and Wald-Z tests. Table S2. List of measured lipid species. Table S3. Weighted gene co-expression network analysis module components. Table S4. Boruta analysis result.
About this article
Cite this article
Zhong, J., Cheung, C.Y.Y., Su, X. et al. Specific triacylglycerol, diacylglycerol, and lyso-phosphatidylcholine species for the prediction of type 2 diabetes: a ~ 16-year prospective study in Chinese. Cardiovasc Diabetol 21, 234 (2022). https://doi.org/10.1186/s12933-022-01677-4