Despite phages’ ubiquitous presence and great importance in shaping microbial communities, little is known about the diversity of specific phages in different ecological niches. Here, we isolated, sequenced, and characterized 38 Escherichia coli-infecting phages (coliphages) from poultry faeces to gain a better understanding of the coliphage diversity in the poultry intestine. All phages belonged to either the Siphoviridae or Myoviridae family and their genomes ranged between 44,324 and 173,384 bp, with a G+C content between 35.5 and 46.4%. Phylogenetic analysis was performed based on single “marker” genes; the terminase large subunit, portal protein, and exonucleases, as well as the full draft genomes. Single gene analysis resulted in six distinct clusters. Only minor differences were observed between the different phylogenetic analyses, including branch lengths and additional duplicate or triplicate subclustering. Cluster formation was according to genome size, G+C content and phage subfamily. Phylogenetic analysis based on the full genomes supported these clusters. Moreover, several of our Siphoviridae phages might represent a novel unclassified phage genus. This study allowed for identification of several novel coliphages and provides new insights to the coliphage diversity in the intestine of poultry. Great diversity was observed amongst the phages, while they were isolated from an otherwise similar ecosystem.

Avian pathogenic Escherichia coli (APEC) is one of the most important bacterial pathogens affecting poultry worldwide. The emergence of multidrug-resistant pathogens has renewed the interest in the therapeutic use of bacteriophages (phages). However, a major concern for the successful implementation of phage therapy is the emergence of phage-resistant mutants. The understanding of the phage-host interactions, as well as underlying mechanisms of resistance, have shown to be essential for the development of a successful phage therapy. Here, we demonstrate that the strictly lytic Escherichia phage vB_EcoM-P10 rapidly selected for resistance in the APEC ST95 O1 strain AM621. Whole-genome sequence analysis of 109 spontaneous phage-resistant mutant strains revealed 41 mutants with single-nucleotide polymorphisms (SNPs) in their core genome. In 32 of these, a single SNP was detected while two SNPs were identified in a total of nine strains. In total, 34 unique SNPs were detected. In 42 strains, including 18 strains with SNP(s), gene losses spanning 17 different genes were detected. Affected by genetic changes were genes known to be involved in phage resistance (outer membrane protein A, lipopolysaccharide-, O- antigen-, or cell wall-related genes) as well as genes not previously linked to phage resistance, including two hypothetical genes. In several strains, we did not detect any genetic changes. Infecting phages were not able to overcome the phage resistance in host strains. However, interestingly the initial infection was shown to have a great fitness cost for several mutant strains, with up to ∼65% decrease in overall growth. In conclusion, this study provides valuable insights into the phage-host interaction and phage resistance in APEC. Although acquired resistance to phages is frequently observed in pathogenic E. coli, it may be associated with loss of fitness, which could be exploited in phage therapy.

The therapeutic use of bacteriophages (phage therapy) represents a promising alternative to antibiotics to control bacterial pathogens. However, the understanding of the phage–bacterium interactions and population dynamics seems essential for successful phage therapy implementation. Here, we investigated the effect of three factors: phage species (18 lytic E. coli-infecting phages); bacterial strain (10 APEC strains); and multiplicity of infection (MOI) (MOI 10, 1, and 0.1) on the bacterial growth dynamics. All factors had a significant effect, but the phage appeared to be the most important. The results showed seven distinct growth patterns. The first pattern corresponded to the normal bacterial growth pattern in the absence of a phage. The second pattern was complete bacterial killing. The remaining patterns were in-between, characterised by delayed growth and/or variable killing of the bacterial cells. In conclusion, this study demonstrates that the phage–host dynamics is an important factor in the capacity of a phage to eliminate bacteria. The classified patterns show that this is an essential factor to consider when developing a phage therapy. This methodology can be used to rapidly screen for novel phage candidates for phage therapy. Accordingly, the most promising candidates were phages found in Group 2, characterised by growth dynamics with high bacterial killing.