ORIGINAL RESEARCH article

Front. Immunol., 19 November 2025

Sec. Molecular Innate Immunity

Volume 16 - 2025 | https://doi.org/10.3389/fimmu.2025.1659425

Plasma lipidomic signatures reveal age-associated patterns of septic shock risk and immune dysregulation in sepsis

  • 1. Department of Infectious Diseases, Shanghai Key Laboratory of Infectious Diseases and Biosafety Emergency Response, National Medical Center for Infectious Diseases, Huashan Hospital, Fudan University, Shanghai, China

  • 2. Shanghai Omicsolution Co., Ltd., Shanghai, China

  • 3. Shanghai Sci-Tech Inno Center for Infection & Immunity, Shanghai, China

Abstract

Background:

Sepsis is associated with significant lipidomic disturbances, but age-associated lipidomic patterns remain poorly characterized. Given the links between aging, immune dysfunction, and metabolic dysregulation, defining age-specific lipid profiles could improve sepsis risk assessment. This study investigates age-stratified lipidomic signatures in sepsis and identifies biomarkers for clinical severity.

​Methods:

We prospectively enrolled 62 sepsis patients (21 <65 years, 41 ≥65 years) and 40 healthy controls. Plasma lipidomics was performed via untargeted LC-MS/MS, identifying 1,277 lipid species across 38 subclasses. Principal component analysis (PCA) and consensus clustering were used to assess lipidomic differences and define patient subtypes. Correlations between lipid subclasses, clinical severity (SOFA scores), and immune cell subsets were analyzed. An age-adjusted risk stratification model was developed to assess septic shock and mortality risk (AUC analysis).

​Results:

Sepsis patients exhibited reduced phosphatidylcholine (PC), cholesteryl ester (CE), and lysophosphatidylcholine (LPC) levels (all P<0.05). Clustering revealed four lipidomic patterns, with Cluster 4 distinguishing two sepsis subtypes (C1, C2). Subtype C2 had higher septic shock incidence (57.1% vs. 14.8%, P = 0.0013) and downregulation of 92 lipids, 35 of which strongly correlated with SOFA scores. A risk stratification model incorporating six key lipids (LPC(19:0), PC(P-19:0), SM 32:3;2O(FA 16:3), PC(P-20:0), PC(O-18:1/20:3), CE(15:0)) and age accurately predicted septic shock (AUC: 0.87 training, 0.82 validation) and mortality risk in elderly patients. PC levels correlated with monocytes, while CE and LPC associated with complement proteins and CD8+ T cells.

​Conclusions:

Our lipid-based model effectively predicts septic shock and mortality, particularly in elderly sepsis patients. Age-associated lipid alterations (PC, LPC, CE reduction) correlate with disease severity and immune dysregulation, suggesting distinct lipid-immune mechanisms in younger vs. elderly patients. These findings support lipidomics as a tool for sepsis risk stratification and personalized therapy.

Introduction

Sepsis is a complex, life-threatening condition characterized by a dysregulated host response to infection, leading to systemic inflammation and potential multiorgan failure (–). It has emerged as a significant global health challenge, with high incidence and mortality rates, particularly among older adults (, ). Plasma lipidomics, which involves comprehensive analysis of lipid molecules in biological samples, holds great potential for elucidating the metabolic and immune regulatory mechanisms underlying the onset and progression of sepsis (, ), thereby reflecting the host’s metabolic status and immune function (, ). However, current research on the lipidomic profiles of sepsis patients remains insufficient, particularly due to the lack of systematic age-stratified analyses, which hinders the identification of age-specific metabolic and immune alterations in these patients (, ).

Previous studies have focused primarily on targeted measurements of a select few lipid species (, ), failing to capture the broader lipidomic landscape. Additionally, while the literature indicates metabolic disturbances in critically ill patients (, ), there is a lack of age-stratified investigations into the overall relationship of the plasma lipidome with disease severity and immune cell dynamics (). Elderly patients, due to physiological decline, immunosenescence, and metabolic changes, typically exhibit heightened susceptibility to infections (). The age-associated deterioration of immune function further complicates disease management and leads to significantly worse clinical outcomes (). Therefore, conducting age-stratified investigations into the lipid metabolic characteristics of sepsis patients is of paramount importance for elucidating disease mechanisms and optimizing clinical management strategies ().

The interplay between lipid metabolism and the immune response is critical to the progression and outcomes of sepsis (). We hypothesized that there are significant age-specific lipid profile changes in patients with sepsis, and that these changes correlate with disease severity and immune cell function. Lipids serve not only as fundamental structural components of cellular membranes but also as critical modulators of inflammatory and immune signaling pathways (, ). Age-specific alterations in lipid metabolic profiles may be intrinsically linked to divergent immune regulatory mechanisms across age groups. Elucidating these changes could yield clinically actionable biomarkers for improved diagnosis and prognostic prediction, thereby advancing the development of personalized therapeutic strategies (). To validate this hypothesis, we performed comprehensive untargeted liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based plasma lipidomic profiling of both sepsis patients and healthy controls, with particular focus on identifying age-stratified lipidomic variations and their clinical correlations.

Methods

Study population and group allocation

This study was conducted in the Department of Infectious Diseases at Huashan Hospital, Fudan University, from September 2022 to September 2023. The research protocol was approved by the Ethics Review Committee of Huashan Hospital. Inpatients with suspected infections were enrolled, and samples were collected on the day of fever onset. Patients meeting the Sepsis-3 criteria were included in the sepsis group. Patients aged 65 years or older were classified into the elderly sepsis group (asepsis), whereas those younger than 65 years were classified into the younger sepsis group (ysepsis). Additionally, we enrolled 20 elderly healthy controls (aHC) and 20 young healthy controls (yHC) from the health examination center.

Among all sepsis patients, those who exhibited persistent hypotension despite adequate fluid resuscitation and required vasopressor support to maintain a mean arterial pressure (MAP) of ≥ 65 mmHg, along with serum lactate levels > 2 mmol/L, were classified into the septic shock group. Healthy controls were sampled from community health screenings. Informed consent was obtained from all participants.

All blood samples were collected at hospital admission, prior to the administration of any antibiotic therapy. Demographic, clinical, microbiological, and laboratory data were collected by dedicated research physicians. Physiological measurements, vital signs, and disease severity scores, including the Sequential Organ Failure Assessment (SOFA) score, were recorded.

Sample preparation

The samples were thawed at 4 °C, and 100 µL of serum was transferred to a new tube. To this mixture, 200 µL of methanol (MeOH) and 200 µL of acetonitrile (ACN) were added. The mixture was vortexed for 30 seconds and then sonicated in an ice-water bath for 10 minutes to ensure thorough mixing. To precipitate the proteins, the samples were incubated for 1 hour at -20°C. Following incubation, the samples were centrifuged at 13,000 rpm for 15 minutes at 4°C, and the supernatant was carefully collected. The supernatant was placed in a rotary evaporator at 4°C and evaporated to dryness. The dry extracts were then resuspended in 100 µL of a 1:1 mixture of acetonitrile and water. Vortexed for 30 seconds and then sonicated in an ice-water bath for an additional 10 minutes. The resulting mixture was centrifuged at 13,000 rpm for 15 minutes at 4°C, and the supernatant was collected and stored at -80°C until further analysis.

LC–MS/MS analysis

LC-MS/MS analysis was performed via a high-performance liquid chromatography (HPLC) system (1290 Series, Agilent Technologies, USA) coupled to a Bruker timsTOF mass spectrometer equipped with an electrospray ionization (ESI) source (Bruker Daltonik, Germany). A Phenomenex Kinetex C18 column (particle size: 1.7 µm; dimensions: 100 mm length × 2.1 mm inner diameter) was utilized for liquid chromatography separation, with the column temperature maintained at 55°C. Mobile phase A consisted of water/acetonitrile (6:4, v/v) with 10 mM ammonium formate, while mobile phase B consisted of isopropanol/acetonitrile (9:1, v/v) with 10 mM ammonium formate. Flow rate: 0.30 mL/min. Injection volume: 2 µL. Run time: 18 minutes.

The elution gradient was programmed as follows: 40% B from 0 to 1.5 minutes, a linear increase from 40% B to 85% B from 1.5 to 10.5 minutes, holding at 85% B from 10.5 to 14.0 minutes, transitioning from 85% B to 100% B from 14.0 to 14.1 minutes, holding at 100% B from 14.1 to 15.0 minutes, decreasing back to 40% B within 0.2 minutes, and maintaining at 40% B for 2.8 minutes for column equilibration. The flow rate was set at 0.3 mL/min, and a 2 µL injection volume was used for biological samples in positive mode.

For mass spectrometry acquisition, the PASEF-DDA scan mode was employed to detect ions with m/z values ranging from 100 Da to 1350 Da and mobility values ranging from 0.55 to 1.90 V·s/cm². The MS parameters were as follows: capillary voltage was set to +4500 V; nebulizer pressure was 2.2 bar; dry gas flow was 10.0 L/min at a temperature of 220 °C; the number of PASEF MS/MS scans was set to 2; the ramping time was 100 ms; the total cycle time was 0.32 s; the charge range was 0–1; the absolute threshold was 100 counts; active exclusion was enabled, with previously targeted ions released after 0.1 min; the isolation window was set to 1.2 Da; and the collision energy was set to 40 eV.

Lipidomic profiling

Lipid molecular annotation and raw data preprocessing of the mass spectrometry data were performed using the online tool Lipid4DAnalyzer (http://lipid4danalyzer.zhulab.cn/). The final lipid metabolite list and data matrix were obtained. To ensure the reliability of subsequent statistical analyses, the raw data were preprocessed as follows before formal analysis: Feature peaks with a non-zero value occurrence rate of >50% across all samples were retained. Missing values in the remaining feature peaks were imputed using the K-Nearest Neighbors (KNN) method. Normalization was performed using quality control (QC) samples, and feature peaks with a relative standard deviation (RSD) >30% in QC samples were excluded. Through this workflow, a total of 1277 lipid molecules were annotated. Based on the LIPID MAPS (Lipid Metabolites and Pathways Strategy) classification system, these 1277 lipid molecules were categorized into 38 lipid subclasses according to their structures.

Statistical analysis

All statistical analyses were performed using R statistical software (version 4.3.2), with data visualization primarily implemented via the ggplot2 package (version 3.5.1).

We conducted PCA using the FactoMineR package (version 2.11) to reduce dimensionality. The Wilcoxon rank-sum test was applied to assess the significance of inter-group distribution differences along PCA principal coordinates, thereby elucidating global lipidomic profile disparities.

Individual lipid species were compared between groups using two-sample t-tests, with statistical significance defined as P < 0.05. To address multiple testing, all P-values were adjusted via the Benjamini-Hochberg false discovery rate (FDR) method.

Linear regression models (ordinary least squares, OLS) were constructed to evaluate relationships between lipid levels and SOFA scores at two levels:

Lipid subclass level: Expression values of all lipids within a subclass were summed, followed by log2-transformation and z-score normalization to mitigate distributional bias.

Individual lipid level: Raw expression values were directly log2-transformed and z-score normalized.

Normalized lipid subclass/individual lipid values served as independent variables, with SOFA scores as the dependent variable. Regression coefficients and their significance (P-values) were computed.

​Lasso regression​ was employed to identify SOFA-associated key lipids.A ​logistic regression model​ integrating selected lipids and patient age was developed.Model performance was evaluated via ​5-fold cross-validation​ (200 repeats) and quantified using receiver operating characteristic (ROC) curves and area under the curve (AUC).

Results

Study population

We enrolled 62 hospitalized patients diagnosed with sepsis and 40 healthy individuals as controls. Among the sepsis patients, 41 were older adults (≥65 years), and 21 were younger adults (<65 years). The control group included 20 older adults and 20 younger adults. (Figure 1) The mean age of young sepsis was 40 (28-46), with a male prevalence of 48%. The average age of older sepsis cases was 77 (69-82), of which 63% were male. 44% of older sepsis cases had hypertension and 49% had diabetes. In young sepsis, 19% had hypertension and 33% had diabetes. In the sepsis cohort, the most common pathogen was Klebsiella pneumoniae (n=25), followed by Acinetobacter baumannii (n=14), Staphylococcus aureus (n=8), and Enterococcus spp. (n=7). Seventeen patients had mixed infections. The incidence of shock in sepsis patients was 38.7% (24/62), and 21 patients had in-hospital mortality (Table 1).Patients with malignancies or autoimmune diseases were excluded.

Figure 1

Table 1

Patient characteristicsThe young-adult cohort (<65y)The elder cohort(≥65y)
ysepsisyHCasepsisaHC
N21204120
Gender, male (%)10 (48)4 (20)26 (63)10 (50)
Age, median(IQR)40 (28-46)44 (32-57)77 (69-82)75 (70-79)
Comorbidities
 Obesity BMI>28,N(%)2 (10)0 (0)0 (0)0 (0)
 Hypertension,N(%)4 (19)0 (0)18 (44)0 (0)
 Diabetes,N(%)7 (33)0 (0)20 (49)0 (0)
Site of infection
 Lung,N(%)12 (57)0 (0)35 (85)0 (0)
Pathogenic bacteria
 K. pneumoniae,N(%)5 (24)0 (0)12 (29)0 (0)
 A. baumannii,N(%)0 (0)0 (0)6 (15)0 (0)
 S. aureus,N(%)2 (10)0 (0)3 (7)0 (0)
 Enterococcus,N(%)3 (14)0 (0)1 (2)0 (0)
 Multipathogen,N(%)5 (24)0 (0)12 (29)0 (0)
Shock,N(%)8 (36)0 (0)16 (40)0 (0)
Death,N(%)4 (19)0 (0)17 (41)0 (0)

Study population characteristics.

The higher proportion of males in the elderly sepsis group (63% vs. 48% in young sepsis) may reflect both the demographic characteristics of severe sepsis in aging populations and the potentially higher susceptibility of aged males to severe infections.

Lipid landscape differences between aHC and yHC

We first examined the lipidome of healthy baseline samples, defined as samples from individuals without any self-reported acute disease, although potentially including asymptomatic chronic conditions such as prediabetes or undiagnosed diseases. We analyzed 40 healthy baseline samples and conducted t-distributed stochastic neighbor embedding (t-SNE) analysis to characterize lipid abundance differences between younger and older healthy adults. The t-SNE visualization showed substantial overlap in lipidomic profiles between aHC and yHC groups (t-SNE1 P = 0.66), suggesting relative stability of the plasma lipidome in healthy aging (Figure 2a). Box plot analysis revealed no significant difference in overall lipid abundance between the older (aHC) and younger (yHC) healthy groups (P = 0.64) (Figure 2b), indicating similar lipid distributions across ages under healthy conditions.

Figure 2

To facilitate interpretation, we classified the 1,277 identified lipid molecules into 38 subcategories on the basis of structure, including PEtOH, ST, ADGGA, AHexCer, CE, Cer, cholesterol, CoQ, DG, DGCC, DGDG, DGGA, FA, GM, Hex3Cer, HexCer, LacCer, LNAPE, LNAPS, LPC, LPE, LPI, MG, NAE, NAGly, NAGlySer, PC, PE, PE-Cer, PG, PI, PMeOH, SE, SL, SM, Sph, TG, and Tocopherol. Assuming that lipid levels might vary by age, we applied a linear mixed model and found a significant negative correlation with age for most lipid categories, particularly PEtOH (coefficients=-0.029, P = 0.026) and ST (coefficients=-0.028, P = 0.028). MG and GM levels were significantly positively correlated with age (MG, coefficient=0.052, P<0.001; GM, coefficient=0.037, P = 0.004) (Figure 2c). The 38 lipid subcategories ranked similarly in terms of abundance between the aHC and yHC, with no major differences in lipid types between the age groups (Figures 2d, e). FA, MG, and GM levels were higher in older healthy adults than in younger adults (P = 0.018, P<0.001, P = 0.045). PC was the most abundant subcategory, comprising 311 lipid species, whereas CE, Sph, and FA had broader dynamic abundance ranges.

Sepsis-associated lipid signatures

We used principal component analysis (PCA) to explore plasma lipidomic differences between sepsis patients and healthy controls. The PCA results indicated clear separation between groups (P = 3.1e−16) (Figure 3a). A total of 603 differentially expressed lipid molecules (|FC| > 1.3, P < 0.05) were identified between the sepsis and healthy control groups, with 35% of lipids reduced, 12% increased, and 53% unchanged (Figure 3b). Detailed lipid species distributions within each major category (for categories with ≥10 species) revealed that, PG and LacCer were elevated during sepsis, whereas PC, CE, and LPC were notably reduced. (Figure 3c).

Figure 3

We further analyzed the correlation between 38 lipid subclasses and age in both healthy controls and sepsis patients. Our findings revealed that the majority of lipid subclasses exhibited varying correlations with age across different health states, with some even showing opposite trends. HexCer correlated positively with age in HC and negatively with age in Sepsis; LPC correlated negatively with age in HC and positively correlated with age in Sepsis. This suggests that the pathological state of sepsis differentially impacts lipid metabolism in patients of different ages(Figure 3d).

All blood samples were collected at hospital admission prior to the initiation of systematic treatment—including antibiotic therapy and other clinical interventions. This sampling approach minimizes acute in-hospital treatment effects on lipidomic profiles; however, potential influences from pre-existing comorbidities and pre-admission medications cannot be entirely excluded.

To evaluate the potential confounding effect of metabolic comorbidities, we conducted subgroup principal component analysis (PCA) comparing lipidomic profiles between patients with and without hypertension or diabetes. No significant separation was observed between these subgroups, either in the full lipidome (P = 0.939) or in the sepsis-specific lipid subset (P = 0.894) (Supplementary Figure S1), indicating that the lipid alterations identified were primarily attributable to sepsis itself rather than to underlying metabolic disorders.

Given the high prevalence of pulmonary infection among sepsis cases in our cohort, we further performed subgroup PCA to compare lipidomic profiles between patients with pulmonary versus non-pulmonary infection sources. The global lipidome (n = 1,277 lipids) showed no significant separation between groups (P = 0.762). Similarly, focused analysis of sepsis-altered lipids (n = 603) revealed overlapping distributions (P = 0.851) (Supplementary Figure S2). The absence of distinct clustering suggests that the site of infection does not substantially influence the core sepsis-associated lipidomic signature, supporting the systemic nature of lipid dysregulation in sepsis irrespective of the primary infection origin.

Age stratification analysis in sepsis patients

We analyzed the age-related differences in sepsis patients by stratifying them into two groups: 41 cases of aging sepsis (asepsis) and 21 cases of young sepsis (ysepsis). No significant differences were observed between the asepsis and ysepsis groups in principal component analysis (PCA) (P > 0.05, Figure 4a). There were 24 differentially expressed lipid molecules between the asepsis and ysepsis groups (|FC| > 1.3, P < 0.05), with 5 upregulated and 19 downregulated in the asepsis compared to the ysepsis (Figure 4b).

Figure 4

To further elucidate the impact of age on plasma lipid profiles during the onset of sepsis, we conducted fuzzy clustering analysis based on the 603 differentially expressed lipid molecules identified between sepsis patients and healthy controls. Clustering trend plots were generated to visualize the distribution patterns of mean lipid levels within each clustering module (Cluster) (Figure 4c). The results showed that in Cluster1, both ysepsis and asepsis were upregulated compared to yHC and aHC. In Cluster2, Cluster3, and Cluster4, both ysepsis and asepsis were downregulated compared to yHC and aHC, but there were differences among the clusters. Specifically, in Cluster2, the lipid abundance in yHC was higher than that in aHC, and the changes in ysepsis were greater than those in asepsis. In Cluster3, the changes were similar in magnitude, while in Cluster4, the changes in ysepsis were smaller than those in asepsis.

Using the lipid molecules in each cluster, we performed consensus clustering on sepsis samples using the ConsensusClusterPlus package based on the Euclidean distance between samples, dividing the samples into two subtypes, C1 and C2. We counted the number of samples with shock and non-shock in each subtype and performed Fisher’s statistical test(Figure 4d). The subtyping based on the lipid molecules in Cluster4 showed a significant difference between shock and non - shock samples (P = 0.0013), with a higher proportion of shock patients in subtype C2 (C1:14.8%, C2:57.1%).

To elucidate the potential mechanisms underlying the lipid-immune axis, we performed KEGG pathway enrichment analysis based on DIA-based untargeted proteomic data from the C1 and C2 patient subgroups. The results revealed significant alterations in several key pathways: the ​​”Complement and coagulation cascades”​​ pathway and the ​​”Cholesterol metabolism”​​ pathway were both markedly dysregulated, suggesting a direct mechanistic link between lipid metabolic disturbance and age-related immune dysfunction. Furthermore, pronounced enrichment of the ​​”Chemical carcinogenesis - reactive oxygen species”​​ pathway indicated that oxidative stress may serve as a critical biological mediator connecting lipidomic alterations to immunosenescence(Supplementary Figure S3).

Age-related Cluster4 and the severity of sepsis

The median age in the C1 group was 67 (46-77), of which asepsis accounted for 55.6% (15/27); in the C2 group, the median age was 71 (54-80), of which asepsis accounted for 71.4% (25/35). Notably, gender distribution was comparable between C1 (66.7% male) and C2 (65.7% male), indicating minor sex-based bias in the lipid-defined subtypes. The clustering of the two subtypes (C1 and C2) in Cluster 4 exhibits significant differences, as demonstrated by the PCA plot (P = 9.7e-11) (Figure 5a). There are 92 differentially expressed lipid molecules (|FC| > 1.3, P < 0.05) between the C1 and C2 groups, all of which are downregulated in C2(Figure 5b). Based on these 92 differentially expressed lipid molecules, linear regression analysis was performed to examine their relationship with SOFA scores, both at the subclass and individual molecule levels. Using standardized lipid subclass or individual molecule expression as the independent variable and SOFA score as the dependent variable, a model was constructed using the Ordinary Least Squares (OLS) method to calculate regression coefficients and their corresponding P-values. The results revealed that LPC, CE, and PC subclasses were significantly negatively correlated with SOFA scores (P < 0.05) (Figure 5c), suggesting that increased expression of these subclasses may be associated with reduced SOFA scores and could play a positive role in improving patient prognosis. At the individual lipid molecule level, 35 lipid molecules showed significant negative correlations with SOFA scores (P < 0.05) (Figure 5d), further indicating their potential role in clinical outcomes.

Figure 5

We propose that C2 represents a more severe subtype of sepsis. Further investigation into the associations of LPC, CE, and PC with clinical indicators revealed differences between C1 and C2. Notably, CE showed a negative correlation with monocytes percentage in C1 but a positive correlation in C2. Both CE and LPC exhibited significant positive correlations with complement C3 and C4 in C2 group(Figure 5e).

We attempted to construct a classification model to predict whether sepsis patients would develop shock using the aforementioned 35 lipid molecules associated with SOFA scores, while also incorporating the age factor of the patients. We screened 35 molecules that showed significant linear correlations with SOFA and divided the samples into training and testing sets in a 6:4 ratio. Using the training set, we performed Lasso regression analysis on the 35 lipids and selected those lipids whose regression coefficients were not penalized to zero, ultimately identifying six lipid molecules: LPC(19:0), PC(P-19:0), SM 32:3;2O(FA 16:3), PC(P-20:0), PC(O-18:1/20:3), and CE(15:0) (Figure 6c). These six lipid molecules, along with age, were used to build the model. The final variables were used to establish a logistic regression model on the training set, employing a K-fold cross-validation method with k=5 and repeats=200. The absolute value of the t-statistic for each model parameter was used to represent the importance of the variables.

Figure 6

The results showed that the model could effectively distinguish between sepsis and septic shock patients, with an AUC of 0.87 for the training set and 0.822 for the testing set(Figure 6a). In the PCA plot, the horizontal coordinate represents the predicted probability of shock, with a threshold of 0.5 (values below 0.5 classified as nonShock and above 0.5 as shock). The actual distribution of shock patients largely aligned with the predictions, with most orange points located in the region above 0.5(Figure 6d). The model demonstrated predictive capability in both elderly and young sepsis patients. The ysepsis group accurately predicted 75%(6/8) of the shocked patients, and the asepsis group accurately predicted 81%(13/16) of the shocked patients.At the level of lipid subclasses to which the 35 molecules belong, PC contributed the most to the model, followed by age, CE, and LPC(Figure 6b). PC was positively correlated with the absolute monocyte count (R = 0.33, P = 0.0088) and monocyte percentage (R = 0.37, P = 0.0035) in sepsis patients, and negatively correlated with commonly used clinical inflammatory markers CRP (R=-0.48, P = 1e-04) and PCT (R=-0.36, P = 0.0056) (Figure 6g).

The distribution of the six lipid molecules across groups is shown in bar charts, and their correlations with clinical indicators are displayed in bubble plots. All six molecules had the highest plasma levels in healthy controls, reduced plasma levels in sepsis and the lowest plasma levels in septic shock(Figure 6e). They were negatively correlated with the inflammatory markers CRP, PCT and IL-6, positively correlated with complement C3, complement C4, percentage of lymphocytes, percentage of monocytes, and inversely correlated with absolute neutrophil count(Figure 6f).

Utility of lipid model for predicting mortality in patients with sepsis

To further explore whether the above lipid model could identify patients’ in-hospital prognosis, we divided sepsis patients into in-hospital survival and in-hospital death groups to observe the classification efficacy of the model. We found that age was the most important predictor of patient prognosis, both at the lipid molecular level and at the lipid subclass level(Figures 7a, b). In young sepsis patients, the model was unable to predict patients’ in-hospital death, but in older sepsis patients, the model correctly identified 65% (11/17) of patients with in-hospital death(Figure 7c). We also analysed the association of the high-contributing lipid subclasses CE and LPC with clinical immune markers in patients, and both CE and LPC were positively correlated with complement C3, C4 and CD8+ cells, which suggests that complement C3, C4 and CD8+ cells correlate with the prognosis of elderly sepsis patients(Figure 7d).

Figure 7

Discussion

This study provides the first systematic evidence of ​age-specific differences​ in plasma lipidomic profiles among sepsis patients. While no significant differences in overall lipid abundance were observed between elderly and young healthy controls, ​sepsis-induced lipidomic remodeling exhibited marked heterogeneity across age groups.Notably, we identified ​severity-associated lipid subclasses​ (PCs, CEs, and LPCs) with distinct age-associated patterns, suggesting that ​lipid metabolism and immune regulation may be differentially dysregulated in sepsis depending on patient age.

Our lipidomic analysis revealed that sepsis patients exhibit a marked reduction in key lipid subclasses—including phosphatidylcholines (PCs), cholesteryl esters (CEs), and lysophosphatidylcholines (LPCs)—which play critical roles in maintaining membrane stability, facilitating cholesterol transport, and modulating immune responses (). Notably, clustering analysis of age-stratified lipid profiles identified Cluster 4 lipids as robust discriminators of sepsis subtypes (C1 vs. C2), with the C2 subtype demonstrating significantly higher shock incidence (57.1% vs. 14.8%). Building upon these findings, we developed a risk stratification model incorporating six lipid species [LPC(19:0), PC(P-19:0), SM 32:3;2O(FA 16:3), PC(P-20:0), PC(O-18:1/20:3), CE(15:0)] and patient age, which effectively stratified shock risk across age groups (AUC: 0.87 training, 0.82 validation) and predicted in-hospital mortality in elderly patients. Mechanistically, we observed that PC levels positively correlated with monocyte abundance, while CE and LPC levels were associated with complement factors (C3, C4) and CD8+ T cell proportions. These results collectively suggest that age-specific lipid remodeling may reflect distinct immunoregulatory pathways in sepsis, positioning lipidomic signatures as both biologically informative markers and clinically actionable predictors of outcomes (Figure 8).

Figure 8

These lipidomic changes may reflect broad disturbances in the metabolic and physiological states of sepsis patients (, ). Although it is challenging to directly attribute these changes to specific biological processes, we hypothesize that the reductions in PC, CE, and LPC are directly linked to the inflammatory responses associated with sepsis and the loss of cellular membrane integrity (, ). The importance of PC and LPC in cellular membranes and their roles in signal transduction suggest that their decreased levels may lead to cellular dysfunction, exacerbating the pathological state of sepsis (, ). Previous studies have indicated that LPC is essential for maintaining T-cell homeostasis (), and its depletion leads to a reduction in the memory T-cell pool, thereby diminishing the secondary response to reinfections (). Additionally, LPC induces chemotaxis and intracellular calcium influx in natural killer cells (). In fact, LPC levels are inversely correlated with mortality in sepsis patients, indicating that the lower LPC levels observed in older sepsis patients in our study may predict poorer outcomes (, ).

Our analysis revealed a significant association between ether-linked (lyso)phosphatidylcholine species and SOFA scores, suggesting their potential role as biomarkers of sepsis severity. Ether phospholipids, particularly plasmalogens, are known to function as endogenous antioxidants that protect cellular membranes against oxidative stress—a hallmark of sepsis pathophysiology (). Their depletion may reflect both increased consumption during radical scavenging and impaired synthesis due to peroxisomal dysfunction in severe sepsis (). The more pronounced reduction of these ether lipids in elderly patients may further indicate an age-associated decline in antioxidant capacity, potentially contributing to the worse outcomes observed in this population (, ). While our study demonstrates strong clinical correlations, further mechanistic studies are needed to elucidate whether ether lipid supplementation could serve as a therapeutic strategy to mitigate oxidative damage in sepsis.

Moreover, a reduction in the CE may impair cholesterol transport and utilization, further weakening the immune response (). CE is the storage form of cholesterol and is involved in cell membrane stability and cholesterol transport. CETP is responsible for the transport of cholesteryl esters in the body, and sepsis can lead to alterations in CETP expression (), subsequently affecting cholesterol transport and CE synthesis (, , ). Venancio et al. proposed a potential protective role for CETP in regulating inflammatory responses, whereas Dusuel et al. suggested that human CETP may worsen outcomes in sepsis (). Therefore, further studies are needed to elucidate the role of CE in sepsis and its relationship with CETP. Notably, the accumulation of CEs is crucial for sustaining T-cell proliferation in humans (), and its decline may directly affect host immune efficacy.

Additionally, we investigated the relationship between lipid changes and immune response. CRP and PCT are widely used inflammatory markers in sepsis, whose levels are typically significantly elevated during infection and inflammatory responses (, ). The negative correlation between PC and CRP/PCT may be attributed to the anti-inflammatory effects of PC in inflammatory responses (). Specifically, PC modulates membrane fluidity, thereby reducing the release of inflammatory mediators (). Additionally, PC metabolites, such as LPC, can inhibit the activation of inflammatory cells, alleviating inflammatory responses (). As a metabolite of PC, LPC exhibits multiple biological activities, including the regulation of inflammatory responses and cellular signaling. Both CE and LPC may influence the activation of the complement system by modulating the lipid microenvironment of cell membranes. For instance, LPC can induce the chemotaxis of neutrophils and monocytes, thereby promoting complement system activation (). Meanwhile, CE metabolites may indirectly affect complement system function by regulating cholesterol transport. Furthermore, LPC can activate CD8+ T cells through its receptors (e.g., G2A) (), while CE metabolites may indirectly modulate CD8+ T cell function by influencing cholesterol transport ().

In our study, patients in subtype C2 exhibited higher severity and a greater proportion of elderly individuals. The lipid subclasses CE, PC, and LPC, which are associated with SOFA scores, were found to be decreased in this subtype. Given the correlations between these lipids and immune cells, we hypothesize that the pathogenesis of severe sepsis in the elderly may be related to reductions in complement C3, complement C4, monocytes, and CD8+ T cells. In the elderly, studies have indicated the presence of complement changes and a low-grade inflammatory state, which is the same trend as in the elderly critically ill patients in our study (). This also suggests the reason for the higher probability of shock in elderly septic patients.These age-related differences in immune response mechanisms suggest that clinical treatment strategies for sepsis may need to be tailored on the basis of the patient’s age to better activate and maintain immune function (, –). Although our study revealed correlations between lipidomic profiles and certain immune parameters, we acknowledge that the immune profiling was limited in scope. Key markers of immune dysfunction in sepsis, such as HLA-DR expression on monocytes, T cell activation/exhaustion markers (e.g., CD28, CD38), were not assessed. Future studies incorporating broader immunophenotyping and cytokine profiling would provide a more comprehensive understanding of lipid-immune interactions in sepsis.

Despite presenting new insights, our study has several limitations. The sample size was relatively small, and there was no multicenter validation. Future studies should involve larger-scale clinical observations and validations while linking lipidomics with detailed immunosenescence profiling (, ).Although we minimized acute treatment effects through admission-time sampling, our study cannot exclude potential influences from long-term medications (e.g., statins, antihypertensive drugs) used for comorbid conditions. The lack of systematic medication data represents a limitation, particularly given that: Metformin and other antidiabetic medications are known to alter circulating lipid profiles (). β-blockers and ACE inhibitors may affect phospholipid metabolism (). Statins directly modulate cholesterol biosynthesis pathways (). Future longitudinal studies with detailed medication records are needed to fully disentangle sepsis-specific lipid changes from pre-existing metabolic influences. While our study identified significant lipid alterations associated with sepsis severity and outcomes, the extent to which these changes are specific to sepsis—rather than general critical illness or systemic inflammation—requires further investigation.

Our study captures lipidomic profiles at a single timepoint during acute sepsis, aiming to identify state biomarkers associated with disease severity rather than to establish causal relationships. Importantly, while the temporal origin of these lipid alterations (cause vs. consequence) remains undetermined, their strong association with clinical outcomes supports their utility as prognostic indicators. To move beyond correlation and establish biological mechanism, future work will focus on functional validation using experimental models such as LPS-stimulated PBMCs or relevant cell lines. Specifically, we will investigate whether supplementation of key depleted lipids (e.g., LPC(19:0), CE(15:0)) can modulate critical immune functions including cytokine production, phagocytosis, and T-cell activation, thereby providing direct evidence for their immunomodulatory roles in sepsis.

In summary, our research highlights the complex relationship between lipidomic changes and clinical indicators in sepsis patients, particularly as influenced by age (). These findings increase our understanding of the pathophysiology of sepsis and provide potential directions for future treatment strategies, especially when age-related factors are considered (). Further investigations should explore immune response mechanisms across different age groups to offer new perspectives for sepsis treatment ().

Conclusion

Our study utilized nontargeted LC-MS/MS lipidomics to reveal significant age-specific differences in plasma lipid profiles among sepsis patients. Furthermore, we found that lipid molecules associated with age during sepsis correlate with high SOFA scores, potentially elucidating the reasons for poorer outcomes in older patients. We developed a robust risk stratification model integrating six lipid biomarkers (LPC(19:0), PC(P-19:0), SM 32:3;2O(FA 16:3), PC(P-20:0), PC(O-18:1/20:3), CE(15:0)) and patient age, accurately forecasting the risk of septic shock across age groups and predicting in-hospital mortality specifically in elderly patients. In sepsis subtypes characterized by older age and greater severity, decreased plasma levels of lipid subclasses CE, PC, and LPC likely reflect impaired immune functions, such as reductions in complement factors (C3, C4), monocytes, and CD8+ T cells, which may contribute mechanistically to the rapid clinical deterioration observed in elderly patients. Our findings highlight the clinical significance of plasma lipidomic profiling, pointing toward novel therapeutic interventions aimed at modulating host lipid metabolism to improve outcomes in sepsis. The stratified analysis underscores the importance of precise clinical diagnosis and treatment, as differing pathophysiological changes may require tailored interventions. Additionally, this research provides a reference for identifying lipid biomarkers associated with immunosenescence in disease states.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by Institutional Review Board of the Ethics Committee at Huashan Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

YS: Formal Analysis, Writing – original draft, Writing – review & editing. KL: Conceptualization, Resources, Writing – review & editing. LW: Data curation, Methodology, Writing – original draft. JA: Conceptualization, Resources, Writing – review & editing. JZ: Software, Writing – review & editing. RS: Investigation, Writing – review & editing. XC: Investigation, Writing – review & editing. YW: Conceptualization, Supervision, Writing – review & editing. PC: Validation, Writing – review & editing. SW: Supervision, Writing – original draft, Writing – review & editing. JW: Conceptualization, Writing – review & editing. JJ: Funding acquisition, Supervision, Visualization, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by the National Key Research and Development Program of China (2022YFC2009802, 2022YFC2009801, 2022YFC2009806), Science and Technology Innovation Project of Shanghai (20Y11900400), Shanghai Hospital Development Center, Clinical Science and Technology Innovation Project (SHDC22020214), the Shanghai Municipal Health Commission (GWVI-11.1-07), and the Shanghai Municipal Science and Technology Major Project (HS2021SHZX001).

Acknowledgments

We thank all the patients who volunteered for this trial and the study site personnel for their contributions.

Conflict of interest

Author LW was employed by the company Shanghai Omicsolution Co., Ltd.

The remaining 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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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/fimmu.2025.1659425/full#supplementary-material

Supplementary Figure 1

Lipidomic profiles stratified by comorbidity status. Principal component analysis (PCA) of lipidomic profiles from sepsis patients with (yellow) and without (blue) metabolic complications (hypertension and/or diabetes). ​PCA using all detected lipid species (n=1277) shows no significant separation between groups (P = 0.939). PCA using only sepsis-associated differential lipids (n=603) similarly demonstrates no significant group separation (P = 0.894). Percentages indicate the proportion of variance explained by each principal component (PC1/PC2). The absence of significant separation suggests that metabolic comorbidities did not substantially alter the sepsis-associated lipidomic signatures identified in this study.

Supplementary Figure 2

Lipidomic profiles stratified by infection site.​​ Principal component analysis (PCA) of plasma lipidomes from sepsis patients with pulmonary infections (red squares) versus others infections (blue circles). Global lipidome analysis (n=1,277 lipids) revealed no significant separation between groups (P = 0.762). Focused analysis of sepsis-altered lipids (n=603) similarly showed overlapping distributions (P = 0.851). The lack of distinct clustering suggests that infection site does not substantially influence the core sepsis-associated lipidomic signature, supporting the systemic nature of lipid dysregulation in sepsis regardless of primary infection origin.

Supplementary Figure 3

KEGG pathway enrichment analysis of proteomic profiles between C1 and C2 sepsis subgroups.​​ Bubble plot showing significantly enriched pathways (Y-axis) based on untargeted proteomics data. Bubble size corresponds to the number of proteins mapped to each pathway (Count), and color intensity represents the Ratio (proportion of significant proteins in the pathway). The X-axis indicates the statistical significance (P-value).

Abbreviations

PC, Phosphatidylcholine; LPC, Lysophosphatidylcholine; CE, Cholesteryl Ester; SM, Sphingomyelin; PE, Phosphatidylethanolamine; TG, Triglyceride; SOFA, Sequential Organ Failure Assessment; CRP, C-Reactive Protein; PCT, Procalcitonin; IL-6, Interleukin-6; WBC, White Blood Cell Count; N%, Neutrophil Percentage; aHC, Aged Healthy Controls; yHC, Young Healthy Controls; Asepsis, Aged Sepsis Patients; Ysepsis, Young Sepsis Patients.

References

Summary

Keywords

sepsis, lipidomics, aging, immunity, shock

Citation

Sun Y, Lin K, Wang L, Ai J, Zhao J, Sun R, Cheng X, Wan Y, Cui P, Wang S, Wu J and Jin J (2025) Plasma lipidomic signatures reveal age-associated patterns of septic shock risk and immune dysregulation in sepsis. Front. Immunol. 16:1659425. doi: 10.3389/fimmu.2025.1659425

Received

04 July 2025

Accepted

24 October 2025

Published

19 November 2025

Volume

16 - 2025

Edited by

Vlad Pădureanu, University of Medicine and Pharmacy of Craiova, Romania

Reviewed by

Sandrine Horman, Université Catholique de Louvain, Belgium

Longxiang Su, Peking Union Medical College Hospital (CAMS), China

Updates

Copyright

*Correspondence: Jialin Jin, ; Jing Wu, ; Sen Wang,

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics