Metapopulation model of phage therapy of an acute Pseudomonas aeruginosa lung infection.

Infections caused by multidrug resistant (MDR) pathogenic bacteria are a global health threat. Bacteriophages ("phage") are increasingly used as alternative or last-resort therapeutics to treat patients infected by MDR bacteria. However, the therapeutic outcomes of phage therapy may be limited by the emergence of phage resistance during treatment and/or by physical constraints that impede phage-bacteria interactions in vivo. In this work, we evaluate the role of lung spatial structure on the efficacy of phage therapy for Pseudomonas aeruginosa infections. To do so, we developed a spatially structured metapopulation network model based on the geometry of the bronchial tree, including host innate immune responses and the emergence of phage-resistant bacterial mutants. We model the ecological interactions between bacteria, phage, and the host innate immune system at the airway (node) level. The model predicts the synergistic elimination of a P. aeruginosa infection due to the combined effects of phage and neutrophils, given the sufficient innate immune activity and efficient phage-induced lysis. The metapopulation model simulations also predict that MDR bacteria are cleared faster at distal nodes of the bronchial tree. Notably, image analysis of lung tissue time series from wild-type and lymphocyte-depleted mice revealed a concordant, statistically significant pattern: infection intensity cleared in the bottom before the top of the lungs. Overall, the combined use of simulations and image analysis of in vivo experiments further supports the use of phage therapy for treating acute lung infections caused by P. aeruginosa, while highlighting potential limits to therapy in a spatially structured environment given impaired innate immune responses and/or inefficient phage-induced lysis. IMPORTANCE: Phage therapy is increasingly employed as a compassionate treatment for severe infections caused by multidrug-resistant (MDR) bacteria. However, the mixed outcomes observed in larger clinical studies highlight a gap in understanding when phage therapy succeeds or fails. Previous research from our team, using in vivo experiments and single-compartment mathematical models, demonstrated the synergistic clearance of acute P. aeruginosa pneumonia by phage and neutrophils despite the emergence of phage-resistant bacteria. In fact, the lung environment is highly structured, prompting the question of whether immunophage synergy explains the curative treatment of P. aeruginosa when incorporating realistic physical connectivity. To address this, we developed a metapopulation network model mimicking the lung branching structure to assess phage therapy efficacy for MDR P. aeruginosa pneumonia. The model predicts the synergistic elimination of P. aeruginosa by phage and neutrophils but emphasizes potential challenges in spatially structured environments, suggesting that higher innate immune levels may be required for successful bacterial clearance. Model simulations reveal a spatial pattern in pathogen clearance where P. aeruginosa are cleared faster at distal nodes of the bronchial tree than in primary nodes. Interestingly, image analysis of infected mice reveals a concordant and statistically significant pattern: infection intensity clears in the bottom before the top of the lungs. The combined use of modeling and image analysis supports the application of phage therapy for acute P. aeruginosa pneumonia while emphasizing potential challenges to curative success in spatially structured in vivo environments, including impaired innate immune responses and reduced phage efficacy.

Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre.

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. We find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

Metapopulation model of phage therapy of an acute Pseudomonas aeruginosa lung infection.

Infections caused by multi-drug resistant (MDR) pathogenic bacteria are a global health threat. Phage therapy, which uses phage to kill bacterial pathogens, is increasingly used to treat patients infected by MDR bacteria. However, the therapeutic outcome of phage therapy may be limited by the emergence of phage resistance during treatment and/or by physical constraints that impede phage-bacteria interactions in vivo. In this work, we evaluate the role of lung spatial structure on the efficacy of phage therapy for Pseudomonas aeruginosa infection. To do so, we developed a spatially structured metapopulation network model based on the geometry of the bronchial tree, and included the emergence of phage-resistant bacterial mutants and host innate immune responses. We model the ecological interactions between bacteria, phage, and the host innate immune system at the airway (node) level. The model predicts the synergistic elimination of a P. aeruginosa infection due to the combined effects of phage and neutrophils given sufficiently active immune states and suitable phage life history traits. Moreover, the metapopulation model simulations predict that local MDR pathogens are cleared faster at distal nodes of the bronchial tree. Notably, image analysis of lung tissue time series from wild-type and lymphocyte-depleted mice (n=13) revealed a concordant, statistically significant pattern: infection intensity cleared in the bottom before the top of the lungs. Overall, the combined use of simulations and image analysis of in vivo experiments further supports the use of phage therapy for treating acute lung infections caused by P. aeruginosa while highlighting potential limits to therapy given a spatially structured environment, such as impaired innate immune responses and low phage efficacy.

Macrophage-induced reduction of bacteriophage density limits the efficacy of in vivo pulmonary phage therapy.

The rise of antimicrobial resistance has led to renewed interest in evaluating phage therapy. In murine models highly effective treatment of acute pneumonia caused by Pseudomonas aeruginosa relies on the synergistic antibacterial activity of bacteriophages with neutrophils. Here, we show that depletion of alveolar macrophages (AM) shortens the survival of mice without boosting the P. aeruginosa load in the lungs. Unexpectedly, upon bacteriophage treatment, pulmonary levels of P. aeruginosa were significantly lower in AM-depleted than in immunocompetent mice. To explore potential mechanisms underlying the benefit of AM-depletion in treated mice, we developed a mathematical model of phage, bacteria, and innate immune system dynamics. Simulations from the model fitted to data suggest that AM reduce bacteriophage density in the lungs. We experimentally confirmed that the in vivo decay of bacteriophage is faster in immunocompetent compared to AM-depleted animals. These findings demonstrate the involvement of feedback between bacteriophage, bacteria, and the immune system in shaping the outcomes of phage therapy in clinical settings.

Complex marine microbial communities partition metabolism of scarce resources over the diel cycle.

Complex assemblages of microbes in the surface ocean are responsible for approximately half of global carbon fixation. The persistence of high taxonomic diversity despite competition for a small suite of relatively homogeneously distributed nutrients, that is, 'the paradox of the plankton', represents a long-standing challenge for ecological theory. Here we find evidence consistent with temporal niche partitioning of nitrogen assimilation processes over a diel cycle in the North Pacific Subtropical Gyre. We jointly analysed transcript abundances, lipids and metabolites and discovered that a small number of diel archetypes can explain pervasive periodic dynamics. Metabolic pathway analysis of identified diel signals revealed asynchronous timing in the transcription of nitrogen uptake and assimilation genes among different microbial groups-cyanobacteria, heterotrophic bacteria and eukaryotes. This temporal niche partitioning of nitrogen uptake emerged despite synchronous transcription of photosynthesis and central carbon metabolism genes and associated macromolecular abundances. Temporal niche partitioning may be a mechanism by which microorganisms in the open ocean mitigate competition for scarce resources, supporting community coexistence.

Quantitative Models of Phage-Antibiotic Combination Therapy.

The spread of multidrug-resistant (MDR) bacteria is a global public health crisis. Bacteriophage therapy (or "phage therapy") constitutes a potential alternative approach to treat MDR infections. However, the effective use of phage therapy may be limited when phage-resistant bacterial mutants evolve and proliferate during treatment. Here, we develop a nonlinear population dynamics model of combination therapy that accounts for the system-level interactions between bacteria, phage, and antibiotics for in vivo application given an immune response against bacteria. We simulate the combination therapy model for two strains of Pseudomonas aeruginosa, one which is phage sensitive (and antibiotic resistant) and one which is antibiotic sensitive (and phage resistant). We find that combination therapy outperforms either phage or antibiotic alone and that therapeutic effectiveness is enhanced given interaction with innate immune responses. Notably, therapeutic success can be achieved even at subinhibitory concentrations of antibiotics, e.g., ciprofloxacin. These in silico findings provide further support to the nascent application of combination therapy to treat MDR bacterial infections, while highlighting the role of innate immunity in shaping therapeutic outcomes.IMPORTANCE This work develops and analyzes a novel model of phage-antibiotic combination therapy, specifically adapted to an in vivo context. The objective is to explore the underlying basis for clinical application of combination therapy utilizing bacteriophage that target antibiotic efflux pumps in Pseudomonas aeruginosa In doing so, the paper addresses three key questions. How robust is combination therapy to variation in the resistance profiles of pathogens? What is the role of immune responses in shaping therapeutic outcomes? What levels of phage and antibiotics are necessary for curative success? As we show, combination therapy outperforms either phage or antibiotic alone, and therapeutic effectiveness is enhanced given interaction with innate immune responses. Notably, therapeutic success can be achieved even at subinhibitory concentrations of antibiotic. These in silico findings provide further support to the nascent application of combination therapy to treat MDR bacterial infections, while highlighting the role of system-level feedbacks in shaping therapeutic outcomes.