Bacteriophages (phages) have coevolved with their bacterial hosts for billions of years. With the rise of antibiotic resistance, the significance of using phages in therapy is increasing. Investigating the dynamics of phage evolution can provide valuable insights for pre-adapting phages to more challenging clones of their hosts that may arise during treatment. Two primary models describe interactions in phage–bacteria systems: arms race dynamics and fluctuating selection dynamics. Numerous factors influence which dynamics dominate the interactions between a phage and its host. These dynamics, in turn, affect the coexistence of phages and bacteria, ultimately determining which organism will adapt more effectively to the other, and whether a stable state will be reached. In this review, we summarize key findings from research on phage–bacteria coevolution, focusing on the different concepts that can describe these interactions, the factors that may contribute to the prevalence of one model over others, and the effects of various dynamics on both phages and bacteria.

Aeromonas popoffii is one of the environmental Aeromonas species. A number of factors of virulence have been described for this species and it has been reported as a causative agent of urinary tract infection. The first A. popoffii bacteriophage AerP_220 along with its host strain A. popoffii CEMTC 4062 were isolated from river water. The phage has a podovirus morphotype, shows a narrow host range and is lytic against the host strain. The AerP_220 genome comprises 45,207 bp and does not contain genes responsible for antibiotic resistance and toxin production. Fifty-nine co-directional putative ORFs were found in the AerP_220 genome. Thirty-three ORFs encoded proteins with predicted functions; the products of 26 ORFs were hypothetical proteins. AerP_220 genome analysis revealed that this phage can be considered a novel species within the Autographiviridae family. Comparative genomic and proteomic analysis revealed that AerP_220 along with the Aeromonas phage vB_AspA_Tola (OM913599) are members of a new putative Tolavirus genus in the family Autographiviridae. The Gajwadongvirus and proposed Tolavirus genera along with Pantoea phage Nufs112 and phage Reminis could form a new Tolavirinae subfamily within the Autographiviridae family.

A novel Enterobacter cloacae phage, EC151, was isolated and characterized. Electron microscopy revealed that EC151 has a siphovirus-like virion morphology. The EC151 nucleotide sequence shows limited similarity to other phage genomes deposited in the NCBI GenBank database. The size of the EC151 genome is 60,753 bp and contains 58 putative genes. Thirty-nine of them encode proteins of predicted function, 18 are defined as hypothetical proteins, and one ORF identifies as the tRNA-Ser-GCT-encoding gene. Six ORFs were predicted to be members of the deazaguanine DNA modification pathway, including the preQ₀ transporter. Comparative proteomic phylogenetic analysis revealed that phage EC151 represents a distinct branch within a group of sequences containing clades formed by members of the Seuratvirus, Nonagvirus, and Vidquintavirus genera. In addition, the EC151 genome showed gene synteny typical of the Seuratvirus, Nonagvirus, and Nipunavirus phages. The average genetic distances of EC151/Seuratvirus, EC151/Nonagvirus, and EC151/Vidquintavirus are approximately equal to those between the Seuratvirus, Nonagvirus, and Vidquintavirus genera (~0.7 substitutions per site). Therefore, EC151 may represent a novel genus within the Siphoviridae family. The origin of the deazaguanine DNA modification pathway in the EC151 genome can be traced to Escherichia phages from the Seuratvirus genus.

Stenotrophomonas maltophilia mainly causes respiratory infections that are associated with a high mortality rate among immunocompromised patients. S. maltophilia exhibits a high level of antibiotic resistance and can form biofilms, which complicates the treatment of patients infected with this bacterium. Phages combined with antibiotics could be a promising treatment option. Currently, ~60 S. maltophilia phages are known, and their effects on biofilm formation and antibiotic sensitivity require further examination. Bacteriophage StM171, which was isolated from hospital wastewater, showed a medium host range, low burst size, and low lytic activity. StM171 has a 44kbp dsDNA genome that encodes 59 open-reading frames. A comparative genomic analysis indicated that StM171, along with the Stenotrophomonas phage Suso (MZ326866) and Xanthomonas phage HXX_Dennis (ON711490), are members of a new putative Nordvirus genus. S. maltophilia strains that developed resistance to StM171 (bacterial-insensitive mutants) showed a changed sensitivity to antibiotics compared to the originally susceptible strains. Some bacterial-insensitive mutants restored sensitivity to cephalosporin and penicillin-like antibiotics and became resistant to erythromycin. StM171 shows strain- and antibiotic-dependent effects on the biofilm formation of S. maltophilia strains.

Anti-phage defense systems are widespread in bacteria due to the latter continuous adaptation to infection by bacteriophages (phages). Stenotrophomonas maltophilia has a high degree of intrinsic antibiotic resistance, which makes phage therapy relevant for the treatment of infections caused by this species. Studying the array of anti-phage defense systems that could be found in S. maltophilia helps in better adapting the phages to the systems present in the pathogenic bacteria. Pangenome analysis of the available S. maltophilia strains with complete genomes that were downloaded from GenBank, including five local genomes, indicated a wide set of 72 defense systems and subsystems that varied between the strains. Seven of these systems were present in more than 20% of the studied genomes and the proteins encoded by the systems were variable in most of the cases. A total of 27 defense islands were revealed where defense systems were found; however, more than 60% of the instances of systems were found in four defense islands. Several elements linked to the transfer of these systems were found. No obvious associations between the pattern of distribution of the anti-phage defense systems of S. maltophilia and the phylogenetic features or the isolation site were found.

Two novel P. protegens bacteriophages PseuP_222 and Pseu_224 and their host P. protegens CEMTC 4060 were isolated from the same sample (Inya river, Siberia). Both phages have siphovirus morphology and belong to lambdoid phages. Comparative genome analysis revealed a low nucleotide and amino acid sequence similarity of PseuP_222 and PseuP_224 between themselves, and between them and other lambdoid phages. Bioinformatics analysis indicated that PseuP_222 and PseuP_224 are members of a genetically diverse group of phages of environmental Pseudomonas spp.; this group is distant from a large group of P. aeruginosa phages. In phylogenetic trees, the positioning of the terminase large subunits, major capsid proteins, tail tape measure proteins, and CI-like repressors of PseuP_222 and PseuP_224 were remote and changed relative to those of the Escherichia lambda phage and lambdoid phages of Pseudomonas spp. However, the nucleoid-associated protein NdpA/YejK and P5-like structural protein from both phages showed high similarity and were not found in lambda phage and other lambdoid phages of Pseudomonas spp. Substantial divergences of the PseuP_222 and PseuP_224 genomes and proteomes indicated that the evolutionary history of these phages was mostly independent and they probably began to use one host only recently.