Email updates

Keep up to date with the latest news and content from BMC Biology and BioMed Central.

Journal App

google play app store

This article is part of the series Beyond Mendel: modeling in biology.

Open Access Research article

Division in Escherichia coli is triggered by a size-sensing rather than a timing mechanism

Lydia Robert12*, Marc Hoffmann3, Nathalie Krell4, Stéphane Aymerich12, Jérôme Robert5 and Marie Doumic67

Author Affiliations

1 INRA, Micalis CNRS-UMR 1319, 78350 Jouy-en-Josas, France

2 AgroParisTech, Micalis CNRS-UMR 1319, 78350 Jouy-en-Josas, France

3 Paris Dauphine University, CNRS-UMR 7534, Place du maréchal De Lattre de Tassigny, 75775 Paris, France

4 University of Rennes 1, CNRS-UMR 6625, Campus de Beaulieu, 35042 Rennes, France

5 CNRS/UPMC University of Paris 6, FRE 3231, Laboratoire Jean Perrin LJP, 75005 Paris, France

6 INRIA Paris-Rocquencourt, Domaine de Voluceau, BP 105, 781153 Le Chesnay, France

7 UPMC University of Paris 6, JL Lions Lab., 4 place Jussieu, 75005 Paris, France

For all author emails, please log on.

BMC Biology 2014, 12:17  doi:10.1186/1741-7007-12-17

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1741-7007/12/17


Received:31 December 2013
Accepted:24 February 2014
Published:28 February 2014

© 2014 Robert et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Many organisms coordinate cell growth and division through size control mechanisms: cells must reach a critical size to trigger a cell cycle event. Bacterial division is often assumed to be controlled in this way, but experimental evidence to support this assumption is still lacking. Theoretical arguments show that size control is required to maintain size homeostasis in the case of exponential growth of individual cells. Nevertheless, if the growth law deviates slightly from exponential for very small cells, homeostasis can be maintained with a simple ‘timer’ triggering division. Therefore, deciding whether division control in bacteria relies on a ‘timer’ or ‘sizer’ mechanism requires quantitative comparisons between models and data.

Results

The timer and sizer hypotheses find a natural expression in models based on partial differential equations. Here we test these models with recent data on single-cell growth of Escherichia coli. We demonstrate that a size-independent timer mechanism for division control, though theoretically possible, is quantitatively incompatible with the data and extremely sensitive to slight variations in the growth law. In contrast, a sizer model is robust and fits the data well. In addition, we tested the effect of variability in individual growth rates and noise in septum positioning and found that size control is robust to this phenotypic noise.

Conclusions

Confrontations between cell cycle models and data usually suffer from a lack of high-quality data and suitable statistical estimation techniques. Here we overcome these limitations by using high precision measurements of tens of thousands of single bacterial cells combined with recent statistical inference methods to estimate the division rate within the models. We therefore provide the first precise quantitative assessment of different cell cycle models.

Keywords:
Cell cycle; Bacteria; Division; Size control; Structured population equations; Numerical simulations; Nonparametric estimation

Background

Coordination between cell growth and division is often carried out by ‘size control’ mechanisms, where the cell size has to reach a certain threshold to trigger some event of the cell cycle, such as DNA replication or cell division [1]. As an example, the fission yeast Schizosaccharomyces pombe exhibits a size threshold at mitosis [2,3]. The budding yeast Saccharomyces cerevisiae also uses a size control mechanism that acts at the G1-S transition [4,5]. In contrast, in some cells such as those of early frog embryos, progression in the cell cycle is size independent and relies on a ‘timer’ mechanism [6].

Bacterial division is often assumed to be under size control but conclusive experimental evidence is still lacking and the wealth of accumulated data presents a complex picture. In 1968, building on the seminal work of Schaechter et al. and Helmstetter and Cooper, Donachie suggested that initiation of DNA replication is triggered when the bacterium reaches a critical size [7-9]. This provided the basis for a long-standing model of size control where cell size triggers replication initiation, which in turn determines the timing of division (see [10] and references therein). However, the coupling of replication initiation to cell mass has been repeatedly challenged [11-13]. In particular, on the basis of recent single-cell analysis, the team headed by N Kleckner proposed that replication initiation is more tightly linked to the time elapsed since birth than to cell mass [13,14]. In addition, the extent to which initiation timing affects division timing is unclear. In particular, variations in initiation timing are known to lead to compensatory changes in the duration of chromosome replication (see [15-17] and references therein). These studies argue against a size control model based on replication initiation. Another model postulates that size control acts directly on septum formation [18,19]. Nevertheless, the nature of the signals triggering the formation of the septal ring and its subsequent constriction are still unknown [17,20] and no molecular mechanism is known to sense cell size and transmit the information to the division machinery in bacteria.

Besides the work of Donachie, the assumption of size control in bacteria originates from a theoretical argument stating that such a control is necessary in exponentially growing cells to ensure cell size homeostasis, i.e. to maintain a constant size distribution through successive cycles. The growth of bacterial populations has long been mathematically described using partial differential equation (PDE) models. These models rely on hypotheses on division control: the division rate of a cell, i.e. the instantaneous probability of its dividing, can be assumed to depend either on cell age (i.e. the time elapsed since birth) or cell size. In the classical ‘sizer’ model, the division rate depends on size and not on age whereas in the ‘timer’ model it depends on age and not on size. Mathematical analysis of these models sheds light on the role of size control in cell size homeostasis. In particular, it has been suggested that for exponentially growing cells, a timer mechanism cannot ensure a stable size distribution [21,22]. Nevertheless, this unrealistic behavior of the timer mechanism is based on a biologically meaningless assumption, namely the exponential growth of cells of infinitely small or large size [23,24]. Cells of size zero or infinity do not exist and particularly small or large cells are likely to exhibit abnormal growth behavior. In conclusion, the mathematical arguments that were previously developed are insufficient to rule out a size-independent, timer model of bacterial division: quantitative comparisons between models and data are needed.

In the present study, we test whether age (i.e. the time elapsed since birth) or size is a determinant of cell division in E. coli. To do so, we analyzed two datasets derived from two major single-cell experimental studies on E. coli growth, performed by Stewart et al.[25] and Wang et al.[26]. Our analysis is based on division rate estimation by state-of-the-art nonparametric inference methods that we recently developed [27,28]. The two datasets correspond to different experimental setups and image analysis methods but lead to similar conclusions. We show that even though a model with a simple timer triggering division is sufficient to maintain cell size homeostasis, such a model is not compatible with the data. In addition, our analysis of the timer model shows that this model is very sensitive to hypotheses regarding the growth law of rare cells of very small or large size. This lack of robustness argues against a timer mechanism for division control in E. coli as well as in other exponentially growing organisms. In contrast, a model where cell size determines the probability of division is in good agreement with experimental data. Unlike the timer model, this sizer model is robust to slight modifications of the growth law of individual cells. In addition, our analysis reveals that the sizer model is very robust to phenotypic variability in individual growth rates or noise in septum positioning.

Results and discussion

Description of the data

Age and size distribution of the bacterial population

The results reported in this study were obtained from the analysis of two different datasets, obtained through microscopic time-lapse imaging of single E. coli cells growing in a rich medium, by Stewart et al. [25] and Wang et al.[26]. Stewart et al. followed single E. coli cells growing into microcolonies on LB-agarose pads at 30°C. The length of each cell in the microcolony was measured every 2 min. Wang et al. grew cells in LB: Luria Bertani medium at 37°C in a microfluidic setup [26] and the length of the cells was measured every minute. Due to the microfluidic device structure, at each division only one daughter cell could be followed (data si: sparse tree), in contrast to the experiment of Stewart et al. where all the individuals of a genealogical tree were followed (data fi: full tree). It is worth noting that the different structures of the data fi and si lead to different PDE models, and the statistical analysis was adapted to each situation (see below and Additional file 1). From each dataset (fi and si) we extracted the results of three experiments (experiments f1,f2 and f3 and s1,s2 and s3). Each experiment fi corresponds to the growth of approximately six microcolonies of up to approximately 600 cells and each experiment si to the growth of bacteria in 100 microchannels for approximately 40 generations.

Additional file 1. Supplementary text and figures.

Format: PDF Size: 1.1MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Given the accuracy of image analysis, we do not take into account variations of cell width within the population, which are negligible compared to cell-cycle-induced length variations. Thus, in the present study we do not distinguish between length, volume and mass and use the term cell size as a catch-all descriptor. Cell age and cell size distributions of a representative experiment from each dataset are shown in Figure 1. These distributions are estimated from the age and size measurements of every cell at every time step of a given experiment fi or si, by using a simple kernel density estimation method (kernel estimation is closely related to histogram construction but gives smooth estimates of distributions, as shown in Figure 1, for instance; for details see the Methods and Additional file 1). As expected for the different data structures (full tree fi or sparse tree si) and different experimental conditions, the distributions for the two datasets are not identical. The age distribution is decreasing with a maximum for age zero and the size distribution is wide and positively skewed, in agreement with previous results using various bacterial models [29-31].

thumbnailFigure 1. Distributions of cell age and cell size. Cell age (A) and cell size (B) distributions for a representative experiment of the fi dataset from Stewart et al. (green) [25] and of the si dataset from Wang et al.[26] (blue).

Testing the timer versus sizer models of division

Age-structured (timer) and size-structured (sizer) models

The timer and sizer hypotheses are easily expressed in mathematical terms: two different PDE models are commonly used to describe bacterial growth, using a division rate (i.e. the instantaneous probability of division) depending either on cell age or cell size. In the age-structured model (Age Model) the division rate Ba is a function only of the age a of the cell. The density n(t,a) of cells of age a at time t is given as a solution to the Mckendrick–Von Foerster equation (see [32] and references therein):

<a onClick="popup('http://www.biomedcentral.com/1741-7007/12/17/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1741-7007/12/17/mathml/M1">View MathML</a>

(1)

with the boundary condition

<a onClick="popup('http://www.biomedcentral.com/1741-7007/12/17/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1741-7007/12/17/mathml/M2">View MathML</a>

In this model, a cell of age a at time t has the probability Ba(a)dt of dividing between time t and t+dt.

In the size-structured model (Size model), the division rate Bs is a function only of the size x of the cell. Assuming that the size of a single cell grows with a rate v(x), the density n(t,x) of cells of size x at time t is given as a solution to the size-structured cell division equation: [32]

<a onClick="popup('http://www.biomedcentral.com/1741-7007/12/17/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1741-7007/12/17/mathml/M3">View MathML</a>

(2)

In the Size Model, a cell of size x at time t has the probability Bs(x)dt of dividing between time t and t+dt. This model is related to the so-called sloppy size control model [33] describing division in S. pombe.

For simplicity, we focused here on a population evolving along a full genealogical tree, accounting for fi data. For data si observed along a single line of descendants, an appropriate modification is made to Equations (1) and (2) (see Additional file 1: Supplementary Text).

Testing the Age Model (timer) and the Size Model (sizer) with experimental data

In this study we tested the hypothesis of an age-dependent versus size-dependent division rate by comparing the ability of the Age Model and Size Model to describe experimental data. The PDE given by Equations (1) and (2) can be embedded into a two-dimensional age-and-size-structured equation (Age & Size Model), describing the temporal evolution of the density n(t,a,x) of cells of age a and size x at time t, with a division rate Ba,sa priori depending on both age and size:

<a onClick="popup('http://www.biomedcentral.com/1741-7007/12/17/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1741-7007/12/17/mathml/M4">View MathML</a>

(3)

with the boundary condition

<a onClick="popup('http://www.biomedcentral.com/1741-7007/12/17/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1741-7007/12/17/mathml/M5">View MathML</a>

In this augmented setting, the Age Model governed by the PDE (1) and the Size Model governed by (2) are restrictions to the hypotheses of an age-dependent or size-dependent division rate, respectively (Ba,s=Ba or Ba,s=Bs).

The density n(t,a,x) of cells having age a and size x at a large time t can be approximated as n(t,a,x)≈eλtN(a,x), where the coefficient λ>0 is called the Malthus coefficient and N(a,x) is the stable age-size distribution. This regime is rapidly reached and time can then be eliminated from Equations (1), (2) and (3), which are thus transformed into equations governing the stable distribution N(a,x). Importantly, in the timer model (i.e. Ba,s=Ba), the existence of this stable distribution requires that growth is sub-exponential around zero and infinity [23,24].

We estimate the division rate Ba of the Age Model using the age measurements of every cell at every time step. Likewise, we estimate the division rate Bs of the Size Model using the size measurements of every cell at every time step. Our estimation procedure is based on mathematical methods we recently developed. Importantly, our estimation procedure does not impose any particular restrictions on the form of the division rate function B, so that any biologically realistic function can be estimated (see Additional file 1: Section 4 and Figure S6). In Additional file 1: Figures S1 and S2, we show the size-dependent and age-dependent division rates Bs(x) and Ba(a) estimated from the experimental data. Once the division rate has been estimated, the stable age and size distribution N(a,x) can be reconstructed through simulation of the Age & Size Model (using the experimentally measured growth rate; for details see the Methods).

We measure the goodness-of-fit of a model (timer or sizer) by estimating the distance between two distributions: the age-size distribution obtained through simulations of the model with the estimated division rate (as explained above), and the experimental age-size distribution. Therefore, a small distance indicates a good fit of the model to the experimental data. To estimate this distance we use a classical metric, which measures the average of the squared difference between the two distributions. As an example, the distance between two bivariate Gaussian distributions with the same mean and a standard deviation difference of 10% is 17%, and a 25% difference in standard deviation leads to a 50% distance between the distributions. The experimental age-size distribution is estimated from the age and size measurements of every cell at every time step of a given experiment fi or si, thanks to a simple kernel density estimation method.

Analysis of single-cell growth

As mentioned above, to avoid unrealistic asymptotic behavior of the Age Model and ensure the existence of a stable size distribution, assumptions have to be made on the growth of very small and large cells, which cannot be exactly exponential. To set realistic assumptions, we first studied the growth of individual cells. As expected, we found that during growth, a cell diameter is roughly constant (see inset in Figure 2A). Figure 2A shows cell length as a function of time for a representative cell, suggesting that growth is exponential rather than linear, in agreement with previous studies [25,26,34-36]. To test this hypothesis further, we performed linear and exponential fits of cell length for each single cell. We then calculated in each case the R2 coefficient of determination, which is classically used to measure how well a regression curve approximates the data (a perfect fit would give R2=1 and lower values indicate a poorer fit). The inset of Figure 2B shows the distribution of the R2 coefficient for all single cells for exponential (red) and linear (green) regressions, demonstrating that the exponential growth model fits the data very well and outperforms the linear growth model. We then investigated whether the growth of cells of particularly small or large size is exponential. If growth is exponential, the increase in length between each measurement should be proportional to the length. Therefore, we averaged the length increase of cells of similar size and tested whether the proportionality was respected for all sizes. As shown in Figure 2B, growth is exponential around the mean cell size but the behavior of very small or large cells may deviate from exponential growth. We therefore determined two size thresholds xmin and xmax below and over which the growth law may not be exponential (e.g. for the experiment f1 shown in Figure 2B, we defined xmin=2.3 µm and xmax=5.3 µm).

thumbnailFigure 2. Analysis of single-cell growth.(A) Cell length vs cell age for a representative cell (black dots); exponential fit (red curve) and linear fit (black line). Inset: Cell width vs cell age for the same cell. (B) Increase in cell length during one time step (i.e. 1 min) as a function of cell length for fi data. During the lifetime of a cell, the cell length is measured at each time step and the increase in cell length between successive time steps is calculated. Black dots are the average length increase for every cell of a given experiment f1, as a function of cell length; error bars are the average +/−2 SEM (standard error of the mean). The red line is a linear fit for lengths between 2.5 µm and 4.5 µm. Inset: For each single cell of f1, the evolution of cell length with age was fitted with a linear or an exponential function (as shown in panel A). We thus obtain a distribution of R2 coefficients corresponding to the linear (green) and exponential (red) fits.

The age-size joint distribution of E. coli corresponds to a size-dependent division rate

We used both the Age Model and Size Model to fit the experimental age-size distributions, following the approach described above. The growth law below xmin and above xmax is unknown. Therefore, to test the Age Model, growth was assumed to be exponential between xmin and xmax and we tested several growth functions v(x) for x<xmin and x>xmax, such as constant (i.e. linear growth) and polynomial functions. Figure 3 shows the best fit we could obtain. Comparing the experimental data f1 shown in Figure 3A (Figure 3B for s1 data) with the reconstructed distribution shown in Figure 3C (Figure 3D for s1 data) we can see that the Age Model fails to reconstruct the experimental age-size distribution and produces a distribution with a different shape. In particular, its localization along the y-axis is very different. For instance, for f1 data (panels A and C), the red area corresponding to the maximum of the experimental distribution is around 2.4 on the y-axis whereas the maximum of the fitted distribution is around 3.9. The y-axis corresponds to cell size. The size distribution produced by the Age Model is thus very different from the size distribution of the experimental data (experimental and fitted size distributions are shown in Additional file 1: Figure S9).

thumbnailFigure 3. Experimental and reconstructed age-size distributions for representative experiments from Stewartet al.[25](f1) and Wanget al.[26](s1).(A,B) Experimental age-size distributions for representative experiments f1 (A) and s1 (B). The frequency of cells of age a and size s in the population is represented by the color at the point of coordinate a on the x-axis and s on the y-axis, according to the scale indicated to the right of the figure. (C,D) Reconstruction of the distributions using the Age Model (C: reconstruction of the data f1 shown in panel A; D: reconstruction of the data s1 shown in panel B). These reconstructed distributions were obtained from simulations with the Age Model using a division rate estimated from the data (C: from f1, D: from s1). The growth functions used for the simulations are detailed in the Methods section. (E,F) Reconstruction of the distributions using the Size Model (E: reconstruction of the data f1 shown in panel A; F: reconstruction of the data s1 shown in panel B). These distributions were obtained from simulations with the Size Model using a division rate estimated from the data (E: from f1, F: from s1) with an exponential growth function (see Methods).

As an additional analysis to strengthen our conclusion, we calculated the correlation between the age at division and the size at birth using the experimental data. If division is triggered by a timer mechanism, these two variables should not be correlated, whereas we found a significant correlation of −0.5 both for si and fi data (P<10−16; see Additional file 1: Figure S7).

We used various growth functions for x<xmin and x>xmax but a satisfying fit could not be obtained with the Age Model. In addition, we found that the results of the Age Model are very sensitive to the assumptions made for the growth law of rare cells of very small and large size (see Additional file 1: Figure S3). This ultra-sensitivity to hypotheses regarding rare cells makes the timer model unrealistic generally for any exponentially growing organisms.

In contrast, the Size Model is in good agreement with the data (Figure 3: A compared to E and B compared to F) and allows a satisfactory reconstruction of the age-size structure of the population. The shape of the experimental and fitted distributions as well as their localization along the y-axis and x-axis are similar (size distributions and age distributions, i.e. projections onto the y-axis and x-axis, are shown in Additional file 1: Figure S8).

The quantitative measure of goodness-of-fit defined above is coherent with the curves’ visual aspects: for the Size Model the distance between the model and the data ranges from 17% to 20% for fi data (16% to 26% for si data) whereas for the Age Model it ranges from 51% to 93% for fi data (45% to 125% for si).

The experimental data has a limited precision. In particular, the division time is difficult to determine precisely by image analysis and the resolution is limited by the time step of image acquisition (for si and fi data, the time step represents respectively 5% and 8% of the average division time). By performing stochastic simulations of the Size Model (detailed in Additional file 1: Section 6), we evaluated the effect of measurement noise on the goodness of fit of the Size Model. We found that noise of 10% in the determination of the division time leads to a distance around 14%, which is of the order of the value obtained with our experimental data. We conclude that the Size Model fits the experimental data well. Moreover, we found that in contrast to the Age Model, the Size Model is robust with respect to the mathematical assumptions for the growth law for small and large sizes: the distance changes by less than 5%.

Size control is robust to phenotypic noise

Noise in the biochemical processes underlying growth and division, such as that created by stochastic gene expression, may perturb the control of size and affect the distribution of cell size. We therefore investigated the robustness of size control to such phenotypic noise. The Size Model describes the growth of a population of cells with variable age and size at division. Nevertheless, it does not take into account potential variability in individual growth rate or the difference in size at birth between two sister cells, i.e. the variability in septum positioning. To do so, we derived two PDE models, which are revised Size Models with either growth rate or septum positioning variability (see Additional file 1: Supplementary Text) and ran these models with different levels of variability.

Variability in individual growth rate has a negligible effect on the size distribution

For each single cell, a growth rate can be defined as the rate of the exponential increase of cell length with time [25,26]. By doing so, we obtained the distribution of the growth rate for the bacterial population (Additional file 1: Figure S4A). In our dataset this distribution is statistically compatible with a Gaussian distribution with a coefficient of variation of approximately 8% (standard deviation/mean =0.08).

We recently extended the Size Model to describe the growth of a population with single-cell growth rate variability (the equation is given in Additional file 1: Section 5) [28]. We simulated this extended Size Model using the growth rate distribution of fi data. The resulting size distribution is virtually identical to the one obtained without growth rate variability (Figure 4A, red and blue lines). Therefore, the naturally occurring variability in individual growth rate does not significantly perturb the size control. To investigate the effect of growth rate variability further, we simulated the model with various levels of noise, using truncated Gaussian growth rate distributions with coefficients of variation from 5 to 60%. We found that to obtain a 10% change in size distribution, a 30% coefficient of variation is necessary, which would represent an extremely high level of noise (Figure 4A, inset).

thumbnailFigure 4. Influence of the variability in individual growth rate and septum positioning on the cell size distribution.(A) Size distributions simulated using the Size Model with the division rate Bs estimated from f1 data and an exponential growth (v(x)=vx). In blue: Simulations with the same growth rate v=0.0274 for every cell. Red dashed line: Simulations with individual growth rates distributed according to the experimentally observed distribution. Green line: Growth rates normally distributed with coefficient of variation CV =60%. Inset: Difference (i.e. normalized integrated squared error) between the size distribution simulated without variability and the distributions simulated with various levels of variability (normally distributed individual growth rates of CV between 10% and 60%). (B) Simulated size distributions using the Size Model with the same division rate Bs as in A and a constant growth rate v=0.0274. In blue: Simulations where division is perfectly symmetrical. Red dashed line: Simulations with a variable septum position distributed according to the experimentally observed distribution. Green line: Simulations with a normally distributed septum position with CV =30%. Inset: Difference between the size distribution simulated without variability in the septum position and the distributions simulated with various levels of variability (normally distributed septum position of CV between 5% and 30%). CV, coefficient of variation.

Variability in septum positioning has a negligible effect on size distribution

The cells divide into two daughter cells of almost identical length. Nevertheless, a slight asymmetry can arise as an effect of noise during septum positioning. We found a 4% variation in the position of the septum (Additional file 1: Figure S4B), which is in agreement with previous measurements [35,37-39]. To test the robustness of size control to noise in septum positioning, we extended the Size Model to allow for different sizes of the two sister cells at birth (the equation is given in Additional file 1: Section 5). We ran this model using the empirical variability in septum positioning (shown in Additional file 1: Figure S4B) and compared the resulting size distribution to the one obtained by simulations without variability. As shown in Figure 4B (comparing the red and blue lines), the effect of natural noise in septum positioning is negligible. We also ran the model with higher levels of noise in septum positioning and found that a three times higher (12%) coefficient of variation is necessary to obtain a 10% change in size distribution (Figure 4B inset, and Additional file 1: Figure S5).

Conclusions

In the present study, we present statistical evidence to support the hypothesis that a size-dependent division rate can be used to reconstruct the experimental age-size distribution of E. coli. In contrast, this distribution cannot be generated by a timer model where the division rate depends solely on age. Even though the timer model can maintain cell size homeostasis, it is quantitatively incompatible with the observed size distribution. Our analysis of two different datasets shows the robustness of our conclusions to changes in experimental setup and image analysis methods. Our results therefore confirm the hypothesis of size control of division in E. coli. In addition, our analysis of the timer model shows that it is very sensitive to mathematical assumptions for the growth law of very rare cells of abnormal size, suggesting that this model is unrealistic for any exponentially growing organisms.

Noise in biochemical processes, in particular gene expression, can have a significant effect on the precision of biological circuits. In particular, it can generate a substantial variability in the cell cycle [5]. Therefore we investigated in bacteria the robustness of size control to noise in the single-cell growth rate and septum positioning, using appropriate extensions of the Size Model. We found that variability of the order of what we estimated from E. coli data does not significantly perturb the distribution of cell size. Therefore, in a natural population exhibiting phenotypic noise, the control of cell size is robust to fluctuations in septum positioning and individual growth rates. From a modeling perspective, this demonstrates that the simple Size Model is appropriate for describing a natural bacterial population showing phenotypic diversity.

Our approach is based on comparisons between PDE models and single-cell data for the cell cycle. Such comparisons were attempted a few decades ago using data from yeasts (e.g. [21,33]). Nevertheless, these interesting studies were hampered by the scarcity and poor quality of single-cell data as well as the lack of appropriate statistical procedures to estimate the division rate within the models. In contrast, we used high-precision measurements of tens of thousands cells in combination with modern statistical inference methods, which allowed us to assess quantitatively the adequacy of different models. We think this approach could prove successful in studying other aspects of the cell cycle, such as the coordination between replication and division or the molecular mechanisms underlying size control of division. Several different mechanisms involved in division control in bacteria have already been unraveled, in particular MinCD inhibition and nucleoid occlusion [40-42]. We believe that a better understanding of the relative roles played by MinCD inhibition and nucleoid occlusion in division control can be gained by analyzing the age-size distributions of minCD and nucleoid occlusion mutants. We are therefore currently performing time-lapse microscopy experiments to record the growth of such mutants.

Methods

Data analysis

The data of Stewart et al. contain the results of several experiments performed on different days, each of them recording the simultaneous growth of several microcolonies of the MG1655 E. coli strain on LB-agar pads at 30°C, with a generation time of approximately 26 min [25]. The first 150 min of growth were discarded to limit the effects of non-steady-state growth (cells undergo a slight plating stress when put on microscopy slides and it takes several generations to recover a stable growth rate). For the dataset obtained by Wang et al., the MG1655 E. coli strain was grown in LB at 37°C in a microfluidic device with a doubling time of approximately 20 min. To avoid any effect of replicative aging such as described in [26], we only kept the first 50 generations of growth. In addition the first ten generations were discarded to ensure steady-state growth. Both datasets were generated by analyzing fluorescent images (the bacteria express the Yellow Fluorescent Protein) using two different software systems. For si data, cell segmentation was based on the localization of brightness minima along the channel direction (see [26]). In the same spirit, for fi data, local minima of fluorescence intensity were used to outline the cells, following by an erosion and dilation step to separate adjacent cells (see [25]). To measure its length, a cell was approximated by a rectangle with the same second moments of pixel intensity and location distribution (for curved cells the measurement was done manually).

For both datasets we extracted data from three experiments done on different days. We did not pool the data together to avoid statistical biases arising from day-to-day differences in experimental conditions. Each analysis was performed in parallel on the data corresponding to each experiment.

Numerical simulations and estimation procedures

All the estimation procedures and simulations were performed using MATLAB. Experimental age-size distributions, such as those shown in Figure 3A,B, were estimated from the size and age measurements of every cell at every time step using the MATLAB kde2D function, which estimates the bivariate kernel density. This estimation was performed on a regular grid composed of 27 equally spaced points on [ 0,Amax] and 27 equally spaced points on [ 0,Xmax], where Amax is the maximal cell age in the data and Xmax the maximal cell size (for instance Amax=60 min and Xmax=10 µm for the experiment f1, as shown in Figure 3A). To estimate the size-dependent division rate Bs for each experiment, the distribution of size at division was first estimated for the cell size grid [ 0,Xmax] using the ksdensity function. This estimated distribution was then used to estimate Bs for the size grid using Equation (20) (for si data) or (22) (for fi data) of Additional file 1. The age-size distributions corresponding to the Size Model (Figure 3E,F) were produced by running the Age & Size Model (Equation (3) in the main text) using the estimated division rate Bs and an exponential growth function (v(x)=vx) with a rate v directly estimated from the data as the average of single-cell growth rates in the population (e.g. v=0.0274 min −1 for the f1 experiment and v=0.0317 min −1 for s1). For the Age & Size Model, we discretized the equation along the grid [ 0,Amax] and [ 0,Xmax], using an upwind finite volume method described in detail in [43]. We used a time step:

<a onClick="popup('http://www.biomedcentral.com/1741-7007/12/17/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1741-7007/12/17/mathml/M11">View MathML</a>

meeting the CFL: Courant-Friedrichs-Lewy stability criterion. We simulated n(t,a,x) iteratively until the age-size distribution reached stability (|(n(t+dt,a,x)−n(t,a,x))|<10−8). To eliminate the Malthusian parameter, the solution n(t,a,x) was renormalized at each time step (for details see [43]).

The age-dependent division rate Ba for each experiment was estimated for the cell age grid [0,Amax] using Equation (14) and (16) of Additional file 1. Using this estimated division rate, the age-size distributions corresponding to the Age Model (Figure 3C,D) were produced by running the Age & Size Model. As explained in the main text, we used various growth functions for small and large cells (i.e. for x<xmin and x>xmax; between xmin and xmax growth is exponential with the same rate as for the Size Model). For instance for the fit of the experiment f1 shown in Figure 3C, for x<2.3 µm and x>5.3 µm, v(x)= max(p(x),0), with p(x)=−0.0033x3+0.036x2−0.094x+0.13. Likewise, for the fit of the experiment s1 shown in Figure 3D, for x<3.5 µm and x>7.2 µm, v(x)= max(p(x),0), with p(x)=−0.0036x3+0.063x2−0.33x+0.67. For each dataset the polynomial p(x) was chosen as an interpolation of the function giving the length increase as a function of length (shown in Figure 2B for f1 data).

Simulations of the extended size models with variability in growth rates or septum positioning (Equations (23) and (24) in Additional file 1) were performed as for the Age & Size Model, with an upwind finite volume scheme. To simulate Equation (23), we used a grid composed of 27 equally spaced points on [ 0,Xmax] and 100 equally spaced points on [ 0.9vmin,1.1vmax], where vmin and vmax are the minimal and maximal individual growth rates in the data.

Abbreviations

PDE: partial differential equation.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

LR and MD conceived and designed the study, performed the analysis and drafted the manuscript. MH participated in the design of the study and helped to draft the manuscript. NK, SA and JR provided analytical tools. All authors read and approved the final manuscript.

Acknowledgments

We thank S Jun and E Stewart for sharing their data, D Chatenay, E Stewart, J Hoffmann, M Elez and G Batt for critical reading of the manuscript and Richard James for English editing. The research of M Doumic was supported by the ERC Starting Grant SKIPPER AD.

References

  1. Turner JJ, Ewald JC, Skotheim J M: Cell size control in yeast.

    Curr Biol 2012, 22:350-359. OpenURL

  2. Fantes PA: Control of cell size and cycle time inSchizosaccharomyces pombe.

    J Cell Sci 1977, 24:51-67. PubMed Abstract | Publisher Full Text OpenURL

  3. Sveiczer A, Novak B, Mitchison JM: The size control of fission yeast revisited.

    J Cell Sci 1996, 109:2947-2957. PubMed Abstract | Publisher Full Text OpenURL

  4. Johnston GC, Pringle JR, Hartwell LH: Coordination of growth with cell division in the yeastSaccharomyces cerevisiae.

    Exp Cell Res 1977, 105:79-98. PubMed Abstract | Publisher Full Text OpenURL

  5. Di Talia S, Skotheim JM, Bean JM, Siggia ED, Cross F R: The effects of molecular noise and size control on variability in the budding yeast cell cycle.

    Nature 2007, 448:947-951. PubMed Abstract | Publisher Full Text OpenURL

  6. Wang P, Hayden S, Masui Y: Transition of the blastomere cell cycle from cell size-independent to size-dependent control at the midblastula stage inXenopus laevis.

    J Exp Zool 2000, 287:128-144. PubMed Abstract | Publisher Full Text OpenURL

  7. Donachie WD: Relationship between cell size and time of initiation of DNA replication.

    Nature 1968, 219:1077-1079. PubMed Abstract | Publisher Full Text OpenURL

  8. Schaechter M, MaalOe O, Kjeldgaard NO: Dependency on medium and temperature of cell size and chemical composition during balanced growth ofSalmonella Typhimurium.

    Microbiology 1958, 19:592-606. OpenURL

  9. Cooper S, Helmstetter C E: Chromosome replication and the division cycle ofEscherichia coli.

    J Mol Biol 1968, 31:519-540. PubMed Abstract | Publisher Full Text OpenURL

  10. Donachie WD, Blakely G W: Coupling the initiation of chromosome replication to cell size inEscherichia coli.

    Curr Opin Microbiol 2003, 6:146-150. PubMed Abstract | Publisher Full Text OpenURL

  11. Boye E, Nordström K: Coupling the cell cycle to cell growth.

    EMBO Rep 2003, 4:757-760. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Wold S, Skarstad K, Steen HB, Stokke T, Boye E: The initiation mass for DNA replication inEscherichia coli k-12 is dependent on growth rate.

    Eur Mol Biol Organ J 1994, 13:2097-2102. OpenURL

  13. Bates D, Kleckner N: Chromosome and replisome dynamics inE. coli: loss of sister cohesion triggers global chromosome movement and mediates chromosome segregation.

    Cell 2005, 121:899-911. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Bates D, Epstein J, Boye E, Fahrner K, Berg H, Kleckner N: TheEscherichia coli baby cell column: a novel cell synchronization method provides new insight into the bacterial cell cycle.

    Mol Microbiol 2005, 57:380-391. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Hill NS, Kadoya R, Chattoraj DK, Levin PA: Cell size and the initiation of DNA replication in bacteria.

    PLoS Genet 2012, 8:1002549. Publisher Full Text OpenURL

  16. Boye E, Stokke T, Kleckner N, Skarstad K: Coordinating DNA replication initiation with cell growth: differential roles for DnaA and SeqA proteins.

    Proc Natl Acad Sci USA 12206., 93 OpenURL

  17. Chien A-C, Hill NS, Levin PA: Cell size control in bacteria.

    Curr Biol 2012, 22:340-349. Publisher Full Text OpenURL

  18. Teather RM, Collins JF, Donachie WD: Quantal behavior of a diffusible factor which initiates septum formation at potential division sites inEscherichia coli.

    J Bacteriol 1974, 118:407-413. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Bi E, Lutkenhaus J: Ftsz regulates frequency of cell division inEscherichia coli.

    J Bacteriol 1990, 172:2765-2768. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Margolin W: Ftsz and the division of prokaryotic cells and organelles.

    Nat Rev Mol Cell Biol 2005, 6:862-871. PubMed Abstract | Publisher Full Text OpenURL

  21. Tyson JJ: The coordination of cell growth and division – intentional or incidental?

    Bio Essays 1985, 2:72-77. OpenURL

  22. Trucco E, Bell GI: A note on the dispersionless growth law for single cells.

    Bull Math Biophys 1970, 32:475. PubMed Abstract | Publisher Full Text OpenURL

  23. Diekmann O, Heijmans HJAM, Thieme HR: On the stability of the cell size distribution.

    J Math Biol 1984, 113:227-248. OpenURL

  24. Metz JAJ, Diekmann O (Eds): The Dynamics of Physiologically Structured Populations. Lecture Notes in Biomathematics. Berlin: Springer; 1986:511.

    Papers from the colloquium held in Amsterdam, 1983

    OpenURL

  25. Stewart EJ, Madden R, Paul G, Taddei F: Aging and death in an organism that reproduces by morphologically symmetric division.

    PLoS Biol 2005, 3:45. Publisher Full Text OpenURL

  26. Wang P, Robert L, Pelletier J, Dang WL, Taddei F, Wright A, Jun S: Robust growth ofEscherichia coli.

    Curr Biol 2010, 20:1099-1103. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Doumic M, Hoffmann M, Reynaud-Bouret P, Rivoirard V: Nonparametric estimation of the division rate of a size-structured population.

    SIAM J Numer Anal 2012, 50:925-950. Publisher Full Text OpenURL

  28. Doumic M, Hoffmann M, Krell N, Robert L: Statistical estimation of a growth-fragmentation model observed on a genealogical tree.

    ArXiv 2012.

  29. Powell EO: Growth rate and generation time of bacteria, with special reference to continuous culture.

    Microbiology 1956, 15:492-511. OpenURL

  30. Kubitschek HE: Growth during the bacterial cell cycle: analysis of cell size distribution.

    Biophys J 1969, 9:792-809. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Collins JF, Richmond MH: Rate of growth ofBacillus cereus between divisions.

    J Gen Microbiol 1962, 28:15-33. PubMed Abstract | Publisher Full Text OpenURL

  32. Perthame B: Transport Equations in Biology. Basel: Birkhäuser Verlag; 2007:198.

    [Frontiers in Mathematics]

    OpenURL

  33. Wheals AE: Size control models ofSaccharomyces cerevisiae cell proliferation.

    Mol Cell Biol 1982, 2:361-368. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Schaechter M, Williamson JP, Hood Jun JR, Koch AL: Growth, cell and nuclear divisions in some bacteria.

    Microbiology 1962, 29:421-434. OpenURL

  35. Marr AG, Harvey RJ, Trentini WC: Growth and division ofEscherichia coli.

    J Bacteriol 1966, 91:2388-2389. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Cooper S: Leucine uptake and protein synthesis are exponential during the division cycle ofEscherichia coli b/r.

    J Bacteriol 1988, 170:436-438. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Yu XC, Margolin W: Ftsz ring clusters inmin and partition mutants: role of both the Min system and the nucleoid in regulating Ftsz ring localization.

    Mol Microbiol 1999, 32:315-326. PubMed Abstract | Publisher Full Text OpenURL

  38. Sun Q, Margolin W: Influence of the nucleoid on placement of FtsZ and MinE rings inEscherichia coli.

    J Bacteriol 2001, 183:1413-1422. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Migocki MD, Freeman MK, Wake RG, Harry EJ: The Min system is not required for precise placement of the midcell Z ring inBacillus subtilis.

    EMBO Rep 2002, 3:1163-1167. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Lutkenhaus J: Assembly dynamics of the bacterial MinCDE system and spatial regulation of the Z ring.

    Annu Rev Biochem 2007, 76:539-562. PubMed Abstract | Publisher Full Text OpenURL

  41. Wu LJ, Errington J: Nucleoid occlusion and bacterial cell division.

    Nat Rev Microbiol 2011, 10:8-12. PubMed Abstract | Publisher Full Text OpenURL

  42. Mulder E, Woldringh CL: Actively replicating nucleoids influence positioning of division sites inEscherichia coli filaments forming cells lacking DNA.

    J Bacteriol 1989, 171:4303-4314. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Doumic M, Perthame B, Zubelli JP: Numerical solution of an inverse problem in size-structured population dynamics.

    Inverse Probl 2009, 25:045008. Publisher Full Text OpenURL