Bistability and phase variation in Salmonella enterica

Cell-to-cell differences in bacterial gene expression can merely reflect the occurrence of noise. In certain cases, however, heterogeneous gene expression is a programmed event that results in bistable expression. If bistability is heritable, bacterial lineages are formed. When programmed bistability is reversible, the phenomenon is known as phase variation. In certain cases, bistability is controlled by genetic mechanisms (e. g., DNA rearrangement). In other cases, bistability has epigenetic origin. A robust epigenetic mechanism for the formation of bacterial lineages is the formation of heritable DNA methylation patterns. However, bistability can also arise upon propagation of gene expression patterns by feedback loops that are stable upon cell division. This review describes examples of bistability and phase variation in Salmonella enterica and discusses their adaptive value, sometimes in a speculative manner.     

A DNA polymerase V homologue encoded by TOL plasmid pWW0 confers evolutionary fitness on Pseudomonas putida under conditions of environmental stress

Plasmids in conjunction with other mobile elements such as transposons are major players in the genetic adaptation of bacteria in response to changes in environment. Here we show that a large catabolic TOL plasmid, pWW0, from Pseudomonas putida carries genes (rulAB genes) encoding an error-prone DNA polymerase Pol V homologue which increase the survival of bacteria under conditions of accumulation of DNA damage. A study of population dynamics in stationary phase revealed that the presence of pWW0-derived rulAB genes in the bacterial genome allows the expression of a strong growth advantage in stationary phase (GASP) phenotype of P. putida. When rulAB-carrying cells from an 8-day-old culture were mixed with Pol V-negative cells from a 1-day-old culture, cells derived from the aged culture out-competed cells from the nonaged culture and overtook the whole culture. At the same time, bacteria from an aged culture lacking the rulAB genes were only partially able to out-compete cells from a fresh overnight culture of the parental P. putida strain. Thus, in addition to conferring resistance to DNA damage, the plasmid-encoded Pol V genes significantly increase the evolutionary fitness of bacteria during prolonged nutritional starvation of a P. putida population. The results of our study indicate that RecA is involved in the control of expression of the pWW0-encoded Pol V.     

Temperature-regulated expression of bacterial virulence genes

Virulence gene expression in most bacteria is a highly regulated phenomenon, affected by a variety of parameters including osmolarity, pH, ion concentration, iron levels, growth phase, and population density. Virulence genes are also regulated by temperature, which acts as an 'on-off' switch in a manner distinct from the more general heat-shock response. Here, we review temperature-responsive expression of virulence genes in four diverse pathogens.     

Positive feedback and noise activate the stringent response regulator rel in mycobacteria

Phenotypic heterogeneity in an isogenic, microbial population enables a subset of the population to persist under stress. In mycobacteria, stresses like nutrient and oxygen deprivation activate the stress response pathway involving the two-component system MprAB and the sigma factor, SigE. SigE in turn activates the expression of the stringent response regulator, rel. The enzyme polyphosphate kinase 1 (PPK1) regulates this pathway by synthesizing polyphosphate required for the activation of MprB. The precise manner in which only a subpopulation of bacterial cells develops persistence, remains unknown. Rel is required for mycobacterial persistence. Here we show that the distribution of rel expression levels in a growing population of mycobacteria is bimodal with two distinct peaks corresponding to low (L) and high (H) expression states, and further establish that a positive feedback loop involving the mprAB operon along with stochastic gene expression are responsible for the phenotypic heterogeneity. Combining single cell analysis by flow cytometry with theoretical modeling, we observe that during growth, noise-driven transitions take a subpopulation of cells from the L to the H state within a "window of opportunity" in time preceding the stationary phase. It is these cells which adapt to nutrient depletion in the stationary phase via the stringent response. We find evidence of hysteresis in the expression of rel in response to changing concentrations of PPK1. Hysteresis promotes robustness in the maintenance of the induced state. Our results provide, for the first time, evidence that bistability and stochastic gene expression could be important for the development of "heterogeneity with an advantage" in mycobacteria and suggest strategies for tackling tuberculosis like targeting transitions from the low to the high rel expression state.     

Bistability versus bimodal distributions in gene regulatory processes from population balance

In recent times, stochastic treatments of gene regulatory processes have appeared in the literature in which a cell exposed to a signaling molecule in its environment triggers the synthesis of a specific protein through a network of intracellular reactions. The stochastic nature of this process leads to a distribution of protein levels in a population of cells as determined by a Fokker-Planck equation. Often instability occurs as a consequence of two (stable) steady state protein levels, one at the low end representing the "off" state, and the other at the high end representing the "on" state for a given concentration of the signaling molecule within a suitable range. A consequence of such bistability has been the appearance of bimodal distributions indicating two different populations, one in the "off" state and the other in the "on" state. The bimodal distribution can come about from stochastic analysis of a single cell. However, the concerted action of the population altering the extracellular concentration in the environment of individual cells and hence their behavior can only be accomplished by an appropriate population balance model which accounts for the reciprocal effects of interaction between the population and its environment. In this study, we show how to formulate a population balance model in which stochastic gene expression in individual cells is incorporated. Interestingly, the simulation of the model shows that bistability is neither sufficient nor necessary for bimodal distributions in a population. The original notion of linking bistability with bimodal distribution from single cell stochastic model is therefore only a special consequence of a population balance model.     

In silico evolved lac operons exhibit bistability for artificial inducers, but not for lactose

Bistability in the lac operon of Escherichia coli has been widely studied, both experimentally and theoretically. Experimentally, bistability has been observed when E. coli is induced by an artificial, nonmetabolizable, inducer. However, if the lac operon is induced with lactose, the natural inducer, bistability has not been demonstrated. We derive an analytical expression that can predict the occurrence of bistability both for artificial inducers and lactose. We find very different conditions for bistability in the two cases. Indeed, for artificial inducers bistability is predicted, but for lactose the condition for bistability is much more difficult to satisfy. Moreover, we demonstrate that in silico evolution of the lac operon generates an operon that avoids bistability with respect to lactose, but does exhibit bistability with respect to artificial inducers. The activity of this evolved operon strikingly resembles the experimentally observed activity of the operon. Thus our computational experiments suggest that the wild-type lac operon, which regulates lactose metabolism, is not a bistable switch. Nevertheless, for engineering purposes, this operon can be used as a bistable switch with artificial inducers.     

Heterogeneity in expression of the <Emphasis Type="Italic">Escherichia coli</Emphasis> colicin K activity gene <Emphasis Type="Italic">cka</Emphasis> is controlled by the SOS system and stochastic factors

Phenotypic diversity provides populations of prokaryotic and eukaryotic organisms with the flexibility required to adapt to and/or survive environmental perturbations. Consequently, there is much interest in unraveling the molecular mechanisms of heterogeneity. A classical example of heterogeneity in Escherichia coli is the subset (3%) of the population that expresses the colicin K activity gene (cka) upon nutrient starvation. Here, we report on the mechanism underlying this variable response. As colicin synthesis is regulated by the LexA protein, the central regulator of the SOS response, we focused on the role of LexA and the SOS system in the variable cka expression. Real-time RT-PCR showed that the SOS system, without exogenous DNA damage, induces moderate levels of cka expression. The use of cka-gfp fusions demonstrated that modification of the conserved LexA boxes in the cka promoter region affected LexA binding affinity and the percentage of cka-gfp expressing cells in the population. A lexA-gfp fusion showed that the lexA gene is highly expressed in a subset of bacteria. Furthermore, cka-gfp fusions cloned into higher copy plasmid vectors increased the percentage of cka-gfp positive bacteria. Together, these results indicate that the bistability in cka expression in the bacterial population is determined by (1) basal SOS activity, (2) stochastic factors and possibly (3) the interplay of LexA dimers at cka operator. Other LexA regulated processes could exhibit similar regulation.     

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.     

The expression of antibiotic resistance genes in antibiotic‐producing bacteria

Antibiotic‐producing bacteria encode antibiotic resistance genes that protect them from the biologically active molecules that they produce. The expression of these genes needs to occur in a timely manner: either in advance of or concomitantly with biosynthesis. It appears that there have been at least two general solutions to this problem. In many cases, the expression of resistance genes is tightly linked to that of antibiotic biosynthetic genes. In others, the resistance genes can be induced by their cognate antibiotics or by intermediate molecules from their biosynthetic pathways. The regulatory mechanisms that couple resistance to antibiotic biosynthesis are mechanistically diverse and potentially relevant to the origins of clinical antibiotic resistance.     

Quorum sensing and communication in bacteria

Bacteria are able “to sense” an increase in the cell population density and to respond to it by the induction of special sets of genes. This type of regulation, called Quorum Sensing (QS), includes the production and excretion of low-molecular-weight signaling molecules (autoinducers, AI), which diffuse readily through the cell wall, from cells into the medium. As the bacterial population reaches the critical level of density, the concentration of these signaling molecules in the medium increases as a function of population density. On reaching the critical threshold concentration, AIs bind to specific receptor regulatory proteins, which induce the expression of target genes. By means of AIs, bacteria accomplish the communication that is the transmission of information between bacteria belonging to the same or different species, genera, and even families: the signaling molecules of some bacteria affect the receptors of others causing a coordinated reply of cells of the bacterial population. Bacteria of different taxonomic groups use the QS systems in regulation of a broad range of physiological activities. These processes include virulence, symbiosis, conjugation, biofilm formation, bioluminescence, synthesis of enzymes, antibiotic substances, etc. Here we review different QS systems of bacteria, the role of QS in bacterial communication, and some applied aspects of QS regulation application.     

Intra-Tumour Signalling Entropy Determines Clinical Outcome in Breast and Lung Cancer

The cancer stem cell hypothesis, that a small population of tumour cells are responsible for tumorigenesis and cancer progression, is becoming widely accepted and recent evidence has suggested a prognostic and predictive role for such cells. Intra-tumour heterogeneity, the diversity of the cancer cell population within the tumour of an individual patient, is related to cancer stem cells and is also considered a potential prognostic indicator in oncology. The measurement of cancer stem cell abundance and intra-tumour heterogeneity in a clinically relevant manner however, currently presents a challenge. Here we propose signalling entropy, a measure of signalling pathway promiscuity derived from a sample’s genome-wide gene expression profile, as an estimate of the stemness of a tumour sample. By considering over 500 mixtures of diverse cellular expression profiles, we reveal that signalling entropy also associates with intra-tumour heterogeneity. By analysing 3668 breast cancer and 1692 lung adenocarcinoma samples, we further demonstrate that signalling entropy correlates negatively with survival, outperforming leading clinical gene expression based prognostic tools. Signalling entropy is found to be a general prognostic measure, valid in different breast cancer clinical subgroups, as well as within stage I lung adenocarcinoma. We find that its prognostic power is driven by genes involved in cancer stem cells and treatment resistance. In summary, by approximating both stemness and intra-tumour heterogeneity, signalling entropy provides a powerful prognostic measure across different epithelial cancers.     

The plant cell defense and Agrobacterium tumefaciens

We previously identified changes in gene expression in Ageratum conyzoides plant cells inoculated with Agrobacterium tumefaciens by using cDNA-AFLP. Here, we show that a subset of defense-related genes is differentially regulated by an Agrobacterium attachment-deficient mutant. The expression pattern triggered by this mutant is similar to that induced by inoculation with non-pathogenic bacteria. We also observed that the expression level of the defense genes was inversely correlated with the efficiency of transformation by Agrobacterium. We propose that the plant defense system has an important role in controlling infection and transformation and that Agrobacterium may dampen some plant defense responses.     

Subpopulations of sensorless bacteria drive fitness in fluctuating environments

Populations of bacteria often undergo a lag in growth when switching conditions. Because growth lags can be large compared to typical doubling times, variations in growth lag are an important but often overlooked component of bacterial fitness in fluctuating environments. We here explore how growth lag variation is determined for the archetypical switch from glucose to lactose as a carbon source in Escherichia coli. First, we show that single-cell lags are bimodally distributed and controlled by a single-molecule trigger. That is, gene expression noise causes the population before the switch to divide into subpopulations with zero and nonzero lac operon expression. While “sensorless” cells with zero preexisting lac expression at the switch have long lags because they are unable to sense the lactose signal, any nonzero lac operon expression suffices to ensure a short lag. Second, we show that the growth lag at the population level depends crucially on the fraction of sensorless cells and that this fraction in turn depends sensitively on the growth condition before the switch. Consequently, even small changes in basal expression can significantly affect the fraction of sensorless cells, thereby population lags and fitness under switching conditions, and may thus be subject to significant natural selection. Indeed, we show that condition-dependent population lags vary across wild E. coli isolates. Since many sensory genes are naturally low expressed in conditions where their inducer is not present, bimodal responses due to subpopulations of sensorless cells may be a general mechanism inducing phenotypic heterogeneity and controlling population lags in switching environments. This mechanism also illustrates how gene expression noise can turn even a simple sensory gene circuit into a bet hedging module and underlines the profound role of gene expression noise in regulatory responses.

Is ignorance bliss for some bacterial cells? Single-cell analysis of the archetypical switch from glucose to lactose as a carbon source in E. coli shows that bacteria can exhibit stochastic bimodal responses to external stimuli because the corresponding sensory circuit is so lowly expressed that some cells are effectively blind to the stimulus.     

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.     

The Burkholderia bcpAIOB Genes Define Unique Classes of Two-Partner Secretion and Contact Dependent Growth Inhibition Systems

Microbes have evolved many strategies to adapt to changes in environmental conditions and population structures, including cooperation and competition. One apparently competitive mechanism is contact dependent growth inhibition (CDI). Identified in Escherichia coli, CDI is mediated by Two-Partner Secretion (TPS) pathway proteins, CdiA and CdiB. Upon cell contact, the toxic C-terminus of the TpsA family member CdiA, called the CdiA-CT, inhibits the growth of CDI(-) bacteria. CDI(+) bacteria are protected from autoinhibition by an immunity protein, CdiI. Bioinformatic analyses indicate that CDI systems are widespread amongst α, β, and γ proteobacteria and that the CdiA-CTs and CdiI proteins are highly variable. CdiI proteins protect against CDI in an allele-specific manner. Here we identify predicted CDI system-encoding loci in species of Burkholderia, Ralstonia and Cupriavidus, named bcpAIOB, that are distinguished from previously-described CDI systems by gene order and the presence of a small ORF, bcpO, located 5' to the gene encoding the TpsB family member. A requirement for bcpO in function of BcpA (the TpsA family member) was demonstrated, indicating that bcpAIOB define a novel class of TPS system. Using fluorescence microscopy and flow cytometry, we show that these genes are expressed in a probabilistic manner during culture of Burkholderia thailandensis in liquid medium. The bcpAIOB genes and extracellular DNA were required for autoaggregation and adherence to an abiotic surface, suggesting that CDI is required for biofilm formation, an activity not previously attributed to CDI. By contrast to what has been observed in E. coli, the B. thailandensis bcpAIOB genes only mediated interbacterial competition on a solid surface. Competition occurred in a defined spatiotemporal manner and was abrogated by allele-specific immunity. Our data indicate that the bcpAIOB genes encode distinct classes of CDI and TPS systems that appear to function in sociomicrobiological community development.