Within-host environments are complex and multidimensional, making it challenging to link the evolutionary responses of colonizing pathogens to causal selective drivers. Loss of quorum sensing (QS) via mutation of the master regulator, lasR, is a common adaptation of Pseudomonas aeruginosa during chronic infections. Here, we use experimental evolution in host-mimicking media to show that loss of QS is constrained by environmental factors associated with host inflammation. Specifically, environments combining oxidative stress and abundant free amino acids limited loss of QS, whereas QS loss was rapid in the absence of oxidative stress, regardless of amino acids. Under oxidative stress, lasR mutations were contingent upon first decoupling regulation of oxidative stress responses from QS via mutations in the promoter region of the primary catalase, katA, or in the oxidative stress regulator, oxyR, enabling maintenance of oxidative stress tolerance. Together, our findings suggest that host inflammatory responses likely delay the loss of QS while colonizers undergo stepwise evolution, first adapting to survive lethal stressors before responding to other nutritional selective drivers that favor loss of QS. IMPORTANCE:Pseudomonas aeruginosa is a common cause of chronic infections characterized by persistent inflammation. Host inflammatory responses alter within-host environments, including by increasing levels of antimicrobial stressors and releasing free amino acids through proteolysis. Here, we show stepwise adaptation of experimental P. aeruginosa populations to inflammation-like environments, first adapting to survive lethal stress by decoupling oxidative stress responses from quorum sensing (QS), before then adapting to the nutritional conditions, delaying the loss of quorum sensing. These results highlight the power of using laboratory evolution experiments to disentangle the multidimensional selective forces driving pathogen adaptation in complex within-host environments.

Pseudomonas aeruginosa is a major nosocomial pathogen that causes severe disease including sepsis. Carbapenem-resistant P. aeruginosa is recognised by the World Health Organisation as a priority 1 pathogen, with urgent need for new therapeutics. As such, there is renewed interest in using bacteriophages as a therapeutic. However, the dynamics of treating pan-resistant P. aeruginosa with phage in vivo are poorly understood. Using a pan-resistant P. aeruginosa in vivo infection model, phage therapy displays strong therapeutic potential, clearing infection from the blood, kidneys, and spleen. Remaining bacteria in the lungs and liver displays phage resistance due to limiting phage adsorption. Yet, resistance to phage results in re-sensitisation to a wide range of antibiotics. In this work, we use phage steering in vivo, pre-exposing a pan resistant P. aeruginosa infection with a phage cocktail to re-sensitise bacteria to antibiotics, clearing the infection from all organs.

Clostridioides difficile is an important human pathogen, for which there are very limited treatment options, primarily the glycopeptide antibiotic vancomycin. In recent years, vancomycin resistance has emerged as a serious problem in several gram-positive pathogens, but high-level resistance has yet to be reported for C. difficile, although it is not known if this is due to constraints upon resistance evolution in this species. Here, we show that resistance to vancomycin can evolve rapidly under ramping selection but is accompanied by fitness costs and pleiotropic trade-offs, including sporulation defects that would be expected to severely impact transmission. We identified 2 distinct pathways to resistance, both of which are predicted to result in changes to the muropeptide terminal D-Ala-D-Ala that is the primary target of vancomycin. One of these pathways involves a previously uncharacterised D,D-carboxypeptidase, expression of which is controlled by a dedicated two-component signal transduction system. Our findings suggest that while C. difficile is capable of evolving high-level vancomycin resistance, this outcome may be limited clinically due to pleiotropic effects on key pathogenicity traits. Moreover, our data identify potential mutational routes to resistance that should be considered in genomic surveillance.

Bacteria have evolved a variety of defence mechanisms to protect against mobile genetic elements, including restrictionmodification systems and CRISPR-Cas. In recent years, dozens of previously unknown defence systems (DSs) have been discovered. Notably, diverse DSs often coexist within the same genome, and some co-occur at frequencies significantly higher than would be expected by chance, implying potential synergistic interactions. Recent studies have provided evidence of defence mechanisms that enhance or complement one another. Here, we review the interactions between DSs at the mechanistic, regulatory, ecological and evolutionary levels.

Plasmids are important vectors of horizontal gene transfer in microbial communities but can impose a burden on the bacteria that carry them. Such plasmid fitness costs are thought to arise principally from conflicts between chromosomal- and plasmid-encoded molecular machineries, and thus can be ameliorated by compensatory mutations (CMs) that reduce or resolve the underlying causes. CMs can arise on plasmids (i.e., plaCM) or on chromosomes (i.e., chrCM), with contrasting predicted effects upon plasmid success and subsequent gene transfer because plaCM can also reduce fitness costs in plasmid recipients, whereas chrCM can potentially ameliorate multiple distinct plasmids. Here, we develop theory and a novel experimental system to directly compare the ecological effects of plaCM and chrCM that arose during evolution experiments between Pseudomonas fluorescens SBW25 and its sympatric mercury resistance megaplasmid pQBR57. We show that while plaCM was predicted to succeed under a broader range of parameters in mathematical models, chrCM dominated in our experiments, including conditions with numerous recipients, due to a more efficacious mechanism of compensation, and advantages arising from transmission of costly plasmids to competitors (plasmid “weaponisation”). We show analytically the presence of a mixed Rock-Paper-Scissors (RPS) regime for CMs, driven by trade-offs with horizontal transmission, that offers one possible explanation for the observed failure of plaCM to dominate even in competition against an uncompensated plasmid. Our results reveal broader implications of plasmid-bacterial evolution for plasmid ecology, demonstrating the importance of specific compensatory mutations for resistance gene spread. One consequence of the superiority of chrCM over plaCM is the likely emergence in microbial communities of compensated bacteria that can act as “hubs” for plasmid accumulation and dissemination.

Plasmids are important vectors of horizontal gene transfer in microbial communities but can impose a burden on the bacteria that carry them. Such plasmid fitness costs are thought to arise principally from conflicts between chromosomal- and plasmid-encoded molecular machineries, and thus can be ameliorated by compensatory mutations (CMs) that reduce or resolve the underlying causes. CMs can arise on plasmids (i.e. plaCM) or on chromosomes (i.e. chrCM), with contrasting predicted effects upon plasmid success and subsequent gene transfer because plaCM can also reduce fitness costs in plasmid recipients, whereas chrCM can potentially ameliorate multiple distinct plasmids. Here, we develop theory and a novel experimental system to directly compare the ecological effects of plaCM and chrCM that arose during evolution experiments between Pseudomonas fluorescens SBW25 and its sympatric mercury resistance megaplasmid pQBR57. We show that while plaCM was predicted to succeed under a broader range of parameters in mathematical models, experimentally chrCM dominated under all conditions, including those with numerous recipients, due to a more efficacious mechanism of compensation, and advantages arising from transmission of costly plasmids to competitors (plasmid ‘weaponisation’). We show analytically the presence of a mixed Rock-Paper-Scissors regime for plaCM, driven by trade-offs with horizontal transmission, that explains the observed failure of plaCM to dominate even in competition against an uncompensated plasmid. Our results reveal broader implications of plasmid-bacterial evolution for plasmid ecology, demonstrating the importance of compensatory mutations for resistance gene spread. One consequence of the superiority of chrCM over plaCM is the likely emergence in microbial communities of compensated bacteria that can act as ‘hubs’ for plasmid accumulation and dissemination.### Competing Interest StatementThe authors have declared no competing interest.

Clostridioides difficile is an important human pathogen, for which there are very limited treatment options, primarily the glycopeptide antibiotic vancomycin. In recent years vancomycin resistance has emerged as a serious problem in several Gram positive pathogens, but high level resistance has yet to be reported for C. difficile , although it is not known if this is due to constraints upon resistance evolution in this species. Here we show that resistance to vancomycin can evolve rapidly under ramping selection but is accompanied by severe fitness costs and pleiotropic trade-offs, including sporulation defects that would be expected to severely impact transmission. We identified two distinct pathways to resistance, both of which are predicted to result in changes to the muropeptide terminal D-Ala-D-Ala that is the primary target of vancomycin. One of these pathways involves a previously uncharacterised D,D-carboxypeptidase, expression of which is controlled by a dedicated two-component signal transduction system. Our findings suggest that while C. difficile is capable of evolving high-level vancomycin resistance, this outcome may be limited clinically due to pleiotropic effects on key pathogenicity trains. Moreover, our data provide a mutational roadmap to inform genomic surveillance.

Competitive bacteria-bacteriophage interactions have resulted in the evolution of a plethora of bacterial defense systems preventing phage propagation. In recent years, computational and bioinformatic approaches have underpinned the discovery of numerous novel bacterial defense systems. Anti-phage systems are frequently encoded together in genomic loci termed defense islands. Here we report the identification and characterisation of a novel anti-phage system, that we have termed Shield, which forms part of the Pseudomonas defensive arsenal. The Shield system comprises the core component ShdA, a membrane-bound protein harboring an RmuC domain. Heterologous production of ShdA alone is sufficient to mediate bacterial immunity against several phages. We demonstrate that Shield and ShdA confer population-level immunity and that they can also decrease transformation efficiency. We further show that ShdA homologues can degrade DNA in vitro and, when expressed in a heterologous host, can alter the organisation of the host chromosomal DNA. Use of comparative genomic approaches identified how Shield can be divided into four subtypes, three of which contain additional components that in some cases can negatively affect the activity of ShdA and/or provide additional lines of phage defense. Collectively, our results identify a new player within the Pseudomonas bacterial immunity arsenal that displays a novel mechanism of protection, and reveals a role for RmuC domains in phage defense.

Plasmids are well known for transferring antibiotic resistance genes between bacteria. A study in the clinic shows that evolutionary dynamics within the gut microbiomes of hospitalized patients lead to rapid adaptive changes, balancing the level of resistance that plasmids provide against the fitness costs that they impose on bacteria.

Bacterial pathogens are confronted with a range of challenges at the site of infection, including exposure to antibiotic treatment and harsh physiological conditions, that can alter the fitness benefits and costs of acquiring antibiotic resistance. Here, we develop an experimental system to recapitulate resistance gene acquisition by Staphylococcus aureus and test how the subsequent evolution of the resistant bacterium is modulated by antibiotic treatment and oxygen levels, both of which are known to vary extensively at sites of infection.

Plasmids play an important role in bacterial genome evolution by transferring genes between lineages. Fitness costs associated with plasmid carriage are expected to be a barrier to gene exchange, but the causes of plasmid fitness costs are poorly understood. Single compensatory mutations are often sufficient to completely ameliorate plasmid fitness costs, suggesting that such costs are caused by specific genetic conflicts rather than generic properties of plasmids, such as their size, metabolic burden, or gene expression level. By combining the results of experimental evolution with genetics and transcriptomics, we show here that fitness costs of 2 divergent large plasmids in Pseudomonas fluorescens are caused by inducing maladaptive expression of a chromosomal tailocin toxin operon. Mutations in single genes unrelated to the toxin operon, and located on either the chromosome or the plasmid, ameliorated the disruption associated with plasmid carriage. We identify one of these compensatory loci, the chromosomal gene PFLU4242 , as the key mediator of the fitness costs of both plasmids, with the other compensatory loci either reducing expression of this gene or mitigating its deleterious effects by up-regulating a putative plasmid-borne ParAB operon. The chromosomal mobile genetic element Tn6291, which uses plasmids for transmission, remained up-regulated even in compensated strains, suggesting that mobile genetic elements communicate through pathways independent of general physiological disruption. Plasmid fitness costs caused by specific genetic conflicts are unlikely to act as a long-term barrier to horizontal gene transfer (HGT) due to their propensity for amelioration by single compensatory mutations, helping to explain why plasmids are so common in bacterial genomes.

In Belgium, the incorporation of phages into magistral preparations for human application has been permitted since 2018. The stability of such preparations is of high importance to guarantee quality and efficacy throughout treatments. We evaluated the ability to preserve infectivity of four different phages active against three different bacterial species in five different buffer and infusion solutions commonly used in medicine and biotechnological manufacturing processes, at two different concentrations (9 and 7 log pfu/mL), stored at 4 °C. DPBS without Ca2+ and Mg2+ was found to be the best option, compared to the other solutions. Suspensions with phage concentrations of 7 log pfu/mL were unsuited as their activity dropped below the effective therapeutic dose (6–9 log pfu/mL), even after one week of storage at 4 °C. Strong variability between phages was observed, with Acinetobacter baumannii phage Acibel004 being stable in four out of five different solutions. We also studied the long term storage of lyophilized staphylococcal phage ISP, and found that the titer could be preserved during a period of almost 8 years when sucrose and trehalose were used as stabilizers. After rehydration of the lyophilized ISP phage in saline, the phage solutions remained stable at 4 °C during a period of 126 days.

Phage therapy is a promising alternative to traditional antibiotics for treating bacterial infections. Such phage-based therapeutics typically contain multiple phages, but how the efficacy of phage combinations scales with phage richness, identity and functional traits is unclear. Here, we experimentally tested the efficacy of 827 unique phage combinations ranging in phage richness from one to 12 phages. The efficacy of phage combinations increased with phage richness. However, complementarity between functionally diverse phages allowed efficacy to be maximized at lower levels of phage richness in functionally diverse combinations. These findings suggest that phage functional diversity is the key property of effective phage combinations, enabling the design of simple but effective phage therapies that overcome the practical and regulatory hurdles that limit development of more diverse phage therapy cocktails.

This document contains all experiment data, including a) population dynamics, b) RBG and c) Growth rate data. (a) Dynamics of the bacteria and phage populations through time. Every three transfers, bacteria and phage population densities were measured. (b) Resistance estimates as 'Reduction in Bacterial Growth' (RBG) for 20 bacterial clones against whole phage populations. Both bacteria and phage were isolated from the final transfer for four replicates of four (out of five) treatments. (c) Growth rate data for ancestral and phage-resistant mutant strains without phage in different media concentrations. Data are OD600nm measurements taken every 30 minutes for 48 hours.

Seasonal environments vary in their amplitude of oscillation but the effects of this temporal heterogeneity for host–parasite coevolution are poorly understood. Here, we combined mathematical modelling and experimental evolution of a coevolving bacteria–phage interaction to show that the intensity of host–parasite coevolution peaked in environments that oscillate in their resource supply with intermediate amplitude. Our experimentally parameterized mathematical model explains that this pattern is primarily driven by the ecological effects of resource oscillations on host growth rates. Our findings suggest that in host–parasite systems where the host's but not the parasite's population growth dynamics are subject to seasonal forcing, the intensity of coevolution will peak at intermediate amplitudes but be constrained at extreme amplitudes of environmental oscillation.

The eco-evolutionary dynamics of microbial communities are predicted to affect both the tempo and trajectory of evolution in constituent species [1]. While community composition determines available niche space, species sorting dynamically alters composition, changing over time the distribution of vacant niches to which species adapt [2], altering evolutionary trajectories [3, 4]. Competition for the same niche can limit evolutionary potential if population size and mutation supply are reduced [5, 6] but, alternatively, could stimulate evolutionary divergence to exploit vacant niches if character displacement results from the coevolution of competitors [7, 8]. Under more complex ecological scenarios, species can create new niches through their exploitation of complex resources, enabling others to adapt to occupy these newly formed niches [9, 10]. Disentangling the drivers of natural selection within such communities is extremely challenging, and it is thus unclear how eco-evolutionary dynamics drive the evolution of constituent taxa. We tracked the metabolic evolution of a focal species during adaptation to wheat straw as a resource both in monoculture and in polycultures wherein on-going eco-evolutionary community dynamics were either permitted or prevented. Species interactions accelerated metabolic evolution. Eco-evolutionary dynamics drove increased use of recalcitrant substrates by the focal species, whereas greater exploitation of readily digested substrate niches created by other species evolved if on-going eco-evolutionary dynamics were prevented. Increased use of recalcitrant substrates was associated with parallel evolution of tctE, encoding a carbon metabolism regulator. Species interactions and species sorting set, respectively, the tempo and trajectory of evolutionary divergence among communities, selecting distinct ecological functions in otherwise equivalent ecosystems.

Phage therapy is a promising alternative to chemotherapeutic antibiotics for the treatment of bacterial infections. However, despite recent clinical uses of combinations of phages to treat multidrug-resistant infections, a mechanistic understanding of how bacteria evolve resistance against multiple phages is lacking, limiting our ability to deploy phage combinations optimally. Here, we show, using Pseudomonas aeruginosa and pairs of phages targeting shared or distinct surface receptors, that the timing and order of phage exposure determine the strength, cost, and mutational basis of resistance. Whereas sequential exposure allowed bacteria to acquire multiple resistance mutations effective against both phages, this evolutionary trajectory was prevented by simultaneous exposure, resulting in quantitatively weaker resistance. The order of phage exposure determined the fitness costs of sequential resistance, such that certain sequential orders imposed much higher fitness costs than the same phage pair in the reverse order. Together, these data suggest that phage combinations can be optimized to limit the strength of evolved resistances while maximizing their associated fitness costs to promote the long-term efficacy of phage therapy.IMPORTANCE Globally rising rates of antibiotic resistance have renewed interest in phage therapy where combinations of phages have been successfully used to treat multidrug-resistant infections. To optimize phage therapy, we first need to understand how bacteria evolve resistance against combinations of multiple phages. Here, we use simple laboratory experiments and genome sequencing to show that the timing and order of phage exposure determine the strength, cost, and mutational basis of resistance evolution in the opportunistic pathogen Pseudomonas aeruginosa These findings suggest that phage combinations can be optimized to limit the emergence and persistence of resistance, thereby promoting the long-term usefulness of phage therapy.

Temperate bacteriophages are a common feature of Pseudomonas aeruginosa genomes, but their role in chronic lung infections is poorly understood. This study was designed to identify the diverse communities of mobile P. aeruginosa phages by employing novel metagenomic methods, to determine cross infectivity, and to demonstrate the influence of phage infection on antimicrobial susceptibility. Mixed temperate phage populations were chemically mobilized from individual P. aeruginosa, isolated from patients with cystic fibrosis (CF) or bronchiectasis (BR). The infectivity phenotype of each temperate phage lysate was evaluated by performing a cross-infection screen against all bacterial isolates and tested for associations with clinical variables. We utilized metagenomic sequencing data generated for each phage lysate and developed a novel bioinformatic approach allowing resolution of individual temperate phage genomes. Finally, we used a subset of the temperate phages to infect P. aeruginosa PAO1 and tested the resulting lysogens for their susceptibility to antibiotics. Here, we resolved 105 temperate phage genomes from 94 lysates that phylogenetically clustered into 8 groups. We observed disease-specific phage infectivity profiles and found that phages induced from bacteria isolated from more advanced disease infected broader ranges of P. aeruginosa isolates. Importantly, when infecting PAO1 in vitro with 20 different phages, 8 influenced antimicrobial susceptibility. This study shows that P. aeruginosa isolated from CF and BR patients harbors diverse communities of inducible phages, with hierarchical infectivity profiles that relate to the progression of the disease. Temperate phage infection altered the antimicrobial susceptibility of PAO1 at subinhibitory concentrations of antibiotics, suggesting they may be precursory to antimicrobial resistance. IMPORTANCE Pseudomonas aeruginosa is a key opportunistic respiratory pathogen in patients with cystic fibrosis and non-cystic fibrosis bronchiectasis. The genomes of these pathogens are enriched with mobile genetic elements including diverse temperate phages. While the temperate phages of the Liverpool epidemic strain have been shown to be active in the human lung and enhance fitness in a rat lung infection model, little is known about their mobilization more broadly across P. aeruginosa in chronic respiratory infection. Using a novel metagenomic approach, we identified eight groups of temperate phages that were mobilized from 94 clinical P. aeruginosa isolates. Temperate phages from P. aeruginosa isolated from more advanced disease showed high infectivity rates across a wide range of P. aeruginosa genotypes. Furthermore, we showed that multiple phages altered the susceptibility of PAO1 to antibiotics at subinhibitory concentrations.

The acquisition of plasmids is often accompanied by fitness costs such that compensatory evolution is required to allow plasmid survival, but it is unclear whether compensatory evolution can be extensive or rapid enough to maintain plasmids when they are very costly. The mercury-resistance plasmid pQBR55 drastically reduced the growth of its host, Pseudomonas fluorescens SBW25, immediately after acquisition, causing a small colony phenotype. However, within 48h of growth on agar plates we observed restoration of the ancestral large colony morphology, suggesting that compensatory mutations had occurred. Relative fitness of these evolved strains, in lab media and in soil microcosms, varied between replicates, indicating different mutational mechanisms. Using genome sequencing we identified that restoration was associated with chromosomal mutations in either a hypothetical DNA-binding protein PFLU4242, RNA polymerase or the GacA/S two-component system. Targeted deletions in PFLU4242, gacA or gacS recapitulated the ameliorated phenotype upon plasmid acquisition, indicating three distinct mutational pathways to compensation. Our data shows that plasmid compensatory evolution is fast enough to allow survival of a plasmid despite it imposing very high fitness costs upon its host, and indeed may regularly occur during the process of isolating and selecting individual plasmid-containing clones.

Phages shape the structure of natural bacterial communities and can be effective therapeutic agents. Bacterial resistance to phage infection, however, limits the usefulness of phage therapies and could destabilise community structures, especially if individual resistance mutations provide cross-resistance against multiple phages. We currently understand very little about the evolution of cross-resistance in bacteria-phage interactions. Here we show that the network structure of cross-resistance among spontaneous resistance mutants of Pseudomonas aeruginosa evolved against each of 27 phages is highly modular. The cross-resistance network contained both symmetric (reciprocal) and asymmetric (nonreciprocal) cross-resistance, forming two cross-resistance modules defined by high within-but low between-module cross-resistance. Mutations conferring cross-resistance within modules targeted either lipopolysaccharide or type IV pilus biosynthesis, suggesting that the modularity of cross-resistance was structured by distinct phage receptors. In contrast, between-module cross-resistance was provided by mutations affecting the alternative sigma factor, RpoN, which controls many lifestyle-associated functions, including motility, biofilm formation, and quorum sensing. Broader cross-resistance range was not associated with higher fitness costs or weaker resistance against the focal phage used to select resistance. However, mutations in rpoN, providing between-module cross-resistance, were associated with higher fitness costs than mutations associated with within-module cross-resistance, i.e., in genes encoding either lipopolysaccharide or type IV pilus biosynthesis. The observed structure of cross-resistance predicted both the frequency of resistance mutations and the ability of phage combinations to suppress bacterial growth. These findings suggest that the evolution of cross-resistance is common, is likely to play an important role in the dynamic structure of bacteria-phage communities, and could inform the design principles for phage therapy treatments.