Phenotypic variation in bacteria: the role of feedback regulation

Evidence for phenotypic bistability resulting from transcriptional interference of bvgAS in Bordetella bronchiseptica

Bordetella species cause respiratory infections in mammals. Their master regulatory system BvgAS controls expression of at least three distinct phenotypic phases in response to environmental cues. The Bvg⁺ phase is necessary and sufficient for respiratory infection while the Bvg⁻ phase is required for survival ex vivo. We obtained large colony variants (LCVs) from the lungs of mice infected with B. bronchiseptica strain RBX9, which contains an in-frame deletion mutation in fhaB, encoding filamentous haemagglutinin. RBX9 also yielded LCVs when switched from Bvg⁻ phase conditions to Bvg⁺ phase conditions in vitro. We determined that LCVs are composed of both Bvg⁺ and Bvg⁻ phase bacteria and that they result from defective bvgAS positive autoregulation. The LCV phenotype was linked to the presence of a divergent promoter 5' to bvgAS, suggesting a previously undescribed mechanism of transcriptional interference that, in this case, leads to feedback-based bistability (FBM). Our results also indicate that a small proportion of RBX9 bacteria modulates to the Bvg⁻ phase in vivo. In addition to providing insight into transcriptional interference and FBM, our data provide an example of an in-frame deletion mutation exerting a 'polar' effect on nearby genes.

Evolutionary tuning of protein expression levels of a positively autoregulated two-component system

Cellular adaptation relies on the development of proper regulatory schemes for accurate control of gene expression levels in response to environmental cues. Over- or under-expression can lead to diminished cell fitness due to increased costs or insufficient benefits. Positive autoregulation is a common regulatory scheme that controls protein expression levels and gives rise to essential features in diverse signaling systems, yet its roles in cell fitness are less understood. It remains largely unknown how much protein expression is 'appropriate' for optimal cell fitness under specific extracellular conditions and how the dynamic environment shapes the regulatory scheme to reach appropriate expression levels. Here, we investigate the correlation of cell fitness and output response with protein expression levels of the E. coli PhoB/PhoR two-component system (TCS). In response to phosphate (Pi)-depletion, the PhoB/PhoR system activates genes involved in phosphorus assimilation as well as genes encoding themselves, similarly to many other positively autoregulated TCSs. We developed a bacteria competition assay in continuous cultures and discovered that different Pi conditions have conflicting requirements of protein expression levels for optimal cell fitness. Pi-replete conditions favored cells with low levels of PhoB/PhoR while Pi-deplete conditions selected for cells with high levels of PhoB/PhoR. These two levels matched PhoB/PhoR concentrations achieved via positive autoregulation in wild-type cells under Pi-replete and -deplete conditions, respectively. The fitness optimum correlates with the wild-type expression level, above which the phosphorylation output saturates, thus further increase in expression presumably provides no additional benefits. Laboratory evolution experiments further indicate that cells with non-ideal protein levels can evolve toward the optimal levels with diverse mutational strategies. Our results suggest that the natural protein expression levels and feedback regulatory schemes of TCSs are evolved to match the phosphorylation output of the system, which is determined by intrinsic activities of TCS proteins.

Noise-induced regime shifts: A quantitative characterization

Diverse complex dynamical systems are known to exhibit abrupt regime shifts at bifurcation points of the saddle-node type. The dynamics of most of these systems, however, have a stochastic component resulting in noise-driven regime shifts even if the system is away from the bifurcation points. In this paper, we propose a new quantitative measure, namely, the propensity transition point as an indicator of stochastic regime shifts. The concepts and the methodology are illustrated for the one-variable May model, a well-known model in ecology and the genetic toggle, a two-variable model of a simple genetic circuit. The general applicability and usefulness of the method for the analysis of regime shifts is further demonstrated in the case of the mycobacterial switch to persistence for which experimental data are available.

Genetic characterization of atypical Citrobacter freundii

The ability of a bacterial population to survive in different niches, as well as in stressful and rapidly changing environmental conditions, depends greatly on its genetic content. To survive such fluctuating conditions, bacteria have evolved different mechanisms to modulate phenotypic variations and related strategies to produce high levels of genetic diversity. Laboratories working in microbiological diagnosis have shown that Citrobacter freundii is very versatile in its colony morphology, as well as in its biochemical, antigenic and pathogenic behaviours. This phenotypic versatility has made C. freundii difficult to identify and it is frequently confused with both Salmonella enterica and Escherichia coli. In order to determine the genomic events and to explain the mechanisms involved in this plasticity, six C. freundii isolates were selected from a phenotypic variation study. An I-CeuI genomic cleavage map was created and eight housekeeping genes, including 16S rRNA, were sequenced. In general, the results showed a range of both phenotypes and genotypes among the isolates with some revealing a greater similarity to C. freundii and some to S. enterica, while others were identified as phenotypic and genotypic intermediary states between the two species. The occurrence of these events in natural populations may have important implications for genomic diversification in bacterial evolution, especially when considering bacterial species boundaries. In addition, such events may have a profound impact on medical science in terms of treatment, course and outcomes of infectious diseases, evading the immune response, and understanding host-pathogen interactions.

Extrinsic noise driven phenotype switching in a self-regulating gene

Analysis of complex gene regulation networks gives rise to a landscape of metastable phenotypic states for cells. Heterogeneity within a population arises due to infrequent noise-driven transitions of individual cells between nearby metastable states. While most previous works have focused on the role of intrinsic fluctuations in driving such transitions, in this Letter we investigate the role of extrinsic fluctuations. First, we develop an analytical framework to study the combined effect of intrinsic and extrinsic noise on a toy model of a Boolean regulated genetic switch. We then extend these ideas to a more biologically relevant model with a Hill-like regulatory function. Employing our theory and Monte Carlo simulations, we show that extrinsic noise can significantly alter the lifetimes of the phenotypic states and may fundamentally change the escape mechanism. Finally, our theory can be readily generalized to more complex decision making networks in biology.

A genome-wide transcriptional profiling of sporulating Bacillus subtilis strain lacking PrpE protein phosphatase

The sporulation process is a complex genetic developmental program leading to profound changes in global gene expression profile. In this work, we have applied genome-wide microarray approach for transcriptional profiling of Bacillus subtilis strain lacking a gene coding for PrpE protein phosphatase. This protein was previously shown to be involved in the regulation of germination of B. subtilis spores. Moreover, the deletion of prpE gene resulted in changing the resistance properties of spores. Our results provide genome-wide insight into the influence of this protein phosphatase on the physiology of B. subtilis cells. Although the precise role of PrpE in shaping the observed phenotype of ΔprpE mutant strain still remains beyond the understanding, our experiments brought observations of possible indirect implication of this protein in the regulation of cell motility and chemotaxis, as well as the development of competence. Surprisingly, prpE-deleted cells showed elevated level of general stress response, which turned out to be growth medium specific.

You are what you talk: quorum sensing induces individual morphologies and cell division modes in Dinoroseobacter shibae

Dinoroseobacter shibae, a member of the Roseobacter clade abundant in marine environments, is characterized by a pronounced pleomorphism. Cell shapes range from variable-sized ovoid rods to long filaments with a high copy number of chromosomes. Time-lapse microscopy shows cells dividing either by binary fission or by budding from the cell poles. Here we demonstrate that this morphological heterogeneity is induced by quorum sensing (QS). D. shibae utilizes three acylated homoserine lactone (AHL) synthases (luxI1-3) to produce AHLs with unsaturated C18 side chains. A ΔluxI1-knockout strain completely lacking AHL biosynthesis was uniform in morphology and divided by binary fission only. Transcriptome analysis revealed that expression of genes responsible for control of cell division was reduced in this strain, providing the link between QS and the observed phenotype. In addition, flagellar biosynthesis and type IV secretion system (T4SS) were downregulated. The wild-type phenotype and gene expression could be restored through addition of synthetic C18-AHLs. Their effectiveness was dependent on the number of double bonds in the acyl side chain and the regulated trait. The wild-type expression level of T4SS genes was fully restored even by an AHL with a saturated C18 side chain that has not been detected in D. shibae. QS induces phenotypic individualization of D. shibae cells rather than coordinating the population. This strategy might be beneficial in unpredictably changing environments, for example, during algal blooms when resource competition and grazing exert fluctuating selective pressures. A specific response towards non-native AHLs might provide D. shibae with the capacity for complex interspecies communication.

New approaches to the problem of generating coherent, reproducible phenotypes

Fundamental, unresolved questions in biology include how a bacterium generates coherent phenotypes, how a population of bacteria generates a coherent set of such phenotypes, how the cell cycle is regulated and how life arose. To try to help answer these questions, we have developed the concepts of hyperstructures, competitive coherence and life on the scales of equilibria. Hyperstructures are large assemblies of macromolecules that perform functions. Competitive coherence describes the way in which organisations such as cells select a subset of their constituents to be active in determining their behaviour; this selection results from a competition between a process that is responsible for a historical coherence and another process responsible for coherence with the current environment. Life on the scales of equilibria describes how bacteria depend on the cell cycle to negotiate phenotype space and, in particular, to satisfy the conflicting constraints of having to grow in favourable conditions so as to reproduce yet not grow in hostile conditions so as to survive. Both competitive coherence and life on the scales deal with the problem of reconciling conflicting constraints. Here, we bring together these concepts in the common framework of hyperstructures and make predictions that may be tested using a learning program, Coco, and secondary ion mass spectrometry.

In situ activity of a dominant Prochlorococcus ecotype (eHL-II) from rRNA content and cell size

In the open ocean genetically diverse clades of the unicellular cyanobacteria Prochlorococcus are biogeographically structured along environmental gradients, yet little is known about their in situ activity. To address this gap, here we use the numerically dominant Prochlorococcus clade eHL-II (eMIT9312) as a model organism to develop and apply a method to examine their in situ activity using rRNA content and cell size as metrics of cellular physiology. For two representative isolates (MIT9312 and MIT9215) rRNA cell(-1) increases linearly with specific growth rate but is anticorrelated with cell size indicated by flow cytometrically measured (SSC). Although each strain has a unique relationship between cellular rRNA (or cell size) and growth rate, both strains have the same strong positive correlation between rRNA cell(-1) SSC(-1) and growth rate. We field test this approach and observe distinct patterns of eHL-II clade specific activity (rRNA cell(-1) SSC(-1)) with depth that are consistent with patterns of photosynthetic rates. This molecular technique provides unique insight into the ecology of Prochlorococcus and could potentially be expanded to include other microbes to unravel the ecological and biogeochemical contributions of genetically distinct marine side scatter microbes.

Population dynamics of bacterial persistence

Persistence is a prime example of phenotypic heterogeneity, where a microbial population splits into two distinct subpopulations with different growth and survival properties as a result of reversible phenotype switching. Specifically, persister cells grow more slowly than normal cells under unstressed growth conditions, but survive longer under stress conditions such as the treatment with bactericidal antibiotics. We analyze the population dynamics of such a population for several typical experimental scenarios, namely a constant environment, shifts between growth and stress conditions, and periodically switching environments. We use an approximation scheme that allows us to map the dynamics to a logistic equation for the subpopulation ratio and derive explicit analytical expressions for observable quantities that can be used to extract underlying dynamic parameters from experimental data. Our results provide a theoretical underpinning for the study of phenotypic switching, in particular for organisms where detailed mechanistic knowledge is scarce.

Stochastic transitions in a bistable reaction system on the membrane

Transitions between steady states of a multi-stable stochastic system in the perfectly mixed chemical reactor are possible only because of stochastic switching. In realistic cellular conditions, where diffusion is limited, transitions between steady states can also follow from the propagation of travelling waves. Here, we study the interplay between the two modes of transition for a prototype bistable system of kinase-phosphatase interactions on the plasma membrane. Within microscopic kinetic Monte Carlo simulations on the hexagonal lattice, we observed that for finite diffusion the behaviour of the spatially extended system differs qualitatively from the behaviour of the same system in the well-mixed regime. Even when a small isolated subcompartment remains mostly inactive, the chemical travelling wave may propagate, leading to the activation of a larger compartment. The activating wave can be induced after a small subdomain is activated as a result of a stochastic fluctuation. Such a spontaneous onset of activity is radically more probable in subdomains characterized by slower diffusion. Our results show that a local immobilization of substrates can lead to the global activation of membrane proteins by the mechanism that involves stochastic fluctuations followed by the propagation of a semi-deterministic travelling wave.

Programmed heterogeneity: epigenetic mechanisms in bacteria

Contrary to the traditional view that bacterial populations are clonal, single-cell analysis reveals that phenotypic heterogeneity is common in bacteria. Formation of distinct bacterial lineages appears to be frequent during adaptation to harsh environments, including the colonization of animals by bacterial pathogens. Formation of bacterial subpopulations is often controlled by epigenetic mechanisms that generate inheritable phenotypic diversity without altering the DNA sequence. Such mechanisms are diverse, ranging from relatively simple feedback loops to complex self-perpetuating DNA methylation patterns.

Agent-based modeling of competence phenotype switching in Bacillus subtilis

Background

It is a fascinating phenomenon that in genetically identical bacteria populations of Bacillus subtilis, a distinct DNA uptake phenotype called the competence phenotype may emerge in 10–20% of the population. Many aspects of the phenomenon are believed to be due to the variable expression of critical genes: a stochastic occurrence termed “noise” which has made the phenomenon difficult to examine directly by lab experimentation.

Methods

To capture and model noise in this system and further understand the emergence of competence both at the intracellular and culture levels in B. subtilis, we developed a novel multi-scale, agent-based model. At the intracellular level, our model recreates the regulatory network involved in the competence phenotype. At the culture level, we simulated growth conditions, with our multi-scale model providing feedback between the two levels.

Results

Our model predicted three potential sources of genetic “noise”. First, the random spatial arrangement of molecules may influence the manifestation of the competence phenotype. In addition, the evidence suggests that there may be a type of epigenetic heritability to the emergence of competence, influenced by the molecular concentrations of key competence molecules inherited through cell division. Finally, the emergence of competence during the stationary phase may in part be due to the dilution effect of cell division upon protein concentrations.

Conclusions

The competence phenotype was easily translated into an agent-based model – one with the ability to illuminate complex cell behavior. Models such as the one described in this paper can simulate cell behavior that is otherwise unobservable in vivo, highlighting their potential usefulness as research tools.

Reverse vaccinology in the 21st century: improvements over the original design

Reverse vaccinology (RV), the first application of genomic technologies in vaccine research, represented a major revolution in the process of discovering novel vaccines. By determining their entire antigenic repertoire, researchers could identify protective targets and design efficacious vaccines for pathogens where conventional approaches had failed. Bexsero, the first vaccine developed using RV, has recently received positive opinion from the European Medicines Agency. The use of RV initiated a cascade of changes that affected the entire vaccine development process, shifting the focus from the identification of a list of vaccine candidates to the definition of a set of high throughput screens to reduce the need for costly and labor intensive tests in animal models. It is now clear that a deep understanding of the epidemiology of vaccine candidates, and their regulation and role in host-pathogen interactions, must become an integral component of the screening workflow. Far from being outdated by technological advancements, RV still represents a paradigm of how high-throughput technologies and scientific insight can be integrated into biotechnology research.

Slow protein fluctuations explain the emergence of growth phenotypes and persistence in clonal bacterial populations

One of the most challenging problems in microbiology is to understand how a small fraction of microbes that resists killing by antibiotics can emerge in a population of genetically identical cells, the phenomenon known as persistence or drug tolerance. Its characteristic signature is the biphasic kill curve, whereby microbes exposed to a bactericidal agent are initially killed very rapidly but then much more slowly. Here we relate this problem to the more general problem of understanding the emergence of distinct growth phenotypes in clonal populations. We address the problem mathematically by adopting the framework of the phenomenon of so-called weak ergodicity breaking, well known in dynamical physical systems, which we extend to the biological context. We show analytically and by direct stochastic simulations that distinct growth phenotypes can emerge as a consequence of slow-down of stochastic fluctuations in the expression of a gene controlling growth rate. In the regime of fast gene transcription, the system is ergodic, the growth rate distribution is unimodal, and accounts for one phenotype only. In contrast, at slow transcription and fast translation, weakly non-ergodic components emerge, the population distribution of growth rates becomes bimodal, and two distinct growth phenotypes are identified. When coupled to the well-established growth rate dependence of antibiotic killing, this model describes the observed fast and slow killing phases, and reproduces much of the phenomenology of bacterial persistence. The model has major implications for efforts to develop control strategies for persistent infections.

Sticking together: building a biofilm the Bacillus subtilis way

Biofilms are ubiquitous communities of tightly associated bacteria encased in an extracellular matrix. Bacillus subtilis has long served as a robust model organism to examine the molecular mechanisms of biofilm formation, and a number of studies have revealed that this process is regulated by several integrated pathways. In this Review, we focus on the molecular mechanisms that control B. subtilis biofilm assembly, and then briefly summarize the current state of knowledge regarding biofilm disassembly. We also discuss recent progress that has expanded our understanding of B. subtilis biofilm formation on plant roots, which are a natural habitat for this soil bacterium.

Stochasticity in colonial growth dynamics of individual bacterial cells

Conventional bacterial growth studies rely on large bacterial populations without considering the individual cells. Individual cells, however, can exhibit marked behavioral heterogeneity. Here, we present experimental observations on the colonial growth of 220 individual cells of Salmonella enterica serotype Typhimurium using time-lapse microscopy videos. We found a highly heterogeneous behavior. Some cells did not grow, showing filamentation or lysis before division. Cells that were able to grow and form microcolonies showed highly diverse growth dynamics. The quality of the videos allowed for counting the cells over time and estimating the kinetic parameters lag time (λ) and maximum specific growth rate (μmax) for each microcolony originating from a single cell. To interpret the observations, the variability of the kinetic parameters was characterized using appropriate probability distributions and introduced to a stochastic model that allows for taking into account heterogeneity using Monte Carlo simulation. The model provides stochastic growth curves demonstrating that growth of single cells or small microbial populations is a pool of events each one of which has its own probability to occur. Simulations of the model illustrated how the apparent variability in population growth gradually decreases with increasing initial population size (N(0)). For bacterial populations with N(0) of >100 cells, the variability is almost eliminated and the system seems to behave deterministically, even though the underlying law is stochastic. We also used the model to demonstrate the effect of the presence and extent of a nongrowing population fraction on the stochastic growth of bacterial populations.

Immediate and heterogeneous response of the LiaFSR two-component system of Bacillus subtilis to the peptide antibiotic bacitracin

BACKGROUND

Two-component signal transduction systems are one means of bacteria to respond to external stimuli. The LiaFSR two-component system of Bacillus subtilis consists of a regular two-component system LiaRS comprising the core Histidine Kinase (HK) LiaS and the Response Regulator (RR) LiaR and additionally the accessory protein LiaF, which acts as a negative regulator of LiaRS-dependent signal transduction. The complete LiaFSR system was shown to respond to various peptide antibiotics interfering with cell wall biosynthesis, including bacitracin.

METHODOLOGY AND PRINCIPAL FINDINGS

Here we study the response of the LiaFSR system to various concentrations of the peptide antibiotic bacitracin. Using quantitative fluorescence microscopy, we performed a whole population study analyzed on the single cell level. We investigated switching from the non-induced 'OFF' state into the bacitracin-induced 'ON' state by monitoring gene expression of a fluorescent reporter from the RR-regulated liaI promoter. We found that switching into the 'ON' state occurred within less than 20 min in a well-defined switching window, independent of the bacitracin concentration. The switching rate and the basal expression rate decreased at low bacitracin concentrations, establishing clear heterogeneity 60 min after bacitracin induction. Finally, we performed time-lapse microscopy of single cells confirming the quantitative response as obtained in the whole population analysis for high bacitracin concentrations.

CONCLUSION

The LiaFSR system exhibits an immediate, heterogeneous and graded response to the inducer bacitracin in the exponential growth phase.

Functional analysis of the ComK protein of Bacillus coagulans

The genes for DNA uptake and recombination in Bacilli are commonly regulated by the transcriptional factor ComK. We have identified a ComK homologue in Bacillus coagulans, an industrial relevant organism that is recalcitrant for transformation. Introduction of B. coagulans comK gene under its own promoter region into Bacillus subtilis comK strain results in low transcriptional induction of the late competence gene comGA, but lacking bistable expression. The promoter regions of B. coagulans comK and the comGA genes are recognized in B. subtilis and expression from these promoters is activated by B. subtilis ComK. Purified ComK protein of B. coagulans showed DNA-binding ability in gel retardation assays with B. subtilis- and B. coagulans-derived probes. These experiments suggest that the function of B. coagulans ComK is similar to that of ComK of B. subtilis. When its own comK is overexpressed in B. coagulans the comGA gene expression increases 40-fold, while the expression of another late competence gene, comC is not elevated and no reproducible DNA-uptake could be observed under these conditions. Our results demonstrate that B. coagulans ComK can recognize several B.subtilis comK-responsive elements, and vice versa, but indicate that the activation of the transcription of complete sets of genes coding for a putative DNA uptake apparatus in B. coagulans might differ from that of B. subtilis.

Expression of SCCmec cassette chromosome recombinases in methicillin-resistant Staphylococcus aureus and Staphylococcus epidermidis

OBJECTIVES

Methicillin resistance in staphylococci is mediated by the mecA gene, which is carried on the staphylococcal cassette chromosome mec (SCCmec). SCCmec is responsible for vertical and horizontal transfer of methicillin resistance. Horizontal transfer implies first SCCmec excision from the chromosome. Site-specific excision is catalysed by the Ccr recombinases, which are encoded by ccrAB genes located on the cassette. The aim of this study is to determine the promoter activity of ccrAB genes in individual cells of methicillin-resistant Staphylococcus aureus (N315, COL and MW2) and Staphylococcus epidermidis (RP62A). One mutant cured of its SCCmec (N315EX) was also used. Exposure to various stresses was included in the study.

METHODS

For each strain, translational promoter-green fluorescent protein (gfp) fusions were used to assess the levels of ccr promoter activity in individual cells. Analyses were performed using epifluorescence microscopy and flow cytometry.

RESULTS

ccr promoter activity was observed in only a small percentage of cell populations. This 'bistable' phenotype was strain dependent (GFP was expressed in N315 and RP62A, but not in COL and MW2) and growth dependent (GFP-expressing cells decreased from approximately 3% to 1% between logarithmic and stationary growth phases). The ccr promoter of strain N315 displayed normal promoter activity when expressed in SCCmec-negative N315EX. Likewise, the ccr promoter of strain COL (which was inactive in COL) showed normal N315-like activity when transformed into N315 and N315EX.

CONCLUSIONS

SCCmec excision operates through bistability, favouring a small fraction of cells to 'sacrifice' their genomic islands for transfer, while the rest of the population remains intact. Determinants responsible for the activity of the ccr promoter were not located on SCCmec, but were elsewhere on the genome. Thus, the staphylococcal chromosome plays a key role in determining SCCmec stability and transferability.

Ultrasensitivity of the Bacillus subtilis sporulation decision

Starving Bacillus subtilis cells execute a gene expression program resulting in the formation of stress-resistant spores. Sporulation master regulator, Spo0A, is activated by a phosphorelay and controls the expression of a multitude of genes, including the forespore-specific sigma factor σ(F) and the mother cell-specific sigma factor σ(E). Identification of the system-level mechanism of the sporulation decision is hindered by a lack of direct control over Spo0A activity. This limitation can be overcome by using a synthetic system in which Spo0A activation is controlled by inducing expression of phosphorelay kinase KinA. This induction results in a switch-like increase in the number of sporulating cells at a threshold of KinA. Using a combination of mathematical modeling and single-cell microscopy, we investigate the origin and physiological significance of this ultrasensitive threshold. The results indicate that the phosphorelay is unable to achieve a sufficiently fast and ultrasensitive response via its positive feedback architecture, suggesting that the sporulation decision is made downstream. In contrast, activation of σ(F) in the forespore and of σ(E) in the mother cell compartments occurs via a cascade of coherent feed-forward loops, and thereby can produce fast and ultrasensitive responses as a result of KinA induction. Unlike σ(F) activation, σ(E) activation in the mother cell compartment only occurs above the KinA threshold, resulting in completion of sporulation. Thus, ultrasensitive σ(E) activation explains the KinA threshold for sporulation induction. We therefore infer that under uncertain conditions, cells initiate sporulation but postpone making the sporulation decision to average stochastic fluctuations and to achieve a robust population response.

A microfluidic platform for rapid, stress-induced antibiotic susceptibility testing of Staphylococcus aureus

The emergence and spread of bacterial resistance to ever increasing classes of antibiotics intensifies the need for fast phenotype-based clinical tests for determining antibiotic susceptibility. Standard susceptibility testing relies on the passive observation of bacterial growth inhibition in the presence of antibiotics. In this paper, we present a novel microfluidic platform for antibiotic susceptibility testing based on stress-activation of biosynthetic pathways that are the primary targets of antibiotics. We chose Staphylococcus aureus (S. aureus) as a model system due to its clinical importance, and we selected bacterial cell wall biosynthesis as the primary target of both stress and antibiotic. Enzymatic and mechanical stresses were used to damage the bacterial cell wall, and a β-lactam antibiotic interfered with the repair process, resulting in rapid cell death of strains that harbor no resistance mechanism. In contrast, resistant bacteria remained viable under the assay conditions. Bacteria, covalently-bound to the bottom of the microfluidic channel, were subjected to mechanical shear stress created by flowing culture media through the microfluidic channel and to enzymatic stress with sub-inhibitory concentrations of the bactericidal agent lysostaphin. Bacterial cell death was monitored via fluorescence using the Sytox Green dead cell stain, and rates of killing were measured for the bacterial samples in the presence and absence of oxacillin. Using model susceptible (Sanger 476) and resistant (MW2) S. aureus strains, a metric was established to separate susceptible and resistant staphylococci based on normalized fluorescence values after 60 min of exposure to stress and antibiotic. Because this ground-breaking approach is not based on standard methodology, it circumvents the need for minimum inhibitory concentration (MIC) measurements and long wait times. We demonstrate the successful development of a rapid microfluidic-based and stress-activated antibiotic susceptibility test by correctly designating the phenotypes of 16 additional clinically relevant S. aureus strains in a blinded study. In addition to future clinical utility, this method has great potential for studying the effects of various stresses on bacteria and their antibiotic susceptibility.

Nanog-dependent feedback loops regulate murine embryonic stem cell heterogeneity

A number of key regulators of mouse embryonic stem (ES) cell identity, including the transcription factor Nanog, show strong expression fluctuations at the single cell level. The molecular basis for these fluctuations is unknown. Here we used a genetic complementation strategy to investigate expression changes during transient periods of Nanog downregulation. Employing an integrated approach, that includes high-throughput single cell transcriptional profiling and mathematical modelling, we found that early molecular changes subsequent to Nanog loss are stochastic and reversible. However, analysis also revealed that Nanog loss severely compromises the self-sustaining feedback structure of the ES cell regulatory network. Consequently, these nascent changes soon become consolidated to committed fate decisions in the prolonged absence of Nanog. Consistent with this, we found that exogenous regulation of Nanog-dependent feedback control mechanisms produced more a homogeneous ES cell population. Taken together our results indicate that Nanog-dependent feedback loops have a role in controlling both ES cell fate decisions and population variability.

Synthetic feedback loop model for increasing microbial biofuel production using a biosensor

Current biofuel production methods use engineered bacteria to break down cellulose and convert it to biofuel. A major challenge in microbial fuel production is that increasing biofuel yields can be limited by the toxicity of the biofuel to the organism that is producing it. Previous research has demonstrated that efflux pumps are effective at increasing tolerance to various biofuels. However, when overexpressed, efflux pumps burden cells, which hinders growth and slows biofuel production. Therefore, the toxicity of the biofuel must be balanced with the toxicity of pump overexpression. We have developed a mathematical model for cell growth and biofuel production that implements a synthetic feedback loop using a biosensor to control efflux pump expression. In this way, the production rate will be maximal when the concentration of biofuel is low because the cell does not expend energy expressing efflux pumps when they are not needed. Additionally, the microbe is able to adapt to toxic conditions by triggering the expression of efflux pumps, which allow it to continue biofuel production. Sensitivity analysis indicates that the feedback sensor model is insensitive to many system parameters, but a few key parameters can influence growth and production. In comparison to systems that express efflux pumps at a constant level, the feedback sensor increases overall biofuel production by delaying pump expression until it is needed. This result is more pronounced when model parameters are variable because the system can use feedback to adjust to the actual rate of biofuel production.

Analysis of two distinct mycelial populations in liquid-grown Streptomyces cultures using a flow cytometry-based proteomics approach

Streptomycetes are proficient producers of enzymes and antibiotics. When grown in bioreactors, these filamentous microorganisms form mycelial pellets that consist of interconnected hyphae. We here employed a flow cytometry approach designed for large particles (COPAS) and demonstrate that liquid-grown Streptomyces cultures consist of two distinct populations of pellets. One population consists of mycelia with a constant mean diameter of approximately 260 μm, whereas the other population contains larger mycelia whose diameter depends on the strain, the age of the culture, and medium composition. Quantitative proteomics analysis revealed that 37 proteins differed in abundance between the two populations of pellets. Stress-related proteins and biosynthetic proteins for production of the calcium-dependent antibiotic were more abundant in the population of large mycelia, while proteins involved in DNA topology, modification, or degradation were overrepresented in the population of small mycelia. Deletion of genes for the cellulose synthase-like protein CslA and the chaplins affected the average size of the population of large pellets but not that of small pellets. Considering the fact that the production of enzymes and metabolites depends on pellet size, these results provide new leads toward rational strain design of Streptomyces strains tailored for industrial fermentations.

Microfluidic study of competence regulation in Streptococcus mutans: environmental inputs modulate bimodal and unimodal expression of comX

Streptococcus mutans regulates genetic competence through a complex network that receives inputs from a number of environmental stimuli, including two signalling peptides designated as CSP and XIP. The response of the downstream competence genes to these inputs shows evidence of stochasticity and bistability and has been difficult to interpret. We have used microfluidic, single-cell methods to study how combinations of extracellular signals shape the response of comX, an alternative sigma factor governing expression of the late competence genes. We find that the composition of the medium determines which extracellular signal (XIP or CSP) can elicit a response from comX and whether that response is unimodal or bimodal across a population of cells. In a chemically defined medium, exogenous CSP does not induce comX, whereas exogenous XIP elicits a comX response from all cells. In complex medium, exogenous XIP does not induce comX, whereas CSP elicits a bimodal comX response from the population. Interestingly, bimodal behaviour required an intact copy of comS, which encodes the precursor of XIP. The comS-dependent capability for both unimodal and bimodal response suggests that a constituent - most likely peptides - of complex medium interacts with a positive feedback loop in the competence regulatory network.

How did metabolism and genetic replication get married?

In addressing the question of the origins of the relationship between metabolism and genetic replication, we consider the implications of a prebiotic, fission-fusion, ecology of composomes. We emphasise the importance of structures and non-specific catalysis on interfaces created by structures. From the assumption that the bells of the metabolism-replication wedding still echo in modern cells, we argue that the functional assemblies of macromolecules that constitute hyperstructures in modern bacteria are the descendants of composomes and that interactions at the hyperstructure level control the cell cycle. A better understanding of the cell cycle should help understand the original metabolism-replication marriage. This understanding requires new concepts such as metabolic signalling, metabolic sensing and Dualism, which entails the cells in a population varying the ratios of equilibrium to non-equilibrium hyperstructures so as to maximise the chances of both survival and growth. A deeper understanding of the coupling between metabolism and replication may also require a new view of cell cycle functions in creating a coherent diversity of phenotypes and in narrowing the combinatorial catalytic space. To take these ideas into account, we propose the Accordion model in which a dynamic interface between lipid domains catalysed monomer to polymer reactions and became decorated with peptides and nucleotides that favoured their own catalysis. In this model, metabolism, replication, differentiation and division all began together at the interface between extended equilibrium structures within protocells or composomes.

Origins of binary gene expression in post-transcriptional regulation by microRNAs

MicroRNA-mediated regulation of gene expression is characterised by some distinctive features that set it apart from unregulated and transcription-factor-regulated gene expression. Recently, a mathematical model has been proposed to describe the dynamics of post-transcriptional regulation by microRNAs. The model explains the observations made in single-cell experiments quite well. In this paper, we introduce some additional features into the model and consider two specific cases. In the first case, a noncooperative positive feedback loop is included in the transcriptional regulation of the target gene expression. In the second case, a stochastic version of the original model is considered in which there are random transitions between the inactive and active expression states of the gene. In the first case we show that bistability is possible in a parameter regime, due to the presence of a nonlinear protein decay term in the gene expression dynamics. In the second case, we derive the conditions for obtaining stochastic binary gene expression. We find that this type of gene expression is more favourable in the case of regulation by microRNAs as compared to the case of unregulated gene expression. The theoretical predictions relating to binary gene expression are experimentally testable.

Culture history and population heterogeneity as determinants of bacterial adaptation: the adaptomics of a single environmental transition

Diversity in adaptive responses is common within species and populations, especially when the heterogeneity of the frequently large populations found in environments is considered. By focusing on events in a single clonal population undergoing a single transition, we discuss how environmental cues and changes in growth rate initiate a multiplicity of adaptive pathways. Adaptation is a comprehensive process, and stochastic, regulatory, epigenetic, and mutational changes can contribute to fitness and overlap in timing and frequency. We identify culture history as a major determinant of both regulatory adaptations and microevolutionary change. Population history before a transition determines heterogeneities due to errors in translation, stochastic differences in regulation, the presence of aged, damaged, cheating, or dormant cells, and variations in intracellular metabolite or regulator concentrations. It matters whether bacteria come from dense, slow-growing, stressed, or structured states. Genotypic adaptations are history dependent due to variations in mutation supply, contingency gene changes, phase variation, lateral gene transfer, and genome amplifications. Phenotypic adaptations underpin genotypic changes in situations such as stress-induced mutagenesis or prophage induction or in biofilms to give a continuum of adaptive possibilities. Evolutionary selection additionally provides diverse adaptive outcomes in a single transition and generally does not result in single fitter types. The totality of heterogeneities in an adapting population increases the chance that at least some individuals meet immediate or future challenges. However, heterogeneity complicates the adaptomics of single transitions, and we propose that subpopulations will need to be integrated into future population biology and systems biology predictions of bacterial behavior.

Optimality principles in the regulation of metabolic networks

One of the challenging tasks in systems biology is to understand how molecular networks give rise to emergent functionality and whether universal design principles apply to molecular networks. To achieve this, the biophysical, evolutionary and physiological constraints that act on those networks need to be identified in addition to the characterisation of the molecular components and interactions. Then, the cellular "task" of the network-its function-should be identified. A network contributes to organismal fitness through its function. The premise is that the same functions are often implemented in different organisms by the same type of network; hence, the concept of design principles. In biology, due to the strong forces of selective pressure and natural selection, network functions can often be understood as the outcome of fitness optimisation. The hypothesis of fitness optimisation to understand the design of a network has proven to be a powerful strategy. Here, we outline the use of several optimisation principles applied to biological networks, with an emphasis on metabolic regulatory networks. We discuss the different objective functions and constraints that are considered and the kind of understanding that they provide.

Surface enhanced Raman spectroscopy (SERS) for the discrimination of Arthrobacter strains based on variations in cell surface composition

Surface enhanced Raman spectroscopy (SERS) is a rapid and highly sensitive spectroscopic technique that has the potential to measure chemical changes in bacterial cell surface in response to environmental changes. The objective of this study was to determine whether SERS had sufficient resolution to differentiate closely related bacteria within a genus grown on solid and liquid medium, and a single Arthrobacter strain grown in multiple chromate concentrations. Fourteen closely related Arthrobacter strains, based on their 16S rRNA gene sequences, were used in this study. After performing principal component analysis in conjunction with Linear Discriminant Analysis, we used a novel, adapted cross-validation method, which more faithfully models the classification of spectra. All fourteen strains could be classified with up to 97% accuracy. The hierarchical trees comparing SERS spectra from the liquid and solid media datasets were different. Additionally, hierarchical trees created from the Raman data were different from those obtained using 16S rRNA gene sequences (a phylogenetic measure). A single bacterial strain grown on solid media culture with three different chromate levels also showed significant spectral distinction at discrete points identified by the new Elastic Net regularized regression method demonstrating the ability of SERS to detect environmentally induced changes in cell surface composition. This study demonstrates that SERS is effective in distinguishing between a large number of very closely related Arthrobacter strains and could be a valuable tool for rapid monitoring and characterization of phenotypic variations in a single population in response to environmental conditions.

Regulatory consequences of gene translocation in bacteria

Gene translocations play an important role in the plasticity and evolution of bacterial genomes. In this study, we investigated the impact on gene regulation of three genome organizational features that can be altered by translocations: (i) chromosome position; (ii) gene orientation; and (iii) the distance between a target gene and its transcription factor gene (‘target-TF distance’). Specifically, we quantified the effect of these features on constitutive expression, transcription factor binding and/or gene expression noise using a synthetic network in Escherichia coli composed of a transcription factor (LacI repressor) and its target gene (yfp). Here we show that gene regulation is generally robust to changes in chromosome position, gene orientation and target-TF distance. The only demonstrable effect was that chromosome position alters constitutive expression, due to changes in gene copy number and local sequence effects, and that this determines maximum and minimum expression levels. The results were incorporated into a mathematical model which was used to quantitatively predict the responses of a simple gene network to gene translocations; the predictions were confirmed experimentally. In summary, gene translocation can modulate constitutive gene expression levels due to changes in chromosome position but it has minimal impact on other facets of gene regulation.

Equation-free analysis of two-component system signalling model reveals the emergence of co-existing phenotypes in the absence of multistationarity

Phenotypic differences of genetically identical cells under the same environmental conditions have been attributed to the inherent stochasticity of biochemical processes. Various mechanisms have been suggested, including the existence of alternative steady states in regulatory networks that are reached by means of stochastic fluctuations, long transient excursions from a stable state to an unstable excited state, and the switching on and off of a reaction network according to the availability of a constituent chemical species. Here we analyse a detailed stochastic kinetic model of two-component system signalling in bacteria, and show that alternative phenotypes emerge in the absence of these features. We perform a bifurcation analysis of deterministic reaction rate equations derived from the model, and find that they cannot reproduce the whole range of qualitative responses to external signals demonstrated by direct stochastic simulations. In particular, the mixed mode, where stochastic switching and a graded response are seen simultaneously, is absent. However, probabilistic and equation-free analyses of the stochastic model that calculate stationary states for the mean of an ensemble of stochastic trajectories reveal that slow transcription of either response regulator or histidine kinase leads to the coexistence of an approximate basal solution and a graded response that combine to produce the mixed mode, thus establishing its essential stochastic nature. The same techniques also show that stochasticity results in the observation of an all-or-none bistable response over a much wider range of external signals than would be expected on deterministic grounds. Thus we demonstrate the application of numerical equation-free methods to a detailed biochemical reaction network model, and show that it can provide new insight into the role of stochasticity in the emergence of phenotypic diversity.

It is a surprising fact that genetically identical bacteria, living in identical conditions, can develop in completely different ways: for example, one subpopulation might grow very fast and another very slowly. These different phenotypes are thought to be one reason why bacteria that cause disease can survive antibiotic treatment or become persistent. This diversity of behaviour is usually attributed to the existence of multiple stable phenotypic states, or to the coexistence of one stable state with another unstable excited state, or finally to the possibility of the whole biochemical system that controls the phenotype being switched on and off. In this paper we describe a different scenario that leads to phenotypic diversity in two-component system signalling, a very common mechanism that bacteria use to sense external signals and control their response to changes in their environment. We use probability theory and equation-free computational analysis to calculate the average number of molecules of each chemical species present in the two-component system and hence show that sporadic production of either of two key chemical components required for signalling can delay the response to the external signal in some bacterial cells and so lead to the emergence of two distinct cell populations.

Prokaryotic regulatory systems biology: Common principles governing the functional architectures of Bacillus subtilis and Escherichia coli unveiled by the natural decomposition approach

Escherichia coli and Bacillus subtilis are two of the best-studied prokaryotic model organisms. Previous analyses of their transcriptional regulatory networks have shown that they exhibit high plasticity during evolution and suggested that both converge to scale-free-like structures. Nevertheless, beyond this suggestion, no analyses have been carried out to identify the common systems-level components and principles governing these organisms. Here we show that these two phylogenetically distant organisms follow a set of common novel biologically consistent systems principles revealed by the mathematically and biologically founded natural decomposition approach. The discovered common functional architecture is a diamond-shaped, matryoshka-like, three-layer (coordination, processing, and integration) hierarchy exhibiting feedback, which is shaped by four systems-level components: global transcription factors (global TFs), locally autonomous modules, basal machinery and intermodular genes. The first mathematical criterion to identify global TFs, the κ-value, was reassessed on B. subtilis and confirmed its high predictive power by identifying all the previously reported, plus three potential, master regulators and eight sigma factors. The functionally conserved cores of modules, basal cell machinery, and a set of non-orthologous common physiological global responses were identified via both orthologous genes and non-orthologous conserved functions. This study reveals novel common systems principles maintained between two phylogenetically distant organisms and provides a comparison of their lifestyle adaptations. Our results shed new light on the systems-level principles and the fundamental functions required by bacteria to sustain life.

A disposable picolitre bioreactor for cultivation and investigation of industrially relevant bacteria on the single cell level

In the continuously growing field of industrial biotechnology the scale-up from lab to industrial scale is still a major hurdle to develop competitive bioprocesses. During scale-up the productivity of single cells might be affected by bioreactor inhomogeneity and population heterogeneity. Currently, these complex interactions are difficult to investigate. In this report, design, fabrication and operation of a disposable picolitre cultivation system is described, in which environmental conditions can be well controlled on a short time scale and bacterial microcolony growth experiments can be observed by time-lapse microscopy. Three exemplary investigations will be discussed emphasizing the applicability and versatility of the device. Growth and analysis of industrially relevant bacteria with single cell resolution (in particular Escherichia coli and Corynebacterium glutamicum) starting from one single mother cell to densely packed cultures is demonstrated. Applying the picolitre bioreactor, 1.5-fold increased growth rates of C. glutamicum wild type cells were observed compared to typical 1 litre lab-scale batch cultivation. Moreover, the device was used to analyse and quantify the morphological changes of an industrially relevant l-lysine producer C. glutamicum after artificially inducing starvation conditions. Instead of a one week lab-scale experiment, only 1 h was sufficient to reveal the same information. Furthermore, time lapse microscopy during 24 h picolitre cultivation of an arginine producing strain containing a genetically encoded fluorescence sensor disclosed time dependent single cell productivity and growth, which was not possible with conventional methods.

SlrA/SinR/SlrR inhibits motility gene expression upstream of a hypersensitive and hysteretic switch at the level of σ(D) in Bacillus subtilis

Exponentially growing Bacillus subtilis cultures are epigenetically differentiated into two subpopulations in which cells are either ON or OFF for σ(d) -dependent gene expression: a pattern suggestive of bistability. The gene encoding σ(D) , sigD, is part of the 31-gene fla/che operon where its location at the 3' end, 25 kb away from the strong P(fla/che) promoter, determines its expression level relative to a threshold. Here we show that addition of a single extra copy of the slrA gene in the chromosome inhibited σ(d) -dependent gene expression. SlrA together with SinR and SlrR reduced sigD transcript by potentiating a distance-dependent decrease in fla/che operon transcript abundance that was not mediated by changes in expression from the P(fla/che) promoter. Consistent with acting upstream of σ(D) , SlrA/SinR/SlrR was bypassed by artificial ectopic expression of sigD and hysteretically maintained for 20 generations by engaging the sigD gene at the native locus. SlrA/SinR/SlrR was also bypassed by increasing fla/che transcription and resulted in a hypersensitive output in flagellin expression. Thus, flagellin gene expression demonstrated hypersensitivity and hysteresis and we conclude that σ(d) -dependent gene expression is bistable.

Microfluidic capture and release of bacteria in a conical nanopore array

We present a simple and inexpensive method for the capture and release of bacteria contained in an array of conical nanopores on a membrane inside a microfluidic device. As an example, we demonstrate that cyanobacteria can be captured, one bacterium per pore, in a defined orientation with over 500 bacteria per membrane with viabilities as high as 100%. The device can also specifically capture cyanobacteria from a mixed suspension of cyanobacteria and chlamydomonas with a selectivity as high as 90%.

Emergent bistability: effects of additive and multiplicative noise

Positive feedback and cooperativity in the regulation of gene expression are generally considered to be necessary for obtaining bistable expression states. Recently, a novel mechanism of bistability termed emergent bistability has been proposed which involves only positive feedback and no cooperativity in the regulation. An additional positive feedback loop is effectively generated due to the inhibition of cellular growth by the synthesized proteins. The mechanism, demonstrated for a synthetic circuit, may be prevalent in natural systems also as some recent experimental results appear to suggest. In this paper, we study the effects of additive and multiplicative noise on the dynamics governing emergent bistability. The calculational scheme employed is based on the Langevin and Fokker-Planck formalisms. The steady state probability distributions of protein levels and the mean first passage times are computed for different noise strengths and system parameters. In the region of bistability, the bimodal probability distribution is shown to be a linear combination of a lognormal and a Gaussian distribution. The variances of the individual distributions and the relative weights of the distributions are further calculated for varying noise strengths and system parameters. The experimental relevance of the model results is also pointed out.

Influence of molecular noise on the growth of single cells and bacterial populations

During the last decades experimental studies have revealed that single cells of a growing bacterial population are significantly exposed to molecular noise. Important sources for noise are low levels of metabolites and enzymes that cause significant statistical variations in the outcome of biochemical reactions. In this way molecular noise affects biological processes such as nutrient uptake, chemotactic tumbling behavior, or gene expression of genetically identical cells. These processes give rise to significant cell-to-cell variations of many directly observable quantities such as protein levels, cell sizes or individual doubling times. In this study we theoretically explore if there are evolutionary benefits of noise for a growing population of bacteria. We analyze different situations where noise is either suppressed or where it affects single cell behavior. We consider two specific examples that have been experimentally observed in wild-type Escherichia coli cells: (i) the precision of division site placement (at which molecular noise is highly suppressed) and (ii) the occurrence of noise-induced phenotypic variations in fluctuating environments. Surprisingly, our analysis reveals that in these specific situations both regulatory schemes [i.e. suppression of noise in example (i) and allowance of noise in example (ii)] do not lead to an increased growth rate of the population. Assuming that the observed regulatory schemes are indeed caused by the presence of noise our findings indicate that the evolutionary benefits of noise are more subtle than a simple growth advantage for a bacterial population in nutrient rich conditions.

Cell individuality: the bistable gene expression of the type III secretion system in Dickeya dadantii 3937

Dickeya dadantii 3937 is a gram-negative phytopathogenic bacterium that expresses genes encoding a type III secretion system (T3SS) in a bistable pattern when cultured in a homogeneous minimal media. In this work, we further characterized the bistable gene expression of T3SS at the single-cell level. We demonstrated that bistable expression of the HrpL-regulon genes, such as hrpA and hrpN, is controlled by the same regulatory mechanism. We also showed that the expression level of the T3SS master regulatory gene hrpL plays an important role in the development of the bistable expression of hrpA. A high expression level of hrpL is required but unable to guarantee the high-state expression of hrpA in a cell. In addition, bistable expression patterns of T3SS genes in other gram-negative pathogens of the Enterobacteriaceae and Pseudomonadaceae families were also described in this study. This suggests that the T3SS bistability might be a conserved population behavior in several gram-negative bacterial pathogens.

Phenotypic variation and host interactions of Xenorhabdus bovienii SS-2004, the entomopathogenic symbiont of Steinernema jollieti nematodes

Xenorhabdus bovienii (SS-2004) bacteria reside in the intestine of the infective-juvenile (IJ) stage of the entomopathogenic nematode, Steinernema jollieti. The recent sequencing of the X. bovienii genome facilitates its use as a model to understand host - symbiont interactions. To provide a biological foundation for such studies, we characterized X. bovienii in vitro and host interaction phenotypes. Within the nematode host X. bovienii was contained within a membrane bound envelope that also enclosed the nematode-derived intravesicular structure. Steinernema jollieti nematodes cultivated on mixed lawns of X. bovienii expressing green or DsRed fluorescent proteins were predominantly colonized by one or the other strain, suggesting the colonizing population is founded by a few cells. Xenorhabdus bovienii exhibits phenotypic variation between orange-pigmented primary form and cream-pigmented secondary form. Each form can colonize IJ nematodes when cultured in vitro on agar. However, IJs did not develop or emerge from Galleria mellonella insects infected with secondary form. Unlike primary-form infected insects that were soft and flexible, secondary-form infected insects retained a rigid exoskeleton structure. Xenorhabdus bovienii primary and secondary form isolates are virulent towards Manduca sexta and several other insects. However, primary form stocks present attenuated virulence, suggesting that X. bovienii, like Xenorhabdus nematophila may undergo virulence modulation.

Regulation of phenotypic heterogeneity permits Salmonella evasion of the host caspase-1 inflammatory response

Sensing and adapting to the environment is one strategy by which bacteria attempt to maximize fitness in an unpredictable world; another is the stochastic generation of phenotypically distinct subgroups within a genetically clonal population. In culture, Salmonella Typhimurium populations are bistable for the expression of flagellin. We report that YdiV controls this expression pattern by preventing transcription of the sigma factor that recruits RNA polymerase to the flagellin promoter. Bistability ensues when the sigma factor is repressed in a subpopulation of cells, resulting in two phenotypes: flagellin expressors and flagellin nonexpressors. Although the ability to swim is presumably a critical survival trait, flagellin activates eukaryotic defense pathways, and Salmonella restrict the production of flagellin during systemic infection. Salmonella mutants lacking YdiV are unable to fully repress flagellin at systemic sites, rendering them vulnerable to caspase-1 mediated colonization restriction. Thus, a regulatory mechanism producing bistability also impacts Salmonella virulence.

Cross-strain quorum sensing inhibition by Staphylococcus aureus. Part 2: A spatially inhomogeneous model

Staphylococcus aureus uses quorum sensing (QS) to enhance its pathogenicity. An intriguing aspect of this is that different strains are capable of inactivating the QS systems of opposing strains. In Part 1 of this study, we presented a model of this phenomenon in a well-mixed environment; here, we incorporate spatial structure. Two competitive strains occupying adjacent habitats with freely diffusing QS signal molecules (QSSMs) are considered. We investigate the effect of the QSSM diffusion coefficient and the relative size of the two populations on the ability of one population to dominate the other. Regarding population size, a larger population is generally at an advantage (initial conditions permitting), while the implications of different diffusivities are more complex and depend upon the sizes of the populations.