The study enrolled 62 patients, of whom 27 died (44%) and 35 survived (56%) after 28 days of mechanical ventilation. Each patient provided two sets of BALF and stool samples, one within 48 hours (Group 01) and the other at one week after mechanical ventilation (Group 02). Table 1 displays the clinical characteristics of the patients. Importantly, there were no significant differences in clinical features, culture results, or antibiotic use between the deceased and survivor groups.
Table 1
Baseline characteristics and outcome.
| Deceased(27) | Survival(35) | P |
---|
Age | 63(53, 69) | 58(45, 75) | 0.121 |
length of ICU stay(d) | 17(11, 23) | 13(8, 23) | 0.286 |
Enteral nutrition time(d) | 2(1, 4) | 2(1, 4) | 0.896 |
SOFA | 14(12, 16) | 13(8, 14) | 0.098 |
Antifungal | 11 | 11 | 0.447 |
Sex(Male) | 15 | 17 | 0.585 |
Emergency Surgery | 11 | 12 | 0.602 |
Smoking history | 7 | 6 | 0.400 |
Alcohol history | 7 | 7 | 0.580 |
VAP | 17 | 19 | 0.492 |
ARDS | 21 | 22 | 0.206 |
Early enteral nutrition (72h) | 20 | 26 | 0.985 |
Diagnosis-abdominal surgery | 3 | 8 | 0.321 |
Diagnosis-hematopathy | 3 | 4 | 1.000 |
Diagnosis-Rheumatic immune disease | 2 | 3 | 1.000 |
Diagnosis-respiratory disease | 9 | 5 | 0.124 |
Diagnosis-neurosurgery | 10 | 10 | 0.480 |
Diagnosis-respiratory disease | 0 | 5 | 0.063 |
complication-hypertension | 4 | 5 | 1.000 |
complication-diabetes | 4 | 3 | 0.698 |
complication-cardiovascular disease | 5 | 4 | 0.485 |
complication-tumour | 4 | 5 | 1.000 |
complication-Neuropathy | 4 | 6 | 1.000 |
Antibiotics | | | |
Enterobacter aerogenes | 6 | 6 | 0.616 |
Carbapenems | 16 | 16 | 0.290 |
Cephalosporins | 9 | 13 | 0.756 |
β- Lactam antibiotics | 12 | 14 | 0.725 |
Tegacyclin | 9 | 6 | 0.140 |
Aminoglycoside | 0 | 2 | 0.500 |
Anti positive bacteria | 16 | 16 | 0.290 |
Quinolone | 3 | 3 | 1.000 |
Polymyxin | 2 | 3 | 1.000 |
Anti anaerobic bacteria | 1 | 0 | 0.436 |
BALF traditional culture | | | |
CRAB | 3 | 7 | 0.491 |
CRKP | 14 | 18 | 0.974 |
CRPA | 6 | 7 | 0.831 |
MRSA | 4 | 10 | 0.199 |
Neisseria | 5 | 5 | 0.735 |
Stenotrophomonas maltophilia | 4 | 4 | 0.719 |
Burkholderia cepacia | 2 | 3 | 1.000 |
Chryseobacterium indologenes | 3 | 3 | 1.000 |
Table 1: Data presented as median with [IQR].SOFA:Sequential Organ Failure Assessment. Definition of abbreviations: VAP:Ventilator associated pneumonia; ARDS:acute respiratory distress syndrome; CRAB:Cabapemne Resistant Acinetobacter Baumannii; CRKP:carbapenem-resistant Klebsiella pneumoniae; CRPA:Carbapenem-resistant Pseudomonas aeruginosa; MASA:methicillin-susceptible Staphylococcus aureus.
Comparison of microbial Alpha diversity in survivors and deceased
We performed 16S rRNA gene sequencing on all samples, with BALF and fecal matter diluted to saturation levels (Fig. 2A,B), indicating sufficient sample coverage and distribution. The α-diversity, estimated by the Shannon index, reflects the diversity and abundance of microbiota. Our results demonstrate a decreasing trend in α-diversity over time in both the BALF and fecal microbiome, with the deceased group exhibiting significantly lower microbial diversity in both the lung (Fig. 2C) and gut (Fig. 2D) compared to the survival group (p < 0.05). To further investigate the association between alpha diversity and clinical outcomes, we analyzed Kaplan-Meier curves for Group 01 patients (Fig. 2E,F). The findings revealed that individuals with low α diversity of lung microbes had a higher likelihood of mortality than those with high α-diversity (p < 0.05).
Comparison of the relative abundance of BALF and fecal microbiota in survivors and deceased
The microbial composition of intestinal and respiratory samples was investigated to identify the prevailing microbial taxa at each site. In respiratory samples, Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria emerged as the dominant phylum, collectively representing almost 70% of the total sequencing reads. Likewise, in the intestinal samples, Firmicutes, Proteobacteria, and Bacteroidetes predominated, accounting for nearly 90% of the total sequencing reads. Prolonged mechanical ventilation was found to accentuate the relative abundance of these phylum. In BALF samples, the abundance of Proteobacteria was higher in the deceased (Group 01: 57.5656% vs. 23.1363%; Group 02: 75.8837% vs. 21.6580%) compared to survivors, whereas the reverse pattern was observed for Firmicutes (Group 01: 18.819% vs. 54.2770%; Group 02: 19.0444% vs. 40.8464%), and this trend intensified over time (Fig. 3A). In feces, the abundance of Bacteroidetes in deceased was lower than that in survivors (Group 01: 19.0477% vs. 17.5287%; Group 02: 21.9145% vs. 21.9145%). While no conspicuous difference was found, Proteobacteria exhibited higher levels in the deceased group (group 01: 18.2066% vs. 32.2428%; group 02: 32.2428% vs. 21.6580%), whereas Firmicutes were more abundant in the survivor group (group 01: 54.2770% vs. 53.1729%; group 02: 40.8464% vs. 64.059%)(Fig. 3B). The Venn diagram reveals a substantial reduction in specific bacterial genera in the deceased group compared to the surviving group, with the exception of shared microbes (Fig. 3C,D). In the respiratory tract of the survival group, Prevotellaceae, Moraxellaceae, and Staphylococcaceae were the predominant families. Interestingly, the respiratory tract of the deceased group displayed distinct characteristics, including the presence of Enterobacteriaceae (p < 0.05), enrichment of Enterococcaceae, and an increase in opportunistic pathogens Xanthomonadaceae (p < 0.05) (Fig. 3E, G). Conversely, the intestines of the survival group exhibited a dominance of symbiotic bacteria, such as Veillonellaceae (p<0.05), Peptostreptococcaceae (p<0.05), Bacteroidaceae, Lachnospiraceae, Lactobacillaceae, among others. However, the intestines of the deceased group were primarily characterized by Enterococcaceae (p < 0.05)(Fig. 3F, H).
Comparison of microbial Beta diversity in survivors and deceased
Beta diversity is a measure of species diversity among ecosystems, according to the absence or presence of species and their abundance. To further clarify the association between respiratory and intestinal microorganisms, we used NMDS and PLS-DA to construct a model for differential analysis of specimen clustering results.It is known that biogeographic barriers exist between respiratory tract and intestine microbial community within individuals. Our results indicate that beta diversity is significantly different between lung and gut in survivors, but not in deceased(p < 0.05)(Fig. 4A,B). There were significant differences between microorganisms in these two regions. However, in the survivor group, Akkermansia, Streptococcus, Prevotella, and Lactobacillus were enriched in the intestinal and respiratory tracts at the beginning of mechanical ventilation, and the commensal flora of both intermingled with each other (Fig. 4C). The respiratory tract was enriched with intestinal beneficial bacteria. However, in the deceased group, the respiratory and intestinal tracts were simultaneously dominated by an increase in pathogenic bacteria and a decrease in commensal bacteria. The aggregation of highly pathogenic Enterobacteriaceae and Enterococcus in the intestine seems to explain the large enrichment of intestinal microbiota in the respiratory tract among the deceased (p < 0.05)(Fig. 4D).
Comparing the fecal SCFA concentrations in survivors and deceased
SCFAs are metabolic by-products of symbiotic intestinal bacteria fermenting dietary fiber, known for their anti-inflammatory, immune-modulatory effects, and distal lung protection. We quantified seven SCFAs, including Caproic acid, Pentanoic acid, Butyric acid, Isobutyric acid, Isovaleric acid, Propionic acid and Acetic acid, in the fecal microbiota of Group 01. Our findings reveal significantly lower concentrations of Pentanoic acid, Butyric acid, Isobutyric acid, and Isovaleric acid in the deceased group compared to the surviving group(P < 0.05)(Fig. 5A).
Associations between lung and gut microbiota and 28-day mortality
We employed Cox proportional hazard regression models to investigate the potential links between changes in intestinal and pulmonary microbiota and 28-day mortality after the initiation of mechanical ventilation. During the follow-up period, a total of 27 patients (44%) succumbed within 28 days. After conducting a preliminary univariate analysis, we specifically focused on lung and gut microbial diversity, BALF Enterobacteriaceae, and fecal Enterococcus, undertaking a multifactor Cox regression analysis. Our results revealed significant associations between lung microbial diversity and fecal Enterococcaceae, both of which independently correlated with 28-day mortality. The results showed that pulmonary microbial diversity and fecal Enterococcaceae were correlated and independently associated with 28-day mortality (adjusted hazard ratio (aHR) = 0.773 ; 95% confidence interval (CI) 0.652 ~ 0.916, p = 0.003; and aHR = 1.022; 95% CI, 1.008–1.037, p = 0.002)(Table 2). These findings underscore the potential prognostic value of analyzing intestinal and pulmonary microbiota in predicting 28-day mortality among mechanically ventilated patients.
Table 2 Definition of abbreviations: HR hazard ratio, aHR adjusted hazard ratio, 95%CI 95% confidence interval.