Study population and baseline characteristics
From an initial screening of 312 ICU admissions, 118 patients met the inclusion criteria and were included in the final analysis cohort. These patients were classified into four pathogen groups: Legionella pneumophila (n = 18), other gram-negative bacteria (n = 68), gram-positive bacteria (n = 24), and fungi (n = 8) (Figure 1). Baseline characteristics were generally comparable across groups (Table 1). Patients with Legionella pneumophila infection showed a higher prevalence of immunosuppression compared with those with other gram-negative infections (22.2% vs 7.4%, p = 0.048). Clinically, Legionella patients more frequently presented with fever (83.3%), cough (100%), and altered mental status (88.9%). Illness severity, assessed by the APACHE II score, was broadly similar across groups, although Legionella cases tended to have slightly higher median scores (28 [IQR 24–36] vs 25 [IQR 18–31] in other gram-negative infections) as shown in Table 1.
Endotoxin levels and diagnostic performance
Infections with Legionella pneumophila were associated with substantially higher endotoxin levels compared with other pathogen groups. The mean endotoxin level was 2.15 ± 1.28 EU/mL in Legionella infections, compared with 0.61 ± 0.74 EU/mL in other gram-negative infections, 0.42 ± 0.58 EU/mL in gram-positive infections, and 0.13 ± 0.15 EU/mL in fungal infections (Figure 2; Table 2). Using a threshold of > 2.5 EU/mL, endotoxin quantification yielded a sensitivity of 77.8% (95% CI: 52.4–93.6) and a specificity of 89.7% (95% CI: 79.9–95.8) for differentiating Legionella infections from other gram-negative infections. The positive predictive value was 66.7% (95% CI: 43.0–85.4) and the negative predictive value was 95.3% (95% CI: 87.1–99.0). Receiver operating characteristic analysis demonstrated good diagnostic performance, with an area under the curve (AUC) of 0.882 (95% CI: 0.809–0.955) (Figure 3; Table 2).
Laboratory biomarkers and clinical correlations
Patients with Legionella pneumophila infections exhibited distinctive laboratory profiles compared with other pathogen groups. Notably, serum sodium levels were significantly lower in Legionella cases (133 ± 9.8 mmol/L) than in other gram-negative infections (142 ± 11.4 mmol/L, p = 0.003). Procalcitonin levels also differed across groups (p = 0.038), with Legionella infections showing a median level of 28.6 ng/mL (IQR 5.2–98.5) compared with 21.3 ng/mL (IQR 0.08–95.6) in other gram-negative infections. White blood cell counts and platelet counts tended to be lower in Legionella infections (9.24 ± 7.15 × 103/µL and 186 ± 78 × 103/µL, respectively) than in other gram-negative infections (14.2 ± 10.3 × 103/µL and 198 ± 92 × 103/µL), although these differences did not reach statistical significance. These laboratory patterns complemented the endotoxin findings described above (Table 3).
Clinical outcomes and endotoxin stratification
Stratification by endotoxin levels revealed that patients with endotoxin >2.5 EU/mL had greater illness severity and distinct clinical characteristics compared with those with lower endotoxin levels. The high-endotoxin group was strongly enriched for Legionella pneumophila infections (66.7% vs 4.1%, p < 0.001) and received Legionella-active antimicrobial therapy more frequently (76.2% vs 18.6%, p < 0.001). Patients in the high-endotoxin group also had higher illness severity, with a median APACHE II score of 28 (IQR 24–36) compared with 25 (IQR 18–31) in the low-endotoxin group (p = 0.042). Mechanical ventilation was more common among patients with endotoxin >2.5 EU/mL (76.2% vs 58.8%, p = 0.042). Regarding clinical outcomes, 28-day mortality was significantly higher in the high-endotoxin group (38.1% vs 24.7%, log-rank p = 0.035), and ICU length of stay was longer (9 vs 6 days, p = 0.024). Cox proportional hazards analysis demonstrated that endotoxin >2.5 EU/mL was independently associated with increased mortality risk (hazard ratio 2.14, 95% CI 1.08–4.23, p = 0.028) (Figure 4; Table 4).
Multivariate analysis and model performance
Univariable logistic regression analysis identified several factors associated with Legionella pneumophila infection, including endotoxin >2.5 EU/mL (odds ratio [OR] 30.2, 95% CI 8.71–104.7, p < 0.001), hyponatremia <135 mmol/L (OR 5.84, 95% CI 1.98–17.2, p = 0.001), presence of cough (OR 18.4, 95% CI 2.31–146.3, p = 0.006), altered mental status (OR 6.73, 95% CI 1.89–23.9, p = 0.003), and immunosuppression (OR 3.56, 95% CI 1.04–12.2, p = 0.043). In multivariable logistic regression, endotoxin >2.5 EU/mL remained the strongest independent predictor of Legionella infection (adjusted OR 15.7, 95% CI 3.84–64.2, p < 0.001), followed by hyponatremia <135 mmol/L (adjusted OR 4.82, 95% CI 1.26–18.4, p = 0.021) and presence of cough (adjusted OR 12.3, 95% CI 1.34–112.6, p = 0.026). Altered mental status showed a trend toward association but did not reach statistical significance after adjustment (adjusted OR 3.41, 95% CI 0.89–13.1, p = 0.073). The combined prediction model integrating endotoxin levels with clinical variables demonstrated excellent discrimination, with an area under the receiver operating characteristic curve (AUC) of 0.938 (95% CI 0.882–0.994). Model sensitivity and specificity were 83.3% and 91.2%, respectively, with good calibration indicated by the Hosmer–Lemeshow test (χ2 = 6.84, p = 0.553). Comparative analysis showed that the combined endotoxin–clinical model outperformed both endotoxin alone (AUC 0.882) and clinical variables alone (AUC 0.784). DeLong testing confirmed significant improvements in discrimination for the combined model compared with endotoxin alone (p = 0.024) and clinical variables alone (p < 0.001) (Figure 5; Figure 6; Table 5).
DATA AVAILABILITY:
All the raw data generated and analyzed during the current study are available in the supplementary materials.

Figure 1: Study flow chart. Flowchart illustrating patient selection for the retrospective cohort study. A total of 312 ICU admissions between January 2020 and December 2023 were screened. 194 patients were excluded, including those with incomplete microbiological data (n = 67), missing endotoxin measurements (n = 89), polymicrobial infections (n = 24), and early deaths within 24 h (n = 14). The final analysis cohort included 118 patients, who were stratified by pathogen group: Legionella pneumophila (n = 18), other gram-negative bacteria (n = 68), gram-positive bacteria (n = 24), and fungi (n = 8). Please click here to view a larger version of this figure.

Figure 2: Serum endotoxin levels by pathogen group. Box plot comparison of serum endotoxin levels across different pathogen groups. Legionella pneumophila infections demonstrated significantly elevated endotoxin levels compared to other gram-negative bacteria, gram-positive bacteria, and fungal infections. The horizontal dashed line indicates the optimal diagnostic cutoff of 2.5 EU/mL. Individual data points are overlaid to show the distribution within each group. Please click here to view a larger version of this figure.

Figure 3: ROC curve analysis for Legionella detection. The area under the curve (AUC) was 0.882 (95% CI: 0.809–0.955), indicating good diagnostic performance. The optimal diagnostic cutoff was 2.5 EU/mL, corresponding to a sensitivity of 77.8% and specificity of 89.7%. The optimism-corrected AUC was 0.874. Please click here to view a larger version of this figure.

Figure 4: Survival analysis by endotoxin levels. Kaplan–Meier survival curves comparing patients with endotoxin levels >2.5 EU/mL (n = 21) and ≤2.5 EU/mL (n = 97). Patients with endotoxin levels >2.5 EU/mL showed significantly lower survival probability, with 28-day mortality rates of 38.1% versus 24.7%. The difference in survival was statistically significant (log-rank p = 0.035). Cox proportional hazards analysis demonstrated that endotoxin levels >2.5 EU/mL were associated with increased mortality risk (hazard ratio 2.14, 95% CI 1.08–4.23, p = 0.028). Numbers at risk are shown below the survival curves. Please click here to view a larger version of this figure.

Figure 5: Multivariate analysis forest plot. Adjusted odds ratios (ORs) with 95% confidence intervals are shown for candidate predictors. Endotoxin levels >2.5 EU/mL were the strongest independent predictor (adjusted OR 15.7, 95% CI 3.84–64.2, p < 0.001). Hyponatremia <135 mmol/L (adjusted OR 4.82, 95% CI 1.26–18.4, p = 0.021) and presence of cough (adjusted OR 12.3, 95% CI 1.34–112.6, p = 0.026) were also independently associated with Legionella infection, whereas altered mental status and immunosuppression were not statistically significant. Please click here to view a larger version of this figure.

Figure 6: Comparative ROC analysis. The combined model incorporating endotoxin levels, hyponatremia, and clinical symptoms demonstrated the highest diagnostic performance (AUC = 0.938), outperforming endotoxin alone (AUC = 0.882) and clinical variables alone (AUC = 0.784). Pairwise comparisons using the DeLong test showed that the combined model significantly improved discrimination compared with endotoxin alone (p = 0.024) and clinical variables (p < 0.001). Please click here to view a larger version of this figure.
| Characteristic | Legionella pneumophila (n = 18) | Other Gram-negative (n = 68) | Gram-positive (n = 24) | Fungi (n = 8) | P-value |
| Demographics | | | | | |
| Age, median (IQR), years | 67 (5275) | 67 (56–81) | 63 (48–74) | 69 (60–83) | 0.678 |
| Male sex, n (%) | 13 (72.2) | 43 (63.2) | 19 (79.2) | 6 (75.0) | 0.512 |
| Comorbidities | | | | | |
| Diabetes mellitus, n (%) | 10 (55.6) | 27 (39.7) | 9 (37.5) | 3 (37.5) | 0.534 |
| Hypertension, n (%) | 8 (44.4) | 33 (48.5) | 10 (41.7) | 4 (50.0) | 0.932 |
| COPD, n (%) | 3 (16.7) | 9 (13.2) | 3 (12.5) | 2 (25.0) | 0.778 |
| Immunosuppression, n (%) | 4 (22.2) | 5 (7.4) | 0 (0) | 1 (12.5) | 0.048 |
| Clinical presentation | | | | | |
| Fever >38 °C, n (%) | 15 (83.3) | 43 (63.2) | 12 (50.0) | 1 (12.5) | 0.012 |
| Cough, n (%) | 18 (100) | 47 (69.1) | 10 (41.7) | 3 (37.5) | <0.001 |
| Altered mental status, n (%) | 16 (88.9) | 38 (55.9) | 15 (62.5) | 3 (37.5) | 0.042 |
| Severity scores | | | | | |
| APACHE II, median (IQR) | 28 (24–36) | 25 (18-31) | 26 (21–32) | 25 (18–34) | 0.168 |
| Outcomes | | | | | |
| Overall mortality, n (%) | 7 (38.9) | 24 (35.3) | 6 (25.0) | 3 (37.5) | 0.689 |
Table 1: Baseline patient characteristics by pathogen group. Abbreviations; APACHE II = Acute Physiology and Chronic Health Evaluation II; COPD = Chronic obstructive pulmonary disease; IQR = Interquartile range. P-values calculated using Kruskal-Wallis test for continuous variables and Fisher's exact test for categorical variables. Bold values indicate statistical significance (p < 0.05).
| Parameter | Legionella pneumophila (n = 18) | Other Gram-negative (n = 68) | Gram-positive (n = 24) | Fungi (n = 8) |
| Endotoxin levels (EU/mL) | | | | |
| Mean ± SD | 2.15 ± 1.28 | 0.61 ± 0.74 | 0.42 ± 0.58 | 0.13 ± 0.15 |
| Median (IQR) | 2.50 (1.35–2.80) | 0.32 (0.15–0.78) | 0.28 (0.12–0.52) | 0.08 (0.04–0.16) |
| Endotoxin >2.5 EU/mL, n (%) | 14 (77.8) | 7 (10.3) | 0 (0) | 0 (0) |
| Diagnostic performance (>2.5 EU/mL threshold)* | | | | |
| Sensitivity, % (95% CI) | 77.8 (52.4–93.6) | — | — | — |
| Specificity, % (95% CI) | 89.7 (79.9–95.8) | — | — | — |
| Positive predictive value, % (95% CI) | 66.7 (43.0–85.4) | — | — | — |
| Negative predictive value, % (95% CI) | 95.3 (87.1–99.0) | — | — | — |
| Positive likelihood ratio (95% CI) | 7.55 (3.84–14.8) | — | — | — |
| Negative likelihood ratio (95% CI) | 0.25 (0.11–0.56) | — | — | — |
| Area under ROC curve (95% CI) | 0.882 (0.809–0.955) | — | — | — |
| Optimism-corrected AUC | 0.874 (0.798–0.947) | — | — | — |
Table 2: Endotoxin levels and diagnostic performance metrics. For differentiating Legionella pneumophila from other gram-negative bacterial infections. Abbreviations; AUC = area under the receiver operating characteristic curve; CI = confidence interval; EU = Endotoxin units; IQR = Interquartile range; ROC: Receiver operating characteristic; SD = Standard deviation. Optimism-corrected AUC calculated using bootstrap resampling with 1000 iterations.
| Biomarker | Legionella (n = 18) | Other Gram-negative (n = 68) | Gram-positive (n = 24) | Fungi (n = 8) | P-value |
| White blood cell count, ×10³/μL | 9.24 ± 7.15 | 14.2 ± 10.3 | 15.1 ± 9.8 | 13.2 ± 8.6 | 0.089 |
| Platelet count, ×10³/μL | 186 ± 78 | 198 ± 92 | 205 ± 88 | 192 ± 76 | 0.823 |
| Procalcitonin, ng/mL | 28.6 (5.2–98.5) | 21.3 (0.08–95.6) | 9.4 (0.12–34.8) | 0.18 (0.06–0.32) | 0.038 |
| C-reactive protein, mg/L | 172.4 ± 61.2 | 135.7 ± 68.9 | 128.3 ± 87.6 | 102.8 ± 42.3 | 0.312 |
| Serum sodium, mmol/L | 133 ± 9.8 | 142 ± 11.4 | 143 ± 7.2 | 141 ± 6.8 | 0.003 |
| Serum creatinine, mg/dL | 1.58 ± 0.92 | 1.32 ± 0.86 | 1.21 ± 0.74 | 1.45 ± 0.88 | 0.142 |
| Lactate, mmol/L | 2.8 ± 1.4 | 2.4 ± 1.6 | 2.2 ± 1.3 | 2.6 ± 1.5 | 0.456 |
| Bilirubin, mg/dL | 1.2 ± 0.8 | 1.1 ± 0.9 | 0.9 ± 0.6 | 1.3 ± 0.7 | 0.678 |
Table 3: Laboratory biomarkers and inflammatory parameters. Values are presented as mean ± standard deviation or median (interquartile range). P-values calculated using analysis of variance or Kruskal-Wallis test as appropriate. Bold values indicate statistical significance (p < 0.05).
| Outcome | Endotoxin >2.5 EU/mL (n = 21) | Endotoxin ≤2.5 EU/mL (n = 97) | P-value |
| Pathogen distribution | | | |
| Legionella pneumophila, n (%) | 14 (66.7) | 4 (4.1) | <0.001 |
| Other gram-negative, n (%) | 7 (33.3) | 61 (62.9) | 0.012 |
| Gram-positive, n (%) | 0 (0) | 24 (24.7) | 0.008 |
| Fungi, n (%) | 0 (0) | 8 (8.2) | 0.198 |
| Severity and management | | | |
| APACHE II score, median (IQR) | 28 (24–36) | 25 (18–31) | 0.042 |
| Legionella-active therapy, n (%) | 16 (76.2) | 18 (18.6) | <0.001 |
| Mechanical ventilation, n (%) | 16 (76.2) | 57 (58.8) | 0.042 |
| Vasopressor requirement, n (%) | 13 (61.9) | 46 (47.4) | 0.128 |
| Clinical outcomes | | | |
| 7-day mortality, n (%) | 3 (14.3) | 8 (8.2) | 0.389 |
| 28-day mortality, n (%) | 8 (38.1) | 24 (24.7) | 0.048 |
| Hospital mortality, n (%) | 9 (42.9) | 31 (32.0) | 0.256 |
| ICU length of stay, days, median (IQR) | 9 (6–14) | 6 (4–10) | 0.024 |
| Hospital length of stay, days, median (IQR) | 18 (12–28) | 15 (10–24) | 0.142 |
| Ventilator-free days at 28 days, median (IQR) | 12 (0–22) | 18 (8–25) | 0.089 |
| Hazard ratio for mortality (adjusted)* | 2.14 (1.08–4.23) | Reference | 0.028 |
Table 4: Clinical outcomes stratified by endotoxin levels. Adjusted for APACHE II score, age, and Charlson Comorbidity Index using Cox proportional hazards regression. Abbreviations; APACHE II = acute physiology and chronic health evaluation II; EU = Endotoxin units; ICU = intensive care unit; IQR = interquartile range. Legionella-active therapy includes fluoroquinolones or macrolides. Bold values indicate statistical significance (p < 0.05).
| Analysis | Value | P-value |
| Univariable logistic regression | | |
| Endotoxin >2.5 EU/mL | OR 30.2 (95% CI: 8.71–104.7) | <0.001 |
| Hyponatremia <135 mmol/L | OR 5.84 (95% CI: 1.98–17.2) | 0.001 |
| Cough present | OR 18.4 (95% CI: 2.31–146.3) | 0.006 |
| Altered mental status | OR 6.73 (95% CI: 1.89–23.9) | 0.003 |
| Immunosuppression | OR 3.56 (95% CI: 1.04–12.2) | 0.043 |
| Multivariable logistic regression* | | |
| Endotoxin >2.5 EU/mL | Adjusted OR 15.7 (95% CI: 3.84–64.2) | <0.001 |
| Hyponatremia <135 mmol/L | Adjusted OR 4.82 (95% CI: 1.26-18.4) | 0.021 |
| Cough present | Adjusted OR 12.3 (95% CI: 1.34–112.6) | 0.026 |
| Altered mental status | Adjusted OR 3.41 (95% CI: 0.89-13.1) | 0.073 |
| Immunosuppression | Adjusted OR 2.18 (95% CI: 0.51–9.32) | 0.291 |
| Model performance metrics | | |
| Combined model AUC (95% CI) | 0.938 (0.882–0.994) | — |
| Combined model sensitivity, % (95% CI) | 83.3 (58.6–96.4) | — |
| Combined model specificity, % (95% CI) | 91.2 (81.8–96.7) | — |
| Hosmer-Lemeshow goodness-of-fit | χ² = 6.84, df = 8 | 0.553 |
| Comparison of biomarker performance (AUC) | | |
| Endotoxin alone | 0.882 (95% CI: 0.809–0.955) | <0.001 |
| Clinical variables alone | 0.784 (95% CI: 0.691–0.877) | <0.001 |
| Combined endotoxin + clinical model | 0.938 (95% CI: 0.882–0.994) | <0.001 |
| Statistical comparisons (DeLong test) | | |
| Combined vs Endotoxin alone | — | 0.024 |
| Combined vs Clinical variables alone | — | <0.001 |
| Endotoxin vs Clinical variables alone | — | 0.018 |
Table 5: Multivariate analysis and ROC curve performance. Variables with p < 0.10 in univariable analysis were included in the multivariable model. Variance inflation factors for all variables < 2.1, indicating acceptable multicollinearity. Abbreviations; AUC = Area under the receiver operating characteristic curve; CI = confidence interval; df = degrees of freedom; EU = Endotoxin units; OR = Odds ratio. Bold values indicate statistical significance (p < 0.05). The combined model includes endotoxin > 2.5 EU/mL, hyponatremia, and cough as independent predictors.