Species-scale genomic analysis of S. aureus genes influencing phage host range and their relationships to virulence and antibiotic resistance genes

Phage therapy has been proposed as a possible alternative treatment for infections caused by the ubiquitous bacterial pathogen Staphylococcus aureus. However, successful phage therapy requires knowing both host and phage genetic factors influencing host range for rational cocktail formulation. To further our understanding of host range, we searched 40,000+ public S. aureus genome sequences for previously identified phage resistance genes. We found that phage adsorption targets and genes that block phage assembly were significantly more conserved than genes targeting phage biosynthesis. Core phage resistance genes had similar nucleotide diversity, ratio of non-synonymous to synonymous substitutions, and functionality (measured by delta-bitscore) to other core genes in a set of 380 non-redundant S. aureus genomes (each from a different MLST sequence type). Non-core phage resistance genes were significantly less consistent with the core genome phylogeny than all non-core genes in this set. Only superinfection immunity genes correlated with empirically determined temperate phage resistance, accessory genome content, and numbers of accessory antibiotic resistance or virulence genes encoded per strain. Taken together, these results suggested that, while phage adsorption genes are heavily conserved in the S. aureus species, they are not undergoing positive selection, arms race dynamics. They also suggested genes classified as involved in assembly are least phylogenetically constrained and superinfection immunity genes best predict both empirical phage resistance and levels of phage-mediated HGT. ImportanceStaphylococcus aureus is a widespread, hospital- and community-acquired pathogen that is commonly antibiotic resistant. It causes diverse diseases affecting both the skin and internal organs. Its ubiquity, antibiotic resistance, and disease burden make new therapies urgent, such as phage therapy, in which viruses specific to infecting bacteria clear infection. S. aureus phage host range not only determines whether phage therapy will be successful by killing bacteria but also horizontal gene transfer through transduction of host genetic material by phages. In this work, we comprehensively reviewed existing literature to build a list of S. aureus phage resistance genes and searched our database of almost 43,000 S. aureus genomes for these genes to understand their patterns of evolution, finding that prophages superinfection immunity correlates best with phage resistance and HGT. These findings improved our understanding of the relationship between known phage resistance genes and phage host range in the species.