According to the specified criteria, we documented 197 episodes of VAP in 129 patients, as depicted in Fig. 1.
The median age of the patients included was 44 years (IQR 35–56), 31% were women (n = 40), and the median Charlson comorbidity index was 1 (IQR 0–2). The median SOFA score was 5 (IQR 3–7), vasopressors were required for 40.60% of the VAP episodes (n = 80), and 26.39% of the episodes met the definition for ARDS (n = 52). The median hospital stay after VAP development was 12 days (IQR 6–21), the 90th percentile of days after VAP development was 29.4 days, and the in-hospital mortality rate for patients with VAP was 34.51% (n = 68). The remaining sociodemographic and clinical characteristics are shown in Table 1.
In terms of local epidemiology and antibiotic resistance, the most frequently isolated organisms were Acinetobacter baumannii (36.04%, n = 71), Enterobacterales (26.39%, n = 52), Pseudomonas aeruginosa (18.25%, n = 36), Staphylococcus aureus (10.15%, n = 20), and Stenotrophomonas maltophilia (8.12%, n = 16). In 41.12% of the episodes, more than one organism was isolated (n = 81). The remaining isolated organisms can be found in Table 1.
According to standardized international terminology for describing acquired resistance profiles, (28) 87.32% of Acinetobacter baumannii isolates were classified as XDR (extensively drug-resistant) Acinetobacter baumannii (XDR-AB) (n = 62). For Enterobacterales, 28.3% of the isolates were MDR (multidrug resistant) Enterobacteriaceae (n = 15), and 7.5% were carbapenem-resistant Enterobacteriaceae (n = 4). For Pseudomonas aeruginosa, 19.44% of the isolates were MDR, and 5.55% were XDR. Only 1% of the Staphylococcus aureus isolates were methicillin-resistant Staphylococcus aureus (MRSA) (n = 2). The rest of the resistance profiles can be found in Table 2.
The most frequently used initial empiric therapy regimens were meropenem (24.87%, n = 49), colistin + tigecycline (18.27%, n = 36), piperacillin/tazobactam (13.7%, n = 27), and meropenem + linezolid (7.61%, n = 15). Other regimens were used in 35.53% (n = 70) of the patients. The remaining antimicrobial regimens are listed in Table S1 of the supplementary material. Initial empiric therapy was inappropriate for 45.18% of VAP episodes (n = 89). However, in 85.79% of the events, patients eventually received an appropriate regimen during the entire VAP episode. The median number of days to correct the initial inappropriate empiric treatment was 3 days (IQR 1-5.75) (Table 1).
Eight antibiotics were tested as monotherapies, and thirteen combined regimens were studied in the WISCA model for the overall pool of pathogens. The median estimated coverages ranged from 4.99% (vancomycin, VA) to 87.23% (cefepime + tigecycline, FEP + TGC). Tigecycline-based combined regimens showed greater median coverage (+ 47.97%) than nontigecycline-based combined regimens. The complete set of medians of the posterior distribution and the associated 95% highest density intervals (HDIs) are shown in Fig. 2.
The algorithm used to sample from the parameters posterior distributions achieved optimal convergence, with R^ index values always near 1 and good mixing for all chains, as shown in Table S2 of the supplementary material. The Monte Carlo Markov chain (MCMC) trace plots of the WISCA model parameters and the density plots of the posterior distributions of the model parameters are depicted in Figures S1 and S2 of the supplementary material.
We found significant differences in the distributions of VAP-causing microorganisms according to the days of invasive mechanical ventilation (IMV) prior to the development of VAP. Specifically, Enterobacterales and Staphylococcus aureus were more frequently isolated within the first 5 days of IMV, whereas Acinetobacter baumannii was more frequently isolated after 5 days of IMV (p < 0.02), as shown in Fig. 3 and Table 2.
Due to these differences in the distributions of the isolated microorganisms, we calculated the posterior coverages of the thirteen most commonly used antibiotic regimens for both time periods, as shown in Fig. 4. The complete set of medians of the posterior distribution and the associated 95% HDIs for these antibiotic regimens during these time periods are shown in Table S3, as well as the R^ index values of the WISCA models parameters in Table S4 of the supplementary material.
According to the univariate analyses, the presence of a low PaO2/FiO2 ratio, inappropriate treatment for VAP episodes, a high SOFA score, the need for vasopressors, the presence of ARDS, and Acinetobacter baumannii infection were positively associated with mortality. The Charlson comorbidity index and days of IMV after VAP development were positively associated with an extended hospital stay. The remaining sociodemographic and clinical variables can be found in Table S5 of the supplementary material.
According to the multivariate regression analysis, inappropriate treatment for the entire VAP event, a higher SOFA score, the need for vasopressors, the presence of ARDS, and Acinetobacter baumannii infection were independently positively associated with mortality. Interestingly, an increase in days of IMV after VAP development was negatively associated with mortality. On the other hand, initial empiric inappropriate treatment and an increase in days of IMV after VAP development were positively associated with an extended hospital stay. The goodness of fit was adequate for both models (Hosmer‒Lemeshow test, p = 0.988 and 0.925, respectively). The complete set of variables for both models is depicted in Table 3.