Serum sphingolipid profiling as a novel biomarker for metabolic syndrome characterization

Background: Sphingolipids are components of cell membrane structure, but also circulate in serum and are essential mediators of many cellular functions. While ceramides have been proposed previously as a useful biomarker for cardiometabolic disease, the involvement of other sphingolipids is still controversial. The aim of this study was to investigate the cross-sectional association between blood sphingolipidomic profiles and metabolic syndrome (MetS) as well as other atherosclerotic risk factors in a large population-based study in the U.S.

Methods: Clinical data and serum sphingolipidomic profiling from 2,063 subjects who participated in the biomarker project of the Midlife in the United States (MIDUS) study were used.

Results: Consistent with previous reports, we found a positive association between most ceramide levels and obesity, atherogenic dyslipidemia, impaired glucose metabolism, and MetS prevalence. In contrast, most simple β-glycosphingolipids (i.e., hexosylceramides and lactosylceramides) were inversely associated with dysmetabolic biomarkers. However, this latter sphingolipid class showed a positive link with inflammatory and vascular damage-associated biomarkers in subjects with MetS. Through metabolic network analysis, we found that the relationship between ceramides and simple β-glycosphingolipids differed significantly not only according to MetS status, but also with respect to the participants' C-reactive protein levels.

Conclusion: Our findings suggest that a comprehensive sphingolipid profile is more informative about MetS than ceramides alone, and it may reveal new insights into the pathophysiology and further diabetic vs. cardiovascular risk in patients with MetS.

Introduction

Atherosclerotic cardiovascular diseases (ASCVD) are still a major leading cause of death worldwide (1). Several risk factors have been well-established, including family history, obesity, high blood pressure, dyslipidemia, type 2 diabetes, proinflammatory processes, and metabolic syndrome (MetS) (1). MetS is typically identified by several disease-related risk factors (i.e., abdominal obesity, high serum triglyceride and glucose, hypertension, and low HDL cholesterol (HDL-c) (2) and meeting these criteria increases risk for ASCVD and type 2 diabetes by two- and five-fold, respectively (3).

Patients with MetS show several lipid abnormalities beyond the well-known proatherogenic lipid profile (increased triglycerides and low HDL-c), including high levels of small, dense and modified-LDL particles (4). In fact, traditional biochemical measures alone do not explain the complexity of metabolic abnormalities associated with MetS or its related atherosclerotic cardiovascular risk (CVR) vs. diabetic risk (4). Thus, there is a lingering need to delineate other lipid biomarkers that may contribute to improve evaluation, prognosis, and treatment of this syndrome and its long-term complications.

Based on technological advances in lipidomics, researchers have identified sphingolipids as novel biomarkers of cardiometabolic risk (5, 6). However, many of the underlying mechanisms involved remain unclear and vary depending on the specific class, or even species.

Sphingolipid metabolism can be characterized as a group of metabolically connected pathways all of which converge toward ceramides (Figure 1). A substantial literature now shows that some species of ceramides may have a causative and pathogenic role in diabetes and other cardiometabolic disorders (7–12). On the other hand, the role of simple glycosylated ceramides, such as hexosylceramides and lactosylceramides, in cardiometabolic disease is unclear. While administration of these glycosphingolipids has shown remarkable beneficial effects on glucose intolerance, metabolic syndrome, and hepatic steatosis in mouse models (13–15), recent studies suggest that high levels of simple β-glycosphingolipids may increase the risk of cardiovascular events in humans (5, 16). Thus, the association of blood sphingolipid levels with MetS features may potentially reveal significant insights into the underlying mechanisms and evolution of cardiometabolic diseases in these subjects (17).

FIGURE 1

Figure 1. Sphingolipid metabolic pathway. A set of enzymes catalyzes the condensation of serine and acyl-CoA (usually palmitate) and de novo synthesis of dihydroceramides (yellow box). Dihydroceramides are then oxidized to ceramides, which are characterized by a long sphingosine carbon backbone (in black) and a fatty acid of variable length (orange box) that defines the species, according to the number of carbon atoms and double bonds. Subsequently, a hydrogen in ceramides can be replaced by a polar group, given rise to more complex sphingolipids. Depending on this polar group (green box) different classes exist, such as sphingomyelins, simple hexosylceramides (i.e., glucosylceramides and galactosylceramides), lactosylceramides, and other complex glycosphingolipids. These complex sphingolipids can be catabolized and recycled to regenerate ceramide through the salvage pathway (red box), while sphingomyelin can be hydrolyzed directly to ceramide (blue box). SPT, serine palmitoyltransferase; CerS, ceramide synthase; DES, dihydroceramide desaturase; SMase, sphingomyelinase; SMS, sphingomyelin synthase; CGT, ceramide galactosyltransferase; GCS, glucosylceramide synthase; LCS, lactosylceramide synthase, CDase, ceramidase.

The current study investigated the relationship between blood sphingolipids and various parameters of MetS and other atherosclerotic risk factors in a large population-based study from U.S. We hypothesized that distinctive sphingolipidome profiles would be evident in adults with MetS and might vary depending on whether they exhibit a predominantly dysmetabolic pattern (i.e., obesity, dyslipidemia, and impaired glucose metabolism) or indicate more pro-inflammatory activity (i.e., high levels of pro-inflammatory proteins and endothelial dysfunction).

Methods

Study sample

Blood samples for the current analysis were collected as part of the Midlife in the United States (MIDUS) study, specifically from middle-aged and older adults evaluated during MIDUS 2 and MIDUS Refresher phases. MIDUS was initiated in 1995/96 to assess the role of behavioral, psychological, and social factors in adult health (18). MIDUS 2 was a longitudinal follow-up conducted approximately 9–10 years after the baseline survey and included a comprehensive biomarker evaluation of a subsample (19). The MIDUS Refresher was intended to broaden the age and racial representation by recruiting a new set of participants after the 2007–2009 great recession (20).

From subjects who participated in the biomarker project (n = 2,118; 1,255 from MIDUS 2 and 863 from MIDUS Refresher), we excluded 55 individuals (2.6%) who did not have lipidomic data or metabolic measures. Thus, the current analysis included 2,063 subjects (54.9% female, 74.8% white of European family backgrounds, age_mean = 55.7 ± 12.5 years) who participated in MIDUS 2 (n = 1,214) and MIDUS Refresher (n = 849) biomarker assessments. Additional descriptive information for the total sample is provided in Supplementary Table 1.

Sociodemographic variables

Sociodemographic variables included age, gender, race, and years of education. Age was treated both as a parametric variable as well as categorical with participants in 3 age groups (<50, 50–65 and >65 years of age). Educational attainment was categorized into 3 levels: (1) high school or less, (2) some college, and (3) postgraduate studies. Race was self-reported and coded as White, African American, and others.

Anthropometric measures and biochemical biomarkers

Health variables included body mass index (BMI, kg/m²), waist circumference (cm), blood pressure (mmHg), LDL-c (mg/dL), HDL-c (mg/dL), triglycerides (mg/dL), oxidized HDL (arb. unit), glycemia (mg/dL), insulinemia (μIU/mL), HbA1c (%), HOMA-IR, cytokine inflammatory markers (pg/mL), and C-reactive protein (μg/mL).

BMI, waist circumference, and blood pressure were measured by a clinical staff in one of the three Clinical and Translation Research Centers, where participants were evaluated during overnight visits (i.e., University of Wisconsin-Madison, UCLA, and Georgetown University). For this analysis, overweight was defined as >25–< 30 kg/m² and BMI ≥30 kg/m² was considered as obesity. In addition, information on relevant medication use was obtained at the time of blood collection.

All biochemical measures were assessed from a fasting blood sample as described elsewhere (21, 22). Briefly, a traditional lipid panel (HDL-c, estimated LDL-c, and triglyceride levels) and measures of glucoregulation (glycemia, insulinemia, HOMA-IR, and HbA1c) were determined by CLIA-certified clinical laboratories (Meriter Labs, Madison, WI and ARUP, Salt Lake City, UT). Oxidized HDL (oxHDL) levels were determined by measuring lipid peroxidation of HDL (23). Plasma fibrinogen and high sensitive C reactive protein (CRP) levels were measured using a BNII nephelometer (Dade Behring Inc., Deerfield, IL). Serum IL-6, soluble ICAM-I (soluble intercellular adhesion molecule-1), and soluble E-selectin levels were measured using high sensitivity ELISAs (R&D Systems, Minneapolis, MN).

Serum sphingolipid profiling

Serum sphingolipid profiling was performed as part of an untargeted overall lipidomic approach by Metabolon, Inc. (Durham, NC), following the protocol described in Supplementary material (24). For the current analysis, only sphingolipid values were used, i.e., dihydroceramides, ceramides, sphingomyelins, hexosylceramides (glucosylceramides and galactosylceramides) and lactosylceramides.

Before conducting the correlational analyses, sphingolipid levels were log-transformed to achieve normal distributions and then further harmonized using z-scoring, which made the effect sizes comparable. Only those lipid species with <20% of missing values (due to being below the lower limit of detection) were included in the statistical analyses.

To calculate sphingolipid ratios, total levels of each class were used. The first ratio divides ceramide levels by lactosylceramide levels, while the second ratio partitions levels of precursor ceramides (DCER and CER) by glycosylated-ceramides (HCER and LCER). When dividing to generate sphingolipid ratios, calculations were done before log transformation.

Metabolic syndrome

MetS was defined by the updated definition of the Adult Treatment Panel III (ATPIII) of the US National Cholesterol Education Program (25). It was based on at least three abnormal findings out of the following criteria: central adiposity (waist circumference >102 cm for men and >88 cm for women), elevated triglycerides (≥150 mg/dL) or presence of pharmacological treatment for hypertriglyceridemia, reduced HDL-c (<40 mg/dL for men and <50 mg/dL for women) or use of prescription medications to raise low HDL-c, elevated blood pressure (systolic ≥130 mmHg or diastolic ≥85 mmHg) or pharmacological treatment for hypertension in patients with a history of this disease, and high fasting glucose (≥100 mg/dL) or pharmacological treatment for hyperglycemia.

Statistics

Participant characteristics were summarized as the mean ± SD for continuous variables or percent occurrence for categorical variables. Statistical analyses were conducted using R-studios Desktop, version 1.2.5001. All analyses were two-tailed and considered statistically significant when p ≤ 0.05.

Volcano plot - For exploratory analysis, a volcano plot was performed to visualize the overall pattern of sphingolipid species in participants with MetS. Fold change threshold was set to higher than 1.15 and statistical significance evaluated with the Wilcoxon test. For these analyses, lipid concentrations were used without any transformation.

Correlation analyses

All biomarker values were log-transformed to achieve graphically normal data distributions and normalized using z-scoring. Associations between sphingolipids and the health biomarkers were evaluated using linear regression models (adjusting for sex, age, race, education attainment, and MetS diagnosis) and graphed as regression coefficient matrices. Color labels represent the regression coefficients. Associations between sphingolipids and MetS components were assessed using binary logistic regression models adjusted for sex, age, race, and education attainment. Regression effects were presented as odd ratios (OR) with 95% confidence interval (CI) or graphed in a regression matrix plot. To assess the relationship between inflammatory CVR biomarkers and sphingolipid levels in subjects with or without MetS, linear regression models were used, adjusting for sociodemographic variables (sex, age, race, and education attainment).

Sphingolipid networks

All participants were categorized with respect to MetS criteria (positive or negative) and blood levels of CRP (below or above the median). Sphingolipid distance network analyses were performed including the sphingolipid species selected in the volcano plot and built based on the absolute Spearman correlation distance (1-|ρ|), setting to 0 those coefficients with p ≥ 0.05 after Benjamin-Hochberg (BH) correction. Therefore, a shorter distance represents a higher correlation magnitude between two lipid species. To identify clusters, the fast greedy modularity optimization algorithm (26) implemented in the igraph R package (27) was used to directly optimize the modularity score by merging clusters pairwise iteratively. The same method was used to analyze the networks according to alterations in glucose metabolism. For this latter analysis, impaired glucoregulation included individuals who had high glucose >100 mg/dL and insulin ≥25 mIU/L. To assess the statistical difference between complete networks, the algorithm proposed by Steiger (28) for comparison of two distance matrixes was used, implemented in the “psych” R package (29). The distance between two sphingolipid classes was calculated as the average of distances between all pairs of lipids with one species of each class. A permutation test with 1,000 repetitions was used to determine the significance of the differences in interclass distances between different groups.

Study approval

All procedures followed were in accordance with ethical standards of the Health Sciences Institutional Review Board at the University of Wisconsin-Madison, as well as the Institutional Review Boards at the University of California-Los Angeles and Georgetown University. Informed consent was obtained from all individual participants included in the study.

Results

Sociodemographic, behavioral and sphingolipid differences between participants with and without MetS

Descriptive summary statistics for sociodemographic and health variables for MIDUS participants with or without MetS are provided in Table 1 (n = 1,292 and 773, respectively), Individuals meeting MetS criteria were older, less educated, and proportionally included more men and individuals of racial minorities. They were less likely to engage in regular exercise and more healthy diet. In addition, more than 95% of the adults with MetS were overweight or obese.

TABLE 1

Table 1. Comparison of sociodemographic and health characteristics between American adults with and without metabolic syndrome (MetS).

To explore the differences in the sphingolipid profile between adults with and without MetS, an exploratory volcano plot analysis of sphingolipid species was performed (Figure 2). Most ceramide and dihydroceramide species were higher in individuals with MetS, and increased levels of CER18 and DCER22:2 were especially prominent. On the other hand, most lactosylceramide and hexosylceramide species were lower in subjects with MetS, with LCER14, LCER24:1, and LCER16 being particularly distinctive. In contrast, there were not any significant differences in sphingomyelin species between the two groups. Based on these findings, subsequent analyses focused on ceramide, dihydroceramide, hexosylceramide, and lactosylceramide levels.

FIGURE 2

Figure 2. Volcano plots of sphingolipid species levels in adults with or without metabolic syndrome. Statistical significance was evaluated by Wilcoxon-test and fold change threshold in sphingolipids was set to higher than 1.15. Positive log₂(fold change) indicates increased abundance in MetS adults and light blue circles denote p < 0.05. For these analyses, blood lipid levels were used without log transformation.

Opposite relationships of ceramides and simple β-glycosphingolipids with MetS criteria

To analyze the association of the sphingolipid classes and MetS prevalence, binary logistic regression models were run, adjusted for sex, age, race, and education attainment.

As shown in Figure 3A, dihydroceramides and ceramides were positively associated with abdominal obesity, high triglycerides, high fasting glucose, blood hypertension, and MetS. Only their association with low HDL-c was not statistically significant. In contrast, hexosylceramides and lactosylceramides had a negative association with most of these indicators. In particular, the levels of lactosylceramides revealed a strong negative association with all the MetS components.

FIGURE 3

Figure 3. Association of blood sphingolipid profile and metabolic syndrome. Binary logistic regression models for the association between MetS and its components with blood sphingolipid levels (log-transformed), adjusted for sex, age, race, and education. (A) Regression coefficient matrix of the adjusted logistic regression of MetS components with blood sphingolipid levels. Color label represents the regression coefficient of those associations with p<0.05 (positive in red, and negative in blue). (B) Odds ratios for MetS (with 95%CI) per increase in blood levels of sphingolipid classes, as well as CER/LCER and (DCER+CER)/(HCER+LCER) ratios. (C) Odds ratios for MetS (with 95% CI) per increase in blood levels of sphingolipid species.

Figure 3B shows the likelihood of meeting MetS criteria (with 95% CI) for incremental increases in blood sphingolipid levels. Based on this analysis, higher dihydroceramide levels were significantly associated with a 19% higher prevalence of MetS (OR 1.19, 95% CI 1.13–1.24). Similarly, higher ceramide levels were associated with a 20% higher prevalence of MetS (OR 1.20, 95% CI 1.13–1.28). In contrast, higher hexosylceramide and lactosylceramide levels were associated with 18% (OR 0.82, 95% CI 1.77–1.87) and 35% lower prevalence of MetS, respectively (OR 0.65, 95% CI 0.60–0.70).

Considering that association patterns of ceramides and simple β-glycosphingolipids with MetS appeared to be consistently inverse, we next examined ceramides/simple β-glycosphingolipids ratios. Based on adjusted logistic regression models, both ratios –including 1 or 2 sphingolipid classes either in the numerator or denominator– were significantly associated with a higher prevalence of MetS (Figure 3B). In particular, an integrative ratio including the four classes showed the strongest association with MetS prevalence (OR 1.58, 95% CI 1.49–1.68).

Because the level of each lipid species may differ within a sphingolipid class (30), the probability of meeting MetS criterion based on incremental increases in each sphingolipid species were evaluated. Figure 3C shows the species that revealed a significant association with MetS. In these analyses, all tested dihydroceramide species and seven out of eight ceramide species were significantly associated with a higher prevalence of MetS. Specifically, DCER22, DCER22:2, DCER24:1, CER18, and CER22:1 stood out as being strongly related to MetS. In parallel, five of seven hexosylceramides and six of seven lactosylceramides were significantly associated with a lower prevalence of MetS. The strongest negative associations were found for LCER14 and LCER16. Therefore, most ceramide and simple β-glycosphingolipid species were well represented by their total class.

Opposite relationships of ceramides and simple β-glycosphingolipids with biomarkers of metabolic disorders

The association of blood sphingolipid levels with several metabolic biomarkers was examined with linear regression modeling adjusted for sociodemographic variables. As shown in Figure 4A, and consistent with the logistic regression, dihydroceramide and ceramide classes evinced a significant positive correlation with almost all risk biomarkers of cardiometabolic disease, including obesity, blood pressure, atherogenic lipids, and impaired glucoregulation. In contrast, hexosylceramide and lactosylceramide classes were negatively associated with most of these biomarkers, indicative of poor metabolic health. Specifically, hexosylceramide and lactosylceramide levels exhibited a significant negative association with obesity and impaired glucoregulation (Figure 4A). Moreover, lactosylceramide levels were associated with a healthier lipid profile (i.e., lower triglycerides and oxHDL, and higher HDL-c). To verify that the separate associations were not just connected to the overall MetS condition, the linear regression models were adjusted for MetS prevalence. As shown in Figure 4A (right panel), most of the individual associations remained significant in the MetS-adjusted model.

FIGURE 4

Figure 4. Adjusted linear regression of blood sphingolipid levels and health biomarkers. (A) Regression coefficient matrix of sphingolipid classes on health biormarkers adjusted for control variables (sex, age, race, and education) and MetS prevalence. (B) Regression coefficient matrix of sphingolipid species on health biormarkers adjusted for sex, age, race, and education. Color label represents the regression coefficient of those associations with p < 0.05 (positive in red, and negative in blue).

We next delineated the relationship between specific species and health biomarkers, after adjusting for sociodemographic variables (Figure 4B) and prevalence of MetS (Supplementary Figure 1). As shown in Figure 4B, the four analyzed dihydroceramide species had a significant positive relationship with obesity, blood pressure, atherogenic lipids, and impaired glucoregulation. Within the ceramides, seven of the eight analyzed species were positively associated with most of the CVR biomarkers, except for CER26 that was associated only with atherogenic lipids (LDL-c and triglycerides). With regard to the hexosylceramide species, only HCER22:1 showed a distinctive relationship with the health biomarkers (Figure 4B). While all other species were negatively associated with obesity and impaired glucose metabolism, HCER22:1 was not correlated with healthier glucoregulation. Moreover, HCER22:1 displayed a positive association with high triglyceride levels, while almost all other species were associated with a healthier lipid profile (i.e., lower triglycerides and oxHDL, and higher HDL-c). A similar pattern was evident with lactosylceramides. Six of seven species displayed a significant negative association with obesity, dyslipidemia, and impaired glucoregulation. Conversely, LCER18:1 was positively associated with some biomarkers, including oxHDL, insulin, and HOMA-IR, but not others (Figure 4B). Thus, just some lipid species of a same class behaved in a distinctive manner, and most of them were well represented by total class levels.

In addition, most associations detected for sphingolipid species remained unchanged after adding MetS as a control variable in the MetS-adjusted models (Supplementary Figure 1).

Positive association between sphingolipid levels and inflammatory biomarkers

Because inflammatory processes can also contribute to atherogenesis and ASCVD, the associations between sphingolipid levels and inflammatory mediators and adhesion proteins were assessed. Simple β-glycosphingolipids were negatively associated with inflammatory measures, including IL-6 and CRP, after adjusting for sociodemographic variables, while dihydroceramides and ceramides were positively related with most of these biomarkers (Figure 5A). Considering potential covariation with the overall MetS diagnosis, the analysis was rerun separately for participants in each subgroup. In adults without MetS, the four sphingolipid classes were still positively associated with at least one inflammatory marker (Figure 5B, left panel). Surprisingly, while no ceramide or dihydroceramide was positively related to inflammatory activity in the participants with MetS, blood levels of lactosylceramide were significantly associated with higher levels of the four inflammatory biomarkers (Figure 5B, right panel).

FIGURE 5

Figure 5. Adjusted linear regression of blood sphingolipid levels and inflammatory cardiovascular risk biomarkers. Regression coefficient matrix between sphingolipid classes and inflammatory mediators and adhesion molecules, adjusted for control variables (sex, age, race, and education) in (A) all subjects or (B) separated by absence or presence of MetS diagnosis. Color label represents the regression coefficient of those associations p < 0.05 (positive in red, and negative in blue).

Sphingolipid network differences according to MetS diagnosis and inflammatory status

To evaluate the alterations of sphingolipid metabolic networks related to MetS and inflammatory status, distance networks were constructed and the relationship between sphingolipid species was assessed (Figure 6). For this analysis, subjects were subdivided regarding MetS criteria (i.e., positive or negative) and their CRP levels (i.e., below or above the median), and then the sphingolipid networks of the 4 subgroups were compared.

FIGURE 6

Figure 6. Sphingolipid distance networks based on metabolic syndrome diagnosis and CRP levels. Network analyses were performed using sphingolipids selected in the volcano plot and built based on significant correlations (BH adjusted-Spearman correlation p < 0.05). Nodes indicate sphingolipid species and connecting lines represent sphingolipid distances. Black lines denote intra-cluster relationships and red lines represent inter-cluster relationships. Various clusters are indicated by different colored cloud clusters. (A) Sphingolipid network in adults without MetS and low levels of CRP (n = 770). (B) Sphingolipid network in adults without MetS and high levels of CRP (n = 517). (C) Sphingolipid network in adults with MetS and low levels of CRP (n = 260). (D) Sphingolipid network in adults with MetS and high levels of CRP (n = 505). Statistical comparisons of whole networks and interclasses distances between groups are shown.

The sphingolipid networks of all subgroups displayed several broad similarities. In agreement with the previous results, ceramides and dihydroceramides formed a different cluster that was distinct from simple β-glycosphingolipids (hexosylceramides and lactosylceramides), a finding that was evident in all 4 groups (Figures 6A–D).

Despite these overall similarities, there were some important differences in the networks across the 4 groups (Figure 6). Participants with low CRP who did not meet MetS criterion were considered as the reference healthier group (MetS⁻CRP^low, Figure 6A). When the distances between sphingolipid classes were compared (Figures 6, 7), adults without MetS and high CRP levels (MetS⁻CRP^high, Figure 6B) exhibited a stronger association between precursor ceramides and hexosylceramide levels as well as a lower correlation between hexosylceramide and lactosylceramide levels. Consequently, the distance between ceramides and hexosylceramides in MetS⁻CRP^high subjects was shorter, while the hexosylceramides were separated from the lactosylceramide cluster (Figure 6B). On the other hand, subjects who exhibited MetS but had low CRP levels (MetS⁺CRP^low, Figure 6C) showed a slight strengthening of the association between dihydroceramides and ceramides, without other differences, when compared with MetS⁻CRP^low subjects (Figure 6A). Interestingly, when participants did have elevated CRP levels in the context of MetS (MetS⁺CRP^high), the association between all classes was increased, which manifested by a shorter distance between them (Figure 6D). Postulated changes along the sphingolipid metabolic pathway based on interclass distances and cluster analyses are presented in Figure 7.

FIGURE 7

Figure 7. Postulated metabolic changes based on sphingolipid networks analyses and comparisons in metabolic syndrome. Different colors indicate different clusters. The length of the arrows represents interclasses distances and the letters–with or without asterisk–indicate the statistical differences obtained between the four networks.

Finally, we determined the impact of poor glucoregulation on sphingolipid networks of participants with or without MetS (Supplementary Table 2). Network patterns were significantly different in adults showing signs of impaired glucoregulation. However, no major differences were found in inter-class distances. Only adults with impaired glucose metabolism and MetS showed a stronger association between hexosylceramides and lactosylceramides, when compared to the healthier subgroup without either condition.

Discussion

The current analyses revealed opposite associations between ceramides and simple β-glycosphingolipid with MetS and its related metabolic, inflammatory, and endothelial dysfunction pathophysiological biomarkers. Specifically, high ceramide but low simple β-glycosphingolipid levels were linked with obesity, atherogenic dyslipidemia, and unhealthy glucoregulation. However, in contrast to their inverse relationship with metabolic markers, high levels of lactosylceramides were positively correlated with inflammatory and endothelial CVR markers in adults with MetS. These differential –even divergent– associations indicate the multifactorial aspects of lipidomic profiles and their association with disease risk factors. It also conveys how sphingolipid profiling may provide new insights into the underlying pathways that mediate metabolic vs. cardiovascular risk in MetS.

Unlike most previous lipidomic studies performed in patients with ASCVD or DM (5, 6, 16, 31), our findings allowed a detailed characterization of sphingolipid levels and the potential pathobiology that lies beneath the context of MetS. For example, if sphingolipid synthesis was globally upregulated with increasing lipid storage observed in MetS (32), all sphingolipids should have a positive correlation with obesity related disorders. The distinctive relationships between the precursor ceramides when compared to the down-stream simple β-glycosphingolipids suggests the regulation of sphingolipid metabolism and its contribution to metabolic disorders are more complex than expected. Moreover, within each sphingolipid class, a variety of species can be defined by structural and chemical features that may have different functions or targets. In this study most species behaved similarly within the same class. The species with the highest positive association to MetS were CER18 and DCER22, and those with the highest negative association were LCER14 and LCER16. These results are consistent with previous studies showing similar profiles in patients with insulin resistance or diabetes (30, 33).

Recent research indicates that visceral obesity and chronic inflammation increase dihydroceramide and ceramide levels, which in turn can worsen MetS-associated pathophysiology, including atherothrombotic status and insulin resistance, which ultimately leads to clinical manifestations of ASCVD and/or type 2 diabetes, respectively (32, 34–37). Indeed, ceramides have been identified as candidate biomarkers of ASCVD risk and mortality (9, 16, 38). However, using different methodologies and populations, three recent studies reported negative correlations between hexosylceramide levels and the biochemical and clinical measures of obesity and diabetes (6, 39, 40). Those results are consistent with our findings from a MetS perspective. Furthermore, one study found a negative association between simple β-glycosphingolipids and the diagnosis of type 2 diabetes during follow-up (6).

Therefore, the metabolic pathway of sphingolipid synthesis is not completely up-regulated in MetS or some other compensatory pathway exists. Considering that simple β-glycosphingolipids were decreased in adults with MetS and diabetes, a rise in the salvage pathway leading to the production of ceramides by catabolism of hexosylceramides and lactosylceramides may be enhanced (Figure 1).

Based on this, we propose to include simple glycosphingolipids as the denominator in a ratio between ceramides and β-glycosphingolipids, to be used as a novel metabolic and diabetic risk indicator with a better predictive value rather than the evaluation of ceramides alone (10). The use of a simple ratio of ceramides to β-glycosphingolipids allows normalization of ceramide values by other sphingolipid classes. This standardization would be particularly useful when making comparisons across racial groups, including African Americans, in whom the overall sphingolipid levels are lower than in European American subjects, who are typically used as reference populations (41, 42). Further longitudinal follow-up in MIDUS as well as validation in additional prospective cohorts should establish the actual predictive value of this sphingolipid ratio, for instance, in forecasting the incidence of type 2 diabetes.

However, this ratio of ceramide and simple β-glycosphingolipid levels would only be valid for glucometabolic alterations present in MetS that are associated with progression to diabetes. It may not necessarily apply as predictor of CVR. Indeed, two large cohort-based lipidomic studies positively associated high levels of hexosylceramides and lactosylceramides with future cardiovascular events and cardiovascular death (5, 16). However, the cross-sectional design of our analysis prevented prospective analysis on the association between simple β-glycosphingolipids and cardiovascular outcomes. Given that we have not completed a detailed follow-up on ASCVD incidence in our cohort, we evaluated the association between sphingolipid levels and pathogenically relevant inflammatory and endothelial markers as surrogates of future CVR. Consistent with previous studies (5, 16), a positive relationship between lactosylceramides with IL-6, CRP, fibrinogen, and sICAM was found in subjects with MetS, strongly suggesting that high levels of these β-glycosphingolipids may be associated with increased risk of future cardiovascular events. Indeed, correlation networks analysis showed significant differences in the sphingolipidome of individuals with or without MetS based on their CRP levels, a key biomarker related with initiation, progression, and ischemic complications of atherosclerosis (43).

Based on these findings, new hypotheses can be postulated regarding changes in sphingolipid metabolism (Figure 7) as cause and/or consequence of metabolic and inflammatory alterations in subjects with MetS. The sphingolipid cluster distance analysis suggests that high levels of CRP generate a greater imbalance in sphingolipid metabolism than the presence of MetS per se. High levels of CRP in the absence of MetS seems to stimulate the flux from dihydroceramides to hexosylceramides. This implies a stronger metabolic association between hexosylceramides and their precursors, rather than lactosylceramide levels downstream. In contrast, subjects with MetS and low CRP levels displayed a strengthened association between dihydroceramides and ceramides, without changes in the relationship with simple β-glycosphingolipids. The salvage pathway for ceramide formation may be activated as a compensatory mechanism, as mentioned above. Finally, the presence of MetS and high levels of CRP seems to affect the entire metabolic pathway of sphingolipids, increasing the association between all 4 classes analyzed. This means that high levels of the precursor forms would also lead to increased production of lactosylceramides, unlike what was observed in other subgroups. Therefore, the sphingolipidomic network profile may allow to differentiate subjects not only by MetS diagnosis, but also by their current inflammatory status.

Even though MetS has previously been associated with higher odds of cardiovascular events, it remains unclear whether MetS should be considered a single clinical entity or a group of different clinical conditions with diverse underlying mechanisms and carrying also distinctive cardiovascular vs. metabolic risks (44). According to our findings, more specific scores for cardiovascular vs. metabolic risk, using novel biomarkers that include both ceramides and glycosphingolipids, should allow a better stratification of MetS-associated risks.

The present study has several strengths, including the availability of lipidomic data and many clinical biomarkers in a large number of American adults. However, the cross-sectional design of this study prevents inferring causal relationships between baseline sphingolipid profiling and the progression of MetS. Moreover, cardiovascular events usually result from long-term and sustained metabolic and inflammatory conditions (45). Thus, further evidence from longitudinal studies and independent cohorts is needed to evaluate if different sphingolipid profiles -beyond ceramides- may improve prediction of progression of MetS into diabetes vs. ASCVD. Additional follow-up studies can also reveal variations in blood sphingolipid patterns along time and whether these sequential changes correlate with MetS progression to chronic disease.

In conclusion, our results showed opposite relationships of ceramides and simple β-glycosphingolipids, specially lactosylceramides, with MetS and its glucometabolic related biomarkers. Interestingly, although simple β-glycosphingolipids were associated inversely with metabolic measures, they had a positive association with inflammatory and endothelial CVR markers in subjects with MetS. Thus, our study suggests that a broader analysis of the blood sphingolipid profile, rather than ceramide levels only, may be applied as a more informative MetS-associated biomarker. In addition, it may provide new and relevant insights on the progression of this increasingly and highly prevalent risk condition that enhances the burden of chronic diseases worldwide.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: https://midus.colectica.org/. Lipidomic data was deposited to OSF database under the https://doi.org/10.17605/OSF.IO/VFR7B.

Ethics statement

The studies involving human participants were reviewed and approved by the Health Sciences Institutional Review Board at the University of Wisconsin-Madison, as well as the Institutional Review Boards at the University of California-Los Angeles and Georgetown University. The patients/participants provided their written informed consent to participate in this study.

Author contributions

LB conducted analyses, interpreted results, and drafted the manuscript. CS assisted with data management, analysis, and manuscript revision. CC and CR participated in conception and design of the study, obtained funding, took care of the project administration, and critically reviewed the manuscript. AR was involved in conception, design of the present study, and substantial revision of the work. All authors contributed to the manuscript and approved the submitted version.

Funding

This research was supported by grants from the National Institute on Aging (P01 AG020166 and U19-AG051426) and by the National Fund of Scientific and Technological Development FONDECYT Chile (Postdoc 3210391). Since 1995 the MIDUS study has been funded by the MacArthur Foundation Research Network on Successful Midlife Development.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2022.1092331/full#supplementary-material

References

1. Benjamin EJ, Virani SS, Callaway CW, Chamberlain AM, Chang AR, Cheng S, et al. Heart disease and stroke statistics−2018 update: a report from the American Heart Association. Circulation. (2018) 137:e67–e492. doi: 10.1161/CIR.0000000000000573

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Grundy SM. Metabolic syndrome update. Trends Cardiovasc Med. (2016) 26:364–73. doi: 10.1016/j.tcm.2015.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Alberti KGMM, Eckel RH, Grundy SM, Zimmet PZ, Cleeman JI, Donato KA, et al. Harmonizing the metabolic syndrome: a joint interim statement of the International Diabetes Federation Task Force on Epidemiology and Prevention; National Heart, Lung, and Blood Institute; American Heart Association; World Heart Federation; International Atherosclerosis Society; and International Association for the Study of Obesity. Circulation. (2009) 120:1640–5. doi: 10.1161/CIRCULATIONAHA.109.192644

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Paredes S, Fonseca L, Ribeiro L, Ramos H, Oliveira JC, Palma I. Novel and traditional lipid profiles in Metabolic Syndrome reveal a high atherogenicity. Sci Rep. (2019) 9:11792. doi: 10.1038/s41598-019-48120-5

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Alshehry ZH, Mundra PA, Barlow CK, Mellett NA, Wong G, McConville MJ, et al. Plasma lipidomic profiles improve on traditional risk factors for the prediction of cardiovascular events in type 2 diabetes mellitus. Circulation. (2016) 134:1637–50. doi: 10.1161/CIRCULATIONAHA.116.023233

PubMed Abstract | CrossRef Full Text | Google Scholar

6. 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) 5:e126925. doi: 10.1172/jci.insight.126925

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Bismuth J, Lin P, Yao Q, Chen C. Ceramide: a common pathway for atherosclerosis? Atherosclerosis. (2008) 196:497–504. doi: 10.1016/j.atherosclerosis.2007.09.018

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Chaurasia B, Summers SA. Ceramides – lipotoxic inducers of metabolic disorders. Trends Endocrinol Metabol. (2018) 29:66–7. doi: 10.1016/j.tem.2017.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

9. de Carvalho LP, Tan SH, Ow GS, Tang Z, Ching J, Kovalik JP, et al. Plasma ceramides as prognostic biomarkers and their arterial and myocardial tissue correlates in acute myocardial infarction. JACC Basic Transl Sci. (2018) 3:163–75. doi: 10.1016/j.jacbts.2017.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Hilvo M, Vasile VC, Donato LJ, Hurme R, Laaksonen R. Ceramides and ceramide scores: clinical applications for cardiometabolic risk stratification. Front Endocrinol. (2020) 11:570628. doi: 10.3389/fendo.2020.570628

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Sokolowska E, Blachnio-Zabielska A. The role of ceramides in insulin resistance. Front Endocrinol. (2019) 10:577. doi: 10.3389/fendo.2019.00577

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yaribeygi H, Bo S, Ruscica M, Sahebkar A. Ceramides and diabetes mellitus: an update on the potential molecular relationships. Diabetic Med. (2020) 37:11–9. doi: 10.1111/dme.13943

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Lalazar G, Zigmond E, Weksler-Zangen S, Ya'acov A, Levy MS, Hemed N, et al. Oral Administration of β-glucosylceramide for the treatment of insulin resistance and nonalcoholic steatohepatitis: results of a double-blind, placebo-controlled trial. J Med Food. (2017) 20:458–64. doi: 10.1089/jmf.2016.3753

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Margalit M, Shalev Z, Pappo O, Sklair-Levy M, Alper R, Gomori M, et al. Glucocerebroside ameliorates the metabolic syndrome in OB/OB mice. J Pharmacol Expe Therap. (2006) 319:105–10. doi: 10.1124/jpet.106.104950

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zigmond E, Zangen SW, Pappo O, Sklair-Levy M, Lalazar G, Zolotaryova L, et al. β-Glycosphingolipids improve glucose intolerance and hepatic steatosis of the Cohen diabetic rat. Am J Physiol Endocrinol Metabol. (2009) 296:E72–8. doi: 10.1152/ajpendo.90634.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Poss AM, Maschek JA, Cox JE, Hauner BJ, Hopkins PN, Hunt SC, et al. Machine learning reveals serum sphingolipids as cholesterol-independent biomarkers of coronary artery disease. J Clin Invest. (2020) 130:1363–76. doi: 10.1172/JCI131838

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Berkowitz L, Cabrera-Reyes F, Salazar C, Ryff CD, Coe C, Rigotti A. Sphingolipid profiling: a promising tool for stratifying the metabolic syndrome-associated risk. Front Cardiovasc Med. (2021) 8:785124. doi: 10.3389/fcvm.2021.785124

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ryff C, Almeida DM, Ayanian J, Carr DS, Cleary PD, Coe C. Midlife in the United States (MIDUS 2), 2004–2006. Inter-university Consortium for Political Social Research (2017). Available online at: https://www.icpsr.umich.edu/web/NACDA/studies/4652

Google Scholar

19. Radler BT, Ryff CD. Who Participates? Accounting for longitudinal retention in the MIDUS national study of health and wellbeing. J Aging Health. (2010) 22:307–31. doi: 10.1177/0898264309358617

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Kirsch JA, Ryff CD. Hardships of the Great Recession and health: understanding varieties of vulnerability. Health Psychol Open. (2016) 3:2055102916652390. doi: 10.1177/2055102916652390

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Amaral WZ, Krueger RF, Ryff CD, Coe CL. Genetic and environmental determinants of population variation in interleukin-6, its soluble receptor and C-reactive protein: Insights from identical and fraternal twins. Brain Behav Immun. (2015) 49:171–81. doi: 10.1016/j.bbi.2015.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Love G, Seeman T, Weinstein M, Ryff C. Bioindicators in the MIDUS national study: protocol, measures, sample, and comparative context. J Aging Health. (2010) 22:1059–80. doi: 10.1177/0898264310374355

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Flaherty SM, Wood EK, Ryff CD, Love GD, Kelesidis T, Berkowitz L, et al. Race and sex differences in HDL peroxide content among American adults with and without type 2 diabetes. Lipids Health Dis. (2022) 21:18. doi: 10.1186/s12944-021-01608-4

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Berkowitz L, Henríquez MP, Salazar C, Rojas E, Echeverría G, Love GD, et al. Association between serum sphingolipids and eudaimonic well-being in white U.S. adults. Sci Rep. (2021) 11:13139. doi: 10.1038/s41598-021-92576-3

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Libby P. Inflammation in atherosclerosis. Arterioscler Thromb Vasc Biol. (2012) 32:2045–51. doi: 10.1161/ATVBAHA.108.179705

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Clauset A, Newman MEJ, Moore C. Finding community structure in very large networks. Phys Rev E. (2004) 70:6. doi: 10.1103/PhysRevE.70.066111

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Csardi G, Nepusz T. The igraph software package for complex network research. Int J Complex Syst. (2006) 1695:1–9. Available online at: http://igraph.sf.net

PubMed Abstract | Google Scholar

28. Steiger JH. Testing pattern hypotheses on correlation matrices: alternative statistics and some empirical results. Multivariate Behav Res. (1980) 15:335–52. doi: 10.1207/s15327906mbr1503_7

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Revelle W. psych: Procedures for Psychological, Psychometric, Personality Research. R Package Version (2021). Available online at: https://cran.r-project.org/web/packages/psych/psych.pdf

Google Scholar

30. Kurz J, Parnham MJ, Geisslinger G, Schiffmann S. Ceramides as novel disease biomarkers. Trends Mol Med. (2019) 25:20–32. doi: 10.1016/j.molmed.2018.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Laaksonen R, Ekroos K, Sysi-Aho M, Hilvo M, Vihervaara T, Kauhanen D, et al. Plasma ceramides predict cardiovascular death in patients with stable coronary artery disease and acute coronary syndromes beyond LDL-cholesterol. Eur Heart J. (2016) 37:1967–76. doi: 10.1093/eurheartj/ehw148

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Castillo RI, Rojo LE, Henriquez-Henriquez M, Silva H, Maturana A, Villar MJ, et al. From molecules to the clinic: Linking schizophrenia and metabolic syndrome through sphingolipids metabolism. Fronti Neurosci. (2016) 10:488. doi: 10.3389/fnins.2016.00488

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Jensen PN, Fretts AM, Yu C, Hoofnagle AN, Umans JG, Howard B, et al. Circulating sphingolipids, fasting glucose, and impaired fasting glucose: the strong heart family study. EBioMedicine. (2019) 41:44–9. doi: 10.1016/j.ebiom.2018.12.046

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Chavez JA, Summers SA. A ceramide-centric view of insulin resistance. Cell Metab. (2012) 15:585–94. doi: 10.1016/j.cmet.2012.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Cowart LA. Sphingolipids: players in the pathology of metabolic disease. Trends Endocrinol Metab. (2009) 20:34–42. doi: 10.1016/j.tem.2008.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Luukkonen PK, Zhou Y, Sädevirta S, Leivonen M, Arola J, Orešič M, et al. Hepatic ceramides dissociate steatosis and insulin resistance in patients with non-alcoholic fatty liver disease. J Hepatol. (2016) 64:1167–75. doi: 10.1016/j.jhep.2016.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Wigger L, Cruciani-Guglielmacci C, Nicolas A, Denom J, Fernandez N, Fumeron F, et al. Plasma dihydroceramides are diabetes susceptibility biomarker candidates in mice and humans. Cell Rep. (2017) 18:2269–79. doi: 10.1016/j.celrep.2017.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Hilvo M, Meikle PJ, Pedersen ER, Tell GS, Dhar I, Brenner H, et al. Development and validation of a ceramide- and phospholipid-based cardiovascular risk estimation score for coronary artery disease patients. Eur Heart J. (2019) 41:371–80. doi: 10.1093/eurheartj/ehz387

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Huynh K, Barlow CK, Jayawardana KS, Weir JM, Mellett NA, Cinel M, et al. High-throughput plasma lipidomics: detailed mapping of the associations with cardiometabolic risk factors. Cell Chem Biol. (2019) 26:71–84.e4. doi: 10.1016/j.chembiol.2018.10.008

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Lemaitre RN, Yu C, Hoofnagle A, Hari N, Jensen PN, Fretts AM, et al. Circulating sphingolipids, insulin, HOMA-IR, and HOMA-B: the strong heart family study. Diabetes. (2018) 67:1663–72. doi: 10.2337/db17-1449

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Buie JNJ, Hammad SM, Nietert PJ, Magwood G, Adams RJ, Bonilha L, et al. Differences in plasma levels of long chain and very long chain ceramides between African Americans and whites: An observational study. PLoS ONE. (2019) 14:e0216213. doi: 10.1371/journal.pone.0216213

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Hu J, Yao J, Deng S, Balasubramanian R, Jiménez MC, Li J, et al. Differences in metabolomic profiles between black and white women and risk of coronary heart disease: an observational study of women from four US cohorts. Circ Res. (2022) 131:601–15. doi: 10.1161/CIRCRESAHA.121.320134

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Amar J, Fauvel J, Drouet L, Ruidavets JB, Perret B, Chamontin B, et al. Interleukin 6 is associated with subclinical atherosclerosis: a link with soluble intercellular adhesion molecule 1. J Hypertens. (2006) 24:1083–8. doi: 10.1097/01.hjh.0000226198.44181.0c

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Scuteri A, Laurent S, Cucca F, Cockcroft J, Cunha PG, Mañas LR, et al. Metabolic syndrome across Europe: different clusters of risk factors. Eur J Prev Cardiol. (2015) 22:486–91. doi: 10.1177/2047487314525529

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Nilsson J, Gonçalves I, Edsfeldt A. Chronic inflammation and atherosclerosis. Early Vasc Aging. (2015) 157–67. doi: 10.1016/B978-0-12-801387-8.00016-8

CrossRef Full Text | Google Scholar