Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol. , 12 March 2025

Sec. Phage Biology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1530819

This article is part of the Research Topic Harnessing Bacteriophages and Phage-Engineered Products for Antibacterial and Anticancer Therapies: Challenges and Opportunities View all 5 articles

Phage-antibiotic combinations against Klebsiella pneumoniae: impact of methodological approaches on effect evaluation

  • 1Department of Biomedicine and Genomics, Lopukhin Federal Research and Clinical Center of Physical-Chemical Medicine of Federal Medical Biological Agency, Moscow, Russia
  • 2Moscow Center for Advanced Studies, Moscow, Russia

Background: The combined use of bacteriophages and antibiotics represents a promising strategy for combating multidrug-resistant bacterial pathogens. However, the lack of uniformity in methods for assessing combination effects and experimental protocols has resulted in inconsistent findings across studies. This study aimed to evaluate the effects of interactions between phages and antibiotics on Klebsiella pneumoniae strains using various statistical approaches to formalize combination effects.

Methods: Effects were assessed for four antibiotics from distinct classes (gentamicin, levofloxacin, meropenem, chloramphenicol), three phages from different genera (Dlv622, Seu621, FRZ284), and a depolymerase (Dep622) on three K. pneumoniae strains of the KL23 capsule type. Antibiotics were used at Cmax concentrations, and phages at sublethal levels. A modified t-test, Bliss independence model, two-way ANOVA, and checkerboard assay were employed to evaluate the results.

Results: Among 48 combinations, 33 effects were statistically significant, including 26 cases of synergy and 7 of antagonism. All statistical methods showed consistency in identifying effects; however, the t-test and Bliss method detected a greater number of effects. The strongest synergy was observed with levofloxacin in combination with Seu621 or Dep622 across all bacterial strains. Checkerboard assays confirmed synergy in selected cases but indicated that combined effects could vary with antimicrobial concentrations.

Conclusion: The choice of analytical method substantially impacts the detection of phage-antibiotic effects. The t-test and Bliss method, due to their simplicity and sensitivity, may be optimal for clinical application, while two-way ANOVA for confirming strong interactions. These results emphasize the need to consider interaction characteristics when designing therapeutic strategies.

1 Introduction

Klebsiella pneumoniae is both a natural component of human microbiota, colonizing the gastrointestinal tract and nasopharynx in healthy individuals (Chang et al., 2021), and a dangerous pathogen, responsible for approximately 11–15% of hospital-acquired pneumonia cases and up to 8% of all healthcare-associated infections (Chen et al., 2023; Grosjean et al., 2024; Rahmat Ullah et al., 2024).

Infections caused by K. pneumoniae represent a significant concern in healthcare settings due to the increasing prevalence of antibiotic-resistant strains. In Europe, over one third of K. pneumoniae isolates demonstrate resistance to at least one major antimicrobial drug class, including fluoroquinolones, third-generation cephalosporins, aminoglycosides, or carbapenems (European Centre for Disease Prevention and Control, 2022). These infections are responsible for an estimated 650,000 deaths annually worldwide (Murray et al., 2022). In light of this significant public health concern, the World Health Organization (WHO) has identified K. pneumoniae as a high-priority pathogen in its 2024 Bacterial Priority Pathogens List (World Health Organization, 2024).

One promising alternative for combating infections caused by antibiotic-resistant strains is bacteriophage (phage) therapy and its derivatives, including polysaccharide depolymerases and endolysins (Danis-Wlodarczyk et al., 2021; Durbas and Machnik, 2022). Bacteriophages—viruses that specifically infect and lyse bacteria—have demonstrated promising results in treatment applications (Pirnay et al., 2024). Although clinical data on recombinant phage proteins in humans are still unavailable, their efficacy has been convincingly demonstrated in animal models (Volozhantsev et al., 2020; Gorodnichev et al., 2021a; Hua et al., 2022; Li et al., 2024).

Phage therapy is commonly implemented under Article 37 of the Declaration of Helsinki, providing “compassionate care” when conventional antimicrobial therapies fail. However, standard treatment protocols generally do not recommend antibiotic withdrawal, meaning that phage therapy is often administered concurrently with antibiotics. This raises important questions about the potential interactions between these antimicrobial agents.

The combined use of phages and antibiotics can yield varying outcomes: synergy, additive effect, or antagonism. Synergy occurs when the bactericidal action of antimicrobial agents together is greater than the sum of their individual effects. Additive effects indicate an impact equal to the sum of the individual effects, while antagonism describes an interaction where the action of one agent inhibits the effect of the other (Liu et al., 2020).

In therapeutic contexts, antagonistic interactions between phages and antibiotics, which can reduce overall treatment efficacy, are considered particularly critical—even more so than synergy among antimicrobial agents. Reflecting this, several phage therapy centers have incorporated in vitro tests for detecting combined effects of phages and antibiotics as part of the standard phage selection process (Onallah et al., 2023).

To date, no cases of antagonism between phages and antibiotics targeting K. pneumoniae have been documented. Instead, both in vitro (Bedi et al., 2009; Verma et al., 2010; Qurat-ul-Ain et al., 2021; Han et al., 2022; Xu et al., 2022; Lew et al., 2023) and in vivo (Chhibber et al., 2008; Pacios et al., 2021; Wang et al., 2021) studies have demonstrated high efficacy of these agents when used together. However, differences in bacterial strains, phages, antibiotics, and methodologies often complicate comparisons across studies.

The aim of this work was to compare various in vitro methods for assessing combination effects of antibiotics and phages. Using three strains of K. pneumoniae, we evaluated the individual and combined actions of four antibiotics (gentamicin, levofloxacin, meropenem, and chloramphenicol), three phages (Dlv622, Seu621, and FRZ284), and a depolymerase (Dep622). To quantify combination effects, we applied the most commonly used methods for such studies: the t-test, Bliss independence model, two-way ANOVA, and checkerboard assay.

2 Materials and methods

2.1 Bacterial strains and characteristics

This study utilized K. pneumoniae strains 9faize, G4-17, and B536-17-2, isolated from patients between 2018 and 2023. Prior to use, the strains were propagated in lysogeny broth (LB, HiMedia, India) at 37°C. Capsule typing of K. pneumoniae was conducted via sequencing of the wzi gene (Brisse et al., 2013), and sequence typing followed the standard MLST protocol (Diancourt et al., 2005). Antibiotic susceptibility was determined using the broth microdilution method according to Clinical & Laboratory Standards Institute (CLSI) guidelines (CLSI, 2015). The antibiotics used included gentamicin (GEN), levofloxacin (LVX), meropenem (MEM), and chloramphenicol (CMP) (HiMedia, India).

2.2 Bacteriophages

Three previously characterized K. pneumoniae phages Dlv622, Seu621, and FRZ284 (GenBank accessions MT939252, MT939253, and MZ602148) were included in the study, classified under the genera Drulisvirus, Mydovirus, and Jiaodavirus, respectively (Gorodnichev et al., 2021a; b). Phage titers and efficiency of plating (EOP) were determined using the spot test and double-layer agar methods (Kutter, 2009; Mazzocco et al., 2009).

Lytic activity of the phages was assessed using Appelman’s titration method (Burrowes et al., 2019). For this, phage suspensions were serially diluted 10-fold in LB broth. A total of 190 μL of each dilution was added to microplate wells, followed by inoculation with 10 μL of bacterial culture in the logarithmic growth phase at a concentration of 105 CFU/well. The plates were incubated overnight at 37°C, and lytic activity was determined as the lowest dilution that inhibited bacterial growth.

2.3 Polysaccharide depolymerase

The recombinant polysaccharide depolymerase Dep622, derived from phage Dlv622, was previously obtained and characterized (Gorodnichev et al., 2021a). The minimum halo-forming concentration (MHF) of depolymerase was assessed using a method similar to phage titration. For the assay, 7 mL of molten, semi-solid LB agar (0.7%) was combined with 100 μL of an overnight culture of the test strain, vortexed, and spread evenly across the surface of a solidified LB agar plate. After solidification, 5 μL of serially diluted recombinant depolymerase (prepared in a 2-fold LB dilution series) was applied to the agar surface using a template. The plate was incubated at 37°C for 18 h, and translucent spots were observed on the bacterial lawn. The MHF was defined as the lowest concentration that produced visible spots.

2.4 Individual and combined effects of antimicrobial agents on planktonic bacterial cells

The individual and combined effects of antibiotics, phages, and depolymerase were assessed using time-kill curves. Antimicrobial agents were added individually or in combination to the wells of a 96-well plate. Antibiotics were applied at concentrations equivalent to their maximum serum concentration (Cmax) (Table 1). Phage lysates were used at a sublethal concentration based on Appelman’s titration (109 PFU/well), while recombinant depolymerase at 100 MHF. A positive bacterial growth control included no antimicrobial agents, while control wells for individual agents contained only the antibiotic, phage, or depolymerase. Wells were inoculated with 10 μL of each strain in the exponential growth phase (OD600 = 0.3), diluted in LB to achieve 105 CFU/well. The final volume in each well was 200 μL. Optical density (OD) was measured every hour for 20 h at 37°C using a FlexA200 microplate reader (Allsheng, China) at 600 nm.

Table 1
www.frontiersin.org

Table 1. MIC breakpoints for K. pneumoniae and maximum concentration of antibiotic in serum.

The individual effects of antibacterial agents on planktonic cells were evaluated by calculating the percent reduction in optical density, according to the formula: (1 - ((OD of treated - ODK−)/(OD of untreated - ODK−))) × 100%. The endpoint and area under the curve (AUC) data were used for the calculation. The AUC was determined using the trapezoidal rule.

Statistical analyses, including t-tests, two-way ANOVA, and Bliss analysis, were used to evaluate combined effects.

Based on endpoint OD600 measurements, a t-test was utilized to assess differences between OD600 values of the antimicrobial combinations and their individual actions. Antagonism was indicated when the combination’s OD600 was significantly higher than that of the most effective agent alone, or when the OD600 of the combination did not differ from the most effective agent, while the second agent showed a significant effect compared to the control. Synergy was defined when the combination’s OD600 was significantly lower than the most effective agent, while the second agent’s OD600 was not different from the control. When the OD600 of the combination was significantly lower than that of the most effective agent, and both individual agents also showed a significant reduction compared to the control, this phenomenon was termed “tentative synergy.”

Bliss analysis followed the approach of Eller et al. (2021). From endpoint OD600 data, a synergy coefficient (S) and p-value were calculated, where S > 0 indicated synergy and S < 0 indicated antagonism. The significance of the S value was assessed using a one-sample t-test with a hypothesized mean of zero, accounting for individual propagated errors.

Two-way ANOVA was conducted following Gordillo Altamirano et al. (2022) for both endpoint OD600 and the area under the curve, calculated using the trapezoidal rule. Data normality was verified using the Shapiro–Wilk test, and multiple comparisons were corrected with Tukey’s test. Interaction plot slopes were used to interpret types of combined effects.

Statistical analyses were performed using GraphPad Prism 8. Any effects not meeting the criteria for statistical significance (p < 0.05) were considered additive.

2.5 Checkerboard assay

A modified checkerboard assay (Nikolic et al., 2022) was performed for specific combinations of antimicrobial agents. Antibiotic concentrations ranged from 2 to 128 μg/mL across columns, while phage titers varied from 103 to 109 PFU/well across rows. The first row contained only phages, and the last row included only antibiotics. The bottom-left well served as a positive control with no antimicrobial agents. Each well was inoculated with 10 μL of bacterial culture in the exponential growth phase (OD600 = 0.3), diluted in LB to a final concentration of 105 CFU/well.

The plates were incubated for 20 h at 37°C, and OD was measured hourly at 600 nm using a FlexA200 microplate reader (Allsheng, China). Each experiment was conducted in triplicate, with results averaged.

To evaluate the effect of phage-antibiotic interactions, the Fractional Inhibitory Concentration Index (FICI) was calculated:

FICI = MIC of phage in combination MIC of phage alone + MIC of antibiotic in combination MIC of antibiotic alone

FICI values were interpreted as follows: FICI <0.5 indicated synergy, FICI between 0.5 and 1 indicated an additive effect, 1 < FICI <2 was considered indifferent, and FICI >2 indicated antagonism. Data visualization was performed using ComplexHeatmap v2.16.0 for R v4.3.0 and SynergyFinder v3.0 (Gu et al., 2016; Ianevski et al., 2022).

3 Results

3.1 Susceptibility of Klebsiella pneumoniae strains to antibiotics, phages, and depolymerase

The K. pneumoniae strains were isolated from clinical samples between 2018 and 2023, belonged to capsule type KL23 and clinically relevant sequence types (Table 2). All strains were resistant to the antibiotics tested, according to CLSI standards.

Table 2
www.frontiersin.org

Table 2. Characterization of sensitivity of K. pneumoniae strains to antimicrobial agents.

Three phages were tested: two capsule-specific phages (Dlv622 and Seu621) and one broad-host-range phage FRZ284 that was not capsule-dependent (Gorodnichev et al., 2021a; b). Phages Dlv622 and Seu621 exhibited similar efficacy on strain B536-17-2 as on their original host strain (EOP = 1), but showed 100-fold lower efficacy on strain 9faiz (EOP = 0.01). On strain G4-17, the EOP values for Dlv622 and Seu621 were 0.1 and 0.001, respectively. Phage FRZ284 demonstrated lytic activity across all three strains. However, when titrated by Appelman’s method, phages did not completely inhibit bacterial growth even at the highest concentration (109 PFU/mL).

The recombinant depolymerase Dep622, obtained in a previous study (Gorodnichev et al., 2021a), formed halos on all three K. pneumoniae strains. The MHF values of Dep622 varied across strains, ranging from 0.68 to 1.37 μg/mL.

3.2 Individual and combined effects of antimicrobial agents on planktonic bacterial cells

Sublethal concentrations of phages (109 PFU/well) and maximum concentration of antibiotic in serum (Table 2) were used to study the effects of antimicrobial agents. Depolymerase Dep622 was added at a concentration equivalent to 100 MHF for each strain.

The individual application of antibiotics had minimal antimicrobial effect, except for CMP on strain 9faiz (Figure 1). Phages exhibited a more substantial effect: phages Dlv622 and Seu621 reduced endpoint OD600 and AUC values by 13–55% and 38–64%, respectively, while phage FRZ284 achieved reductions of up to 90%. Depolymerase treatment reduced endpoint OD600 and AUC by no more than 25%.

Figure 1
www.frontiersin.org

Figure 1. Individual effects (endpoint and AUC) of antibiotics, bacteriophages, and depolymerase on planktonic cells of K. pneumoniae. The color gradient shows the degree of reduction. GEN, Gentamicin; LVX, Levofloxacin; MEM, Meropenem; CMP, Chloramphenicol.

Combined effects were assessed using t-test, Bliss, and two-way ANOVA methods, with ANOVA applied to both endpoint and AUC data (Figure 2; Supplementary Figure S1). The t-test identified 25 significant cases of phage-antibiotic effects: 13 cases of synergy, 7 cases of tentative synergy, and 5 cases of antagonism. The Bliss method detected 23 synergistic and 6 antagonistic interactions. Two-way ANOVA identified the fewest significant effects: 9 and 8 cases for endpoint OD600 and AUC, respectively, with 5 cases of synergy based on endpoint OD600 and 3 based on AUC.

Figure 2
www.frontiersin.org

Figure 2. Cases of synergy, tentative synergy, additive effect and antagonism predicted by different methods. p < 0.05.

3.3 Comparison of methods

The t-test and Bliss methods demonstrated the highest convergence and number of statistically significant results. Even when statistical significance was not reached by one method, the effect direction (S-value magnitude or OD value for the agent combination relative to individual effects) remained consistent across approaches (Figure 2). The term “tentative synergy,” introduced for the t-test, showed a strong correlation with Bliss results: of the seven cases of tentative synergy, five aligned with Bliss-defined synergy, one displayed an additive effect, and one showed antagonism (9faiz + FRZ284 + CMP).

Analysis of endpoint OD600 data by two-way ANOVA revealed nine effects consistent with those identified by the t-test and Bliss. Using AUC data, two-way ANOVA identified eight effects: five combinations (two cases of synergy and three of antagonism) aligned with results based on endpoint OD600, while the remaining three cases (one synergy and two antagonism) indicated a more complex nature of phage-antibiotic interactions. In the case of synergy (G4-17 + Seu621 + MEM), the result was consistent with Bliss but showed no significant difference from the most potent agent alone by t-test. In the first antagonism case (9faiz + Dlv622 + MEM), only two-way ANOVA based on AUC values showed a significant effect, whereas in the second (9faiz + FRZ284 + GEN), two-way ANOVA indicated antagonism, contrasting with Bliss-defined synergy and tentative synergy identified by the t-test.

Discrepancies in effects identified by different methods (9faiz + FRZ284 + CMP and 9faiz + FRZ284 + GEN) may be explained by the fact that in these cases both antimicrobial agents reduced optical density and the combined effect of their individual actions exceeded 100%. Consequently, synergy in these combinations would only be observed if the combined optical density were negative—a technically unachievable condition. Thus, accurately determining the effect of these combinations is challenging at the phage and antibiotic concentrations used.

3.4 Checkerboard assay

A checkerboard assay was conducted for three phage-antibiotic combinations against K. pneumoniae strains. For the combination G4-17 + Dlv622 + LVX, synergy was confirmed by all three analytical methods, while for G4-17 + Seu621 + LVX and B536-17-2 + Dlv622 + LVX, synergy was confirmed by both the t-test and Bliss analysis (Figure 2).

The FICI values for G4-17 + Seu621 + LVX and G4-17 + Dlv622 + LVX were 0.5, indicating an additive effect. In contrast, the B536-17-2 + Dlv622 + LVX combination displayed a synergistic effect with a FICI of 0.25 (Figure 3A).

Figure 3
www.frontiersin.org

Figure 3. (A) Heat plots OD600 of checkerboard assays (the color key shows the optical density of the culture at the time of 20 h incubation), (B) checkerboard assays analyzed by the Bliss algorithm. Synergy is shown in red, antagonism in green.

Analysis of these data with Bliss scoring in SynergyFinder v.3.0 showed that the G4-17 + Seu621 + LVX combination exhibited strong synergy at most phage and antibiotic concentrations (Figure 3B). Conversely, G4-17 + Dlv622 + LVX showed weak synergy, with additive and antagonistic effects predominating; however, at concentrations used in previous experiments, Bliss analysis still indicated synergy. For the final combination, B536-17-2 + Dlv622 + LVX, synergy was primarily observed, though antagonism emerged at a phage concentration of 108 PFU/well and antibiotic levels of 4–8 μg/mL, comparable to previous experimental conditions.

4 Discussion

In recent years, the scientific community has actively explored the potential effects of combined phage-antibiotic therapy. Numerous in vitro studies have demonstrated both synergistic and antagonistic interactions (Schindler, 1964; Tagliaferri et al., 2019; Jiang et al., 2020; Liu et al., 2020; Nicholls et al., 2023; Shein et al., 2024), and these effects have been further validated in in vivo models (Oechslin et al., 2017; Grygorcewicz et al., 2020; Santamaría-Corral et al., 2023). These findings highlight the potential for translating such combined approaches into clinical practice, offering promising avenues for improved therapeutic outcomes.

One of the key challenges in studying phage-antibiotic combinations is the difficulty in comparing results across studies due to significant differences in experimental methods, limited characterization of bacterial strains and phages, and a restricted number of tested combinations. These limitations often result in a tendency to describe synergistic effects as unique to specific phage-antibiotic pairs (Bedi et al., 2009; Xu et al., 2022; Zhao et al., 2023). In our study, we aimed to carefully standardize and structure each phase of evaluating the combined effects of phages and antibiotics on K. pneumoniae.

In selecting K. pneumoniae strains, we focused on isolates that were likely to require combination therapy for eradication. We included multidrug-resistant strains with reduced susceptibility to phages, specifically with capsular type KL23 and clinically significant sequence types ST39 and ST1869 (Table 2). The KL23 capsular type is found in approximately 4% of the general population but accounts for 9–17% of carbapenemase-producing strains (Satlin et al., 2017; Sonda et al., 2018). ST39 is associated with several outbreaks of carbapenem-resistant K. pneumoniae in Russia and Greece (Kuzina et al., 2023; Tryfinopoulou et al., 2023). ST1869, although rare in clinical practice, exhibits carbapenem resistance (Li et al., 2019) and differs by only one allele from ST11, a prevalent sequence type in China known for high rates of multidrug resistance and hypervirulence (Gu et al., 2018).

To maximize the diversity of observed effects, we selected phages from different taxonomic families: Dlv622 (Autographiviridae, Slopekvirinae, Drulisvirus), Seu621 (Vequintavirinae, Mydovirus), and FRZ284 (Straboviridae, Tevenvirinae, Jiaodavirus). The lytic phage FRZ284, capable only of lysis from without, was included to demonstrate that even phages unlikely to be chosen for therapy based on spot tests can exhibit synergy with antibiotics. This phage lysed K. pneumoniae strains regardless of capsular type and has a receptor-binding protein of unknown specificity (Gorodnichev et al., 2021b).

The phages Dlv622 and Seu621 possess homologous receptor-binding proteins (query coverage 98%, identity 65.94%), characterized as polysaccharide depolymerases (Gorodnichev et al., 2021a). One of these proteins, Dep622, the receptor-binding protein of phage Dlv622, was used in the current study. Since depolymerase itself does not lyse bacterial cells, investigating patterns in the combined effects of depolymerases and antibiotics offers a novel perspective on the therapeutic application of depolymerases.

We also selected antibiotics from different classes (aminoglycosides, fluoroquinolones, carbapenems, and phenicols) to maximize the diversity of observed effects. These antibiotics are widely used in clinical practice: meropenem is generally considered a first-line treatment for serious K. pneumoniae infections, while gentamicin is frequently used as part of combination therapy for severe infections. The use of chloramphenicol and levofloxacin is typically guided by susceptibility testing results to ensure their effectiveness against resistant strains (CLSI, 2021).

The concentrations of antimicrobial agents used in phage-antibiotic combination studies are pivotal for interpreting results and assessing clinical relevance. Most studies on phage-antibiotic interactions employ sublethal doses of antimicrobials, which may be above or below therapeutic doses depending on the bacterial strain’s susceptibility. In this study, we based antibiotic concentrations on Cmax—the maximum serum concentration of an antibiotic reached after a single dose (Table 1). We believe this approach enhances clinical relevance, providing a more clinically relevant evaluation of combination efficacy. While achieving an equivalent concentration for phages would be ideal, this is likely unfeasible due to phage replication at infection sites.

To assess and compare combined antimicrobial effects, we used several established methods: the t-test, Bliss independence model, two-way ANOVA, and checkerboard assay (Liu et al., 2020; Eller et al., 2021; Gordillo Altamirano et al., 2022; Nikolic et al., 2022; Xu et al., 2022). Among these, the t-test, commonly applied in many studies, only compares the combined effect of agents with the effect of the most effective agent alone. This approach does not account for the contribution of the second agent, thus limiting its ability to determine the overall interaction effect. To address this limitation, we refined the criteria for t-test interpretation, defining conditions for synergy, antagonism, and additive effects. In specific cases where the individual effects of both agents were statistically significantly different from the control and the inhibitory effect of their combination was significantly stronger, we introduced the term “tentative synergy.” This term reflects a potential positive interaction, yet the exact contribution of the second agent remains uncertain, precluding a definitive classification as synergy.

Another frequently used method is the Bliss independence model, which assumes that two agents act independently through stochastic processes (Eller et al., 2021). Given that antibiotics and bacteriophages target bacteria via distinct biochemical mechanisms, the Bliss model is applicable for analyzing their combined effects. However, a notable drawback of this method is its high sensitivity—it does not account for the absence of an effect. Consequently, if the effects of both agents and their combination do not significantly differ from the control, the Bliss model will classify the interaction as additive, even when a true interaction might be absent.

Another commonly used approach for assessing interactions is two-way ANOVA, which specifically evaluates the interaction term, making it relevant for distinguishing between synergy and antagonism. However, the interaction p-value alone does not provide information on the direction of the effect. This interpretation can be inferred from the slope of interaction plots (Liu et al., 2020), though standard error may limit the accuracy of this approach.

Our results indicate that the t-test and Bliss models identified the highest number of statistically significant effects, likely due to the lack of adjustments for multiple comparisons, allowing for the detection of less pronounced synergy and antagonism effects (Figure 2). In contrast, the more conservative two-way ANOVA, which requires a strong and consistent effect to reach statistical significance, identified considerably fewer effects.

The broad range of results highlights the need to select methods based on specific research objectives. For selecting potential phage-antibiotic combinations for therapeutic application, using the t-test or Bliss (or a combination of these) may be optimal due to their simplicity in calculation and interpretation. Additionally, for therapeutic purposes, identifying potential antagonistic interactions can be more critical than seeking combinations with the highest synergy. In turn, two-way ANOVA could be an effective tool for exploratory research, particularly when the objective is to pinpoint combinations with strong interactions for deeper investigation into bacterial response.

The checkerboard assay, with its evaluation of interactions across a broad concentration range, has unique advantages. While this approach appears more informative, the FICI values calculated for two out of three cases conflicted with the effects observed in t-test, Bliss, and two-way ANOVA analyses. Notably, the checkerboard assay operates at antimicrobial agent concentrations unlikely to be achieved in real therapeutic scenarios. However, checkerboard results, when supplemented with Bliss analysis, can illustrate how the combined effect varies with agent concentrations. For example, when using linezolid near Cmax (4–8 μg/mL), the combinations G4-17 + Seu621 + LVX and B536-2-17 + Dlv622 + LVX consistently demonstrate synergy at any phage concentration, while a similar effect for G4-17 + Dlv622 + LVX is observed only at the highest phage concentrations. Such observations may be valuable for understanding the potential variability in phage titers at infection sites (Schooley et al., 2017).

Overall, while these methods cannot be directly compared due to differences in statistical assumptions and sensitivity, their complementary use provides a more robust assessment of phage-antibiotic interactions. The consistency of results across multiple methods strengthens confidence in observed effects, whereas discrepancies highlight potential context-dependent factors such as dose–response relationships. These findings emphasize the importance of selecting the appropriate analytical framework based on the specific goals of a study rather than relying on a single method.

In total, 35 of the 48 phage-antibiotic combinations tested demonstrated statistically significant effects based on at least one method (Figure 4). Two previously described cases (9faiz + FRZ284 + CMP and 9faiz + FRZ284 + GEN) were excluded due to difficulties in accurately interpreting their effects; these combinations were thus classified as statistically insignificant.

Figure 4
www.frontiersin.org

Figure 4. Cumulative effects of combinations of antibiotics and phages or depolymerase on all three strains from a practical point of view. Synergy is shown in blue, antagonism in orange, additive effect in gray.

Phage and depolymerase combinations with levofloxacin exhibited the highest proportion of synergistic effects. However, only the combinations of levofloxacin with Seu621 phage or Dep622 depolymerase consistently produced effects across all three K. pneumoniae strains tested. Previous studies have shown synergy between fluoroquinolone antibiotics and phages in disrupting K. pneumoniae biofilms, specifically with a ciprofloxacin and T7-like phage combination (Verma et al., 2009); however, this is the first report of synergy involving a Mydovirus phage or depolymerase. Nevertheless, expanding the diversity of bacterial strains tested may impact the consistent synergistic interaction observed between levofloxacin and either Seu621 phage or Dep622 depolymerase. In fact, as shown in the checkerboard assay, certain effects may reverse depending on the antimicrobial concentration (Figure 3).

Gentamicin ranked second in terms of the proportion of synergistic effects observed. Although synergistic effects between aminoglycoside antibiotics and K. pneumoniae phages have been reported in animal models (Wang et al., 2021; Zhao et al., 2023), numerous studies suggest that aminoglycosides may inhibit phages (Schindler, 1964; Bowman, 1967; Jiang et al., 2020; Kever et al., 2022). In our study, all phage and depolymerase combinations with gentamicin showed trends toward synergy or additive effects, with antagonism detected in only one case.

Phage combinations with meropenem and chloramphenicol yielded the highest proportion of statistically insignificant results and a higher frequency of antagonistic effects. Beta-lactam antibiotics, particularly meropenem, are frequently explored in studies examining phage synergy. For example, meropenem combined with a K. pneumoniae Webervirus phage showed synergy in biofilm degradation (Xu et al., 2022), while amoxicillin reduced biofilm biomass when combined with an unclassified phage (Bedi et al., 2009). Synergy between K. pneumoniae phages and chloramphenicol has not been previously documented; however, similar effects have been observed with a P. aeruginosa N4-like phage KPP21 (Uchiyama et al., 2018) and Inovirus phages (Hagens et al., 2006).

Notably, combinations involving the FRZ284 phage and the Dep622 depolymerase revealed diverse effects. In plating assays, phage FRZ284 only displayed lysis from without, indicating that while it could attach to bacteria, it could not proliferate within the strain. Despite this, cases of synergy were observed with three out of four antibiotics (gentamicin, levofloxacin, and meropenem), alongside antagonistic effects with all four antibiotics. The Dep622 depolymerase itself did not directly kill bacteria but sensitized them to antimicrobial agents. In this context, synergistic effects with antibiotics were expected, as demonstrated with gentamicin, levofloxacin, and meropenem. However, the antagonistic effect observed with chloramphenicol on strain 9faiz—validated across all testing methods—remains challenging to interpret.

Combinations of phages Dlv622 or FRZ284 with meropenem, as well as FRZ284 with gentamicin, demonstrated all three types of interactions (synergy, additive effect, and antagonism) across the three K. pneumoniae strains. This unpredictability stems not only from methodological inconsistencies but also from the genetic diversity of bacterial strains and the variability of phage lifecycles, which significantly influence outcomes.

The numerous studies have described the mechanisms underlying synergy and antagonism in phage-antibiotic interactions (Liu et al., 2020; Zuo et al., 2021; Bulssico et al., 2023; Luo et al., 2023; Pons et al., 2023; Kunz Coyne et al., 2024; Qin et al., 2024; Yu et al., 2025). However, it is important to recognize that these mechanisms have been described for specific bacteria-phage-antibiotic combinations, and given the vast diversity of bacterial strains and phages, they may not represent the only possible explanations. The variability in effects observed in our study further supports the idea that multiple, context-dependent factors influence phage-antibiotic interactions. Given this complexity, predicting the clinical applicability of specific phage-antibiotic combinations requires careful experimental validation, but their translation into clinical settings remains uncertain without standardized experimental models.

To advance the field, standardizing experimental approaches for assessing combination effects is crucial. The goal of our study is to initiate this discussion by demonstrating how methodological variability can lead to inconsistent interpretations.

Ultimately, the development of effective phage-antibiotic therapies requires a shift from studying isolated mechanisms to a systematic, standardized framework that accounts for the complexity of microbial interactions. Establishing such approaches will improve reproducibility and enable the rational design of phage-based treatments in clinical practice.

5 Conclusion

The combined application of phages and antibiotics in therapy requires reliable methods to assess their interactions; however, the lack of standardized approaches complicates data interpretation and clinical implementation. In this study, we applied the most commonly used methods (t-test, Bliss method, two-way ANOVA, and checkerboard assay), complemented by the standardization of antibiotic concentrations at Cmax. The t-test proved valuable for rapid interaction assessment, especially after expanding the criteria and introducing “tentative synergy,” allowing consideration of complex cases where the contribution of the second agent is less apparent. The Bliss method aligned well with the t-test results and was particularly effective when combined with checkerboard assay data, enabling visualization of the interaction effect based on agent concentration. In contrast, the FICI was less informative and poorly correlated with other methods. Two-way ANOVA identified only the most consistent and substantial effects, highlighting strong interactions and making it useful for comprehensive investigations.

Overall, our approach identified only two combinations with consistent effects across all three K. pneumoniae strains, suggesting that synergy or antagonism is not solely determined by the nature of phage or antibiotic but depends on a complex interplay of factors. These findings emphasize the need for standardized interaction assessment methods to improve comparability of data and enhance the efficacy of combination therapies.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

RG: Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing. AK: Formal analysis, Investigation, Methodology, Writing – original draft. MK: Investigation, Methodology, Writing – original draft. NA: Investigation, Methodology, Software, Writing – original draft. MM: Investigation, Writing – original draft. MZ: Investigation, Writing – original draft. DB: Formal analysis, Visualization, Writing – original draft. VM: Investigation, Writing – original draft. ES: Conceptualization, Supervision, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was supported by the State Assignment “Bacteriophage-2” R&D state registration number 122022800139-0.

Acknowledgments

Part of this work was performed using the core facilities of the Lopukhin FRCC PCM “Genomics, proteomics, metabolomics” (http://rcpcm.org/?p=2806).

Conflict of interest

The authors declare that the research 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 authors declare that no Gen AI was used in the creation of this manuscript.

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/fmicb.2025.1530819/full#supplementary-material

References

Bedi, M. S., Verma, V., and Chhibber, S. (2009). Amoxicillin and specific bacteriophage can be used together for eradication of biofilm of Klebsiella pneumoniae B5055. World J. Microbiol. Biotechnol. 25, 1145–1151. doi: 10.1007/s11274-009-9991-8

Crossref Full Text | Google Scholar

Bowman, B. U. (1967). Biological activity of φX DNA: I. Inhibition of infectivity by streptomycin. J. Mol. Biol. 25, 559–561. doi: 10.1016/0022-2836(67)90207-0

PubMed Abstract | Crossref Full Text | Google Scholar

Brisse, S., Passet, V., Haugaard, A. B., Babosan, A., Kassis-Chikhani, N., Struve, C., et al. (2013). Wzi gene sequencing, a rapid method for determination of capsulartype for klebsiella strains. J. Clin. Microbiol. 51, 4073–4078. doi: 10.1128/JCM.01924-13

PubMed Abstract | Crossref Full Text | Google Scholar

Bulssico, J., Papukashvili, I., Espinosa, L., Gandon, S., and Ansaldi, M. (2023). Phage-antibiotic synergy: cell filamentation is a key driver of successful phage predation. PLoS Pathog. 19:e1011602. doi: 10.1371/journal.ppat.1011602

PubMed Abstract | Crossref Full Text | Google Scholar

Burrowes, B. H., Molineux, I. J., and Fralick, J. A. (2019). Directed in vitro evolution of therapeutic bacteriophages: the Appelmans protocol. Viruses 11:241. doi: 10.3390/v11030241

PubMed Abstract | Crossref Full Text | Google Scholar

Chang, D., Sharma, L., Dela Cruz, C. S., and Zhang, D. (2021). Clinical epidemiology, risk factors, and control strategies of Klebsiella pneumoniae infection. Front. Microbiol. 12:750662. doi: 10.3389/FMICB.2021.750662

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, J., Li, J., Huang, F., Fang, J., Cao, Y., Zhang, K., et al. (2023). Clinical characteristics, risk factors and outcomes of Klebsiella pneumoniae pneumonia developing secondary Klebsiella pneumoniae bloodstream infection. BMC Pulm. Med. 23:102. doi: 10.1186/s12890-023-02394-8

PubMed Abstract | Crossref Full Text | Google Scholar

Chhibber, S., Kaur, S., and Kumari, S. (2008). Therapeutic potential of bacteriophage in treating Klebsiella pneumoniae B5055-mediated lobar pneumonia in mice. J. Med. Microbiol. 57, 1508–1513. doi: 10.1099/jmm.0.2008/002873-0

PubMed Abstract | Crossref Full Text | Google Scholar

CLSI (2015). M07-A10: Methods for dilution antimicrobial susceptibility tests for bacteria that grow aerobically. Clinical and Laboratory Standards Institute. Available online at: https://www.clsi.org.%0Awww.clsi.org.

Google Scholar

CLSI (2021). Performance standards for antimicrobial susceptibility testing, M100, 31st ed. Wayne, PA: Clinical and Laboratory Standards Institute.

Google Scholar

Danis-Wlodarczyk, K. M., Wozniak, D. J., and Abedon, S. T. (2021). Treating bacterial infections with bacteriophage-based enzybiotics: in vitro, in vivo and clinical application. Antibiotics 10, 1–36. doi: 10.3390/antibiotics10121497

PubMed Abstract | Crossref Full Text | Google Scholar

Diancourt, L., Passet, V., Verhoef, J., Grimont, P. A. D., and Brisse, S. (2005). Multilocus sequence typing of Klebsiella pneumoniae nosocomial isolates. J. Clin. Microbiol. 43, 4178–4182. doi: 10.1128/JCM.43.8.4178-4182.2005

PubMed Abstract | Crossref Full Text | Google Scholar

Durbas, I., and Machnik, G. (2022). Phage therapy: an old concept with new perspectives. J. Appl. Pharmaceut. Sci. 12, 27–38. doi: 10.7324/JAPS.2022.120502

Crossref Full Text | Google Scholar

Eller, K. A., Aunins, T. R., Courtney, C. M., Campos, J. K., Otoupal, P. B., Erickson, K. E., et al. (2021). Facile accelerated specific therapeutic (FAST) platform develops antisense therapies to counter multidrug-resistant bacteria. Commun. Biol. 4:331. doi: 10.1038/s42003-021-01856-1

PubMed Abstract | Crossref Full Text | Google Scholar

European Centre for Disease Prevention and Control (2022). Antimicrobial resistance in the EU/EEA (EARS-net) AER for 2021. Stockholm. Available online at: https://atlas.ecdc.europa.eu/ (Accessed December 5, 2022).

Google Scholar

Gonçalves-Pereira, J., Martins, A., and Póvoa, P. (2010). Pharmacokinetics of gentamicin in critically ill patients: pilot study evaluating the first dose. Clin. Microbiol. Infect. 16, 1258–1263. doi: 10.1111/j.1469-0691.2009.03074.x

PubMed Abstract | Crossref Full Text | Google Scholar

Gordillo Altamirano, F. L., Kostoulias, X., Subedi, D., Korneev, D., Peleg, A. Y., and Barr, J. J. (2022). Phage-antibiotic combination is a superior treatment against Acinetobacter baumannii in a preclinical study. EBioMedicine 80:104045. doi: 10.1016/j.ebiom.2022.104045

PubMed Abstract | Crossref Full Text | Google Scholar

Gorodnichev, R., Kornienko, M., Kuptsov, N., Malakhova, M., Bespiatykh, D., Veselovsky, V., et al. (2021b). Molecular genetic characterization of three new Klebsiella pneumoniae bacteriophages suitable for phage therapy. Med. Extrem. Situat. 3, 90–97. doi: 10.47183/mes.2021.035

Crossref Full Text | Google Scholar

Gorodnichev, R. B., Volozhantsev, N. V., Krasilnikova, V. M., Bodoev, I. N., Kornienko, M. A., Kuptsov, N. S., et al. (2021a). Novel Klebsiella pneumoniae K23-specific bacteriophages from different families: similarity of Depolymerases and their therapeutic potential. Front. Microbiol. 12:669618. doi: 10.3389/FMICB.2021.669618

PubMed Abstract | Crossref Full Text | Google Scholar

Grosjean, V., Gressens, S. B., Pham, T., Gaudry, S., Ait-Oufella, H., De Prost, N., et al. (2024). Community-acquired Klebsiella pneumoniae pneumonia in ICU: a multicenter retrospective study. Ann. Intensive Care 14:69. doi: 10.1186/s13613-024-01269-3

PubMed Abstract | Crossref Full Text | Google Scholar

Grygorcewicz, B., Roszak, M., Golec, P., Śleboda-taront, D., Łubowska, N., Górska, M., et al. (2020). Antibiotics act with vb_abap_agc01 phage against acinetobacter baumannii in human heat-inactivated plasma blood and galleria mellonella models. Int. J. Mol. Sci. 21, 1–14. doi: 10.3390/ijms21124390

PubMed Abstract | Crossref Full Text | Google Scholar

Gu, D., Dong, N., Zheng, Z., Lin, D., Huang, M., Wang, L., et al. (2018). A fatal outbreak of ST11 carbapenem-resistant hypervirulent Klebsiella pneumoniae in a Chinese hospital: a molecular epidemiological study. Lancet Infect. Dis. 18, 37–46. doi: 10.1016/S1473-3099(17)30489-9

PubMed Abstract | Crossref Full Text | Google Scholar

Gu, Z., Eils, R., and Schlesner, M. (2016). Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 32, 2847–2849. doi: 10.1093/bioinformatics/btw313

PubMed Abstract | Crossref Full Text | Google Scholar

Hagens, S., Habel, A., and Bläsi, U. (2006). Augmentation of the antimicrobial efficacy of antibiotics by filamentous phage. Microb. Drug Resist. 12, 164–168. doi: 10.1089/mdr.2006.12.164

PubMed Abstract | Crossref Full Text | Google Scholar

Han, M.-L., Nang, S. C., Lin, Y.-W., Zhu, Y., Yu, H. H., Wickremasinghe, H., et al. (2022). Comparative metabolomics revealed key pathways associated with the synergistic killing of multidrug-resistant Klebsiella pneumoniae by a bacteriophage-polymyxin combination. Comput. Struct. Biotechnol. J. 20, 485–495. doi: 10.1016/j.csbj.2021.12.039

PubMed Abstract | Crossref Full Text | Google Scholar

Hua, Y., Wu, Y., Guo, M., Ma, R., Li, Q., Hu, Z., et al. (2022). Characterization and functional studies of a novel Depolymerase against K19-type Klebsiella pneumoniae. Front. Microbiol. 13:878800. doi: 10.3389/fmicb.2022.878800

PubMed Abstract | Crossref Full Text | Google Scholar

Ianevski, A., Giri, A. K., and Aittokallio, T. (2022). SynergyFinder 3.0: an interactive analysis and consensus interpretation of multi-drug synergies across multiple samples. Nucleic Acids Res. 50, W739–W743. doi: 10.1093/nar/gkac382

PubMed Abstract | Crossref Full Text | Google Scholar

Jiang, Z., Wei, J., Liang, Y., Peng, N., and Li, Y. (2020). Aminoglycoside antibiotics inhibit Mycobacteriophage infection. Antibiotics 9, 1–6. doi: 10.3390/antibiotics9100714

PubMed Abstract | Crossref Full Text | Google Scholar

Kever, L., Hardy, A., Luthe, T., Hünnefeld, M., Gätgens, C., Milke, L., et al. (2022). Aminoglycoside antibiotics inhibit phage infection by blocking an early step of the infection cycle. mBio 13:e0078322. doi: 10.1128/mbio.00783-22

PubMed Abstract | Crossref Full Text | Google Scholar

Kunz Coyne, A. J., Eshaya, M., Bleick, C., Vader, S., Biswas, B., Wilson, M., et al. (2024). Exploring synergistic and antagonistic interactions in phage-antibiotic combinations against ESKAPE pathogens. Microbiol. Spectr. 12:e0042724. doi: 10.1128/spectrum.00427-24

PubMed Abstract | Crossref Full Text | Google Scholar

Kutter, E. (2009). Phage host range and efficiency of plating. Methods Mol. Biol. (Clifton, N.J.) 501, 141–149. doi: 10.1007/978-1-60327-164-6_14/FIGURES/14_2_978-1-60327-164-6

Crossref Full Text | Google Scholar

Kuzina, E. S., Kislichkina, A. A., Sizova, A. A., Skryabin, Y. P., Novikova, T. S., Ershova, O. N., et al. (2023). High-molecular-weight plasmids carrying Carbapenemase genes blaNDM-1, blaKPC-2, and blaOXA-48 coexisting in clinical Klebsiella pneumoniae strains of ST39. Microorganisms 11:459. doi: 10.3390/microorganisms11020459

PubMed Abstract | Crossref Full Text | Google Scholar

Lew, B. Y. X., Njondimackal, N. L., Ravisankar, V., and Norman, N. A. (2023). “Enhanced antibacterial activity of a novel phage-antibiotic combination against Klebsiella pneumoniae isolates” in Proceedings of the 9th IRC conference on science, engineering, and technology. IRC-SET 2023 (Singapore: Springer), 244–256.

Google Scholar

Li, P., Guo, G., Zheng, X., Xu, S., Zhou, Y., Qin, X., et al. (2024). Therapeutic efficacy of a K5-specific phage and depolymerase against Klebsiella pneumoniae in a mouse model of infection. Vet. Res. 55:59. doi: 10.1186/s13567-024-01311-z

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Y., Shen, H., Zhu, C., and Yu, Y. (2019). Carbapenem-resistant Klebsiella pneumoniae infections among ICU admission patients in Central China: prevalence and prediction model. Biomed. Res. Int. 2019:9767313. doi: 10.1155/2019/9767313

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, C. G., Green, S. I., Min, L., Clark, J. R., Salazar, K. C., Terwilliger, A. L., et al. (2020). Phage-antibiotic synergy is driven by a unique combination of antibacterial mechanism of action and stoichiometry. MBio 11, 1–19. doi: 10.1128/mBio.01462-20

PubMed Abstract | Crossref Full Text | Google Scholar

Luo, J., Xie, L., Yang, M., Liu, M., Li, Q., Wang, P., et al. (2023). Synergistic antibacterial effect of phage pB3074 in combination with antibiotics targeting Cell Wall against multidrug-resistant Acinetobacter baumannii in vitro and ex vivo. Microbiol. Spectr. 11:e0034123. doi: 10.1128/spectrum.00341-23

PubMed Abstract | Crossref Full Text | Google Scholar

Mazzocco, A., Waddell, T. E., Lingohr, E., and Johnson, R. P. (2009). Enumeration of bacteriophages using the small drop plaque assay system. Bacteriophages 501, 81–85. doi: 10.1007/978-1-60327-164-6

Crossref Full Text | Google Scholar

Murray, C. J., Ikuta, K. S., Sharara, F., Swetschinski, L., Robles Aguilar, G., Gray, A., et al. (2022). Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet 399, 629–655. doi: 10.1016/S0140-6736(21)02724-0

PubMed Abstract | Crossref Full Text | Google Scholar

Nicholls, P., Clark, J. R., Gu Liu, C., Terwilliger, A., and Maresso, A. W. (2023). Class-driven synergy and antagonism between a Pseudomonas phage and antibiotics. Infect. Immun. 91:e0006523. doi: 10.1128/iai.00065-23

PubMed Abstract | Crossref Full Text | Google Scholar

Nikolic, I., Vukovic, D., Gavric, D., Cvetanovic, J., Aleksic Sabo, V., Gostimirovic, S., et al. (2022). An optimized checkerboard method for phage-antibiotic synergy detection. Viruses 14:1542. doi: 10.3390/v14071542

PubMed Abstract | Crossref Full Text | Google Scholar

Oechslin, F., Piccardi, P., Mancini, S., Gabard, J., Moreillon, P., Entenza, J. M., et al. (2017). Synergistic interaction between phage therapy and antibiotics clears Pseudomonas Aeruginosa infection in endocarditis and reduces virulence. J. Infect. Dis. 215, 703–712. doi: 10.1093/infdis/jiw632

PubMed Abstract | Crossref Full Text | Google Scholar

Onallah, H., Hazan, R., and Nir-Paz, R.Israeli Phage Therapy Center (IPTC) Study Team (2023). Compassionate use of bacteriophages for failed persistent infections during the first 5 years of the Israeli phage therapy center. Open Forum Infect. Dis. 10:ofad221. doi: 10.1093/ofid/ofad221

PubMed Abstract | Crossref Full Text | Google Scholar

Pacios, O., Fernández-García, L., Bleriot, I., Blasco, L., González-Bardanca, M., López, M., et al. (2021). Enhanced antibacterial activity of repurposed mitomycin C and imipenem in combination with the lytic phage vB_KpnMVAC13 against clinical isolates of klebsiella pneumoniae. Antimicrob. Agents Chemother. 65:e0090021. doi: 10.1128/AAC.00900-21

PubMed Abstract | Crossref Full Text | Google Scholar

Pirnay, J.-P., Djebara, S., Steurs, G., Griselain, J., Cochez, C., De Soir, S., et al. (2024). Personalized bacteriophage therapy outcomes for 100 consecutive cases: a multicentre, multinational, retrospective observational study. Nat. Microbiol. 9, 1434–1453. doi: 10.1038/s41564-024-01705-x

PubMed Abstract | Crossref Full Text | Google Scholar

Pons, B. J., Dimitriu, T., Westra, E. R., and van Houte, S. (2023). Antibiotics that affect translation can antagonize phage infectivity by interfering with the deployment of counter-defenses. Proc. Natl. Acad. Sci. USA 120:e2216084120. doi: 10.1073/pnas.2216084120

PubMed Abstract | Crossref Full Text | Google Scholar

Qin, K., Shi, X., Yang, K., Xu, Q., Wang, F., Chen, S., et al. (2024). Phage-antibiotic synergy suppresses resistance emergence of Klebsiella pneumoniae by altering the evolutionary fitness. MBio 15:e0139324. doi: 10.1128/mbio.01393-24

PubMed Abstract | Crossref Full Text | Google Scholar

Qurat-ul-Ain, H., Ijaz, M., Siddique, A. B., Muzammil, S., Shafique, M., Rasool, M. H., et al. (2021). Efficacy of phage-antibiotic combinations against multidrug-resistant Klebsiella pneumoniae clinical isolates. Jundishapur J. Microbiol. 14:e111926. doi: 10.5812/jjm.111926

Crossref Full Text | Google Scholar

Rahmat Ullah, S., Jamal, M., Rahman, A., and Andleeb, S. (2024). Comprehensive insights into Klebsiella pneumoniae: unravelling clinical impact, epidemiological trends and antibiotic-resistance challenges. J. Antimicrob. Chemother. 79, 1484–1492. doi: 10.1093/jac/dkae184

PubMed Abstract | Crossref Full Text | Google Scholar

Santamaría-Corral, G., Senhaji-Kacha, A., Broncano-Lavado, A., Esteban, J., and García-Quintanilla, M. (2023). Bacteriophage-antibiotic combination therapy against Pseudomonas aeruginosa. Antibiotics (Basel) 12:1089. doi: 10.3390/antibiotics12071089

PubMed Abstract | Crossref Full Text | Google Scholar

Satlin, M. J., Chen, L., Patel, G., Gomez-Simmonds, A., Weston, G., Kim, A. C., et al. (2017). Multicenter clinical and molecular epidemiological analysis of bacteremia due to Carbapenem-resistant Enterobacteriaceae (CRE) in the CRE epicenter of the United States. Antimicrob. Agents Chemother. 61, 1–13. doi: 10.1128/AAC.02349-16

PubMed Abstract | Crossref Full Text | Google Scholar

Schindler, J. (1964). Inhibition of reproduction of the f2 bacteriophage by streptomycin. Folia Microbiol. 9, 269–276. doi: 10.1007/BF02873305

PubMed Abstract | Crossref Full Text | Google Scholar

Schooley, R. T., Biswas, B., Gill, J. J., Hernandez-Morales, A., Lancaster, J., Lessor, L., et al. (2017). Development and use of personalized bacteriophage-based therapeutic cocktails to treat a patient with a disseminated resistant Acinetobacter baumannii infection. Antimicrob. Agents Chemother. 61, 1–15. doi: 10.1128/AAC.00954-17

PubMed Abstract | Crossref Full Text | Google Scholar

Shein, A. M. S., Wannigama, D. L., Hurst, C., Monk, P. N., Amarasiri, M., Wongsurawat, T., et al. (2024). Phage cocktail amikacin combination as a potential therapy for bacteremia associated with carbapenemase producing colistin resistant Klebsiella pneumoniae. Sci. Rep. 14:28992. doi: 10.1038/s41598-024-79924-9

Crossref Full Text | Google Scholar

Sonda, T., Kumburu, H., van Zwetselaar, M., Alifrangis, M., Mmbaga, B. T., Lund, O., et al. (2018). Molecular epidemiology of virulence and antimicrobial resistance determinants in Klebsiella pneumoniae from hospitalised patients in Kilimanjaro, Tanzania. Eur. J. Clin. Microbiol. Infect. Dis. 37, 1901–1914. doi: 10.1007/s10096-018-3324-5

Crossref Full Text | Google Scholar

Tagliaferri, T. L., Jansen, M., and Horz, H.-P. (2019). Fighting pathogenic Bacteria on two fronts: phages and antibiotics as combined strategy. Front. Cell. Infect. Microbiol. 9:22. doi: 10.3389/fcimb.2019.00022

PubMed Abstract | Crossref Full Text | Google Scholar

Thalhammer, F., Schenk, P., Burgmann, H., El Menyawi, I., Hollenstein, U. M., Rosenkranz, A. R., et al. (1998). Single-dose pharmacokinetics of meropenem during continuous venovenous hemofiltration. Antimicrob. Agents Chemother. 42, 2417–2420. doi: 10.1128/AAC.42.9.2417

PubMed Abstract | Crossref Full Text | Google Scholar

Tryfinopoulou, K., Linkevicius, M., Pappa, O., Alm, E., Karadimas, K., Svartström, O., et al. (2023). Emergence and persistent spread of carbapenemase-producing Klebsiella pneumoniae high-risk clones in Greek hospitals, 2013 to 2022. Euro Surveill. 28:2300571. doi: 10.2807/1560-7917.ES.2023.28.47.2300571

PubMed Abstract | Crossref Full Text | Google Scholar

Uchiyama, J., Shigehisa, R., Nasukawa, T., Mizukami, K., Takemura-Uchiyama, I., Ujihara, T., et al. (2018). Piperacillin and ceftazidime produce the strongest synergistic phage-antibiotic effect in Pseudomonas aeruginosa. Arch. Virol. 163, 1941–1948. doi: 10.1007/s00705-018-3811-0

PubMed Abstract | Crossref Full Text | Google Scholar

Verma, V., Harjai, K., and Chhibber, S. (2009). Restricting ciprofloxacin-induced resistant variant formation in biofilm of Klebsiella pneumoniae B5055 by complementary bacteriophage treatment. J. Antimicrob. Chemother. 64, 1212–1218. doi: 10.1093/jac/dkp360

Crossref Full Text | Google Scholar

Verma, V., Harjai, K., and Chhibber, S. (2010). Structural changes induced by a lytic bacteriophage make ciprofloxacin effective against older biofilm of Klebsiella pneumoniae. Biofouling 26, 729–737. doi: 10.1080/08927014.2010.511196

PubMed Abstract | Crossref Full Text | Google Scholar

Volozhantsev, N. V., Shpirt, A. M., Borzilov, A. I., Komisarova, E. V., Krasilnikova, V. M., Shashkov, A. S., et al. (2020). Characterization and therapeutic potential of bacteriophage-encoded polysaccharide depolymerases with β galactosidase activity against klebsiella pneumoniae K57 capsular type. Antibiotics 9, 1–16. doi: 10.3390/antibiotics9110732

PubMed Abstract | Crossref Full Text | Google Scholar

Wagenlehner, F. M. E., Kinzig-Schippers, M., Sörgel, F., Weidner, W., and Naber, K. G. (2006). Concentrations in plasma, urinary excretion and bactericidal activity of levofloxacin (500 mg) versus ciprofloxacin (500 mg) in healthy volunteers receiving a single oral dose. Int. J. Antimicrob. Agents 28, 551–559. doi: 10.1016/j.ijantimicag.2006.07.026

Crossref Full Text | Google Scholar

Wang, Z., Cai, R., Wang, G., Guo, Z., Liu, X., Guan, Y., et al. (2021). Combination therapy of phage vB_KpnM_P-KP2 and gentamicin combats acute pneumonia caused by K47 serotype Klebsiella pneumoniae. Front. Microbiol. 12:674068. doi: 10.3389/fmicb.2021.674068

PubMed Abstract | Crossref Full Text | Google Scholar

Wenk, M., Vozeh, S., and Follath, F. (1984). Serum level monitoring of antibacterial drugs. A review. Clin. Pharmacokinet 9, 475–492. doi: 10.2165/00003088-198409060-00001

Crossref Full Text | Google Scholar

World Health Organization (2024). WHO Bacterial Priority Pathogens List, 2024: bacterial pathogens of public health importance to guide research, development and strategies to prevent and control antimicrobial resistance : World Health Organization.

Google Scholar

Xu, W., Zhao, Y., Qian, C., Yao, Z., Chen, T., Wang, L., et al. (2022). The identification of phage vB_1086 of multidrug-resistant Klebsiella pneumoniae and its synergistic effects with ceftriaxone. Microb. Pathog. 171:105722. doi: 10.1016/j.micpath.2022.105722

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, Y., Wang, M., Ju, L., Li, M., Zhao, M., Deng, H., et al. (2025). Phage-mediated virulence loss and antimicrobial susceptibility in carbapenem-resistant Klebsiella pneumoniae. MBio 16:e0295724. doi: 10.1128/mbio.02957-24

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, Y., Feng, L., Zhou, B., Zhang, X., Yao, Z., Wang, L., et al. (2023). A newly isolated bacteriophage vB8388 and its synergistic effect with aminoglycosides against multi-drug resistant Klebsiella oxytoca strain FK-8388. Microb. Pathog. 174:105906. doi: 10.1016/j.micpath.2022.105906

PubMed Abstract | Crossref Full Text | Google Scholar

Zuo, P., Yu, P., and Alvarez, P. J. J. (2021). Aminoglycosides antagonize bacteriophage proliferation, attenuating phage suppression of bacterial growth, biofilm formation, and antibiotic resistance. Appl. Environ. Microbiol. 87:e0046821. doi: 10.1128/AEM.00468-21

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: phage, antibiotic, capsule depolymerase, synergy, antagonism, Klebsiella pneumoniae

Citation: Gorodnichev RB, Krivulia AO, Kornienko MA, Abdraimova NK, Malakhova MV, Zaychikova MV, Bespiatykh DA, Manuvera VA and Shitikov EA (2025) Phage-antibiotic combinations against Klebsiella pneumoniae: impact of methodological approaches on effect evaluation. Front. Microbiol. 16:1530819. doi: 10.3389/fmicb.2025.1530819

Received: 19 November 2024; Accepted: 17 February 2025;
Published: 12 March 2025.

Edited by:

Tamas Feher, HUN-REN Biological Research Centre, Hungary

Reviewed by:

Swapnil Ganesh Sanmukh, Université Clermont Auvergne, France
Hosoon Choi, United States Department of Veterans Affairs, United States

Copyright © 2025 Gorodnichev, Krivulia, Kornienko, Abdraimova, Malakhova, Zaychikova, Bespiatykh, Manuvera and Shitikov. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Roman B. Gorodnichev, Z29yb2RuaWNoZXYuci5iQGdtYWlsLmNvbQ==

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.

Research integrity at Frontiers

Man ultramarathon runner in the mountains he trains at sunset

94% of researchers rate our articles as excellent or good

Learn more about the work of our research integrity team to safeguard the quality of each article we publish.


Find out more