Patient characteristics
The study included a total of 95 PDS patients, with 60 patients in the mNGS group, and 35 patients in the non-mNGS group. Table 1 documented the patient characteristics and baselines information. Statistical analyses were conducted to assess factors like age, sex, SOFA score, comorbidities, and patient outcomes. All factors, except for the SOFA score, showed no significant differences, suggesting a relatively similar baseline between the two groups. Despite the longer ICU length of stay in the mNGS group compared to the non-mNGS group (20 vs 8, P = 0.115), the mNGS group exhibited a significantly lower mortality rate (35.0% vs 57.1%, P = 0.034) (Table 1).
Table 1
Clinical characteristics of participants.
Characteristic | mNGS (n = 60) | Non-mNGS (n = 35) | P-value |
Age, median (IQR) | 68 (60–77) | 73 (62–85) | 0.177 |
Sex, n (%) | | | |
Male | 47 (78.3%) | 26 (74.3%) | 0.652 |
Female | 13 (21.7%) | 9 (25.7%) | |
SOFA score, median (IQR) | 11 (7–13) | 8 (5–11) | 0.039* |
Comorbidities, n (%) | | | |
Respiratory failure | 42 (70.0%) | 27 (77.1%) | 0.451 |
Hypertension | 39 (65.0%) | 19 (54.3%) | 0.302 |
Diabetes | 23 (38.3%) | 8 (22.9%) | 0.121 |
Hypoproteinemia | 14 (23.3%) | 13 (37.1%) | 0.150 |
CAD | 6 (10.0%) | 8 (22.9%) | 0.088 |
MODS | 30 (50.0%) | 15 (42.9%) | 0.501 |
ARDS | 20 (33.3%) | 17 (48.6%) | 0.142 |
CKD | 14 (23.3%) | 7 (20%) | 0.706 |
Outcomes, median (IQR) | | | |
Hospital length of stay(days) | 31 (22–53) | 22 (9–46) | 0.111 |
ICU length of stay(days) | 20 (11–28) | 8 (5–14) | 0.115 |
Mortality(%) | 21 (35.0%) | 20 (57.1%) | 0.034* |
IQR, interquartile range; SOFA, sepsis-related organ failure assessment; CAD, Coronary atherosclerotic heart disease; MODS, Multiple organ dysfunction syndrome; ARDS, Acute respiratory distress syndrome; CKD, Chronic kidney disease.*,P-values ≦ 0.05; **,P-values ≦ 0.01; ***,P-values ≦ 0.001 |
Detection performance of plasma and BALF mNGS
The performance of mNGS was compared to CMT in all patients (Fig. 2A). The pathogen positivity rate was higher in plasma and BALF mNGS compared to blood and BALF culture, respectively (61.7% vs 6.7%, P < 0.001; 100% vs 40%, P < 0.001). Specifically, BALF mNGS exhibited a 100% positivity rate. No significant difference was observed in the positivity rates of different CMT methods between the mNGS and non-mNGS groups.
Subsequently, we analyzed the quantity of microorganisms identified by various methods (Fig. 2B). Within the mNGS group, 8.3% of the patients exhibited a singular microorganism through plasma mNGS, whereas 55% identified more than 4 microorganisms through BALF mNGS. The CMT method identified 1 or 2 microorganisms in 78.3% of patients in the mNGS group and 65.7% of patients in the non-mNGS group (Additional file 1: Figure S1).
Then we assessed the consistency between mNGS and CMT methods. Among the mNGS group, a total of 33 patients were concurrently detected with microorganisms in plasma mNGS, BALF mNGS, and CMT (Fig. 2C). Out of the 33 patients who tested positively for both plasma mNGS and CMT, only 33% had consistent microbial detection (Fig. 2D). Of the 52 patients who tested positive for both BALF mNGS and CMT, 75% displayed at least one consistent microbial detection (Fig. 2E). When combing plasma and BALF sample mNGS results, out of the 52 patients who tested positively for both mNGS and CMT, 79% had consistent microbial detection (Fig. 2F).
Clinical adjudication of plasma and BALF mNGS results
In the mNGS group, 37 patients had both plasma and BALF mNGS test results positive, and among them, 28 patients were identified to have at least one consistent microorganism (Fig. 3A). A total of 284 microorganisms were identified through mNGS. Out of these, 60 microorganisms were identified in the plasma samples, while 260 microorganisms were identified in BALF samples, and 36 of them were consistently identified in both sample types (Fig. 3B). Out of the 36 detections, 55.6% (20/36) was classified as definite causes of the sepsis alert (Fig. 3C). Among the microorganisms identified exclusively in either plasma or BALF mNGS, only 20.8% (5/24 in plasma) and 18.8% (42/224 in BALF) were classified as definite causes of the sepsis alert (Fig. 3C). The proportion of definite causes increased by about 35% when microorganisms were identified in both samples compared to those identified in only one sample, the proportion of microorganisms classified as "Definite + Probable" in cases where microorganisms were simultaneously detected in both samples (88.9%) increased by 9 to 26% compared to microorganisms detected in a single sample alone (62.5% in plasma and 79.9% in BALF). This means that detecting consistent microorganisms can better clarify the clinical significance of pathogens.
Microorganisms identified in patients with PDS
A total of 26 kinds of pathogens detected by plasma and BALF mNGS were classified as definite causes of the sepsis alert, including 15 bacteria, 5 fungi, 4 viruses, and 2 special pathogens (Fig. 4). Acinetobacter baumannii, Stenotrophomonas maltophilia, Candida albicans, and Human mastadenovirus B were the most common pathogens in PDS patients. Noting the particular sensitivity of plasma in identifying special pathogens such as Orientia tsutsugamushi and Mycobacterium tuberculosis complex. A total of 43 pathogens were classified as probable causes of sepsis alert. The most frequently detected bacteria were Corynebacterium striatum, Staphylococcus haemolyticus, Enterococcus faecium, and Burkholderia cepacia, as well as Epstein-Barr virus and cytomegalovirus. Similar to the “definite” pathogens, viruses were frequently detected in both plasma and BALF, whereas fungi were exclusively detected in BALF (Additional file 2:Figure S2). The CMT method identified a total of 19 pathogens in both the mNGS and non-mNGS groups. The most common pathogens in both the mNGS groups and non-mNGS groups were Candida app (Additional file 3:Figure S3).
Then we analyzed the pathogen burden based on RPTM among microorganisms identified by plasma and BALF mNGS. There were no significant differences in pathogen burden among microorganisms classified as definite, probable, possible, and unlikely causes of the sepsis alert in plasma (Fig. 5A). The pathogen burden of microorganisms classified as definite causes of the sepsis alert was significantly higher than that of microorganisms classified as probable, possible, and unlikely causes of the sepsis alert in BALF (Fig. 5B). Moreover, the pathogen burden of microorganisms identified in both plasma and BALF samples was significantly higher than that of microorganisms identified only in plasma or BALF samples (Fig. 5C and D).
The associations between inflammatory markers and mNGS results
To investigate potential factors associated with the matched results between plasma and BALF mNGS findings, we compared differences in inflammatory markers, including WBC, Hb, PLT, NE%, LY%, NE, HCT%, CRP, and PCT. The inflammatory markers between the mNGS and non-mNGS groups were generally similar, except for the Hb index, which was lower in the mNGS group (Additional file 4: Table S1). Subsequently, we compared the differences between the matched or partially matched cases and unmatched cases. There were no significant differences in Hb, PLT, NE%, LY%, HCT%, CRP, and PCT between these two groups. However, the unmatched cases exhibited significantly higher WBC indices (11.03 ± 5.42 versus 14.30 ± 6.50 10^9/L, P = 0.041) and NE indices (8.75 ± 5.27 vs 12.17 ± 5.78 10^9/L, P = 0.021) compared to the matched or partially matched cases (Table 2).
Table 2
Comparison of inflammatory markers between matched or partially matched cases and unmatched cases of plasma and BALF mNGS detection.
Laboratory tests | Plasma and BALF mNGS matched or partially matched cases (n = 28) | Plasma and BALF mNGS unmatched cases (n = 32) | P-value |
WBC (10^9/L) | 11.03 ± 5.42 | 14.30 ± 6.50 | 0.041* |
Hb (g/L) | 90.64 ± 22.60 | 88.42 ± 22.95 | 0.707 |
PLT (10^9/L) | 187.54 ± 120.41 | 242.75 ± 138.20 | 0.107 |
NE% | 82.14 ± 15.64 | 83.44 ± 8.87 | 0.689 |
LY% | 11.22 ± 12.10 | 11.89 ± 17.83 | 0.866 |
NE (10^9/L) | 8.75 ± 5.27 | 12.17 ± 5.78 | 0.021* |
HCT% | 30.22 ± 22.46 | 27.41 ± 7.60 | 0.507 |
CRP (mg/L) | 108.03 ± 88.30 | 116.61 ± 94.32 | 0.736 |
PCT (ng/ml) | 1.32 ± 2.32 | 8.84 ± 22.37 | 0.153 |
WBC, white blood cell; Hb, hemoglobin; PLT, blood platelet; NE%, Percentage of neutrophils in blood; LY%, Percentage of lymphocytes and total white blood cell; NE, neutrophil; HCT%, hematocrit; CRP, C-reactive protein; PCT, procalcitonin. |
Application of mNGS for antibiotic adjustment
All hospitalized patients received antibiotic treatment. Among the patients in the mNGS group, 56 individuals received antibiotic treatment prior to undergoing mNGS testing (Table 3). Amon the 56 patients, 85.71% (48/56) had adjustments made to their antibiotic treatment, a higher percentage compared to the non-mNGS group (33.33%, 11/33) (Table 3). Based on the mNGS results, 19.64% (11/56) of the patients had a reduction in the number or spectrum of initial empiric antimicrobial drugs, while no patients in the non-mNGS group experienced de-escalation. Antibiotic escalation was observed in 42.86% (24/56) of patients in the mNGS group, primarily due to the addition of targeted antimicrobial drugs. Additionally, one patient in the mNGS group showed symptom improvement and was discharged from the hospital before the mNGS results were obtained (Table 3).
Table 3
The adjustment of antibiotics in the mNGS group and the non-mNGS group.
Antibiotic adjustments | mNGS (n = 56) | Non-mNGS (n = 33) | P value |
Adjustments | 48 (85.71%) | 11 (33.33%) | < 0.001*** |
De-escalation | 11 (19.64%) | 0 (0.00%) | 0.007** |
Reduce spectrum of agents | 3 (5.36%) | 0 (0.00%) | 0.176 |
Reduce number of agents | 8 (14.29%) | 0 (0.00%) | 0.023* |
Escalation | 24 (42.86%) | 4 (12.12%) | 0.003** |
Increase number of agents | 20 (35.71%) | 4 (12.12%) | 0.015* |
Increase spectrum of agents | 4 (7.14%) | 0 (0.00%) | 0.116 |
Same level replacement | 13 (23.21%) | 7 (21.21%) | 0.827 |
Unadjustments | 7 (12.5%) | 22 (66.67%) | < 0.001*** |
No change | 7 (12.5%) | 22 (66.67%) | < 0.001*** |
Discharged before mNGS | 1(1.79%) | | |
Next, we compared the antibiotic adjustment between the matched or partially matched cases and unmatched cases from dual sample mNGS tests. Among the 27 patients with matched or partially matched results, 81.48% (22/27) had adjustments in antibiotic treatment based on the mNGS results, while 89.66% (26/29) had adjustments in antibiotic treatment in the 29 patients with unmatched results (Additional file 5: Table S2). Compared to patients with unmatched results, more patients in the group with matched or partially matched results increased spectrum of agents (14.81% vs 0.00%, P = 0.032) (Table S2).
Factors influencing prognosis in patients with PDS
Compared to the non-mNGS group, the mNGS group exhibited a prolonged length of stay in the ICU and a reduced mortality rate (Table 1). Further analysis was conducted to compare the clinical outcomes among cases with matched or partially matched plasma and BALF mNGS results and those with unmatched results. The mortality rate was lower in the plasma and BALF consistent mNGS group compared to the inconsistent group, but the difference was not statistically significant (35.71% vs 40.61%, P = 0.696) (Table 4).
Table 4
Clinical outcomes between matched or partially matched cases and unmatched cases of plasma and BALF mNGS detection.
Outcomes | Blood and BALF mNGS matched or partially matched cases (n = 28) | Blood and BALF mNGS unmatched cases (n = 32) | P-value |
Hospital length of stay (days) | 31 (22–52) | 30 (22–51) | 0.510 |
ICU length of stay (days) | 20 (11–28) | 19 (11–28) | 0.625 |
Mortality (%) | 10 (35.71%) | 13 (40.61%) | 0.696 |
The initially univariate regression analysis identified potential factors influencing sepsis prognosis. The univariate Cox regression analysis indicated a significant association between age (adjusted HR = 0.98) and CMT positive (adjusted HR = 2.36) with adverse outcomes (P ≦ 0.05) (Table 5). To conduct a comprehensive assessment of independent risk factors and determine the impact of matched dual mNGS results on patient prognosis, we incorporated ‘dual mNGS consistency’ and variables with P < 0.2 in the univariate analysis into the multivariate analysis. The multivariate analysis revealed that NE% emerged as an independent risk factor for poor prognosis among sepsis patients (adjusted HR = 0.94, P = 0.002) (Table 5). The consistency of dual-mNGS exhibited a trend as a risk factor for mortality in sepsis patients (adjusted HR = 1.35, P = 0.451).
Table 5
Multivariate regression analysis of predictors of clinical outcomes.
Variables | Univariate analyses | Multivariate analysis |
HR (95% CI) | P-value | HR (95% CI) | P-value |
Gender | 0.77 (0.38–1.54) | 0.457 | —— | —— |
Age | 0.98 (0.97-1.00) | 0.048* | 0.99 (0.96–1.01) | 0.256 |
SOFA | 1.00 (0.94–1.08) | 0.914 | —— | —— |
RF | 0.70 (0.36–1.38) | 0.307 | —— | —— |
Dual mNGS consistent | 0.94 (0.49–1.82) | 0.860 | 1.35 (0.62–2.93) | 0.451 |
CMT positive | 2.36 (1.02–5.46) | 0.044* | 0.54 (0.15-2.00) | 0.358 |
WBC | 1.02 (0.98–1.08) | 0.332 | —— | —— |
Hb | 1.01 (1.00-1.02) | 0.156 | 1.01 (0.99–1.03) | 0.194 |
PLT | 1.00 (1.00–1.00) | 0.356 | —— | —— |
NE% | 0.98 (0.96-1.00) | 0.075 | 0.94 (0.91–0.98) | 0.002** |
HCT | 0.99 (0.97–1.01) | 0.270 | —— | —— |
CRP | 1.00 (1.00–1.00) | 0.747 | —— | —— |
PCT | 1.01 (1.00-1.02) | 0.072 | 1.01 (1.00-1.02) | 0.064 |
HR, hazard ratio; CI, confidence interval; SOFA, sequential organ failure assessment; CMT, conventional microbiological tests. |