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.

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.

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, usingPseudomonas aeruginosaand 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.IMPORTANCEGlobally 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 pathogenPseudomonas 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.