Vibrio cholerae is a Gram-negative bacterial pathogen that consists of over 200 serogroups with differing pathogenic potential. Only strains that express the virulence factors cholera toxin (CT) and toxin-coregulated pilus (TCP) are capable of pandemic spread of cholera diarrhoea. Regardless, all V. cholerae strains sequenced to date harbour genes for the type VI secretion system (T6SS) that translocates effectors into neighbouring eukaryotic and prokaryotic cells. Here we report that the effectors encoded within these conserved gene clusters differ widely among V. cholerae strains, and that immunity proteins encoded immediately downstream from the effector genes protect their host from neighbouring bacteria producing corresponding effectors. As a consequence, strains with matching effector-immunity gene sets can coexist, while strains with different sets compete against each other. Thus, the V. cholerae T6SS contributes to the competitive behaviour of this species.
Vibrio cholerae is a Gram-negative bacterium predominantly found in marine environments1. V. cholerae comprises more than 200 highly diverse serogroups, but only toxigenic O1 serogroup strains that encode the virulence factors cholera toxin (CT) and toxin-coregulated pilus (TCP) have caused pandemic spread of the diarrhoeal disease cholera2,3,4. TCP biosynthesis genes are encoded within Vibrio pathogenicity island-1 (VPI-1), and TCP protruding from the bacterium serves as a receptor for a filamentous bacteriophage (CTX-ф) that transduces CT genes into the chromosomes of pandemic V. cholerae strains5. The contribution of CTX-ф is crucial to pathogenesis as only V. cholerae strains encoding CT are capable of pandemic spread5. When humans consume contaminated water, toxigenic V. cholerae passes through the gastric acid barrier, expresses TCP, secretes CT and colonizes the small intestine. Next, V. cholerae multiplies rapidly and exits the human host during diarrhoeal purges6. If left untreated, subsequent dehydration of the host can lead to death. During the presently ongoing seventh pandemic, approximately five million infections and 120,000 deaths due to cholera are estimated to occur each year7. Rehydration therapy is the preferred treatment, whereas the use of antibiotics is limited because of the emergence of antibiotic-resistant strains8.
In addition to the principal virulence factors TCP and CT, V. cholerae genomes encode the type VI secretion system (T6SS), a contact-dependent toxin delivery system9,10. Genes encoding the T6SS are present in ~25% of all sequenced Gram-negative bacteria, and are found in all sequenced V. cholerae strains11. The T6SS is a molecular puncturing device employed by bacteria to translocate effectors into eukaryotic and prokaryotic cells10,12. The hallmark of the secretion system is an inner tube composed of haemolysin-coregulated protein (Hcp) polymers13. Similar to a mechanism found in tailed bacteriophages, contraction of an outer sheath around the inner Hcp tube causes the inner tube to be ejected14,15. Three V. cholerae T6SS effectors (the lipase TseL16,17, the membrane-disrupting VasX18 and the peptidoglycan-degrading VgrG-3 (refs 16, 19)) believed to be situated at the distal end of the ejected Hcp tube target prokaryotic cells. On T6SS-mediated effector transfer from one bacterium into another, TseL, VasX and VgrG-3 are inactivated in target cells that express corresponding immunity genes tsiV1, tsiV2 and tsiV3, respectively. The immunity genes are encoded directly downstream of their corresponding effector gene as part of an effector-immunity pair in each of the three T6SS gene clusters, the large and the two small auxiliary gene clusters (Fig. 1 and refs 16, 19).
We first noticed that V. cholerae strains compete against each other in a T6SS-dependent manner when we discovered that V. cholerae strains from the Rio Grande river delta (Southwestern United States) are able to kill each other20; however, the molecular mechanism and the biological significance of this intraspecies competition remained unclear. Here we apply a bioinformatics approach to identify the molecular basis of T6SS-mediated (in)compatibility. This analysis allows us to subdivide this bacterial species into strains that are able to coexist (compatible strains) and strains that outcompete each other in a T6SS-dependent manner (incompatible strains). Phylogenetic analysis reveals that toxigenic (as defined by the presence of VPI-1 and CTX-Φ) V. cholerae strains cluster within a single group of compatible strains.
Conserved T6SS clusters harbour diverse effector modules
We hypothesized that a given strain of V. cholerae deploys effectors to kill strains lacking cognate immunity proteins. Therefore, compatible strains do not kill each other because they express the same immunity proteins, whereas incompatible strains are unable to coexist because they kill each other as a result of having different effector and immunity gene pairs21. To test this hypothesis, the nucleotide sequences of the three T6SS clusters of 37 V. cholerae strains collected from a variety of sources over the last 77 years (Supplementary Table 1) were aligned and compared. As shown in Fig. 1, auxiliary clusters 1, 2 and the large cluster were found to be highly conserved among the available V. cholerae genomes except for regions of high diversity and low GC content (shaded blue) that encode the bacterial effectors TseL, VasX and VgrG-3, and their immunity proteins TsiV1, TsiV2 and TsiV3, respectively. Conserved T6SS genes contained an average GC content of ~51%, whereas the GC content in the region coding the effector-immunity pairs was an average of ~9% lower (Fig. 1), suggesting that these have a different origin than the remainder of the T6SS gene clusters and may have been acquired independently. We named these regions ‘effector modules’, each encoding a T6SS effector and its cognate immunity protein. The combination of all three effector modules in one strain comprises an effector module set.
T6SS effector modules determine compatibility
Effector modules at each of the three loci in the V. cholerae chromosomes were grouped into 15 families based on the amino-acid sequence identity of their immunity proteins TsiV1, TsiV2 and TsiV3; immunity proteins with >30% identity were grouped as a family because proteins sharing more than 30% identity are predicted to fold similarly21. As shown in Fig. 2, distinct families of effector modules at each locus were labelled alphabetically: three families were found in auxiliary cluster 1 (A–C), five in auxiliary cluster 2 (A–E) and seven in the large cluster (A–G), allowing the hypothetical existence of 105 different effector module sets. This is a conservative estimate, because our analysis does not account for the existence of additional, yet unknown effector modules.
Effector module families also comprise multiple subfamilies, each consisting of identical immunity proteins and effectors that can display varying degrees of polymorphisms. Subfamilies are indicated by a subscript (for example, A2) as shown in Fig. 3 and Supplementary Fig. 1. Effector module sets of V. cholerae strains were labelled with a three-letter code (for example, A1A1A1) representing distinct families/subfamilies in T6SS auxiliary cluster 1, auxiliary cluster 2 and large cluster, respectively. Thus, we hypothesized that immunity proteins form the basis of (in)compatibility, allowing strains with polymorphic effectors to be compatible as long as these effectors are paired with identical immunity proteins.
To test whether the effector module sets determine compatibility between strains, we chose V52 and C6706 as A1A1A1 module set representatives. Under laboratory conditions, V52 expresses T6SS effector and immunity genes, whereas C6706 expresses only immunity genes16,22. Therefore, although C6706 is T6SS-effector silent, C6706 exhibits immunity towards V52. An immunity deletion mutant (ΔtsiV1, ΔtsiV2 and ΔtsiV3) of T6SS-silent C6706 is viable under standard laboratory growth conditions, allowing us to study the genetic basis of compatibility between strains with the same effector module set. C6706 strains encoding one, none or all immunity genes (tsiV1, tsiV2 and tsiV3) were challenged with V52 that express one, none or all T6SS effectors (tseL, vasX or vgrG-3). Each V. cholerae V52 mutant expressing any one of the three T6SS effectors killed a C6706 immunity mutant lacking all three immunity genes 10- to 100-fold (Fig. 4, triangles). In contrast, C6706 mutants encoding only one immunity gene survived in the presence of V52 equipped with only the single corresponding T6SS effector (Fig. 4, arrows). These results indicate that cognate effector and immunity gene pairs function independently of one another, and that all three immunity proteins are required for complete protection. V. cholerae strains that encode the same effector module sets are able to coexist and are compatible. Strains that lack cognate immunity genes are outcompeted and are thus incompatible.
Strains with different effector module sets are incompatible
To challenge our compatibility model, we coincubated V52 (A1A1A1 effector module set) with T6SS-active strains (Supplementary Figs 2A and 3) harbouring the following effector module sets: A1A1A1 (C6706), A3A1B1 (MZO-2), C5D3A1 (V51), C1E1A1 (MZO-3), C4E1C2 (DL4215), C3D2C1 (1587), C1D4F1 (AM-19226) or C6E2E1 (DL4211) (Fig. 5a). In agreement with our hypothesis, all A1A1A1 strains tested were able to coexist (Fig. 5b), whereas non-A1A1A1 V. cholerae were effectively outcompeted by wild-type V52 when inoculated at a 1:1 ratio (Fig. 5a). Control experiments carried out with a V52 vasK in-frame deletion mutant unable to translocate effectors confirmed that killing was T6SS dependent (Fig. 5a,b).
To further examine the competitive interactions between A1A1A1 V. cholerae and strains that harbour a different effector module set, V52 (A1A1A1) and AM-19226 (C1D4F1) were inoculated at different ratios. Both strains exhibit constitutive T6SS activity under laboratory conditions as shown by secretion of the inner tube component Hcp and T6SS-dependent killing of Escherichia coli (Supplementary Figs 2B,C and 4). The two strains were coincubated at equal numbers on agar plates, and enumeration of survivors after 4 h showed that wild-type V52 outcompeted wild-type AM-19226 in a T6SS-dependent manner (Fig. 6a). When mixed at different input ratios, V52 remained competitive even when inoculated in the minority (Fig. 6b). In conclusion, V. cholerae strains with different effector module sets compete with each other and are thus incompatible.
Toxigenic strains carry compatible T6SS effector module sets
To determine the genetic relationship between compatible and incompatible strains, a phylogenetic tree for the 37 V. cholerae strains (Supplementary Table 1) used in this study was built based on six polymorphic housekeeping genes23 (Fig. 7). TCP was found exclusively in AAA strains. CTX-Φ was present in all but two AAA strains (2740-80, M66-2, both belonging to the O1 serogroup) and in one non-AAA strain, V51. Except V51—a clinical O141 serotype strain from the United States, all non-AAA strains analysed were VPI-1- and CTX-Φ-negative. Toxigenic strains (including C6706, MO10, V52, N16961 and O395) have the AAA effector module set, while V. cholerae strains with any other effector module set belong to phylogenetic lineages that exclude toxigenic strains.
This report describes a diverse set of genetic elements within the T6SS gene clusters of V. cholerae strains that we named effector modules. These modules differ in their GC content from the remainder of the clusters and encode T6SS effector and immunity genes. The three effector modules comprise an effector module set characteristic for each V. cholerae strain. Our data suggest that V. cholerae strains with the same effector module set can coexist and are therefore compatible, whereas when strains with different effector module sets engage their T6SS, one or both will experience cell death resulting in incompatibility. Compatibility between strains based on effector modules of the same family needs to be further explored as immunity proteins can differ within a family. We speculate that compatibility between strains decreases as the diversity between their immunity proteins increases. Conversely, incompatible strains might coexist under conditions that do not support production of all T6SS effectors, or as a result of cross-immunity24 provided by homologous immunity proteins.
Our study identified a wide variety of effector modules encoded within V. cholerae T6SS gene clusters. The differences in GC content between genes of the effector modules and genes that encode structural T6SS components suggest independent acquisition and integration events into the genome. This hypothesis is also supported by phylogenetic evidence (Fig. 7). For example, despite having undergone 162 genetic changes in the housekeeping genes since diverging from their last common ancestor (Fig. 7, Supplementary Fig. 5A), strains V51 and VL426 carry an identical copy of the immunity gene tsiV3A1 (the family of the effector module that encodes tsiV3 is indicated in superscript; Supplementary Fig. 5B); sequence identity indicates that V51 and VL426 recently acquired the tsiV3A1 allele horizontally. Similarly, strains MZO-2 and 12129(1) are distantly related to O1 strains (Fig. 7), but carry an identical tsiV2A1 allele (Supplementary Fig. 6), strongly supporting the notion that tsiV2A1 was recently acquired by horizontal gene transfer. Evidence for these and other immunity gene transfer events can be found for any of the three effector modules, suggesting that horizontal acquisition of effector modules has been a recurrent theme in V. cholerae evolution. Furthermore, a phylogenetic tree constructed from six housekeeping genes (Fig. 7) revealed that strains derived from the same recent common ancestors frequently carry distinct effector module sets, suggesting that horizontal movement of effector modules, at least in some cases, occurred sometime after differentiation into distinct strains. For example, DL4215 (C4E1C2) and TMA21 (C2B1F2) share a recent common ancestor but differ in two of the three effector modules. Analogously, AAA strains share a common ancestor with MZO-3 (C1E1A1), yet differ in two effector modules. In these lineages, the effector module sets of closely related strains are not a result of recent common ancestry, but probably a result of subsequent horizontal effector module acquisition. In these cases, the divergence of newly evolving strains from a common ancestor might have been accelerated by the acquisition of different effector modules. T6SS-mediated competition interferes with cell–cell contact and therefore possibly affects contact-dependent exchange of genetic elements. Subsequently, such events may be restricted to bacteria that encode the same effector module sets and represent one way for a newly evolving strain to separate from its ancestor.
The cluster of strains with the AAA effector module sets (blue-shaded portion of Fig. 7) consists of strains that differ in their VPI-1/CTX content (+/+, +/− and −/−) suggesting that toxigenic and by extension pandemic V. cholerae acquired VPI-1 (and in some cases CTX-Φ) after AAA effector modules. As shown in Fig. 7, this holds true regardless of classical (responsible for the sixth pandemic and probably for the first five pandemics25,26,27) or El Tor biotype (responsible for the current seventh pandemic28). This indicates that AAA acquisition also preceded biotype differentiation and that toxigenic strains encode the AAA effector module set as a result of common ancestry. The conservation of the AAA effector module set over the evolutionary distance separating classical and El Tor biotypes suggests that natural selection could have pressured toxigenic strains to retain the contributions of the AAA module set because it is advantageous to human colonization (as they are likely to come in close contact with V. cholerae employing other effectors in the small intestine). Other strains such as MZO-2 (A3A1B1) and 623-39 (C5D1C1) separated by similar evolutionary distances have experienced changes in all three effector modules (Fig. 7). These observations, including that none of the toxigenic strains included in this study have acquired non-A modules, suggest that natural selection favours the AAA effector module set in toxigenic strains. T6SS-mediated competition of V. cholerae occurs in vivo29; thus, intraspecies competition might allow AAA strains to outcompete incompatible V. cholerae strains on host entry to establish a colonization niche. In such a scenario, AAA strains similar to 2740-80 lacking TCP and CT might be viewed as cheaters30,31,32 because when gaining access to a host along with a toxigenic AAA strain, they could coexist and colonize successfully without producing virulence factors.
The discovery of diverse effector modules has translational implications because inhibition of T6SS immunity would not affect bacterial viability per se, and be detrimental only when the target bacterium comes in close contact with a T6SS-active competitor. Current antibiotic compounds target bacterial metabolism, cell wall and growth, applying selective pressure towards essential functions, thereby favouring the emergence of adaptive mechanisms8. As a result, mutations that give rise to antibiotic resistance give rise to the emergence of strains that complicate the treatment of bacterial infections. We suggest that antibacterial strategies against a contact-dependent killing mechanism conserved among toxigenic strains would apply less selective pressure than currently employed antibiotics that target essential cellular housekeeping functions. AAA immunity genes might represent particularly attractive targets because all toxigenic (TCP+/CT+) strains capable of causing pandemic cholera that we evaluated in our study carry the AAA effector module set. Finally, AAA effector module sets could be targeted as a biomarker for toxigenic strain identification for epidemic risk assessment. Classification of the V. cholerae strains based on effector module sets may affect future outbreak management, diagnostics and therapeutic approaches.
Bacterial strains and cell culture
V. cholerae was grown in Luria Bertani (LB) broth supplemented with either rifampicin, 50 μg ml−1; ampicillin, 100 μg ml−1; or streptomycin, 100 μg ml−1 as appropriate. An endonuclease-minus mutant of AM-19226 was used as a wild type33 and as the parent strain for the vasK deletion mutant. A complete list of strains used in this study and their sources is shown in Supplementary Table 2.
This assay was performed as previously described10. In brief, predator and prey bacteria were mixed at a 1:1 ratio unless otherwise stated. A total number of 2 × 108 bacteria were spotted on LB plates and incubated at 37 °C for 4 h. Bacteria were resuspended, diluted and plated on LB plates containing antibiotics that allowed selection of predator or prey. Each experiment was done twice in duplicates. Colony-forming units per ml of V52 were divided by the colony-forming units per ml of the competing strain at time points t=0 h and t=4 h; the V52/competing strain ratio at t=4 h was divided by the V52/competing strain ratio at t=0 h to calculate the competitive index.
In-frame deletion of effector genes, immunity genes and vasK was performed using pWM91-based knockout constructs34. A list of primer sequences used is shown in Supplementary Table 3.
Analysed genome sequences not generated by in-house Illumina NGS sequencing were obtained from NCBI or The Broad Institute databases (Supplementary Table 4). The GC content of V. cholerae V52 was analysed using Geneious (version 6.1.5). The T6SS gene clusters of the 37 V. cholerae strains were compared using the MAUVE alignment tool35 and visualized using Geneious. Prediction of effector functions was based on Phyre2 homology searches36. To group the strains by effector modules (Fig. 2), the amino-acid sequences of the immunity proteins were aligned using MUSCLE37. Immunity proteins with a percentage amino-acid identity of at least 30% were grouped into the same family21. Amino-acid sequences of immunity proteins within a subfamily are identical. The phylogenetic tree shown in Fig. 7 was assembled based on the concatenated nucleotide sequences of adk (adenylate kinase, locus tag VC0986), gyrB (DNA gyrase subunit, locus tag VC0015), mdh (malate dehydrogenase, locus tag VC0432), recA (recombinase A, locus tag VC0543), pgi (glucose-6-phosphate isomerase, locus tag VC0374) and rpoB (DNA-dependent RNA polymerase, locus tag VC0328) to establish the relatedness of V. cholerae to the corresponding genes from V. mimicus as the outgroup23. The data set containing these homologues had 38 sequences and 11,046 positions. Alignments were done using ClustalW38. All positions were determined to be homologous and were included in the alignment. jModelTest2 (refs 39, 40) was used to find the best model of nucleotide evolution for the sequences, incorporating corrections for invariable sites as well as a four-category gamma correction for rate variation when appropriate. For the tree shown in Fig. 7, a substitution model that accounted for both invariant sites and gamma-distributed rate heterogeneity (GTR+I+G) was chosen.
PhyML v. 2.4.4 (ref. 40) was used for maximum likelihood analysis, and to generate ML bootstrap values based on 100 pseudoreplicates of each data set. The phylogenetic tree shown in Fig. 7 is the PhyML consensus topology.
To determine the phylogeny of the tsiVA nucleotide sequences in Supplementary Fig. 6, the nucleotide sequences of the indicated genes were aligned using MUSCLE. RAxML41 (boostrap=100) was used for the phylogenetic analysis.
Protein secretion profile
Western blot samples were prepared and analysed as described previously20. In brief, bacteria were grown to mid-logarithmic phase, centrifuged, and supernatants and pellets were subject to SDS–PAGE followed by immunodetection with polyclonal rabbit anti-Hcp antiserum42 (diluted 1:500) and monoclonal mouse anti-DnaK antibody (Stressgen, diluted 1:15,000). Goat anti-mouse horseradish peroxidase and goat anti-rabbit horseradish peroxidase (both Santa Cruz, diluted 1:3,000) were used as secondary antibodies.
The two-tailed Student’s t-test was used to determine statistical significance between indicated groups of data points. Before performing the t-test, variances within two groups were tested for equality using the F-test. In case of significantly different variances, the logarithmic values of the data points were used for the t-test. Additional information to the individual statistical analysis performed in this study can be found in Supplementary Table 5.
Accession codes: Nucleotide sequences of V. cholerae T6SS gene clusters have been deposited in the GenBank nucleotide database under accession codes KF228943, KF228947, KF228950, KC955251, KF228941, KF228945, KF228942, KF228946, KF228949, KF228944, KF228948 and KF228951.
How to cite this article: Unterweger, D. et al. The Vibrio cholerae type VI secretion system employs diverse effector modules for intraspecific competition. Nat. Commun. 5:3549 doi: 10.1038/ncomms4549 (2014).
Kaper, J., Lockman, H., Colwell, R. R. & Joseph, S. W. Ecology, serology, and enterotoxin production of Vibrio cholerae in Chesapeake Bay. Appl. Environ. Microbiol. 37, 91–103 (1979).
Holmgren, J. Actions of cholera toxin and the prevention and treatment of cholera. Nature 292, 413–417 (1981).
Taylor, R. K., Miller, V. L., Furlong, D. B. & Mekalanos, J. J. Use of phoA gene fusions to identify a pilus colonization factor coordinately regulated with cholera toxin. Proc. Natl Acad. Sci. USA 84, 2833–2837 (1987).
Herrington, D. A. et al. Toxin, toxin-coregulated pili, and the toxR regulon are essential for Vibrio cholerae pathogenesis in humans. J. Exp. Med. 168, 1487–1492 (1988).
Waldor, M. K. & Mekalanos, J. J. Lysogenic conversion by a filamentous phage encoding cholera toxin. Science 272, 1910–1914 (1996).
Sack, D. A., Sack, R. B., Nair, G. B. & Siddique, A. K. Cholera. Lancet 363, 223–233 (2004).
WHO. Cholera. Fact Sheet N° 107, (2012).
Kitaoka, M., Miyata, S. T., Unterweger, D. & Pukatzki, S. Antibiotic resistance mechanisms of Vibrio cholerae. J. Med. Microbiol. 60, 397–407 (2011).
Pukatzki, S. et al. Identification of a conserved bacterial protein secretion system in Vibrio cholerae using the Dictyostelium host model system. Proc. Natl Acad. Sci. USA 103, 1528–1533 (2006).
MacIntyre, D. L., Miyata, S. T., Kitaoka, M. & Pukatzki, S. The Vibrio cholerae type VI secretion system displays antimicrobial properties. Proc. Natl Acad. Sci. USA 107, 19520–19524 (2010).
Boyer, F., Fichant, G., Berthod, J., Vandenbrouck, Y. & Attree, I. Dissecting the bacterial type VI secretion system by a genome wide in silico analysis: what can be learned from available microbial genomic resources? BMC Genomics 10, 104 (2009).
Pukatzki, S., Ma, A. T., Revel, A. T., Sturtevant, D. & Mekalanos, J. J. Type VI secretion system translocates a phage tail spike-like protein into target cells where it cross-links actin. Proc. Natl Acad. Sci. USA 104, 15508–15513 (2007).
Mougous, J. D. et al. A virulence locus of Pseudomonas aeruginosa encodes a protein secretion apparatus. Science 312, 1526–1530 (2006).
Leiman, P. G. et al. Type VI secretion apparatus and phage tail-associated protein complexes share a common evolutionary origin. Proc. Natl Acad. Sci. USA 106, 4154–4159 (2009).
Basler, M., Pilhofer, M., Henderson, G. P., Jensen, G. J. & Mekalanos, J. J. Type VI secretion requires a dynamic contractile phage tail-like structure. Nature 483, 182–186 (2012).
Dong, T. G., Ho, B. T., Yoder-Himes, D. R. & Mekalanos, J. J. Identification of T6SS-dependent effector and immunity proteins by Tn-seq in Vibrio cholerae. Proc. Natl Acad. Sci. USA 110, 2623–2628 (2013).
Russell, A. B. et al. Diverse type VI secretion phospholipases are functionally plastic antibacterial effectors. Nature 496, 508–512 (2013).
Miyata, S. T., Kitaoka, M., Brooks, T. M., McAuley, S. B. & Pukatzki, S. Vibrio cholerae requires the type VI secretion system virulence factor VasX to kill Dictyostelium discoideum. Infect. Immun. 79, 2941–2949 (2011).
Brooks, T. M., 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. 288, 7618–7625 (2013).
Unterweger, D. et al. Constitutive type VI secretion system expression gives Vibrio cholerae intra- and interspecific competitive advantages. PLoS One 7, e48320 (2012).
Rost, B. Twilight zone of protein sequence alignments. Protein Eng. 12, 85–94 (1999).
Miyata, S. T., Unterweger, D., Rudko, S. P. & Pukatzki, S. Dual expression profile of type VI secretioin system immunity genes protects pandemic Vibrio cholerae. PLoS Pathog. 9, e1003752 (2013).
Boucher, Y. et al. Local mobile gene pools rapidly cross species boundaries to create endemicity within global Vibrio cholerae populations. mBio 2, e00335–10 (2011).
Russell, A. B. et al. A widespread bacterial type VI secretion effector superfamily identified using a heuristic approach. Cell Host Microbe 11, 538–549 (2012).
Colwell, R. R. Global climate and infectious disease: the cholera paradigm. Science 274, 2025–2031 (1996).
Safa, A., Nair, G. B. & Kong, R. Y. Evolution of new variants of Vibrio cholerae O1. Trends Microbiol. 18, 46–54 (2010).
Devault, A. M. et al. Second-pandemic strain of Vibrio cholerae from the Philadelphia cholera outbreak of 1849. New Engl. J. Med. 370, 334–340 (2014).
Cvjetanovic, B. & Barua, D. The seventh pandemic of cholera. Nature 239, 137–138 (1972).
Fu, Y., Waldor, M. K. & Mekalanos, J. J. Tn-Seq analysis of Vibrio cholerae intestinal colonization reveals a role for T6SS-mediated antibacterial activity in the host. Cell Host Microbe 14, 652–663 (2013).
West, S. A., Diggle, S. P., Buckling, A., Gardner, A. & Griffins, A. S. The social lives of microbes. Annu. Rev. Ecol. Evol. Syst. 38, 53–77 (2007).
Smith, J. The social evolution of bacterial pathogenesis. Proc. Biol. Sci. 268, 61–69 (2001).
Ennis, H. L., Dao, D. N., Pukatzki, S. U. & Kessin, R. H. Dictyostelium amoebae lacking an F-box protein form spares rather than stalk in chimeras with wild type. Proc. Natl Acad. Sci. USA 97, 3292–3297 (2000).
Tam, V. C., Serruto, D., Dziejman, M., Brieher, W. & Mekalanos, J. J. A type III secretion system in Vibrio cholerae translocates a formin/spire hybrid-like actin nucleator to promote intestinal colonization. Cell Host Microbe 1, 95–107 (2007).
Metcalf, W. W. et al. Conditionally replicative and conjugative plasmids carrying lacZα for cloning, mutagenesis, and allele replacement in bacteria. Plasmid 35, 1–13 (1996).
Darling, A. C., Mau, B., Blattner, F. R. & Perna, N. T. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 14, 1394–1403 (2004).
Kelley, L. A. & Sternberg, M. J. Protein structure prediction on the Web: a case study using the Phyre server. Nat. Protoc. 4, 363–371 (2009).
Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797 (2004).
Larkin, M. A. et al. Clustal W and Clustal X version 2.0. Bioinformatics 23, 2947–2948 (2007).
Darriba, D., Taboada, G. L., Doallo, R. & Posada, D. jModelTest 2: more models, new heuristics and parallel computing. Nat. Methods 9, 772 (2012).
Guindon, S. & Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 52, 696–704 (2003).
Stamatakis, A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22, 2688–2690 (2006).
Ma, A. T., McAuley, S., Pukatzki, S. & Mekalanos, J. J. Translocation of a Vibrio cholerae type VI secretion effector requires bacterial endocytosis by host cells. Cell Host Microbe 5, 234–243 (2009).
We thank Yan Boucher, Fabini Orata, Joel Dacks, Emily Herman and Bart Hazes (University of Alberta) for helpful discussions and Marcia Craig for critically reading the manuscript. We thank John Mekalanos (Harvard Medical School), Michelle Dziejman (University of Rochester), Jun (Jay) Zhu (University of Pennsylvania) and Shelley Payne (The University of Texas at Austin) for providing strains. Work in the laboratory of S.P. is supported by the Canadian Institute of Health Research Operating Grant MOP-84473 and by Alberta Innovates Health Solutions. D.U. and S.T.M. are supported by Graduate Student Scholarships from Alberta Innovates Health Solutions. D.P. was supported by NIH Grants MD001091-01 and GM068855-02.
The authors declare no competing financial interests.
Supplementary Figures 1-6, Supplementary Tables 1-5 and Supplementary References (PDF 1283 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 3.0 Unported License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/3.0/
About this article
Cite this article
Unterweger, D., Miyata, S., Bachmann, V. et al. The Vibrio cholerae type VI secretion system employs diverse effector modules for intraspecific competition. Nat Commun 5, 3549 (2014). https://doi.org/10.1038/ncomms4549
This article is cited by
Bacterial type VI secretion system (T6SS): an evolved molecular weapon with diverse functionality
Biotechnology Letters (2023)
Single nucleotide polymorphism determines constitutive versus inducible type VI secretion in Vibrio cholerae
The ISME Journal (2022)
Dynamic Subspecies Population Structure of Vibrio cholerae in Dhaka, Bangladesh
Microbial Ecology (2022)
Type VI secretion system mutations reduced competitive fitness of classical Vibrio cholerae biotype
Nature Communications (2021)
Glucose confers protection to Escherichia coli against contact killing by Vibrio cholerae
Scientific Reports (2021)
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.