Spo0A∼P Imposes a Temporal Gate for the Bimodal Expression of Competence in Bacillus subtilis

ComK transcriptionally controls competence for the uptake of transforming DNA in Bacillus subtilis. Only 10%–20% of the cells in a clonal population are randomly selected for competence. Because ComK activates its own promoter, cells exceeding a threshold amount of ComK trigger a positive feedback loop, transitioning to the competence ON state. The transition rate increases to a maximum during the approach to stationary phase and then decreases, with most cells remaining OFF. The average basal rate of comK transcription increases transiently, defining a window of opportunity for transitions and accounting for the heterogeneity of competent populations. We show that as the concentration of the response regulator Spo0A∼P increases during the entry to stationary phase it first induces comK promoter activity and then represses it by direct binding. Spo0A∼P activates by antagonizing the repressor, Rok. This amplifies an inherent increase in basal level comK promoter activity that takes place during the approach to stationary phase and is a general feature of core promoters, serving to couple the probability of competence transitions to growth rate. Competence transitions are thus regulated by growth rate and temporally controlled by the complex mechanisms that govern the formation of Spo0A∼P. On the level of individual cells, the fate-determining noise for competence is intrinsic to the comK promoter. This overall mechanism has been stochastically simulated and shown to be plausible. Thus, a deterministic mechanism modulates an inherently stochastic process

Improved statistical analysis of low abundance phenomena in bimodal bacterial populations.

Accurate detection of subpopulation size determinations in bimodal populations remains problematic yet it represents a powerful way by which cellular heterogeneity under different environmental conditions can be compared. So far, most studies have relied on qualitative descriptions of population distribution patterns, on population-independent descriptors, or on arbitrary placement of thresholds distinguishing biological ON from OFF states. We found that all these methods fall short of accurately describing small population sizes in bimodal populations. Here we propose a simple, statistics-based method for the analysis of small subpopulation sizes for use in the free software environment R and test this method on real as well as simulated data. Four so-called population splitting methods were designed with different algorithms that can estimate subpopulation sizes from bimodal populations. All four methods proved more precise than previously used methods when analyzing subpopulation sizes of transfer competent cells arising in populations of the bacterium Pseudomonas knackmussii B13. The methods' resolving powers were further explored by bootstrapping and simulations. Two of the methods were not severely limited by the proportions of subpopulations they could estimate correctly, but the two others only allowed accurate subpopulation quantification when this amounted to less than 25% of the total population. In contrast, only one method was still sufficiently accurate with subpopulations smaller than 1% of the total population. This study proposes a number of rational approximations to quantifying small subpopulations and offers an easy-to-use protocol for their implementation in the open source statistical software environment R