Genome-scale reconstruction of Gcn4/ATF4 networks driving a growth program.

Growth and starvation are considered opposite ends of a spectrum. To sustain growth, cells use coordinated gene expression programs and manage biomolecule supply in order to match the demands of metabolism and translation. Global growth programs complement increased ribosomal biogenesis with sufficient carbon metabolism, amino acid and nucleotide biosynthesis. How these resources are collectively managed is a fundamental question. The role of the Gcn4/ATF4 transcription factor has been best studied in contexts where cells encounter amino acid starvation. However, high Gcn4 activity has been observed in contexts of rapid cell proliferation, and the roles of Gcn4 in such growth contexts are unclear. Here, using a methionine-induced growth program in yeast, we show that Gcn4/ATF4 is the fulcrum that maintains metabolic supply in order to sustain translation outputs. By integrating matched transcriptome and ChIP-Seq analysis, we decipher genome-wide direct and indirect roles for Gcn4 in this growth program. Genes that enable metabolic precursor biosynthesis indispensably require Gcn4; contrastingly ribosomal genes are partly repressed by Gcn4. Gcn4 directly binds promoter-regions and transcribes a subset of metabolic genes, particularly driving lysine and arginine biosynthesis. Gcn4 also globally represses lysine and arginine enriched transcripts, which include genes encoding the translation machinery. The Gcn4 dependent lysine and arginine supply thereby maintains the synthesis of the translation machinery. This is required to maintain translation capacity. Gcn4 consequently enables metabolic-precursor supply to bolster protein synthesis, and drive a growth program. Thus, we illustrate how growth and starvation outcomes are both controlled using the same Gcn4 transcriptional outputs that function in distinct contexts.

Dynamics of genetic variation in transcription factors and its implications for the evolution of regulatory networks in Bacteria.

The evolution of regulatory networks in Bacteria has largely been explained at macroevolutionary scales through lateral gene transfer and gene duplication. Transcription factors (TF) have been found to be less conserved across species than their target genes (TG). This would be expected if TFs accumulate mutations faster than TGs. This hypothesis is supported by several lab evolution studies which found TFs, especially global regulators, to be frequently mutated. Despite these studies, the contribution of point mutations in TFs to the evolution of regulatory network is poorly understood. We tested if TFs show greater genetic variation than their TGs using whole-genome sequencing data from a large collection of Escherichia coli isolates. TFs were less diverse than their TGs across natural isolates, with TFs of large regulons being more conserved. In contrast, TFs showed higher mutation frequency in adaptive laboratory evolution experiments. However, over long-term laboratory evolution spanning 60 000 generations, mutation frequency in TFs gradually declined after a rapid initial burst. Extrapolating the dynamics of genetic variation from long-term laboratory evolution to natural populations, we propose that point mutations, conferring large-scale gene expression changes, may drive the early stages of adaptation but gene regulation is subjected to stronger purifying selection post adaptation.

Early fate of exogenous promoters in E. coli.

Gene gain by horizontal gene transfer is a major pathway of genome innovation in bacteria. The current view posits that acquired genes initially need to be silenced and that a bacterial chromatin protein, H-NS, plays a role in this silencing. However, we lack direct observation of the early fate of a horizontally transferred gene to prove this theory. We combine sequencing, flow cytometry and sorting, followed by microscopy to monitor gene expression and its variability after large-scale random insertions of a reporter gene in a population of Escherichia coli bacteria. We find that inserted promoters have a wide range of gene-expression variability related to their location. We find that high-expression clones carry insertions that are not correlated with H-NS binding. Conversely, binding of H-NS correlates with silencing. Finally, while most promoters show a common level of extrinsic noise, some insertions show higher noise levels. Analysis of these high-noise clones supports a scenario of switching due to transcriptional interference from divergent ribosomal promoters. Altogether, our findings point to evolutionary pathways where newly-acquired genes are not necessarily silenced, but may immediately explore a wide range of expression levels to probe the optimal ones.

Evidence of non-random mutation rates suggests an evolutionary risk management strategy.

A central tenet in evolutionary theory is that mutations occur randomly with respect to their value to an organism; selection then governs whether they are fixed in a population. This principle has been challenged by long-standing theoretical models predicting that selection could modulate the rate of mutation itself. However, our understanding of how the mutation rate varies between different sites within a genome has been hindered by technical difficulties in measuring it. Here we present a study that overcomes previous limitations by combining phylogenetic and population genetic techniques. Upon comparing 34 Escherichia coli genomes, we observe that the neutral mutation rate varies by more than an order of magnitude across 2,659 genes, with mutational hot and cold spots spanning several kilobases. Importantly, the variation is not random: we detect a lower rate in highly expressed genes and in those undergoing stronger purifying selection. Our observations suggest that the mutation rate has been evolutionarily optimized to reduce the risk of deleterious mutations. Current knowledge of factors influencing the mutation rate­including transcription-coupled repair and context-dependent mutagenesis­do not explain these observations, indicating that additional mechanisms must be involved. The findings have important implications for our understanding of evolution and the control of mutations.

Using stochastic cell division and death to probe minimal units of cellular replication.

The invariant cell initiation mass measured in bacterial growth experiments has been interpreted as a minimal unit of cellular replication. Here we argue that the existence of such minimal units induces a coupling between the rates of stochastic cell division and death. To probe this coupling we tracked live and dead cells in Escherichia coli populations treated with a ribosome-targeting antibiotic. We find that the growth exponent from macroscopic cell growth or decay measurements can be represented as the difference of microscopic first-order cell division and death rates. The boundary between cell growth and decay, at which the number of live cells remains constant over time, occurs at the minimal inhibitory concentration (MIC) of the antibiotic. This state appears macroscopically static but is microscopically dynamic: division and death rates exactly cancel at MIC but each is remarkably high, reaching 60% of the antibiotic-free division rate. A stochastic model of cells as collections of minimal replicating units we term 'widgets' reproduces both steady-state and transient features of our experiments. Sub-cellular fluctuations of widget numbers stochastically drive each new daughter cell to one of two alternate fates, division or death. First-order division or death rates emerge as eigenvalues of a stationary Markov process, and can be expressed in terms of the widget's molecular properties. High division and death rates at MIC arise due to low mean and high relative fluctuations of widget number. Isolating cells at the threshold of irreversible death might allow molecular characterization of this minimal replication unit.

The genome-scale interplay amongst xenogene silencing, stress response and chromosome architecture in Escherichia coli.

The gene expression state of exponentially growing Escherichia coli cells is manifested by high expression of essential and growth-associated genes and low levels of stress-related and horizontally acquired genes. An important player in maintaining this homeostasis is the H-NS-StpA gene silencing system. A Δhns-stpA deletion mutant results in high expression of otherwise-silent horizontally acquired genes, many located in the terminus-half of the chromosome, and an indirect downregulation of many highly expressed genes. The Δhns-stpA double mutant displays slow growth. Using laboratory evolution we address the evolutionary strategies that E. coli would adopt to redress this gene expression imbalance. We show that two global gene regulatory mutations-(i) point mutations inactivating the stress-responsive sigma factor RpoS or σ38 and (ii) an amplification of ∼40% of the chromosome centred around the origin of replication-converge in partially reversing the global gene expression imbalance caused by Δhns-stpA. Transcriptome data of these mutants further show a three-way link amongst the global gene regulatory networks of H-NS and σ38, as well as chromosome architecture. Increasing gene expression around the terminus of replication results in a decrease in the expression of genes around the origin and vice versa; this appears to be a persistent phenomenon observed as an association across ∼300 publicly-available gene expression data sets for E. coli. These global suppressor effects are transient and rapidly give way to more specific mutations, whose roles in reversing the growth defect of H-NS mutations remain to be understood.

Identification of neoblast- and regeneration-specific miRNAs in the planarian Schmidtea mediterranea.

In recent years, the planarian Schmidtea mediterranea has emerged as a tractable model system to study stem cell biology and regeneration. MicroRNAs are small RNA species that control gene expression by modulating translational repression and mRNA stability and have been implicated in the regulation of various cellular processes. Though recent studies have identified several miRNAs in S. mediterranea, their expression in neoblast subpopulations and during regeneration has not been examined. Here, we identify several miRNAs whose expression is enriched in different neoblast subpopulations and in regenerating tissue at different time points in S. mediterranea. Some of these miRNAs were enriched within 3 h post-amputation and may, therefore, play a role in wound healing and/or neoblast migration. Our results also revealed miRNAs, such as sme-miR-2d-3p and the sme-miR-124 family, whose expression is enriched in the cephalic ganglia, are also expressed in the brain primordium during CNS regeneration. These results provide new insight into the potential biological functions of miRNAs in neoblasts and regeneration in planarians.

Context-dependent conservation of DNA methyltransferases in bacteria.

DNA methytransferases (MTs) in bacteria are best understood in the context of restriction-modification (R-M) systems, which act as bacterial immune systems against incoming DNA including phages, but have also been described as selfish elements. But several orphan MTs, which are not associated with any restriction enzyme, have also been characterized and may protect against parasitism by R-M systems. The occurrence of MTs in these two contexts, namely as part of R-M systems or as orphans, is poorly understood. Here we report the results of a comparative genomic survey of DNA MTs across ∼1000 bacterial genomes. We show that orphan MTs overwhelm R-M systems in their occurrence. In general, R-M MTs are poorly conserved, whereas orphans are nearly as conserved within a genus as any average gene. However, oligonucleotide usage and conservation patterns across genera suggest that both forms of MTs might have been horizontally acquired. We suggest that many orphan MTs might be 'degradation' products of R-M systems, based on the properties of orphan MTs encoded adjacent to highly diverged REs. In addition, several fully degraded R-M systems exist in which both the MT and the RE are highly divergent from their corresponding reference R-M pair. Despite their sporadic occurrence, conserved R-M systems are present in strength in two highly transformable genera, in which they may contribute to selection against integration of foreign DNA.

Genomic analysis of DNA binding and gene regulation by homologous nucleoid-associated proteins IHF and HU in Escherichia coli K12.

IHF and HU are two heterodimeric nucleoid-associated proteins (NAP) that belong to the same protein family but interact differently with the DNA. IHF is a sequence-specific DNA-binding protein that bends the DNA by over 160°. HU is the most conserved NAP, which binds non-specifically to duplex DNA with a particular preference for targeting nicked and bent DNA. Despite their importance, the in vivo interactions of the two proteins to the DNA remain to be described at a high resolution and on a genome-wide scale. Further, the effects of these proteins on gene expression on a global scale remain contentious. Finally, the contrast between the functions of the homo- and heterodimeric forms of proteins deserves the attention of further study. Here we present a genome-scale study of HU- and IHF binding to the Escherichia coli K12 chromosome using ChIP-seq. We also perform microarray analysis of gene expression in single- and double-deletion mutants of each protein to identify their regulons. The sequence-specific binding profile of IHF encompasses ∼30% of all operons, though the expression of <10% of these is affected by its deletion suggesting combinatorial control or a molecular backup. The binding profile for HU is reflective of relatively non-specific binding to the chromosome, however, with a preference for A/T-rich DNA. The HU regulon comprises highly conserved genes including those that are essential and possibly supercoiling sensitive. Finally, by performing ChIP-seq experiments, where possible, of each subunit of IHF and HU in the absence of the other subunit, we define genome-wide maps of DNA binding of the proteins in their hetero- and homodimeric forms.

Intracellular peroxynitrite perturbs redox balance, bioenergetics, and Fe-S cluster homeostasis in Mycobacterium tuberculosis.

The ability of Mycobacterium tuberculosis (Mtb) to tolerate nitric oxide (•NO) and superoxide (O2•-) produced by phagocytes contributes to its success as a human pathogen. Recombination of •NO and O2•- generates peroxynitrite (ONOO-), a potent oxidant produced inside activated macrophages causing lethality in diverse organisms. While the response of Mtb toward •NO and O2•- is well established, how Mtb responds to ONOO- remains unclear. Filling this knowledge gap is important to understand the persistence mechanisms of Mtb during infection. We synthesized a series of compounds that generate both •NO and O2•-, which should combine to produce ONOO-. From this library, we identified CJ067 that permeates Mtb to reliably enhance intracellular ONOO- levels. CJ067-exposed Mtb strains, including multidrug-resistant (MDR) and extensively drug-resistant (XDR) clinical isolates, exhibited dose-dependent, long-lasting oxidative stress and growth inhibition. In contrast, Mycobacterium smegmatis (Msm), a fast-growing, non-pathogenic mycobacterial species, maintained redox balance and growth in response to intracellular ONOO-. RNA-sequencing with Mtb revealed that CJ067 induces antioxidant machinery, sulphur metabolism, metal homeostasis, and a 4Fe-4S cluster repair pathway (suf operon). CJ067 impaired the activity of the 4Fe-4S cluster-containing TCA cycle enzyme, aconitase, and diminished bioenergetics of Mtb. Work with Mtb strains defective in SUF and IscS involved in Fe-S cluster biogenesis pathways showed that both systems cooperatively protect Mtb from intracellular ONOO- in vitro and inducible nitric oxide synthase (iNOS)-dependent growth inhibition during macrophage infection. Thus, Mtb is uniquely sensitive to intracellular ONOO- and targeting Fe-S cluster homeostasis is expected to promote iNOS-dependent host immunity against tuberculosis (TB).

Direct and indirect effects of H-NS and Fis on global gene expression control in Escherichia coli.

Nucleoid-associated proteins (NAPs) are global regulators of gene expression in Escherichia coli, which affect DNA conformation by bending, wrapping and bridging the DNA. Two of these--H-NS and Fis--bind to specific DNA sequences and structures. Because of their importance to global gene expression, the binding of these NAPs to the DNA was previously investigated on a genome-wide scale using ChIP-chip. However, variation in their binding profiles across the growth phase and the genome-scale nature of their impact on gene expression remain poorly understood. Here, we present a genome-scale investigation of H-NS and Fis binding to the E. coli chromosome using chromatin immunoprecipitation combined with high-throughput sequencing (ChIP-seq). By performing our experiments under multiple time-points during growth in rich media, we show that the binding regions of the two proteins are mutually exclusive under our experimental conditions. H-NS binds to significantly longer tracts of DNA than Fis, consistent with the linear spread of H-NS binding from high- to surrounding lower-affinity sites; the length of binding regions is associated with the degree of transcriptional repression imposed by H-NS. For Fis, a majority of binding events do not lead to differential expression of the proximal gene; however, it has a significant indirect effect on gene expression partly through its effects on the expression of other transcription factors. We propose that direct transcriptional regulation by Fis is associated with the interaction of tandem arrays of Fis molecules to the DNA and possible DNA bending, particularly at operon-upstream regions. Our study serves as a proof-of-principle for the use of ChIP-seq for global DNA-binding proteins in bacteria, which should become significantly more economical and feasible with the development of multiplexing techniques.

Cysteine desulfurase (IscS)-mediated fine-tuning of bioenergetics and SUF expression prevents Mycobacterium tuberculosis hypervirulence.

Iron-sulfur (Fe-S) biogenesis requires multiprotein assembly systems, SUF and ISC, in most prokaryotes. M. tuberculosis (Mtb) encodes a complete SUF system, the depletion of which was bactericidal. The ISC operon is truncated to a single gene iscS (cysteine desulfurase), whose function remains uncertain. Here, we show that MtbΔiscS is bioenergetically deficient and hypersensitive to oxidative stress, antibiotics, and hypoxia. MtbΔiscS resisted killing by nitric oxide (NO). RNA sequencing indicates that IscS is important for expressing regulons of DosR and Fe-S-containing transcription factors, WhiB3 and SufR. Unlike wild-type Mtb, MtbΔiscS could not enter a stable persistent state, continued replicating in mice, and showed hypervirulence. The suf operon was overexpressed in MtbΔiscS during infection in a NO-dependent manner. Suppressing suf expression in MtbΔiscS either by CRISPR interference or upon infection in inducible NO-deficient mice arrests hypervirulence. Together, Mtb redesigned the ISC system to "fine-tune" the expression of SUF machinery for establishing persistence without causing detrimental disease in the host.

Biosensor-integrated transposon mutagenesis reveals rv0158 as a coordinator of redox homeostasis in Mycobacterium tuberculosis.

Mycobacterium tuberculosis (Mtb) is evolutionarily equipped to resist exogenous reactive oxygen species (ROS) but shows vulnerability to an increase in endogenous ROS (eROS). Since eROS is an unavoidable consequence of aerobic metabolism, understanding how Mtb manages eROS levels is essential yet needs to be characterized. By combining the Mrx1-roGFP2 redox biosensor with transposon mutagenesis, we identified 368 genes (redoxosome) responsible for maintaining homeostatic levels of eROS in Mtb. Integrating redoxosome with a global network of transcriptional regulators revealed a hypothetical protein (Rv0158) as a critical node managing eROS in Mtb. Disruption of rv0158 (rv0158 KO) impaired growth, redox balance, respiration, and metabolism of Mtb on glucose but not on fatty acids. Importantly, rv0158 KO exhibited enhanced growth on propionate, and the Rv0158 protein directly binds to methylmalonyl-CoA, a key intermediate in propionate catabolism. Metabolite profiling, ChIP-Seq, and gene-expression analyses indicate that Rv0158 manages metabolic neutralization of propionate toxicity by regulating the methylcitrate cycle. Disruption of rv0158 enhanced the sensitivity of Mtb to oxidative stress, nitric oxide, and anti-TB drugs. Lastly, rv0158 KO showed poor survival in macrophages and persistence defect in mice. Our results suggest that Rv0158 is a metabolic integrator for carbon metabolism and redox balance in Mtb.

Draft genome sequence of Staphylococcus aureus ST672, an emerging disease clone from India.

We report the draft genome sequence of methicillin-resistant Staphylococcus aureus (MRSA) strain ST672, an emerging disease clone in India, from a septicemia patient. The genome size is about 2.82 Mb with 2,485 open reading frames (ORFs). The staphylococcal cassette chromosome mec (SCCmec) element (type V) and immune evasion cluster appear to be different from those of strain ST772 on preliminary examination.

Draft genome sequence of Staphylococcus aureus 118 (ST772), a major disease clone from India.

We report the draft genome sequence of an ST772 Staphylococcus aureus disease isolate carrying staphylococcal cassette chromosome mec (SCCmec) type V from a pyomyositis patient. Our de novo short read assembly is ∼2.8 Mb and encodes a unique Panton-Valentine leukocidin (PVL) phage with structural genes similar to those of ϕ7247PVL and novel lysogenic genes at the N termini.

Comparative genomics of cyclic-di-GMP signalling in bacteria: post-translational regulation and catalytic activity.

Cyclic-di-GMP is a bacterial second messenger that controls the switch between motile and sessile states. It is synthesized by proteins containing the enzymatic GGDEF domain and degraded by the EAL domain. Many bacterial genomes encode several copies of proteins containing these domains, raising questions on how the activities of parallel c-di-GMP signalling systems are segregated to avoid potentially deleterious cross-talk. Moreover, many 'hybrid' proteins contain both GGDEF and EAL domains; the relationship between the two apparently opposing enzymatic activities has been termed a 'biochemical conundrum'. Here, we present a computational analysis of 11 248 GGDEF- and EAL-containing proteins in 867 prokaryotic genomes to address these two outstanding questions. Over half of these proteins contain a signal for cell-surface localization, and a majority accommodate a signal-sensing partner domain; these indicate widespread prevalence of post-translational regulation that may segregate the activities of proteins that are co-expressed. By examining the conservation of amino acid residues in the GGDEF and EAL catalytic sites, we show that there are predominantly two types of hybrid proteins. In the first, both sites are intact; an additional regulatory partner domain, present in most of these proteins, might determine the balance between the two enzymatic activities. In the second type, only the EAL catalytic site is intact; these--unlike EAL-only proteins--generally contain a signal-sensing partner domain, suggesting distinct modes of regulation for EAL activity under different sequence contexts. Finally, we discuss the role of proteins that have lost GGDEF and EAL catalytic sites as potential c-di-GMP-binding effectors. Our findings will serve as a genomic framework for interpreting ongoing molecular investigations of these proteins.

An overview of prokaryotic transcription factors : a summary of function and occurrence in bacterial genomes.

Transcriptional initiation is arguably the most important control point for gene expression. It is regulated by a combination of factors, including DNA sequence and its three-dimensional topology, proteins and small molecules. In this chapter, we focus on the trans-acting factors of bacterial regulation. Initiation begins with the recruitment of the RNA polymerase holoenzyme to a specific locus upstream of the gene known as its promoter. The sigma factor, which is a component of the holoenzyme, provides the most fundamental mechanisms for orchestrating broad changes in gene expression state. It is responsible for promoter recognition as well as recruiting the holoenzyme to the promoter. Distinct sigma factors compete with for binding to a common pool of RNA polymerases, thus achieving condition-dependent differential expression. Another important class of bacterial regulators is transcription factors, which activate or repress transcription of target genes typically in response to an environmental or cellular trigger. These factors may be global or local depending on the number of genes and range of cellular functions that they target. The activities of both global and local transcription factors may be regulated either at a post-transcriptional level via signal-sensing protein domains or at the level of their own expression. In addition to modulating polymerase recruitment to promoters, several global factors are considered as "nucleoid-associated proteins" that impose structural constraints on the chromosome by altering the conformation of the bound DNA, thus influencing other processes involving DNA such as replication and recombination. This chapter concludes with a discussion of how regulatory interactions between transcription factors and their target genes can be represented as a network.

Replication-Dependent Organization Constrains Positioning of Long DNA Repeats in Bacterial Genomes.

Bacterial genome organization is primarily driven by chromosomal replication from a single origin of replication. However, chromosomal rearrangements, which can disrupt such organization, are inevitable in nature. Long DNA repeats are major players mediating rearrangements, large and small, via homologous recombination. Since changes to genome organization affect bacterial fitness-and more so in fast-growing than slow-growing bacteria-and are under selection, it is reasonable to expect that genomic positioning of long DNA repeats is also under selection. To test this, we identified identical DNA repeats of at least 100 base pairs across ∼6,000 bacterial genomes and compared their distribution in fast- and slow-growing bacteria. We found that long identical DNA repeats are distributed in a non-random manner across bacterial genomes. Their distribution differs in the overall number, orientation, and proximity to the origin of replication, between fast- and slow-growing bacteria. We show that their positioning-which might arise from a combination of the processes that produce repeats and selection on rearrangements that recombination between repeat elements might cause-permits less disruption to the replication-dependent genome organization of bacteria compared with random suggesting it as a major constraint to positioning of long DNA repeats.

Understanding the Genome-Wide Transcription Response To Various cAMP Levels in Bacteria Using Phenomenological Models.

Attempts to understand gene regulation by global transcription factors have largely been limited to expression studies under binary conditions of presence and absence of the transcription factor. Studies addressing genome-wide transcriptional responses to changing transcription factor concentration at high resolution are lacking. Here, we create a data set containing the entire Escherichia coli transcriptome in Luria-Bertani (LB) broth as it responds to 10 different cAMP concentrations spanning the biological range. We use the Hill's model to accurately summarize individual gene responses into three intuitively understandable parameters, Emax, n, and k, reflecting the sensitivity, nonlinearity, and midpoint of the dynamic range. Our data show that most cAMP-regulated genes have an n of >2, with their k values centered around the wild-type concentration of cAMP. Additionally, cAMP receptor protein (CRP) affinity to a promoter is correlated with Emax but not k, hinting that a high-affinity CRP promoter need not ensure transcriptional activation at lower cAMP concentrations and instead affects the magnitude of the response. Finally, genes belonging to different functional classes are tuned to have different k, n, and Emax values. We demonstrate that phenomenological models are a better alternative for studying gene expression trends than classical clustering methods, with the phenomenological constants providing greater insights into how genes are tuned in a regulatory network. IMPORTANCE Different genes may follow different trends in response to various transcription factor concentrations. In this study, we ask two questions: (i) what are the trends that different genes follow in response to changing transcription factor concentrations and (ii) what methods can be used to extract information from the gene trends so obtained. We demonstrate a method to analyze transcription factor concentration-dependent genome-wide expression data using phenomenological models. Conventional clustering methods and principal-component analysis (PCA) can be used to summarize trends in data but have limited interpretability. The use of phenomenological models greatly enhances the interpretability and thus utility of conventional clustering. Transformation of dose-response data into phenomenological constants opens up avenues to ask and answer many different kinds of question. We show that the phenomenological constants obtained from the model fits can be used to generate insights about network topology and allows integration of other experimental data such as chromatin immunoprecipitation sequencing (ChIP-seq) to understand the system in greater detail.

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

Restriction-modification (RM) systems are the most ubiquitous bacterial defence systems against bacteriophages. Using genome sequence data, we showed that RM systems are often shared among bacterial strains in a structured way. Examining the network of interconnections between bacterial strains within genera, we found that many strains share more RM systems than expected compared with a suitable null model. We also found that many genera have a larger than expected number of bacterial strains with unique RM systems. We used population dynamics models of closed and open phage-bacteria ecosystems to qualitatively understand the selection pressures that could lead to such network structures with enhanced overlap or uniqueness. In our models, we found that the phages impose a selection pressure that favours bacteria with greater number of RM systems, and higher overlap of RM systems with other strains, but in bacteria-dominated states, this is opposed by the increased cost-to-growth rate of these bacteria. Similar to what we observed in the genome data, we found that two distinct bacterial strategies emerge - strains either have a greater overlap than expected, or, at the other extreme, have unique RM systems. 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.