Host Range Breadth Correlates with Genic Diversity in Honeybee Phages.

Bacteriophages can evolve rapidly. Mutation and recombination via horizontal gene transfer allow them to counter adaptive responses by microbial hosts. However, little is known about the genomic processes underlying phage evolution within an ecological context-especially within natural microbial communities. This is due in part to the difficulty in resolving aspects of phage ecology, such as host range. To better understand the interplay of phage ecology and evolution within natural microbial communities, we combined measures of phage host range in vivo with measures of genome evolution in order to infer the evolutionary pressures acting on phage genomes within individual honeybee worker microbiomes. We show that near-identical phage genomes, cooccurring across multiple honeybee colonies, exhibit large variation with respect to gene modules, despite retaining a highly similar core genome. Estimates of genic diversity suggest deviations from neutral evolutionary models and identify loci under putative diversifying selection. We then use HiC-resolved metagenomics and show that the honeybee gut contains a dense phage community that exhibits a wide degree of host range variation. This variation differed across individual metagenomes in both the number and phylogenetic distance of potential hosts. We show that common measures of genetic variation positively correlate with host range in bee-associated phages and that functional targets of diversifying selection are partitioned differently between broad or narrow host range phages. Our work underscores the high host range variation associated with phages within host-associated microbial communities and provides evidence that this variation impacts rates of phage evolution.