Emergence of networks of shared restriction-modification systems in phage-bacteria ecosystems

1Restriction-modification (RM) systems are the most ubiquitous bacterial defense system against bacteriophages and an important part of controlling phage predation. Using genomic sequence data, we show that RM systems are often shared among bacterial strains in a structured way. Examining the network of interconnections between bacterial strains within each genus, we find that in many genera strains share more RM systems than expected from a random network. We also find that many genera have a larger than expected number of bacterial strains with unique RM systems. We use population dynamics models of closed and open phage-bacteria ecosystems to qualitatively understand the selection pressures that could lead to these non-random network structures with enhanced overlap or uniqueness. In our models we find that the phages impose a pressure that favours bacteria with more RM systems, and more overlap of RM systems with other strains, but in bacteria dominated states this is opposed by the increased cost to the growth rate of these bacteria. Similar to what we observe in the genome data, we find that two distinct bacterial strategies emerge - strains either have a larger overlap than expected, or they have more unique RM systems than one expects from a null model. The former strategy appears to dominate when the repertoire of available RM systems is smaller but the average number of RM systems per strain is larger.