Identification of circulating sphingosine kinase-related metabolites for prediction of type 2 diabetes

Sphingosine Kinase (SphK) that catalyzes sphingosine (Sph) to sphingosine 1-phosphate (S1P), plays a key role in both sphingolipid metabolism and cellular signaling. While SphK has been implicated in type 2 diabetes mellitus (T2DM), it is unexplored in humans. Herein, we investigated whether circulating SphK-related metabolites are associated with T2DM incidence in an established prospective cohort. Levels of SphK-related sphingolipid metabolites, including Sph, S1P, dihydrosphingosine (dhSph) and dihydro-S1P (dhS1P) in serum were measured by targeted-lipidomic analyses. By accessing to an established prospective cohort that involves a total of 2486 non-diabetic adults at baseline, 100 subjects who developed T2DM after a mean follow-up of 4.2-years, along with 100 control subjects matched strictly with age, sex, BMI and fasting glucose, were randomly enrolled for the present study. Comparison with the control group, medians of serum dhS1P and dhS1P/dhSph ratio at baseline were elevated significantly prior to the onset of T2DM. Each SD increment of dhS1P and dhS1P/dhSph ratio was associated with 53.5% and 54.1% increased risk of incident diabetes, respectively. The predictive effect of circulating dhS1P and dhS1P/dhSph ratio on T2DM incidence was independent of conventional risk factors in multivariate regression models. Furthermore, combination of serum dhS1P and dhS1P/dhSph ratio with conventional clinical indices significantly improved the accuracy of T2DM prediction (AUROC, 0.726), especially for normoglycemic subjects (AUROC, 0.859). Circulating levels of dhS1P and dhS1P/dhSph ratio are strongly associated with increased risk of T2DM, and could serve as a useful biomarker for prediction of incident T2DM in normoglycemic populations.


Introduction
Type 2 diabetes mellitus (T2DM) is a major global health problem that currently affects 463 million people over the world [1]. It is estimated that about 4.2 million deaths were attributable to diabetes in 20-79 years old adults, which accounted for 11.3% of total deaths globally in 2019 [2]. The most effective way to reduce such a global health burden is to prevent diabetes itself, because prevention of incident diabetes could reduce the all-cause mortality and the death from cardiovascular disease by 38.8%, and 44.2%, respectively [3]. However, the current strategy of prevention is still insufficient to curb the escalating epidemics of T2DM [1]. It is widely believed that both energy-dense diet and sedentary lifestyle are the primary causes of the epidemics of T2DM, however a substantial rise of incident diabetes was also seen in people with normal body weight or from the countries that have no major changes in the category of food consumption in recent decades [4]. Therefore, new approaches are urgently needed for better prediction and prevention of diabetes in general populations.
Sphingolipids are a class of essential lipids that function as both cell membrane constituents and signaling messengers [5]. Sphingosine kinase (SphK) is a key enzyme in sphingolipid metabolic network, which catalyzes the conversion of ceramides/sphingosine (Sph) to sphingosine 1-phosphate (S1P). Because the upstream substrates and the products of this enzyme often exhibit opposing biological functions, SphK is regarded as a "switch" of the sphingolipid rheostat [6]. There are two isoforms of SphK in human, i.e., SphK1 and SphK2, encoded by two different genes [6,7]. Both SphK1 and SphK2 catalyzes the conversion of Sph to S1P, while SphK2 can preferably catalyze the formation of dihydro-S1P (dhS1P) from dihydrosphingosine (dhSph) in mammals [7]. Besides, SphK1 and SphK2 often possess distinct tissue distribution, subcellular localization and biochemical properties, and thus exhibit different and even opposite functions in a context-dependent manner [7,8]. For instance, we have recently reported that global deletion of Sphk1 in mice resulted in a significant loss of pancreatic β-cell mass and promoted the onset of diabetes [9], whereas Sphk2 deficiency ameliorated diabetic phenotype by protecting β-cells against lipoapoptosis [10]. These findings indicate that both SphK1 and SphK2 are critically involved in the regulation of β-cell viability and function in opposite directions. Moreover, the implications of SphK1 and SphK2 in insulin resistance have recently emerged. For examples, Anderson et al. reported that adipocyte-specific deletion of Sphk1 in mice impaired glucose tolerance and elevated fasting blood glucose, which is believed attributable to the suppressed lipolysis in adipocytes [11]. More recently, we have reported that hepatocytespecific deficiency of SphK2 in mice increased hepatic glucose production, impaired glucose homeostasis and caused insulin resistance [12]. Thus, it appears that SphK is critically involved in both insulin resistance and β-cell failure, suggesting a crucial role in diabetes. However, the relationships of SphK and its related sphingolipid metabolites with diabetes in humans remain unexplored. Therefore, the present study aimed to investigate the association of SphK-related metabolites with the risk of T2DM in a community-based prospective cohort.

Study samples
For the animal study, the colonies of Sphk1 −/− , Sphk2 −/− and wild-type littermates were derived on the same C57BL/6 background as detailed in our previous studies [9,10]. All mice were housed in a temperature-controlled, pathogen-free environment on a 12-h light/dark cycle and had ad libitum access to food and water. Animal studies were approved by the Animal Use and Care Committee of Fudan University and conformed with the U.S. National Institutes of Health (NIH) U.S. Department of Health and Human Services Public Health Service Policy on Humane Care and Use of Laboratory Animals (NIH publication No. 15-8013; Bethesda, MD, USA).
For the cohort study, the participants were enrolled from the Shanghai Changfeng Study, a populationbased prospective cohort study on metabolic diseases in Chinese community residents aged over 45 years as previously reported in details [13,14]. A retrospective nested case-control study was conducted. Among 2486 non-diabetic participants at baseline, a total of 100 participants with newly diagnosis of T2DM after a mean follow-up of 4.2 years (3.5-5 years) were randomly selected and each case was matched with an individual without incident diabetes from the same cohort. All cases and control subjects were strictly matched with respects to age, sex, body mass index (BMI; ± 1 kg/ m 2 ) and fasting plasma glucose (FPG; ± 0.1 mmol/L) at baseline. All study participants provided written informed consent, and the study was approved by the Ethical Committee of Zhongshan Hospital affiliated to Fudan University (No. 2008-119).

Anthropometric and biochemical profiling
Weight, height and waist circumference were measured by trained interviewers. BMI was calculated by dividing the weight (kg) by the square of the height (m 2 ). The mean value of three resting blood pressure (BP) measurements was used for the analyses. Venous blood samples were collected for the biochemical examinations after an overnight fast of at least 12 h. Serum total cholesterol, highdensity lipoprotein (HDL) cholesterol, triglycerides (TG) levels and uric acid were measured by an oxidase method. Low-density lipoprotein (LDL) cholesterol was calculated using the Friedewald equation. All participants underwent a 75 g oral glucose tolerance test (OGTT). Fasting and 2-h post-load plasma glucose (PPG) were measured using a glucose oxidase method. Serum insulin concentrations were measured using electrochemiluminescence immunoassay. Homeostatic Model Assessment of Insulin Resistance (HOMA-IR) was calculated by multiplying the FPG (mmol/L) by fasting insulin (μU/mL) and dividing by 22.5 [15]. For both the baseline and follow-up evaluations, diabetes was diagnosed if FPG ≥ 7.0 mmol/l or OGTT 2 h PPG ≥ 11.1 mmol/l, and normal glucose tolerance defined as FPG < 6.1 mmol/L and OGTT 2 h PPG < 7.8 mmol/L, according to the WHO criteria [16].

Sphingolipid analyses by HPLC-MS/MS
For the measurement of Sph, dhSph, S1P and dhS1P, the serum samples collected from the study subjects were removed from − 80 °C storage just before the measurement. Sphingolipids were measured in duplicate of each sample (100 μL each), using a High Performance Liquid Chromatography-Tandem Mass Spectrometry (HPLC-MS/MS) methodology as previously described with some modifications [9,17]. Briefly, the samples were mixed with 10 μL of internal standards (S1P-d7; Avanti, Alabama, USA) and 500 μL of extraction solution consisting of isopropanol, methanol and formic acid (45:45:10 v/v). Then, the mixtures were vortexed, sonicated and centrifuged at 1500×g for 10 min at 4 °C. Supernatants were diluted with 30% methanol before analysis. The samples were analyzed by positive-ion Ultrahigh-Performance Liquid Chromatography-Electrospray Ionization-Tandem Mass Spectrometry (UPLC-ESI-MS) in a multiple reaction monitoring mode performed on a Shimadzu LC-20AD liquid chromatography instrument (Shimadzu, Japan) and coupled to an QTrap5500 mass spectrometer (Applied Biosystems Sciex, Canada). Chromatographic separations were performed under a gradient elution on an Eclipse XDB-C8 (2.1 × 100 mm, 1.8 μm) column (Agilent Technologies) using a binary solvent system at a flow rate of 300 μl/min. Prior to sample injection, the column was equilibrated for 2 min with a solvent mixture of 70% mobile phase A (methanol/H2O/formic acid, 69/29/2/, v/v/v, with 2 mmol ammonium formate) and 30% phase B (methanol/tetrahydrofuran/formic acid, 97/2/1, v/v/v, with 2 mmol ammonium formate) as previously described [18]. After injection, the mobile phase A/B ratio was maintained at 70/30 for 2 min, followed by a linear gradient to 100% B over 5 min. The flow was then maintained for 10 min at 100% B, followed by a wash of the column with 70:30 A/B for 0.5 min before the next run. The mass transitions were m/z 300. S1P-d7, respectively. Levels of Sph, dhSph, S1P and dhS1P were quantified using external standard curves ranging from 0.1 to 1000 nmol, and normalized to each the synthetic internal standard (Avanti, Alabama, USA). All of the data were acquired and processed using Analyst 1.6 software from Sciex Inc., Concord, ON, Canada.

Statistical analysis
The estimation of sample size was based on our preliminary experiment in 30 pairs of subjects. According to the preliminary data, subjects with incident diabetes had significantly higher log 10 dhS1P/dhSph ratio levels of 0.04 with a SD of 0.11. Therefore, we determined that a sample size of 82 achieves 90% power to detect a mean of paired differences of 0.04 with an estimated standard deviation of differences of 0.11 and a significance level (alpha) of 0.05 using PASS tests for paired means. Assuming a drop-out rate of 20%, approximately 100 pairs subjects would be required.
All statistical analyses were performed using STATIS-TICAL ANALYSIS SYSTEM (SAS) 9.2 and R 3.4.2 software (SAS Institute). The continuous data were compared using the Student's t-test, whereas the χ2 test was used for comparisons of categorical variables. Skewed variables were presented as the median with the interquartile range. The sphingolipids were non-normally distributed, so the sphingolipid parameters were log-transformed before analysis. Since the Shanghai Changfeng Cohort Study is a fixed cohort study with an average of 4.2-year (3.5-5 years) follow-up, multivariate logistic regression models were used to estimate the odds ratios (ORs) and 95%CIs of incident diabetes per-SD increase of each individual sphingolipid parameter at baseline after multiple adjustment. Model 1 was adjusted for age, sex and BMI changes as the basic model. Model 2 was adjusted for family history of diabetes, FPG and OGTT 2 h PPG in addition to factors included in Model 1. Model 3 was adjusted for serum total cholesterol, triglycerides, systolic blood pressure and uric acid in addition to factors included in Model 2. The performance of each serum sphingolipid component at baseline in prediction of incident diabetes was evaluated using receiver operating characteristic (ROC) curve analyses, and was compared with conventional risk factors of diabetes. A combination model of serum sphingolipids with conventional risk factors was also established and assessed. Subgroup analysis was performed in the participants with normal glucose tolerance. All tests were two-sided and a P value of less than 0.05 was considered statistically significant.

Circulating levels of the SphK-related sphingolipid metabolites in mice
Because no established methodology currently available for measuring SphK activity in vivo, we applied Sphk1-deficient (Sphk1 −/− ) and Sphk2-deficient (Sphk2 −/− ) mice to measure circulating levels of the key sphingolipid metabolites that are mostly relevant to SphK activity, including Sph, dhSph, S1P and dhS1P, using HPLC-MS/MS-based targeted-lipidomic analyses. In line with the enzymatic property of SphK, i.e., converting Sph to S1P, Sphk1 −/− mice resulted in a significant decrease in serum S1P but increase in Sph levels ( Fig. 1). Correspondingly, the ratio of circulating S1P/Sph was dramatically decreased. However, the S1P/Sph ratio was not significantly altered in Sphk2 −/− mice, as S1P levels were significantly increased in Sphk2 −/− mice. The latter was consistent with the previous reports [19,20]. By contrast, Sphk2 −/− significantly decreased dhS1P and slightly increased dhSph levels, and thus a remarkable reduction in serum dhS1P/dhSph ratio was detected. However, the circulating dhS1P/dhSph ratio was not significant altered by Sphk1 −/− , although a decrease trend existed (Fig. 1). Taken together, these results indicate that the in vivo activity of SphK1 and SphK2 is likely reflected in the circulating ratio of S1P/Sph and dhS1P/dhSph, respectively.

Circulating levels of the SphK-related sphingolipid metabolites in human
We have recently reported that global deletion of either Sphk1 or Sphk2 exhibited a significant impact on diabetes in the animal models, i.e., loss of Sphk1 promotes, but Sphk2 deficiency protects, the disease [9,10]. We thus attempted to apply these experimental findings into a clinically relevant setting and investigate whether SphK is associated with T2DM incidence in humans. To this end, we conducted targeted-lipidomic analyses of the SphK-related sphingolipids in serum collected from an established community-based prospective cohort with a Fig. 1 Circulating levels of SphK-related sphingolipids in Sphk1 −/− and Sphk2 −/− mice. Levels of circulating SphK-related sphingolipids, including sphingosine (Sph), sphingosine 1-phosphate (S1P), dihydro-Sph (dhSph) and dhS1P were measured in Sphk1 −/− , Sphk2 −/− and wild-type (WT) mice (n = 7 per group). Ratios of S1P/Sph and dhS1P/dhSph are shown in the right panels. Data are expressed as means ± SD. *P < 0.05, **P < 0.01, ***P < 0.001 nested case-control design. We randomly selected 100 cases who developed T2DM (incident T2DM) after a mean follow-up of 4.2 years, which were matched with 100 subjects who did not develop diabetes (serving as a control group) from a total of 2486 participants without diabetes at baseline in the cohort. The two groups were strictly matched in subjects' age, sex, baseline BMI and FPG levels. Baseline characteristics of the study population are shown in Table 1. In addition to the matched parameters, there were no difference in waist circumference, blood pressure, serum LDL-C, HDL-C, uric acid, 2-h PPG and HOMA-IR between the two groups. However, levels of total cholesterol in the incident diabetic group were higher than the control group  (Table 2 and Fig. 2). However, there were no significant difference in serum levels of Sph, S1P and dhSph and ratio of S1P/ Sph between the control and incident diabetic group either at baseline or after follow-up.

Association of the SphK-related sphingolipids with incident T2DM
Having shown the elevated level of dhS1P and dhS1P/ dhSph ratio occurred prior to diabetes onset in the case group, we then further evaluated the independent association of each sphingolipid component at baseline with the risk of incident diabetes by multiple logistic regression analyses. As shown in Table 3, univariate regression analysis revealed that levels of serum dhS1P and dhS1P/ dhSph ratio, but not others, were significantly and positively associated with increased risk of incident diabetes, and each SD increment of dhS1P and dhS1P/dhSph ratio was associated with 53.5% (OR, 1.54; 95%CI, 1.12-2.10) and 54.1% (OR, 1.54; 95%CI, 1.10-2.17) increased risk of incident T2DM, respectively. Parameter estimates of the correlation between baseline dhS1P level and dhS1P/ dhSph ratio and risk of incident diabetes remained significant after successively adjusted for age, sex, BMI changes, family history of diabetes, FPG, 2-h PPG, as well as the addition of total cholesterol, triglycerides, systolic blood pressure and uric acid levels (dhS1P: OR, 1.50; 95%CI, 1.05-2.14; dhS1P/dhSph: OR, 1.70; 95%CI, 1.14-2.53).

Performance of the SphK-related sphingolipids for prediction of incident T2DM
In attempt to assess the potential value of serum dhS1P and dhS1P/dhSph ratio at baseline for prediction of diabetes, we performed receiver-operating-characteristic curves (ROC) analyses. The data indicate that the prediction performance of serum dhS1P and dhS1P/dhSph ratio for the incidence of diabetes was comparable to that of the combined conventional risk models, including age, BMI, LDL, HDL, triglycerides, FPG, OGTT 2 h PPG, HOMA-IR and systolic blood pressure. A combination of serum dhS1P and dhS1P/dhSph with the conventional diabetes risk factors could further improve the accuracy for prediction of incident diabetes (AUROC, 0.726; 95%CI, 0.647-0.805), in comparison with the conventional prediction models (AUROC, 0.654; 95%CI, 0.568-0.739) (Fig. 3a). Furthermore, in subgroup analysis when dhS1P and dhS1P/dhSph ratio were added to established risk prediction model of T2DM (AUROC,0.699; 95%CI, 0.590-0.793) in subjects with normal glucose tolerance at baseline, the prediction was significantly improved to an AUROC of 0.859 (95%CI, 0.767-0.924) (P = 0.002 compared with conventional diabetes prediction models (Fig. 3b).

Discussion
Early prediction of T2DM is regarded as one of the most effective ways to reduce its incidence in the general population and improve the long-term prognosis of the disease. The current study based on an established prospective cohort with a nested case-control design has for the first revealed time that elevated circulating levels of dhS1P and its ratio to dhSph were able to effectively predict diabetes incidence in otherwise healthy, normoglycemic individuals. The levels of dhS1P and dhS1P/dhSph ratio were indeed elevated before the changes of blood glucose and insulin levels in the subjects who developed T2DM after an average of 4.2-year follow-up, suggesting that these SphK-related sphingolipids could serve as a potential biomarker for early prediction of T2DM. Previous human studies have reported that several sphingolipids, especially ceramides, were associated with Table 2 General characteristics and levels of SphK-related sphingolipids at follow-up in the control and incident T2DM groups Data are presented as mean ± SD, except for skewed variables, which are presented as the median with the interquartile range given in parentheses   insulin resistance (reviewed in [21,22]). Based on a large family-based cohort study, the Strong Heart Family Study of American Indians, it was found that higher levels of serum ceramides with saturated fatty acids are associated with higher fasting serum insulin and HOMA-IR, both cross-sectionally and prospectively [23]. More recently, by using large-scale lipidomic analysis, studies on 2 ethnically distinct cohorts have shown that ceramides correlated positively whereas hexosylceramides correlated negatively with clinical and biochemical characteristics of obesity and diabetes [24,25]. Notably, ceramides in serum are actually present in a broad variation of subspecies, which consist of either saturated or unsaturated sphingoid bases (i.e., dhSph and Sph) with different acylchain lengths of fatty acid (C14:0-C30:0) [26]. However, only a certain subspecies of ceramides can be found associated with insulin resistance in both experimental animals and humans [27,28]. This strongly suggests the notion that the balance of sphingolipid metabolism, rather than ceramide accumulation, in blood is associated with insulin resistance and diabetes [24,25,27,28]. While the regulation of circulating level of sphingolipids has yet been identified, recent studies have suggested that the liver is a major contributor of circulating ceramide species [29]. More recently, by taking the advantage of studies on hepatocyte-specific Sphk2-deficient mice, we found that the liver is also a key organ that regulates circulating levels of sphingolipids, including Sph, dhSph, S1P and dhS1P [12]. Remarkably, the liver-specific deficiency of Sphk2 exhibits pronounced insulin resistance and glucose intolerance, demonstrating a key role of hepatic SphK2 in the regulation of insulin sensitivity and glucose homeostasis. Similar to SphK1, SphK2 is also responsible for the conversion of Sph to S1P in cells. However, circulating S1P appears chiefly attributable to SphK1 activity [30,31]. Indeed, we found that global deletion of Sphk1 in mice results in a more than 50% reduction of S1P in serum, whereas loss of Sphk2 unexpectedly increases circulating S1P (Fig. 1), which is consistent with previous reports [19,20]. The circulating S1P is highly dynamic and its balance is regulated by a process that involves de-phosphorylation of S1P by S1P-phosphatase, re-phosphorylation of Sph by SphK2 and irreversible degradation of the newly formed S1P by its specific lyase [32]. Unlike SphK1, SphK2 can also specifically catalyze the phosphorylation of dhSph to form dhS1P in mammals [33]. Correspondingly, while deletion of Sphk2 resulted in a significant decrease in dhS1P and increased dhSph in serum, neither dhS1P nor dhSph was influenced by Sphk1 deficiency in mice (Fig. 1). In unstimulated cells, SphK1 is a cytosolic enzyme, whereas SphK2 can be found in the nucleus, the cytoplasm, and other intracellular compartments. Although it was reported that a minor fraction of SphK1 is constitutively exported from endothelial cells, possibly contributing to circulating S1P [34], we were unable to detect the enzymatic activity of either SphK1 or SphK2 in murine and human serum. Thus, the findings that circulating S1P/Sph and dhS1P/ dhSph ratios were dramatically decreased by deficiency of Sphk1 and Sphk2, respectively, indicate that the ratios could represent a favorable readout of in vivo activity for these two isoenzymes. With nearly 4000 sphingolipid metabolites known to be present in mammalian cells (http:// www. lipid maps. org), many of which are bioactive and play important roles in biological function. Despite recent progress in mass spectrometry-based technologies that allow for better discrimination of individual sphingolipid species than was previously possible, only a relatively small portion of sphingolipids can be detected in blood by non-targeted sphingolipidomic or large-scale lipidomic analyses [35,36]. The current study was driven primary by our previous experimental data that suggest an important role of SphK in the pathogenesis of diabetes. Thus, we aimed to explore the relationship between sphingolipids and T2DM with a focus on SphK-related metabolites, including Sph, dhSph, S1P and dhS1P. The most significance of the current study is the finding that high levels of circulating dhS1P and its ratio to dhSph are prospectively associated with T2DM incidence, revealing a novel biomarker for susceptibility to develop diabetes, up to 5 years before the disease diagnosed. The prediction performance of circulating dhS1P/dhSph ratio is better than dhS1P alone. In marked contrast, neither S1P levels nor S1P/Sph ratio is associated with the risk of incident diabetes. To the best of our knowledge, this is the first study showing the association of T2DM incidence with circulating levels of SphK-related sphingolipids in a prospective epidemiologic setting. Early identification of individuals with high risk of incident diabetes is of critical importance, because diabetes and its complications have been fully proven to be preventable [37,38]. Although conventional risk factors, such as FPG, BMI and dyslipidemia, could predict incident diabetes to a certain degree [39], their prediction accuracy and sensitivity are limited. A better clinical prediction model for incident diabetes remains urgently needed. Previous efforts on combining the genetic risk factors of T2DM with conventional predictors add modestly values to the risk assessment [40]. In the present study, the application of serum dhS1P/dhSph ratio in early prediction of incident diabetes showed comparable accuracy to the conventional risk model before the changes of blood glucose. Combination of dhS1P levels and dhS1P/dhSph ratio with those conventional risk factors might provide a promising tool for early prediction of T2DM clinically.
The finding that dhS1P and dhS1P/dhSph ratio are elevated in serum up to 5 years before diabetes onset suggests that SphK2 activity could be stably associated with T2DM predisposition. Indeed, in recent years, an increasing body of evidence has suggested an important role of SphK2 in both insulin resistance and pancreatic β-cell dysfunction, the two key pathogenic factors of diabetes [10,12,41]. While SphK1 has been identified as protective factor against β-cell death [9], we have recently uncovered that SphK2 is an endogenous pro-apoptotic factor promoting β-cell lipotoxicity and dysfunction [10]. Global ablation of Sphk2 significantly reserves β-cell mass and insulin production, resulting in amelioration of diabetic phenotype in a mouse model of T2DM [10]. Intriguingly, SphK2 is likely a metabolically protective factor in the liver, as Sphk2 −/− predisposes to non-alcoholic and alcoholic fatty liver diseases [42,43], whereas adenoviral overexpression of SphK2 primarily in the liver improves insulin resistance [41]. More recently, we found that hepatocyte-specific Sphk2-deficent mice exhibited evident insulin resistance and pre-diabetes [12]. Thus, it appears that SphK2 can exert both protective and pathogenic effects, leading to different metabolic outcomes in a highly context-dependent manner. Although we are currently unable to assess the contribution of these different functions of SphK2 to the overall metabolic homeostasis, the association of circulating dhS1P/dhSph ratio with T2DM incidence has shed light on an important role of SphK2 in diabetes in a clinically relevant setting.
In this study, all cases and control subjects were strictly matched with respects to age, sex, BMI, FPG at baseline. Intriguingly, levels of total cholesterol in the incident diabetic group were higher than the control group. It is well acknowledged that total cholesterol is an established risk factor for diabetes. To exclude the potential confounding factor of cholesterol, we first performed a correlation analysis and found that dhS1P and dhS1P/dhSph ratio were not related to total cholesterol (P > 0.10). Moreover, in the logistic regression, after adjustment for total cholesterol, parameter estimates of the correlation between baseline dhS1P level and dhS1P/dhSph ratio and risk of incident diabetes remained significant. Hence, it is more likely to be a happenchance but not a real phenomenon associated with the changes in SphK-related sphingolipid metabolite levels.
It is interesting to note that high levels of serum dihydroceramides have recently been identified as a potential biomarker of T2DM susceptibility in both mice and humans [28]. Dihydroceramides, along with dhSph and dhSph, represent a key subspecies of sphingolipids, namely dihydrosphingolipids, which does not contain a C4-trans-double bond in the sphingoid backbone. They are produced in early steps of the de novo ceramide synthetic pathway (reviewed in [44]). Dihydroceramides are yielded from the conversion of dhSph by (dihydro-) ceramide synthases, and the latter can be phosphorylated by SphK2 to form dhS1P. Theoretically, loss of Sphk2 leads to dhSph accumulation, which would increase dihydroceramide formation. Likewise, elevated SphK2 activity and increased dhS1P would result in a reduction of dihydroceramides. However, as dihydroceramides were not measured in the present study, we are unable to ascertain whether the individuals who had high levels of dhS1P with incident T2DM also had changes in dihydroceramides. This is worth to be tested in the future.
There are several limitations in the present study. First, while the nested case-control design has provided us a powerful tool to identify serum dhS1P/dhSph ratio as a biomarker for T2DM prediction in a given cohort, we did not evaluate its performance in additional large-scale population of cohorts due to funding constraints. Thus, the predictive value of serum dhS1P and the ratio of dhS1P/dhSph for incident T2DM in general population remains to be further validated. Second, our study participants were only of Han Chinese, and thus the further studies on a variety of races and ethnicities are required. Third, all participants in this study were middle-aged and elderly adults, and further work is needed to determine whether our findings can be extrapolated to all-ages population. Finally, in the present study we only measured circulating levels of the sphingolipid metabolites mostly related to SphK activity, including Sph, S1P, dhSph and dhS1P. Further work is needed to examine whether combining the current findings with more sphingolipid measurements, such as ceramides, dihydroceramides and glucosylceramides, could further improve the predictive value for incident T2DM.
In summary, we have shown that high levels of circulating dhS1P and its ratio to dhSph are prospectively associated with T2DM incidence in a prospective cohort, revealing a novel biomarker for susceptibility to develop diabetes, up to 5 years before the disease is diagnosed. Further studies are warranted to validate its prediction power in a general population and further understand the mechanisms underlying the effects of SphK-related sphingolipids in the pathogenesis of T2DM.