M. xanthus inter-strain antagonism related to the presence of Mx-alpha was originally observed as ‘swarm inhibition’, during which a nonmotile ancestor strain (Mx-alpha⁺) inhibited the outward swarming of a motile strain (missing two of three Mx-alpha units) during co-culture on agar. This phenomenon is traA-dependent and therefore is likely an outcome of OME (Pathak et al., 2013) (see Figure 1). Swarm inhibition was further demonstrated to be caused by cell death of the motile strain (Dey et al., 2016). We first sought to identify the specific genetic determinant on Mx-alpha that was required for antagonism and cell death of the susceptible motile strain. Sequence analysis of the two Mx-alpha units retained in the ancestral DK101 strain (a.k.a. DZF1) (Müller et al., 2013) but lost in DK1622 revealed a candidate toxin gene (MXF1DRAFT_07513), which we have designated sitA1 for swarm inhibition toxin. This gene encodes a predicted lipoprotein that contains a C-terminal nuclease domain with a WHH motif (Zhang et al., 2011) (Figure 1A). The absence of lysine at the +2 position or alanine at the +7 position relative to the +1 cysteine in the N-terminal lipobox suggests that this protein is localized to the OM (Bhat et al., 2011) (see Supplementary file 1). Given that OME efficiently transfers OM lipoproteins between cells (Wei et al., 2011), this open reading frame (ORF) represents a promising candidate for the antagonistic determinant. Immediately downstream of sitA1 is a gene (sitI1) that shows homology to the SUKH-family of immunity proteins commonly found in polymorphic toxin systems (Zhang et al., 2011) (Figure 1A). Upstream of sitA1 is a hypothetical gene (sitB1) of unknown function. The sitB1 ORF overlaps with sitA1 by 11 base pairs, suggesting that the genes form an operon and function together (Figure 1A).

Figure 1

Download asset Open asset

To test if sitA1 is the swarm inhibition toxin, we used the swarm inhibition assay as a readout for the contribution of sitA or sitI toward cell death of the susceptible motile strain during co-culture with the antagonistic, nonmotile ancestral strain. In this assay, cell death of the susceptible strain results in no cells visibly escaping the mixed culture spot. In contrast, abrogation of cell death results in the appearance of the motile strain moving outward from the colony co-culture. As shown previously (Pathak et al., 2013; Dey et al., 2016; Dey and Wall, 2014), nonmotile ancestral cells inhibited the motility of DK1622, but not the DK1622 ∆traA strain (Figure 1B, rows 1–4). Importantly, nonmotile ancestors carrying a sitA1 mutation did not inhibit DK1622 (Figure 1B, row 5). Further, expression of sitI1 in motile DK1622 cells also prevented antagonism (Figure 1B, row 6), consistent with the prediction that sitI encodes an immunity protein that neutralizes SitA1. To test whether sitBAI1 is sufficient to convert non-antagonistic cells into killers, we expressed the gene cassette in a DK1622-derived nonmotile strain, which lacks two Mx-alpha units (see Figure 2A), and does not cause swarm inhibition. As predicted, ectopic sitBAI1 expression allowed the nonmotile ∆Mx-alpha cells to inhibit DK1622, thus recapitulating the antagonistic phenotype exhibited by the nonmotile ancestor strain (Figure 1B, rows 7–8). These combined results suggest that sitBAI1 may function as a toxin/immunity system responsible for the antagonistic behaviors previously described (Pathak et al., 2013; Dey et al., 2016; Dey and Wall, 2014). Therefore, the loss of two Mx-alpha units, and thus the sitBAI1 operon (Figure 2A), during the construction of DK1622 from an ancestral DK101 strain explains why the latter strain antagonizes the former.

Figure 2 with 3 supplements see all

Download asset Open asset

The three tandem Mx-alpha units in the ancestor strain are related and contain different alleles of >80 genes. This region therefore represents a rare bacterial polyploid element – that is, it contains three Mx-alpha prophage genomes with divergent gene allele sets (Dey et al., 2016). Inspection of the other two Mx-alpha units revealed additional putative sitBAI operons. The second Mx-alpha unit (absent from DK1622) carries sitA2 (MXF1DRAFT_07313), which encodes a putative lipoprotein with clear homology to the N-terminal region of SitA1, though the C-terminal domains are unrelated (Figure 2A). The sitA2 gene is flanked by a sitB1 homolog, sitB2, and a downstream candidate immunity gene, sitI2 (Figure 2A). The third Mx-alpha unit, which is shared between the ancestral strain and DK1622, also appears to contain a sitBAI operon. Although the sitA3 gene (MXF1DRAFT_05864 or MXAN_1899) has low sequence homology with sitA1 and sitA2, the three genes nevertheless share several key features: (1) sitA3 is preceded by sitB3, which is homologous to sitB1 and sitB2, (2) sitA3 occupies a similar position within its Mx-alpha unit as the other sitA genes, (3) sitA3 encodes an OM lipoprotein signal sequence, (4) sitA3 encodes a predicted C-terminal tRNase toxin domain, and (5) sitA3 is adjacent to a downstream putative immunity gene, sitI3 (Figure 2A). This analysis suggests that the three Mx-alpha units each contain distinct sitBAI toxin/immunity operons.

To determine whether SitA lipoproteins function as toxins, we expressed each sitBAI cassette in DK1622 and tested the competitive fitness of the resulting inhibitor strains against parental DK1622 target cells that lack the corresponding sitI gene. Target strains were labeled with fluorescent markers and co-cultured with inhibitor strains on agar for 24 hr. Competition outcomes were assessed by competitive index, which is the ratio of target cells to toxin-producing inhibitor cells at 24 hr relative to the starting ratio. For example, a competitive index of 0.01 indicates that the ratio of target cells to inhibitor cells decreased 100-fold. In all instances, sitBAI-expressing inhibitor cells significantly outcompeted target cells, whereas the mock-inhibitor did not (Figure 2B). Delivery of SitA1 and SitA2 over a 24 hr period induced filamentation and lysis of target cells, whereas SitA3 induced rounding and lysis of target cells (Figure 2—figure supplement 1). Furthermore, a ∆traA mutation in either the target or inhibitor strain abolished the inhibitor’s competitive advantage (Figure 2B). We note that assessing competitive index by microscopy gives a quick and reproducible metric of one strain’s ability to outcompete another, but does not capture the full dynamic range of competition because many enumerated target cells have severe morphological abnormalities and are likely not viable at the 24 hr time point. However, microscopy allows competitive indices to be determined for these otherwise WT strains, which are not easily amenable to enumeration as colony forming units (CFU) because they form extraordinarily cohesive biofilms in isogenic co-cultures (dependent on type IV pili). These results show that SitA lipoproteins provide a competitive advantage, conferring the ability to kill and/or inhibit the growth of competitors in a TraA-dependent manner.

To examine SitA localization, we generated an mCherry reporter that carries the N-terminal lipobox from SitA1. Cells expressing this fusion have membrane-localized fluorescence as expected for a lipoprotein (Figure 2C). The TraA-dependent function of SitA shown in Figure 2B suggests that the protein is delivered by OME. Therefore, we tested whether the reporter fusion is transferred between cells. We co-cultured the reporter strain with a target strain expressing cytoplasmic GFP (which is not exchanged [Wei et al., 2011]) and microscopically assayed for transfer of the reporter. At 6 hr of co-culture, we observed the mCherry signal present in the cell envelope of the GFP target strain, indicating cell-to-cell transfer of the SitA1-mCherry reporter (Figure 2C, upper panel). Deletion of traA in the GFP target strain prevented the acquisition of mCherry signal (Figure 2C, lower panel), recapitulating our earlier findings that OM-localized reporters are exchanged between cells in a TraA/B-dependent manner (Nudleman et al., 2005; Wei et al., 2011; Pathak et al., 2012). Because inner membrane lipoproteins are not transferred during OME (Wei et al., 2011), these data also suggest that the lipobox directs SitA1 to the OM. We confirmed that full-length SitA1 localizes to the cell envelope using immunofluorescence microscopy to detect FLAG epitope-tagged SitA1 in formaldehyde-fixed cells (Figure 2D). Taken together, these results demonstrate that SitA1 resides in the OM and is transferred cell-to-cell by OME.

The fact that sitI1 expression protects WT DK1622 cells from swarm inhibition suggests that this gene encodes an immunity protein that neutralizes SitA1 toxicity. To determine whether SitI proteins block SitA-mediated growth inhibition, we expressed each sitI allele individually in DK1622 ∆sitBAI3 cells and co-cultured the resulting strains with strains that express each of the three sitBAI cassettes. For each competition, only strains that express the cognate sitI were protected from growth inhibition (Figure 2E), consistent with immunity function.

Immunity proteins typically interact with the C-terminal domain of polymorphic toxins (Zhang et al., 2012; Poole et al., 2011). To test whether this was true for SitA, we expressed the predicted C-terminal toxin domains (CTD) of each SitA toxin in E. coli MG1655 under the inducible P_BAD promoter. Expression of SitA1-CTD and SitA3-CTD blocked growth, but co-expression of cognate sitI from a second plasmid restored cell growth (Figure 2F). These results confirm that SitI proteins specifically neutralize cognate SitA toxins. In addition, because the SitA-CTD constructs lack secretion signal sequences, these data show that the domains exert their toxic effects in the cytoplasm. We also tested SitA2-CTD expression constructs, but none inhibited E. coli MG1655 growth. Because SitA2-mediated inhibition is obvious in M. xanthus competition co-cultures (Figure 2B and E), we tested the SitA-CTD expression constructs in M. xanthus and found that each inhibited cell growth (Figure 2—figure supplement 2A). Therefore, SitA2-CTD is indeed toxic when expressed in the cytoplasm of M. xanthus.

Given that sitBAI2 expression confers a significant advantage in competition co-culture (see Figure 2B), it is unclear why swarm inhibition was not observed with the sitA1^– nonmotile ancestral strain (see Figure 1B, row 5), considering these cells should still deploy SitA2 and that DK1622 lacks the SitI2 immunity protein. Therefore, we tested whether nonmotile ∆Mx-alpha cells that ectopically express sitBAI2 are able to inhibit DK1622 swarming. In agreement with the prior result, we found that DK1622 motility was only partially inhibited by the sitBAI2-expressing strain (Figure 1C, row 1). This result confirms that SitA2 contributes to the swarm inhibition phenotype, but is not sufficiently potent by itself to block outward swarming of DK1622. Together, these results indicate that SitA1 is the major swarm-inhibition toxin. We note that SitA3 does not contribute to the originally observed swarm inhibition phenotype because both ancestor and DK1622 strains contain the sitBAI3 operon (see Figure 2A). However, we found that a nonmotile strain expressing sitBAI3 fully inhibits the motility of a DK1622 ∆sitBAI3 strain that lacks the sitI3 immunity gene (Figure 1C, row 2).

Homologous CTDs are often associated with different toxin delivery systems from phylogenetically distant bacteria (Zhang et al., 2012). SitA3-CTD is homologous to a previously characterized tRNase domain found at the C-terminus of CdiA from Burkholderia pseudomallei 1026b (Morse et al., 2012; Nikolakakis et al., 2012) and an orphan CdiA-CTD encoded by Yersinia pseudotuberculosis YPIII (Figure 3—figure supplement 1). To determine whether SitA3-CTD also has tRNase activity, we expressed the toxin in E. coli and compared its activity to the CDI toxins. Induction of SitA3-CTD expression inhibited cell growth in the same manner as the CdiA-CTD toxins (Figure 3A). Examination of tRNA from SitA3-CTD intoxicated cells revealed cleavage of tRNA_UGC^Ala, similar to the specific tRNase activity of the B. pseudomallei toxin (Figure 3A).

Figure 3 with 1 supplement see all

Download asset Open asset

Next, we investigated the toxic activities of SitA1-CTD and SitA2-CTD. We had previously observed that SitA1 induces cell filamentation and loss of DAPI staining in M. xanthus target cells, which is consistent with DNase activity mediated by the predicted Colicin-DNase domain (Pfam 12639, E = 6.7 e-21) containing the WHH motif (Zhang et al., 2011). To test this, we induced SitA1-CTD expression in E. coli and found that cells became filamentous and had reduced DAPI stain signal (Figure 3B). By contrast, E. coli cells that were intoxicated by SitA3-CTD retained DAPI staining, though their nucleoids became more compact (Figure 3B). Together, these results suggest that SitA1-CTD has DNase activity. HMM-HMM comparison (HHpred [Söding et al., 2005]) of C-terminal residues 699–783 of SitA2 revealed distant homology to another CdiA-CTD from Y. pseudotuberculosis YPIII (locus tag, Ga0077885_11586), which was previously characterized as a DNase (Morse et al., 2015). To examine this possibility, we expressed each SitA-CTD in M. xanthus under the control of an IPTG-inducible promoter. Expression of SitA2-CTD in M. xanthus resulted in cell filamentation and reduced DAPI staining (Figure 2—figure supplement 2B), suggesting that SitA2-CTD degrades DNA. Expression of SitA1-CTD and SitA3-CTD in the cytoplasm of M. xanthus cells yielded similar results as when they were expressed in E. coli (Figure 2—figure supplement 2B), although the DAPI signal from SitA3-CTD expressing cells was brighter than the control and many cells contained two distinct nucleoids (Figure 2—figure supplement 2B), suggesting a block in cell division.

To determine the phylogenetic distribution of SitA toxins, we conducted a BLAST search using the N-terminal domains of SitA1/2 (which are homologous) and SitA3 as query sequences. This search recovered >100 sitA orthologs that are common in the Myxococcales (Supplementary file 1). More sensitive search algorithms such as HMMER (Finn et al., 2011) failed to identify significant homologs outside of the Myxococcales. Consistent with the finding that SitA is delivered through OME, all orthologs contain lipoboxes and are only found in species that contain traAB. Moreover, SitA C-terminal domains are variable and typically show homology to nuclease domains when subjected to HMM-HMM comparison using HHpred (Supplementary file 3). Interestingly, many sitA genes are not linked to upstream sitB orthologs, particularly in species that are distantly related to M. xanthus (Supplementary file 1). This suggests that SitB may not be required for SitA function, or perhaps that SitB proteins, encoded at unlinked loci, function promiscuously between multiple SitA proteins. Notably, some sitA loci are found outside of Mx-alpha-like elements, however, these genes are typically adjacent to other mobile genetic elements. To test cross-genotype compatibility of SitA orthologs, we cloned two sitA gene cassettes from Myxococcus fulvus HW1 for heterologous expression in M. xanthus. One of these operons (sitAI3^Mf1) does not contain a sitB gene. As predicted, M. xanthus cells that express heterologous sitAI3^Mf1 or sitBAI1^Mf1 outcompeted the parental strain in a traA-dependent manner (Figure 2—figure supplement 3). These results indicate that SitA toxin delivery is not limited by its specific species/strain of origin, and that the systems are functional after horizontal gene transfer (HGT) of a minimal set of components (sitAI).

Our results show that sitBAI-expressing cells inhibit OME-compatible strains that lack a cognate SitI immunity protein. We hypothesized that this antagonism should be sufficient to mediate territorial exclusion. Territorial exclusion promotes physical segregation of nonself organisms (Gibbs and Greenberg, 2011), which in turn drives both diversification and maintenance of cooperation (Papke and Ward, 2004; Velicer and Vos, 2009). Territorial exclusion between wild M. xanthus genotypes, including those that are closely related and that are isolated in close proximity to one another has been well studied, but the specific determinants underlying this behavior are unknown (Vos and Velicer, 2009; Rendueles et al., 2015; Vos and Velicer, 2006; Wielgoss et al., 2016). We tested whether SitA is sufficient for territorial exclusion by conducting colony merger assays with every combination of DK1622 strains expressing one of the five described sitAI alleles. In these assays, two liquid cultures are spotted next to one another on agar, and the colonies are allowed to swarm toward each other. If the converging swarms merge, then the strains are considered compatible. For each combination, the expression of different sitAI cassettes resulted in dramatic lines of demarcation between the two strains (Figure 4A), and the formation of these demarcation zones was traA-dependent (Figure 4B). These results show that otherwise isogenic strains are rendered socially incompatible and geographically isolated by the acquisition of a single sitAI cassette.

Figure 4

Download asset Open asset

We previously found that swarm inhibition occurs efficiently even when motile cells outnumber the antagonistic nonmotile strain 40 to 1 (Dey and Wall, 2014). This observation suggests that an individual SitA-expressing cell can inhibit many targets. To further explore this phenomenon, we quantified viable target cells in a series of competition co-cultures in which the ratio of SitA producing cells to targets was progressively decreased by factors of 10. To facilitate CFU enumeration in these experiments, we used ∆pilA cells, which are unable form type IV pili-dependent biofilms. We found that target cell CFU were reduced approximately 10⁶-fold when the strains were mixed at a 1:1 ratio (Figure 5A). The higher degree of killing reported here (compared to competitive index in Figure 2B) provides a clearer understanding of the killing efficiency because the CFU assay measures a broader dynamic range of viable cell number. Remarkably, the SitA-producing strain still reduced target cell viability >10⁴-fold in co-cultures seeded at a 1:1000 ratio of inhibitors to target cells (Figure 5A). These observations imply that each inhibitor cell intoxicates several thousand target cells during co-culture. In one explanation we hypothesized that OME delivery allows a series of SitA transfer events from one target cell to other cells. We consider this serial transfer mechanism plausible if translocation of all toxin molecules from the target cell OM to the cytoplasm is not completed before subsequent OME events occur. This model predicts that SitA toxins could spread through the population like an infectious agent, intoxicating target cells that never made direct contact with the producer. To test this hypothesis, we conducted three-strain competitions with (1) a sitBAI1 inhibitor strain that contains M. fulvus traAB (traAB^Mf) as its only traAB alleles, (2) a susceptible target strain that contains M. xanthus traAB^DK1622 alleles and thus is incompatible for OME with the inhibitor, and (3) a susceptible intermediary strain that carries both traAB^Mf and traAB^DK1622 (Figure 5B). If serial toxin transfer occurs, the traAB merodiploid strain should act as an intermediary carrier/conduit to deliver toxin to traAB^DK1622 targets (Figure 5B). As a control, we first showed that traAB^Mf inhibitors do not inhibit traAB^DK1622 targets (Figure 5C), consistent with the incompatibility of their traAB alleles. Importantly, inclusion of intermediary cells, which are inhibited (Figure 5D), also resulted in the inhibition of traAB^DK1622 target cells (Figure 5E). As expected, the intermediary strain (which lacks sitBAI1) did not inhibit traAB^DK1622 targets in co-cultures containing only those two strains (Figure 5F). To exclude a SitA-independent mechanism of target cell inhibition, we conducted the same three-strain competition, but provided traAB^DK1622 targets with the sitI1 immunity gene. In this latter co-culture, the intermediary strain was inhibited, but traAB^DK1622 targets were not (Figure 5G). Finally, we tested an intermediary strain that lacks traA and found that neither intermediary nor target cells were inhibited during co-culture (Figure 5H).

Figure 5 with 1 supplement see all

Download asset Open asset

To directly visualize serial transfer, we co-cultured traAB^Mf cells that express the SitA1-mCherry fusion (described in Figure 2) with traAB merodiploid intermediary cells that express cytoplasmic tdTomato, and GFP-labeled, traAB^DK1622 target cells. Microscopic examination of target cells at 6 hr revealed that all had acquired mCherry fluorescence, with the signal localized to the cell envelope (Figure 5I, left panels). By contrast, no SitA1-mCherry transfer was detected when the intermediary was absent from the co-culture (Figure 5I, right panels). Although these results are consistent with serial transfer of SitA, it is also possible that target cells acquire TraA^Mf and/or inhibitor cells acquire TraA^DK1622 by OME-dependent exchange of TraA with the intermediary strain, which would then allow direct transfer of SitA between inhibitor and target cells. To test if TraA is transferred during OME, we used a TraA-mCherry fusion to monitor transfer of TraA. The TraA fusion promoted efficient transfer of an OM sfGFP reporter (Figure 5—figure supplement 1A), demonstrating that it is functional to catalyze OME. However, TraA-mCherry itself was not transferred (Figure 5—figure supplement 1B). One explanation for why TraA does not transfer is that it interacts with TraB, which contains an OmpA domain known to bind the cell wall. For this reason, we suspect TraA is anchored to the cell envelope and is unable to transfer (Cao and Wall, 2017). Taken together, these results indicate that SitA1 can be transferred from the initial target to secondary recipients, supporting a model in which SitA acts like an infectious agent that disseminates through a population by OME.

In M. xanthus and its close relatives, sitA is typically accompanied by an overlapping sitB cistron. In the case of DK101 the genes overlap by 11 bp in all three sitBAI cassettes. SitB shows no significant homology to other proteins or domains using HMMER or HHpred, though it does contain a type I signal sequence (SignalP 4.1 [Petersen et al., 2011]). I-TASSER (Zhang, 2008) predicts that SitB adopts a transmembrane β-barrel structure characteristic of OM proteins. To examine the role of SitB1, we tested the activity of inhibitor cells that express either sitBAI1 or sitAI1 (cells lack sitB1) against a susceptible target strain. The inhibitors in these experiments also carried a ∆sitB3 mutation to eliminate the possibility of promiscuous interactions between SitB3 and SitA1. At a 1:1 (inhibitor to target) ratio, sitAI inhibitors had less of an advantage against targets than sitBAI inhibitors, but still retained activity compared to the mock inhibitor control (Figure 6A, left). The sitAI inhibitors were less effective at a 1:10 ratio, and at 1:100 were indistinguishable from mock inhibitors (Figure 6A). In contrast, sitBAI1 inhibitors were equally effective at outcompeting the target strain at each of the three ratios (Figure 6A). Thus, SitB1 contributes significantly to SitA1-mediated inhibition.

Figure 6

Download asset Open asset

Progressive loss of function at increasing target to inhibitor ratios with the sitAI inhibitors could indicate defects in serial toxin transfer compared to sitBAI inhibitors. Therefore, we tested serial transfer using three-strain co-cultures as described in Figure 5B. To improve the sensitivity of this assay, we increased the inhibitor to intermediary to target strain ratio to 10:1:1 to compensate for the inhibition defect of sitAI1 cells. As observed in Figure 5, sitBAI1 cells outcompeted both the intermediary and target strains (Figure 6B). Because of the high inhibitor cell ratio, the sitAI1 cells outcompeted the intermediary strain to a similar extent as sitBAI cells, but importantly, sitAI cells had little to no effect on target cells (Figure 6B). We extended the experiment to 48 hr and also performed experiments with SitA-resistant intermediary cells, to increase the number of conduit cells, but again we did not observe SitAI-mediated antagonism of target cells. These results support the hypothesis that SitB promotes SitA transfer, including the serial transfer from primary to secondary target cells. Although SitB1 clearly contributes to SitA1-mediated inhibition, it is not strictly required, which may explain why many myxobacterial sitA genes are not linked to sitB. Finally, these results are congruent with our above conclusion that TraA is not transferred, because if it was, then direct transfer of SitA1 would occur between sitAI inhibitors and the target strain. However, this did not occur because the target strain was not inhibited.

M. xanthus uses multiple inhibitory mechanisms to antagonize non-kin. However, because SitA toxins are serially transferred between cells, we hypothesized that they should be powerful determinants of competitive outcomes during inter-genotype conflict. To test this hypothesis, we quantified the contribution of TraA and SitA to competitive outcomes in co-cultures of DK1622 with wild M. xanthus soil isolates. Isolates A66 and A88 (from Tübingen, Germany) (Vos and Velicer, 2006) and DK801 (from California, USA) (Martin et al., 1978) contain traA alleles in the same recognition group as DK1622 (originally isolated from Iowa, USA) (Pathak et al., 2013; Dey et al., 2016). We compared the fitness outcomes of WT, ∆traA, and ∆sitBAI3 genotypes when co-cultured with these environmental isolates by monitoring the ratio of fluorescently labeled DK1622-derived cells to isolate cells at 4, 8 and 24 hr. As a second metric, we enumerated CFU of the DK1622-derived strains at the 24 hr time point. In every case, mutant strains that cannot deploy SitA3 had dramatically decreased competitive fitness outcomes and viability compared to WT (Figure 7A). Remarkably, against all three isolates, the presence of traA was the determining factor in which strain prevailed, demonstrating up to a 10⁶-fold swing in strain ratio (vs. A66) and a near 10⁷-fold difference in CFU (vs. DK801) between WT and ∆traA strains (Figure 7A). These results indicate that the ability to deliver SitA3 is a dominant determinant of competitive fitness under these conditions. Not surprisingly, the finding that environmental isolates outcompeted and killed ∆traA strains (Figure 7A), confirms the existence of OME-independent killing mechanisms at play. Interestingly, these experiments revealed that ∆traA cells had less competitive fitness than ∆sitBAI3, indicating a competitive fitness defect in ∆traA cells beyond just the inability to deploy SitA3. As a control, we competed the DK1622 genotypes against isolates A23 and A47 (from Tübingen [Vos and Velicer, 2006]), which are outside the DK1622 TraA recognition group (Pathak et al., 2013). ∆traA and ∆sitBAI3 genotypes resulted in similar competitive outcomes as WT as would be expected when inter-strain OME does not occur and SitA3 cannot be deployed (Figure 7B). These results show that SitA contributes significantly to fitness during competition with TraA-compatible non-kin genotypes, and likely plays a key role in competition and survival in nature.

Figure 7

Download asset Open asset

Here, we describe a novel family of proteins that carry polymorphic toxin domains and are delivered between myxobacteria by OME. SitA proteins are similar to other polymorphic toxins in that they carry diverse C-terminal domains, are neutralized by cognate immunity proteins, and are delivered in a cell contact-dependent manner. However, they are unique in their N-terminal domains, and in that they are lipoproteins transferred with other cargo during OME. To our knowledge, this is the first example of a polymorphic toxin system in which the toxin itself is a lipoprotein. Unlike CDI, T4SS, T6SS and T7SS toxins, there appears to be no requirement for a specialized apparatus to export the toxins. Instead, toxins are exchanged bi-directionally and simultaneously during OME. Therefore, SitA delivery likely only requires OM localization of the toxin and compatible TraA receptors. This discovery highlights the diversity of mechanisms used by bacteria to deliver polymorphic toxins.

Importantly, the SitA toxin family constitutes a second identity constraint upon OME with partner cells. For two cells to engage in OME, not only must they present compatible TraA receptors, but they must also contain immunity proteins to each other’s toxins. In this ‘recognize and verify’ system (Wall, 2016), if the latter constraint is not met, then the recipient of the toxin is poisoned. TraA homotypic interactions alone are considered kind or greenbeard recognition, in which social interactions are based on a single gene locus. The finding that myxobacteria verify relatedness with sitAI confirms the notion that myxobacteria apply a bona fide kin discrimination mechanism during OME by requiring identity verification at multiple polymorphic loci. Interestingly, this system allows TraA interactions to promote contrasting behaviors – cooperative or antagonistic – depending on relatedness. Either outcome makes OME potentially beneficial regardless of the partner by conferring the ability to both share goods with clonemates and poison non-kin.

SitA delivery range is restricted to within a single TraA recognition group. Considering traA allele diversity is high (Pathak et al., 2013), this significantly limits the use of SitA to related but nonself individuals. This suggests that one of the primary functions of SitA is the discrimination of exchange partners, consistent with the notion that sharing large amounts of goods with non-kin is costly. Within TraA recognition groups, OM material is a shared good; a resource to be guarded from exploitation by OME compatible, yet nonself populations. Myxobacteria achieve this safeguard by inextricably linking the delivery of these goods with the delivery of SitA toxins. Another example is the Burkholderia thailandensis CDI system, which couples a communication signal with polymorphic toxin delivery during biofilm formation (Anderson et al., 2014). Similarly, the CDI system of E. coli mediates both antagonism and cooperative intercellular adhesion to related cells (Ruhe et al., 2013, 2015). In Proteus mirabilis, IdsD/IdsE interactions communicate identity and may promote cooperative behaviors (Cardarelli et al., 2015). This intercellular communication is also coupled to toxic T6SS effector delivery (Wenren et al., 2013). In these examples, organisms link goods, signals, and/or cooperative behaviors to polymorphic toxin delivery, which ensures that potential cooperators are related.

Differential acquisition of antagonistic systems can affect cooperation compatibility between originally identical genotypes. We demonstrated that when two otherwise isogenic colonies express different SitA toxins, they are no longer able to merge swarms. The acquisition of a single sitAI operon would therefore alter strain identity and population structure between previously clonal cells. However, the impact of SitA on strain identity and population structure in natural soil habitats depends on two factors: (1) strains that belong to the same traA recognition group must exist in proximity, and (2) sitA loci must be sufficiently diverse to ensure different toxin/immunity types are represented at fine geographic scales. Velicer and colleagues have examined the compatibility of natural M. xanthus isolates obtained from a centimeter-scale plot of soil (Vos and Velicer, 2009; Wielgoss et al., 2016). Colony merger assays between these geographically proximal strains reveal compatibilities among the most closely related isolates, but also strong incompatibilities between strains that differ by only several dozens of mutations outside of the Mx-alpha region (Wielgoss et al., 2016). Indeed, these incompatibilities correlate with gene variation at hyper-variable Mx-alpha loci, where sitBAI genes commonly reside. Our analysis of their published sequences reveals 69 total and 15 unique sitA alleles distributed over 22 isolates (see Materials and methods for search criteria). Between these strains there are two traA alleles that are known to be incompatible (Pathak et al., 2013). Within each TraA recognition group we have observed an apparently high degree of correlation between the published colony merger compatibility of the strains (Wielgoss et al., 2016) and the sitA genes they possess. This suggests that sitBAI polymorphisms contribute to swarm incompatibility and inter-strain competition among natural soil isolates. We are currently investigating this possibility. Moreover, because Mx-alpha produces defective phage particles that promote specialized transduction (Starich and Zissler, 1989; Starich et al., 1985), these elements are apparent hotspots for HGT, perhaps explaining the high degree of Mx-alpha variation discovered between otherwise related isolates (Wielgoss et al., 2016). These variations likely contribute to the emergence of new compatibility types which underlie complex population structures, and explain the observation of rapidly evolving social antagonisms in M. xanthus (Velicer and Vos, 2009). Thus, TraA and SitA may act as powerful evolutionary drivers of myxobacteria diversification.

sitAI genes are often associated with prophage-like elements or other mobile elements and thus are likely acquired by HGT. By conferring a fitness advantage to their host, they may play an important role in transmission and retention of mobile DNA. For example, a HGT event into one cell in a population endowing it with a novel sitBAI operon allows that cell to infect its clonemates with toxins, thus ensuring the propagation of that element within the population. Similarly, the loss of the element would be lethal because the surrounding cells harbor this toxin-immunity pair would kill susceptible cells. Importantly, this model explains why many lab strains have stably maintained three large tandem repeats of Mx-alpha, which is expected to be genetically unstable (Starich and Zissler, 1989; Roth et al., 1996). In cases in which strains have spontaneously lost Mx-alpha units (Dey et al., 2016), those events likely occurred during propagation in liquid media, where OME cannot occur. Mx-alpha has the attributes of a selfish or addictive element that exploits the social nature of myxobacteria and OME. However, the origin of sitBAI loci is unclear because these genes also reside outside of selfish elements in myxobacterial genomes.

SitA toxins are uniquely powerful determinants of identity, likely because they are transmitted as infectious agents between recipient cells. The infectious model is consistent with the observations that individual M. xanthus cells typically make contact with multiple cells simultaneously within a swarm, OME is constitutively active, and that prodigious amounts of material are transferred during OME (Nudleman et al., 2005). Furthermore, SitA entry into the cytoplasm occurs by a secondary and uncoupled pathway to OME (Dey et al., 2016; Dey and Wall, 2014). Thus, it is possible that SitA lingers in the OM of the primary target long enough to allow transfer to secondary target cells through subsequent OME events (Figure 8A). Perhaps a cellular protein is required for SitA cell entry (Dey et al., 2016; Dey and Wall, 2014), but this protein is outnumbered by SitA molecules in the OM, making cytoplasmic entry a rate-limiting step. Based on these inferences, we propose two non-exclusive models that explain serial transfer: (1) Following transfer of SitA to a primary infected cell, OME with a secondary cell occurs before the full complement of SitA enters the cytoplasm of the initial recipient (Figure 8B); or (2) that three or more cells are engaged in OME simultaneously (Figure 8C). Our results suggest that SitB functions to promote serial transfer. Given that ∆sitB inhibitors have defects in direct transfer, serial transfer could be blocked simply by a decrease in number of SitA molecules delivered or a decrease in rate of delivery. Alternatively, SitB could stabilize SitA in the OM of the inhibitor and/or target cell, thereby increasing OM dwell times to promote serial transfer. The finding that SitB is an accessory protein is consistent with bioinformatics analysis in which many of the sitA genes that reside outside of the M. xanthus species are not linked with a sitB gene. In the case of M. xanthus isolates (Wielgoss et al., 2016), we found only 3 of 69 sitA genes were not associated with sitB. The mechanism of serial transfer, the function of SitB, and the ability of SitA to traverse the cell envelope are topics for future study. With respect to the latter point, in prior work we found that a mutation that disrupts the inner membrane protein, OmrA, renders target cells resistant to SitA1 (Dey et al., 2016; Dey and Wall, 2014). Thus, as proposed for CDI toxins, one possibility is that SitA toxins exploit inner-membrane proteins to gain access to the cytoplasm (Willett et al., 2015).

Figure 8

Download asset Open asset

Remarkably, SitA appears to determine the competitive outcome of two-strain co-cultures between DK1622 and wild-isolates, despite the T6SS and a host of other antagonistic machinery at play (Konovalova et al., 2010; Smith and Dworkin, 1994). We hypothesize that the rapid spread of SitA by serial transfer may disrupt the antagonistic capabilities of competitors. When deployed between converging swarms, serial transfer provides a mechanism to inhibit cells behind the front-lines. Thus, SitA mediated antagonism results in the formation of distinct territorial boundaries (see Figure 4), which in turn minimize OME, resource sharing, and social interactions between non-clonal swarms.

The selective pressure from SitA antagonism within a TraA recognition group may help drive the generation and fixation of traA polymorphisms that determine recognition specificity. For instance, we recently showed that simply substituting a single amino acid residue in TraA can alter recognition specificity while retaining OME function (Cao and Wall, 2017). More broadly, it is a puzzle how organisms select and maintain genetic variation in social genes such as traA. Diversification of beneficial greenbeard genes is theorized to be selected against because social groups with more common alleles receive benefit more often than those with less common alleles. This pressure to possess the allele that gains the most benefit is thought to ultimately erode allele diversity that originally allowed discrimination. This problem is known as Crozier’s paradox (Strassmann et al., 2011; Crozier, 1986). Our results provide one solution to this paradox in that variation at a second locus (sitA, a ‘harming greenbeard' [Gardner and West, 2010]) exerts selective pressure to diversify TraA, a helping greenbeard (Wall, 2016). For example, if one strain is killed by another via SitA, a TraA mutation that alters specificity within the losing population would be immune and retain OME function, and would thus be selected for. We suggest that other greenbeard systems could involve a similar balance between antagonism and cooperation that promotes maintenance of diversity for beneficial greenbeard genes.

Our results provide a description of a novel polymorphic toxin system that helps direct cellular goods to clonemates to promote multicellular cooperation. Polymorphic toxin systems are widely prevalent in bacteria and their role in population structure, ecology and evolution of microbes is only beginning to be understood. Importantly, this study highlights the diversity of delivery mechanisms for these domains and how they have adapted to the lifestyle of their host genomes.

All strains are listed in Supplementary file 2A. M. xanthus was routinely grown in CTT medium [1% casitone; 10 mM Tris⋅HCl (pH 7.6); 8 mM MgSO₄; 1 mM KH₂PO₄] in the dark at 33°C. E. coli TOP10 and MG1655 were grown in LB media at 37°C. As needed for selection or induction, media were supplemented with kanamycin (50 µg/mL), oxytetracycline (10 µg/mL), ampicillin (100 µg/mL), streptomycin (50 µg/mL), arabinose (0.2%), or IPTG (1 mM). TPM buffer (CTT without casitone) or PBS was used to wash cells. CTT or LB agar was used as a solid growth medium for routine strain maintenance. For all assays, strains were grown to logarithmic growth phase, washed, and re-suspended to the appropriate density.

All plasmids and primers are listed in Supplementary file 2B,C. Plasmids were constructed and maintained in E. coli and subsequently electroporated into M. xanthus. In the case of cloning IPTG-dependent sitA-CTDs, plasmids were maintained in XL1-Blue, which overexpresses LacI and reduces clone toxicity. Insertion mutations were created by amplifying an approximately 500 bp fragment of the gene of interest by PCR and cloning the fragment into the pCR-TOPO XL or pCR-TOPO 2.1 vectors (Invitrogen, Carlsbad, CA). For gene expression in M. xanthus, we cloned the appropriate gene(s) into pMR3487 using XbaI and NdeI restriction sites and T4 DNA ligase. If the gene(s) of interest contained these restriction sites, we used Gibson Assembly (Gibson et al., 2009) (New England Biolabs, Ipswich, MA) for plasmid construction. pMR3487 recombines at a specific site in the M. xanthus chromosome and expression is induced with IPTG (Iniesta et al., 2012). Traditional restriction endonuclease cloning was used to create pBAD30, pCH450 and pKSAT derived plasmids in which the sitA-CTD fragments had an ATG start codon engineered into the insert. Expression of sitI genes in pKSAT is constitutive, driven by the Km^r promoter. pCV10 for GFP expression was created by ligating tandem rRNA promoters, used for expressing lacI from the pMR3487 plasmid, with EGFP and the pSWU19 plasmid backbone. This plasmid recombines into the M. xanthus chromosome at the Mx8 phage attachment site. The deletion of sitBAI3 was constructed by cloning in-frame regions flanking and partially overlapping the start and stop codons of sitB3 and sitI3, respectively, into pBJ114 using Gibson Assembly. After recombination in M. xanthus, mutants were grown in CTT for 24 hr and plated on CTT containing 2% galactose to select for spontaneous loss of the galK marker. Deletion mutants were distinguished from WT by PCR with primers that flanked the deletion site.

Competition experiments, unless otherwise noted, were done using 1:1 strain mixtures of 3 × 10⁸ cells per mL spotted (20 µL) on agar plates containing 0.5× CTT with 2 mM CaCl₂ and 1 mM IPTG (competition media). Culture spots were harvested at 24 hr and observed on glass slides by microscopy to quantify strain ratios based on fluorescent labels. Typically, between 200 and 800 cells were counted. Competitive index in all assays was quantified by calculating the change in ratio of target to toxin-producing inhibitor cells over 24 hr. For example, if the 0 hr ratio was 1 to 1 and the 24 hr ratio was 1 to 100, the competitive index was. 01, indicating that the target strain was outcompeted. Swarm inhibition experiments were done identically but cells were not collected and instead were imaged after 72 hr. To determine the potency of killing, competition assays were conducted where the target cell volume and density were held constant (50 μL, 3 × 10⁸ cells per mL) while the number of toxin producing cells were titrated 1 to 10 for each sample. Cells were harvested at 48 hr, serially diluted and plated on CTT containing Km to enumerate viable target cells.

For the SitA1 serial transfer assay, competition was done at 1:5 or 1:5:5 mixtures of inhibitors to target(s) using a culture density of 3 × 10⁹ cells per mL. For competition with environmental isolates, liquid cultures of the DK1622-derived strain and the environmental isolate were adjusted to 3 × 10⁹ cells per mL liquid culture, mixed (100 µL DK1622 to 50 µL isolate), and spread onto competition media to ensure there were enough non-lysed cells to count by fluorescent microscopy. At the indicated time points, 2 mL of TPM was added to the agar plate, agitated with a plate spreader, and collected by pipette. Cells were centrifuged at low speed to help prevent clumping and either cell ratio was quantified by microscopy (as above) or the cells were resuspended in 1 mL TPM for CFU determination. Clumping was not an issue due to a non-isogenic mix of strains and massive cell lysis interfering with cell-cell adhesion. CFU of DK1622-derived strains were enumerated by 1 to 10 serial dilution and plating on CTT with oxytetracycline. All figures that contain error bars indicate the experiments were done in triplicate on different days. All statistical tests comparing two results are unpaired, two-tailed t-tests.

To test for growth inhibition by toxin expression, overnight cultures of each strain grown in LB with appropriate antibiotics were adjusted to OD₆₀₀ = 1.0 and back diluted 1 to 10 into fresh media containing 0.2% arabinose (to induce expression from pBAD), ampicillin 100 μg/mL and streptomycin 50 μg/mL. Cultures were then incubated in a shaker at 37°C for 5 hr. 1 to 10 serial dilutions of each culture were plated on LB with 0.2% arabinose and 100 μg/mL ampicillin and imaged after overnight growth at 37°C.

For DAPI staining, cells were grown as described above. At 0 hr and 6 hr post arabinose induction cells were adjusted to OD₆₀₀ = 0.4 and 1 mL was collected by centrifugation. Cells were fixed in freshly prepared 4% (vol/vol) formaldehyde in PBS for 15 min, rotating at room temperature (RT). The reaction was quenched by addition of an equal volume of 250 mM glycine (pH 7.5) in PBS. Cells were collected and washed 3x and resuspended to 100 μL. Cells were then spotted on poly-L-lysine coated slides and incubated for 10 min. Excess liquid was removed and cells were rinsed with water. 5 μL of fluorogel-II with DAPI (Electron Microscopy Sciences, Hatfield, PA) was placed on cells, a coverslip was applied and cells were imaged at 100× magnification with a Nikon E800 microscope coupled to a digital imaging system (Wei et al., 2011).

To assess tRNA processing, total RNA was isolated from E. coli clones using guanidine isothiocyanate-phenol extraction as described (Garza-Sánchez et al., 2006). Here, cultures were grown to mid-log phase and then diluted to OD₆₀₀ = 0.05 and grown for 30 min before 0.4% arabinose was added to induce toxin expression. tRNAs were analyzed by Northern blot hybridization using the following 5´-radiolabeled oligonucleotide probes: tRNA_UGC^Ala (5´ - TCC TGC GTG CAA AGC AG), tRNA_ICG^Arg (5´ - CCT CCG ACC GCT CGG TTC G), tRNA_CGA^Ser (5´ - GTA GAG TTG CCC CTA CTC CGG), and tRNA_GCA^Cys (5´ - GGA CTA GAC GGA TTT GCA A).

These assays used a modified competition media (replacing 1.5% agar with 1.0% agarose for imaging and addition of 5 μg/mL oxytetracycline). A multichannel pipette was used to simultaneously pipette competing strain suspensions (1.5 × 10⁹ cells per mL) until culture spots were nearly touching. Aliquots were air dried and plates were then incubated in a humid chamber at 33°C for 3 days. Spots were viewed on an Olympus SZX10 stereomicroscope coupled to a digital imaging system.

To discover SitA homologs, we performed BLAST (Altschul et al., 1990) analysis against the IMG database (Markowitz et al., 2012), using the conserved N-terminal sequence (first 507 amino acids of SitA1 and the first 441 amino acids of SitA3) as queries. HMMER (Finn et al., 2011), HHPred (Söding et al., 2005) and I-TASSER (Zhang, 2008) analysis were performed as described above using default parameters.

To detect any correlation between colony merger and sitA genes from the referenced study (Wielgoss et al., 2016), we first identified SitA homologs by performing local BLAST analysis of the sequences published at http://www.odose.nl/u/michiel/h/22-myxo-genomes-w-annotation (Wielgoss et al., 2016) using SitB1 as the query. Any gene located downstream of a SitB homolog that contained the described features of SitA was considered a SitA homolog. We also used SitA homolog sequences as queries to find any sitA genes without an accompanying sitB. Redundant sequences were removed using NCBI FASTA Tools Unique Sequences webpage. The sequences were then clustered according to >96% pairwise amino acid identity and each cluster was considered a unique polymorphic toxin group. The two closest toxin groups were 90% identical. This analysis was done blind with respect to the published colony merger compatibility types (Wielgoss et al., 2016).

Log-phase reporter and target cell liquid cultures were adjusted to 1.5^–3 × 10⁹ cells per mL, mixed at the indicated ratios, and plated on ½ CTT agar with 2 mM CaCl₂, +/- IPTG, as needed. Strains were incubated at 33°C for the indicated time-period, collected, and visualized by fluorescent microscopy as described (Wei et al., 2011).

Log phase cultures were resuspended to 6 × 10⁸ cells per mL in reduced osmolarity PBS solution (“mPBS” = 7.2 mM NaCl, 5.4 mM KCl, 10 mM Na₂HPO₄, 1.8 mM KH₂PO₄). 100 μL of 10% paraformaldehyde solution (in mPBS) and 1 μL of 5% glutaraldehyde solution (in H₂O) were added to 400 μL of cell suspension. Each mixture was spotted on a poly-lysine coated slide. Fixation proceeded for 30 min at RT. After rinsing with mPBS, cells were permeabilized with 0.025% Triton X-100 for 10 min and washed. Cells were blocked for 1 hr at RT with 4% BSA in mPBS, then probed with a 1:1000 final concentration of anti-FLAG antibody (in 4% BSA, Sigma, St. Louis, MO) for 1 hr at RT and subsequently washed 2× for 5 min and 1× for 10 min with mPBS. Cells were then probed with a 1:2000 final concentration of secondary antibody (in 4% BSA, Alexa Fluor 488-conjugated donkey anti-rabbit IgG; Jackson ImmunoResearch, Westgrove, PA) for 30 min at RT and subsequently washed as before. SlowFade Gold antifade reagent (Invitrogen) was added to the slide and the cells were visualized with a 100× objective lens. Western blot was performed according to standard protocols, using anti-FLAG antibody described above, and horseradish peroxidase-conjugated goat anti-rabbit secondary antibody (Thermo Scientific, Waltham, MA).

Each SitA-CTD was expressed from an IPTG-inducible promoter. For liquid growth-inhibition, cells were grown to log-phase and diluted to 5 × 10⁷ cells per mL in two separate flasks containing fresh media (CTT, 2.5 μg/mL oxytetracycline). To one of two flasks, IPTG was added to 1.0 mM. Cells were grown in the dark with shaking at room temperature. Culture growth was monitored at the indicated time points by measuring turbidity using a Klett meter. For DAPI staining, cells were grown as above for 30 hr, washed and resuspended in mPBS to a concentration of 1.5 × 10⁹ cells per mL. 1 μL of a 50 μg/mL DAPI solution (Life Technologies, Carlsbad, CA) was added for each mL of culture. Cells were incubated for 20 min in the dark with rotation, washed, concentrated, and visualized by fluorescent microscopy.

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

We would like to thank Dr. Greg Velicer and Dr. Sébastien Wielgoss for constructive comments.

© 2017, Vassallo et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.