Skip to main content
Download PDF

Metabolomics and machine learning approaches for diagnostic and prognostic biomarkers screening in sepsis

BMC Anesthesiology volume 23, Article number: 367 (2023) Cite this article

Abstract

Background

Sepsis is a life-threatening disease with a poor prognosis, and metabolic disorders play a crucial role in its development. This study aims to identify key metabolites that may be associated with the accurate diagnosis and prognosis of sepsis.

Methods

Septic patients and healthy individuals were enrolled to investigate metabolic changes using non-targeted liquid chromatography-high-resolution mass spectrometry metabolomics. Machine learning algorithms were subsequently employed to identify key differentially expressed metabolites (DEMs). Prognostic-related DEMs were then identified using univariate and multivariate Cox regression analyses. The septic rat model was established to verify the effect of phenylalanine metabolism-related gene MAOA on survival and mean arterial pressure after sepsis.

Results

A total of 532 DEMs were identified between healthy control and septic patients using metabolomics. The main pathways affected by these DEMs were amino acid biosynthesis, phenylalanine metabolism, tyrosine metabolism, glycine, serine and threonine metabolism, and arginine and proline metabolism. To identify sepsis diagnosis-related biomarkers, support vector machine (SVM) and random forest (RF) algorithms were employed, leading to the identification of four biomarkers. Additionally, analysis of transcriptome data from sepsis patients in the GEO database revealed a significant up-regulation of the phenylalanine metabolism-related gene MAOA in sepsis. Further investigation showed that inhibition of MAOA using the inhibitor RS-8359 reduced phenylalanine levels and improved mean arterial pressure and survival rate in septic rats. Finally, using univariate and multivariate cox regression analysis, six DEMs were identified as prognostic markers for sepsis.

Conclusions

This study employed metabolomics and machine learning algorithms to identify differential metabolites that are associated with the diagnosis and prognosis of sepsis patients. Unraveling the relationship between metabolic characteristics and sepsis provides new insights into the underlying biological mechanisms, which could potentially assist in the diagnosis and treatment of sepsis.

Trial registration

This human study was approved by the Ethics Committee of the Research Institute of Surgery (2021–179) and was registered by the Chinese Clinical Trial Registry (Date: 09/12/2021, ChiCTR2200055772).

Peer Review reports

Introduction

Sepsis is a significant public health concern characterized by a dysregulated host response to infection, resulting in life-threatening organ dysfunction [1]. The global burden of sepsis is increasing, with approximately 45 million new cases and 11 million deaths reported in 2017 [2]. Accurate definition, diagnosis, and early recognition of sepsis are crucial for improving treatment effects of patients. In addition to the early diagnostic symptoms and signs, the identification of available biomarkers is essential for monitoring and treating sepsis. Several molecules, including lactic acid, c-reactive protein (CRP), procalcitonin (PCT), and B-type natriuretic peptide (BNP), have been proposed as candidate sepsis biomarkers. However, the number of useful predictive biomarkers for assessing the severity and prognosis of sepsis in clinical practice remains limited [3,4,5]. Therefore, it is urgent to find specific biomarkers for sepsis.

Metabolomics is a novel approach developed following genomics, transcriptomics, proteomics and lipidomics [6, 7]. It involves the systematic identification and quantification of a wide range of small molecule metabolites in biological samples. Since the research focus is on endogenous molecule compounds at the end of metabolic pathways, this approach directly reflects the body's state, thereby offering a new means of identifying potential biomarkers [8]. Metabolomics has been extensively applied to investigate the pathologic mechanism and biomarker of sepsis. Mickiewicz et al. employed targeted metabolomics based on nuclear magnetic resonance (NMR) to identify 186 metabolites in the serum of patients in the Intensive Care Unit, suggesting that metabolomics holds promise for predicting mortality in septic shock [9]. Another study found that two metabolites could distinguish severe sepsis from systemic inflammatory response syndrome [10]. Liquid chromatography and tandem mass spectrometry (LC–MS/MS) combines high performance liquid chromatography with electrospray ionization mass spectrometry (ESI–MS) metabolomics technology to achieve comprehensive detection of different types of metabolites in samples [11]. In comparison to targeted metabolomics using NMR, LC–MS/MS offers the advantages of heightened sensitivity and a broader dynamic range [12].

Machine learning is a key area within the field of artificial intelligence [13, 14], representing an algorithmic framework that enables intelligent processing of data. Through feature extraction, machine learning continually enhances its performance by automatically learning internal data patterns. The main algorithms include artificial neural network (ANN), support vector machine (SVM), random forest (RF), and decision tree (DT) [15, 16]. Utilizing machine learning opens up new possibilities for enhancing diagnostic efficiency and achieving more objective and personalized patient assessments, deviating from traditional diagnostic and treatment methods [17]. Despite its increasing prominence in the diagnostic field, it is rarely used to identify potential diagnostic and prognostic goals for sepsis.

The present study employed liquid chromatography quadrupole time of flight mass spectrometry (LC-QTOF/MS) to analyze the serum metabolites of sepsis patients. Machine learning techniques were then utilized to identify potential biomarkers for the diagnosis and prognosis of sepsis.

Materials and methods

Ethical review

This human study was approved by the Ethics Committee of the Research Institute of Surgery (2021–179) and was registered by the Chinese Clinical Trial Registry (Date: 09/12/2021, ChiCTR2200055772). All the participants provided written informed consent before inclusion in this study. All experimental procedures were approved by the Animal Ethics Committee of Army Medical University, China, and were performed in accordance with the National Institutes of Health “Guidelines for the Care and Use of Laboratory Animals”. The study was carried out in compliance with the ARRIVE guidelines. All methods are reported in accordance with ARRIVE guidelines (https://arriveguidelines.org) for the reporting of animal experiments.

Study design and population recruitment

From December 2021 to April 2022, a total of 30 septic patients were recruited from the Intensive Care Unit (ICU) of Daping Hospital, Army Medical University, based on the Sepsis-3 criteria for sepsis and septic shock. Additionally, 15 age-matched healthy volunteers were enrolled from the State Key Laboratory of Trauma, Burns, and Combined Injury. The inclusion and exclusion criteria for septic patients are presented in Supplementary Table 1. Healthy participants also follow the similar inclusion and exclusion criteria. Blood samples were collected from all septic patients within 24 h of admission, while blood samples from healthy volunteers were obtained during enrollment.

Blood sample collection and metabolites extraction

A 10 mL peripheral venous blood sample was collected from all participants in both groups. Subsequently, the samples were centrifuged at 3000 rpm for 10 min, and the resulting supernatant was aliquoted into 1.5 mL centrifuge tubes and stored at -80 °C. For analysis, 100 μL of each sample was transferred to an EP tube. To this, 400 μL of extract solution (methanol containing isotopically-labelled internal standard mixture) was added, followed by 30 s of vortexing, 10 min of sonication in an ice-water bath, and 1 h of incubation at -40 °C to precipitate proteins. The sample was then centrifuged at 12,000 rpm (RCF = 13,800(× g), R = 8.6 cm) for 15 min at 4 °C. The resulting supernatant was transferred to a fresh glass vial for analysis. A quality control (QC) sample was prepared by combining equal aliquots of the supernatants from all the individual samples.

LC–MS/MS analysis

LC–MS/MS analyses were conducted using a UHPLC system (Vanquish, Thermo Fisher Scientific) coupled to an Orbitrap Exploris 120 mass spectrometer (Orbitrap MS, Thermo) with a UPLC HSS T3 column (2.1 mm × 100 mm, 1.8 μm). The mobile phase consisted of 5 mmol/L ammonium acetate and 5 mmol/L acetic acid in water (A) and acetonitrile (B). The auto-sampler temperature was maintained at 4 °C, and the injection volume was 2 μL. The Orbitrap Exploris 120 mass spectrometer was selected for its capability to acquire MS/MS spectra in information-dependent acquisition (IDA) mode, controlled by Xcalibur software (Thermo). In this mode, the acquisition software continuously evaluates the full scan MS spectrum. The ESI source conditions were set as follows: sheath gas flow rate at 50 Arb, aux gas flow rate at 15 Arb, capillary temperature at 320 °C, full MS resolution at 60,000, MS/MS resolution at 15,000, collision energy at 10/30/60 in NCE mode, and spray voltage at 3.8 kV (positive) or -3.4 kV (negative), respectively.

Data preprocessing and annotation

The raw data were converted to the mzXML format using ProteoWizard and processed with an in-house program developed in R. This program utilized XCMS for peak detection, extraction, alignment, and integration. Metabolite annotation was performed using an in-house MS2 database called BiotreeDB, with a cutoff set at 0.3 for annotation.

Identification of differentially expressed metabolites via RF and SVM Algorithms

In this study, various machine learning algorithms were employed to identify differentially expressed metabolites (DEMs). The random forest (RF) algorithm was utilized to determine the optimum DEMs between the sepsis and healthy control groups. Using the "e1071" R package, the support vector machine (SVM) generated eigenvectors were eliminated to extract the optimal variables for identifying diagnostic DEMs in sepsis. Additionally, a SVM classifier was constructed with tenfold cross-validation. Finally, a Venn diagram was employed to identify the shared DEMs.

Dataset collection

A total of 18 genes related to phenylalanine metabolism were obtained from the Molecular Signatures Database (MSigDB) (https://www.gsea-msigdb.org/gsea/) for further investigation. The mRNA expression matrix was downloaded from the Gene Expression Omnibus (GEO) dataset (https://www.ncbi.nlm.nih.gov/geo/). In this study, the GSE57065 dataset was utilized for screening differentially expressed genes (DEGs). The R package "SVA" was employed to normalize the gene expression matrices and eliminate batch effects. Perl scripts were used to convert the dataset probes into corresponding gene symbols.

Sepsis model establishment

The ethics and protocols for the animal experiments were approved by the La-boratory Animal Welfare and Ethics Committee of the Army Medical University (Approval No. AMUWEC20224867). Adult male and female Sprague–Dawley rats weighing between 200-220 g were bred in a facility with filtered positive-pressure ventilation, under a 12:12-h dark/light cycle, with ad libitum access to food and water. The rats were randomly assigned to three groups: Control (n = 30), Sepsis (n = 30), and RS8359-treated Sepsis (n = 30). The RS8359 compound (HY-14260, Med-ChemExpress) was administered at a dosage of 5 mg/kg, dissolved in ddH2O, through the tail vein 30 min prior to inducing sepsis in the RS8359-treated group. The Sepsis group received an equivalent volume of ddH2O. The rats were anesthetized with sodium pentobarbital (45 mg/kg, i.p.) for the surgical procedure. Cecal ligation and puncture (CLP) were conducted to create the sepsis model following established protocols [8]. Myocardial tissues were harvested at 12 h after CLP and stored at − 80 °C until further processing.

Statistical analysis

All statistical analyses were conducted using R software (version 4.1.2, http://www.r-project.org). The DEGs were identified using the "Limma" package with a significance threshold of P < 0.05. Receiver operating characteristic (ROC) analysis was performed, and the area under the curve (AUC) was calculated to assess the predictive performance of the classifier.

Results

Clinical characteristics of the studied population

Our study included 30 septic patients and 15 healthy controls. The patients in the sepsis group had a mean age of 60.1 ± 14.5, while the control group had a mean age of 50.4 ± 17.2 (p > 0.05). There were no significant differences in gender and BMI among the two groups (p > 0.05). However, we observed significant differences in temperature (p = 0.01), heart rate (p = 0.001), and white blood cell count (p = 0.003) when comparing septic patients to the healthy controls (Table 1).

Table 1 Baseline characteristics of the septic patients and healthy controls in this study
Full size table

Metabolic reprogramming related differentially expressed metabolites participated in the occurrence of sepsis

The ionization source of Orbitrap is electrospray ionization, including positive and negative ion modes (POS and NEG, respectively), which can increase the metabolite coverage and improve the detection effect. Subsequent data analysis involved separate analyses of the two groups. Principal component analysis (PCA) revealed significant differences in metabolites between sepsis patients (SP) and healthy controls (HC) (Fig. 1A, B). Further assessment using Orthogonal Partial Least Squares-Discriminant Analysis (OPLS-DA) demonstrated the robustness of the original model without any overfitting observed (Fig. 1C, D). Metabolomics analysis detected a total of 532 differentially expressed metabolites (DEMs) (VIP > 1, P < 0.05) between the sepsis and control groups, with 327 identified in the positive ion mode and 205 in the negative ion mode (Fig. 1G). The pie plots based on super class classification revealed that the DEMs primarily belonged to Lipids and lipid-like molecules (37.69% in POS and 42.975% in NEG) and Organic acids and derivatives (11.294% in POS and 21.074% in NEG) (Fig. 1 E, F). Identified through radar plots, the top 10 DEMs with the highest fold change between sepsis patients and healthy controls in POS were (15a,20R)-Dihydroxypregn-4-en-3-one 20-[glucosyl-(1- > 4)-6-acetyl-glucoside], Acar (14:1), N1-Methyl-2-pyridone-5-carboxamide, Kanzonol K, Ricinine, Zapotin, L-trans-4-Methyl-2-pyrrolidinecarboxylic acid, 2-Ethyl-4,5-dimethyloxazole, L-glycyl-L-hydroxyproline, and Phosphoric acid. Meanwhile, in NEG, the top 10 DEMs were 16-Methylheptadecanoic acid, Eicosadienoic acid, Oleic acid, FA (18:1), FA (18:0), FA (19:0), FA (19:1), FA (17:0), FA (22:4), and FA (19:2) (Fig. 1H).

Fig. 1
figure 1

Identification of differentially expressed metabolites (DEMs) between control and sepsis. A, B Principal component analysis (PCA) scores plot for metabolomics analysis in the sepsis and control groups. PC1 and PC2 in the figure represent the scores of the first and second principal components respectively. Each scatter represents a sample. The red circle represents septic patients, and the blue circle represents the healthy controls. C, D OPLS-DA replacement test of NEG and POS. The abscissa represents the displacement retention degree of the displacement test, the ordinate represents the value of R2Y or Q2, the green dot represents the R2Y value obtained from the displacement test, the blue square dot represents the Q2 value obtained from the displacement test, and the two dotted lines represent the regression lines of R2Y and Q2 respectively. Pie charts of super class of DEMs in NEG E and POS (F) models. (G) Heat maps and (H) Radar plots analyzed by TBtools showing the significantly changed metabolites in sepsis and control

Full size image

Figure 2A and B illustrate that the differentially expressed metabolites (DEMs) pathways primarily involved Phenylalanine metabolism, Pyruvate metabolism, and Cysteine and methionine metabolism in the NEG model. On the other hand, in the POS model, the DEMs were enriched in Linoleic acid metabolism and Arginine and proline metabolism. To further assess the functional implications of these metabolites, we conducted Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. In the NEG model, the most significant KEGG pathways included biosynthesis of amino acids, phenylalanine metabolism, cysteine and methionine metabolism, as well as glycine, serine, and threonine metabolism. Meanwhile, in the POS model, KEGG analysis revealed significant enrichment in pathways related to tyrosine metabolism, arginine and proline metabolism, and linoleic acid metabolism (Fig. 2C, D). These findings indicate a close association between sepsis progression and metabolic reprogramming.

Fig. 2
figure 2

Correlation and function enrichment analysis of DEMs. Pathway analysis of DEMs in (A) NEG and (B) POS models. Each square in the rectangular tree represents a metabolic pathway. The size of the square represents the impact degree of the pathway, and the color represents the P value. C, D Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEMs in NEG and POS. The vertical axis represents the names of the enriched KEGG metabolic pathways, and the horizontal axis represents the percentage of annotated metabolites in each pathway relative to the total number of metabolites

Full size image

Identification of key DEMs in sepsis via various machine learning algorithms

We constructed a Random Forest (RF) model using minimum error regression trees for screening differentially expressed metabolites (DEMs) (Fig. 3A). Through this approach, we identified 26 DEMs in both the NEG and POS models. Additionally, we established a support vector machine (SVM) for DEM selection (Fig. 3B). By comparing the results from RF and SVM, we determined 4 key DEMs (3,4-Dihydrocoumarin, Phenol, Benzaldehyde, and DL-Phenylalanine) based on the intersection of these two methods (Fig. 3C). Furthermore, we conducted Pearson's correlation analysis to examine the relationship between these key DEMs and conventional indicators, which is represented in the correlation thermogram (Fig. 3D). The analysis revealed positive correlations between the 4 DEMs and high sensitivity cardiac troponin (hs-CTn), temperature, heart rate, high density lipoprotein (HDL), white blood count, aspartate aminotransferase/alanine aminotransferase (AST/ALT), and triglyceride (TG). Conversely, a negative correlation was observed between the 4 DEMs and creatine, potassium (K+), total cholesterol (TC), and B-type natriuretic peptide (BNP).

Fig. 3
figure 3

Screening of DEMs via the comprehensive strategy. A DEMs screened by random forest (RF) algorithm. B Support vector machine (SVM) was applied for DEMs. C VENN diagram of hub DEMs. D Person’s correlation analysis between hub DEMs and routine detection indexes. *P < 0.05, **P < 0.01, ***P < 0.001

Full size image

The 4 DEMs were diagnostic biomarkers of sepsis

We further screened metabolites by ROC analysis and calculated the area under the curve (AUC). The results showed that the AUC of the 4 key DEMs analyzed by machine learning were both 1.0, which showed an excellent prediction ability for sepsis (Fig. 4A-D). The expressions of the 4 key DEMs (3,4-Dihydrocoumarin, Phenol, Benzaladehyde, and DL-Phenylalanine) were higher in the sepsis group than in the control group (Fig. 4E-H).

Fig. 4
figure 4

Receiver operating characteristic curve (ROC) analysis of the 4 DEMs. ROC curve of (A) 3,4-Dihydrocoumarin, (B) Phenol, (C) Benzaladehyde, and (D) DL-Phenylalanine. The expression of (E) 3,4-Dihydrocoumarin, (F) Phenol, (G) Benzaladehyde, and (H) DL-Phenylalanine between sepsis and control groups

Full size image

The phenylalanine metabolism-related key gene MAOA played important role in sepsis

Due to the enrichment analysis of differential metabolites and machine learning screening showing the significant role of phenylalanine in the progression of sepsis, we further investigated the role and regulatory mechanism of phenylalanine metabolism in sepsis. We analyzed transcriptome data from the GEO database. This dataset included 25 blood samples from healthy controls and 82 blood samples from septic patients. We identified differentially expressed genes (DEGs) between the sepsis and healthy control groups using a threshold of |log2FC|> 1 and p < 0.001. As depicted in Fig. 5A (volcano plot), Monoamine Oxidase A (MAOA), the sole phenylalanine metabolism-related gene, exhibited a significant increase post-sepsis. To demonstrate the impact of MAOA in sepsis, we employed a sepsis model in rats using CLP (cecal ligation and puncture) and observed the effect of MAOA inhibition through its antagonist, RS-8359. RS-8359 was administered through the tail vein at a dosage of 5 mg/kg during sepsis induction. Firstly, we confirmed the inhibitory effect of RS-8359 on MAOA in myocardial tissue via Western blot analysis (Fig. 5B, C). Subsequently, our study revealed that mean arterial pressure (MAP) decreased after CLP, while administration of RS-8359 increased the MAP compared to the sepsis group (Fig. 5D). None of the septic rats survived beyond 24 h, with an average survival time of 7.25 ± 3.49 h. Conversely, in the RS-8359 group, 25% (4/16) of the rats survived for 24 h, with an average survival time of 13.16 ± 7.07 h, which was significantly longer than the sepsis group (Fig. 5E, F).

Fig. 5
figure 5

The role of MAOA in sepsis rats. A Volcano plot of differentially-expressed genes (DEGs) of GSE57065. Statistically significant DEGs were identified as those with a student’s t-test (|logFC|> 1, P < 0.001). B Representative Western blotting images and (C) quantification of Western blotting results of MAOA between the control group, the sepsis group and the RS-8359 group in heart tissues (n = 3 independent experiments). ***P < 0.001 vs control, ###P < 0.001 vs sepsis. D The level of phenylalanine in heart tissues detected by the phenylalanine assay kit (ab241000, Abcam) (n = 6 each group). ***P < 0.001 vs control, #P < 0.05 vs sepsis. (E) MAP within 6 h detecting starting at 12 h after CLP between the control group, the sepsis group and the RS-8359 group (n = 8 each group). ***P < 0.001 vs control, ###P < 0.001 vs sepsis. F The survival rate of rats between the control group, the sepsis group and the RS-8359 group (n = 16 each group). **P < 0.01 vs control, #P < 0.05 vs sepsis. G The survival time of rats in between the control group, the sepsis group and the RS-8359 group. **P < 0.01 vs control, #P < 0.05 vs sepsis

Full size image

DEMs could indicate the prognosis for sepsis

In the study, we followed up with patients for 28 days from the diagnosis of sepsis in the ICU, the 28 days mortality of septic patients was 36.7% (11/30), and the ICU mortality was 13.3% (4/30). Univariate Cox regression and multivariate Cox regression analysis were used to identify the prognostic related DEMs. After observing the survival outcome of septic patients within 28 days, DEMs associated with survival were identified by univariate Cox regression analysis in NEG and POS (Tables 2 and 3). Subsequently, a total 6 prognostic related DEMs, including FA (22:5) (HR = 2.98, p = 0.007), L-Threonic (HR = 0.16, p = 0.006), Cer/AS(d17:3/16:2) (HR = 2.38, p = 0.022), Diethylphosphate (HR = 8.78, p = 0.003), Hydroxycotinine (HR = 2.17, p = 0.009) and Syringol (HR = 0.17, P = 0.004) were identified by multivariate Cox regression analysis (Fig. 6A, B).

Table 2 Univariate Cox regression analysis of DEMs in NEG
Full size table
Table 3 Univariate Cox regression analysis of DEMs in POS
Full size table
Fig. 6
figure 6

Prognostic analysis of DEMs in sepsis. Multivariate Cox regression analysis of DEMs between septic patients and healthy controls in (A) NEG and (B) POS models

Full size image

Discussion

In this study, metabolomics analysis identified a total of 532 differentially expressed metabolites (DEMs). These DEMs were primarily involved in various metabolic pathways, including amino acid biosynthesis, phenylalanine metabolism, tyrosine metabolism, glycine, serine, and threonine metabolism, as well as arginine and proline metabolism. Machine learning algorithms SVM and RF were utilized to identify four metabolic biomarkers (3,4-Dihydrocoumarin, Phenol, Benzaladehyde, and DL-phenylalanine) that are associated with the diagnosis of sepsis. Additionally, an analysis of transcriptome data from sepsis patients in the GEO database revealed a significant up-regulation of the phenylalanine-metabolism-related gene MAOA in the sepsis group. Inhibition of MAOA with the inhibitor RS-8359 resulted in improved mean arterial pressure and survival rate in septic rats. Finally, univariate and multivariate Cox regression analysis identified six DEMs (FA (22:5), L-Thyronic, Cer/AS (d17:3/16:2), Diethylphophase, Hydroxyco-tinine, and Syringol) that are associated with the prognosis of sepsis.

Machine learning algorithms have played a crucial role in the mining of metabolomic data [18]. A lot of new biomarkers can be found by analyzing metabolomics based on machine learning. Fortino et al. [19] utilized comprehensive transcriptome analysis coupled with machine learning to identify 89 biomarkers and decipher disease-related characteristic genes. In another study, Bi-farin et al. [20] employed machine learning algorithms to analyze urine metabolomics of renal cell carcinoma (RCC) patients and successfully identified seven metabolites for the diagnosis of RCC, achieving an impressive AUC of up to 0.98. Similarly, in a study involving serum samples from 109 gout patients, 102 asymptomatic hyperuricemia patients, and 119 normal uric acid control groups, three machine learning algorithms (random forest, support vector machine, and logistic regression) were employed, leading to the identification of 13 metabolites as potential biomarkers to distinguish hyperuricemia, gout, and normal uric acid conditions [21]. Given the substantial amount of metabolomics information in sepsis, the use of intelligent algorithms to mine variables becomes crucial. In this particular study, the SVM and RF algorithms were employed to screen for four distinctive DEMs in the serum of sepsis patients. Subsequently, the diagnostic efficacy of these four metabolites was evaluated using ROC curve analysis, ultimately identifying them as potential new markers for sepsis diagnosis.

Previous studies have shown that bacterial infection can lead to an imbalance of amino acid metabolism [22]. In sepsis patients, several metabolic characteristics have been observed, including hypermetabolism, catabolism, negative nitrogen balance, muscle and visceral protein decomposition for energy, liver uptake of amino acids for gluconeogenesis, and significant metabolic changes in various amino acids [23]. Amino acid metabolism disorder has been found to play a vital role in sepsis progression. Chen et al. conducted a non-targeted metabolomics analysis on serum samples from sepsis patients and identified extensive amino acid metabolism abnormalities. Notably, phenylalanine and histidine metabolism exhibited the most significant alterations in sepsis [24]. Aromatic amino acids (phenylalanine) can compete with branched-chain amino acids to penetrate the blood–brain barrier. The increase of aromatic amino acids causes the formation of pseudo-neurotransmitters, thus inhibiting the central nervous system [25]. Another study also indicated the association between phenylalanine metabolism and sepsis-related acute renal injury [26]. Amino acids like arginine and glycine possess antioxidant and immunomodulatory effects [27]. Glycine, in particular, can enhance myocardial function and diminish the production of free radicals, thereby mitigating endotoxin-induced myocardial damage [28]. Furthermore, by acting on the glycine receptor (GlyR) on cell membranes, glycine can prevent the excessive activation of Kupffer cells and play a significant role in liver protection [29]. Thus, glycine serves as a crucial metabolic regulator in human cells. Citrulline also plays a vital role in cell metabolism and organ function regulation [30]. This study not only observed abnormal amino acid metabolism, including phenyl-alanine metabolism, tyrosine metabolism, glycine, serine, and threonine metabolism, and arginine and proline metabolism in sepsis patients, but also examined the expression of genes related to phenylalanine metabolism in the sepsis patient transcriptome and the protective effect of MAOA inhibitors on sepsis in rats.

This study utilized metabolomics, machine learning algorithms, and prognosis analysis to screen potential metabolic biomarkers associated with the diagnosis and prognosis of sepsis. Moreover, it specifically examined the alterations in phenylalanine metabolism following sepsis and subsequently analyzed the regulatory mechanism of MAOA. However, since this study is considered exploratory and based on a small sample size of very new biomarkers, we opted to use completely healthy individuals as the control group. This decision could potentially lead to an overestimation of the sensitivity and specificity of the biomarkers. In future large sample studies, we will select case-controls as the control group to authenticate the diagnostic efficacy of these biomarkers.

Conclusions

The current study utilized metabolomics and machine learning algorithms to identify distinct metabolites associated with the diagnosis and prognosis of sepsis patients. Investigating the association between metabolic characteristics and sepsis offers novel insights into the underlying biological mechanisms, which may potentially enhance the diagnosis and treatment of sepsis.

Availability of data and materials

The datasets analyzed in this study are available from the corresponding authors upon reasonable request. Additionally, the datasets presented in this study can be found in an online repository, specifically at https://www.ncbi.nlm.nih.gov/geo/ with the accession number GSE57065.

References

  1. Bauer M, Gerlach H, Vogelmann T, Preissing F, Stiefel J, Adam D. Mortality in sepsis and septic shock in Europe, North America and Australia between 2009 and 2019- results from a systematic review and meta-analysis. Crit Care (London, England). 2020;24(1):239.

    Article  Google Scholar 

  2. Rudd KE, Kissoon N, Limmathurotsakul D, Bory S, Mutahunga B, Seymour CW, Angus DC, West TE. The global burden of sepsis: barriers and potential solutions. Crit Care (London, England). 2018;22(1):232.

    Article  Google Scholar 

  3. Phua J, Koay ES, Lee KH. Lactate, procalcitonin, and amino-terminal pro-B-type natriuretic peptide versus cytokine measurements and clinical severity scores for prognostication in septic shock. Shock (Augusta, Ga). 2008;29(3):328–33.

    Article  CAS  PubMed  Google Scholar 

  4. She H, Tan L, Du Y, Zhou Y, Guo N, Zhang J, Du Y, Wang Y, Wu Z, Ma C, et al. VDAC2 malonylation participates in sepsis-induced myocardial dysfunction via mitochondrial-related ferroptosis. Int J Biol Sci. 2023;19(10):3143–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Hamasaki MY, Barbeiro HV, de Souza HP, Machado MC, da Silva FP. sRAGE in septic shock: a potential biomarker of mortality. Rev Bras Ter Intensiva. 2014;26(4):392–6.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Pinu FR, Beale DJ, Paten AM, Kouremenos K, Swarup S, Schirra HJ, Wishart D. Systems biology and multi-omics integration: viewpoints from the metabolomics research community. Metabolites. 2019;9(4):76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. She H, Tan L, Wang Y, Du Y, Zhou Y, Zhang J, Du Y, Guo N, Wu Z, Li Q, et al. Integrative single-cell RNA sequencing and metabolomics decipher the imbalanced lipid-metabolism in maladaptive immune responses during sepsis. Front Immunol. 2023;14:1181697.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. She H, Tan L, Zhou Y, Zhu Y, Ma C, Wu Y, Du Y, Liu L, Hu Y, Mao Q, et al. The Landscape of featured metabolism-related genes and imbalanced immune cell subsets in sepsis. Front Genet. 2022;13:821275.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Mickiewicz B, Vogel HJ, Wong HR, Winston BW. Metabolomics as a novel approach for early diagnosis of pediatric septic shock and its mortality. Am J Respir Crit Care Med. 2013;187(9):967–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Schmerler D, Neugebauer S, Ludewig K, Bremer-Streck S, Brunkhorst FM, Kiehntopf M. Targeted metabolomics for discrimination of systemic inflammatory disorders in critically ill patients. J Lipid Res. 2012;53(7):1369–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Theodoridis GA, Gika HG, Want EJ, Wilson ID. Liquid chromatography-mass spectrometry based global metabolite profiling: a review. Anal Chim Acta. 2012;711:7–16.

    Article  CAS  PubMed  Google Scholar 

  12. Zhang P, Wang Z, Qiu H, Zhou W, Wang M, Cheng G. Machine learning applied to serum and cerebrospinal fluid metabolomes revealed altered arginine metabolism in neonatal sepsis with meningoencephalitis. Comput Struct Biotechnol J. 2021;19:3284–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Reel PS, Reel S, Pearson E, Trucco E, Jefferson E. Using machine learning approaches for multi-omics data analysis: A review. Biotechnol Adv. 2021;49:107739.

    Article  CAS  PubMed  Google Scholar 

  14. Hogan CA, Rajpurkar P, Sowrirajan H, Phillips NA, Le AT, Wu M, Garamani N, Sahoo MK, Wood ML, Huang C, et al. Nasopharyngeal metabolomics and machine learning approach for the diagnosis of influenza. EBioMedicine. 2021;71:103546.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Zhu Y, Zhang RH, Moum JN, Wang F, Li X, Li D. Physics-informed deep-learning parameterization of ocean vertical mixing improves climate simulations. Natl Sci Rev. 2022;9(8):nwac044.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Amiri M, Fisher PM, Raimondo F, Sidaros A, Hribljan MC, Othman MH, Zibrandtsen I, Albrechtsen SA, Bergdal O, Hansen AE et al: Multimodal prediction of residual consciousness in the intensive care unit: the CONNECT-ME study. Brain. 2022.

  17. Mukherjee S, Patra A, Khasawneh H, Korfiatis P, Rajamohan N, Suman G, Majumder S, Panda A, Johnson MP, Larson NB, et al. Radiomics-based machine-learning models can detect pancreatic cancer on Prediagnostic computed tomography scans at a substantial lead time before clinical diagnosis. Gastroenterology. 2022;163(5):1435-1446.e1433.

    Article  PubMed  Google Scholar 

  18. Liu Y, Wang Y, Ni Y, Cheung CKY, Lam KSL, Wang Y, Xia Z, Ye D, Guo J, Tse MA, et al. Gut Microbiome Fermentation Determines the Efficacy of Exercise for Diabetes Prevention. Cell Metab. 2020;31(1):77-91.e75.

    Article  CAS  PubMed  Google Scholar 

  19. Fortino V, Wisgrill L, Werner P, Suomela S, Linder N, Jalonen E, Suomalainen A, Marwah V, Kero M, Pesonen M, et al. Machine-learning-driven biomarker discovery for the discrimination between allergic and irritant contact dermatitis. Proc Natl Acad Sci USA. 2020;117(52):33474–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Bifarin OO, Gaul DA, Sah S, Arnold RS, Ogan K, Master VA, Roberts DL, Bergquist SH, Petros JA, Fernández FM, et al. Machine Learning-Enabled Renal Cell Carcinoma Status Prediction Using Multiplatform Urine-Based Metabolomics. J Proteome Res. 2021;20(7):3629–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Shen X, Wang C, Liang N, Liu Z, Li X, Zhu ZJ, Merriman TR, Dalbeth N, Terkeltaub R, Li C, et al. Serum Metabolomics Identifies Dysregulated Pathways and Potential Metabolic Biomarkers for Hyperuricemia and Gout. Arthritis Rheumatol (Hoboken, NJ). 2021;73(9):1738–48.

    Article  Google Scholar 

  22. Xiong L, Teng JL, Botelho MG, Lo RC, Lau SK, Woo PC. Arginine metabolism in bacterial pathogenesis and cancer therapy. Int J Mol Sci. 2016;17(3):363.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Seymour CW, Yende S, Scott MJ, Pribis J, Mohney RP, Bell LN, Chen YF, Zuckerbraun BS, Bigbee WL, Yealy DM, et al. Metabolomics in pneumonia and sepsis: an analysis of the GenIMS cohort study. Intensive Care Med. 2013;39(8):1423–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Chen Q, Liang X, Wu T, Jiang J, Jiang Y, Zhang S, Ruan Y, Zhang H, Zhang C, Chen P, et al. Integrative analysis of metabolomics and proteomics reveals amino acid metabolism disorder in sepsis. J Transl Med. 2022;20(1):123.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Basler T, Meier-Hellmann A, Bredle D, Reinhart K. Amino acid imbalance early in septic encephalopathy. Intensive Care Med. 2002;28(3):293–8.

    Article  PubMed  Google Scholar 

  26. Ping F, Guo Y, Cao Y, Shang J, Yao S, Zhang J, Li Y. Metabolomics analysis of the renal cortex in rats with acute kidney injury induced by sepsis. Front Mol Biosci. 2019;6:152.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Chen YL, Lin MT, Wang WH, Yeh SL, Yeh CL. Intravenous Arginine Administration Attenuates the Inflammatory Response and Improves Metabolic Profiles in Diet-Induced Obese Mice after Sleeve Gastrectomy. Metabolites. 2022;12(2):153.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Zhang Y, Ma X, Jiang D, Chen J, Jia H, Wu Z, Kim IH, Yang Y. Glycine Attenuates Lipopolysaccharide-Induced Acute Lung Injury by Regulating NLRP3 Inflammasome and NRF2 Signaling. Nutrients. 2020;12(3):611.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Al-Saeedi M, Nickkholgh A, Schultze D, Flechtenmacher C, Zorn M, Liang R, Gutt CN, Schemmer P. Glycine Protects the Liver from Reperfusion Injury following Pneumoperitoneum. Eur Surg Res. 2018;59(1–2):91–9.

    Article  CAS  PubMed  Google Scholar 

  30. Asgeirsson T, Zhang S, Nunoo R, Mascarenas C, Dujovny N, Luchtefeld M, Cavey GS, Senagore A. Citrulline: a potential immunomodulator in sepsis. Surgery. 2011;150(4):744–51.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the study participants.

Funding

This study was funded by the Natural Science Foundation of Chongqing (No. CSTB2023NSCQ-MSX0713), and National Natural Science Foundation of China (No. 82300561, No. 82270523, and No.81830065).

Author information

Author notes

  1. Han She, Yuanlin Du and Yunxia Du contributed equally to this work.

Authors and Affiliations

  1. Department of Anesthesiology, Daping Hospital, Army Medical University, Chongqing, 400042, China

    Han She, Yuanlin Du, Yunxia Du, Lei Tan, Shunxin Yang, Xi Luo & Yi Hu

  2. State Key Laboratory of Trauma, Burns and Combined Injury, Shock and Transfusion Department, Daping Hospital, Army Medical University, Chongqing, 400042, China

    Han She, Yuanlin Du, Yunxia Du, Lei Tan, Shunxin Yang, Xi Luo, Qinghui Li, Xinming Xiang, Liangming Liu & Tao Li

  3. Department of Intensive Care Unit, Daping Hospital, Army Medical University, Chongqing, 400042, China

    Haibin Lu

Authors
  1. Han She
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yuanlin Du
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yunxia Du
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Lei Tan
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Shunxin Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xi Luo
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Qinghui Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Xinming Xiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Haibin Lu
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Yi Hu
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Liangming Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Tao Li
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

TL, LML, and YH conceived and designed the study. HS, YLD, and LT analyzed the data. HS and YLD drafted the manuscript. TL, LML, and YH revised the manuscript. TL, LML, and HS acquired the financial support. All authors performed the experimental procedures.

Corresponding authors

Correspondence to Yi Hu, Liangming Liu or Tao Li.

Ethics declarations

Ethics approval and consent to participate

This human study was approved by the Ethics Committee of the Research Institute of Surgery (2021–179) and was registered by the Chinese Clinical Trial Registry (Date: 09/12/2021, ChiCTR2200055772). All the participants provided written informed consent before inclusion in this study. All experimental procedures were approved by the Animal Ethics Committee of Army Medical University, China, and were performed in accordance with the National Institutes of Health “Guidelines for the Care and Use of Laboratory Animals”. The study was carried out in compliance with the ARRIVE guidelines. All methods are reported in accordance with ARRIVE guidelines (https://arriveguidelines.org) for the reporting of animal experiments.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table 1.

The inclusion and exclusion criteria of septic patients.

Additional file 2.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

She, H., Du, Y., Du, Y. et al. Metabolomics and machine learning approaches for diagnostic and prognostic biomarkers screening in sepsis. BMC Anesthesiol 23, 367 (2023). https://doi.org/10.1186/s12871-023-02317-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12871-023-02317-4

Keywords

BMC Anesthesiology

ISSN: 1471-2253

Contact us