The Cost of Virulence: Retarded Growth of Salmonella Typhimurium Cells Expressing Type III Secretion System 1

1The ability to infect a host and elicit disease is dictated by the virulence factors expressed by a given pathogen. This may include, but is not limited to, protective factors neutralizing antibacterial defenses, enzymes involved in nutrient acquisition within the host, regulators of virulence factor expression and toxins or secretion systems for subverting host cell signal transduction. The coordinated expression of such virulence factors enhances colonization, growth/survival within the host and transmission. However, most studies of virulence factor function and pathogen fitness have been performed in bulk assays, averaging the phenotypes over large pathogen populations of genetically identical cells. In contrast, little is known about the potential advantages, costs or burdens arising from virulence factor expression by an individual cell of the pathogen population. Therefore, single cell analyses might be of significant interest, in particular if virulence factors, which are expressed in a bistable fashion by some but not all members of a pathogen population, e.g. the ttss-1 system of S. Typhimurium,,,,, as described in this paper.

2Bistable gene expression is genetically encoded. In most cases, one particular genotype expresses one predictable phenotype in a given environment. However, in some cases, two different phenotypes are expressed by isogenic organisms living in the same environment. This is termed phenotypic variation, bimodal gene expression or bistability and represents a special case of gene expression. The importance of bistability for pathogenic bacterial fitness and evolution is just beginning to be understood.

3Like other cases of gene expression, bistability is generally observed in response to particular environmental cues. The response is driven by a dedicated (set of) regulator(s), which responds to environmental signals (operon model of Jacob ). This response is subject to stochastic fluctuations. In particular in the case of regulators expressed in a few copies per cell, this can significantly affect the active regulator concentration thus randomizing the corresponding phenotype in a population,. In combination with non-linear responses (e.g. regulator multimerization, feedback loops), this can lead to formation of phenotypically distinct and stable subpopulations of isogenic bacteria,,,,. In terms of evolution, two models may explain the advantage of bistability: i. in “bet hedging”, the optimally adapted phenotype will prevail and ensure the survival of the shared genotype in a changing environment. ii. in “division of labor”, both phenotypes cooperate to ensure survival of the shared genotype. In either way, the bistable expression of certain genes is thought to promote the survival of the genotype. However, it has remained poorly understood whether/how bistability may affect the lifestyle of pathogenic bacteria.

4Salmonella enterica subspecies 1 serovar Typhimurium (S. Tm) is a pathogenic Gram-negative bacterium causing numerous cases of diarrhea, worldwide. Its' type III secretion system 1 (TTSS-1) was recently identified as an example for bistable gene expression,,,. TTSS-1 is a well-known virulence determinant of S. Tm required for eliciting diarrheal disease,,. The needle like TTSS-1 apparatus injects effector proteins into host epithelial cells, thus triggering host cell invasion and pro-inflammatory responses,,. TTSS-1 is encoded on a genomic island (Salmonella pathogenicity island 1 (SPI-1)), which also harbors genes for effector proteins and for several regulators of ttss-1 expression, e.g. hilA, hilC and hilD,.

5The bistable ttss-1 expression is controlled by a complex regulatory network, which includes coupled positive feedback loops, controls the threshold for ttss-1 induction and amplifies ttss-1 expression,. Bistable ttss-1 expression is observed in “ttss-1 inducing” environments, i.e. the gut lumen of infected mice or in non-host environments, e.g. when S. Tm is grown to late logarithmic phase in LB,,,. This yields mixed populations of isogenic S. Tm cells that express ttss-1 (TTSS-1+ phenotype), or do not (TTSS-1− phenotype), in a bimodal fashion. In the mouse gut, only the TTSS-1+ cells can actively invade the mucosal tissue and efficiently trigger inflammation,. This inflammatory response may help to overcome the commensal microflora, thus enhancing Salmonella growth and transmission,,,,,,. Experimental data indicate that bistable ttss-1 expression might represent an example of “division of labor”, but further data is required to settle this point. At any rate, ttss-1 expression seems to be instrumental for eliciting diarrheal disease and enhancing pathogen transmission. But the functional properties of the TTSS-1+ phenotype are not well understood.

6The complex setting of the infected animal gut has hampered the analysis of the TTSS-1+ phenotype. In vitro experiments are essential for gaining detailed mechanistic insights. Here, we have analyzed the induction of ttss-1 expression and its effects on the growth rate of the TTSS-1+ phenotype by single cell reporter assays, competitive growth experiments and mathematical modeling. In such non-host environments, expression of the ttss-1 virulence system expression imposed a growth penalty on the TTSS-1+ cells. This may have important implications with respect to compensatory mechanisms during the infection of animal hosts.

7We started our analysis of the TTSS-1+ phenotype by probing ttss-1 expression at the single cell level. For this purpose, we chose the sicA promoter (PsicA), which controls expression of the chromosomal sicAsipBCDA operon (Fig. S1C). This operon encodes key parts of the TTSS-1 virulence system. On the one hand, we employed a transcriptional sipA-tsrvenus reporter gene cassette placing the reporter downstream of the sicAsipBCDA operon (Fig. S1;, ). Due to its localization at the bacterial poles, the tsrvenus reporter allows detecting <10 proteins per cell. Thus, sipA-tsrvenus provides a highly sensitive reporter for the TTSS-1+ phenotype.

8Next, we verified the performance of the sipA-tsrvenus reporter. sipA-tsrvenus expression was bistable and TTSS-1− and TTSS-1+ individuals were distinguishable by the presence/absence of Tsrvenus spots at the bacterial poles (; Fig. 1A; Fig. S1D). TTSS-1 expression and virulence were not compromised (Fig. 1B). The accurate response of sipA-tsrvenus to Salmonella signaling cascades was established by disturbing known elements of the TTSS-1 gene regulation network and FACS analysis of sipA-tsrvenus expression (Fig. 1C, D). In line with the published work on ttss-1 regulation (Fig. 1D): i. Over-expression of positive TTSS-1 regulators increased the abundance of tsrvenus-expressing individuals (Fig. 1C; Fig. S1D). In particular, hilA, hilC and hilD over-expression increased the fraction of sipA-tsrvenus expressing individuals from ∼20% to 80–100%. ii. The median signal intensity per sipA-tsrvenus expressing cell increased when positive regulators were over-expressed (philA: 3.8±0.3-fold; philC: 4.0±0.1-fold; philD: 4±0.1-fold; median ± s.d.). iii. Control experiments in a ΔhilA mutant verified that expression of the TTSS-1+ phenotype depended on the ttss-1 master-regulator, HilA (Fig. 1C; open bars) and iv. The average HilA protein levels of the analyzed strains correlated positively with the fraction of tsrvenus-expressing individuals (r2 = 0.78; quantitative Western blot; Fig. 1E). These data verified the accurate performance of the sipA-tsrvenus reporter and demonstrated that hilA-dependent regulation affects both, the fraction of TTSS-1+ individuals and the level of ttss-1 expression per cell.

9In addition, we employed psicA-gfp, a reporter plasmid expressing gfp under control of the sicA promoter. This construct yielded brighter fluorescence than the chromosomal sipA-tsrvenus and was better suited for FACS analysis. Again, this reporter yielded a bistable expression pattern (Fig. 1F). Using wt S. Tm psicA-gfp we separated TTSS-1+ and TTSS-1- subpopulations by FACS. Western blot analysis of the FACS-sorted subpopulations verified coincident expression of psicA-gfp and the TTSS-1 protein SipC (Fig. 1F, G). This indicated that our fluorescent reporter constructs are faithful reporters of the bistable expression of the TTSS-1+ phenotype.

10During our experiments, we observed that hilA, hilC and hilD over-expression led to reduced culture densities (e.g. OD600 for wt sipA-tsrvenus: 3.4±0.3 vs. wt sipA-tsrvenus philA: 2.0±0.3; mean ± s.d.). This was a first hint suggesting that retarded growth might be a general feature of the TTSS-1+ phenotype. However, it remained to be shown whether growth retardation occurs in wild type cells expressing normal levels of hilA, hilC and hilD.

11The growth rate of the TTSS-1+ individuals was analyzed by time-lapse microscopy. Wild type S. Tm harboring gfp- or tsrvenus-reporters for ttss-1 expression were placed on an agar pad (LB, 1.5% agarose), the TTSS-1+ individuals were identified by fluorescence microscopy and growth was analyzed by time-lapse phase contrast microscopy (1 frame/30 min; Fig. 2A). Imaging did not impose detectable photo damage to the bacteria, as indicated by the unaltered growth rate (Fig. S2). Strikingly, TTSS-1+ individuals grew slower than TTSS-1− individuals (wt S. Tm sipA-tsrvenus (M2001); µT1+ = 0.90 h−1 vs. µT1 − = 1.30 h−1; p = 0.027 for the factor ‘phenotype’ in a two-way ANOVA; Fig. 2B). The negative control strain ΔhilA sipA-tsrvenus yielded only TTSS-1− individuals, which grew at the “fast” rate (µT1 − = 1.16 h−1; Fig. 2B). Thus, TTSS-1+ individuals seemed to grow at a reduced rate.

12To exclude potential artifacts attributable to the sipA-tsrvenus reporter, we analyzed unmodified wild type S. Tm not harboring any reporter (Fig. 2C; Fig. S2). Using a maximum likelihood approach, we identified two populations with distinct growth rates (likelihood ratio test for two populations versus one population, p<0.001, µslow = 0.66 h−1 vs. µfast = 1.27 h−1; Fig. 2C), very similar to the ones described above (Fig. 2B). Furthermore, unmarked mutants lacking the entire SPI-1 region (Δspi-1) or the positive ttss-1 regulator hilD yielded exclusively fast growing cells, while deletion of the negative ttss-1 regulator hilE yielded only slow growing cells (Fig. 2C). Finally, wild type S. Tm harboring psicA-gfp or a chromosomal gfp-reporter for the TTSS-1 gene prgH yielded slow growing TTSS-1+ and fast growing TTSS-1− cells (µT1+ = 0.51 h−1 vs. µT1 − = 1.2 h−1; p = 0.006 for the factor ‘phenotype’ in a two-way ANOVA; Fig. 2D). Bacteria expressing the psicA-gfp or prgH-gfp reporters grew even slower than the TTSS-1− sipA-tsrvenus bacteria or the slow-growing wt S. Tm subpopulation (Fig. 2BC). Presumably, this was attributable to the additional “burden” conferred by the GFP expression, as described, before.

13Thus, the time-lapse microscopy experiments verified bistable ttss-1 expression and revealed that the TTSS-1− phenotype has a reduced growth rate, even at wild type HilA and TTSS-1 levels (µT1+ in the range of 0.7 h−1 vs. µT1 − in the range of 1.3 h−1). This was confirmed in a dye dilution assay (Fig. S3).

14Our data suggested that ttss-1 expression represents a “cost” to the bacterial cell. However the mechanism explaining this growth retardation had remained unclear. We speculated that expression of the TTS apparatus itself or the sheer load of the proteins transported by the TTSS-1 (effectors, translocon proteins) might play a role. To test these hypotheses, we analyzed two additional S. Tm mutants. In the first mutant, termed Δprg-orgΔinv-spa, we deleted most apparatus-encoding genes (Table S1). This mutant formed two populations with distinct growth rates (likelihood ratio test for two populations versus one population, p<0.001, µslow = 0.72 h−1 vs. µfast = 1.36 h−1; Fig. 2E), very similar to those described for wild type S. Tm (Fig. 2C). The second mutant, termed Δ8Δsip, was lacking the genes for most TTSS-1 effector proteins and the secreted translocon components including sipB, sipC, sipD, sipA, sptP, sopE, sopE2, sopB and sopA (Tab. S1). In contrast to wild type S. Tm, we could not distinguish two subpopulations in this mutant (likelihood ratio test for two populations versus one population, p = 0.73; Fig. 2E). Instead, this mutant displayed a median growth rate of µ = 1.10 h−1, similar to the fast growing subpopulation of S. Tm wt and the mutants Δspi-1 and hilD (Fig. 2C). This data suggests, that expression of the effector proteins and translocon components is “costly” and provides at least in part a mechanistic explanation for the growth retardation of wild type S. Tm cells of the TTSS-1+ phenotype.

15When monitoring growth and bistable ttss-1 expression in a wt S. Tm (psicA-gfp) culture, the fraction of TTSS-1+ individuals began to rise after 2.5 h as soon as the culture entered the late logarithmic phase, increased in a linear fashion, and reached approx. 60% after 7 h once the culture entered the stationary phase (Fig. 3A).

16Our results implied that two different parameters affect the fraction of TTSS-1+ individuals and the overall growth progression in the late logarithmic phase: i. Competitive growth. TTSS-1+ individuals are steadily outgrown by the fast-growing TTSS-1− individuals (µT1+<µT1 −; Fig. 2); this constantly reduces the size of the TTSS-1+ subpopulation. ii. ttss-1 induction. Presumably, initiation of ttss-1 expression in TTSS-1− individuals compensates the “TTSS-1+ losses” attributable to competitive growth and explains the increasing fractions of TTSS-1+ individuals during the late logarithmic phase.

17To infer the dynamic initiation rate ri of ttss-1 expression in the late logarithmic phase from our experimental data, we devised a mathematical model describing the growth of the TTSS-1+ (NT1+; growth rate µT1+) and the TTSS-1− population (NT1 −; growth rate µT1 −) as a function of time (t):

18It should be noted that the model does not include a term for “switching off” ttss-1 expression. This was justified by our failure to observe “off switching” in the experiments shown in Fig. 2 and further supported by other data (Fig. S2 and data shown below). During the late logarithmic phase, the relative abundance of the TTSS-1+ individuals increased, and the fraction α of TTSS-1− individuals (NT1-) decreased in a linear fashion (Fig. 3A):

19Equation (2) can be rearranged to calculate ri(t) (see Text S1 for details):

20With the data from Fig. 3A and by using equation (3) we could determine NT1 − (t) and, after fitting an empirical function to NT1 − (t), also dNT1 −/dt. Using equation (4), this allowed calculating ri(t) during the late logarithmic phase (see Text S1 for details). We found that the mean initiation rate (ri) of ttss-1 expression increased continuously during the late logarithmic phase, e.g. from 0.28 h−1 at 3.5 h to 0.54 h−1 at 5.5 h (SEM = 0.03 h−1; Fig. 3B).

21The initiation rate of ttss-1 expression seemed to increase upon entry into the late logarithmic growth phase (Fig. 3A). Therefore, it might be induced by growth-related environmental signals (e.g. oxygen depletion, quorum signals, nutrient depletion, metabolite accumulation). To address this, we analyzed the partial oxygen pressure (pO2) during growth. As expected, pO2 declined to <30% relative aeration during the first three hours (Fig. 3C). After approximately 3.5 h, we detected a transient rebound of the oxygen pressure followed by a steady decline to <3% relative aeration during the next hour. This undulation of oxygen pressure is indicative of a change in the growth physiology at 3.5 h and was in line with the reduced growth rate (Fig. 3A, shaded area).

22The data suggested that altered metabolism, nutrient availability, waste product accumulation, the reduced growth rate or the low oxygen pressure might represent cues inducing ttss-1 expression. As a first approach to test the role of pO2, we performed batch culture growth experiments in identical 250 ml culture flasks filled with the indicated volumes of media (wt S. Tm psicA gfp grown in 5, 10, 25, 50 or 100 ml LB; Fig. 3D). This setup allowed analyzing the effect of reduced pO2 (i.e. in larger, poorly aerated culture volumes) at equivalent growth rates. We observed that the fraction of ttss-1 expressing cells increased in larger culture volumes. Therefore, low oxygen tension might represent one environmental cue directly or indirectly inducing bistable ttss-1 expression. However, the evidence is merely circumstantial at this moment and other cues might well be involved. Identification of these cues will benefit from the strategies for determining ri as described above.

23In liquid culture, the initiation of ttss-1 expression occurred in the late logarithmic phase. However, our initial time lapse microscopy data for bacteria sampled from this growth phase did not show initiation of ttss-1 expression (Fig. 2). We reasoned that this might be attributable to the lack of inducing environmental signals, as these experiments had been performed on agar pads soaked with fresh LB medium. To test this hypothesis, we modified the time lapse microscopy experiment and imaged bacteria (S. Tm psicA-gfp) placed on agar pads soaked with filter-sterilized spent medium taken from a culture at the same growth phase (OD600 = 0.9, see Materials and Methods). We analyzed growth of 191 micro colonies. At the beginning, 135 did not express ttss-1. But remarkably, we observed 15 of 135 initially TTSS-1− micro colonies, in which individual bacteria induced ttss-1 expression during the course of our imaging experiment (e.g. Fig. 4A, Fig. S4; Video S1). After induction, the TTSS-1+ cells grew at a slower rate than their TTSS-1− siblings. In addition, we observed numerous TTSS-1+ bacteria (56 micro colonies) and TTSS-1− bacteria (120 micro colonies) which did not “switch” their ttss-1 expression status. In line with the results above, ttss-1 expression and the interval between two cell divisions was negatively correlated (Fig. 4A,B,C, Spearman's rho = −0.747, p<0.0001, N = 29).

24These experiments support the stochastic initiation of ttss-1 expression. But the initiation rate of ttss-1 expression (<0.04 h−1) was lower than that predicted from the batch culture experiment shown in Fig. 3 (ri = 0.18−0.45 h−1). This might be attributable to the lack of some environmental cue, e.g. low oxygen pressure, as time lapse microscopy was performed at ambient atmosphere. Only two micro colonies showed a decrease in fluorescence as expected for “off-switching”. Hence, the rate of off-switching is not substantial. This indicated that our mathematical model, which assumed that “switching off” the ttss-1 expression would be negligible, was justified (equation (1) did not include ri(t)NT1-(t)). These experiments verified that ttss-1 expression is initiated in a stochastic fashion under “inducing” environmental conditions and that the TTSS-1+ phenotype exhibits a growth defect.

25Finally, we confirmed the growth penalty attributable to ttss-1 expression in the late logarithmic phase in competition experiments. Wt S. Tm expresses ttss-1 in a bistable fashion and forms a significant fraction of slow-growing TTSS-1+ cells during the late logarithmic phase (Fig. 3). This slows down the apparent growth of the total wild type population (see above). In contrast, hilA or hilD mutants, which do not express ttss-1, yield a pure population of fast-growing TTSS-1− cells (Figs. 1 and 2). Thus, in a mixed culture, hilA or hilD mutants should outgrow wt S. Tm. Indeed, both mutants out-competed the wt strain during the late logarithmic phase of the mixed culture (ΔhilA, ΔhilD; Fig. 5A,B). In contrast, a hilE mutant, which forms a larger fraction of TTSS-1+ cells than wt S. Tm (Fig. 2), was outcompeted by wt S. Tm in this type of assay (ΔhilE, Fig. 5C). This verified the growth penalty of TTSS-1+ cells in LB batch cultures.

26The effect of virulence factor expression on the fitness of an individual pathogen cell has remained unclear. We have analyzed the fitness costs associated with the expression of ttss-1, which encodes a key virulence function of S. Tm. An in vitro system was chosen for a detailed analysis of the growth phenotype of TTSS-1+ cells. We found that these cells have a reduced growth rate. This established that ttss-1 expression represents a burden (and not an advantage) at the level of the individual cell, at least in the non-host environment of our assay system. The growth penalty affects the fraction of TTSS-1+ individuals and the overall growth progression in a S. Tm culture. Mathematical modeling and experimental data demonstrated that this growth penalty and an increasing initiation rate of ttss-1 expression during the late logarithmic growth phase were sufficient to explain the dynamic abundance of TTSS-1+ and TTSS-1− individuals in a clonal S. Tm batch culture.

27Evidence for bistability of ttss-1 expression has only recently been accumulated. Under inducing conditions, single cell reporters for expression of ttss-1 or effector proteins yielded cells in the “on” and cells in the “off” state,,,,. The regulatory network controlling ttss-1 expression includes at least three positive feedback loops and this architecture is thought to set the threshold for initiating ttss-1 expression and to amplify the level of expression,,. The TTSS-1+ phenotype can persist for several hours, even if the bacteria are shifted into environments normally not inducing ttss-1 expression (histeresis; shift to fresh LB, Fig. 2; Fig. S2). However, it should also be noted that it has not been possible to define unequivocally where stochasticity is introduced. In fact, stochastic initiation of ttss-1 expression might hinge on different regulators in different environments.

28TTSS-1+ cells have at least two important characteristics. First, they express the virulence factors enabling host manipulation and elicitation of disease,,. Second, as we have found here, they grow at a reduced rate. ttss-1 expression may represent a “burden” in itself. The mechanism explaining the growth defect of TTSS-1+ cells is of significant interest. A partial disruption of the proton gradient by “leaky” TTSS assembly-intermediates and/or the metabolic energy required for biosynthesis of the TTSS may offer plausible explanations. Typical TTSS-1+ cells are estimated to express 20–200 TTS apparatuses and approx. 3−10×104 effector proteins, amounting to a significant fraction of the total cellular protein,. Indeed, deleting the translocon and most effector proteins significantly increased the growth rate of the TTSS-1+ cells (Δ8Δsip; Fig. 2E), indicating that these proteins account at least in part for the cost of ttss-1 expression. However, the growth rate of Δ8Δsip (µ = 1.10 h−1) was still lower than that of the TTSS-1− subpopulation of wt S. Tm (µfast = 1.27 h−1), suggesting that other factors do also contribute to growth retardation.

29An alternative explanation for the reduced growth rate of TTSS-1+ cells might reside in coordinated expression of a complex regulon. This might be reminiscent of the prf virulence regulon of Listeria monocytogenes, which coordinates metabolism and virulence gene expression thus controlling environment-specific fitness phenotypes in vitro and in vivo. Several global regulators (e.g. crp, mlc, fur;,, ) and silencing proteins (hns, hha;, ) can control ttss-1 expression. Moreover, HilA may control multiple loci apart from ttss-1 (25). And we have observed co-expression of ttss-1 and of fliC, which encodes a key structural component of the flagella, in the late logarithmic phase (Fig. S5). Accordingly, ttss-1 expression might be one feature of a “differentiated” state which also includes adaptations reducing the growth rate. It is tempting to speculate that this state might be particularly adapted for mucosal tissue invasion. This would be an important topic for future research.

30Interestingly, similar phenomena have been observed in other ttss-expressing pathogens. In Pseudomonas aeruginosa, growth in suboptimal media was shown to result in bistable ttss expression. But it remained unclear whether growth might be affected. In contrast, the plasmid-encoded TTSS of Yersinia spp. is well known to cause growth retardation in response to host cell contact or low calcium environments,. However, in this case, ttss induction seems to be uniform even in suboptimal media. Thus, bistability and growth retardation do occur in other ttss expressing bacteria, but specific adaptations may exist for each pathogen.

31Which environmental cues induce ttss-1 expression in S. Tm? ttss-1 is expressed in the lumen of the host's intestine and in the late logarithmic phase in LB-batch culture. Low oxygen pressure is common to both environments and may represent an inducing signal (see Fig. 3C). In line with this hypothesis, Shigella flexneri, a closely related gut pathogen, can modulate the activity of its TTSS in response to low oxygen pressures typically observed at the gut wall. Similarly, HilA-mediated ttss-1 expression is known to respond to oxygen pressure,. In addition, numerous other internal and external cues are known to affect ttss-1 expression, including osmolarity, pH, growth rate, or the presence of short chain fatty acids like acetate,,,,,,. The sum of these environmental cues seems to determine the level of ttss-1 induction. This might explain our observation of a low, but detectable initiation rate of ttss-1 expression on agar pads soaked with spent medium (Fig. 4). This environment should harbor most cues present in the late log culture medium, but lacks low oxygen pressure, which could not be established in the real time microscopy setup.

32In summary, our findings indicate that the TTSS-1+ phenotype is more complex than previously anticipated. Currently, we can only speculate how this affects the real infection and transmission in vivo. Our results suggest that the TTSS-1+ subpopulation is constantly drained by the burdens inflicted by immune defenses within the infected gut mucosa and by the reduced growth rate (this work). The latter should represent a competitive disadvantage against all other bacteria (commensals and TTSS-1− S. Tm cells) present in the gut lumen. Moreover, this burden should materialize even before invading the gut tissue and may explain why ttss-1 defective mutants are sometimes (though rarely) found in infected animal flocks and isolated in one case of a human outbreak,. In order to explain the evolution and mainentance of bistable ttss-1 expression and the successful propagation of the ttss-1 genotype, one has to predict that the TTSS-1+ phenotype must confer some type of advantage. According to the “division of labor” model, the advantage might emanate from a “public good”, i.e. the TTSS-1 induced gut inflammation fostering Salmonella growth in the gut lumen and enhancing transmission. Alternatively, the TTSS-1+ phenotype might include (unidentified) features enhancing the survival and growth of the ttss-1 expressing bacteria themselves, e.g. in permissive niches of the host's intestine or by enhancing the chances of chronic infection and long-term shedding. Identifying these mechanisms will represent an important step for understanding the evolution of bistable ttss-1 expression.

33All strains were derivatives of Salmonella Typhimurium SL1344 or ATCC14028 (see Tab. S1 and Text S2 for references). All plasmids and primers are shown in Tab. S2 and S3. Bacteria were inoculated (1∶100 in LB) from 12 h overnight cultures (LB, supplemented with the appropriate antibiotics) and grown under mild aeration for 4 h at 37°C, if not stated otherwise. In Fig. 1C,E, the medium included 0.01% arabinose.

34The mutants were constructed using the lambda red recombination system. The chloramphenicol or kanamycin resistance cassette of pKD3 (cat) resp. pKD4 (aphT) were amplified by PCR using the primer pairs ÄhilA::kan-fw and ÄhilA::kan-rev, ÄhilD::kan-fw and ÄhilD::kan-rev, ÄhilE::cat-fw and ÄhilE::cat-rev and electroporated into SL1344 harboring pKD46 to generate the regulator mutants M2005 (ÄhilA::aphT), M2007 (ÄhilD::aphT) and M2008 (ÄhilE::cat). Mutants were selected by plating on LB-Agar (50 µg/ml kanamycin or 30 µg/ml chloramphenicol). M2072 (termed Δprg-orgΔinv-spa in this paper) was also generated using the lambda red system using the primers invG-fw and spaS-rev as well as prgH-fw and orgC–rev and the plasmids pKD3 and pKD4 to generate prgHIJKorgABC::aphT, invGEABCIJspaOPQRS::cat, a mutant lacking most genes of the TTS apparatus. For construction of strain M2532 (termed Δ8Δsip in this paper), we transduced the ÄsipBCDA-sptP::aphT allele from SB245 (SL1344, ÄsipBCDA-sptP::aphT fliGHI::Tn10; K. Kaniga and J. E. Galan, unpublished data) via P22 into M2400 (SL1344, ÄsopE, ÄsopE2, ÄsopB, ÄsipA, ÄsptP, ÄsopA, ÄspvB, ÄspvC), which has been previously described. M2532 fails to express most TTSS-1 effector proteins and the translocon components.

35To create the suicide plasmid pM2002, pVS152Tsr was digested with the restriction endonucleases Eco47III and XmaI. The tsrvenus encoding fragment was ligated into pM1300 (digested with MslI and XmaI, ) downstream of a truncated sipA fragment (nt 1156–2058 of the orf), to finally create pM2002 and introduced by homologous recombination into the genome of ATCC14028 to generate the reporter strain M2001. To obtain the tsrvenus reporter for hilA (M2076), the c-terminal region of hilA (nt 114 to 1661 of the orf) was amplified using the primer pair hilA-fw-XmaI-NcoI and hilA-rev-NheI-XbaI and cloned into pBluescriptII (Invitrogen) using the restriction endonucleases XmaI and XbaI, yielding pM2090. This plasmid was digested with NheI and NotI to introduce the tsrvenus encoding PCR fragment (template pM2002, primers: venus-NheI-fw and venus-NotI-rev, digested with NheI and NotI) to obtain pM2095. The entire region ranging from hilA to tsrvenus was cloned into pSB377 using the restriction enzymes NotI and XmaI yielding the suicide plasmid pM2080. This plasmid was used to generate the hilA reporter strain M2076 by homologous recombination into the genome of ATCC14028. To obtain the tsrvenus reporter for fliC, tsrvenus was amplified by PCR (primers: tsr-XmaI-fw and venus-XbaI-rev) and cloned into pBluescriptII using XmaI and XbaI thus yielding pM2533. After amplification of fliC by PCR using SL1344 chromosomal DNA as template and primers fliC-XhoI-fw and fliC-HindIII-rev, the fliC encoding fragment was cloned via XhoI and HindIII upstream of the tsrvenus gene into pM2533, thus yielding pM2539. Subsequently, the construct was moved via XhoI and XbaI into the suicide plasmid pGP704, thus yielding pM2819. This plasmid was used to create the fliC-tsrvenus reporter strain M2821 by homologous recombination into the genome of SL1344.

36All over-expression plasmids from pM2010 to pM2042 were obtained by digesting the indicated PCR fragments (Table S2 and S3 for plasmids and primers) with EcoRI and XbaI into pBAD24.

37All mutations were verified by PCR or DNA sequencing.

38HilA expression was analyzed by quantitative Western blot using an affinity-purified rabbit α-HilA antiserum (Fig. 1E). Recombinant HilA was used for normalization. SipC was detected using an α-SipC serum (Fig. 1G).

39For invasion, MDCK cells were grown in MEM (Invitrogen), infected for 30 min (MOI = 5;, washed and incubated in MEM (400 µg/ml gentamicin; 1 h). Intracellular bacteria were enumerated by plating.

40Prior to analysis, fluorophore formation was ensured (2 h, RT, 30 µg/ml chloramphenicol). Tsrvenus and Gfp emission was analyzed at 530 nm (supplement; FACSCalibur 4-color, Becton Dickinson). Bacteria were identified by side scatter (SSC). Data were analyzed with FlowJo software (Tree Star, Inc.). For Tsrvenus (Fig. 1), ln-transformed fluorescence values for 40000 events were median-normalized (subtraction) and compared to the similarly normalized data from the reporterless control strain, thus yielding the fraction of TTSS-1+ individuals. For sorting bacterial cells, S. Tm (psicA-gfp) cells were sorted by FACS (Aria Becton Dickinson, FACSDiva Software).

41Bacteria were placed on a 1.5% agarose pad equilibrated with LB, sealed under a glass coverslip and mounted (37°C temp. control; Axioplan2; Plan-APOCHROMAT 63x/1.4 oil; Zeiss or IX81, UPlanFLN 100x/1.3 Oil, Olympus). Reporter fluorescence (Exc. 470/20 nm; BP 495 nm; Em. 505–530 nm) and micro colony growth (phase contrast) were monitored and evaluated using Axiovision software (Zeiss). The slope of the ln-tranformed bacterial numbers (t), as determined from the logarithmic growth phase, yielded the growth rate µ. For sipA-tsrvenus and prgH-gfp, the micro colonies were scored visually as TTSS-1+ or TTSS-1−. To analyze differences in growth rates between TTSS-1+ and TTSS-1− micro colonies, we performed a full-factorial analysis of variance with the two factors phenotype (fixed) and experiment (random). Variance was analyzed in SPSS 17.0 (SPSS Inc. - Chicago, IL).

42Growth rates w/o reporter were analyzed via a maximum likelihood approach to test for two subpopulations with different growth rates. The growth rate measurements from five independent experiments (87 micro colonies) were combined. Using maximum likelihood, we fitted a bi-modal distribution (the sum of two normal probability density functions) and a unimodal (normal) distribution, and compared the two fits with a likelihood ratio test using R software.

43In Fig. 4, cell growth and ttss-1 expression were analyzed using a modified version of the cell tracking software described in. The first cell in each micro colony that could be observed over a whole division was used to analyze the statistical association between ttss-1 expression and the interval between two divisions (by non-parametric correlation analysis using PASW Statistics 18.0.0). 157 micro colonies were analyzed to estimate the fraction of micro colonies in which all cells, none of the cells, and a fraction of the cells expressed ttss-1. These groupings were based on visual inspection of each micro colony.