Isolation and characterization of 24 phages infecting the plant growth-promoting rhizobacterium Klebsiella sp. M5al.

Bacteriophages largely impact bacterial communities via lysis, gene transfer, and metabolic reprogramming and thus are increasingly thought to alter nutrient and energy cycling across many of Earth's ecosystems. However, there are few model systems to mechanistically and quantitatively study phage-bacteria interactions, especially in soil systems. Here, we isolated, sequenced, and genomically characterized 24 novel phages infecting Klebsiella sp. M5al, a plant growth-promoting, nonencapsulated rhizosphere-associated bacterium, and compared many of their features against all 565 sequenced, dsDNA Klebsiella phage genomes. Taxonomic analyses revealed that these Klebsiella phages belong to three known phage families (Autographiviridae, Drexlerviridae, and Straboviridae) and two newly proposed phage families (Candidatus Mavericviridae and Ca. Rivulusviridae). At the phage family level, we found that core genes were often phage-centric proteins, such as structural proteins for the phage head and tail and DNA packaging proteins. In contrast, genes involved in transcription, translation, or hypothetical proteins were commonly not shared or flexible genes. Ecologically, we assessed the phages' ubiquity in recent large-scale metagenomic datasets, which revealed they were not widespread, as well as a possible direct role in reprogramming specific metabolisms during infection by screening their genomes for phage-encoded auxiliary metabolic genes (AMGs). Even though AMGs are common in the environmental literature, only one of our phage families, Straboviridae, contained AMGs, and the types of AMGs were correlated at the genus level. Host range phenotyping revealed the phages had a wide range of infectivity, infecting between 1-14 of our 22 bacterial strain panel that included pathogenic Klebsiella and Raoultella strains. This indicates that not all capsule-independent Klebsiella phages have broad host ranges. Together, these isolates, with corresponding genome, AMG, and host range analyses, help build the Klebsiella model system for studying phage-host interactions of rhizosphere-associated bacteria.