Phase variation in the phage growth limitation system of Streptomyces coelicolor A3(2)

Identification and characterization of the WYL BrxR protein and its gene as separable regulatory elements of a BREX phage restriction system

Bacteriophage exclusion ('BREX') phage restriction systems are found in a wide range of bacteria. Various BREX systems encode unique combinations of proteins that usually include a site-specific methyltransferase; none appear to contain a nuclease. Here we describe the identification and characterization of a Type I BREX system from Acinetobacter and the effect of deleting each BREX ORF on growth, methylation, and restriction. We identified a previously uncharacterized gene in the BREX operon that is dispensable for methylation but involved in restriction. Biochemical and crystallographic analyses of this factor, which we term BrxR ('BREX Regulator'), demonstrate that it forms a homodimer and specifically binds a DNA target site upstream of its transcription start site. Deletion of the BrxR gene causes cell toxicity, reduces restriction, and significantly increases the expression of BrxC. In contrast, the introduction of a premature stop codon into the BrxR gene, or a point mutation blocking its DNA binding ability, has little effect on restriction, implying that the BrxR coding sequence and BrxR protein play independent functional roles. We speculate that elements within the BrxR coding sequence are involved in cis regulation of anti-phage activity, while the BrxR protein itself plays an additional regulatory role, perhaps during horizontal transfer.

Reassembling a cannon in the DNA defense arsenal: Genetics of StySA, a BREX phage exclusion system in Salmonella lab strains

Understanding mechanisms that shape horizontal exchange in prokaryotes is a key problem in biology. A major limit on DNA entry is imposed by restriction-modification (RM) processes that depend on the pattern of DNA modification at host-specified sites. In classical RM, endonucleolytic DNA cleavage follows detection of unprotected sites on entering DNA. Recent investigation has uncovered BREX (BacteRiophage EXclusion) systems. These RM-like activities employ host protection by DNA modification, but immediate replication arrest occurs without evident of nuclease action on unmodified phage DNA. Here we show that the historical stySA RM locus of Salmonella enterica sv Typhimurium is a variant BREX system. A laboratory strain disabled for both the restriction and methylation activity of StySA nevertheless has wild type sequence in pglX, the modification gene homolog. Instead, flanking genes pglZ and brxC each carry multiple mutations (μ) in their C-terminal domains. We further investigate this system in situ, replacing the mutated pglZμ and brxCμ genes with the WT counterpart. PglZ-WT supports methylation in the presence of either BrxCμ or BrxC-WT but not in the presence of a deletion/insertion allele, ΔbrxC::cat. Restriction requires both BrxC-WT and PglZ-WT, implicating the BrxC C-terminus specifically in restriction activity. These results suggests that while BrxC, PglZ and PglX are principal components of the BREX modification activity, BrxL is required for restriction only. Furthermore, we show that a partial disruption of brxL disrupts transcription globally.

Horizontal gene transfer is a major driver of evolution and adaptation in bacteria. Genes from outside may be beneficial or dangerous to the receiving cell. Benefits include new food sources such as sugars, or new homes by adhesion, or new resistances, as to antibiotics. Dangers are posed by bacteriophages--viruses that take over the cell machinery, multiply, and release progeny to kill sister cells. Host-dependent restriction-modification systems enable defense that distinguishes relatives from strangers: using a modification pattern (M) carried by DNA bases added by the host cell to prevent restriction (R). Sisters and cousin cells will have the same protective pattern on DNA, while DNA of foreign origin will have the wrong M pattern and be restricted (R, rejected). Typically, restriction involves nuclease digestion. Here we address the enigmatic StySA RM system, one of the earliest to be genetically characterized. It is a variant of the newly recognized defense mechanism, BREX. BREX systems also track DNA history via modification pattern, but restrict by a novel, uncharacterized mechanism. Like other BREX family systems, StySA-BREX modification requires multiple components. When StySA-BREX transcription is unbalanced, we find global disruption of gene transcription. The disruption pattern does not suggest SOS-inducing damage to DNA.

Phage-Resistant Phase-Variant Sub-populations Mediate Herd Immunity Against Bacteriophage Invasion of Bacterial Meta-Populations

Hypermutable loci are widespread in bacteria as mechanisms for rapid generation of phenotypic diversity within a population that enables survival of fluctuating, often antagonistic, selection pressures. Localized hypermutation can mediate phase variation and enable survival of bacteriophage predation due to high frequency, reversible alterations in the expression of phage receptors. As phase variation can also generate population-to-population heterogeneity, we hypothesized that this phenomenon may facilitate survival of spatially-separated bacterial populations from phage invasion in a manner analogous to herd immunity to infectious diseases in human populations. The lic2A gene of Haemophilus influenzae is subject to “ON” and “OFF” switches in expression mediated by mutations in a 5′CAAT repeat tract present within the reading frame. The enzyme encoded by lic2A mediates addition of a galactose moiety of the lipopolysaccharide. This moiety is required for attachment of the HP1C1 phage such that the ON state of the lic2A gene is associated with HP1c1 susceptibility while the OFF state is resistant to infection. We developed an “oscillating prey assay” to examine phage spread through a series of sub-populations of Haemophilus influenzae whose phage receptor is in an ON or OFF state. Phage extinction was frequently observed when the proportion of phage-resistant sub-populations exceeded 34%. In silico modeling indicated that phage extinction was interdependent on phage loss during transfer between sub-populations and the frequency of resistant sub-populations. In a fixed-area oscillating prey assay, heterogeneity in phage resistance was observed to generate vast differences in phage densities across a meta-population of multiple bacterial sub-populations resulting in protective quarantining of some sub-populations from phage attack. We conclude that phase-variable hypermutable loci produce bacterial “herd immunity” with resistant intermediary-populations acting as a barricade to reduce the viral load faced by phage-susceptible sub-populations. This paradigm of meta-population protection is applicable to evolution of hypermutable loci in multiple bacteria-phage and host-pathogen interactions.

The phage growth limitation system in Streptomyces coelicolor A(3)2 is a toxin/antitoxin system, comprising enzymes with DNA methyltransferase, protein kinase and ATPase activity

The phage growth limitation system of Streptomyces coelicolor A3(2) is an unusual bacteriophage defence mechanism. Progeny ϕC31 phage from an initial infection are thought to be modified such that subsequent infections are attenuated in a Pgl(+) host but normal in a Pgl(-) strain. Earlier work identified four genes required for phage resistance by Pgl. Here we demonstrate that Pgl is an elaborate and novel phage restriction system that, in part, comprises a toxin/antitoxin system where PglX, a DNA methyltransferase is toxic in the absence of a functional PglZ. In addition, the ATPase activity of PglY and a protein kinase activity in PglW are shown to be essential for phage resistance by Pgl. We conclude that on infection of a Pgl(+) cell by bacteriophage ϕC31, PglW transduces a signal, probably via phosphorylation, to other Pgl proteins resulting in the activation of the DNA methyltransferase, PglX and this leads to phage restriction.

BREX is a novel phage resistance system widespread in microbial genomes

The perpetual arms race between bacteria and phage has resulted in the evolution of efficient resistance systems that protect bacteria from phage infection. Such systems, which include the CRISPR-Cas and restriction-modification systems, have proven to be invaluable in the biotechnology and dairy industries. Here, we report on a six-gene cassette in Bacillus cereus which, when integrated into the Bacillus subtilis genome, confers resistance to a broad range of phages, including both virulent and temperate ones. This cassette includes a putative Lon-like protease, an alkaline phosphatase domain protein, a putative RNA-binding protein, a DNA methylase, an ATPase-domain protein, and a protein of unknown function. We denote this novel defense system BREX (Bacteriophage Exclusion) and show that it allows phage adsorption but blocks phage DNA replication. Furthermore, our results suggest that methylation on non-palindromic TAGGAG motifs in the bacterial genome guides self/non-self discrimination and is essential for the defensive function of the BREX system. However, unlike restriction-modification systems, phage DNA does not appear to be cleaved or degraded by BREX, suggesting a novel mechanism of defense. Pan genomic analysis revealed that BREX and BREX-like systems, including the distantly related Pgl system described in Streptomyces coelicolor, are widely distributed in ~10% of all sequenced microbial genomes and can be divided into six coherent subtypes in which the gene composition and order is conserved. Finally, we detected a phage family that evades the BREX defense, implying that anti-BREX mechanisms may have evolved in some phages as part of their arms race with bacteria.

The origin of the bacterial immune response

Bacteriophages are probably the oldest viruses, having appeared early during bacterial evolution. Therefore, bacteria and bacteriophages have a long history of co-evolution in which bacteria have developed multiple resistance mechanisms against bacteriophages. These mechanisms, that are very diverse and are in constant evolution, allow the survival of the bacteria. Bacteriophages have adapted to bacterial defense systems, devised strategies to evade these anti-phage mechanisms and restored their infective capacity. In this chapter, we review the bacterial strategies that hinder the phage infection as well as the counter-defense mechanisms developed.

A new group of cosmopolitan bacteriophages induce a carrier state in the pandemic strain of Vibrio parahaemolyticus

A clonal population of pathogenic Vibrio parahaemolyticus O3 : K6 serovar has spread in coastal waters, causing outbreaks worldwide since 1996. Bacteriophage infection is one of the main factors affecting bacterial strain concentration in the ocean. We studied the occurrence and properties of phages infecting this V. parahaemolyticus pandemic strain in coastal waters. Analysing 143 samples, phages were found in 13. All isolates clustered in a closely related group of podophages with at least 90% nucleotide sequence identity in three essential genes, despite distant geographical origins. These bacteriophages were able to multiply on the V. parahaemolyticus pandemic strain, but the impact on host concentration and subsequent growth was negligible. Infected bacteria continued producing the phage but were not lysogenized. The phage genome of prototype strain VP93 is 43 931 nucleotides and contains 337 bp direct terminal repeats at both ends. VP93 is the first non-Pseudomonas phage related to the PhiKMV-like subgroup of the T7 supergroup. The lack of a major effect on host growth suggests that these phages exert little control on the propagation of the pandemic strain in the environment. This form of phage growth can be modelled if phage-sensitive and -resistant cells that convert to each other with a high frequency are present in clonal cultures of pandemic V. parahaemolyticus.

Quantitative models of in vitro bacteriophage-host dynamics and their application to phage therapy

Phage therapy is the use of bacteriophages as antimicrobial agents for the control of pathogenic and other problem bacteria. It has previously been argued that successful application of phage therapy requires a good understanding of the non-linear kinetics of phage–bacteria interactions. Here we combine experimental and modelling approaches to make a detailed examination of such kinetics for the important food-borne pathogen Campylobacter jejuni and a suitable virulent phage in an in vitro system. Phage-insensitive populations of C. jejuni arise readily, and as far as we are aware this is the first phage therapy study to test, against in vitro data, models for phage–bacteria interactions incorporating phage-insensitive or resistant bacteria. We find that even an apparently simplistic model fits the data surprisingly well, and we confirm that the so-called inundation and proliferation thresholds are likely to be of considerable practical importance to phage therapy. We fit the model to time series data in order to estimate thresholds and rate constants directly. A comparison of the fit for each culture reveals density-dependent features of phage infectivity that are worthy of further investigation. Our results illustrate how insight from empirical studies can be greatly enhanced by the use of kinetic models: such combined studies of in vitro systems are likely to be an essential precursor to building a meaningful picture of the kinetic properties of in vivo phage therapy.

Phage therapy is an antimicrobial treatment based on specific viruses which are natural predators of bacteria. This approach is being promoted as a possible alternative treatment for use against antibiotic-resistant strains of bacteria. Despite its long history and many potential benefits, adoption of phage therapy has been retarded by a variety of factors, including a poor understanding of the therapeutic consequences of the phage–bacteria relationship. In our work we bring together theory and data by testing kinetic models of phage–bacteria interactions against data for an important agent of human food poisoning, Campylobacter jejuni. Our model explicitly allows for resistant bacteria because these have not been properly accounted for in previous phage therapy theory but will be relevant to practical applications. The excellent fit of our model to the data confirms the value of such combined approaches and supports an interpretative viewpoint based on critical density-dependent thresholds that are not part of standard pharmacology. We also find that phage activity appears to be dose-dependent, and we speculate on possible causes for this. Our work illustrates how mathematical models can considerably enhance insights from empirical studies, as an important step in advancing the understanding of phage therapy.

The evolution of bet-hedging adaptations to rare scenarios

When faced with a variable environment, organisms may switch between different strategies according to some probabilistic rule. In an infinite population, evolution is expected to favor the rule that maximizes geometric mean fitness. If some environments are encountered only rarely, selection may not be strong enough for optimal switching probabilities to evolve. Here we calculate the evolution of switching probabilities in a finite population by analyzing fixation probabilities of alleles specifying switching rules. We calculate the conditions required for the evolution of phenotypic switching as a form of bet-hedging as a function of the population size N, the rate theta at which a rare environment is encountered, and the selective advantage s associated with switching in the rare environment. We consider a simplified model in which environmental switching and phenotypic switching are one-way processes, and mutation is symmetric and rare with respect to the timescale of fixation events. In this case, the approximate requirements for bet-hedging to be favored by a ratio of at least R are that sN>log(R) and thetaN>square root R .

Hypervariation and phase variation in the bacteriophage 'resistome'

Most bacteria encode proteins for defence against infection by bacteriophages. The mechanisms that bring about phage defence are extremely diverse, suggesting frequent independent evolution of novel processes. Phage defence determinants are often plasmid or phage-encoded and many that are chromosomal show evidence of lateral transfer. Recent studies on restriction-modification (R-M) systems show that these genes are amongst the most rapidly evolving. Some bacteria have contingency genes that encode alternative target specificity determinants for Type I or Type III R-M systems, thus expanding the range of phages against which the host population is immune. The most counter-intuitive observation, however, is the prevalence of phase variation in many restriction systems, but recent arguments suggest that switching off expression of R-M systems can aid phage defence.

Comparative genomic hybridizations reveal absence of large Streptomyces coelicolor genomic islands in Streptomyces lividans

Background

The genomes of Streptomyces coelicolor and Streptomyces lividans bear a considerable degree of synteny. While S. coelicolor is the model streptomycete for studying antibiotic synthesis and differentiation, S. lividans is almost exclusively considered as the preferred host, among actinomycetes, for cloning and expression of exogenous DNA. We used whole genome microarrays as a comparative genomics tool for identifying the subtle differences between these two chromosomes.

Results

We identified five large S. coelicolor genomic islands (larger than 25 kb) and 18 smaller islets absent in S. lividans chromosome. Many of these regions show anomalous GC bias and codon usage patterns. Six of them are in close vicinity of tRNA genes while nine are flanked with near perfect repeat sequences indicating that these are probable recent evolutionary acquisitions into S. coelicolor. Embedded within these segments are at least four DNA methylases and two probable methyl-sensing restriction endonucleases. Comparison with S. coelicolor transcriptome and proteome data revealed that some of the missing genes are active during the course of growth and differentiation in S. coelicolor. In particular, a pair of methylmalonyl CoA mutase (mcm) genes involved in polyketide precursor biosynthesis, an acyl-CoA dehydrogenase implicated in timing of actinorhodin synthesis and bldB, a developmentally significant regulator whose mutation causes complete abrogation of antibiotic synthesis belong to this category.

Conclusion

Our findings provide tangible hints for elucidating the genetic basis of important phenotypic differences between these two streptomycetes. Importantly, absence of certain genes in S. lividans identified here could potentially explain the relative ease of DNA transformations and the conditional lack of actinorhodin synthesis in S. lividans.

The loss of adaptive plasticity during long periods of environmental stasis

Adaptive plasticity allows populations to adjust rapidly to environmental change. If this is useful only rarely, plasticity may undergo mutational degradation and be lost from a population. We consider a population of constant size N undergoing loss of plasticity at functional mutation rate m and with selective advantage s associated with loss. Environmental change events occur at rate theta per generation, killing all individuals that lack plasticity. The expected time until loss of plasticity in a fluctuating environment is always at least tau, the expected time until loss of plasticity in a static environment. When mN > 1 and N theta >> 1, we find that plasticity will be maintained for an average of at least 10(8) generations in a single population, provided tau > 18/theta. In a metapopulation, plasticity is retained under the more lenient condition tau > 1.3/theta, irrespective of mN, for a modest number of demes. We calculate both exact and approximate solutions for tau and find that it is linearly dependent only on the logarithm of N, and so, surprisingly, both the population size and the number of demes in the metapopulation make little difference to the retention of plasticity. Instead, tau is dominated by the term 1/(m+s/2).

A microbial modified prisoner's dilemma game: how frequency-dependent selection can lead to random phase variation

Random phase variation (RPV) is a control strategy in which the expression of a cell state or phenotype randomly alternates between discrete 'on' and 'off' states. Though this mode of control is common for bacterial virulence factors like pili and toxins, precise conditions under which RPV confers an advantage have not been well defined. In Part I of this study, we predicted that fluctuating environments select for RPV if transitions between different selective environments cannot be reliably sensed (J. Theor. Biol. (2005)). However, selective forces both inside and outside of human hosts are also likely to be frequency dependent in the sense that the fitnesses of some bacterial states are greatest when rare. Here we show that RPV at slow rates can provide a survival advantage in such a frequency-dependent environment by generating population heterogeneity, essentially mimicking a polymorphism. More surprisingly, RPV at a faster 'optimal' rate can shift the population composition toward an optimal growth rate that exceeds that possible for polymorphic populations, but this optimal strategy is not evolutionarily stable. The population would be most fit if all cells randomly phase varied at the optimal rate, but individual cells have a growth-rate incentive to defect (mutate) to other switching rates or non-phase variable phenotype expression, leading to an overall loss of fitness of the individual and the population. This scenario describes a modified Prisoner's Dilemma game (Evolution and the Theory of Games, Cambridge University Press, Cambridge, New York, 1982, viii, 224pp.; Nature 398 (6726) (1999) 367), with random phase variation at optimal switching rates serving as the cooperation strategy.

Diversity in times of adversity: probabilistic strategies in microbial survival games

Population diversification strategies are ubiquitous among microbes, encompassing random phase-variation (RPV) of pathogenic bacteria, viral latency as observed in some bacteriophage and HIV, and the non-genetic diversity of bacterial stress responses. Precise conditions under which these diversification strategies confer an advantage have not been well defined. We develop a model of population growth conditioned on dynamical environmental and cellular states. Transitions among cellular states, in turn, may be biased by possibly noisy readings of the environment from cellular sensors. For various types of environmental dynamics and cellular sensor capability, we apply game-theoretic analysis to derive the evolutionarily stable strategy (ESS) for an organism and determine when that strategy is diversification. We find that: (1) RPV, effecting a sort of Parrondo paradox wherein random alternations between losing strategies produce a winning strategy, is selected when transitions between different selective environments cannot be sensed, (2) optimal RPV cell switching rates are a function of environmental lifecycle asymmetries and environmental autocorrelation, (3) probabilistic diversification upon entering a new environment is selected when sensors can detect environmental transitions but have poor precision in identifying new environments, and (4) in the presence of excess additive noise, low-pass filtering is required for evolutionary stability. We show that even when RPV is not the ESS, it may minimize growth rate variance and the risk of extinction due to 'unlucky' environmental dynamics.

Phase and antigenic variation in bacteria

Phase and antigenic variation result in a heterogenic phenotype of a clonal bacterial population, in which individual cells either express the phase-variable protein(s) or not, or express one of multiple antigenic forms of the protein, respectively. This form of regulation has been identified mainly, but by no means exclusively, for a wide variety of surface structures in animal pathogens and is implicated as a virulence strategy. This review provides an overview of the many bacterial proteins and structures that are under the control of phase or antigenic variation. The context is mainly within the role of the proteins and variation for pathogenesis, which reflects the main body of literature. The occurrence of phase variation in expression of genes not readily recognizable as virulence factors is highlighted as well, to illustrate that our current knowledge is incomplete. From recent genome sequence analysis, it has become clear that phase variation may be more widespread than is currently recognized, and a brief discussion is included to show how genome sequence analysis can provide novel information, as well as its limitations. The current state of knowledge of the molecular mechanisms leading to phase variation and antigenic variation are reviewed, and the way in which these mechanisms form part of the general regulatory network of the cell is addressed. Arguments both for and against a role of phase and antigenic variation in immune evasion are presented and put into new perspective by distinguishing between a role in bacterial persistence in a host and a role in facilitating evasion of cross-immunity. Finally, examples are presented to illustrate that phase-variable gene expression should be taken into account in the development of diagnostic assays and in the interpretation of experimental results and epidemiological studies.