The tripartite immunity system of phages P1 and P7

Insight into the Mechanism of Lysogeny Control of phiCDKH01 Bacteriophage Infecting Clinical Isolate of Clostridioides difficile

Clostridioides difficile is a causative agent of antibiotic-associated diarrhea as well as pseudomembranous colitis. So far, all known bacteriophages infecting these bacteria are temperate, which means that instead of prompt lysis of host cells, they can integrate into the host genome or replicate episomally. While C. difficile phages are capable of spontaneous induction and entering the lytic pathway, very little is known about the regulation of their maintenance in the state of lysogeny. In this study, we investigated the properties of a putative major repressor of the recently characterized C. difficile phiCDKH01 bacteriophage. A candidate protein belongs to the XRE family and controls the transcription of genes encoding putative phage antirepressors, known to be involved in the regulation of lytic development. Hence, the putative major phage repressor is likely to be responsible for maintenance of the lysogeny.

Isolation and characterization of novel Staphylococcus aureus bacteriophage Hesat from dairy origin

A novel temperate phage, named Hesat, was isolated by the incubation of a dairy strain of Staphylococcus aureus belonging to spa-type t127 with either bovine or ovine milk. Hesat represents a new species of temperate phage within the Phietavirus genus of the Azeredovirinae subfamily. Its genome has a length of 43,129 bp and a GC content of 35.11% and contains 75 predicted ORFs, some of which linked to virulence. This includes (i) a pathogenicity island (SaPln2), homologous to the type II toxin-antitoxin system PemK/MazF family toxin; (ii) a DUF3113 protein (gp30) that is putatively involved in the derepression of the global repressor Stl; and (iii) a cluster coding for a PVL. Genomic analysis of the host strain indicates Hesat is a resident prophage. Interestingly, its induction was obtained by exposing the bacterium to milk, while the conventional mitomycin C-based approach failed. The host range of phage Hesat appears to be broad, as it was able to lyse 24 out of 30 tested S. aureus isolates. Furthermore, when tested at high titer (108 PFU/ml), Hesat phage was also able to lyse a Staphylococcus muscae isolate, a coagulase-negative staphylococcal strain. KEY POINTS: • A new phage species was isolated from a Staphylococcus aureus bovine strain. • Pathogenicity island and PVL genes are encoded within phage genome. • The phage is active against most of S. aureus strains from both animal and human origins.

Identification and genomic analysis of temperate Halomonas bacteriophage vB_HmeY_H4907 from the surface sediment of the Mariana Trench at a depth of 8,900 m

Viruses play crucial roles in the ecosystem by modulating the host community structure, mediating biogeochemical cycles, and compensating for the metabolism of host cells. Mariana Trench, the world's deepest hadal habitat, harbors a variety of unique microorganisms that have adapted to its extreme conditions of low temperatures, high pressure, and nutrient scarcity. However, our knowledge about isolated hadal phage strains in the hadal trench is still limited. This study reported the discovery of a temperate phage, vB_HmeY_H4907, infecting Halomonas meridiana H4907, isolated from surface sediment from the Mariana Trench at a depth of 8,900 m. To our best knowledge, it is the deepest isolated siphovirus from the ocean. Its 40,452 bp linear dsDNA genome has 57.64% GC content and 55 open reading frames, and it is highly homologous to its host. Phylogenetic analysis and average nucleotide sequence identification reveal that vB_HmeY_H4907 is separated from the isolated phages and represents a new family, Suviridae, with eight predicted proviruses and six uncultured viral genomes. They are widely distributed in the ocean, suggesting a prevalence of this viral family in the deep sea. These findings expand our understanding of the phylogenetic diversity and genomic features of hadal lysogenic phages, provide essential information for further studies of phage-host interactions and evolution, and may reveal new insights into the lysogenic lifestyles of viruses inhabiting the hadal ocean. IMPORTANCE Halomonas phage vB_HmeY_H4907 is the deepest isolated siphovirus from the ocean, and it represents a novel abundant viral family in the ocean. This study provides insights into the genomic, phylogenetic, and ecological characteristics of the new viral family, namely, Suviridae.

An adenine/thymidine-rich region is integral to RepL-mediated DNA replication

The lytic replication of bacteriophage P1 requires RepL expression and the lytic stage origin, oriL, which is postulated to be located within repL gene sequence. The exact sequence of P1 oriL and the mechanism(s) of RepL-mediated DNA replication, however, are not fully understood. By using repL gene expression to induce DNA replication of a gfp and a rfp reporter plasmids, we demonstrated that synonymous base substitution in an adenine/thymidine-rich region of repL gene sequence, termed AT2, significantly inhibited the RepL-mediated signal amplification. Contrastingly, mutations in an IHF and two DnaA binding sites did not affect the RepL-mediated signal amplification significantly. A truncated repL sequence with the AT2 region allowed RepL-mediated signal amplification in trans therefore verifying a significant role of the AT2 region in RepL-mediated DNA replication. A combination of repL gene expression and a non-protein-coding copy of repL gene sequence (termed nc-repL) was able to amplify the output of an arsenic biosensor. Furthermore, mutation(s) at single or multiple positions within the AT2 region produced varying levels of RepL-mediated signal amplification. Overall, our results provide novel insights into the identity and location of P1 oriL as well as demonstrating the potential of using repL constructs to amplify and modulate the output of genetic biosensors.

Interactions between Viral Regulatory Proteins Ensure an MOI-Independent Probability of Lysogeny during Infection by Bacteriophage P1

Phage P1 is a temperate phage which makes the lytic or lysogenic decision upon infecting bacteria. During the lytic cycle, progeny phages are produced and the cell lyses, and in the lysogenic cycle, P1 DNA exists as a low-copy-number plasmid and replicates autonomously. Previous studies at the bulk level showed that P1 lysogenization was independent of multiplicity of infection (MOI; the number of phages infecting a cell), whereas lysogenization probability of the paradigmatic phage λ increases with MOI. However, the mechanism underlying the P1 behavior is unclear. In this work, using a fluorescent reporter system, we demonstrated this P1 MOI-independent lysogenic response at the single-cell level. We further observed that the activity of the major repressor of lytic functions (C1) is a determining factor for the final cell fate. Specifically, the repression activity of P1, which arises from a combination of C1, the anti-repressor Coi, and the corepressor Lxc, remains constant for different MOI, which results in the MOI-independent lysogenic response. Additionally, by increasing the distance between phages that infect a single cell, we were able to engineer a λ-like, MOI-dependent lysogenization upon P1 infection. This suggests that the large separation of coinfecting phages attenuates the effective communication between them, allowing them to make decisions independently of each other. Our work establishes a highly quantitative framework to describe P1 lysogeny establishment. This system plays an important role in disseminating antibiotic resistance by P1-like plasmids and provides an alternative to the lifestyle of phage λ.

Virus-Host Interaction Gets Curiouser and Curiouser. PART I: Phage P1vir Enhanced Development in an E. coli DksA-Deficient Cell

Bacteriophage P1 is among the best described bacterial viruses used in molecular biology. Here, we report that deficiency in the host cell DksA protein, an E. coli global transcription regulator, improves P1 lytic development. Using genetic and microbiological approaches, we investigated several aspects of P1vir biology in an attempt to understand the basis of this phenomenon. We found several minor improvements in phage development in the dksA mutant host, including more efficient adsorption to bacterial cell and phage DNA replication. In addition, gene expression of the main repressor of lysogeny C1, the late promoter activator Lpa, and lysozyme are downregulated in the dksA mutant. We also found nucleotide substitutions located in the phage immunity region immI, which may be responsible for permanent virulence of phage P1vir. We suggest that downregulation of C1 may lead to a less effective repression of lysogeny maintaining genes and that P1vir may be balancing between lysis and lysogeny, although finally it is able to enter the lytic pathway only. The mentioned improvements, such as more efficient replication and more “gentle” cell lysis, while considered minor individually, together may account for the phenomenon of a more efficient P1 phage development in a DksA-deficient host.

Repression of the lysogenic PR promoter in bacteriophage TP901-1 through binding of a CI-MOR complex to a composite OM-OR operator

A functional genetic switch from the lactococcal bacteriophage TP901-1, deciding which of two divergently transcribing promoters becomes most active and allows this bi-stable decision to be inherited in future generations requires a DNA region of less than 1 kb. The fragment encodes two repressors, CI and MOR, transcribed from the P_R and P_L promoters respectively. CI can repress the transcription of the mor gene at three operator sites (O_R, O_L, and O_D), leading to the immune state. Repression of the cI gene, leading to the lytic (anti-immune) state, requires interaction between CI and MOR by an unknown mechanism, but involving a CI:MOR complex. A consensus for putative MOR binding sites (O_M sites), and a common topology of three O_M sites adjacent to the O_R motif was here identified in diverse phage switches that encode CI and MOR homologs, in a search for DNA sequences similar to the TP901-1 switch. The O_R site and all putative O_M sites are important for establishment of the anti-immune repression of P_R, and a putative DNA binding motif in MOR is needed for establishment of the anti-immune state. Direct evidence for binding between CI and MOR is here shown by pull-down experiments, chemical crosslinking, and size exclusion chromatography. The results are consistent with two possible models for establishment of the anti-immune repression of cI expression at the P_R promoter.

Diversity of P1 phage-like elements in multidrug resistant Escherichia coli

The spread of multidrug resistance via mobile genetic elements is a major clinical and veterinary concern. Pathogenic Escherichia coli harbour antibiotic resistance and virulence genes mainly on plasmids, but also bacteriophages and hybrid phage-like plasmids. In this study, the genomes of three E. coli phage-like plasmids, pJIE250-3 from a human E. coli clinical isolate, pSvP1 from a porcine ETEC O157 isolate, and pTZ20_1P from a porcine commensal E. coli, were sequenced (PacBio RSII), annotated and compared. All three elements are coliphage P1 variants, each with unique adaptations. pJIE250-3 is a P1-derivative that has lost lytic functions and contains no accessory genes. In pTZ20_1P and pSvP1, a core P1-like genome is associated with insertion sequence-mediated acquisition of plasmid modules encoding multidrug resistance and virulence, respectively. The transfer ability of pTZ20_1P, carrying antibiotic resistance markers, was also tested and, although this element was not able to transfer by conjugation, it was able to lysogenize a commensal E. coli strain with consequent transfer of resistance. The incidence of P1-like plasmids (~7%) in our E. coli collections correlated well with that in public databases. This study highlights the need to investigate the contribution of phage-like plasmids to the successful spread of antibiotic resistant pathotypes.

Evolution of Superinfection Immunity in Cluster A Mycobacteriophages

Many aspects regarding superinfection, immunity, virulence, and the evolution of immune specificities are poorly understood due to the lack of large collections of isolated and sequenced phages with a spectrum of genetic diversity. Using a genetically diverse collection of Cluster A phages, we show that the classical and relatively straightforward patterns of homoimmunity, heteroimmunity, and virulence result from interactions between homotypic and heterotypic phages at the extreme edges of an evolutionary continuum of immune specificities. Genetic interactions between mesotypic phages result in more complex mesoimmunity phenotypes and virulence profiles. These results highlight that the evolution of immune specificities can be shaped by homotypic and mesotypic interactions and may be more dynamic than previously considered.

Temperate phages encode an immunity system to control lytic gene expression during lysogeny. This gene regulatory circuit consists of multiple interacting genetic elements, and although it is essential for controlling phage growth, it is subject to conflicting evolutionary pressures. During superinfection of a lysogen, the prophage’s circuit interacts with the superinfecting phage’s circuit and prevents lytic growth if the two circuits are closely related. The circuitry is advantageous since it provides the prophage with a defense mechanism, but the circuitry is also disadvantageous since it limits the phage’s host range during superinfection. Evolutionarily related phages have divergent, orthogonal immunity systems that no longer interact and are heteroimmune, but we do not understand how immunity systems evolve new specificities. Here, we use a group of Cluster A mycobacteriophages that exhibit a spectrum of genetic diversity to examine how immunity system evolution impacts superinfection immunity. We show that phages with mesotypic (i.e., genetically related but distinct) immunity systems exhibit asymmetric and incomplete superinfection phenotypes. They form complex immunity networks instead of well-defined immunity groups, and mutations conferring escape (i.e., virulence) from homotypic or mesotypic immunity have various escape specificities. Thus, virulence and the evolution of new immune specificities are shaped by interactions with homotypic and mesotypic immunity systems.

1st German Phage Symposium-Conference Report

In Germany, phage research and application can be traced back to the beginning of the 20th century. However, with the triumphal march of antibiotics around the world, the significance of bacteriophages faded in most countries, and respective research mainly focused on fundamental questions and niche applications. After a century, we pay tribute to the overuse of antibiotics that led to multidrug resistance and calls for new strategies to combat pathogenic microbes. Against this background, bacteriophages came into the spotlight of researchers and practitioners again resulting in a fast growing “phage community”. In October 2017, part of this community met at the 1st German Phage Symposium to share their knowledge and experiences. The participants discussed open questions and challenges related to phage therapy and the application of phages in general. This report summarizes the presentations given, highlights the main points of the round table discussion and concludes with an outlook for the different aspects of phage application.

The genome sequence of Escherichia coli tailed phage D6 and the diversity of Enterobacteriales circular plasmid prophages

The temperate Escherichia coli bacteriophage D6 can exist as a circular plasmid prophage, and we report here its 91,159bp complete genome sequence. It is a distant relative of the well-studied phage P1, but it is sufficiently different that it typifies a previously undescribed tailed phage type or cluster. Examination of the database of bacterial genome sequences revealed that phage P1 and D6 prophage plasmids are common in the Enterobacteriales, and in addition, previously described Salmonella phage SSU5 represents a different type of temperate tailed phage with a circular plasmid prophage that is also very common in this host order. This analysis also discovered additional divergent clusters of putative circular plasmid prophages within the two larger P1 and SSU5 groups (superclusters) that inhabit the Enterobacteriales as well as bacteria in several other orders in the Gamma-proteobacteria class. Very few of these sequences are annotated as putative prophages.

Towards the complete small RNome of Acinetobacter baumannii

In recent years, the Gram-negative bacterium Acinetobacter baumannii has garnered considerable attention for its unprecedented capacity to rapidly develop resistance to antibacterial therapeutics. This is coupled with the seemingly epidemic emergence of new hyper-virulent strains. Although strain-specific differences for A. baumannii isolates have been well described, these studies have primarily focused on proteinaceous factors. At present, only limited publications have investigated the presence and role of small regulatory RNA (sRNA) transcripts. Herein, we perform such an analysis, describing the RNA-seq-based identification of 78 A. baumannii sRNAs in the AB5075 background. Together with six previously identified elements, we include each of these in a new genome annotation file, which will serve as a tool to investigate regulatory events in this organism. Our work reveals that the sRNAs display high expression, accounting for >50 % of the 20 most strongly expressed genes. Through conservation analysis we identified six classes of similar sRNAs, with one found to be particularly abundant and homologous to regulatory, C4 antisense RNAs found in bacteriophages. These elements appear to be processed from larger transcripts in an analogous manner to the phage C4 molecule and are putatively controlled by two further sRNAs that are strongly antisense to them. Collectively, this study offers a detailed view of the sRNA content of A. baumannii, exposing sequence and structural conservation amongst these elements, and provides novel insight into the potential evolution, and role, of these understudied regulatory molecules.

P1 Ref Endonuclease: A Molecular Mechanism for Phage-Enhanced Antibiotic Lethality

Ref is an HNH superfamily endonuclease that only cleaves DNA to which RecA protein is bound. The enigmatic physiological function of this unusual enzyme is defined here. Lysogenization by bacteriophage P1 renders E. coli more sensitive to the DNA-damaging antibiotic ciprofloxacin, an example of a phenomenon termed phage-antibiotic synergy (PAS). The complementary effect of phage P1 is uniquely traced to the P1-encoded gene ref. Ref is a P1 function that amplifies the lytic cycle under conditions when the bacterial SOS response is induced due to DNA damage. The effect of Ref is multifaceted. DNA binding by Ref interferes with normal DNA metabolism, and the nuclease activity of Ref enhances genome degradation. Ref also inhibits cell division independently of the SOS response. Ref gene expression is toxic to E. coli in the absence of other P1 functions, both alone and in combination with antibiotics. The RecA proteins of human pathogens Neisseria gonorrhoeae and Staphylococcus aureus serve as cofactors for Ref-mediated DNA cleavage. Ref is especially toxic during the bacterial SOS response and the limited growth of stationary phase cultures, targeting aspects of bacterial physiology that are closely associated with the development of bacterial pathogen persistence.

Function of the N-terminal segment of the RecA-dependent nuclease Ref

The bacteriophage P1 Ref (recombination enhancement function) protein is a RecA-dependent, HNH endonuclease. It can be directed to create targeted double-strand breaks within a displacement loop formed by RecA. The 76 amino acid N-terminal region of Ref is positively charged (25/76 amino acid residues) and inherently unstructured in solution. Our investigation of N-terminal truncation variants shows this region is required for DNA binding, contains a Cys involved in incidental dimerization and is necessary for efficient Ref-mediated DNA cleavage. Specifically, Ref N-terminal truncation variants lacking between 21 and 47 amino acids are more effective RecA-mediated targeting nucleases. We propose a more refined set of options for the Ref-mediated cleavage mechanism, featuring the N-terminal region as an anchor for at least one of the DNA strand cleavage events.

The N-terminal domain of the repressor of Staphylococcus aureus phage Φ11 possesses an unusual dimerization ability and DNA binding affinity

Bacteriophage Φ11 uses Staphylococcus aureus as its host and, like lambdoid phages, harbors three homologous operators in between its two divergently oriented repressor genes. None of the repressors of Φ11, however, showed binding to all three operators, even at high concentrations. To understand why the DNA binding mechanism of Φ11 repressors does not match that of lambdoid phage repressors, we studied the N-terminal domain of the Φ11 lysogenic repressor, as it harbors a putative helix-turn-helix motif. Our data revealed that the secondary and tertiary structures of the N-terminal domain were different from those of the full-length repressor. Nonetheless, the N-terminal domain was able to dimerize and bind to the operators similar to the intact repressor. In addition, the operator base specificity, binding stoichiometry, and binding mechanism of this domain were nearly identical to those of the whole repressor. The binding affinities of the repressor and its N-terminal domain were reduced to a similar extent when the temperature was increased to 42°C. Both proteins also adequately dislodged a RNA polymerase from a Φ11 DNA fragment carrying two operators and a promoter. Unlike the intact repressor, the binding of the N-terminal domain to two adjacent operator sites was not cooperative in nature. Taken together, we suggest that the dimerization and DNA binding abilities of the N-terminal domain of the Φ11 repressor are distinct from those of the DNA binding domains of other phage repressors.

Toxin-antitoxin genes of the Gram-positive pathogen Streptococcus pneumoniae: so few and yet so many

Pneumococcal infections cause up to 2 million deaths annually and raise a large economic burden and thus constitute an important threat to mankind. Because of the increase in the antibiotic resistance of Streptococcus pneumoniae clinical isolates, there is an urgent need to find new antimicrobial approaches to triumph over pneumococcal infections. Toxin-antitoxin (TA) systems (TAS), which are present in most living bacteria but not in eukaryotes, have been proposed as an effective strategy to combat bacterial infections. Type II TAS comprise a stable toxin and a labile antitoxin that form an innocuous TA complex under normal conditions. Under stress conditions, TA synthesis will be triggered, resulting in the degradation of the labile antitoxin and the release of the toxin protein, which would poison the host cells. The three functional chromosomal TAS from S. pneumoniae that have been studied as well as their molecular characteristics are discussed in detail in this review. Furthermore, a meticulous bioinformatics search has been performed for 48 pneumococcal genomes that are found in public databases, and more putative TAS, homologous to well-characterized ones, have been revealed. Strikingly, several unusual putative TAS, in terms of components and genetic organizations previously not envisaged, have been discovered and are further discussed. Previously, we reported a novel finding in which a unique pneumococcal DNA signature, the BOX element, affected the regulation of the pneumococcal yefM-yoeB TAS. This BOX element has also been found in some of the other pneumococcal TAS. In this review, we also discuss possible relationships between some of the pneumococcal TAS with pathogenicity, competence, biofilm formation, persistence, and an interesting phenomenon called bistability.

Bacteriophage 434 Hex protein prevents recA-mediated repressor autocleavage

In a λ^imm434 lysogen, two proteins are expressed from the integrated prophage. Both are encoded by the same mRNA whose transcription initiates at the P_RM promoter. One protein is the 434 repressor, needed for the establishment and maintenance of lysogeny. The other is Hex which is translated from an open reading frame that apparently partially overlaps the 434 repressor coding region. In the wild type host, disruption of the gene encoding Hex destabilizes λ^imm434 lysogens. However, the hex mutation has no effect on lysogen stability in a recA⁻ host. These observations suggest that Hex functions by modulating the ability of RecA to stimulate 434 repressor autocleavage. We tested this hypothesis by identifying and purifying Hex to determine if this protein inhibited RecA‑stimulated autocleavage of 434 repressor in vitro. Our results show that in vitro a fragment of Hex prevents RecA-stimulated autocleavage of 434 repressor, as well as the repressors of the closely related phage P22. Surprisingly, Hex does not prevent RecA‑stimulated autocleavage of phage lambda repressor, nor the E. coli LexA repressor.

Inhibition of superinfection and the evolution of viral latency

Latent viruses generally defend their host cell against superinfection by nonlatent virulent mutants that could destroy the host cell. Superinfection inhibition thus seems to be a prerequisite for the maintenance of viral latency. Yet viral latency can break down when resistance to superinfection inhibition, known as ultravirulence, occurs. To understand the evolution of viral latency, we have developed a model that analyzes the epidemiology of latent infection in the face of ultravirulence. We show that latency can be maintained when superinfection inhibition and resistance against it coevolve in an arms race, which can result in large fluctuations in virulence. An example is the coevolution of the virulence and superinfection repressor protein of phage lambda (cI) and its binding target, the lambda oLoR operator. We show that this repressor/operator coevolution is the driving force for the evolution of superinfection immunity groups. Beyond latent phages, we predict analogous dynamics for any latent virus that uses a single repressor for the simultaneous control of virulence and superinfection.

Regions and residues of an asymmetric operator DNA interacting with the monomeric repressor of temperate mycobacteriophage L1

Previously, the repressor protein of mycobacteriophage L1 bound to two operator DNAs with dissimilar affinity. Surprisingly, the putative operator consensus sequence, 5'GGTGGa/cTGTCAAG, lacks the dyad symmetry reported for the repressor binding operators of lambda and related phages. To gain insight into the structure of the L1 repressor-asymmetric operator DNA complex, we have performed various in vitro experiments. A dimethyl sulfate protection assay revealed that five guanine bases, mostly distributed in the two adjacent major grooves of the 13 bp operator DNA helix, participate in repressor binding. Hydroxyl radical footprinting demonstrated that interaction between the repressor and operator DNA is asymmetric in nature and occurs primarily through one face of the DNA helix. Genetic studies not only confirmed the results of the dimethyl sulfate protection assay but also indicated that other bases in the 13 bp operator DNA are critical for repressor binding. Interestingly, repressor that weakly induced bending in the asymmetric operator DNA interacted with this operator as a monomer. The tertiary structure of the L1 repressor-operator DNA complex therefore appears to be distinct from those of the lambdoid phages even though the number of repressor molecules per operator site closely matched that of the lambda phage system.

Isolation and genomic characterization of the first phage infecting Iodobacteria: ϕPLPE, a myovirus having a novel set of features

The aquatic phage ϕPLPE infects a bacterium of the genus Iodobacter that are common inhabitants of rivers, streams and canals that produce violacein-like pigments. Our characterization of ϕPLPE reveals it to be a small, contractile-tailed phage whose 47.5 kb genome sequence is phylogenetically distant from all previously characterized phages. The genome has a generally modular organization (e.g. replication/recombination, structure/morphogenesis, lysis/lysogeny) and approximately half of its 84 open reading frames have no known homologues. It behaves as a virulent phage under the host growth conditions we have employed and, with the exception of an anti-repressor (ant) homologue, the genome lacks all the genes associated with a lysogenic lifestyle. Thus, either ϕPLPE was once a temperate phage that has lost most of its lysogeny cassette or it is a virulent phage that acquired an ant-like gene presumably for some function other than the control of lysogeny. The ϕPLPE genome has few bacterial gene homologues with the interesting exception of a putative acylhydrolase (acylase). This function has been implicated in bacterial quorum sensing since it degrades homoserine-lactone signalling molecules and can disrupt or modulate quorum signalling from either the emitter or its competitors. ϕPLPE may be an example of a phage co-opting components of the bacterial quorum-sensing apparatus to its own advantage.

Comprehensive comparative-genomic analysis of type 2 toxin-antitoxin systems and related mobile stress response systems in prokaryotes

Background

The prokaryotic toxin-antitoxin systems (TAS, also referred to as TA loci) are widespread, mobile two-gene modules that can be viewed as selfish genetic elements because they evolved mechanisms to become addictive for replicons and cells in which they reside, but also possess "normal" cellular functions in various forms of stress response and management of prokaryotic population. Several distinct TAS of type 1, where the toxin is a protein and the antitoxin is an antisense RNA, and numerous, unrelated TAS of type 2, in which both the toxin and the antitoxin are proteins, have been experimentally characterized, and it is suspected that many more remain to be identified.

Results

We report a comprehensive comparative-genomic analysis of Type 2 toxin-antitoxin systems in prokaryotes. Using sensitive methods for distant sequence similarity search, genome context analysis and a new approach for the identification of mobile two-component systems, we identified numerous, previously unnoticed protein families that are homologous to toxins and antitoxins of known type 2 TAS. In addition, we predict 12 new families of toxins and 13 families of antitoxins, and also, predict a TAS or TAS-like activity for several gene modules that were not previously suspected to function in that capacity. In particular, we present indications that the two-gene module that encodes a minimal nucleotidyl transferase and the accompanying HEPN protein, and is extremely abundant in many archaea and bacteria, especially, thermophiles might comprise a novel TAS. We present a survey of previously known and newly predicted TAS in 750 complete genomes of archaea and bacteria, quantitatively demonstrate the exceptional mobility of the TAS, and explore the network of toxin-antitoxin pairings that combines plasticity with selectivity.

Conclusion

The defining properties of the TAS, namely, the typically small size of the toxin and antitoxin genes, fast evolution, and extensive horizontal mobility, make the task of comprehensive identification of these systems particularly challenging. However, these same properties can be exploited to develop context-based computational approaches which, combined with exhaustive analysis of subtle sequence similarities were employed in this work to substantially expand the current collection of TAS by predicting both previously unnoticed, derived versions of known toxins and antitoxins, and putative novel TAS-like systems. In a broader context, the TAS belong to the resistome domain of the prokaryotic mobilome which includes partially selfish, addictive gene cassettes involved in various aspects of stress response and organized under the same general principles as the TAS. The "selfish altruism", or "responsible selfishness", of TAS-like systems appears to be a defining feature of the resistome and an important characteristic of the entire prokaryotic pan-genome given that in the prokaryotic world the mobilome and the "stable" chromosomes form a dynamic continuum.

Reviewers

This paper was reviewed by Kenn Gerdes (nominated by Arcady Mushegian), Daniel Haft, Arcady Mushegian, and Andrei Osterman. For full reviews, go to the Reviewers' Reports section.

Two different Panton-Valentine leukocidin phage lineages predominate in Japan

We determined the entire nucleotide sequence of phiSa2958-carrying Panton-Valentine leukocidin (PVL) gene, which was lysogenized in a sequence type 5 staphylococcal cassette chromosome mec (SCCmec) type II strain of methicillin-resistant Staphylococcus aureus (MRSA). Based on the nucleotide sequences of PVL phages, we developed PCRs to discriminate among five PVL phages, with a preliminary classification into two morphological groups (elongated-head type and icosahedral-head type) with four PCRs, including two PCRs for identifying the gene lineage between lukS-PV and the tail gene. The phages were then classified into five types by four PCRs identifying each phage-specific structure. With these PCRs, we examined the PVL phage types of 67 MRSA strains isolated in Japan from 1979 through 1985 and since 2000 and found that two morphologically distinct phages were predominant in Japan. The icosahedral-head-type phage, represented by the phi108PVL type, was identified for 39 of 53 strains isolated from 1979 through 1985. Of 26 other Japanese isolates, 25 belonged either definitively or presumably to elongated-head types as follows: 3 belonged to the phiSa2958 type; 8 were determined to belong to an elongated-head type, but a determination of greater specificity was not made; and 14 belonged to a phiSa2958-like phage of unknown type. We induced prophages by treatment with mitomycin C from six strains of the phiSa2958 type or of phiSa2958-like unknown-type phages; five of six strains carried intact PVL-carrying phages, which can infect other S. aureus strains and might generate novel PVL-positive strains of S. aureus. That various SCCmec elements were carried by different strains of the same phage type suggests that S. aureus strains might independently acquire PVL phages before they acquire various SCCmec elements.

Transcription-termination-mediated immunity and its prevention in bacteriophage SfV of Shigella flexneri

The temperate phage SfV encodes the genes responsible for the serotype conversion of Shigella flexneri strains from serotype Y to 5a. Bacteriophages often encode proteins that prevent subsequent infection by homologous phages; the mechanism by which this is accomplished is referred to as superinfection immunity. The serotype conversion mediated following lysogenization of SfV is one such mechanism. Another mechanism is the putative lambda-like CI protein within SfV. This study reports the characterization of a third superinfection mechanism, transcription termination, in SfV. The presence of a small immunity-mediating RNA molecule, called CI RNA, and its essential role in the establishment of immunity, is shown. The novel role of the gene orf77, located immediately downstream from the transcription termination region, in inhibiting the establishment of CI RNA-mediated immunity is also presented.

Repressor of temperate mycobacteriophage L1 harbors a stable C-terminal domain and binds to different asymmetric operator DNAs with variable affinity

BACKGROUND

Lysogenic mode of life cycle of a temperate bacteriophage is generally maintained by a protein called 'repressor'. Repressor proteins of temperate lambdoid phages bind to a few symmetric operator DNAs in order to regulate their gene expression. In contrast, repressor molecules of temperate mycobacteriophages and some other phages bind to multiple asymmetric operator DNAs. Very little is known at present about the structure-function relationship of any mycobacteriophage repressor.

RESULTS

Using highly purified repressor (CI) of temperate mycobacteriophage L1, we have demonstrated here that L1 CI harbors an N-terminal domain (NTD) and a C-terminal domain (CTD) which are separated by a small hinge region. Interestingly, CTD is more compact than NTD at 25 degrees C. Both CTD and CI contain significant amount of alpha-helix at 30 degrees C but unfold partly at 42 degrees C. At nearly 200 nM concentration, both proteins form appreciable amount of dimers in solution. Additional studies reveal that CI binds to O64 and OL types of asymmetric operators of L1 with variable affinity at 25 degrees C. Interestingly, repressor-operator interaction is affected drastically at 42 degrees C. The conformational change of CI is most possibly responsible for its reduced operator binding affinity at 42 degrees C.

CONCLUSION

Repressors encoded by mycobacteriophages differ significantly from the repressor proteins of lambda and related phages at functional level but at structural level they are nearly similar.

DNA recombination with a heterospecific Cre homolog identified from comparison of the pac-c1 regions of P1-related phages

Sequencing of the 7 kb immC region from four P1-related phages identified a novel DNA recombinase that exhibits many Cre-like characteristics, including recombination in mammalian cells, but which has a distinctly different DNA specificity. DNA sequence comparison to the P1 immC region showed that all phages had related DNA terminase, C1 repressor and DNA recombinase genes. Although these genes from phages P7, phi(w39) and p15B were highly similar to those from P1, those of phage D6 showed significant divergence. Moreover, the D6 sequence showed evidence of DNA deletion and substitution in this region relative to the other phages. Characterization of the D6 site-specific DNA recombinase (Dre) showed that it was a tyrosine recombinase closely related to the P1 Cre recombinase, but that it had a distinct DNA specificity for a 32 bp DNA site (rox). Cre and Dre are heterospecific: Cre did not catalyze recombination at rox sites and Dre did not catalyze recombination at lox sites. Like Cre, Dre catalyzed both integrative and excisive recombination and required no other phage-encoded proteins for recombination. Dre-mediated recombination in mammalian cells showed that, like Cre, no host bacterial proteins are required for efficient Dre-mediated site-specific DNA recombination.

Genome of bacteriophage P1

P1 is a bacteriophage of Escherichia coli and other enteric bacteria. It lysogenizes its hosts as a circular, low-copy-number plasmid. We have determined the complete nucleotide sequences of two strains of a P1 thermoinducible mutant, P1 c1-100. The P1 genome (93,601 bp) contains at least 117 genes, of which almost two-thirds had not been sequenced previously and 49 have no homologs in other organisms. Protein-coding genes occupy 92% of the genome and are organized in 45 operons, of which four are decisive for the choice between lysis and lysogeny. Four others ensure plasmid maintenance. The majority of the remaining 37 operons are involved in lytic development. Seventeen operons are transcribed from sigma(70) promoters directly controlled by the master phage repressor C1. Late operons are transcribed from promoters recognized by the E. coli RNA polymerase holoenzyme in the presence of the Lpa protein, the product of a C1-controlled P1 gene. Three species of P1-encoded tRNAs provide differential controls of translation, and a P1-encoded DNA methyltransferase with putative bifunctionality influences transcription, replication, and DNA packaging. The genome is particularly rich in Chi recombinogenic sites. The base content and distribution in P1 DNA indicate that replication of P1 from its plasmid origin had more impact on the base compositional asymmetries of the P1 genome than replication from the lytic origin of replication.

immX immunity region of rhizobium phage 16-3: two overlapping cistrons of repressor function

16-3 is a temperate phage of the symbiotic nitrogen-fixing bacterium Rhizobium meliloti 41. Its prophage state and immunity against superinfection by homoimmune phages are governed by a complex set of controls: the immC and immX repressor systems and the avirT element are all located in well-separated, distinct regions which span 25 kb on the bacteriophage chromosome. The anatomy and function of the immC region are well documented; however, fewer analyses have addressed the immX and avirT regions. We focused in this paper on the immX region and dissected it into two major parts: X(U/L) and X(V). The X(U/L) part (0.6 kb) contained two overlapping cistrons, X(U) and X(L), coding for proteins pXU and pXL, respectively. Inactivation of either gene inactivated the repressor function of the immX region. Loss-of-function mutants of X(U) and X(L) complemented each other in trans in double lysogens. The X(V) part (1 kb) contained a target for X(U/L) repressor action. Mutations at three sites in X(V) led to various degree of ImmX insensitivity in a hierarchic manner. Two sites (X(V1) and X(V3)) exhibited the inverted-repeat structures characteristic of many repressor binding sites. However, X(V1) could also be folded into a transcription terminator. Of the two immunity regions of 16-3, immX seems to be unique both in its complex genetic anatomy and in its sequence. To date, no DNA or peptide sequence homologous to that of ImmX has been found in the data banks. In contrast, immC shares properties of a number of immunity systems commonly found in temperate phages.

Prophage genomics

The majority of the bacterial genome sequences deposited in the National Center for Biotechnology Information database contain prophage sequences. Analysis of the prophages suggested that after being integrated into bacterial genomes, they undergo a complex decay process consisting of inactivating point mutations, genome rearrangements, modular exchanges, invasion by further mobile DNA elements, and massive DNA deletion. We review the technical difficulties in defining such altered prophage sequences in bacterial genomes and discuss theoretical frameworks for the phage-bacterium interaction at the genomic level. The published genome sequences from three groups of eubacteria (low- and high-G+C gram-positive bacteria and gamma-proteobacteria) were screened for prophage sequences. The prophages from Streptococcus pyogenes served as test case for theoretical predictions of the role of prophages in the evolution of pathogenic bacteria. The genomes from further human, animal, and plant pathogens, as well as commensal and free-living bacteria, were included in the analysis to see whether the same principles of prophage genomics apply for bacteria living in different ecological niches and coming from distinct phylogenetical affinities. The effect of selection pressure on the host bacterium is apparently an important force shaping the prophage genomes in low-G+C gram-positive bacteria and gamma-proteobacteria.

Escherichia coli SspA is a transcription activator for bacteriophage P1 late genes

The stringent starvation protein A (SspA), an Escherichia coli RNA polymerase (RNAP)-associated protein, has been reported to be essential for lytic growth of bacteriophage P1. Unlike P1 early promoters, P1 late promoters are not recognized by RNAP alone. A phage-encoded early protein, Lpa (late promoter activator protein, formerly called gp10), has been shown to be required for P1 late transcription in vivo. Here, we demonstrate that SspA is a transcription activator for P1 late genes. Our results indicated that Lpa is not limiting in an sspA mutant. However, the transcription of P1 late genes was deficient in an sspA mutant in vivo. We demonstrated that SspA/Lpa are required for transcription activation of the P1 late promoter Ps in vitro. In addition, SspA and Lpa were shown to facilitate the binding of RNAP to Ps late promoter DNA. Activation of late transcription by SspA/Lpa was dependent on holoenzyme containing sigma70 but not sigmaS, indicating that the two activators discriminate between the two forms of the holoenzyme. Furthermore, P1 early gene expression was downregulated in the wild-type background, whereas it persisted in the sspA mutant background, indicating that SspA/Lpa mediate the transcriptional switch from the early to the late genes during P1 lytic growth. Thus, this work provides the first evidence for a function of the E. coli RNAP-associated protein SspA.

Development of a thermally regulated broad-spectrum promoter system for use in pathogenic gram-positive species

Selectively regulating gene expression is an essential molecular tool that is lacking for many pathogenic gram-positive bacteria. In this report, we describe the evaluation of a series of promoters regulated by the bacteriophage P1 temperature-sensitive C1 repressor in Enterococcus faecium, Enterococcus faecalis, and Staphylococcus aureus. Using the lacZ gene to monitor gene expression, we examined the strength, basal expression, and induced expression of synthetic promoters carrying C1 operator sites. The promoters exhibited extremely low basal expression and, under inducing conditions, gave high levels of expression (100- to 1,000-fold induction). We demonstrate that the promoter system could be modulated by temperature and showed rapid induction and that the mechanism of regulation occurred at the level of transcription. Controlled expression with the same constructs was also demonstrated in the gram-negative bacterium Escherichia coli. However, low basal expression and the ability to achieve derepression were dependent on both the number of mismatches in the C1 operator sites and the promoter driving c1 expression. Since the promoters were designed to contain conserved promoter elements from gram-positive species and were constructed in a broad-host-range plasmid, this system will provide a new opportunity for controlled gene expression in a variety of gram-positive bacteria.

Complete genomic sequence of SfV, a serotype-converting temperate bacteriophage of Shigella flexneri

Bacteriophage SfV is a temperate serotype-converting phage of Shigella flexneri. SfV encodes the factors involved in type V O-antigen modification, and the serotype conversion and integration-excision modules of the phage have been isolated and characterized. We now report on the complete sequence of the SfV genome (37,074 bp). A total of 53 open reading frames were predicted from the nucleotide sequence, and analysis of the corresponding proteins was used to construct a functional map. The general organization of the genes in the SfV genome is similar to that of bacteriophage lambda, and numerous features of the sequence are described. The superinfection immunity system of SfV includes a lambda-like repression system and a P4-like transcription termination mechanism. Sequence analysis also suggests that SfV encodes multiple DNA methylases, and experiments confirmed that orf-41 encodes a Dam methylase. Studies conducted to determine if the phage-encoded methylase confers host DNA methylation showed that the two S. flexneri strains analyzed encode their own Dam methylase. Restriction mapping and sequence analysis revealed that the phage genome has cos sites at the termini. The tail assembly and structural genes of SfV show homology to those of phage Mu and Mu-like prophages in the genome of Escherichia coli O157:H7 and Haemophilus influenzae. Significant homology (30% of the genome in total) between sections of the early, regulatory, and structural regions of the SfV genome and the e14 and KpLE1 prophages in the E. coli K-12 genome were noted, suggesting that these three phages have common evolutionary origins.

Controlled expression in Klebsiella pneumoniae and Shigella flexneri using a bacteriophage P1-derived C1-regulated promoter system

The utility of promoters regulated by the bacteriophage P1 temperature-sensitive C1 repressor was examined in Shigella flexneri and Klebsiella pneumoniae. Promoters carrying C1 operator sites driving LacZ expression had induction/repression ratios of up to 240-fold in S. flexneri and up to 50-fold in K. pneumoniae. The promoters exhibited remarkably low basal expression, demonstrated modulation by temperature, and showed rapid induction. This system will provide a new opportunity for controlled gene expression in enteric gram-negative bacteria.

Dual regulatory control of a particle maturation function of bacteriophage P1

A unique arrangement of promoter elements was found upstream of the bacteriophage P1 particle maturation gene (mat). A P1-specific late-promoter sequence with conserved elements located at positions -22 and -10 was expected from the function of the gene in phage morphogenesis. In addition to a late-promoter sequence, a -35 element and an operator sequence for the major repressor protein, C1, were found. The -35 and -10 elements constituted an active Escherichia coli sigma(70) consensus promoter, which was converted into a P1-regulated early promoter by the superimposition of a C1 operator. This combination of early- and late-promoter elements regulates and fine-tunes the expression of the particle maturation gene. During lysogenic growth the gene is turned off by P1 immunity functions. Upon induction of lytic growth, the expression of mat starts simultaneously with the expression of other C1-regulated P1 early functions. However, while most of the latter functions are downregulated during late stages of lytic growth the expression of mat continues throughout the entire lytic growth cycle of bacteriophage P1. Thus, the maturation function has a head start on the structural components of the phage particle.

Characterization of the major control region of Vibrio cholerae bacteriophage K139: immunity, exclusion, and integration

The temperate bacteriophage K139 is highly associated with pathogenic O1 Vibrio cholerae strains. The nucleotide sequence of the major control region of K139 was determined. The sequences of four (cox, cII, cI, and int) of the six deduced open reading frames and their gene order indicated that K139 is related to the P2 bacteriophage family. Two genes of the lysogenic transcript from the mapped promoter PL encode homologs to the proteins CI and Int, with deduced functions in prophage formation and maintenance. Between the cI and int genes, two additional genes were identified: orf2, which has no significant similarity to any other gene, and the formerly characterized gene glo. Further analysis revealed that Orf2 is involved in preventing superinfection. In a previous report, we described that mutations in glo cause an attenuation effect in the cholera mouse model (J. Reidl and J. J. Mekalanos, Mol. Microbiol. 18:685-701, 1995). In this report, we present strong evidence that Glo participates in phage exclusion. Glo was characterized to encode a 13.6-kDa periplasmic protein which inhibits phage infection at an early step, hence preventing reinfection of vibriophage K139 into K139 lysogenic cells. Immediately downstream of gene int, the attP site was identified. Upon analysis of the corresponding attB site within the V. cholerae chromosome, it became evident that phage K139 is integrated between the flagellin genes flaA and flaC of O1 El Tor and O139 V. cholerae lysogenic strains.

The primary immunity determinant in modulating the lysogenic immunity of the filamentous bacteriophage cf

Bacteriophage cf is the first single-stranded DNA phage that has been shown to set up a stable lysogenic state with its genome integrated into the host chromosome. From the isolation and characterization of a virulent mutant, cf-tv2, we report the first investigation into the mechanisms of the immunity established by the filamentous bacteriophage. The mutation in cf-tv2 enables the phage to produce plaques on lawns of cf lysogenic cells. The mutation was defined as a 49-nucleotide deletion located in a 0.59 kb NcoI/KpnI fragment of cf replicative form DNA. Two messages, cM1 and cM2, transcribed from the immunity region of wild-type cf but in opposite directions, were detected. In cf-tv2, the 49-nucleotide deletion abolishes cM2 transcription. The primer extension assay suggests a possible RNA-RNA interaction directed by base-pairing of the cM1 and cM2 RNAs. A frameshift mutation of the open reading frame ORF 165, encoded by cM2, resulted in a 10(5) plating efficiency on the cf lysogen. These observations suggest that both RNA-RNA interaction and repressor protein inhibition are involved in the mechanism of cf immunity. A model is proposed for the regulation of cf immunity.

Transcriptional silencing by the mycobacteriophage L5 repressor

The success of a temperate bacteriophage is dependent upon its ability to completely shut down expression of its lytic genes during lysogenic growth. Mycobacteriophage L5 accomplishes this by an atypical phage repressor, gp71, which binds to multiple asymmetric DNA sites. L5 gp71 regulates transcription initiation at an early lytic promoter, Pleft, but also affects downstream gene expression at 'stoperator' sites in the phage genome. The L5 genome is replete with stoperator sites located within short intergenic spaces in both the early and late lytic operons and oriented specifically with respect to transcription. Binding of gp71 to these sites results in a strong orientation-dependent polar effect on downstream gene expression and global silencing of prophage gene expression.

Multiple regulatory mechanisms controlling phage-plasmid P4 propagation

Bacteriophage P4 autonomous replication may result in the lytic cycle or in plasmid maintenance, depending, respectively, on the presence or absence of the helper phage P2 genome in the Escherichia coli host cell. Alternatively, P4 may lysogenize the bacterial host and be maintained in an immune-integrated condition. A key step in the choice between the lytic/plasmid vs. the lysogenic condition is the regulation of P4 alpha operon. This operon may be transcribed from two promoters, PLE and PLL, and encodes both immunity (promoter proximal) and replication (promoter distal) functions. PLE is a constitutive promoter and transcription of the downstream replication genes is regulated by transcription termination. The trans-acting immunity factor that controls premature transcription termination is a short RNA encoded in the PLE proximal part of the operon. Expression of the replication functions in the lytic/plasmid condition is achieved by activation of the PLL promoter. Transcription from PLL is insensitive to the termination mechanism that acts on transcription starting from PLE.PLL is also negatively regulated by P4 orf88, the first gene downstream of PLL. An additional control on P4 DNA replication is exerted by the P4 cnr gene product.

Multiple initiation mechanisms adapt phage T4 DNA replication to physiological changes during T4's development

We summarize the evidence for multiple pathways to initiate phage T4 DNA replication. In any infecting chromosome, leading DNA strands can be primed from pre-replicative transcripts, independent of primase activity, at one of several origins. Within each origin region, there are multiple RNA-DNA transition sites. However, the priming potential at each single site is very low. Our results suggest that origin transcripts can become primers for leading strand DNA synthesis without being processed, but that a promoter-proximal segment of each origin transcript plays an important structural role, as a proposed wedge, in the transition from RNA to DNA synthesis. Two recombination-dependent pathways render subsequent phage T4 DNA replication independent of transcription. The first of these requires proteins that are synthesized during the pre-replicative phase of infection. It is active as soon as the first growing points, initiated at origins, have reached a chromosomal end. The other one requires at least one late protein: endonuclease VII, a resolvase that cuts recombinational junctions. The latter pathway can bypass primase deficiencies by allowing retrograde DNA synthesis without Okazaki pieces. We discuss the integration of these multiple and redundant pathways into the developmental program of T4. Competition between these initiation mechanisms and with other DNA transactions allows for integration of replication controls with transcription, recombination and packaging of the DNA.