ORIGINAL RESEARCH article

Front. Endocrinol., 16 February 2023

Sec. Translational and Clinical Endocrinology

Volume 14 - 2023 | https://doi.org/10.3389/fendo.2023.1060470

Untargeted metabolomic profiling of sepsis-induced cardiac dysfunction

  • 1. Department of Emergency, Hunan Provincial People’s Hospital, The First Affiliated Hospital of Hunan Normal University, Changsha, Hunan, China

  • 2. Clinical Research Center for Emergency and Critical Care in Hunan Province, Hunan Provincial People’s Hospital, The First Affiliated Hospital of Hunan Normal University, Changsha, China

  • 3. Department of Cardiology, Hunan Provincial People’s Hospital, The First Affiliated Hospital of Hunan Normal University, Changsha, China

  • 4. Clinical Research Center for Heart Failure of Hunan Province, Hunan Provincial People’s Hospital, The First Affiliated Hospital of Hunan Normal University, Changsha, China

  • 5. Emergency Department, The First Affiliated Hospital of Hunan Normal University, Changsha, China

  • 6. Institute of Emergency Medicine, Hunan Provincial People’s Hospital, The First Affiliated Hospital of Hunan Normal University, Changsha, China

Abstract

Objective:

Sepsis is a life-threatening condition secondary to infection that evolves into a dysregulated host response and is associated with acute organ dysfunction. Sepsis-induced cardiac dysfunction is one of the most complex organ failures to characterize. This study performed comprehensive metabolomic profiling that distinguished between septic patients with and without cardiac dysfunction.

Method:

Plasma samples collected from 80 septic patients were analysed by untargeted liquid chromatography-mass spectrometry (LC-MS) metabolomics. Principal component analysis (PCA), partial least squares discrimination analysis (PLS-DA), and orthogonal partial least square discriminant analysis (OPLS-DA) were applied to analyse the metabolic model between septic patients with and without cardiac dysfunction. The screening criteria for potential candidate metabolites were as follows: variable importance in the projection (VIP) >1, P < 0.05, and fold change (FC) > 1.5 or < 0.7. Pathway enrichment analysis further revealed associated metabolic pathways. In addition, we constructed a subgroup metabolic analysis between the survivors and non-survivors according to 28-day mortality in the cardiac dysfunction group.

Results:

Two metabolite markers, kynurenic acid and gluconolactone, could distinguish the cardiac dysfunction group from the normal cardiac function group. Two metabolites, kynurenic acid and galactitol, could distinguish survivors and non-survivors in the subgroup analysis. Kynurenic acid is a common differential metabolite that could be used as a candidate for both diagnosis and prognosis for septic patients with cardiac dysfunction. The main associated pathways were amino acid metabolism, glucose metabolism and bile acid metabolism.

Conclusion:

Metabolomic technology could be a promising approach for identifying diagnostic and prognostic biomarkers of sepsis-induced cardiac dysfunction.

Introduction

In critically ill patients, sepsis remains a common condition that is associated with high mortality and substantial morbidity (). Cardiac dysfunction caused by sepsis, referred to as sepsis-induced cardiomyopathy, has long been a subject of interest because mortality can be greatly increased if sepsis is combined with this complication (). Thus, it is critical to assess the probability of septic cardiomyopathy and predict clinical outcome in the early stage of patient admission. However, the prognostic and therapeutic importance of physiological changes observed in sepsis-induced cardiac dysfunction remains poorly understood. In addition, although early application of ultrasonography is helpful for the diagnosis of septic cardiomyopathy, it lacks the basis for prognostic prediction. There is also no robust evidence to identify the association of troponin and brain natriuretic peptide (BNP) with septic cardiomyopathy (). Troponins, including troponin T and troponin I, are widely used markers for myocardial injury. The proportion of patients with elevated troponin T is up to 80% in sepsis (), but no evidence has been found to identify the relationship between left ventricular (LV) dysfunction and troponin T (). Another study reported that troponin I level was also not an independent predictor of mortality rate in sepsis (). BNP, secreted by cardiomyocytes of the ventricle, is dependent on LV filling pressures and myocardial wall stretch (). Julien et al. reported that BNP plasma levels are possibly useful for detecting myocardial dysfunction, but persistence of high BNP levels is associated with age and acute kidney injury (). These confounders should also be considered.

Metabolomics technology is an emerging omics science developed after genomics and proteomics. As an important part of systems biology, this technology qualitatively and quantitatively analyses a broad spectrum of small molecule metabolites, especially energy metabolites, in organisms (). The metabolome influences cellular physiology by regulating other levels of “omics”, including the genome, epigenome, transcriptome and proteome (). Recent advances in metabolomic technology have led to its increasing application in biomedicine. In particular, the application of metabolomics provides a strategic advantage for elucidating the mechanism of disease, discovering biomarkers, and innovating new therapeutics ().

Emerging studies highlight the potential application of metabolomics technology in sepsis (). Feng and colleagues reported that sepsis-induced acute kidney injury is accompanied by an increased oxygen consumption, systemic aerobic and anaerobic metabolism, and abnormal fatty acid metabolism (). A preliminary study demonstrated that 2-ethyl-2-hydroxybutyric acid regulates the expression of programmed cell death protein-1 on the surface of CD4+ T cells through the action of interleukin-2 or lactate, thereby affecting the prognosis of septic patients (). However, to our knowledge, no study has been reported to clarify metabolomics changes in sepsis-induced cardiac dysfunction.

In this study, nontargeted metabolomic profiling was utilized to investigate metabolomic alterations in patients with sepsis-induced cardiac dysfunction and to provide evidence for early metabolic biomarkers for diagnosis and prognosis. Moreover, related metabolic pathways were also explored.

Materials and methods

Study design and participants

Between November 2020 and March 2022, adults who had been admitted to the emergency department of Hunan Provincial People’s Hospital were selected for the study. Patients with sepsis were included in the study within the first 24 hours after admission. The inclusion criteria of patients were based on Sepsis 3.0 in the Third International Consensus Conference. Sepsis was defined as life-threatening organ dysfunction induced by a dysregulated host response to infection and a Sequential Organ Failure Assessment (SOFA) score ≥ 2 (). Sepsis-induced cardiac dysfunction was defined as impaired but reversible cardiac dysfunction under echocardiography, including LV systolic dysfunction, LV diastolic dysfunction, and right ventricle (RV) systolic dysfunction (). Patients who met the exclusion criteria were as follows: (1) age < 18 years old; (2) pregnancy; and (3) history of heart disease, such as acute coronary ischaemia, LV insufficiency, dilated cardiomyopathy, hypertrophic cardiomyopathy, valvular heart disease, or recurrent arrhythmia (). The study complied with the guidelines of the Declaration of Helsinki and the Conference for Coordination of Clinical Practice and was approved by the Ethics Committee of Hunan Provincial People’s Hospital. Informed consent was obtained from all participants.

Blood sampling

Fasting venous blood was obtained from patients diagnosed with sepsis within 24 hours of admission. Blood samples were collected in ethylenediaminetetraacetic acid (EDTA) anticoagulant tubes and centrifuged for 10 min at 4°C and 3,000 rpm to obtain plasma. The fasting plasma samples were stored in a -80°C freezer and kept frozen until metabolomic analysis.

Metabolite extraction and data processing

One hundred microlitres of the plasma sample was transferred to a 1.5 mL Eppendorf tube, and 20 μL of L-2-chlorophenylalanine (0.3 mg/mL) was dissolved in methanol as an internal standard. The tube was vortexed for 10 seconds. Subsequently, 300 μL of protein precipitant (methanol and acetonitrile, 2:1, vol/vol) was added, and the mixture was vortexed for 1 min. Then, the whole sample was sonicated for 10 min in an ice-water bath and stored at -20°C for 30 min. The extract was centrifuged at 4°C (13,000 rpm) for 10 min. The supernatant of each sample was collected. Quality control samples (QCs) were prepared by mixing equal volumes of extracts from all samples.

An ACQUITY UPLC I-Class system (Waters Corporation, Milford, USA) and a VION IMS QTOF mass spectrometer (Waters Corporation, Milford, USA) were used to analyse the metabolic profiles in both ESI positive and ESI negative ion modes. The target compounds were separated by an ACQUITY UPLC BEH C18 column (1.7 μm, 2.1 × 100 mm) at 45°C with 2 μL sample injection. Water and acetonitrile, both containing 0.1% formic acid, were used as mobile phases A and B, respectively. The gradient was set as follows: 0–2 min, 5% B; 4 min, 30% B; 8 min, 50% B; 10 min, 80% B; 14 min, 100% B; 15 min, 100% B; 15.1 min, 5% B and 16 min, 5% B. The flow rate was 0.35 mL/min, and the column temperature was 45 °C.

Primary and secondary mass spectrometry data were collected by a VION IMS QTOF mass spectrometer. The parameters of mass spectrometry were set as follows: a low-energy scan (CE 4 eV) and a high-energy scan (CE ramp 20–45 eV) to fragment the ions. Argon was used as the collision-induced dissociation gas; scan time: 0.2 s; interscan delay: 0.02 s; capillary voltage, 1 kV (negative mode) or 2 kV (positive mode); capillary voltage: 2.5 kV; cone voltage: 40 V; source temperature: 115°C; desolvation gas temperature: 450°C; and desolvation gas flow, 900 L/h.

The original data were analysed by Progenesis QI V2.3 (Nonlinear Dynamics, Newcastle, UK) software for baseline filtering, peak recognition, integration, retention time correction, peak alignment and normalization. Compound identification was based on the precise mass-to-charge ratio (m/z), MS2 fragments, and isotopic distribution using the Human Metabolome Database (HMDB) or Metlin, and self-built databases to do qualitative analysis. The metabolites in the self-built database are all established by standards, which contain retention time, first-level accurate mass information and second-level mass spectrum fragment information. The module pathway analysis was based on the KEGG database.

Statistical analysis

Clinical data were analysed using SPSS Statistics 25 software (IBM®, Armonk, NY, USA). Data are represented as the mean ± standard deviation or median and interquartile range. Continuous variables were compared using Student’s t test or the Mann−Whitney U test. Student’s t test is used when two samples are small and meet the conditions of normal distribution and homogeneity of variance. The Mann−Whitney U test was used when the samples did not meet the conditions of normal distribution and homogeneity of variance. Categorical variables between the two groups were compared by Fisher’s exact probability method. P < 0.05 was considered statistically significant.

The metabolic profiles were imported into R for principal component analysis (PCA) to observe the overall distribution among the samples and the stability of the entire analysis process. Partial least squares discriminant analysis (PLS-DA) and orthogonal partial least squares discriminant analysis (OPLS-DA) were used to distinguish differential metabolites between groups. To prevent overfitting, 7-fold cross-validation and 200 response permutation tests were utilized to evaluate the quality of the model. Variable importance of projection (VIP) values obtained from the OPLS-DA model were used to rank the overall contribution of each variable to group discrimination. A two-tailed Student’s t test was further used to verify whether the differences in metabolites between groups were significant. Differential metabolites were selected with VIP >1.0, P < 0.05, and fold change (FC) >1.5 or <0.7. Binary logistic regression analysis was constructed to screen independent risk factors. Receiver operating characteristic (ROC) curves were constructed to evaluate the diagnostic ability of differential metabolites between the tested groups.

Results

Study design and clinical synopsis

Initially, a total of 132 patients were enrolled. Based on the exclusion criteria, 37 patients with normal cardiac function and 43 patients with cardiac dysfunction were eventually included in our study. The flow chart is shown in Figure 1.

Figure 1

The clinical characteristics are shown in Table 1. Cardiac dysfunction is common in sepsis, affecting 54% (n=43) of patients. No significant differences in sex, age, comorbidities, or site of infection were noted between the two groups. Lactate levels and troponin I levels were higher in the cardiac dysfunction group (P=0.037 and P=0.023, respectively). N-terminal pro-brain natriuretic peptide (NT-proBNP) was elevated in the cardiac dysfunction group, but the difference was not statistically significant (P=0.051). The 28-day mortality was increased in the cardiac dysfunction group, although the difference was not statistically significant (P=0.178). ROC curves were constructed to determine the predictive value of variables in diagnosing sepsis-induced cardiac dysfunction. The AUC for lactate level (AUC 0.636, SE 95% confidence interval 0.513 to 0.759, P = 0.037) was equivalent to troponin I (AUC 0.647, SE 95% confidence interval 0.519 to 0.775, P = 0.024) and NT-proBNP (AUC 0.627, SE, 95% confidence interval 0.505 to 0.750, P = 0.051). The three AUCs were less than 0.75, indicating that the three variables were not robust enough for diagnosing sepsis-induced cardiac dysfunction.

Table 1

VariablesOverall (N=80)Septic patients with normal cardiac function (N=37)Septic patients with cardiac dysfunction (N=43)P value
Characteristic
Sex (Male%)47 (58.8)25 (67.6)22 (51.2)0.208
Age (years)59.2 ± 13.959.7 ± 13.658.7 ± 14.40.731
Combidities
Hypertension26 (32.5)12 (32.4)14 (32.6)0.990
Diabetes19 (23.8)8 (21.6)11 (25.6)0.678
Chronic respiratory disease20 (25.0)9 (24.3)11 (25.6)0.897
Chronic kidney disease8 (10.0)2 (5.4)6 (14.0)0.275
Neurologic disease7 (8.8)2 (5.4)5 (11.6)0.316
Site of infection
Lung44 (55.0)22 (59.5)22 (51.2)0.457
Gentiourinary tract30 (37.5)11 (29.7)19 (44.2)0.248
Hepatobiliary22 (27.5)12 (32.4)10 (23.3)0.359
Gastrointestine14 (17.5)6 (16.2)8 (18.6)0.779
Blood15 (18.8)5 (13.5)10 (23.3)0.390
Neural system2 (2.5)1 (2.7)1 (2.3)0.914
Soft tissue6 (7.5)3 (8.1)3 (7.0)0.848
Lactate(mmol/L)2.0 (1.3-3.6)1.7 (1.3-2.4)2.2 (1.6-4.9)0.037
Troponin I (ng/ml)0.07 (0.02-0.84)0.03 (0.01-0.93)0.10 (0.04-0.92)0.023
NT-proBNP (pg/ml)1405.0 (526.3-5972.5)925.0 (479.8-2596.8)2740.0 (621.0-6910.0)0.051
28-day mortality16 (20.0)5 (13.5)11 (25.6)0.178

Clinical characteristics of subjects enrolled in the study.

Data represented as mean ± standard deviation or median and interquartile range. Continuous variables were compared using Student’s t-test or the Mann-Whitney U test.Categorical variable were compared using the chi-square test or Fisher’s exact test as appropriate.Patients with sepsis was enrolled into the study within 24 hours from admission to emergency department.

Metabolite analysis and model validation in the normal cardiac function group and cardiac dysfunction group

The PCA model in the positive and negative ion modes showed that the instrument was stable throughout the experiment. The red squares, green triangles, and green circles in the figure represent the cardiac dysfunction group, normal cardiac function group, and QC group, respectively (Figures 2A, B).

Figure 2

Then, the PLS-DA method was applied to analyse the metabolite profile: in the positive ion mode, there were significant differences between the normal cardiac function group and the cardiac dysfunction group (CV-ANOVA,F=17.796, P<0.001]; in the negative ion mode, the two groups also exhibited significant differences (CV-ANOVA, F=19.066, P<0.001) (Figures 2C, D). These findings indicate that the PLS-DA model could be used to distinguish septic patients with cardiac dysfunction from those with normal cardiac function.

To achieve the greatest separation of differential metabolites between the two groups, OPLS-DA was performed based on values of VIP > 1. OPLS-DA demonstrated notable separation in metabolic profiles between the cardiac dysfunction group and the normal cardiac function group (Figures 2E, F). The permutation chart verified the validity of the model. R2 and Q2 were generated by permutation test. When R2 > 0, Q2 < 0, it indicates that the model is reliable (). The results of the permutation chart showed that the models were reliable (Figures 2G, H).

Screening of differential metabolites between the normal cardiac function group and cardiac dysfunction group

Univariate statistical analysis was applied to screen differential metabolites. A volcano plot was utilized to visualize the P value, VIP and fold change value. The red dots represent the upregulated metabolites, while the blue dots represent the downregulated metabolites (Figures 3A, C). The heatmap intuitively displays differential metabolites in different samples. In the heatmap of positive and negative ion modes, the Z scores of potential biomarkers were labelled in terms of the types of metabolites (Figures 3B, D). The most differential metabolites, which can be candidate for biomarkers, are presented in Table 2. The screening criteria for potential biomarkers were as follows: VIP >1, P < 0.05, and FC > 1.5 or < 0.7 (). Seventeen metabolites differed significantly between the normal cardiac function group and the cardiac dysfunction group.

Figure 3

Table 2

Ionization mode MetabolitesFormulaCompound ID
HMDB (or Metlin)
Normal cardiac function group vs.
cardiac dysfunction group
VIP  FCP value
PosL-AcetylcarnitineC9H17NO4HMDB000020111.3591.901<0.001
CreatineC4H9N3O275.1591.7080.022
5-Aminopentanoic acidC5H11NO2HMDB00033552.7400.591<0.001
N1-Methyl-4-pyridone-3-carboxamideC7H8N2O2HMDB00041941.6381.6890.042
3-Hydroxy-N6,N6,N6-
4-trimethyl-L-lysine
C9H20N2O362341.2271.7310.013
Quinolinic acidC7H5NO4HMDB00002321.1592.7200.017
5-MethoxyindoleacetateC11H11NO370161.0491.6250.010
25-HydroxycholesterolC27H46O2HMDB00062471.2200.6770.004
NegGluconic acidC6H12O73459.5903.6220.009
Glycochenodeoxycholic acid 3-glucuronideC32H51NO11HMDB00025794.7150.3670.033
3alpha-Androstanediol glucuronideC25H40O8HMDB02462524.2011.7480.016
D-Xylonic acidC5H4N4O3HMDB00597503.2331.899<0.001
Pregnanediol-3-
glucuronide
C27H44O8HMDB00103182.8212.0680.008
UridineC9H12N2O6HMDB00002962.7821.6770.001
Kynurenic acidC10H7NO3HMDB00007151.3894.6980.003
5’-Phosphoribosyl-N-formylglycinamideC8H15N2O9PHMDB00013081.3471.5830.018
GluconolactoneC6H10O6HMDB00001501.1802.299<0.001

In two ionization modes, differential metabolites between normal cardiac function group and cardiac dysfunction group.

Given VIP >1.0, P < 0.05, and FC >1.5 or < 0.7, 17 metabolites were identified.

The differential metabolites that enriched into pathways were listed in a decreasing order according to VIP.

Subsequently, we constructed binary logistic regression analysis to screen risk factors and ROC curve analysis to evaluate the predictive ability in the differential metabolites. Kynurenic acid, gluconolactone, 3-hydroxy-N6,N6,N6-trimethyl-L-lysine and 25-hydroxycholesterol were indentified as independent risk factors for sepsis-induced cardiac dysfunction, and the AUC values of kynurenic acid and gluconolactone were > 0.75 (Supplymentary Tables 1, 2). The AUCs for kynurenic acid and gluconolactone were 0.801, 0.754, respectively. We also constructed ROC analysis combining kynurenic acid and gluconolactone. However, the AUC for combining kynurenic acid and gluconolactone was 0.789, which was not superior to kynurenic acid (Figures 3E–G).

Differential metabolic pathways between the normal cardiac function group and cardiac dysfunction group

KEGG and HMDB were applied to analyse cardiac dysfunction-associated metabolites, and the results were submitted to MetaboAnalyst. In the top 20 metabolic pathway enrichment maps, significance is indicated by logarithm of the reciprocal of P value based on 10 (Figure 4A). In the top 20 bubble chart, the colour and size of each circle are determined by P values. (Figure 4B). The following three important metabolic pathways were identified: tryptophan metabolism, the pentose phosphate pathway, and lysine degradation. The results of pathway enrichment analysis for these differential metabolites are summarized in Supplementary Table 3.

Figure 4

Metabolite analysis and model validation in the survivors and non-survivors of cardiac dysfunction group

It has been reported that mortality can be increased if sepsis is combined with cardiac dysfunction (). In our study, 28-day mortality was increased in the cardiac dysfunction group, although the difference was not statistically significant due to the small sample size. It is crucial for physicians to predict clinical outcome exactly in the early stage of patient admission; thus, clarifying the metabolic changes in non-survivors is of great importance. A comparison between the survivors and non-survivors in the cardiac dysfunction group was performed. The PCA, PLS-DA and OPLS-DA models distinguished the metabolic profiles of the survivors from those of the non-survivors. The blue squares, red triangles and green circles in the figure represent the non-survivors, survivors and QC samples, respectively. PLS-DA analysis of the metabolite profile of the plasma sample: results in the positive ion mode suggested that there were significant differences between survivors and non-survivors (CV-ANOVA, F=23.301, P<0.001); in the negative ion mode, the two groups also exhibited significant differences (CV-ANOVA, F=26.217, P<0.001). In the OPLS-DA model, permutation tests showed that the model is reliable. (Figures 5A–H).

Figure 5

Screening of differential metabolites between the survivors and non-survivors of cardiac dysfunction group

Volcano maps and heatmaps were generated (Figures 6A–D). The most differential metabolites are presented in Table 3 following these criteria: VIP >1, P < 0.05, and FC > 1.5 or < 0.7. Twenty-five metabolites differed significantly between the survivors and non-survivors.

Figure 6

Table 3

Ionization mode Metabolites  FormulaCompound ID
HMDB (or Metlin)
Survivors vs.Non-survivors
VIPFCP value
PosL-AcetylcarnitineC9H17NO4HMDB00002019.6571.7120.014
Glycocholic acidC26H43NO6HMDB00001386.9432.5660.026
Taurochenodesoxycholic acidC12H22O111375.9938.9270.041
SucroseC7H8N2O2HMDB00041941.6384.9430.042
Taurocholic acidC26H45NO7S345425.5056.2740.029
CreatineC4H9N3O275.0021.9350.034
PhytosphingosineC18H39NO370663.5346.7730.023
SorbitolC6H14O6HMDB00002472.7402.7420.002
N-AcetylserotoninC12H14N2O2HMDB00012382.5426.847<0.001
UDP-D-galacturonic acidC15H22N2O18P2HMDB00123022.4493.1180.025
N1-Methyl-4-pyridone-3-carboxamideC7H8N2O2HMDB00041942.2841.8990.039
N1-(alpha-D-ribosyl)-5,6-dimethyl-benzimidazoleC14H18N2O4HMDB00111121.6514.0260.003
S-Adenosylmethioninam-ineC14H23N6O3SHMDB00009881.65012.2080.048
L-TryptophanC11H12N2O2331.6253.565<0.001
beta-D-ribosylnicotinateC11H13NO6HMDB03045401.4185.1410.017
L-HistidineC6H9N3O2211.3381.4010.044
S-AdenosylhomocysteineC14H20N6O5SHMDB00009391.0133.5150.006
NegPhenylacetylglutamineC13H16N2O4HMDB00063449.0235.6730.043
Hippuric acidC9H9NO3HMDB00007142.4693.3210.035
D-MaltoseC12H22O11HMDB00001632.1095.6190.017
UridineC9H12N2O6HMDB00002961.6461.5080.046
Salicylic acidC7H6O3HMDB00018951.4852.6140.034
GalactitolC6H14O6HMDB00001071.3103.6400.001
Kynurenic acidC10H7NO3HMDB00007151.2112.8500.022
Prostaglandin I2C20H32O5HMDB00013351.2061.7640.036

In two ionization modes, differential metabolites between survivors and non-survivors.

Given VIP >1.0, P < 0.05, and FC >1.5 or < 0.7, 25 metabolites were identified.

The differential metabolites that enriched into pathways were listed in a decreasing order according to VIP.

Then, we constructed binary logistic regression analysis to screen risk factors and ROC curve analysis to evaluate the predictive ability. Only kynurenic acid and galactitol were indentified as independent risk factors for 28-day mortality in sepsis-induced cardiac dysfunction and the AUC values were > 0.75 (Supplementary Tables 4, 5).The AUCs for kynurenic acid and galactitol were 0.821, 0.773, respectively. We also constructed ROC analysis combining kynurenic acid and galactitol. Whereas, the AUC for combining kynurenic acid and alactitol was 0.796, which was not superior to kynurenic acid (Figures 6E–G).

Differential metabolic pathways between survivors and non-survivors in the cardiac dysfunction group

In the top 20 metabolic pathway enrichment maps and bubble chart (Figures 7A, B), the following three important metabolic pathways were identified as differential metabolic pathways between the survivors and non-survivors: galactose metabolism, primary bile acid biosynthesis, and phenylalanine metabolism. The results of pathway enrichment analysis for these differential metabolites are summarized in Supplementary Table 6.

Figure 7

Discussion

Metabolomics is a rapidly developing technology that has been applied in many fields, such as biomarker investigation. To the best of our knowledge, this is the first study to perform metabolic analysis for septic patients with cardiac dysfunction. We first identified that key metabolites were changed during the development of sepsis-induced cardiac dysfunction. Two metabolite markers, kynurenic acid, and gluconolactone could distinguish the cardiac dysfunction group from the normal cardiac function group. Accordingly, we also discovered three important metabolic pathways related to sepsis-induced cardiac dysfunction, including tryptophan metabolism, pentose phosphate pathway, and lysine degradation. Next, we performed a comparison between the survivors and non-survivors in the cardiac dysfunction group. Two metabolites, including kynurenic acid and galactitol, could distinguish survivors and non-survivors in patients with sepsis-induced cardiac dysfunction. The crucial metabolic pathways are galactose metabolism, primary bile acid biosynthesis, and phenylalanine metabolism.

The most important finding was the identification of kynurenic acid as a potential metabolite biomarker for diagnosing sepsis-induced cardiac dysfunction and predicting poor outcomes. Kynurenic acid is a key degradation metabolite of tryptophan through the kynurenine pathway. Compared with kynurenine, kynurenic acid is not metabolized and stable in solution (). Hence, kynurenic acid is reported to be a potential candidate laboratory biomarker of metabolic diseases as atherosclerosis and diabetes (). Tryptophan is an essential amino acid that can only be provided from nutrition intake. Its metabolites play key roles in a variety of physiological processes, from cell growth and maintenance to coordinating the organism’s response to the environment (). Tryptophan catabolism by the kynurenine pathway is a pivotal metabolic pathway closely associated with the innate immune system (). In the early stage of sepsis, indoleamine 2,3-dioxygenase 1 (IDO1), a key enzyme that converts tryptophan to kynurenine, is activated and enhances proinflammatory effects that may result in severe tissue damage and septic shock (). Afterwards, the accumulation of kynurenine and its metabolites, such as kynurenic acid, can induce immunotolerance or even immunoparalysis through the expansion of regulatory T cells (Tregs) and T helper 2 (Th2) cells and with the mitigation of inflammatory Th1 and Th17-mediated responses (). Immunosuppression may contribute to the inability to overcome secondary infection and late mortality. These findings may explain why kynurenic acid levels were remarkably elevated in the cardiac dysfunction group and the non-survivors in this group. To our knowledge, this is the first study to find increased kynurenic acid in sepsis-induced cardiac injury. However, the underlying mechanism needs further investigation.

Our data also revealed that another important metabolite, gluconolactone, in the pentose phosphate pathway could be used to distinguish the cardiac dysfunction group and normal cardiac function group. Gluconolactone is considered to be a free radical scavenger. It relieves ischaemia/perfusion-induced cardiac injury although activation of extracellular signal-regulated kinase signalling (). However, the levels of gluconolactone were not greatly increased in the non-survivors of the cardiac dysfunction group. According to relevant literature and our results, we speculate that the alteration of gluconolactone might be a protective adaptive response to overcome the cardiotoxic effects of sepsis.

Another important finding is that the level of galactitol was markedly increased in the non-survivors of the cardiac dysfunction group. ROC analysis demonstrated that it is an independent risk factor to predict 28-day mortality in cardiac dysfunction patients. Galactitol, a downstream metabolite of galactose metabolism, is generated by the hydrogenation of galactose via aldose reductase and cannot be further metabolized (). A study also reported that galactitol emerged as the most obvious differential product in sepsis-related liver injury, suggesting that galactose metabolites are related to liver injury in the endotoxaemic state (). However, the underlying mechanisms of galactitol and sepsis-induced organ dysfunction need further discussion. Under physiological conditions, aldose reductase has relatively low affinity for galactose. In the case of galactose accumulation, the production of galactitol increases, and multiple-aspect damage is induced. Galactitol leads to the exhaustion of nicotinamide adenine dinucleotide phosphate and the reduction of glutathione reductase activity and then acts as a metabolic toxin to the body and induces the accumulation of free radicals (). Galactitol is incapable of diffusing across cellular membranes because of poor liposolubility, resulting in an increase in intracellular osmotic pressure (). The free radicals together with the osmotic effect of galactitol eventually result in rupture of the cellular membrane and impairment of mitochondrial DNA and proteins (). Further studies are needed to verify these potential mechanisms in sepsis.

Particular changes in the metabolic pathway occurred in this study. In the cardiac dysfunction group of septic patients, we found that the altered pathways were mainly involved in tryptophan metabolism, lysine degradation, and pentose phosphate pathway. In the non-survivors of the cardiac dysfunction group, we found that the changed pathways were mainly involved in galactose metabolism, primary bile acid biosynthesis and phenylalanine metabolism. Tryptophan metabolism, lysine degradation and phenylalanine metabolism belong to amino acid metabolism pathways. Both pentose phosphate pathway and galactose metabolism belong to the glucose metabolism pathways. Primary bile acid biosynthesis is associated with bile acid metabolism. Among these, tryptophan metabolism deserves special attention. Kynurenic acid, a candidate biomarker for diagnosing sepsis-induced cardiac dysfunction and predicting outcome, is an important metabolite of tryptophan metabolism. An increasing number of studies have revealed that tryptophan and its metabolites play key roles in inflammation-associated processes (). In addition to influencing T-cell immunity and leading to immune tolerance, emerging studies have highlighted the pivotal role of tryptophan metabolism in modulating B-cell functions and humorall immunity (). Aryl hydrocarbon receptor, a receptor that responds to tryptophan metabolites, affects the proliferation and switching of the immunoglobulin isotype in B cells (). It has been reported that B cells and neutrophils regulate each other in bone marrow, and B cells modulate neutrophils’ tissue-damaging properties by influencing neutrophils in sepsis (). Given the above, tryptophan metabolism has a major influence on clinical outcomes in critically ill septic patients.

Some limitations of our study should be considered. This study is a small, single-centre trial. Blood samples were analysed at only one time point. Further research with more participants will be performed to comprehensively evaluate the timing and dynamic changes in differential metabolites and metabolic pathways, providing a more precise view of changes during convalescence or deterioration. Our study also has some advantages. Currently, this is the first study to evaluate differential metabolites in septic patients with cardiac dysfunction. Patients were prospectively observed until 28 days after emergency admission. Mitochondrial function has a close relationship with cardiovascular disease. Future studies should include correlation analysis between the differential metabolites and mitochondrial bioenergetics and homeostasis to explore the underlying mechanism.

Conclusion

We used metabolomics technique to demonstrate that the metabolites of patients with sepsis-induced cardiac dysfunction change substantially and are mainly associated with the metabolic pathways of amino acid metabolism, glucose metabolism and bile acid metabolism. We also clearly distinguished septic patients with and without cardiac dysfunction using metabolites. As a consequence, kynurenic acid and gluconolactone are candidate biomarkers for diagnosing sepsis-induced cardiac dysfunction; kynurenic acid and galactitol are candidate biomarkers for predicting 28-day mortality. Kynurenic acid is a common differential metabolite in the two analyses and could be used as a diagnostic and prognostic biomarker. The differential metabolites and pathways may also be useful as targets for the development of new therapies for septic patients with cardiac dysfunction.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Hunan Provincial People’s Hospital. The patients/participants provided their written informed consent to participate in this study.

Author contributions

Conceptualization, YC, WX and XH. Methodology, YC and ZL. Formal analysis, YC and ZL. Investigation and data curation, YJ and WX. Writing—original draft preparation, YC, YP and WM. Writing—review and editing, CF and JH. Funding acquisition, YC, XH. WX and XH are the master investigators of the manuscript, take responsibility for guiding the whole work, from start to publication article. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by Clinical Medical Technology Innovation Guidance Program of Hunan Province [2020SK50909], Traditional Chinese Medicine Research Program of Hunan Province [C2023020], Natural Science Foundation of Hunan Province [2020JJ8082], Program of Clinical Research Center for Emergency and Critical Care in Hunan Province [2021SK4011], Major Program of Hunan Province Science and Technology Innovation [2020SK1010]

Acknowledgments

All the authors would like to express appreciation to all the patients for their cooperation. Furthermore, our thanks would go to the foundation from Clinical Medical Technology Innovation Guidance Program of Hunan Province [2020SK50909], Traditional Chinese Medicine Research Program of Hunan Province [C2023020], Natural Science Foundation of Hunan Province [2020JJ8082], Program of Clinical Research Center for Emergency and Critical Care in Hunan Province [2021SK4011], Major Program of Hunan Province Science and Technology Innovation [2020SK1010]

Conflict of interest

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

Publisher’s note

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

Supplementary material

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

References

  • 1

    CecconiMEvansLLevyMRhodesA. Sepsis and septic shock. LANCET (2018) 392(10141):7587. doi: 10.1016/S0140-6736(18)30696-2

  • 2

    BeesleySJWeberGSargeTNikravanSGrissomCKLanspaMJet al. Septic cardiomyopathy. Crit Care Med (2018) 46(4):625–34. doi: 10.1097/CCM.0000000000002851

  • 3

    EhrmanRRSullivanANFavotMJSherwinRLReynoldsCAAbidovAet al. Pathophysiology, echocardiographic evaluation, biomarker findings, and prognostic implications of septic cardiomyopathy: A review of the literature. Crit Care (2018) 22(1):112. doi: 10.1186/s13054-018-2043-8

  • 4

    VasileVCChaiH-SAbdeldayemDAfessaBJaffeAS. Elevated cardiac troponin T levels in critically ill patients with sepsis. Am J OF Med (2013) 126(12):1114–21. doi: 10.1016/j.amjmed.2013.06.029

  • 5

    PulidoJNAfessaBMasakiMYuasaTGillespieSHerasevichVet al. Clinical spectrum, frequency, and significance of myocardial dysfunction in severe sepsis and septic shock. MAYO CLINIC Proc (2012) 87(7):620–8. doi: 10.1016/j.mayocp.2012.01.018

  • 6

    InnocentiFPalmieriVStefanoneVTD'ArgenzioFCiganaMMontuoriMet al. Comparison of troponin I levels versus myocardial dysfunction on prognosis in sepsis. Internal Emergency Med (2022) 17(1):223–31. doi: 10.1007/s11739-021-02701-3

  • 7

    HollenbergSMSingerM. Pathophysiology of sepsis-induced cardiomyopathy. Nat Rev Cardiol (2021) 18(6):424–34. doi: 10.1038/s41569-020-00492-2

  • 8

    CharpentierJLuytC-EFullaYVinsonneauCCariouAGrabarSet al. Brain natriuretic peptide: A marker of myocardial dysfunction and prognosis during severe sepsis. Crit Care Med (2004) 32(3):660–5. doi: 10.1097/01.ccm.0000114827.93410.d8

  • 9

    GuijasCMontenegro-BurkeJRWarthBSpilkerMESiuzdakG. Metabolomics activity screening for identifying metabolites that modulate phenotype. Nat Biotechnol (2018) 36(4):316–20. doi: 10.1038/nbt.4101

  • 10

    RinschenMMIvanisevicJGieraMSiuzdakG. Identification of bioactive metabolites using activity metabolomics. Nat Rev Mol Cell Biol (2019) 20(6):353–67. doi: 10.1038/s41580-019-0108-4

  • 11

    van der HooftJJJMohimaniHBauermeisterADorresteinPCDuncanKRMedemaMH. Linking genomics and metabolomics to chart specialized metabolic diversity. Chem Soc Rev (2020) 49(11):3297–314. doi: 10.1039/d0cs00162g

  • 12

    TrongtrakulKThonusinCPothiratCChattipakornSCChattipakornN. Past experiences for future applications of metabolomics in critically ill patients with sepsis and septic shocks. Metabolites (2021) 12(1):1–32. doi: 10.3390/metabo12010001

  • 13

    PingFLiYCaoYShangJZhangZYuanZet al. Metabolomics analysis of the development of sepsis and potential biomarkers of sepsis-induced acute kidney injury. Oxid Med Cell Longev (2021) 2021:6628847. doi: 10.1155/2021/6628847

  • 14

    BuYWangHMaXHanCJiaXZhangJet al. Untargeted metabolomic profiling of the correlation between prognosis differences and PD-1 expression in sepsis: A preliminary study. Front Immunol (2021) 12:594270. doi: 10.3389/fimmu.2021.594270

  • 15

    SingerMDeutschmanCSSeymourCWShankar-HariMAnnaneDBauerMet al. The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA (2016) 315(8):801–10. doi: 10.1001/jama.2016.0287

  • 16

    HeureuxMSternbergMBrathLTurlingtonJKashiourisMG. Sepsis-induced cardiomyopathy: A comprehensive review. Curr Cardiol Rep (2020) 22(5):35. doi: 10.1007/s11886-020-01277-2.

  • 17

    DuYFanPZouLJiangYGuXYuJet al. Serum metabolomics study of papillary thyroid carcinoma based on HPLC-Q-TOF-MS/MS. Front Cell Dev Biol (2021) 9:593510. doi: 10.3389/fcell.2021.593510

  • 18

    CervenkaIAgudeloLZRuasJL. Kynurenines: Tryptophan's metabolites in exercise, inflammation, and mental health. SCIENCE (2017) 357(6349):eaaf9794. doi: 10.1126/science.aaf9794

  • 19

    ZhenDLiuJZhangXDSongZ. Kynurenic acid acts as a signaling molecule regulating energy expenditure and is closely associated with metabolic diseases. Front Endocrinol (2022) 13:847611. doi: 10.3389/fendo.2022.847611

  • 20

    DehhaghiMKazemiSPHGuilleminG. Microorganisms, tryptophan metabolism, and kynurenine pathway: A complex interconnected loop influencing human health status. Int J tryptophan Res IJTR (2019) 12:178646919852996. doi: 10.1177/1178646919852996

  • 21

    KrupaAKowalskaI. The kynurenine pathway-new linkage between innate and adaptive immunity in autoimmune endocrinopathies. Int J Mol Sci (2021) 22(18):9879. doi: 10.3390/ijms22189879

  • 22

    BadawyAA-B. Kynurenine pathway of tryptophan metabolism: Regulatory and functional aspects. Int J tryptophan Res IJTR (2017) 10:1178646917691938. doi: 10.1177/1178646917691938

  • 23

    KrupaAKrupaMMPawlakK. Kynurenine pathway-an underestimated factor modulating innate immunity in sepsis-induced acute kidney injury? Cells (2022) 11(16):2604. doi: 10.3390/cells11162604

  • 24

    QinXLiuBGaoFHuYChenZXuJet al. Gluconolactone alleviates myocardial Ischemia/Reperfusion injury and arrhythmias via activating PKCϵ/Extracellular signal-regulated kinase signaling. Front Physiol (2022) 13:856699. doi: 10.3389/fphys.2022.856699

  • 25

    ConteFvan BuuringenNVoermansNCLefeberDJ. Galactose in human metabolism, glycosylation and congenital metabolic diseases: Time for a closer look. Biochim ET Biophys ACTA-GENERAL Subj (2021) 1865(8):129898. doi: 10.1016/j.bbagen.2021.129898

  • 26

    LiTBaiJDuYTanPZhengTChenYet al. Thiamine pretreatment improves endotoxemia-related liver injury and cholestatic complications by regulating galactose metabolism and inhibiting macrophage activation. Int Immunopharmacol (2022) 108:108892. doi: 10.1016/j.intimp.2022.108892

  • 27

    Coelho AnaIBerryGTRubio-GozalboME. Galactose metabolism and health. Curr Opin IN Clin Nutr AND Metab Care (2015) 18(4):422–7. doi: 10.1097/MCO.0000000000000189

  • 28

    Bo-HtayCPaleeSApaijaiNChattipakornSCChattipakornN. Effects of d-galactose-induced ageing on the heart and its potential interventions. J OF Cell AND Mol Med (2018) 22(3):1392–410. doi: 10.1111/jcmm.13472

  • 29

    KanovaMKohoutP. Tryptophan: A unique role in the critically ill. Int J Mol Sci (2021) 22(21):11714. doi: 10.3390/ijms222111714

  • 30

    Dagenais-LussierXLoucifHBejiCTelittchenkoRRoutyJPvan GrevenyngheJ. Latest developments in tryptophan metabolism: Understanding its role in b cell immunity. Cytokine Growth FACTOR Rev (2021) 59:111–7. doi: 10.1016/j.cytogfr.2021.02.003

  • 31

    VaidyanathanBChaudhryAYewdellWTAngelettiDYenWFWheatleyAKet al. The aryl hydrocarbon receptor controls cell-fate decisions in b cells. J OF Exp Med (2017) 214(1):197208. doi: 10.1084/jem.20160789

  • 32

    GawishRMaierBObermayerGWatzenboeckMLGorkiADQuattroneFet al. A neutrophil-b-cell axis impacts tissue damage control in a mouse model of intra-abdominal bacterial infection via Cxcr4. eLife (2022) 11:e78291. doi: 10.7554/eLife.78291

Summary

Keywords

sepsis, cardiac dysfunction, metabolomics, diagnosis, prognosis

Citation

Cao Y, Liu Z, Ma W, Fang C, Pei Y, Jing Y, Huang J, Han X and Xiao W (2023) Untargeted metabolomic profiling of sepsis-induced cardiac dysfunction. Front. Endocrinol. 14:1060470. doi: 10.3389/fendo.2023.1060470

Received

03 October 2022

Accepted

23 January 2023

Published

16 February 2023

Volume

14 - 2023

Edited by

Yao Lu, Third Xiangya Hospital, Central South University, China

Reviewed by

Yaping Shao, Dalian Medical University, China; Weidong Dai, Tea Research Institute (CAAS), China

Updates

Copyright

*Correspondence: Xiaotong Han, ; Weiwei Xiao,

†These authors have contributed equally to this work and share first authorship

This article was submitted to Translational Endocrinology, a section of the journal Frontiers in Endocrinology

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