ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 10 July 2026

Sec. Intestinal Microbiome

Volume 16 - 2026 | https://doi.org/10.3389/fcimb.2026.1811654

Gut microbiota and serum metabolic profiles in patients with sepsis-induced cardiomyopathy and their association with the disease

  • 1. First Clinical Medical College, Ningxia Medical University, Yinchuan, China

  • 2. Department of Critical Care Medicine, General Hospital of Ningxia Medical University, Yinchuan, China

Abstract

Background:

The pathogenesis of sepsis-induced cardiomyopathy (SIC) remains unclear, and the lack of effective treatment options results in high mortality rates. There may be interactions between gut microbiota dysbiosis, serum metabolic changes, and SIC; however, research in this area is limited. This study aims to provide a theoretical basis for targeted interventions to improve the gut microbiome, thereby enhancing cardiac function and prognosis in patients with SIC.

Methods:

This prospective cohort study enrolled 48 septic patients admitted to the intensive care unit (ICU) of a tertiary general hospital in Ningxia, China, between March 2024 and March 2025. Patients were stratified into SIC (n=28) and SEPSIS (n=20) groups. Age- and sex-matched healthy controls (HC, n=10) were also recruited. Fecal samples from all three groups were analyzed via 16S rRNA gene amplicon sequencing, and serum samples from the SIC and SEPSIS groups were analyzed via LC–MS-based untargeted metabolomics. Spearman’s correlations linked microbial and metabolic alterations to echocardiographic parameters and myocardial injury markers.

Results:

Compared with the SEPSIS and HC groups, the SIC group presented significant differences in both the α and β diversity of the gut microbiota (all P<0.05). Compared with that in the SEPSIS group, the Bacillota/Bacteroidota ratio was significantly lower in the SIC group (P<0.05); compared with the SEPSIS and HC groups, the SIC group presented significantly increased abundances of genera such as Parabacteroides and Clostridium innocuum group and significantly decreased abundances of genera such as Oscillospiraceae_UCG-003 and Lachnospiraceae_UCG-004 (all P<0.05). Forty-nine differential metabolites were identified (VIP≥1, P<0.05): hippuric acid, estrone glucuronide, and TMAO were elevated, while 15-deoxy-Δ12,14-prostaglandin J2 was reduced in the SIC group (all P<0.05). Differential genera/metabolites correlated strongly with SIC diagnostic indicators: e.g., Parabacteroides positively correlated with estrone glucuronide, hippuric acid, and GLS, and negatively with LVEF.

Conclusion:

Gut microbiota dysbiosis differs between patients with Sepsis-induced cardiomyopathy and healthy individuals or patients with sepsis. Additionally, the serum metabolic profile changes in SIC patients differ from those in patients with sepsis, and these characteristic changes are significantly associated with clinical indicators related to the diagnosis of SIC.

Introduction

Sepsis is a life-threatening organ dysfunction caused by a dysregulated host immune response to infection (), and the heart is one of the most vulnerable target organs involved in sepsis (). Sepsis-induced cardiac dysfunction is defined as sepsis-induced cardiomyopathy (SIC), which is characterized by diffuse and reversible myocardial dysfunction. SIC complicates 10%–70% of sepsis patients, with a mortality rate of up to 70%–90%, severely impairing the long-term prognosis of septic patients (; ). The underlying mechanism of SIC is complex and not fully understood; the absence of specific targeted therapies hinders prognostic improvement (), making SIC a research hotspot in critical care medicine.

As the largest reservoir of intestinal bacteria and endotoxins in the human body, the gut is recognized as the “engine” of sepsis-associated multiple organ dysfunction syndrome (MODS) (). Accumulating evidence has demonstrated that gut microbiota dysbiosis plays a pivotal role in the pathogenesis of sepsis-related MODS (; ). Whether SIC is accompanied by characteristic gut flora dysbiosis remains to be clarified. Previous studies have confirmed that the structure and composition of the gut microbiota differ significantly between SIC patients and non-SIC sepsis patients (). In addition to bacterial flora, intestinal viral and fungal communities also undergo structural remodeling in SIC patients (). Previous studies have verified the interactions among sepsis, sepsis-related MODS and gut microbiota dysbiosis. Given that the heart is a major target organ of sepsis injury, the unique gut microbial characteristics of SIC patients (compared with those of healthy individuals and simple sepsis patients) may reveal the dynamic correlation between disease progression (from healthy individuals to simple sepsis patients and further to SIC patients) and gut dysbiosis, providing a reference for interventions in the gut microbiota to block disease progression and optimize the clinical prognosis.

Gut microbiota dysbiosis is closely implicated in the pathogenesis of multiple cardiovascular diseases. It participates in the pathological progression of hypertension, atherosclerosis, heart failure and atrial fibrillation by increasing intestinal permeability, promoting lipopolysaccharide (LPS) translocation and triggering chronic systemic inflammation (; ; ; ). Moreover, gut microbial metabolites exert regulatory effects on cardiac function. Short-chain fatty acids (SCFAs) exert cardioprotective effects through anti-inflammatory activity and intestinal barrier maintenance (; ). Elevated circulating TMAO and abnormal bile acid metabolism accelerate the progression of multiple cardiovascular diseases (; ). Therefore, exploring the characteristics of gut microbiota dysbiosis, serum metabolic alterations and their disease relevance in SIC can deepen the understanding of SIC pathogenesis and provide a novel strategy for improving SIC prognosis via gut microbiome regulation, which has important clinical translational value.

On the basis of the above background, this study adopted 16S rRNA gene amplicon sequencing to compare the gut microbial composition among SIC patients, non-SIC sepsis patients and healthy controls. Moreover, serum untargeted metabolomics via LC–MS was used to characterize metabolic profile changes in SIC patients. Furthermore, we conducted correlation analyses of differential gut bacterial genera, serum metabolites, echocardiographic parameters and myocardial injury biomarkers to elucidate the intrinsic associations among gut dysbiosis, serum metabolic disorders and myocardial damage in SIC patients.

Study design and methods

Study design and participants

This was a prospective observational cohort study recruiting sepsis patients hospitalized in the ICU of the General Hospital of Ningxia Medical University between March 2024 and March 2025.

The inclusion criteria were as follows: ① met the Sepsis 3.0 diagnostic criteria () and ②were aged ≥18 years. Exclusion criteria: ① Complicated with acute coronary syndrome or receiving cardiac interventional surgery during hospitalization; ② history of chronic heart failure, long-term dialysis-dependent chronic renal insufficiency, or Child–Pugh grade C chronic liver disease; ③ preexisting hypertrophic, dilated or restrictive cardiomyopathy, as well as rheumatic heart disease; ④ history of life-threatening severe arrhythmias (e.g., ventricular tachycardia, ventricular fibrillation); ⑤ unavailability of bedside echocardiography during the ICU; and ⑥ refusal of enrollment or voluntary withdrawal by patients or their legal guardians.

This study was approved by the Ethics Committee of the General Hospital of Ningxia Medical University (Approval No. KYLL-2025-0038) and was registered in the Chinese Clinical Trial Registry (Registration No. CTR2500097664). All participants or their guardians signed written informed consent prior to enrollment.

Sample size calculation was performed via beta diversity PERMANOVA (α=0.05, statistical power=0.8, 1:1 matching design) and the method proposed by . With an expected effect size R2 = 0.3, the estimated sample size per group was 18–20 cases. Owing to limited clinical case enrollment and strict inclusion/exclusion criteria, the final enrollment met the statistical requirements. No prior sample size calculation was performed for metabolomics analysis because of the lack of reference data in the same population.

The diagnostic criteria for SIC are as follows (; ; ): ① acute onset with potential functional reversal within 7–10 days; ② global biventricular systolic and/or diastolic dysfunction accompanied by decreased myocardial contractility; ③ left ventricular dilatation; ④ poor hemodynamic response to fluid resuscitation and catecholamine therapy; and ⑤ exclusion of acute coronary syndrome as the primary cause of cardiac dysfunction.

Serial echocardiography was performed dynamically to evaluate changes in cardiac function. Professional sonographers completed echocardiographic examinations on days 1, 3 and 5 (if necessary) after ICU admission. Two senior critical care physicians jointly made the SIC diagnosis on the basis of the patients’ clinical history and dynamic echocardiographic parameters.

The enrolled sepsis patients were divided into an SIC group (n=28) and a non-SIC sepsis group (n=20). The HC group (n=10) consisted of age- and sex-matched healthy volunteers with normal echocardiographic results and no underlying cardiovascular, metabolic or gastrointestinal diseases.

Data collection

Echocardiographic parameters: A GE Vivid E9 ultrasound system was used for cardiac ultrasound examination. The detected parameters included the following: left ventricular global longitudinal strain (GLS), left ventricular ejection fraction (LVEF), and left ventricular outflow tract time–velocity integral (VTI); diastolic function indicators: mitral early diastolic flow velocity (E), mitral late diastolic flow velocity (A), E/A ratio, mitral annular early diastolic tissue velocity (e’), E/e’ ratio, left atrial volume index (LAVI), and peak tricuspid regurgitation velocity (TRV); and right ventricular systolic function indicators: tricuspid annular plane systolic excursion (TAPSE) and tricuspid systolic flow velocity (S). All ultrasound images were acquired by dedicated cardiac sonographers and independently reviewed by two senior echocardiographers in strict accordance with the American Society of Echocardiography guidelines ().

Clinical baseline data: Demographic characteristics, underlying comorbidities, infection sites, and APACHE II and SOFA scores were recorded. Within 24 hours after sepsis diagnosis, laboratory indicators, including the mean arterial pressure (MAP), white blood cell (WBC) count, lymphocyte percentage (LY%), platelet count (PLT), C-reactive protein (CRP), procalcitonin (PCT), blood lactate (Lac), oxygenation index (OI), serum creatinine (CREA), total bilirubin (TBIL), cardiac troponin I (cTnI) and N-terminal pro-B-type natriuretic peptide (NT-proBNP), were collected.

Fecal sample collection and 16S rRNA gene sequencing: Fecal samples were collected upon ICU admission synchronously with echocardiography. For SIC patients, samples within 24 hours after confirmed SIC diagnosis were used for detection; for simple sepsis patients, the first fecal sample after enrollment was analyzed. Genomic DNA was extracted via an E.Z.N.A.® Stool DNA Kit (Omega Biotek, USA) following the manufacturer’s instructions. The DNA concentration and purity were determined via a QUBIT 3.0 fluorometer (Invitrogen, USA). The V3-V4 hypervariable region of the 16S rRNA gene was amplified with the primers 341F (5’-CCTACGGGNGGCWGCAG-3’) and 805R (5’-GACTACHVGGGTATCTAATCC-3’). PCR amplification was performed on an Easy Cycler 96 PCR instrument (Analytik Jena, Germany). The purified PCR products were mixed equally and sequenced on an Illumina NextSeq 2000 platform.

Serum sample collection and LC–MS untargeted metabolomics: Peripheral venous blood samples were collected synchronously with echocardiography upon ICU admission. Serum was isolated by centrifugation at 3000 rpm for 10 min and stored at -80 °C for subsequent detection. Serum samples were thawed slowly at 4 °C, mixed with precooled methanol:acetonitrile:water mixture (2:2:1), vortexed, subjected to low-temperature ultrasonic extraction for 30 min, and then incubated at -20 °C for 10 min. After centrifugation at 14000 × g for 20 min at 4 °C, the supernatant was collected and vacuum-dried. For mass spectrometry detection, the dried sample was redissolved in 100 μL of acetonitrile-water mixture (1:1), vortexed and centrifuged at 14000 × g for 15 min at 4 °C, and the resulting mixture was injected into the sample for analysis. Chromatographic separation was performed on a Vanquish ultrahigh-performance liquid chromatography system coupled with an Orbitrap Exploris™ 480 mass spectrometer (Thermo Fisher, USA). Electrospray ionization was applied in both positive and negative ion modes.

Data quality control

16S rRNA sequencing data: Raw paired-end sequencing data were sorted by sample and sequencing batch. Cutadapt was used to trim primer sequences and filter abnormal fragments. FastQC and MultiQC were applied to evaluate sequencing quality and base distribution. The quality-filtered sequences were imported into the QIIME 2 (v2024.2) platform and processed via the DADA2 algorithm, including quality trimming, error correction, denoising, paired-end read assembly and chimera removal. Only sequences with overlapping lengths ≥20 bp were retained to generate amplicon sequence variants (ASVs) and abundance profiles. Species annotation was performed via the SILVA 138.2 database, which is based on the naive Bayes classifier. Low-abundance ASVs (sequence count <5 or detected in only one sample) and off-target sequences (mitochondria, chloroplasts, archaea) were removed. Samples with total sequencing reads <4000 were excluded. The average sequencing depth of all the samples was 38528 reads. All samples were normalized to the minimum sequencing depth of 19654 reads for subsequent alpha/beta diversity analysis.

Serum metabolomics data: Raw mass spectrometry data were converted to mzXML format by ProteoWizard. XCMS software was used for peak alignment, retention time correction and peak area extraction. Ion peaks with missing value rates >50% were removed, and missing values were imputed via the K-nearest neighbor algorithm. The feature peaks with a relative standard deviation (RSD) >50% were filtered out. Quality control (QC) sample total ion chromatogram overlap analysis revealed the high stability of instrument detection. Principal component analysis confirmed good clustering repeatability of the QC samples. When the proportion of ion peaks with an RSD ≤30% exceeded 70%, the data were considered suitable for subsequent analysis. Metabolite annotation was performed by matching the exact mass–charge ratio (mass error ≤5 ppm), characteristic fragment ions and retention time with the HMDB and KEGG databases. Metabolite identification confidence strictly followed the four-level criteria of the Metabolomics Standards Initiative (MSI). Targeted MS/MS verification with standard compounds was conducted for key differentially abundant metabolites ().

Bioinformatics and statistical analysis

Bioinformatics analysis: Microbial sequencing data were processed via the DADA2 algorithm to obtain ASV abundance tables. Alpha diversity was calculated to evaluate gut microbial richness and diversity. Principal coordinate analysis (PCoA) based on Bray–Curtis, unweighted and weighted UniFrac distances was used to assess intergroup microbial structural differences. PERMANOVA (Adonis) analysis was used to test the statistical significance of community differences. Intergroup differences in microbial relative abundance were analyzed via the Mann–Whitney U test (two groups) and the Kruskal–Wallis H test (three groups). Linear discriminant analysis effect size (LEfSe) was used to screen for biomarker bacterial genera among the groups.

Serum LC/MS non-targeted metabolomics analysis: For the metabolomics data, PCA and orthogonal partial least squares discriminant analysis (OPLS-DA) were performed after total peak area normalization and UV scaling. OPLS-DA model validation was conducted via 200 bootstrap permutations. Differentially abundant metabolites with VIP>1 and P<0.05 (Student’s t test) were identified. KEGG pathway enrichment analysis was performed via the MetaboAnalyst online platform.

Correlation analysis: Spearman’s correlation analysis was used to analyze the correlation of differential bacterial genera and metabolites with echocardiographic parameters and myocardial injury biomarkers. CCA analysis was performed using the `cca` function in the `vegan` package of R. Key differential metabolite profiles were integrated with clinical and physicochemical indicators, such as those associated with SIC diagnosis, to construct a combined metabolic-clinical-environmental factor matrix; microbial community composition at the genus level served as the species matrix. All environmental variables were standardized using UV normalization (Standard Z-score normalization) to eliminate dimensional differences and numerical weighting biases. This approach revealed covariation patterns between microbial data and host multi-omics indicators from a holistic perspective.

Statistical methods: All the statistical analyses were performed via R (v4.5.1) and SPSS 27.0. The Shapiro–Wilk test was used for normality detection. Normally distributed data are presented as the means ± standard deviations and were compared via independent sample t tests; nonnormally distributed data are presented as medians (interquartile ranges) and were analyzed via nonparametric tests. Categorical variables are expressed as case numbers (percentages) and were compared via the chi-square test or Fisher’s exact test. A two-tailed P<0.05 was considered statistically significant.

Results

Patient characteristics

This study included a total of 48 patients with sepsis, comprising the SIC group (n=28, 58%) and the SEPSIS group (n=20, 42%). The mean age of the enrolled patients was 54.87 ± 14.66 years, with 30 (62.5%) being male. There were no statistically significant differences between the two groups in terms of sex, age, body surface area, comorbidities, site of infection, APACHE II score, or SOFA score (all P>0.05) (Table 1). The HCs (n=10) had a mean age of 55 ± 4.82 years, were male (n=6, 60%), and were fully matched with the enrolled patients. To avoid the influence of confounding factors, intergroup comparisons were performed for several confounding variables; there were no statistically significant differences between the two groups in BMI, endotracheal intubation, or types of antibiotics used (all P>0.05, Supplementary Table 8).

Table 1

VariableHC (n=10)SEPSIS (n=20)SIC (n=28)F/χ2/Z/tP
Age (years, ± s)55.00 ± 4.8254.80± 15.7055.00 ± 14.2000.88b
Gender0.040.76b
 Male n (%)6 (60%)12 (60%)18 (64%)
 Female n (%)4 (40%)8 (40%)10 (36%)
BMI (kg/m2, ± s)22.79 ± 2.3224.44 ± 4.253.790.09a
Body surface area (m2, ± s)1.60 ± 0.201.70 ± 0.201.80 ± 0.20-1.000.39b
Comorbidities
 Diabetes n (%)5 (25%)5 (18%)0.080.72a
 Hypertension n (%)6 (30%)8 (29%)0.010.91 a
 COPD n (%)1 (5.0%)2 (7.1%)0.040.76 a
 Coronary heart disease n (%)3 (15%)2 (7.4%)0.120.38 a
 Cancer n (%)1 (5.0%)5 (18%)0.190.19 a
Source of infection
 Hematogenous infection n (%)4 (20%)7 (25%)0.050.68 a
 Pulmonary infection n (%)17 (85%)20 (71%)0.150.27 a
 Intra-abdominal infection n (%)7 (35%)15 (53.6%)0.180.20 a
 Urinary tract infection n (%)1 (5.0%)3 (10.7%)0.100.48 a
 Infection in other sites n (%)7 (35%)5 (17.9%)0.190.17 a
Infectious pathogens
 Gram-positive bacteria n (%)12 (60%)15 (53.6%)0.060.65 a
 Gram-negative bacteria n (%)6 (30%)14 (50%)0.200.16 a
Disease Severity Scores
 APACHE II Score (points, IQR)22.00 (16.50,26.50)23.00 (18.00,26.50)-0.380.69 a
 SOFA score (points, ± s)12.10 ± 3.4012.10 ± 2.80-0.100.62 a

Comparison of baseline characteristics among patients with sepsis-induced cardiomyopathy, sepsis, and healthy controls.

SIC, Sepsis-induced cardiomyopathy group; SEPSIS, SEPSIS group; HC, Healthy control group; COPD, Chronic Obstructive Pulmonary Disease; APACHE II Score, Acute Physiology and Chronic Health Evaluation; SOFA Score, Sequential Organ Failure Assessment; a: Comparison between the sepsis-associated cardiomyopathy and sepsis groups; b: Comparison among the sepsis-associated cardiomyopathy, sepsis, and healthy control groups.

Analysis of immune function and inflammation-related indicators in both groups revealed that, except for the absolute total T lymphocyte count, which was significantly greater in the SIC group than in the SEPSIS group, there were no significant differences in the other indicators (all P>0.05). Comparisons of organ function-related indicators, such as creatinine and total bilirubin, between the two groups also revealed no statistically significant differences (all P>0.05, Table 2)

Table 2

VariableSEPSIS (n=20)SIC (n=28)F/χ2/Z/tP
Clinical Variables Related to Infection and Immunity
 White Blood Cells (109/L,IQR)10.20 (6.10,18.90)9.00 (6,13.6)-0.840.40
 Procalcitonin (ng/ml,IQR)9.00 (0.80,24.00)15.00 (3.20,37.00)0.780.43
 Interleukin-6 (pg/ml,IQR)130.50 (55.00,606.00)686.00 (51.90-2108.00)0.630.52
 Total T-lymphocyte percentage (%, ± s)61.80 ± 15.8070.00 ± 11.90-1.930.06
 Total B-lymphocyte percentage (%, IQR)18.10 (13.70,27.20)14.90 (11.70,18.70)-1.680.09
 NK cell percentage (%, IQR)7.70 (4.60,12.00)8.31 (4.70,10.60)0.001.00
 Absolute lymphocyte count (cells/μL, IQR)566.80 (326.70,807.81)804.90 (414.71,1051.80)1.320.18
 Absolute total T-lymphocyte count (cell/μL, IQR)192.70 (134.60,354.20)441.10 (156.70,691.20)-2.090.03*
 Absolute total B-cell count (cells/μL, IQR)104.40 (69.12,147.00)100.71 (45.21,160.10]-0.300.76
 Absolute NK-cell count (cells/μL, IQR)39.10 (19.70,88.80)54.40 (31.71,84.00)0.920.35
 Serum Immunoglobulin G (g/L, ± s)8.10 ± 2.518.50 ± 3.20-0.520.60
 Serum Immunoglobulin M (g/L, IQR)1.43 (0.91,2.42)1.40 (1.00,1.60)-0.420.66
 Serum Immunoglobulin A (g/L, IQR)1.10 (0.51,1.61)0.90 (0.60,1.20)-0.120.90
 Complement C3 (g/L, IQR)0.71 (0.61,0.95)0.60 (0.50,0.80)-1.440.14
 Complement C4 (g/L, IQR)0.20 (0.22,0.32)0.20 (0.10,0.30)-1.220.22
Clinical variables for organ function
 Creatinine (umol/L, ± s)95.91 ± 35.12120.80 ± 64.60-1.710.09
 Total Bilirubin (umol/L, IQR)23.63 (19.25,40.76)36.21 (19.42,55.92)0.880.37
 Lactate (mmol/L, IQR)1.80 (1.10,2.50)2.00 (1.41,2.80)0.940.34
 Hemoglobin (g/L, ± s)112.82 ± 31.61107.80 ± 18.600.630.53
 Platelets (109/L, ± s)117.41 ± 68.62120.70 ± 64.00-0.160.86

Comparison of laboratory data between patients with sepsis-induced cardiomyopathy and patients with sepsis.

SIC, Sepsis-induced cardiomyopathy group; SEPSIS, SEPSIS group; *P < 0.05.

In the SIC group, echocardiographic variables such as LVEF, absolute GLS values, E, and TAPSE were significantly lower than those in the SEPSIS group (all P<0.05). Compared with those in the SEPSIS group, the levels of myocardial damage markers such as cTnI were significantly greater (P<0.05), and the level of nt-proBNP in the SEPSIS group tended to increase, suggesting impaired cardiac function in SIC patients (Table 3).

Table 3

VariableSEPSIS (n=20)SIC (n=28)F/χ2/Z/tP
Echocardiographic Parameters 
GLS ( ± s)-16.68 ± 2.04-10.03 ± 2.630.68<0.001***
LVEF (%, ± s)67.77 ± 4.4952.60 ± 10.2710.23<0.001***
VTI (cm/s, ± s)19.24 ± 2.7713.28 ± 3.511.05<0.001***
SV (mL, ± s)59.548 ± 11.1246.5 0 ± 15.782.400.003**
E (m/s, ± s)0.86 ± 0.240.70 ± 0.220.240.026*
A (m/s, ± s)0.80 ± 0.240.66 ± 0.162.610.031*
E/A ( ± s)1.10 ± 0.371.00 ± 0.490.230.234
e’-S (cm/s, ± s)8.85 ± 3.297.00 ± 2.681.240.038*
e’-L (cm/s,  ± s)11.35 ± 3.119.32 ± 2.720.020.021*
E/e’ ( ± s)9.03 ± 3.048.98 ± 2.260.510.951
TVR (m/s,  ± s)2.44 ± 0.412.35 ± 0.400.020.459
TAPSE (mm,  ± s)22.72 ± 3.9714.72 ± 2.507.92<0.001***
IVC (%, ± s)0.19 ± 0.120.12 ± 0.120.180.053
LAVI (ml/m2, ± s)17.79 ± 5.9817.04 ± 4.780.690.634
S (cm/s, ± s)16.85 ± 2.7510.92 ± 3.440.98<0.001***
Myocardial injury markers
cTnI (ng/mL, IQR)0.00 (0-0.22)0.30 (0.00,0.50)1.990.04*
NT-proBNP (pg/mL, IQR)1640.00 (873.80,3870.90)4745.00 (1255.00,12750.00)1.700.08

Comparison of echocardiographic parameters and cardiac injury markers between patients with sepsis-induced cardiomyopathy and patients with sepsis.

SIC group, Sepsis-induced cardiomyopathy group; SEPSIS group, SEPSIS group; LVEF, Left ventricular ejection fraction; TVR, Tricuspid regurgitation velocity; LAVI, Left atrial volume index; e’-S, Interventricular septum e’ velocity; e’-L, Lateral wall e’ velocity; E, early diastolic mitral flow velocity; A, late diastolic mitral flow velocity; E/e’, E/e’ ratio; E/A, E/A ratio; TAPSE, tricuspid annular systolic excursion; S, systolic tricuspid annular velocity; IVC-CI, inferior vena cava collapse index; GLS, global longitudinal strain; VTI, time-integrated left ventricular blood flow velocity; cTnI, cardiac troponin I; nt-proBNP, N-terminal pro-B-type natriuretic peptide; *P < 0.05, **P < 0.01, ***P < 0.001.

Gut microbiome diversity in patients with sepsis-induced cardiomyopathy

The sequencing dilution curves for all the samples were flat (Supplementary Figure 1), indicating that the current sequencing depth was sufficient to cover the vast majority of the microbial species in the samples and that the data were suitable for subsequent analysis. We assessed the internal structure of the microbial communities and intergroup differences via both α diversity and β diversity.

In terms of α diversity, the Shannon and Simpson indices were significantly different among the three groups (P<0.01). The Shannon index in the SIC and SEPSIS groups was significantly lower than that in the HC group (P<0.05), the Simpson index in the SIC and SEPSIS groups was significantly greater than that in the HC group (P<0.05), and the ACE and Chao indices were significantly different between the SEPSIS and HC groups (P = 0.040, P = 0.03). (Figures 1A–C; Supplementary Figure 2; Supplementary Table 1).

Figure 1

In the β diversity analysis, principal coordinate analysis (PCoA) based on unweighted UniFrac distances revealed that the sample points of the SIC and SEPSIS groups were completely separated from those of the HC group. The sample points of the SIC and SEPSIS groups both overlapped and were partially separated, indicating that the gut microbiota of the SIC and SEPSIS groups differed significantly from that of the HC group overall (P<0.05), whereas the SIC and SEPSIS groups exhibited relatively smaller differences in their gut microbiota overall. PERMANOVA confirmed significant differences in the gut community structure among the three groups (R²=0.122, P = 0.001). Although the differences in the gut microbiota between the SIC and SEPSIS groups were not as pronounced as those between the SIC and HC groups were, significant differences still existed overall (P = 0.038, Figures 1D, E).

Changes in the gut microbiota structure in patients with S sepsis-induced cardiomyopathy

Compared with the SEPSIS group, the SIC group presented a downward trend in the relative abundances of Bacillota and Actinomycetota at the phylum level, while the relative abundance of Bacteroidota was significantly greater (0.34 vs. 0.18, P = 0.0198), and the Bacillota/Bacteroidota ratio was significantly lower (P<0.05). Compared with the HC group, the SIC group presented significant decreases in the relative abundances of Bacillota and Pseudomonadota (0.37 vs 0.58, P = 0.01; 0.21 vs 0.05, P = 0.023), whereas the relative abundance of Actinomycetota tended to increase, and the Bacillota/Bacteroidota ratio tended to decrease (Figures 2A, B).

Figure 2

At the genus level, the relative abundances of 71 bacterial genera differed significantly among the three groups (all P<0.05) (Supplementary Figure 3). We selected bacterial genera that presented significant differences and consistent trends between the SIC group and the other two groups for further analysis. Among these genera, the abundances of genera such as Clostridium innocuum group, Morganella, and Thomasclavelia were significantly elevated in the SIC group, whereas the relative abundances of Oscillospiraceae_UCG-003, Lachnospiraceae_UCG-004, Paraprevotella, and Eubacterium ruminantium group were significantly reduced in the SIC group (all P<0.05) (Figure 2C).

LEfSe analysis results indicated that Parabacteroides had the highest LDA value in the SIC group, making it the primary characteristic enriched genus in this group; therefore, this genus was selected for further study. Enterococcus was the primary enriched genus in the SEPSIS group, whereas Faecalibacterium and other genera were the primary enriched genera in the HC group (Figure 2D).

Changes in serum metabolites in patients with sepsis-induced cardiomyopathy

A total of 1,303 metabolites were detected in positive and negative ion scanning modes (Supplementary Table 2). Principal component analysis (PCA) based on the relative abundances of all the metabolites revealed that the metabolic distributions of the samples from the SIC and SEPSIS groups both partially overlapped and distinctly separated. (Figures 3A, B).

Figure 3

Using VIP≥1 and P<0.05 (t test) as screening criteria, a total of 49 differentially expressed metabolites were identified between the two groups. A volcano plot illustrates the differential expression of metabolites between groups, with 37 significantly upregulated and 12 significantly downregulated metabolites in the SIC group (Figure 3C; Supplementary Table 3). In SIC patients, the serum levels of differentially expressed metabolites such as hippuric acid, estrone glucuronide, trimethylamine N-oxide, 3-methoxy-4-hydroxyphenylglycol sulfate, 2,6-diaminopimelic acid, and N-acetyl-L-tyrosine were significantly elevated (all P<0.05), whereas the serum levels of differentially expressed metabolites such as 15-deoxy-Δ12,14-prostaglandin J2 were significantly decreased (P<0.05) (Figures 4A–G).

Figure 4

KEGG pathway enrichment analysis of these significantly differentially expressed metabolites revealed that they primarily participate in functional pathways, including phenylalanine metabolism, tyrosine metabolism, and arachidonic acid metabolism (Figure 4H).

Correlation analysis of the gut microbiota, serum metabolites, and clinical indicators in patients with sepsis-induced cardiomyopathy

Differential bacterial genera and metabolites in SIC patients were strongly and significantly correlated with echocardiographic parameters and myocardial injury biomarkers. (Supplementary Figure 4; Supplementary Tables 4–6). Among these genera, the 8 key gut bacterial genera significantly differentially expressed in the SIC group presented particularly close associations with key serum metabolites and clinical indicators, leading to the construction of a correlation network diagram. The results revealed significant correlations between the key differentially expressed bacterial genera and key serum metabolites identified in the SIC group and cardiac echocardiographic variables and myocardial injury markers. The relative abundance of the Parabacteroides genus was significantly positively correlated with estrone glucuronide and hippuric acid levels (r=0.290, P = 0.045 and r=0.369, P<0.01). Additionally, this genus showed a significant positive correlation with GLS (r=0.358, P = 0.012) and a significant negative correlation with LVEF (r=−0.420, P<0.01); the relative abundance of the Oscillospiraceae_UCG-003 genus was significantly negatively correlated with N-acetyl-L-tyrosine (r=−0.295, P = 0.038) and was also significantly negatively correlated with GLS and cTnI (r=−0.287, P = 0.047 and r=−0.298, P = 0.039) and was significantly positively correlated with TAPSE (r=0.307, P = 0.033); the relative abundance of the Lachnospiraceae_UCG-004 genus was significantly negatively correlated with hippuric acid and 2,6-diaminopimelic acid (r=−0.352, P = 0.013 and r=−0.333, P = 0.020), was negatively correlated with GLS (r=−0.304, P = 0.035), and was significantly positively correlated with TAPSE (r=0.309, P = 0.032) (Figure 5). To integrate the holistic relationships across multiple omics, we performed CCA. The results revealed that samples from the SEPSIS and SIC groups were clearly divided into two clusters in the ordination space. The metabolite 15-Deoxy-.Δ.12,14-prostaglandin J2 and the clinical indicator NT-proBNP, indicated by the longest arrow lengths, were the primary environmental factors driving the differentiation of the microbial community structure. Genus groups such as Morganella, Parabacteroides, Clostridium innocuum group, and Thomasclavelia were positively correlated with the distribution in the SIC group, whereas Eubacterium ruminantium group, Oscillospiraceae_UCG-003, Lachnospiraceae_UCG-004, and Paraprevotella were negatively correlated. Furthermore, there was a broad correlation between key bacterial genera, serum differentially abundant metabolites, and clinical indicators related to SIC diagnosis (Figure 6).

Figure 5

Figure 6

Discussion

Our data agree with earlier findings on dysbiosis in SIC, which is consistent with the findings of previous limited studies (; ) reporting dysbiosis in SIC patients, thereby providing preliminary support for our results. However, these studies focused primarily on comparisons between SIC patients and non-SIC patients and did not elucidate the characteristics of microbial evolution during the progression from a healthy state to sepsis and subsequently to SIC. Our study analyzed and compared the gut microbiota of SIC patients with that of non-SIC sepsis patients and healthy individuals, revealing that SIC patients exhibit unique gut microbiota characteristics. This is manifested as follows: at the phylum level, during the progression from a healthy state to sepsis and subsequently to SIC, the relative abundance of Bacillota gradually decreased; compared with that in patients in the SEPSIS group, the Bacillota/Bacteroidota ratio was significantly lower in SIC patients; at the genus level, compared with those in patients in the SEPSIS and HC groups, the relative abundances of opportunistic pathogens such as Parabacteroides, Morganella, Clostridium innocuum group, Thomasclavelia and Paraprevotella were significantly greater, whereas the relative abundances of probiotic genera such as Oscillospiraceae_UCG-003 and Lachnospiraceae_UCG-004 were significantly lower.

The Bacillota/Bacteroidota ratio is often regarded as a marker of gut microbiota dysbiosis (), suggesting that gut microbiota dysbiosis progresses as the condition evolves from a healthy state to sepsis and ultimately to SIC. The Parabacteroides genus plays a crucial role in maintaining intestinal homeostasis and regulating immunity and metabolism; however, under specific conditions, it can also transform into an opportunistic pathogen (; ). Studies have shown that Parabacteroides distasonis is positively correlated with cholesterol, low-density lipoprotein, and cardiac creatine kinase isoenzymes (), and its relative abundance is elevated in patients with hypertension (), suggesting that it may be a potential pathogen in cardiovascular disease. As opportunistic pathogens, Clostridium innocuum group is associated with a high 30-day mortality rate in infected patients and is frequently accompanied by sepsis, shock, and infective endocarditis (; ). Its pathogenicity may depend on lipid rafts, which may activate inflammatory pathways—including NOD2/NF-κB—thereby exacerbating inflammatory responses and causing damage to distant organs (). Among patients with bloodstream infections caused by the genus Morganella, 12.6% developed myocardial infarction, and 25.3% developed congestive heart failure; these cardiovascular complications significantly increased patient mortality (). The high mortality rate may be related to the unique antibiotic resistance mechanisms and virulence factors of bacteria (). The genera Lachnospiraceae_UCG-004 and Oscillospiraceae_UCG-003 are considered gut probiotics, and many of their members are short-chain fatty acid (SCFA) producers (; ). This study revealed that the abundance of these genera was significantly reduced in SIC patients and that Oscillospiraceae_UCG-003 and Lachnospiraceae_UCG-004 may play a role in reducing SCFA production. No short-chain fatty acids were detected in the serum metabolites in this study, which may be related to the reduction in a large number of beneficial bacterial genera and impaired short-chain fatty acid production. SCFAs are important alternative energy substrates for the myocardium; sepsis leads to massive depletion of gut probiotics, resulting in SCFA deficiency, which can exacerbate myocardial energy metabolism disorders and induce cardiac dysfunction (; ). SCFAs also act as important immunomodulators. Through G protein-coupled receptors such as GPR41, GPR43, and GPR109A, as well as epigenetic regulation via the inhibition of histone deacetylases, they promote the development and functional activation of regulatory T cells and suppress the secretion of proinflammatory cytokines (), suggesting that they may exert a protective effect on myocardial function in sepsis via similar pathways.

This study revealed that echocardiographic parameters were significantly more abnormal in patients with SIC than in patients with sepsis, which is consistent with the characteristics of this disease. Our previous research identified the GLS, E-wave, and TAPSE as echocardiographic predictors of SIC, in addition to LVEF (). Tricuspid annular systolic displacement (TAPSE) is a reliable indicator for assessing right ventricular systolic function, and a TAPSE value <16 mm is the standard for diagnosing right ventricular dysfunction (). Global longitudinal strain (GLS) is a specific indicator for assessing left ventricular systolic function; it is normally negative, and a larger negative value indicates better left ventricular systolic function (). This study revealed that in SIC patients, the relative abundances of the intestinal genera Lachnospiraceae_UCG-004, Oscillospiraceae_UCG-003 and Paraprevotella were negatively correlated with GLS and positively correlated with TAPSE, indicating that a significant reduction in the levels of these genera is associated with impaired left and right ventricular systolic function. The relative abundances of genera such as Parabacteroides and Morganella were positively correlated with GLS, whereas those of genera such as Thomasclavelia and Morganella were negatively correlated with TAPSE, suggesting that significantly elevated levels of these genera are associated with reduced left and right ventricular systolic function. cTnI and nt-proBNP are the most widely used biomarkers for the clinical assessment of SIC; they play complementary and critical roles in identifying myocardial injury, evaluating cardiac function, and predicting prognosis (; ). Oscillospiraceae_UCG-003 was significantly negatively correlated with cTnI, suggesting that a significant reduction in this probiotic genus may affect myocardial function. These findings indicate a broad correlation between gut microbiota dysbiosis and myocardial injury in SIC patients.

SIC is a metabolic disease (), and its gut microbiota dysbiosis may affect cardiac function through serum metabolites. In the context of sepsis, damage to the intestinal barrier and dysbiosis trigger drastic changes in circulating metabolites, leading to transient, severe fluctuations in metabolite concentrations. This exposes the myocardium to acute inflammation and energy stress, inducing the reversible cardiac dysfunction characteristic of SIC (). Microbial metabolites exert a bidirectional effect on cardiomyocytes: on the one hand, short-chain fatty acids and certain indole derivatives have anti-inflammatory and endothelial-protective effects; on the other hand, compounds such as trimethylamine oxide disrupt myocardial energy metabolism (). Metabolomic studies have confirmed that sepsis is accompanied by widespread metabolic disturbances, which can modulate organ dysfunction (). Consistent with these advances in metabolomics research, the nontargeted serum metabolomics analysis in this study further revealed characteristic metabolic disturbances in patients with SIC. Among these, changes in several key differentially expressed metabolites may have significant clinical implications. Hippuric acid is a cometabolite produced by the host and gut microbiota. The results of the present study revealed that the abundance of this genus increased with increasing relative abundances of the genera Clostridium innocuum group and Parabacteroides. Additionally, hippuric acid was significantly positively correlated with GLS, an indicator of left ventricular systolic function, suggesting an adverse effect on cardiac contractility. In a sepsis model induced by Escherichia coli infection, HA enhances M1-like polarization of macrophages and amplifies MyD88-dependent Toll-like receptor ligand-mediated inflammatory responses; moreover, elevated HA levels in vivo are significantly associated with mortality in sepsis patients (). Studies have shown that in patients with chronic kidney disease (CKD), high HA levels are associated with atherosclerosis and left ventricular hypertrophy (); the mechanism may involve the induction of ROS production by mitochondria, thereby exacerbating the atherosclerotic process (). Gut bacteria produce β-glucuronidase (GUS), which converts estrone glucuronide into active estrone, one of the active forms of estrogen, which is believed to have a protective effect on the cardiovascular system (). This study revealed that as the relative abundance of the Parabacteroides genus increased, the serum levels of estrone glucuronide increased, and this increase was significantly positively correlated with the cardiac injury marker NT-proBNP, suggesting that myocardial damage may occur with increasing levels of this metabolite. The mechanism may be related to the decrease in free estrogen levels caused by elevated estrone glucuronide levels. Trimethylamine (TMA) is a purely microbial metabolite whose production depends directly on the metabolism of the gut microbiota; TMA is converted in the liver to trimethylamine N-oxide (TMAO) (), and its serum levels increase with increasing abundance of Thomasclavelia. Trimethylamine N-oxide is a key factor that promotes the onset and progression of cardiovascular disease and is also a risk factor for cardiovascular disease (). Studies have shown that TMAO can activate proinflammatory signaling pathways such as the NF-κB and MAPK pathways in vascular endothelial cells and smooth muscle cells, impairing vascular endothelial function, inducing vascular inflammation, and exacerbating cardiac dysfunction (). 15-Deoxy-Δ12,14-prostaglandin J2 is a metabolite of the arachidonic acid epoxidase pathway (). Its serum levels are significantly positively correlated with LVEF and TAPSE and negatively correlated with GLS, suggesting that reduced serum levels may be associated with impaired left and right ventricular systolic function. Studies have confirmed that a decrease in the plasma level of postpercutaneous coronary intervention (PCI) is positively correlated with the extent of myocardial injury and that it reduces myocardial cell damage by inhibiting the excessive activation of Ca2+/calmodulin-dependent kinase II (CaMKII) ().

To explore the relationships among SIC-associated gut microbiota dysbiosis, serum metabolic changes, and their associations with disease, we performed a Spearman correlation analysis, which revealed a statistically significant correlation among the three factors. To further integrate the evolution of the gut microbiota structure and holistic interrelationships across multiple omics, we performed a CCA analysis (). The results revealed that the gut bacterial genera in the SEPSIS and SIC groups were clearly divided into two distinct clusters in the ordination space, suggesting significant differences in the gut microbial community structure between the two groups. Furthermore, gut microbiota dysbiosis, serum metabolic alterations, and clinical indicators in SIC patients exhibit multiple interrelated patterns. These findings, together with the results of the Spearman correlation analysis, confirm from different perspectives that SIC patients exhibit a characteristic gut microbiota–metabolism–cardiac function association pattern distinct from that of sepsis patients, suggesting an interaction between gut microbiota dysbiosis and altered serum metabolism with SIC. The above results further indicate that there is a broad association among characteristic gut dysbiosis, differential serum metabolites, and cardiac dysfunction in SIC patients.

This study has certain limitations: the study sample size was relatively small, and it was a single-center study. The metabolites identified were not subjected to FDR multiple correction. The interactions and associations between the differentially expressed gut bacterial genera and metabolites in SIC and the clinical indicators related to the diagnosis of SIC reflect only statistical correlations. Larger-scale, multicenter studies and animal experiments may be needed in the future to further validate the biological significance of the “gut microbiota–serum metabolites–clinical indicators” association in SIC.

Conclusion

This study describes the evolution of gut microbiota phylum/genus-level characteristics during the progression from healthy individuals to sepsis and subsequently to SIC, as well as changes in the serum metabolome of SIC patients. Key bacterial genera, differentially expressed metabolites, and clinical indicators associated with SIC diagnosis were all found to be correlated, suggesting that the gut microbiota may interact with SIC either directly or by regulating serum metabolites. These findings provide a theoretical basis for targeting improvements in the gut microbiome to enhance cardiac function and prognosis in SIC patients.

Statements

Data availability statement

The data presented in this study have been deposited in the NCBI Data Repository under accession numbers PRJNA1477312 (SIC group and SEPSIS group) and PRJNA691455 (healthy control group).

Ethics statement

This study was conducted in accordance with the Declaration of Helsinki (revised 2024) and was approved by the Institutional Review Board of the Ethics Committee of the General Hospital of Ningxia Medical University (Approval No: KYLL-2025-0038). Informed consent was obtained from all participants. 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

KC: Formal analysis, Writing – original draft, Data curation. JL: Writing – original draft, Visualization, Software, Validation. YZ: Writing – original draft, Visualization, Software. YL: Writing – original draft, Validation, Visualization. XY: Funding acquisition, Project administration, Writing – review & editing, Resources. XY: Conceptualization, Writing – review & editing, Funding acquisition, Formal analysis.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the National Key Research and Development Program of China (Grant No. 2023ZD0506503) and the Ningxia Key Research and Development Program (Grant No. 2023BEG03040). The National Key Research and Development Program (2023ZD0506503) was listed as the first source of funding.

Acknowledgments

The authors appreciate the contributions of all the ICU teams and participants involved in the study.

Conflict of interest

The author(s) declared that this work 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) declared that generative AI was not 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/fcimb.2026.1811654/full#supplementary-material

References

Summary

Keywords

gut microbiota, intensive care unit, metabolomics, sepsis, sepsis-induced cardiomyopathy

Citation

Chen K, Lin J, Zhang Y, Liu Y, Yang X and Yang X (2026) Gut microbiota and serum metabolic profiles in patients with sepsis-induced cardiomyopathy and their association with the disease. Front. Cell. Infect. Microbiol. 16:1811654. doi: 10.3389/fcimb.2026.1811654

Received

15 February 2026

Revised

01 June 2026

Accepted

02 June 2026

Published

10 July 2026

Volume

16 - 2026

Edited by

Valeriy Poroyko, Laboratory Corporation of America Holdings (LabCorp), United States

Reviewed by

Johid Malik, University of Nebraska Medical Center, United States

Swarnima Pandey, University of Maryland, United States

Updates

Copyright

*Correspondence: Xiaojuan Yang, ; Xiaojun Yang,

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