Overall descriptive characteristics of the study sample
The study sample included 967 middle-aged and older adults; over half were women (54.4%) and had some college education (Table 1). More than 50% of the participants were overweight or obese with a median BMI of 28.1 kg/m2. With regard to medication use, 34.6% of the participants took anti-hypertensive drugs, 29.6% hypolipidemic agents, 4.2% corticosteroids, and 15.8% were using anti-depressants.
Table 1
Descriptive statistics for study variables in participants from MIDUS (n = 967)
Variables
|
%
|
Median
|
Range
|
Age, years
|
|
55
|
34–84
|
Sex, % males
|
45.6
|
|
|
Education, score points
|
|
8
|
1–12
|
Anti-hypertensive meds, % yes
|
34.6
|
|
|
Cholesterol lowering meds, % yes
|
29.6
|
|
|
Corticosteroids, % yes
|
4.2
|
|
|
Anti-depressants, % yes
|
15.8
|
|
|
BMI, kg/m2
|
|
28.1
|
15.0-63.2
|
As an initial approach to visualize the relationship between well-being dimensions and blood sphingolipid classes, a hierarchical log-heatmap clustering analysis was performed (Fig. 2). Overall, the heatmap analysis showed that subjects with higher levels of well-being had lower serum levels of most sphingolipids and vice versa (Fig. 2).
Inverse associations between eudaimonic well-being dimensions and sphingolipid classes
To further dissect the relationship between serum levels of sphingolipids and eudaimonic well-being dimensions, bivariate analyses were conducted for each of the major variables. As shown in Model 1 (Table 2), higher levels of environmental mastery, purpose in life, and self-acceptance were significantly correlated with lower levels of total ceramides, dihydroceramides, and hexosylceramides. Significant associations were not found for other components of well-being or for other classes of sphingolipids (Table S1).
Table 2
Association between blood levels of sphingolipid classes and well-being dimension scores.
|
Dihydroceramides
|
Ceramides
|
Hexosylceramides
|
|
Coef.
|
p-value
|
Coef.
|
p-value
|
Coef.
|
p-value
|
Model 1: Bivariate analysis
|
Environmental mastery
|
-0.088
|
0.006
|
-0.122
|
0.000
|
-0.080
|
0.013
|
Purpose in life
|
-0.067
|
0.038
|
-0.101
|
0.002
|
-0.088
|
0.006
|
Self-acceptance
|
-0.066
|
0.040
|
-0.106
|
0.001
|
-0.081
|
0.011
|
Model 2: Multivariate analysis adjusting for demographic variables, medication usage and depression.
|
Environmental mastery
|
-0.728
|
0.059
|
-0.118
|
0.002
|
-0.051
|
0,152
|
Purpose in life
|
-0.023
|
0.521
|
-0.531
|
0.126
|
-0.059
|
0.071
|
Self-acceptance
|
-0.017
|
0.657
|
-0.066
|
0.079
|
-0.056
|
0.109
|
* p-values < 0.05 and their standardized regression coefficients are denoted in bold.
|
The summary statistics for Model 2 (Table 2 and S1) indicate that the significance of some aspects of the associations found between sphingolipids levels and well-being dimensions were modified after adjusting for demographic factors, medication usage, and depression, which should be considered potential covariates. As conveyed by this model, associations between dihydroceramides and hexosylceramides with environmental mastery, purpose in life, and self-acceptance were fully explained by the covariates, particularly by sex and cholesterol medications (Table S2), suggesting that these factors account for the bivariate association (Table 2, Model 2). However, the relationship between ceramides and environmental mastery remained significant even when considering these covariates, including after adjusting for depressive symptoms (assessed by CES-D Scale) (Table 2, Model 2). Therefore, the link between ceramides and environmental mastery did not appear to be dependent on sociodemographic variables, medication use, or statistical biases attributable to negative affect.
Association between ceramides and environmental mastery: effect of health behaviors
Hierarchical linear regression models were used to further examine the associations between environmental mastery and blood ceramides, and to discern whether health behaviors might be potential mediators of the relationships between PWB and these novel lipids.
Models 3, 4, 5 and 6 examined the role of each health behavior as a possible mediator of the negative association between a sense of more environmental mastery and ceramides (Table 3). Only BMI -as a distal integrative indicator of overall lifestyle- and sleep quality were found to be significantly correlated with ceramide levels. Moreover, the correlation between environmental mastery and ceramide levels remained significant even after adjustments that included the four separate measures of health behavior, suggesting an association that was at least partially independent of these factors (Table 3).
Table 3
Multivariate association analysis between environmental mastery and ceramides considering the effect of health behaviors.
|
Ceramides vs
|
|
Coef.
|
p-value
|
Coef.
|
p-value
|
Coef.
|
p-value
|
Coef.
|
p-value
|
Environmental mastery
|
-0.115
|
0.003
|
-0.112
|
0.003
|
-0.116
|
0.002
|
-0.118
|
0.002
|
Diet
|
-0.041
|
0.221
|
-0.035
|
0.295
|
-0.030
|
0.366
|
-0.032
|
0.353
|
Exercise
|
|
|
-0.035
|
0.285
|
-0.020
|
0.549
|
-0.019
|
0.573
|
BMI
|
|
|
|
|
0.113
|
0.001
|
0.116
|
0.001
|
Sleep quality
|
|
|
|
|
|
|
0.095
|
0.011
|
* p-values < 0.05 and their standardized regression coefficients are denoted in bold. |
Models 3, 4, 5 and 6 were also analyzed for other dimensions of well-being and SPL classes. However, no statistically significant associations were found (Table S1). These follow-up analyses also included a seventh model that adjusted for CVD, which could potentially influence the SPL levels. Nevertheless, observed relationships did not lose significance after this additional measure was considered (Table S1).
Subsequently, mediation bootstrapping analyses were conducted to statistically confirm the mediating effect of BMI and sleep quality on the relationship between environmental mastery and total ceramides. As observed in the mediation model (Fig. 3) and concurring with the hierarchical analysis (Table 3), both BMI and quality of sleep significantly affected blood ceramide levels. However, based on the bootstrapping method, only sleep quality was a statistically significant mediator of the relationship between environmental mastery and total ceramides, when considering both parameters as independent mediators (Table 4). Nevertheless, when the joint effect of both mediators was considered as representative overall health behaviors, the mediation was statistically significant (Table 4).
Table 4
Bootstrapping to test indirect and total effects of environmental mastery on blood ceramide levels that were mediated by sleep quality (PSQI) and BMI, separately or combined.
Effects
|
Coef, ± SE
|
95% CI
|
p-value
|
Indirect effect
|
|
|
Sleep (PSQI)
|
-0.024 ± 0.011
|
-0.046, -0.004
|
0.022
|
BMI
|
-0.005 ± 0.005
|
-0.016, 0.003
|
ns
|
Both
|
-0.030 ± 0.011
|
-0.053, -0.008
|
0.008
|
Total effect
|
|
|
|
Total
|
-0.152 ± 0.038
|
-0.229, -0.078
|
< 0.001
|
Association between environmental mastery and sphingolipid species
Within each sphingolipid class, there are many component species defined by structural and chemical features (i.e., carbon chain length, double bonds) of the attached fatty acids (Fig. 1A) that correlate with different functions and targets17. Because the levels of each species can also differ17, the same hierarchical models were applied to more specifically delineate the relationship between each molecular species of sphingolipids and the different dimensions of eudaimonic well-being (Table 5).
Table 5. Multivariate association analysis of the relationship between psychological well-being domains and individual species of sphingolipids.

*Model 1: bivariate analysis. Model 2: adjusted for demographic factors, medication usage, and depression. Models 3, 4, 5 and 6 determined the influence of each health behavior as a possible mediators, adding them sequentially: food intake quality (MIDUS-HEI) in Model 3, exercise (METs) in Model 4, sleep quality (PSQI) in Model 5, and Body Mass Index (BMI, kg/m2) in Model 6. Model 7 includes adjustments for CVD and diabetes.
*Gray boxes represent statistically significant associations (all with p<0.05 and regression coefficients <0).
Based on the hierarchical analysis shown for Model 1, all three dihydroceramides included in this analysis were negatively correlated with environmental mastery. DCER24 and DCER24:1 remained significantly associated with environmental mastery after adjustments by the covariates in Model 2. Models 3, 4, 5, and 6 indicated these relationships were partially attenuated by BMI (based on β reduction and p < 0.05 for BMI), but remained significant after full adjustments (Table 5).
Similarly, six out of seven species of ceramides exhibited negative and statistically significant correlations with environmental mastery (Table 5). Considering these six ceramide species, only the correlation observed with CER16 lost significance after being adjusted for the covariates in Model 2. Interestingly, Models 3, 4, 5, and 6 indicated that all health behaviors were significantly associated with ceramides at the species level of analysis (p < 0.05 for each health behavior). Nonetheless, the fully adjusted model demonstrated that the association between these ceramide species and environmental mastery remained significant (Table 5).
Some hexosylceramide, lactosylceramide and sphingomyelin species showed significant inverse correlations with environmental mastery (Table 5). However, only the associations between HCER20, LCER18:1 and SM20:1 with environmental mastery remained significant after adjustments for covariates (Table 5). The health behaviors included in Models 3, 4, 5, and 6 did not appear to exert any significant influence on these latter associations (Table 5).
Association between purpose in life and self-acceptance with sphingolipid species
Six species of ceramides (CER16, 18, 20, 22, 24, and 24:1) were negatively correlated with purpose in life and self-acceptance (Table 5). Only three (CER18, CER20, and CER22) were significantly associated with purpose in life after the covariate adjustments. In contrast, 4 ceramides (CER16, 18, 20, and 24:1) evidenced a significant relationship with self-acceptance after adjustments (Model 2). In all cases, the associations between these two PWB domains were attenuated to non-significant levels when looking at the species level of ceramides after adjusting for health behaviors (Models 3, 4, 5, and 6). This outcome suggested that the observed correlations were more appropriately explained by health behaviors. Similarly, most of the dihydroceramide and hexosylceramide species appeared initially to have a negative relationship with purpose in life and self-acceptance. However, these correlations no longer attained statistical significance after the covariate adjustments (Table 5). The only exception was HCER22, which remained significantly associated with purpose in life in the hierarchical models after inclusion of health behaviors. Only a few individual species of the sphingomyelins and lactosylceramides had significant associations with these two PWB domains, and all became nonsignificant after adjustments (Table 5). The same models were performed for other PWB domains (autonomy, personal growth and positive relations). However, no species showed significant associations after adjusting for control variables.