Combinatorial phenotypic landscape enables bacterial resistance to phage infection

Success of phage therapies is limited by bacterial defenses against phages. While a large variety of anti- phage defense mechanisms has been characterized, how expression of these systems is distributed across individual cells and how their combined activities translate into protection from phages has not been studied. Using bacterial single-cell RNA sequencing, we profiled the transcriptomes of ∼50,000 cells from cultures of a human pathobiont, Bacteroides fragilis, infected with a lytic bacteriophage. We quantified the asynchronous progression of phage infection in single bacterial cells and reconstructed the infection timeline, characterizing both host and phage transcriptomic changes as infection unfolded. We discovered a subpopulation of bacteria that remained uninfected and determined the heterogeneously expressed host factors associated with protection. Each cell’s vulnerability to phage infection was defined by combinatorial phase-variable expression of multiple genetic loci, including capsular polysaccharide (CPS) biosynthesis pathways, restriction-modification systems (RM), and a previously uncharacterized operon likely encoding fimbrial genes. By acting together, these heterogeneously expressed phase-variable systems and anti-phage defense mechanisms create a phenotypic landscape where distinct protective combinations enable the survival and re-growth of bacteria expressing these phenotypes without acquiring additional mutations. The emerging model of complementary action of multiple protective mechanisms heterogeneously expressed across an isogenic bacterial population showcases the potent role of phase variation and stochasticity in bacterial anti-phage defenses.