A widespread family of polymorphic toxins encoded by temperate phages
© Jamet et al. 2017
Received: 29 May 2017
Accepted: 7 August 2017
Published: 29 August 2017
Polymorphic toxins (PTs) are multi-domain bacterial exotoxins belonging to distinct families that share common features in terms of domain organization. PTs are found in all major bacterial clades, including many toxic effectors of type V and type VI secretion systems. PTs modulate the dynamics of microbial communities by killing or inhibiting the growth of bacterial competitors lacking protective immunity proteins.
In this work, we identified a novel widespread family of PTs, named MuF toxins, which were exclusively encoded within temperate phages and their prophages. By analyzing the predicted proteomes of 1845 bacteriophages and 2464 bacterial genomes, we found that MuF-containing proteins were frequently part of the DNA packaging module of tailed phages. Interestingly, MuF toxins were abundant in the human gut microbiome.
Our results uncovered the presence of the MuF toxin family in the temperate phages of Firmicutes. The MuF toxin family is likely to play an important role in the ecology of the human microbiota where pathogens and commensal species belonging to the Firmicutes are abundant. We propose that MuF toxins could be delivered by phages into host bacteria and either influence the lysogeny decision or serve as bacterial weapons by inhibiting the growth of competing bacteria.
KeywordsBacterial genomics Phages Bacterial toxins
Polymorphic toxins (PTs) are multi-domain proteins involved in competition between bacteria and in pathogenesis . Most lineages of bacteria encode at least one PT system in their genome . PTs encompass colicins, soluble pyocins, toxic effectors of type V secretion systems (T5SSs), some toxic effectors of type VI secretion systems (T6SSs), and MafB toxins [1, 3–5]. A PT family is defined by an N-terminal region harboring one or more conserved domains. For instance, a domain of unknown function named DUF1020 (PF06255) is found at the N-termini of all MafB toxins , whereas phage late control gene D protein (Phage_GPD; PF05954) and phage base plate assembly protein (Phage_base_V; PF04717) domains are found at the N-termini of VgrG toxins . In contrast, the C-terminus of PTs harbors a set of diverse C-terminal toxin domains, which can have homologs in other distinct PT families . This shared pool of toxin domains between families is a hallmark of PT systems, with more than 150 distinct toxin domains identified so far . Most PTs have RNase, DNase, peptidase, or other protein-modifying activities . PTs can be secreted by various secretion systems. CdiA and BcpA toxins are secreted by CdiB and BcpB, respectively, and belong to the two-partner secretion protein family (i.e., T5SS). CdiA and BcpA toxins require a direct contact for inhibition of bacterial competitors, hence the name contact-dependent inhibition (CDI) system [7, 8]. Certain PT families, including VgrG, Hcp, and PAAR proteins, have structural roles in the T6SS machinery in addition to their toxic activity . In these families, there are two main domain architectures: canonical proteins without a C-terminal extension and PTs bearing a C-terminal extension with toxic activity [6, 10–12]. Besides, it has been demonstrated that a C-terminal extension of a VgrG protein of an enteroaggregative Escherichia coli strain did not harbor a toxic domain but mediated binding and transport of a toxic effector .
When a bacterium produces an antibacterial toxin, it needs to protect itself from autointoxication and prevent self-inhibition. In most cases a small open reading frame (ORF) encoding a specific immunity protein is located immediately downstream of the toxin gene . PTs modulate the dynamics of microbial communities by killing or inhibiting the growth of competitors lacking the cognate immunity protein. For instance, the predominance of E. coli strain EC93 in the intestine of some commercial rats has been linked to CdiA toxin production . Indeed, CdiA of EC93 was shown to inhibit the growth of E. coli K12 strains . In Burkholderia thailandensis, BcpA is also involved in inhibiting growth of neighboring bacteria through contact [8, 16]. In addition, the Bcp system allows kin discrimination through the immunity proteins produced by bacteria . Furthermore, delivery of BcpA to immune bacteria mediates a contact-dependent signaling that promotes cooperative behavior such as biofilm formation . CdiA and BcpA toxins have only been shown to modulate competitive or cooperative relations between bacteria of the same species. In contrast, effectors of T6SS are involved in both intra- and inter-species competition and can also be delivered into eukaryotic cells [12, 18]. Indeed, a VgrG protein of Vibrio cholerae is responsible for remodeling of the actin cytoskeleton when injected into eukaryotic host cells , whereas a VgrG toxin of another Vibrio strain hydrolyzes the cell wall of Gram-negative competitors .
In comparison with the wealth of studies aiming to characterize PT systems in various genera of Proteobacteria (e.g ., Escherichia, Pseudomonas, Burkholderia, Vibrio, and Neisseria), only the WXG/LXG and Rhs PT families have been studied in monoderms [19, 20]. Monoderms encompass Firmicutes and Actinobacteria and are pivotal in human health. Indeed, they are major constituents of human microbiota accounting for > 50% of the species recovered from skin, nose, stomach, and vagina , and they also include major human pathogens . Several WXG/LXG toxins have been described in Bacillus subtilis [19, 23], in Staphylococcus aureus , and in Streptococcus intermedius . The role of these toxins in inter-bacterial competition has been demonstrated in the last two species. In addition, an Rhs toxin named WapA confers a competitive advantage in competition assays  and is involved in kin discrimination in B. subtilis .
In this study, we provide an in-depth description of a new family of PTs harboring a domain of the MuF superfamily in their N-terminal region; this family has the unusual feature of being associated with temperate phages. Viruses that infect bacteria (hereafter designated “phages”) may be “virulent,” and thus restricted to act through the lytic cycle, or “temperate.” The latter may either behave like virulent phages or integrate the bacterial chromosomes as prophages. Bacteria harboring prophages are called lysogens and account for nearly half of the sequenced bacteria . Most of the known phages are tailed phages belonging to the Caudovirales order . Overall, we found that 35% of the 1753 sequenced tailed phages and 30% of the 2622 prophages harbored a muf gene. Among 1515 muf genes, 13% encode toxin domains. The presence of a PT system in phages could have important implications for phage biology and microbial population dynamics.
The MuF domain-containing proteins constitute a novel polymorphic toxin family
Definition of the MuF superfamily of proteins
To detail preliminary indications of a strong association of muf genes with phages, we searched for them in bacteriophage and bacterial genomes. We used protein profiles of the four MuF domain families (Additional file 2: Table S3) to retrieve all proteins containing a MuF domain from the predicted proteomes of 1845 bacteriophages and 2464 bacterial genomes (Additional file 2: Tables S1, S2). Altogether, we identified 614 and 901 MuF proteins in bacteriophage genomes and bacterial chromosomes, respectively (Additional file 2: Tables S1, S2, S4, S5). With rare exceptions, phages had only one muf gene (614 muf genes were retrieved from 611 distinct phage genomes).
Domain architecture of MuF domain-containing proteins
As expected for PTs, most toxin domains found among MuF proteins had homologs in other PT families (75%, e.g., cd13442, Ntox50, and EndoU_bacteria domains) (Fig. 2d, Additional file 2: Table S7). More than half of the 191 putative MuF toxins had a nuclease domain, and 20% presented a metallopeptidase domain. ADP-ribosyl transferase and RelA-like domains were also frequently identified (Additional file 2: Tables S4, S5, S7).
Intriguingly, phages of Proteobacteria had many MuF proteins harboring a C-terminal extension without known domains (Fig. 2c and Additional file 3: Figure S2). We defined a new domain (termed Ct_MAD for C-terminal MuF Associated Domain) present in one third of these C-terminal extensions (Additional file 2: Tables S3–S5). However, it remains to be determined if the Ct_MAD domain and the remaining C-terminal extensions without known domains correspond to novel toxin domains.
Arguments in favor of MuF toxicity
If a bacterium produces a MuF protein with a toxic activity targeting bacterial cytosolic compounds (i.e., a nuclease), it should also produce a protective immunity protein to prevent self-intoxication. ORFs encoding immunity proteins are difficult to identify because they are often very short (less than 150 aa) and do not contain known domains [2, 14]. Most (89.4% of the 191 muf toxin genes) of the genes encoding MuF toxins identified in our study are followed by a small ORF potentially encoding a polypeptide of less than 150 aa (Additional file 4: Figure S3). In contrast, only 25.5% of the 952 genes encoding short MuF proteins are followed by a similarly small ORF (Additional file 4: Figure S3). MuF toxin genes were more likely than short MuF genes to be associated with small ORFs (p < 0.0001, two-tailed Student’s t test for the comparison of the lengths of ORFs downstream of muf genes) and supports the hypothesis that these small ORFs would encode immunity proteins, necessary only if a toxin domain is present in the MuF protein.
Furthermore, we identified one instance where a MuF toxin C-terminal region together with the small ORF downstream are homologous to a toxin-immunity module of a MafB toxin . As expected, since MuF and Maf are distinct PT systems, the N-terminal regions of the MuF and MafB toxins are unrelated (Additional file 5: Figure S4). Interestingly, this MuF toxin is encoded by a Streptococcus mitis phage , and the MafB toxin is encoded on a genomic island of a Neisseria meningitidis strain. Since both species share the same niche in the human nasopharynx, the possibility of DNA exchange between them may be advocated.
Distribution of MuF proteins and toxins
MuF proteins and toxins are encoded by bacteriophages and bacterial genomes
The presence of muf genes in our bacteriophage dataset led us to investigate in a systematic manner the association of muf genes harbored by bacterial genomes with chromosomally integrated phages (i.e., prophages). First, we identified prophages in all bacterial chromosomes (see Methods) and found that 98% of genomes encoding MuF proteins were lysogens. Among the 901 muf genes identified in these bacterial chromosomes, most (90%) were located within prophages and 3% within putative prophage remnants (elements smaller than 18 kb; see Methods). Among the remaining, nearly half were found to be located close to genes encoding proteins with similarity to portal or terminase proteins, and could also correspond to prophage remnants (see Methods). Altogether, almost all (96%) of the muf genes found in bacterial genomes were associated with either complete prophages or prophage remnants (Additional file 2: Table S6). Reciprocally, around 30% of the 2622 identified prophages (>18 kb) encoded MuF proteins.
MuF toxins are associated with temperate tailed phages
We predicted the lifestyle (virulent versus temperate) of 979 Caudovirales phages (see Methods, Additional file 2: Table S1), which allowed us to determine that virulent phages encoded very few MuF proteins and contained no single MuF protein with a known toxin domain (Fig. 5). In contrast, many temperate phages infecting Firmicutes encoded MuF proteins with toxin domains (Additional file 3: Figure S2C). In Proteobacteria, MuF proteins with a Ct_MAD domain were also restricted to temperate phages (Additional file 2: Table S4). Since MuF toxins are restricted to temperate phages, and the genes found in bacterial genomes are in prophages, we hypothesized that there is most likely a link between lysogeny and MuF proteins carrying toxin domains.
Among the prophages harboring MuF toxin genes, several have been shown to be inducible, including SF370.1 of Streptococcus pyogenes , SM1 of Streptococcus mitis , StB12 of Staphylococcus hominis , and several prophages of Enterococcus faecalis . Hence, MuF toxins are present in fully functional prophages. Furthermore, in a proteomic analysis of the StB12 phage, the MuF toxin was detected .
Genes encoding MuF toxins are present in the human gut microbiota
Our results on bacteriophages and bacterial genomes demonstrated a strong association of MuF toxins with both temperate phages and prophages of Firmicutes. Hence, these toxins could play an important role in gut microbiota, especially since more than half of the bacteria found in the gut belong to the Firmicutes phylum . We thus searched for MuF proteins in human gut microbiomes, using a non-redundant catalog of 10 million proteins . We found that 51% of the 6406 MuF proteins present in the human gut catalog have a C-terminal extension and 67% of these have known toxin domains (Additional file 2: Table S8). These results strongly favor a role for MuF toxins in the bacterial population ecology of the human gut microbiome by influencing bacteria-phage interactions.
In this work, we discovered many PTs associated with MuF domains on temperate phages. The most studied MuF protein is encoded by Bacillus subtilis phage SPP1 (Gp7) [31–33], where it is present in one to two copies per capsid in a complex with portal proteins . Gp7 harbors no toxin domain. To our knowledge, the only MuF protein with a toxin domain studied so far is the EFV toxin (Q838U8) encoded by a lysogenic phage of the Enterococcus faecalis strain V583 . EFV toxin is a MuF1 protein with an ADP-ribosyl transferase activity, which is toxic when expressed in yeast .
We hypothesize that the toxin domain of MuF proteins could be delivered into host bacteria during phage DNA injection. In contrast to the widespread distribution of muf genes, those encoding a toxin domain were present only in temperate phages and were vastly over-represented in prophages of Firmicutes. These findings suggest that MuF toxins could influence the lysogenic decision. Alternatively, MuF toxins could act as molecular weapons in inter-bacterial competition.
There are many examples of prophages encoding toxins with anti-eukaryotic activities involved in bacterial virulence . These include E. coli prophages encoding Shiga toxin, filamentous phage CTXφ encoding the cholera toxin, or S. aureus prophages encoding the Panton-Valentine leukocidin . Besides, several prophage toxin-antitoxin (TA) genes were reported [44–47], e.g., in extra-chromosomal prophages P1 and N15, where they may stabilize the presence of prophages though post-segregational killing, as is often the case in plasmids . In addition, multi-protein structures termed “tailocins” have been described in the genomes of Pseudomonas. Tailocins are morphologically similar to phage tails and exhibit a bacteriotoxic activity via direct perforation of the cell envelope. Bacteriocins are often encoded close to the tail cassette .
The predicted targets of MuF toxins encompass RNA molecules, for putative ribonucleases, and (p)ppGpp metabolism, for putative RelA-like MuF toxins (Additional file 2: Table S7). Both substrates are likely to be primarily encountered in the host cytoplasm. We hypothesize that MuF toxins are loaded into the viral head to be delivered to bacteria concomitantly with the injection of phage DNA into the host cytoplasm. Infecting phages are known to inject proteins along with their DNA in the bacterial cytoplasm, as demonstrated for several ADP-ribosyl transferase proteins of virulent T4 phage . In support of this hypothesis, 14% of the 191 MuF toxins identified in our study are also putative ADP-ribosyl transferases (Additional file 2: Tables S4, S5).
Some PTs with nuclease activity are known to be highly toxic and serve as antibacterial weapons. For instance, it has been demonstrated that the cd13442 toxin domain of the CdiA-2 toxin (I1WVY3) of Burkholderia pseudomallei 1026b has tRNA nuclease activity and is required for contact-dependent growth inhibition of neighboring competitors [51, 52]. Such highly toxic MuFs could enhance the role of temperate phages in bacterial competition (Fig. 6b). It has been proposed that bacteria use their prophages to remove competitors from a niche , and that this could impact bacterial pathogenesis . This mechanism would work as follows: a lysogenic population would produce phages by lysis of a subpopulation of cells and thus kill sensitive competitors (siblings are protected by their own prophages encoding an immunity protein). However, the effect of this mechanism is short-lived, since the phage will eventually lysogenize some competitors, rendering them immune to superinfection . This is where the MuF toxin may play a role. If infection by the phage results in lysogeny, the incoming toxin could, at least temporarily, inhibit bacterial growth and provide a competitive advantage to the population carrying the prophage (Fig. 6b). Hence, phages would either kill or inhibit the growth of the competing population.
The abundance of MuF PTs suggests the existence of a toxin phage-delivery system and raises intriguing possibilities for their function in the context of phage-bacteria and inter-bacterial interactions. The latter could be important in environments such as human-associated ecosystems, where multiple strains of the same species or closely related species compete for the same niche . Strikingly, MuF toxins are particularly abundant in phages and prophages of Firmicutes that are significant members of the human microbiota in various niches  and include major human pathogens. Indeed, we found MuF toxins in important human pathogens including Streptococcus pyogenes and Enterococcus faecalis.
The sequences and corresponding annotations of 1845 complete bacteriophage genomes and 2462 complete bacterial genomes were retrieved from GenBank Refseq (last accessed September 2016) .
The phage dataset contained 1753 Caudovirales (95%), of which 53%, 27%, and 19% were Siphoviridae, Myoviridae, and Podoviridae, respectively. The order, the family, the type of nucleic acid, and the bacterial host of bacteriophages were extracted from the GenBank files. The lifestyle of Caudovirales phages was predicted using Phage Classification Tool Set (PHACTS) v0.3  (as in ). PHACTS predicts the lifestyle of a phage (i.e., virulent or temperate) from genomic data using both a similarity and a supervised random forest algorithm. The classification is based on the similarity to phages with known lifestyles in a manually curated database. Predictions were considered as confident only if the averaged probability score of the predicted lifestyle is two standard deviations (SDs) away from the averaged probability score of the other lifestyle, as recommended by the authors (who claim a precision rate of 99% with this parameter). We classified with such high confidence 55% of the dataset (i.e., 410 virulent and 600 temperate phages). The general characteristics and the result of lifestyle prediction for each phage are listed in Additional file 2: Table S1.
Prophages were detected in bacterial chromosomes as in  using Phage Finder v4.6 . The elements larger than 18 kb were considered as prophages. The smaller elements identified by Phage Finder may be prophage remnants or erroneous assignments. Thus, we identified 2622 large prophages (>18 kb). We found that 50% of the strains were lysogenic (i.e., contained at least one large prophage), consistent with our previous analysis . The characteristics of the bacterial genomes and their prophages are reported in Additional file 2: Table S2.
Amino acid sequences of the integrated reference catalog of the human gut microbiome were retrieved from .
Construction of MuF and Ct_MAD profiles
The HMM profiles PF04233 corresponding to MuF1 and PF06152 corresponding to MuF2 were retrieved from the Pfam 28.0 database  (for details see Additional file 2: Table S3). We built profiles for MuF3 and MuF4. For this, we used N-terminal regions (250 amino acids) of the Gp35 protein from Streptococcus mitis phage SM1 (NC_004996) and of the Gp4 protein from mycobacteriophage Angel (NC_012788) as seed sequences for MuF3 and MuF4, respectively. Seed sequences were used as queries for the Position-Specific Iterative (PSI)-Basic Local Alignment Search Tool (BLAST) algorithm of the BLASTP 2.3.1+ program  that was used to search homologies in the NCBI's non-redundant (nr) protein database with default parameters. Redundancy of the sequence set retrieved from BLASTP was reduced using CD-HIT v4.6  with a 90% identity threshold. The longest representative sequence of each cluster was then aligned with MUSCLE v3.8.31 . HMM profiles were built from the multiple sequence alignments using hmmbuild of HMMER v3.1b1  with default parameters.
To build the HMM profile Ct_MAD, we used the C-terminal extensions of MuF proteins without known domains and clustered them with Single Linkage Clustering of Sequences (SiLiX, sequence identity ≥ 30% and overlap ≥ 50%) . We obtained eight families with more than five sequences. Domain analysis with CDvist tool  showed that four of these families have hits within the C-terminal region of COG2369, which is a cluster of orthologous groups including proteins with a MuF1 domain. The profile Ct_MAD was built from the multiple sequence alignment of the sequences of the four families following the same procedure as for MuF3 and MuF4. MuF3, MuF4, and Ct_MAD HMM profiles are available (Additional files 7, 8, and 9).
HMM profiles of the four MuF domains were compared using the HHsearch program of HH-suite v2.0.15 with default parameters . HMM profiles used for comparison were built with the HHmake program of HH-suite from the same seed alignments previously used to build the profiles with HMMER. Cytoscape v3.4.0 was used to visualize the resulting domain association network (Additional file 1: Figure S1) .
Genetic context of genes encoding a MuF domain in bacteriophages and bacterial chromosomes
A muf gene was regarded as part of a prophage if it was within the boundaries of the large prophage coordinates identified within bacterial chromosomes by Phage Finder (see above). It was considered as part of a remnant prophage if it was located within prophages smaller than 18 kb or if it was near portal and/or terminase genes.
We retrieved all the HMM profiles known to be associated with portal (P), terminase small (TSS), and large subunit (TLS) proteins from Pfam and TIGRFam (Additional file 2: Table S3). We performed a search of these profiles in bacteriophages and bacterial chromosomes using HMMER3 v3.1b1 with the –cut_ga option. We only selected the best e-value profile for each hit and considered, for each of the three categories, i.e., P, TSS, and TSL, the closest gene to the muf gene. Then, we computed the minimal distance between each muf gene and P, TLS, and TSS genes (10 genes around each muf gene). This allowed us to clearly identify the genetic context in the vicinity of muf genes in both bacteriophages and bacterial chromosomes.
Since more than 60% of muf genes were located in +1 of a P gene, we analyzed proteins encoded by genes immediately downstream of P genes (+1) in the phage dataset in order to search for possible additional MuF domains. We did not identify proteins belonging to additional MuF families among the proteins encoded by these genes.
Identification of MuF proteins
We used the HMM profiles of MuF1, MuF2, MuF3, and MuF4 domains to scan our datasets of bacteriophages and bacterial genomes using the hmmsearch program of HMMER v3.1b1 with the –cut_ga option. The results of the detection and the RefSeq accession numbers of all the MuF proteins are reported for phages and bacterial genomes in Additional file 2: Tables S4 and S5, respectively. The list of the MuF proteins retrieved from the integrated reference catalog of the human gut microbiome is provided in Additional file 2: Table S8.
Architecture of MuF proteins
Domain architectures of MuF proteins were analyzed with the CDvist tool  against the Pfam 28.0 domain database  and the Conserved Domain Database (CDD) v3.12 . To determine the minimal length required to define the presence of a C-terminal extension, we searched for the shortest C-terminal extension containing a toxin domain in our dataset. We found a MuF4 toxin with an endoU domain from Bifidobacterium bifidum carrying a C-terminal extension of 113 amino acids. Consequently, we set a cut-off equal to or greater than 100 amino acids to define the presence of a C-terminal extension.
The occurrence of MuF-associated toxin domains in other PT families was investigated using the Conserved Domain Architecture Retrieval Tool (CDART) of NCBI . Cytoscape v3.4.0 was used to visualize the resulting domain association network .
All statistical analyses were performed in R. A Fisher’s exact test (fisher.test in R) was used to determine whether proportions of one variable between two groups are significantly different. A two-sample unpaired Student’s t test (t.test in R) was used to determine statistically significant differences between the mean values on data with two categories.
The authors are grateful to Paulo Tavares for his comments on the manuscript and thank Ogün Adebali for helping with the use of the CDvist tool and Kan Nishida for helping with the Exploratory tool.
The laboratory of Xavier Nassif is supported by La Fondation pour la Recherche Médicale, l’Institut National de la Santé et de la Recherche Médicale, le Centre National de la Recherche Scientifique, and l’Université Paris Descartes. The laboratory of Eduardo Rocha is supported by the Centre National de la Recherche Scientifique and the Institut Pasteur.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its additional files.
AJ, MT, XN, AC, MR, and EPCR conceived and designed the study. AJ and MT performed the experiments and analysis. JAC and BRG contributed to the data analysis. AJ, MT, AC, and EPCR wrote the manuscript, with contributions and comments from all authors. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Jamet A, Nassif X. New players in the toxin field: polymorphic toxin systems in bacteria. MBio. 2015;6(3):e00285–00215.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang D, de Souza RF, Anantharaman V, Iyer LM, Aravind L. Polymorphic toxin systems: comprehensive characterization of trafficking modes, processing, mechanisms of action, immunity and ecology using comparative genomics. Biol Direct. 2012;7:18.View ArticlePubMedPubMed CentralGoogle Scholar
- Poole SJ, Diner EJ, Aoki SK, Braaten BA, t'Kint de Roodenbeke C, Low DA, Hayes CS. Identification of functional toxin/immunity genes linked to contact-dependent growth inhibition (CDI) and rearrangement hotspot (Rhs) systems. PLoS Genet. 2011;7(8):e1002217.View ArticlePubMedPubMed CentralGoogle Scholar
- Ghequire MG, De Mot R. Ribosomally encoded antibacterial proteins and peptides from Pseudomonas. FEMS Microbiol Rev. 2014;38(4):523–68.View ArticlePubMedGoogle Scholar
- Jamet A, Jousset AB, Euphrasie D, Mukorako P, Boucharlat A, Ducousso A, Charbit A, Nassif X. A new family of secreted toxins in pathogenic Neisseria species. PLoS Pathog. 2015;11(1):e1004592.View ArticlePubMedPubMed CentralGoogle Scholar
- Brooks TM, Unterweger D, Bachmann V, Kostiuk B, Pukatzki S. Lytic activity of the Vibrio cholerae type VI secretion toxin VgrG-3 is inhibited by the antitoxin TsaB. J Biol Chem. 2013;288(11):7618–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Aoki SK, Pamma R, Hernday AD, Bickham JE, Braaten BA, Low DA. Contact-dependent inhibition of growth in Escherichia coli. Science. 2005;309(5738):1245–8.View ArticlePubMedGoogle Scholar
- Anderson MS, Garcia EC, Cotter PA. The Burkholderia bcpAIOB genes define unique classes of two-partner secretion and contact dependent growth inhibition systems. PLoS Genet. 2012;8(8):e1002877.View ArticlePubMedPubMed CentralGoogle Scholar
- Durand E, Cambillau C, Cascales E, Journet L. VgrG, Tae, Tle, and beyond: the versatile arsenal of Type VI secretion effectors. Trends Microbiol. 2014;22(9):498–507.View ArticlePubMedGoogle Scholar
- Shneider MM, Buth SA, Ho BT, Basler M, Mekalanos JJ, Leiman PG. PAAR-repeat proteins sharpen and diversify the type VI secretion system spike. Nature. 2013;500(7462):350–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Ma J, Pan Z, Huang J, Sun M, Lu C, Yao H. The Hcp proteins fused with diverse extended-toxin domains represent a novel pattern of antibacterial effectors in type VI secretion systems. Virulence. 2017;Jan 6:1–14. doi:https://doi.org/10.1080/21505594.2017.1279374.
- Pukatzki S, Ma AT, Revel AT, Sturtevant D, Mekalanos JJ. Type VI secretion system translocates a phage tail spike-like protein into target cells where it cross-links actin. Proc Natl Acad Sci U S A. 2007;104(39):15508–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Flaugnatti N, Le TT, Canaan S, Aschtgen MS, Nguyen VS, Blangy S, Kellenberger C, Roussel A, Cambillau C, Cascales E, et al. A phospholipase A1 antibacterial Type VI secretion effector interacts directly with the C-terminal domain of the VgrG spike protein for delivery. Mol Microbiol. 2016;99(6):1099–118.View ArticlePubMedGoogle Scholar
- Zhang D, Iyer LM, Aravind L. A novel immunity system for bacterial nucleic acid degrading toxins and its recruitment in various eukaryotic and DNA viral systems. Nucleic Acids Res. 2011;39(11):4532–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Aoki SK, Poole SJ, Hayes CS, Low DA. Toxin on a stick: modular CDI toxin delivery systems play roles in bacterial competition. Virulence. 2011;2(4):356–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Anderson MS, Garcia EC, Cotter PA. Kind discrimination and competitive exclusion mediated by contact-dependent growth inhibition systems shape biofilm community structure. PLoS Pathog. 2014;10(4):e1004076.View ArticlePubMedPubMed CentralGoogle Scholar
- Garcia EC, Perault AI, Marlatt SA, Cotter PA. Interbacterial signaling via Burkholderia contact-dependent growth inhibition system proteins. Proc Natl Acad Sci U S A. 2016;113(29):8296–301.View ArticlePubMedPubMed CentralGoogle Scholar
- Jiang F, Waterfield NR, Yang J, Yang G, Jin Q. A Pseudomonas aeruginosa type VI secretion phospholipase D effector targets both prokaryotic and eukaryotic cells. Cell Host Microbe. 2014;15(5):600–10.View ArticlePubMedGoogle Scholar
- Holberger LE, Garza-Sanchez F, Lamoureux J, Low DA, Hayes CS. A novel family of toxin/antitoxin proteins in Bacillus species. FEBS Lett. 2012;586(2):132–6.View ArticlePubMedGoogle Scholar
- Koskiniemi S, Lamoureux JG, Nikolakakis KC, t'Kint de Roodenbeke C, Kaplan MD, Low DA, Hayes CS. Rhs proteins from diverse bacteria mediate intercellular competition. Proc Natl Acad Sci U S A. 2013;110(17):7032–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Cho I, Blaser MJ. The human microbiome: at the interface of health and disease. Nat Rev Genet. 2012;13(4):260–70.PubMedPubMed CentralGoogle Scholar
- Woodford N, Livermore DM. Infections caused by Gram-positive bacteria: a review of the global challenge. J Infect. 2009;59 Suppl 1:S4–S16.View ArticlePubMedGoogle Scholar
- Elbaz M, Ben-Yehuda S. Following the fate of bacterial cells experiencing sudden chromosome loss. MBio. 2015;6(3):e00092–00015.View ArticlePubMedPubMed CentralGoogle Scholar
- Cao Z, Casabona MG, Kneuper H, Chalmers JD, Palmer T. The type VII secretion system of Staphylococcus aureus secretes a nuclease toxin that targets competitor bacteria. Nat Microbiol. 2016;2:16183.View ArticlePubMedPubMed CentralGoogle Scholar
- Whitney JC, Peterson SB, Kim J, Pazos M, Verster AJ, Radey MC, Kulasekara HD, Ching MQ, Bullen NP, Bryant D, et al. A broadly distributed toxin family mediates contact-dependent antagonism between gram-positive bacteria. Elife. 2017;6.Google Scholar
- Lyons NA, Kraigher B, Stefanic P, Mandic-Mulec I, Kolter R. A combinatorial kin discrimination system in Bacillus subtilis. Curr Biol. 2016;26(6):733–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Touchon M, Bernheim A, Rocha EP. Genetic and life-history traits associated with the distribution of prophages in bacteria. ISME J. 2016;10(11):2744–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Fokine A, Rossmann MG. Molecular architecture of tailed double-stranded DNA phages. Bacteriophage. 2014;4(1):e28281.View ArticlePubMedPubMed CentralGoogle Scholar
- Harshey RM. The Mu story: how a maverick phage moved the field forward. Mob DNA. 2012;3(1):21.View ArticlePubMedPubMed CentralGoogle Scholar
- Giphart-Gassler M, Wijffelman C, Reeve J. Structural polypeptides and products of late genes of bacteriophage Mu: characterization and functional aspects. J Mol Biol. 1981;145(1):139–63.View ArticlePubMedGoogle Scholar
- Droge A, Santos MA, Stiege AC, Alonso JC, Lurz R, Trautner TA, Tavares P. Shape and DNA packaging activity of bacteriophage SPP1 procapsid: protein components and interactions during assembly. J Mol Biol. 2000;296(1):117–32.View ArticlePubMedGoogle Scholar
- Stiege AC, Isidro A, Droge A, Tavares P. Specific targeting of a DNA-binding protein to the SPP1 procapsid by interaction with the portal oligomer. Mol Microbiol. 2003;49(5):1201–12.View ArticlePubMedGoogle Scholar
- Vinga I, Droge A, Stiege AC, Lurz R, Santos MA, Daugelavicius R, Tavares P. The minor capsid protein gp7 of bacteriophage SPP1 is required for efficient infection of Bacillus subtilis. Mol Microbiol. 2006;61(6):1609–21.View ArticlePubMedGoogle Scholar
- Soding J. Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005;21(7):951–60.View ArticlePubMedGoogle Scholar
- Siboo IR, Bensing BA, Sullam PM. Genomic organization and molecular characterization of SM1, a temperate bacteriophage of Streptococcus mitis. J Bacteriol. 2003;185(23):6968–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Overbeek R, Fonstein M, D'Souza M, Pusch GD, Maltsev N. The use of gene clusters to infer functional coupling. Proc Natl Acad Sci U S A. 1999;96(6):2896–901.View ArticlePubMedPubMed CentralGoogle Scholar
- Huynen M, Snel B, Lathe 3rd W, Bork P. Predicting protein function by genomic context: quantitative evaluation and qualitative inferences. Genome Res. 2000;10(8):1204–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Canchaya C, Desiere F, McShan WM, Ferretti JJ, Parkhill J, Brussow H. Genome analysis of an inducible prophage and prophage remnants integrated in the Streptococcus pyogenes strain SF370. Virology. 2002;302(2):245–58.View ArticlePubMedGoogle Scholar
- Deghorain M, Bobay LM, Smeesters PR, Bousbata S, Vermeersch M, Perez-Morga D, Dreze PA, Rocha EP, Touchon M, Van Melderen L. Characterization of novel phages isolated in coagulase-negative staphylococci reveals evolutionary relationships with Staphylococcus aureus phages. J Bacteriol. 2012;194(21):5829–39.View ArticlePubMedPubMed CentralGoogle Scholar
- Matos RC, Lapaque N, Rigottier-Gois L, Debarbieux L, Meylheuc T, Gonzalez-Zorn B, Repoila F, Lopes Mde F, Serror P. Enterococcus faecalis prophage dynamics and contributions to pathogenic traits. PLoS Genet. 2013;9(6):e1003539.View ArticlePubMedPubMed CentralGoogle Scholar
- Li J, Jia H, Cai X, Zhong H, Feng Q, Sunagawa S, Arumugam M, Kultima JR, Prifti E, Nielsen T, et al. An integrated catalog of reference genes in the human gut microbiome. Nat Biotechnol. 2014;32(8):834–41.View ArticlePubMedGoogle Scholar
- Fieldhouse RJ, Turgeon Z, White D, Merrill AR. Cholera- and anthrax-like toxins are among several new ADP-ribosyltransferases. PLoS Comput Biol. 2010;6(12):e1001029.View ArticlePubMedPubMed CentralGoogle Scholar
- Canchaya C, Proux C, Fournous G, Bruttin A, Brussow H. Prophage genomics. Microbiol Mol Biol Rev. 2003;67(2):238–276, table of contents.Google Scholar
- Guo Y, Quiroga C, Chen Q, McAnulty MJ, Benedik MJ, Wood TK, Wang X. RalR (a DNase) and RalA (a small RNA) form a type I toxin-antitoxin system in Escherichia coli. Nucleic Acids Res. 2014;42(10):6448–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Romero P, Croucher NJ, Hiller NL, Hu FZ, Ehrlich GD, Bentley SD, Garcia E, Mitchell TJ. Comparative genomic analysis of ten Streptococcus pneumoniae temperate bacteriophages. J Bacteriol. 2009;191(15):4854–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Goeders N, Van Melderen L. Toxin-antitoxin systems as multilevel interaction systems. Toxins (Basel). 2014;6(1):304–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Chan WT, Yeo CC, Sadowy E, Espinosa M. Functional validation of putative toxin-antitoxin genes from the Gram-positive pathogen Streptococcus pneumoniae: phd-doc is the fourth bona-fide operon. Front Microbiol. 2014;5:677.PubMedPubMed CentralGoogle Scholar
- Van Melderen L, Saavedra De Bast M. Bacterial toxin-antitoxin systems: more than selfish entities? PLoS Genet. 2009;5(3):e1000437.View ArticlePubMedPubMed CentralGoogle Scholar
- Ghequire MG, Dillen Y, Lambrichts I, Proost P, Wattiez R, De Mot R. Different ancestries of R tailocins in rhizospheric Pseudomonas isolates. Genome Biol Evol. 2015;7(10):2810–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Corda D, Di Girolamo M. Functional aspects of protein mono-ADP-ribosylation. EMBO J. 2003;22(9):1953–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Nikolakakis K, Amber S, Wilbur JS, Diner EJ, Aoki SK, Poole SJ, Tuanyok A, Keim PS, Peacock S, Hayes CS, et al. The toxin/immunity network of Burkholderia pseudomallei contact-dependent growth inhibition (CDI) systems. Mol Microbiol. 2012;84(3):516–29.View ArticlePubMedPubMed CentralGoogle Scholar
- Morse RP, Nikolakakis KC, Willett JL, Gerrick E, Low DA, Hayes CS, Goulding CW. Structural basis of toxicity and immunity in contact-dependent growth inhibition (CDI) systems. Proc Natl Acad Sci U S A. 2012;109(52):21480–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Bossi L, Fuentes JA, Mora G, Figueroa-Bossi N. Prophage contribution to bacterial population dynamics. J Bacteriol. 2003;185(21):6467–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Davies EV, James CE, Kukavica-Ibrulj I, Levesque RC, Brockhurst MA, Winstanley C. Temperate phages enhance pathogen fitness in chronic lung infection. ISME J. 2016;10(10):2553–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Gama JA, Reis AM, Domingues I, Mendes-Soares H, Matos AM, Dionisio F. Temperate bacterial viruses as double-edged swords in bacterial warfare. PLoS One. 2013;8(3):e59043.View ArticlePubMedPubMed CentralGoogle Scholar
- Sassone-Corsi M, Nuccio SP, Liu H, Hernandez D, Vu CT, Takahashi AA, Edwards RA, Raffatellu M. Microcins mediate competition among Enterobacteriaceae in the inflamed gut. Nature. 2016;540(7632):280–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Resource Coordinators NCBI. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2017;45(D1):D12–7.View ArticleGoogle Scholar
- McNair K, Bailey BA, Edwards RA. PHACTS, a computational approach to classifying the lifestyle of phages. Bioinformatics. 2012;28(5):614–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Fouts DE. Phage_Finder: automated identification and classification of prophage regions in complete bacterial genome sequences. Nucleic Acids Res. 2006;34(20):5839–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, Potter SC, Punta M, Qureshi M, Sangrador-Vegas A, et al. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 2016;44(D1):D279–85.View ArticlePubMedGoogle Scholar
- Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28(23):3150–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Mistry J, Finn RD, Eddy SR, Bateman A, Punta M. Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 2013;41(12):e121.View ArticlePubMedPubMed CentralGoogle Scholar
- Miele V, Penel S, Duret L. Ultra-fast sequence clustering from similarity networks with SiLiX. BMC Bioinformatics. 2011;12:116.View ArticlePubMedPubMed CentralGoogle Scholar
- Adebali O, Ortega DR, Zhulin IB. CDvist: a webserver for identification and visualization of conserved domains in protein sequences. Bioinformatics. 2015;31(9):1475–7.View ArticlePubMedGoogle Scholar
- Remmert M, Biegert A, Hauser A, Soding J. HHblits: lightning-fast iterative protein sequence searching by HMM-HMM alignment. Nat Methods. 2011;9(2):173–5.View ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.View ArticlePubMedPubMed CentralGoogle Scholar
- Geer LY, Domrachev M, Lipman DJ, Bryant SH. CDART: protein homology by domain architecture. Genome Res. 2002;12(10):1619–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Sampson T, Broussard GW, Marinelli LJ, Jacobs-Sera D, Ray M, Ko CC, Russell D, Hendrix RW, Hatfull GF. Mycobacteriophages BPs, Angel and Halo: comparative genomics reveals a novel class of ultra-small mobile genetic elements. Microbiology. 2009;155(Pt 9):2962–77.View ArticlePubMedPubMed CentralGoogle Scholar