Distinguishing different modes of growth using singlecell data
Abstract
Collection of highthroughput data has become prevalent in biology. Large datasets allow the use of statistical constructs such as binning and linear regression to quantify relationships between variables and hypothesize underlying biological mechanisms based on it. We discuss several such examples in relation to singlecell data and cellular growth. In particular, we show instances where what appears to be ordinary use of these statistical methods leads to incorrect conclusions such as growth being nonexponential as opposed to exponential and vice versa. We propose that the data analysis and its interpretation should be done in the context of a generative model, if possible. In this way, the statistical methods can be validated either analytically or against synthetic data generated via the use of the model, leading to a consistent method for inferring biological mechanisms from data. On applying the validated methods of data analysis to infer cellular growth on our experimental data, we find the growth of length in E. coli to be nonexponential. Our analysis shows that in the later stages of the cell cycle the growth rate is faster than exponential.
Editor's evaluation
In this manuscript, the authors describe a generative modelbased framework to better analyze stochastic growth data, including bacterial cell growth. They show how this framework can be applied to gain insight into the processes underlying these phenomena. This work is wellsupported by simulations and data analysis and will likely be of interest to those trying to understand the processes governing bacterial growth, as well as those studying stochastic growth processes in biology more broadly.
https://doi.org/10.7554/eLife.72565.sa0eLife digest
All cells – from bacteria to humans – tightly control their size as they grow and divide. Cells can also change the speed at which they grow, and the pattern of how fast a cell grows with time is called ‘mode of growth’. Mode of growth can be ‘linear’, when cells increase their size at a constant rate, or ‘exponential’, when cells increase their size at a rate proportional to their current size. A cell’s mode of growth influences its inner workings, so identifying how a cell grows can reveal information about how a cell will behave.
Scientists can measure the size of cells as they age and identify their mode of growth using single cell imaging techniques. Unfortunately, the statistical methods available to analyze the large amounts of data generated in these experiments can lead to incorrect conclusions. Specifically, Kar et al. found that scientists had been using specific types of plots to analyze growth data that were prone to these errors, and may lead to misinterpreting exponential growth as linear and vice versa.
This discrepancy can be resolved by ensuring that the plots used to determine the mode of growth are adequate for this analysis. But how can the adequacy of a plot be tested? One way to do this is to generate synthetic data from a known model, which can have a specific and known mode of growth, and using this data to test the different plots. Kar et al. developed such a ‘generative model’ to produce synthetic data similar to the experimental data, and used these data to determine which plots are best suited to determine growth mode. Once they had validated the best statistical methods for studying mode of growth, Kar et al. applied these methods to growth data from the bacterium Escherichia coli. This showed that these cells have a form of growth called ‘superexponential growth’.
These findings identify a strategy to validate statistical methods used to analyze cell growth data. Furthermore, this strategy – the use of generative models to produce synthetic data to test the accuracy of statistical methods – could be used in other areas of biology to validate statistical approaches.
Introduction
The last decade has seen a tremendous increase in the availability of highquality large datasets in biology, in particular in the context of singlecell level measurements. Such data are complementary to ‘bulk’ measurements made over a population of cells. They have led to new biological paradigms and motivated the development of quantitative models (Osella et al., 2017; Facchetti et al., 2017; Ho et al., 2018; Soifer et al., 2016; Jun et al., 2018; Amir and Balaban, 2018; Kohram et al., 2021). Nevertheless, they have also led to new challenges in data analysis, and here we will point out some of the pitfalls that exist in handling such data. In particular, we will show that the commonly used procedure of binning data and linear regression may hint at specific functional relations between the two variables plotted that are inconsistent with the true functional relations. As we shall show, this may come about due to the ‘hidden’ noise sources that affect the binning procedure and the phenomenon of ‘inspection bias’ where certain bins have biased contributions. One of our main take home messages is the significance of having an underlying model (or models) to guide/test/validate data analysis methods. The underlying model is referred to as a generative model in the sense that it leads to similar data to that observed in the experiments. The importance of a socalled generative model has been beautifully advocated in the context of astrophysical data analysis (Hogg et al., 2010), yet biology brings in a plethora of exciting differences: while in physics noise from measurement instruments often dominates, in the biological examples we will dwell on here it is the intrinsic biological noise that can obscure the mathematical relation between variables when not handled properly. In the following, we will illustrate this rather philosophical introduction on a concrete and fundamental example, albeit e pluribus unum. We will focus on the analysis of the Escherichia coli growth curves obtained via high throughput optical microscopy. Nevertheless we anticipate the conceptual points made here – and demonstrated on a particular example of interest – will translate to other types of measurements, which make use of microscopy but also beyond.
Binning corresponds to grouping data based on the value of the xaxis variable, and finding the mean of the fluctuating yaxis variable for this group. By removing the fluctuations of the yvariable, the binning process often aims to expose the ‘true’ functional relation between the two variables which can be used to infer the underlying biological mechanism. While binning may provide a smooth nonlinear relation between variables, linear regression is used to find a linear relationship between the variables. In addition to binning, we use the ordinary least squares regression where the slope and the intercept of the best linear fit line are obtained by minimizing the squared sum of the difference between the dependent variable raw data and the predicted value. Here, the best fit/the best linear fit is obtained using the raw data and not the binned data. Similar to binning, the assumption underlying linear regression is that our knowledge of xaxis variable is precise while the noise is in the yaxis variable.
It is important to discuss the sources of fluctuations in the yaxis variable before we proceed. In biology, fluctuations in the variables arise inevitably from the intrinsic variability within a cell population. Cells growing in the same medium and environment have different characteristics (e.g. growth rate) due to the stochastic nature of biochemical reactions in the cell (Kiviet et al., 2014). For example, the division event is controlled by stochastic reactions, whose variability leads to cell dividing at a size smaller or larger than the mean. In this paper, when modeling the data, we will consider the intrinsic noise as the only source of variability and assume that the measurement error is much smaller than the intrinsic variation in the population.
One example of the use of binning and linear regression is shown in Figure 1A where size at division (${L}_{d}$) vs size at birth (${L}_{b}$) is plotted using experimental data obtained by Tanouchi et al. for E. coli growing at 25 °C (Tanouchi et al., 2017). In Figure 1A, the functional relation between length at division and length at birth for E. coli is observed to be linear and close to ${L}_{d}={L}_{b}+\mathrm{\Delta}L$ (see the Experimental data section for details). The relation obtained allows us to hypothesize a coarsegrained biological model known as the adder model as shown in Figure 1B in which the length at division is set by addition of length $\mathrm{\Delta}L$ from birth (Soifer et al., 2016; Harris and Theriot, 2016; Si et al., 2019; Amir, 2014; Campos et al., 2014; TaheriAraghi et al., 2015; Eun et al., 2018). This previously discussed example demonstrates and reiterates the use of statistical analysis on singlecell data to understand the underlying cell regulation mechanisms. Using statistical methods such as binning and linear regression, other phenomenological models apart from adder have also been proposed in E. coli where the division length (${L}_{d}$) is not directly ‘set’ by that at birth (Ho and Amir, 2015; Micali et al., 2018; Witz et al., 2019). The phenomenological models, in turn, can be related to mechanistic (molecularlevel) models of cell size and cell cycle regulation (Barber et al., 2017). Recent work has shed light on the subtleties involved in interpreting the linear regression results for the ${L}_{d}$ vs ${L}_{b}$ plot where seemingly adder behavior in length can be obtained from a sizer model (division occurring on reaching a critical size) due to the interplay of multiple sources of variability (Facchetti et al., 2019). This issue is similar in spirit to those we highlight here.
The volume growth of single bacterial cells has been typically assumed to be exponential (Godin et al., 2010; Wang et al., 2010; Campos et al., 2014; Cermak et al., 2016; Soifer et al., 2016; IyerBiswas et al., 2014). Assuming ribosomes to be the limiting component in translation, growth is predicted to be exponential and growth rate depends on the active ribosome content in the cell (Scott et al., 2010; Lin and Amir, 2018; MetzlRaz et al., 2017). Under the assumption of exponential growth, the size at birth (${L}_{b}$), the size at division (${L}_{d}$), and the generation time (${T}_{d}$) are related to each other by,
where $\lambda $ is the growth rate. Understanding the mode of growth is important for example, due to its potential effects on cell size homeostasis. Exponentially growing cells cannot employ a mechanism where they control division by timing a constant duration from birth but such a mechanism is possible in case of linear growth (Amir, 2014; Kafri et al., 2016; Ho et al., 2018). Linear regression performed on $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ plot, where $\u27e8\lambda \u27e9$ is the mean growth rate, was used to infer the mode of growth in the archaeon H. salinarum (Eun et al., 2018), and in the bacteria M. smegmatis (Logsdon et al., 2017) and C. glutamicum (Messelink et al., 2020), for example. If the best linear fit follows the y = x trend, the resulting functional relation might point to growth being exponential. A corollary to this is the rejection of exponential growth when the slope and intercept of the best linear fit deviate from one and zero, respectively (Messelink et al., 2020). Thus, binning and linear regression applied on singlecell data appear to provide information about the underlying biology, in this case, the mode of cellular growth. We will test the validity of such inference by analyzing synthetic data generated using generative models. We find that linear regression performed on the plot $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$, surprisingly, does not provide information about the mode of growth. Nonetheless, we show that other methods of statistical analysis such as binning growth rate vs age plots are adequate in addressing the problem. Using these validated methods on experimental data, we find that E. coli grows nonexponentially. In later stages of the cell cycle, the growth rate is higher than that in early stages.
Results
Statistical methods like binning and linear regression should be interpreted based on a model
To illustrate the pitfalls associated with binning, we use data from recent experiments on E. coli where the length at birth, the length at division and the generation time were obtained for multiple cells (see Experimental methods and [TiruvadiKrishnan et al., 2021]). Phasecontrast microscopy was used to obtain cell length at equal intervals of time. Note that we consider length to reflect cell size in this paper rather than other cell geometry characteristics such as surface area and volume. The length growth rate that we elucidate in the paper can be different from the cell volume growth rate as shown in Appendix 1 assuming a simple cell morphology and exponential growth. Using the same cell morphology, we also find the length growth rate to be identical to cell surface growth rate. To investigate if the cell growth was exponential, we plotted $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ for cells growing in M9 alanine minimal medium at 28 °C ($\u27e8{T}_{d}\u27e9$ = 214 min). The linear regression of these data yields a slope of 0.3 and an intercept of 0.4 as shown in Figure 2A. The binned data and the best linear fit deviate significantly from the y = x line (see Supplementary file 1). Additionally, the binned data follows a nonlinear trend and flattens out at longer generation times. We also found similar deviations in the binned data and best linear fit in glycerol medium ($\u27e8{T}_{d}\u27e9$ = 164 min) shown in Figure 2—figure supplement 1A, and glucosecas medium ($\u27e8{T}_{d}\u27e9$ = 65 min) shown in Figure 2—figure supplement 1B. Qualitatively similar results have been recently obtained for another bacterium, C. glutamicum, in Messelink et al., 2020. These results might point to growth being nonexponential.
Next, we will approach the same problem but with a generative model. We will first show that the $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ binned plot could not distinguish exponential growth from nonexponential growth. For that purpose, we use a previously studied model (Eun et al., 2018) which considers growth to be exponential with the growth rate distributed normally and independently between cell cycles with mean growth rate $\u27e8\lambda \u27e9$ and standard deviation $C{V}_{\lambda}\u27e8\lambda \u27e9$. $C{V}_{\lambda}$ is thus the coefficient of variation (CV) of the growth rate and is assumed to be small. To maintain a narrow distribution of cell size, cells must employ regulatory mechanisms. In our model, we assume that, barring the noise due to stochastic biochemical reactions, cells attempt to divide at a particular size ${L}_{d}$ given size at birth ${L}_{b}$. Keeping the model as generic as possible, we can write ${L}_{d}$ as a function of ${L}_{b}$, f(${L}_{b}$) which can be thought of as a coarsegrained model for the regulatory mechanism. Amir, 2014 provides a framework to capture the regulatory mechanisms by choosing f(${L}_{b}$) = 2 ${L}_{b}^{1\alpha}{L}_{0}^{\alpha}$. L_{0} is the typical size at birth and $\alpha $, which can take values between 0 and 2, reflects the strength of regulation strategy. $\alpha $ = 0 corresponds to the timer model where division occurs on average after a constant time from birth, and $\alpha $ = 1 is the sizer model where a cell divides upon reaching a critical size. $\alpha $ = 1/2 can be shown to be equivalent to the adder model where division is controlled by addition of constant size from birth (Amir, 2014). In addition to the deterministic function (f) specifying division, the size at division is affected by noise ($\frac{\zeta}{\u27e8\lambda \u27e9}$) in division timing. We assume it has a Gaussian distribution with mean zero and standard deviation $\frac{{\sigma}_{n}}{\u27e8\lambda \u27e9}$ and that it is independent of the growth rate. Thus, the generation time (${T}_{d}$) can be mathematically written as ${T}_{d}=\frac{1}{\lambda}\mathrm{ln}(\frac{\mathrm{f}({L}_{b})}{{L}_{b}})+\frac{\zeta}{\u27e8\lambda \u27e9}$ and is influenced by growth rate noise and division timing noise. Note that replacing the time additive division timing noise with a size additive division timing noise will not affect the results qualitatively (see ‘Model’ and ‘Exponential growth’ sections for details and Supplementary file 1 for variable definitions).
For perfectly symmetrically dividing cells whose sizes are narrowly distributed, we find the trend in the binned data for $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ plot to be (see section ‘Predicting the results of statistical constructs applied on $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ and $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$’),
Fixing $C{V}_{\lambda}$ = ${\sigma}_{n}$ = 0.15, we show using simulations in Figure 2C the nonlinear trend in the binned data even though we assumed exponential growth. Similarly, on performing linear regression on the raw data of $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ plot, we find that the slope of the best linear fit is not equal to one and the intercept is nonzero (see Equation 27 and 28 and Figure 2C). Equation 2 shows that the trend in the binned data depends on the ratio of growth rate noise and division timing noise. The slope is equal to one and intercept is zero only if the noise in growth rate is negligible as compared to the division timing noise. In experiments that is rarely the case, hence, the binned data trend and the best linear fit deviate from the y = x line even though growth might be exponential. Thus, we cannot rule out exponential growth in the E. coli experiments despite the binned data trend being nonlinear and the bestfit line deviating from the y = x line.
Why does a nonlinear relationship in the binned data for the plot $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ arise even for exponential growth? According to the model, ${L}_{d}$ is determined by a deterministic strategy, f(${L}_{b}$) and a time/size additive division timing noise. The noise component which affects ${L}_{d}$ and subsequently the quantity $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ is thus the noise in division timing and not the growth rate. The generation time (${T}_{d}$) plotted on the xaxis is influenced by the noise in division timing as well as the noise in growth rate. Binning assumes that for a fixed value of the xaxis variable, the noise from other sources affects only the yaxis variable (the binned variable). Similarly for linear regression, the underlying assumption is that the independent variable on xaxis is precisely known while the dependent variable on the yaxis is influenced by the independent variable and from external factors other than the independent variable. In this case, only $\u27e8\lambda \u27e9{T}_{d}$ plotted on xaxis is influenced by growth rate noise while both $\u27e8\lambda \u27e9{T}_{d}$ and $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ are influenced by noise in division time. This does not fit the assumption for binning and linear regression and hence, the best linear fit for $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ plot might deviate from the y = x line even in the case of exponential growth.
Another way of explaining the deviation from the linear y = x trend is by inspection bias, which arises when certain data is overrepresented (Stein and Dattero, 2018). Cells which have a longer generation time than the mean will most likely have a slower growth rate. Thus, in Figure 2A and C, at larger values of $\u27e8\lambda \u27e9{T}_{d}$ or ${T}_{d}$, the bin averages are biased by slower growing cells, thus making $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ or $\lambda {T}_{d}$ to be lower than expected. This provides an explanation for the flattening of the trend.
It follows from the previous discussion that if one bins data by $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ then the assumption for binning is met. Both of the variables $\u27e8\lambda \u27e9{T}_{d}$ and $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ are influenced by the noise in division time but $\u27e8\lambda \u27e9{T}_{d}$ plotted on the yaxis is also influenced by the growth rate noise. Thus, the yaxis variable, $\u27e8\lambda \u27e9{T}_{d}$ is determined by the xaxis variable, $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$, and an external source of noise, in this case, the growth rate noise. Thus, based on our model, we expect the trend in binned data and linear regression performed on the interchanged axes to follow the y = x trend for exponentially growing cells (see section ‘Predicting the results of statistical constructs applied on $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ and $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$’). Indeed, on interchanging the axis and plotting $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ for synthetic data, we find that the trend in the binned data and the best linear fit closely follows the y = x line (Figure 2D). We also find that the best linear fit follows the y = x line in the case of alanine (Figure 2B), glycerol (Figure 2—figure supplement 1A) and glucosecas (Figure 2—figure supplement 1B). A change from nonlinear behavior to that of linear on interchanging the axes is also observed in a related problem where growth rate ($\lambda $) and inverse generation time ($\frac{1}{{T}_{d}}$) are considered (Figure 2—figure supplement 2 and Section ‘Interchanging axes in growth rate vs inverse generation time plot might lead to different interpretations’).
Thus far, we showed for a range of models where birth controls division that the binned data trend for $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ as function of $\u27e8\lambda \u27e9{T}_{d}$ is nonlinear and dependent on the noise ratio $\frac{{\sigma}_{n}}{C{V}_{\lambda}}$ in the case of exponential growth. On interchanging the axes the binned data trend agrees with the y = x line independent of the growth rate and division time noise. However, we will show next that this agreement with the y = x trend cannot be used as a ‘smoking gun’ for inferring exponential growth from the data.
To investigate this further, let us consider linear growth, which has also been suggested to be followed by E. coli cells (Mitchison, 2005; Abner et al., 2014). The underlying equation for linear growth is,
where ${\lambda}^{\prime}$ is the the elongation speed that is, $\frac{dL}{dt}$. For cells growing linearly, the best linear fit for the plot $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ is expected to deviate from the y = x line. As before, we fix $\u27e8\lambda \u27e9$ to be the mean of $\frac{1}{{T}_{d}}\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$, agnostic of the linear mode of growth. Surprisingly, we found that for the class of models where birth controls division by a strategy f(${L}_{b}$) and cells grow linearly, the best linear fit for $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ agrees closely with the y = x trend. On carrying out analytical calculations based on this model, we obtain the slope and the intercept of the $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ plot to be $\frac{3}{2}\mathrm{ln}(2)\approx $ 1.04 and –0.03 respectively, which is very close to that for exponential growth (see section ‘Differentiating linear from exponential growth’). This is shown for simulations of linear growth with cells following an adder model in Figure 3A. Given no information about the underlying model, Figure 3A could be interpreted as cells undergoing exponential growth contrary to the assumption of linear growth in simulations. Thus, when handling experimental data, cells undergoing either exponential or linear growth might seem to agree closely with the y = x trend. Deforet et al., 2015 used the linear binned data trend in case of $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ plot to infer exponential growth but as we showed in this section, the linear trend does not rule out linear growth. This again reiterates our message of having a generative model to guide the data analysis methods such as binning and linear regression. For completeness, we also test the utility of $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8{T}_{d}\u27e9\lambda $ and its interchanged axes plots to elucidate the mode of growth (Appendix 2). We find that binning and linear regression applied on these plots can not differentiate between exponential and linear growth.
To conclude the discussion of linear growth, we note that the natural plot for this growth regime is $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ vs ${l}_{d}{l}_{b}$ and the plot obtained on interchanging the axes (see the Linear growth section and Figure 3—figure supplement 1A and B). Here l_{b}, l_{d} and ${\lambda}_{lin}$ are defined to be quantities ${L}_{b}$, ${L}_{d}$ and ${\lambda}^{\prime}$, respectively, normalized by the mean length at birth. For cells growing exponentially, the best linear fit for the $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ vs ${l}_{d}{l}_{b}$ plot is expected to deviate from the y = x line. This is indeed what is observed in Figure 3—figure supplement 1 where simulations of exponentially growing cells following the adder model are presented (see ‘Differentiating linear from exponential growth’ for extended discussion).
In all the cases above, the problem at hand deals with distilling the biologically relevant functional relation between two variables. However, the data is assumed to be subjected to fluctuations of various sources, and it is important to ensure that the statistical construct we are using (e.g. binning) is robust to these. How can we know a priori whether the statistical method is appropriate and a ‘smoking gun’ for the functional relation we are conjecturing? The examples shown above suggest that performing statistical tests on synthetic data obtained using a generative model is a convenient and powerful approach. Note that in cases such as the ones studied here where analytical calculations may be performed, one may not even need to perform any numerical simulations to test the validity of the methods.
Growth rate vs age plots are consistent with the underlying growth mode
In the last section, we showed that the plots $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ and $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ are not decisive in identifying the mode of growth. Recent works on B. subtilis (Nordholt et al., 2020) and fission yeast (Knapp et al., 2019) have used differential methods of quantifying growth namely growth rate ( = $\frac{1}{L}\frac{dL}{dt}$) vs age plots and elongation speed (=$\frac{dL}{dt}$) vs age plots to probe the mode of growth within a cell cycle. Here, $L$ denotes the size of the cell after time $t$ from birth in the cell cycle and age denotes the ratio of time $t$ to ${T}_{d}$ within a cell cycle (hence it ranges from 0 to 1 by construction within a cell cycle). In this section, using various models of cell growth and cell cycle, we test the growth rate vs age method. Note that the growth rate vs age and the elongation speed vs age plots are not dimensionless unlike the previous plots. Using the growth rate vs age and elongation speed vs age plots, we aim to quantify the growth rate changes within a cell cycle. For cells assumed to be growing exponentially, growth rate is constant throughout the cell cycle. On averaging over multiple cell cycles, the trend of binned data is expected to be a horizontal line with value equal to mean growth rate which is indeed what we find in the numerical simulations of the adder and the adder per origin model (Ho and Amir, 2015), as shown in Figure 3B. The binned data trend in each of the models matches the theoretical predictions of growth rate (shown as dotted lines). In contrast, for linearly growing cells, the elongation speed is expected to remain constant. We show this constancy using numerical simulations of linearly growing cells following the adder model (Figure 3—figure supplement 3A). In accordance with this result, the growth rate is expected to decrease with cell age as $\lambda \propto \frac{1}{1+age}$. This is verified in Figure 3B by again using the numerical simulations of linear growth with cells following the adder model. The binned data trend for linear growth (green squares) matches the theoretical predictions of $\lambda \propto \frac{1}{1+age}$ (green dotted line).
Thus, the two growth modes (exponential and linear) could be differentiated using the growth rate vs age plot (for details see ‘Growth rate vs age and elongation speed vs age plots’ section). However, the growth rate vs age plots can be used to infer the mode of growth beyond the two discussed above. We show this by using simulations of cells following the adder model and undergoing faster than exponential or superexponential growth (see the Simulations section for details). In such a case, the growth rate is expected to increase. This increase in growth rate is shown in Figure 3B using simulations. The binned data trend (red triangles) again matches the growth rate mode used in the simulations (red dotted line). Thus, the growth rate vs age plots are a consistent method to distinguish linear from exponential and superexponential growths.
Using the validated growth rate vs age plots, we obtained the growth rate trend for experimental data on E. coli for the three growth conditions studied in this paper (Figure 4AC). We found an increase in growth rate in all growth conditions during the course of the cell cycle. One may wonder whether such an increase may be explained by the E. coli morphology alone, due to the presence of hemispherical poles. For exponentially growing cell volume and considering a geometry of E. coli with spherical caps at the poles, the percentage increase in the growth rate of length over a cell cycle is around 3 % which is significantly smaller than that observed in our experimental data. Considering cell size trajectories (cell size, $L$ at time, $t$ data) where cell lengths were tracked beyond the cell division event (by considering cell size in both daughter cells), we also found that the growth rate decreases close to division (age ≈ 1) and returns to a value nearly equal to that observed at the beginning of cell cycle (age ≈ 0) as shown in Figure 4—figure supplement 1 (see ‘Growth rate vs age and elongation speed vs age plots’ section for extended discussion).
The above question of mode of growth within a cell cycle can also be analyzed in relation to a specific event. Several studies have pointed to a change in growth rate at the onset of constriction (Reshes et al., 2008; Banerjee et al., 2017). This change in growth rate can be probed using growth rate vs time plots where time is taken relative to the onset of constriction as shown in Figure 4—figure supplement 2. These plots show a decrease in growth rates at the two extremes of the plot. These decreases are due to inspection bias, where the growth rate trend is affected by the biased contribution of cells with a higher than average generation time or equivalently slower growth rate (see ‘Growth rate vs time from specific event plots are affected by inspection bias’ section for extended discussion). Inspection bias is also observed when timing is considered relative to other cell events such as cell birth (see ‘Growth rate vs time from specific event plots are affected by inspection bias’ section and Figure 3—figure supplement 2).
It might not always be possible to obtain growth rate trajectories as a function of time/cell age. Godin et al. instead obtained the instantaneous biomass growth speed ($\frac{dM}{dt}$) as a function of its buoyant mass ($M$) (Godin et al., 2010). On applying linear regression for instantaneous mass growth speed vs mass, we expect the slope of the best linear fit obtained to provide the average growth rate ($\u27e8\frac{1}{M}\frac{dM}{dt}\u27e9$) under the assumption of exponential growth while for linear growth the intercept provides the average growth speed. Using this method, biomass was suggested to be growing exponentially. This method can be applied to study the length growth rate within the cell cycle by plotting elongation speed as a function of length (Cadart et al., 2019). We find that the binned data trend and the best linear fit of this plot follow the expected trend for linear and exponential growth as shown in Figure 3—figure supplement 3B and Figure 3—figure supplement 3D, respectively, for a cell cycle model where division is controlled via an adder mechanism from birth. However, the trend obtained appears to be modeldependent as shown in Figure 3—figure supplement 3F where the underlying cell cycle model used in the simulations is the adder per origin model. For this model, the binned data trend is found to be nonlinear with the growth rate speeding up at large sizes, despite the synthetic data being generated for perfectly exponential growth. This nonlinear trend can lead to growth rate being misinterpreted as nonexponential within the cell cycle (see ‘Results of elongation speed vs size plots are modeldependent’ section for details). Thus, an analysis using the elongation speed vs size plot must be accompanied with an underlying cell cycle model.
In summary, we found that the growth rate vs age plot was a consistent method to determine the changes in growth rate within a cell cycle. Unlike the growth rate vs age plots, the inference from the growth rate vs size plots was found to be modeldependent. Using the growth rate vs age plots, we show that the length growth of E. coli can be faster than exponential.
Discussion
Statistical methods such as binning and linear regression are useful for interpreting data and generating hypotheses for biological models. However, we show in this paper that predicting the relationships between experimentally measured quantities based on these methods might lead to misinterpretations. Constructing a generic model and verifying the statistical analysis on the synthetic data generated by this model provides a more rigorous way to mitigate these risks.
In the paper, we provide examples in which $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ and $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ plots fail as a method to infer the mode of growth. The binned data trend and the best linear fit for the $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ plot was found to be dependent upon the noise parameters in the class of models where birth controlled division (Equation 2). We also show that $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ plot could not differentiate between exponential and linear modes of growth (Figures 2D and 3A). Thus, we conclude that the best linear fit for the above plots might not be a suitable method to infer the mode of growth but they are just one of the many correlations which the correct cell cycle model should be able to predict.
We found growth rate vs age and elongation speed vs age plots to be consistent methods to probe growth within a cell cycle. The method was validated using simulations of various cell cycle models (such as the adder, and adder per origin model, where in the latter, control over division is coupled to DNA replication) and the binned growth rate trend agreed closely with the underlying mode of growth for the wide range of models considered (Figure 3B). In the case of growth rate vs time plots, it was important to take into consideration the effects of inspection bias. We used cell cycle models to show the time regimes where inspection bias could be observed (Figure 3—figure supplement 2). In the regime with negligible inspection bias, we could reconcile the growth rate trend obtained using growth rate vs age (Figure 4AC) and growth rate vs time plots (Figure 4—figure supplement 2). The authors in Messelink et al., 2020 circumvent inspection bias in the elongation speed vs time from birth plots by focusing their analysis on the time period from cell birth to the generation time of the fastest dividing cell. The authors of Panlilio et al., 2021, while investigating the division behavior in the cells undergoing nutrient shift within their cell cycle, use both models and experimental data from steadystate conditions to identify inspection bias. These serve as good examples of using models to aid data analysis.
Statistics obtained from linear regression such as in Figure 1A help narrow down the landscape of cell cycle models, but many have potential pitfalls lurking which might lead to misinterpretations (Figures 2C and 3A). There are additional issues beyond those concerning linear regression and binning discussed here. For example, Willis and Huang, 2017 discusses Simpson’s paradox (Simpson, 1951) where distinct cellular subpopulations might lead to erroneous interpretation of cell cycle mechanisms. Examples of such distinct subpopulations are found in asymmetrically dividing bacteria such as M. smegmatis (Aldridge et al., 2012; Logsdon et al., 2017). Another source of misinterpretation could arise from presence of measurement errors. Throughout this work, we deal with intrinsic noise and neglect measurement error. However, when measurement noise affects both xaxis and yaxis variables, the slope of the best linear fit is biased towards zero. This can lead to potentially related variables being misinterpreted as uncorrelated. Measurement errors can, however, be handled based on a model. Using a model which includes measurement error as a source of noise, we can guide the binning analysis. Using this methodology, we verified that typical measurement errors ($\approx 0.02{L}_{b}$) Messelink et al., 2020; Kaiser et al., 2018 have negligible effects on the growth rate trends obtained from the experimental data used in our work.
Single cell size in E. coli has been reported to grow exponentially (Campos et al., 2014; Wang et al., 2010; Cermak et al., 2016; Soifer et al., 2016; IyerBiswas et al., 2014; Godin et al., 2010), linearly (Mitchison, 2005), bilinearly (Kubitschek, 1981) or trilinearly (Reshes et al., 2008). These are inconsistent with our observations in Figure 4AC where we find that growth can be superexponential. The nonmonotonic behavior in the fastestgrowth condition is reminiscent of the results reported in Nordholt et al., 2020 for B. subtilis. The authors of Nordholt et al., 2020 attribute the increase in growth rate to a multitude of cell cycle processes such as initiation of DNA replication, divisome assembly, septum formation. In the two slower growth conditions (Figure 4AB), we find that the growth rate increase starts before the time when the septal cell wall synthesis starts i.e., the constriction event. However, in the fastest growth condition (Figure 4C), the timing of growth rate increase seems to coincide with the onset of constriction which is in agreement with previous findings (Reshes et al., 2008; Banerjee et al., 2017).
It is important to distinguish between length growth and biomass growth. Oldewurtel et al., 2021 measures biomass and cell volume and finds the massdensity variations within the cellcycle to be small. In this paper, since we observe the length growth to be nonexponential (Figure 4), it remains to be seen whether biomass growth also follows a similar nonexponential behavior or if it is exponential as previously suggested (Godin et al., 2010; Oldewurtel et al., 2021).
In conclusion, the paper draws the attention of the readers to the careful use of statistical methods such as linear regression and binning. Although shown in relation to cell growth, this approach to data analysis seems ubiquitous. The general framework of carrying out data analysis is presented in Figure 5. It proposes the construction of a generative model based on the experimental data collected. Of course, we do not always know whether the model used is an adequate description of the system. What is the fate of the methodology described here in such cases? First, we should be reminded of Box’s famous quote ‘all models are wrong, some are useful’. The goal of a model is not to provide as accurate a description of a system as possible, but rather to capture the essence of the phenomena we are interested in and stimulate further ideas and understanding. In our context, the goal of the model is to provide a rigorous framework in which data analysis tools can be critically tested. If verified within the model, it is by no means proof of the success of the model and the method itself, and further comparisons with the data may falsify it leading to the usual (and productive) cycle of model rejection and improvement via comparison with experiments. However, if the best model we have at hand shows that the data analysis method is noninformative, as we have shown here on several methods used to identify the mode of growth, then clearly we should revise the analysis as it provides us with a nonconsistent framework, where our modeling is at odds with our data analysis. Furthermore, testing the methods on a simplified model is still advantageous compared with the option of using the methods without any validation. To mitigate the risk of using irrelevant models, in some cases it may be desirable to test the analysis methods on as broad a class of models as possible as we have done in the paper, for example by our use of a general value of $\alpha $ to describe the sizecontrol strategy within our models. Thus, guided by the model, the data analysis methods can be ultimately applied to experimental data and underlying functional relationships can be inferred. Reiterating the message of the authors in Hogg et al., 2010, the data analysis using this framework aims to justify the methods being used, thus, reducing arbitrariness and promoting consensus among the scientists working in the field.
Materials and methods
Experimental methods
Request a detailed protocolStrain engineering: STK13 strain (ΔftsN::frtYpetFtsN, ΔdnaN::frtmCherrydnaN) is derivative of E. coli K12 BW27783 (CGSC#: 12119) constructed by $\lambda $Red engineering (Datsenko and Wanner, 2000) and by P1 transduction (Thom, 2007). For chromosomal replacement of ftsN with fluorescence derivative, we used primers carrying 40nt tails with identical sequence to the ftsN chromosomal locus and a plasmid carrying a copy of ypet preceded by a kanamycin resistance cassette flanked by frt sites (frtkan^{R}frtYpetlinker) as PCR template (a kind gift from R. ReyesLamothe McGill University, Canada; ReyesLamothe et al., 2010). The resulting PCR product was transformed by electroporation into a strain carrying the $\lambda $Redexpressing plasmid pKD46. Colonies were selected by kanamycin resistance, verified by fluorescence microscopy and by PCR using primers annealing to regions flanking ftsN gene. After removal of kanamycin resistance by expressing the Flp recombinase from plasmid pCP20 (Cherepanov and Wackernagel, 1995), we transferred the mCherrydnaN gene fusion (BN1682 strain; a kind gift from Nynke Dekker from TUDelft, The Netherlands, Moolman et al., 2014) into the strain by P1 transduction. To minimize the effect of the insertion on the expression levels of the gene we removed the kanamycin cassette using Flp recombinase expressing plasmid pCP20.
Cells growth, preparation, and culturing E. coli in mother machine microfluidic devices: All cells were grown and imaged in M9 minimal medium (Teknova) supplemented with 2 mM magnesium sulfate (Sigma) and corresponding carbon sources at 28 °C. Three different carbon sources were used: 0.5 % glucose supplemented by 0.2 % casamino acids (Cas) (Sigma), 0.3 % glycerol (Fisher), and 0.3 % alanine (Fisher) supplemented with 1 x trace elements (Teknova).
For microscopy, we used mother machine microfluidic devices made of PDMS (polydimethylsiloxane). These were fabricated following to previously described procedure (Yang et al., 2018). To grow and image cells in microfluidic device, we pipetted 2–3 µl of resuspended concentrated overnight culture of OD_{600}∼ 0.1 into main flow channel of the device and let cells to populate the deadend channels. Once these channels were sufficiently populated (about 1 hr), tubing was connected to the device, and the flow of fresh M9 medium with BSA (0.75 µg/ml) was started. The flow was maintained at 5 µl/min during the entire experiment by an NE1000 Syringe Pump (New Era Pump Systems, NY). To ensure steadystate growth, the cells were left to grow in channels for at least 14 hr before imaging started.
Microscopy: A Nikon TiE inverted epifluorescence microscope (Nikon Instruments, Japan) with a 100 X (NA = 1.45) oil immersion phase contrast objective (Nikon Instruments, Japan), was used for imaging the bacteria. Images were captured on an iXon DU897 EMCCD camera (Andor Technology, Ireland) and recorded using NISElements software (Nikon Instruments, Japan). Fluorophores were excited by a 200 W Hg lamp through an ND8 neutral density filter. A Chroma 41,004 filtercube was used for capturing mCherry images, and a Chroma 41,001 (Chroma Technology Corp., VT) for Ypet images. A motorized stage and a perfect focus system were utilized throughout timelapse imaging. Images in all growth conditions were obtained at 4 min frame rate.
Image analysis: Image analysis was carried out using Matlab (MathWorks, MA) scripts based on Matlab Image Analysis Toolbox, Optimization Toolbox, and DipImage Toolbox (https://www.diplib.org/). Cell lengths were determined based on segmented phase contrast images. Dissociation of YpetFtsN label from cell middle was used to determine the exact timing of cell divisions.
Further experimental details can also be found in TiruvadiKrishnan et al., 2021.
Model
Request a detailed protocolConsider a model of cell cycle characterized by two events: cell birth and division. In our model, we assume that, barring the noise, cells tend to divide at a particular size v_{d} given size at birth v_{b}, via some regulatory mechanism. Hence, we can write v_{d} as a function of v_{b}, f(v_{b}). Amir, 2014 provides a framework to capture the regulatory mechanisms by choosing f(v_{b}) = 2 ${v}_{b}^{1\alpha}{v}_{0}^{\alpha}$. v_{0} is the typical size at birth and $\alpha $ captures the strength of regulation strategy. $\alpha $ = 0 corresponds to the timer model where division occurs after a constant time from birth, and $\alpha $ = 1 is the sizer where a cell divides on reaching a critical size. $\alpha $ = 1/2 can be shown to be equivalent to an adder where division is controlled by addition of constant size from birth (Amir, 2014). From here on, we would be using the length of the cell (${L}_{b}$, ${L}_{d}$, etc.) as a proxy for size (v_{b}, v_{d}, etc.). To reiterate, the length growth is not the same as cell volume growth as shown in Appendix 1. All of the variable definitions are summarized in Supplementary file 1. We also define ${l}_{b}=\frac{{L}_{b}}{\u27e8{L}_{b}\u27e9}$ and ${l}_{d}=\frac{{L}_{d}}{\u27e8{L}_{b}\u27e9}$. Using this, we can write the division strategy f(l_{b}) to be l_{d} = f(l_{b}) = 2 ${l}_{b}^{1\alpha}$. The total division size obtained will be a combination of f(l_{b}) and noise in the division timing, the source of which could be the stochasticity in biochemical reactions controlling division.
We will assume that division is perfectly symmetric i.e., size at birth in the ${(n+1)}^{th}$ generation (${l}_{b}^{n+1}$) is half of size at division in the ${n}^{th}$ generation (${l}_{d}^{n}$). Using the size additive division timing noise (${\zeta}_{s}(0,{\sigma}_{bd})$) and f(l_{b}) specified above, we obtain,
where x_{n} = $\mathrm{ln}({l}_{b}^{n})$. Size at birth (${L}_{b}$) is narrowly distributed, hence ${l}_{b}\approx 1$ and we can write $x$ = $\mathrm{ln}({l}_{b})$ = $\mathrm{ln}(1+\delta )$ where $\delta $ is a small number. We obtain $x\ll 1$ and,
The size additive noise, ${\zeta}_{s}(0,{\sigma}_{bd})$ is assumed to be small and has a normal distribution with mean 0 and standard deviation ${\sigma}_{bd}$. Note that ${\sigma}_{bd}$ is a dimensionless quantity. Since ${\zeta}_{s}(0,{\sigma}_{bd})$ is assumed to be small and ${x}_{n}\ll 1$, we can Taylor expand the last term of Equation 4 to first order,
Equation 6 shows a recursive relation for cell size and it is agnostic of the mode of growth. We will show later for exponential growth that replacing the size additive noise with time additive noise does not change the structure of Equation 6.
Exponential growth
Next, we will try to obtain the generation time (${T}_{d}$) in the case of exponentially growing cells. For exponential growth, the time at division ${T}_{d}$ is given by,
For simplicity, we will assume a constant growth rate ($\lambda $) within the cellcycle. Growth rate is fixed at the start of the cellcycle and is given by $\lambda $ = $\u27e8\lambda \u27e9$ + $\u27e8\lambda \u27e9\xi (0,C{V}_{\lambda})$, where $\u27e8\lambda \u27e9$ is the mean growth rate and $\xi (0,C{V}_{\lambda})$ is assumed to be small with a normal distribution that has mean 0 and standard deviation $C{V}_{\lambda}$ . $C{V}_{\lambda}$ denotes the coefficient of variation (CV) of the growth rate. This captures the variability in growth rate within cells arising from the stochastic nature of biochemical reactions occurring within the cell.
Size additive noise
Request a detailed protocolHere we will calculate the generation time using the division strategy f(l_{b}) and a size additive division timing noise (${\zeta}_{s}(0,{\sigma}_{bd})$) as described previously. On substituting ${L}_{d}=(f({l}_{b})+{\zeta}_{s})\u27e8{L}_{b}\u27e9$ into Equation 7 we obtain,
where the size additive noise (${\zeta}_{s}(0,{\sigma}_{bd})$) is Gaussian with mean 0 and standard deviation ${\sigma}_{bd}$.
The noise ${\zeta}_{s}(0,{\sigma}_{bd})$ is assumed to be small, and we obtain to first order,
Since ${x}_{n}\ll $ 0, on Taylor expanding $\frac{1}{{(1+{x}_{n})}^{1\alpha}}$ to first order,
Assuming noise in growth rate to be small and expanding to first order, we obtain,
Equation 11 gives the generation time for the class of models where birth controls division under the assumption that growth is exponential.
Time additive noise
Request a detailed protocolNext, we ensure that the recursive relation for size at birth and the expression for the generation time given by Equations 6 and 11, respectively, are robust to the nature of noise assumed. In this section, the generation time is obtained using the division strategy f(l_{b}) as described previously along with a time additive division timing noise ($\frac{\zeta}{\u27e8\lambda \u27e9}$). In such a case, ${T}_{d}$ is obtained to be,
The time additive noise, $\frac{\zeta (0,{\sigma}_{n})}{\u27e8\lambda \u27e9}$, is assumed to be small and has a normal distribution with mean 0 and standard deviation $\frac{{\sigma}_{n}}{\u27e8\lambda \u27e9}$. Note that ${\sigma}_{n}$ is a dimensionless quantity.
Assuming noise in growth rate to be small, we find ${T}_{d}$ to first order to be,
Equation 13 is same as Equation 11, if the time additive noise term, $\zeta (0,{\sigma}_{n})$, in Equation 12 is replaced by ${\zeta}_{s}(0,{\sigma}_{bd})/2$. Using Equation 13, the variance in ${T}_{d}$ (${\sigma}_{t}^{2}$) is,
For exponential growth, we also find,
On substituting Equation 12 into Equation 15 we obtain to first order,
On replacing the time additive noise term, $\zeta (0,{\sigma}_{n})$, in Equation 16 with ${\zeta}_{s}(0,{\sigma}_{bd})/2$, we recover the recursive relation for size at birth obtained in the case of size additive noise shown in Equation 6. Hence, the model is insensitive to noise being size additive or time additive with a simple mapping for going from one noise type to another in the small noise limit.
At steady state, $x$ has a normal distribution with mean 0 and variance ${\sigma}_{x}^{2}$ whose value is given by,
We note that some of the derivations above have also been presented in Eun et al., 2018, but are provided here for completeness.
Predicting the results of statistical constructs applied on $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ and $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$
Request a detailed protocolObtaining the best linear fit
Request a detailed protocolNext, we calculate the equation for the best linear fit for the choice of $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ as yaxis and $\u27e8\lambda \u27e9{T}_{d}$ as xaxis and vice versa. For simplicity, in this section, we will consider time additive division timing noise. However, the results obtained here will hold for size additive noise as well because the model is robust to the type of noise added as shown in the previous section.
First, we calculate the correlation coefficient (${\rho}_{exp}$) for $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ and time of division ${T}_{d}$,
where ${\sigma}_{l}$ is the standard deviation in $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$. Using Equations 15 and 16 we obtain,
Substituting Equations 13 and 19 into the numerator of Equation 18,
As the terms $\zeta (0,{\sigma}_{n})$, $\xi (0,C{V}_{\lambda})$ and x_{n} are independent of each other, $\u27e8\xi (0,C{V}_{\lambda})\zeta (0,{\sigma}_{n})\u27e9$ = 0, $\u27e8\xi (0,C{V}_{\lambda}){x}_{n}\u27e9$ = 0 and $\u27e8{x}_{n}\zeta (0,{\sigma}_{n})\u27e9$ = 0. Equation 20 simplifies to,
The variance of $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ obtained using Equation 19 is,
Inserting Equations 14, 21 and 22 into Equation 18, we get,
The slope of a linear regression line is given by,
where ${\sigma}_{x}$, ${\sigma}_{y}$, and $\rho $ are the standard deviation of the xvariable, the standard deviation of the yvariable and the correlation coefficient of the (x,y) pair, respectively. The intercept is,
On the xaxis, we plot $\u27e8\lambda \u27e9{T}_{d}$ and the yaxis is chosen as $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$. The slope for this choice (${m}_{tl}$) can be calculated by,
On substituting the values we get,
Only for $C{V}_{\lambda}\ll {\sigma}_{n}$ we would expect a slope close to 1.
The intercept ($c}_{tl$) for the $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ plot is given by,
However, if we choose the xaxis as $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ and the yaxis is chosen as $\u27e8\lambda \u27e9{T}_{d}$, we obtain the slope ${m}_{lt}$,
On substituting the values we obtain ${m}_{lt}$ = 1 independent of the noise parameters and find that the intercept is zero.
Nonlinearity in binned data
Request a detailed protocolIn the Main text, for the plot $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$, we find the binned data to be nonlinear (see Figure 2C of the Main text). In this section, we explain the nonlinearity observed using the model developed in the previous sections.
Binning data based on the xaxis means taking an average of the yvariable conditioned on the value of the xvariable. Mathematically, this amounts to calculating $\mathbb{E}[y\mid x]$ i.e., the conditional expectation of the yvariable given that x is fixed. In our case, we need to calculate $\mathbb{E}[\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})\mid \u27e8\lambda \u27e9{T}_{d}]$ . $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$= $\lambda {T}_{d}$ by definition of exponential growth, hence,
Since ${T}_{d}$ is fixed, this is equivalent to calculating $\mathbb{E}[\lambda \mid {T}_{d}]$. Using Equation 13,
$p(x,\xi ,\zeta )$ is the joint probability distribution of $x$ and noise parameters $\xi $ and $\zeta $. Since, they are independent of each other, the joint distribution is product of the individual distributions ${f}_{1}(x)$, ${f}_{2}(\xi )$ and ${f}_{3}(\zeta )$, the distributions being Gaussian with mean 0 and standard deviation ${\sigma}_{x}$, $C{V}_{\lambda}$ and ${\sigma}_{n}$, respectively. ${\sigma}_{x}$, ${\sigma}_{n}$ are related by Equation 17. Since $x$, $\xi $, and $\zeta $ are narrowly distributed around zero, the contribution from large positive or negative values is extremely small. This ensures that ${T}_{d}$ is also close to its mean and nonnegative despite the limits of the integral being $\mathrm{\infty}$ to $\mathrm{\infty}$. Using $\lambda $ = $\u27e8\lambda \u27e9+\u27e8\lambda \u27e9\xi (0,C{V}_{\lambda})$ in Equation 31,
On evaluating the integrals, we obtain,
Thus, the trend of binned data is found to be,
In the regime $C{V}_{\lambda}\ll {\sigma}_{n}$ , the last two terms on the RHS of Equation 34 vanish and the binned data follows the trend y = x.
For the $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ plot, we need to calculate $\mathbb{E}[\u27e8\lambda \u27e9{T}_{d}\mid \mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})]$. Using Equations 13 and 19, we obtain,
$\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ is independent of $\xi (0,C{V}_{\lambda})$. Using this, we can write $\mathbb{E}[\u27e8\lambda \u27e9{T}_{d}\mid \mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})]$ as,
Note that the integral over $\u27e8\lambda \u27e9{T}_{d}$ goes from $\mathrm{\infty}$ to $\mathrm{\infty}$ although $\u27e8\lambda \u27e9{T}_{d}$ cannot be negative. As before, this is not an issue because we assume $\u27e8\lambda \u27e9{T}_{d}$ to be tightly regulated around $\mathrm{ln}(2)$ and the contribution to the integral from $\mathrm{\infty}$ to 0 is negligible. ${f}_{4}(\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}}))$ denotes the probability distribution for $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$, the distribution being Gaussian with mean $\mathrm{ln}(2)$, and standard deviation ${\sigma}_{l}$ which is calculated in Equation 22. Putting the Gaussian form of ${f}_{2}(\xi )$ into the integral and simplifying we get,
The trend of binned data to first order in noise and $x$ is $\mathbb{E}[\u27e8\lambda \u27e9{T}_{d}\mid \mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})]$ = $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$. This is shown in Figure 2D of the Main text where the binned data follows the y = x line.
Linear growth
Request a detailed protocolIn this section, we will focus on finding the equation of the best linear fit for relevant plots in the case of linear growth. The time at division for linear growth is given by,
Note that ${\lambda}^{\prime}$ has units of [length/time] and is defined as the elongation speed. This is different from the exponential growth rate which has units [1/time]. Here, we will work with the normalized length at birth (l_{b}) and division (l_{d}),
Consider the normalized elongation speed to be ${\lambda}_{lin}$ = $\u27e8{\lambda}_{lin}\u27e9+\u27e8{\lambda}_{lin}\u27e9{\xi}_{lin}(0,C{V}_{\lambda ,lin})$, where $\u27e8{\lambda}_{lin}\u27e9$ is the mean normalized elongation speed for a lineage of cells and ${\xi}_{lin}(0,C{V}_{\lambda ,lin})$ is normally distributed with mean 0 and standard deviation $C{V}_{\lambda ,lin}$. Thus, the CV of elongation speed is $C{V}_{\lambda ,lin}$. The regulation strategy which the cell undertakes is equivalent to that in previous sections and is given by g(l_{b}) = $2+2(1\alpha )({l}_{b}1)$. Note that we can obtain g(l_{b}) by Taylor expanding f(l_{b}) around l_{b} = 1. Using the regulation strategy g(l_{b}) and adding a size additive noise ${\zeta}_{s}(0,{\sigma}_{bd})$ which is independent of l_{b}, we find,
Note that we chose size additive division timing noise (${\zeta}_{s}(0,{\sigma}_{bd})$) for convenience in this section. However, it can be shown as done previously that the model is robust to the noise in division timing being size additive or time additive. Assuming that the noise terms ${\xi}_{lin}(0,C{V}_{\lambda ,lin})$ and ${\zeta}_{s}(0,{\sigma}_{bd})$ are small, we obtain to first order,
The terms l_{b}, ${\zeta}_{s}(0,{\sigma}_{bd})$ and ${\xi}_{lin}(0,C{V}_{\lambda ,lin})$ are independent of each other. The standard deviation of ${T}_{d}$ (${\sigma}_{t}$) can be calculated to be,
Assuming perfectly symmetric division and using ${l}_{d}^{n}=g({l}_{b}^{n})+{\zeta}_{s}(0,{\sigma}_{bd})$, we find the recursive relation for ${l}_{b}^{n}$ to be,
Note that Equation 43 is the same as Equation 6 under the approximation ${x}_{n}={l}_{b}^{n}1$. At steady state, the standard deviation of l_{b} is denoted by ${\sigma}_{b}$ and using Equation 43 its value is obtained to be,
Similarly, the standard deviation of l_{d}l_{b}, or equivalently ${\lambda}_{lin}{T}_{d}$, denoted by ${\sigma}_{l,lin}$, is calculated to be,
For linear growth, a natural plot is l_{d}l_{b} vs $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ (reminiscent of the $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ plot for exponential growth). To calculate the slope of the best linear fit, we have to calculate the correlation coefficient ${\rho}_{lin}$ given by,
Again using the independence of terms l_{b}, ${\zeta}_{s}(0,{\sigma}_{bd})$ and ${\xi}_{lin}(0,C{V}_{\lambda ,lin})$ from each other, we get,
The slope of best linear fit for the plot ${l}_{d}{l}_{b}$ vs $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ is given by,
The intercept ${c}_{tl,lin}$ is found to be,
On flipping the axis, the slope (${m}_{lt,lin}$) for the plot $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ vs ${l}_{d}{l}_{b}$ is obtained to be,
The intercept ${c}_{lt,lin}$ is found to be,
The best linear fit for the $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ vs ${l}_{d}{l}_{b}$ plot follows the trend y = x.
Simulations of the adder model for linearly growing cells were carried out. The deviation of the best linear fit for the ${l}_{d}{l}_{b}$ vs $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ plot from the y = x line is shown in Figure 3—figure supplement 1A, while in Figure 3—figure supplement 1B, the best linear fit for the plot $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ vs ${l}_{d}{l}_{b}$ is shown to agree with the y = x line.
Differentiating linear from exponential growth
Request a detailed protocolIn this section, we explore the equation for the best linear fit of $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ vs ${l}_{d}{l}_{b}$ plot in the case of exponential growth and $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ plot for linear growth. Intuitively, we expect the best linear fit in both cases to deviate from the y = x line. In this section, we will calculate the best linear fit explicitly. Surprisingly, we will find that, in the case of linear growth, the best linear fit for the $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ plot follows the y = x line closely.
Let us begin with exponential growth with growth rate, $\lambda $ = $\u27e8\lambda \u27e9$ + $\u27e8\lambda \u27e9\xi (0,C{V}_{\lambda})$ as defined previously. Again, $\xi (0,C{V}_{\lambda})$ has a normal distribution with mean 0 and standard deviation $C{V}_{\lambda}$, it being the CV of the growth rate. The time at division is given by Equation 7. The average growth rate $\u27e8\lambda \u27e9$ = $\u27e8\frac{\mathrm{ln}(2)}{{T}_{d}}\u27e9\approx \frac{\mathrm{ln}(2)}{\u27e8{T}_{d}\u27e9}$. For exponential growth, we will plot $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ vs ${l}_{d}{l}_{b}$. As previously defined, $\u27e8{\lambda}_{lin}\u27e9$ is the mean normalized elongation speed and $\u27e8{\lambda}_{lin}\u27e9=\u27e8\frac{1}{{T}_{d}}\u27e9\approx \frac{1}{\u27e8{T}_{d}\u27e9}$. $\u27e8{\lambda}_{lin}\u27e9$ is related to $\u27e8\lambda \u27e9$ by,
${l}_{d}{l}_{b}$can be calculated by using the regulation strategy f(l_{b}) introduced in the Model section and a normally distributed size additive noise ${\zeta}_{s}(0,{\sigma}_{bd})$. Note that we have chosen the noise in division timing to be size additive. However, the model is robust to the choice of type of noise as we showed in the Exponential growth section. Using Equations 5 and 6 we obtain,
Using Equation 11, $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ is obtained to be,
To calculate the expression for ${m}_{lt,lin}$, the slope of the best linear fit for $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ vs ${l}_{d}{l}_{b}$ plot, we first calculate ${\rho}_{lin}$ given by Equation 46. The expression for ${\sigma}_{l,lin}$ (standard deviation of ${l}_{d}{l}_{b}$) and ${\sigma}_{t}$ (standard deviation of ${T}_{d}$) are found to be,
${\sigma}_{x}$is related to ${\sigma}_{n}$ via Equation 17. In Exponential growth section, we also showed that ${\sigma}_{n}$ = $\frac{{\sigma}_{bd}}{2}$. Using these, we can write,
Now using the expressions for ${\sigma}_{t}$, ${\sigma}_{l,lin}$ and the fact that $x$, $\xi (0,C{V}_{\lambda})$ and ${\zeta}_{s}(0,{\sigma}_{bd})$ are independent of each other, we get,
For the plot $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ vs ${l}_{d}{l}_{b}$, the slope ${m}_{lt,lin}$ is given by,
Inserting Equation 55 into Equation 59 and substituting ${\sigma}_{x}^{2}$ given by Equation 57, we obtain,
The intercept ${c}_{lt,lin}$ is found to be,
For the adder model ($\alpha $ = $\frac{1}{2}$), we get the value of slope ${m}_{lin,lt}$ = $\frac{1}{2\mathrm{ln}(2)}\approx 0.7213$ and intercept ${c}_{lin,lt}$ = $1\frac{1}{2\mathrm{ln}(2)}\approx 0.279$. This is different from the best linear fit obtained for same regulatory mechanism controlling division in linearly growing cells where we found that the best linear fit follows the y = x line. Intuitively, we expect the best linear fit of $\u27e8{\lambda}_{lin}\u27e9{T}_{d}$ vs ${l}_{d}{l}_{b}$ plot to deviate from y = x line in the case of exponential growth. We showed analytically that for a class of models where birth controls division, it is indeed the case. This is also shown using simulations of the adder model in Figure 3—figure supplement 1C.
In the ‘ Obtaining the best linear fit’ section, we found the best linear fit for $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ plot to follow the y = x line for exponentially growing cells where division is regulated by birth event via regulation strategy f(l_{b}). Next, we calculate the equation for the best linear fit of $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ plot given that growth is linear. The model for division control will be same as that in the Linear growth section that is, the regulation strategy for division is given by g(l_{b}) = $2+2(1\alpha )({l}_{b}1)$ which is also equivalent to f(l_{b}). The linearly growing cells grow with elongation speed ${\lambda}_{lin}$ = $\u27e8{\lambda}_{lin}\u27e9(1+{\xi}_{lin}(0,C{V}_{\lambda ,lin}))$. As discussed before, ${\xi}_{lin}(0,C{V}_{\lambda ,lin})$ has a normal distribution with mean 0 and standard deviation $C{V}_{\lambda ,lin}$, it being the CV of the elongation speed. Using Equations 5 and 6, we get,
Using Equations 5 and 52, we obtain from Equation 41,
Since $x$, ${\xi}_{lin}(0,C{V}_{\lambda ,lin})$ and ${\zeta}_{s}(0,{\sigma}_{bd})$ are uncorrelated, the standard deviation of $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ and ${T}_{d}$ denoted by ${\sigma}_{l}$ and ${\sigma}_{t}$ respectively are calculated to be,
We calculate the correlation coefficient for the pair ($\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$, $\u27e8\lambda \u27e9{T}_{d}$). Since the correlation coefficient is unaffected by multiplying one of the variables with a positive constant, we can calculate the correlation coefficient for the pair ($\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$, ${T}_{d}$) or ${\rho}_{exp}$ as given by Equation 18. Using the independence of terms $x$, ${\xi}_{lin}(0,C{V}_{\lambda ,lin})$ and ${\zeta}_{s}(0,{\sigma}_{bd})$,
For the plot $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$, the slope ${m}_{lt}$ of the best linear fit is given by,
Inserting Equation 64 into Equation 67 and using Equation 57, we get,
Similarly the intercept (${c}_{lt}$) for the plot $\u27e8\lambda \u27e9{T}_{d}$ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ is found to be,
This is very close to y = x trend obtained for the same regulatory mechanism controlling division in exponentially growing cells (Figure 3A).
Growth rate vs age and elongation speed vs age plots
Request a detailed protocolIn the previous sections, we found that binning and linear regression on the plot $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$, and the plot obtained by interchanging the axes, were inadequate to identify the mode of growth. In this section, we try to validate the growth rate vs age plot as a method to elucidate the mode of growth.
In addition to cell size at birth and division and the generation time, cell size trajectories (cell size, $L$ vs time from birth, $t$) were obtained for multiple cell cycles. In our case, the cell size trajectories were collected either via simulations (in Figure 3B) or from experiments (for Figure 4AC) at intervals of 4 min. Note that if the measurements were to be carried out at equal length intervals instead of time, the results discussed in the paper would still remain unchanged. For each trajectory, growth rate at time $t$ or age $\frac{t}{{T}_{d}}$ is calculated as $\frac{1}{L(t)}\frac{L(t+\mathrm{\Delta}t)L(t)}{\mathrm{\Delta}t}$ where $\mathrm{\Delta}t$ is the time between consecutive measurements. To obtain elongation speed vs age plots, the formula before needs to be replaced with $\frac{L(t+\mathrm{\Delta}t)L(t)}{\mathrm{\Delta}t}$. The growth rate is interpolated to contain 200 points at equal intervals of time for each cell trajectory. The growth rate trends appear to be robust with regards to a different number of interpolated points (from 100 to 500 points). To obtain the growth rate trend as a function of cell age, we use the method previously applied in Nordholt et al., 2020. In this method, growth rate is binned based on age for each individual trajectory (50 bins) and the average growth rate is obtained in each of the bins. The binned data trend for growth rate vs age is then found by taking the average of the growth rate in each bin over all trajectories. Binning the growth rate for each trajectory ensures that each trajectory has an equal contribution to the final growth rate trend so as to avoid inspection bias. This step is especially important when data collected at equal intervals of time is analyzed. In such a case, cells with larger generation times have a greater number of measurements than cells with smaller generation times. Obtaining the growth rate trend without binning growth rate for each trajectory would have biased the binned data trend for the growth rate vs age plot to a smaller value because of overrepresentation by slowergrowing cells (or equivalently cells with longer generation time). This bias toward lower growth rate values in the growth rate vs age plots is an instance of inspection bias.
In Figure 4AC, we find the growth rate obtained from E. coli experiments to change within the cell cycle. In the two slower growth media (Figure 4A and B), the growth rate is found to increase with cell age while for the fastest growth media (Figure 4C) the growth rate follows a nonmonotonic behaviour similar to that observed in Nordholt et al., 2020 for B. subtilis. Abrupt changes in growth rate are reported at constriction in Reshes et al., 2008; Banerjee et al., 2017. We find that the growth rate changes start before constriction in the two slower growth conditions considered. One possibility is that this increase is due to preseptal cell wall synthesis (Pazos et al., 2018). Preseptal cell wall synthesis does not require activity of PBP3 (FtsI) but instead relies on bifunctional glycosyltransferases PBP1A and PBP1B that link to FtsZ via ZipA. One hypothesis that can be tested in future works is that at the onset of constriction, activity from PBP1A and PBP1B starts to gradually shift to the PBP3/FtsW complex and therefore no abrupt change in growth rate is observed. In the fastest growth condition (glucosecas medium), we find that the increase in growth rate approximately coincides with onset of constriction, in agreement with the previous findings (Reshes et al., 2008; Banerjee et al., 2017).
In Figure 4AC, the growth rate trends are not obtained for age close to one. This is because growth rate at age = 1 is given by $\frac{1}{L({T}_{d})}\frac{L({T}_{d}+\mathrm{\Delta}t)L({T}_{d})}{\mathrm{\Delta}t}$ and this requires knowing the cell lengths beyond the division event ($L({T}_{d}+\mathrm{\Delta}t)$). To estimate growth rates at age close to one, we approximate $L({T}_{d}+\mathrm{\Delta}t)$ to be the sum of cell sizes of the two daughter cells. In order to minimize inspection bias, we considered only those cell size trajectories which had $L(t)$ data for 12 min after division (corresponding to an age of approximately 1.1). However, the growth rate trends in all three growth media were robust with regards to a different time for which $L(t)$ was considered (4 min to 20 min after division). We use the binning procedure discussed before in this section. To validate this method, we applied it on synthetic data obtained from the simulations of exponentially growing cells following the adder and the adder per origin model. Cells were assumed to divide in a perfectly symmetric manner and both of the daughter cells were assumed to grow with the same growth rate, independent of the growth rate in the mother cell. The growth rate trends for the two models considered (adder and adder per origin) are expected to be constant even for cell age >1. We found that the growth rate trends were indeed approximately constant as shown in Figure 4—figure supplement 1D. We also considered linear growth with division controlled via an adder model. The daughter cells were assumed to grow with the same elongation speed, independent of the elongation speed in the mother cell. In this case, we expect the elongation speed trend to be constant for cell age >1. This is indeed what we observed as shown in the inset of Figure 4—figure supplement 1D. We used this method on E. coli experimental data and found that the growth rate trends obtained for the three growth conditions (Figure 4—figure supplement 1A–1C) were consistent with that shown in Figure 4AC in the relevant age ranges. For cell age close to one, we found that the growth rate decreased to a value close to the growth rate near cell birth (age ≈ 0) for all three growth conditions considered.
In summary, we find that the growth rate vs age plots are a consistent method to probe the mode of cell growth within a cell cycle.
Growth rate vs time from specific event plots are affected by inspection bias
Request a detailed protocolTo probe the growth rate trend in relation to a specific cell cycle event, for example cell birth, growth rate vs time from birth plots are obtained for simulations of exponentially growing cells following the adder model. In the growth rate vs time from birth plot, the rate is found to stay constant and then decrease at longer times (Figure 3—figure supplement 2C) even though cells are exponentially growing. Because of inspection bias (or survivor bias), at later times, only the cells with larger generation times (or slower growth rates) ‘survive’. The average generation time of the cells averaged upon in each bin of Figure 3—figure supplement 2C is shown in Figure 3—figure supplement 2D. The decrease in growth rate in Figure 3—figure supplement 2C occurs around the same time when an increase in generation time is observed in Figure 3—figure supplement 2D. Thus, the trend in growth rate is biased toward lower values at longer times. The problem might be circumvented by restricting the time on the xaxis to the smallest generation time of all the cell cycles considered (Messelink et al., 2020).
To check for growth rate changes at constriction, we used plots of growth rate vs time from constriction ($t{T}_{n}$). Growth rate trends obtained from E. coli experimental data show a decrease at the edges of the plots (Figure 4—figure supplement 2A and C, and 2E). These deviate from the trends obtained using the growth rate vs age plots (Figure 4AC). To investigate this discrepancy, we use a model which takes into account the constriction and the division event. Currently, it is unknown how constriction is related to division. For the purpose of methods validation, we use a model where cells grow exponentially, constriction occurs after a constant size addition from birth, and division occurs after a constant size addition from constriction. Note that other models where constriction occurs after a constant size addition from birth while division occurs after a constant time from constriction, as well as a mixed timeradder model proposed in Banerjee et al., 2017, lead to similar results. We expect the growth rate trend to be constant for exponentially growing cells. However, we find using numerical simulations that it decreases at the plot edges both before and after the constriction event (Figure 3—figure supplement 2A). This decrease can be attributed to inspection bias. The average growth rate in time bins at the extremes are biased by cells with smaller growth rates. This is shown in Figure 3—figure supplement 2B where the average generation time for the cells contributing in each of the bins of Figure 3—figure supplement 2A is plotted. The time at which the growth rate decreases on both sides of the constriction event is close to the time at which the average generation time increases. For example, in alanine medium, the generation time for each of the bins is plotted in Figure 4—figure supplement 2B. The average generation time for the cells contributing to each of the bins is almost constant for the timings between –80 min and 20 min. Thus, for this time range the changes in growth rate are not because of inspection bias but are a real biological effect. The behavior of growth rate within this time range in Figure 4—figure supplement 2A is in agreement with the trend in growth rate vs age plot of Figure 4A. On accounting for inspection bias, the growth rate vs age plots agree with the growth rate vs time from constriction plots in other growth media as well (Figure 4—figure supplement 2C, Figure 4—figure supplement 2E). Thus, growth rate vs time plots are also a consistent method to probe growth rate modulation in the time range when avoiding the regimes prone to inspection bias.
Results of elongation speed vs size plots are modeldependent
Request a detailed protocolCells assumed to undergo exponential growth have elongation speed proportional to their size. In the case of exponential growth, the binned data trend of the plot elongation speed vs size is expected to be linear with the slope of the best linear fit providing the value of growth rate and intercept being zero. In this section, we use the simulations to test if binning and linear regression on the elongation speed vs size plots are suitable methods to differentiate exponential growth from linear growth (Cadart et al., 2019).
To test the method, we generate cell size trajectories using simulations of the adder model with a size additive division timing noise and assuming exponential growth. Elongation speed at size $L(t)$ is calculated for each trajectory as $\frac{L(t+\mathrm{\Delta}t)L(t)}{\mathrm{\Delta}t}$ where $\mathrm{\Delta}t$ is the time between consecutive measurements ( = 4 min in our case). Each trajectory is binned into 10 equally sized bins based on their cell sizes and the average elongation speed is obtained for each bin. The final trend of elongation speed as a function of size is then obtained by binning (based on size) the pooled average elongation speed data of all the cell cycles.
We find that the binned data trend is linear with the slope of the best linear fit close to the average growth rate considered in the simulations (Figure 3—figure supplement 3D). This is in agreement with our expectations for exponential growth. In order to check if this method could differentiate between exponential growth and linear growth, we used simulations of the adder model undergoing linear growth to generate cell size trajectories for multiple cell cycles. For linear growth, elongation speed is expected to be constant, independent of its cell size. The binned data trend for the elongation speed vs size plot is also obtained to be constant for the simulations of linearly growing cells (Figure 3—figure supplement 3B). The intercept of the best linear fit obtained is close to the average elongation speed considered in the simulations. The binned data trend for linear and exponential growth are clearly different as shown in Figure 3—figure supplement 3B and Figure 3—figure supplement 3D, respectively, and this result holds for a broad class of models where the division event is controlled by birth and the growth rate (for exponential growth)/elongation speed (for linear growth) is distributed normally and independently between cellcycles.
Next, we consider the adder per origin cell cycle model for exponentially growing cells (Ho and Amir, 2015). In this model space, the cell initiates DNA replication by adding a constant size per origin from the previous initiation size. The division occurs on average after a constant time from initiation. For exponentially growing cells, the binned data trend is still expected to be linear as before. Instead, we find using simulations that the trend is nonlinear and it might be misinterpreted as nonexponential growth (Figure 3—figure supplement 3F).
Thus, the results of binning and linear regression for the plot elongation speed vs size is modeldependent.
Interchanging axes in growth rate vs inverse generation time plot might lead to different interpretations
Request a detailed protocolSo far, our discussion was focused on the question of mode of singlecell growth. A related problem regards the relation between growth rate ($\lambda $) and the inverse generation time ($\frac{1}{{T}_{d}}$). On a population level, the two are clearly proportional to each other. However, singlecell studies based on binning showed an intriguing nonlinear dependence between the two, with the two variables becoming uncorrelated in the fastergrowth media (Kennard et al., 2016; IyerBiswas et al., 2014). Within the same medium, the binned data curve for the plot $\lambda $ vs $\frac{1}{{T}_{d}}$ flattened out for faster dividing cells. The trend in the binned data was different from the trend of y = $\mathrm{ln}(2)$ x line as observed for the population means. A priori one might speculate that the flattening in faster dividing cells could be because the faster dividing cells might have less time to adapt their division rate to transient fluctuations in the environment. Kennard et al., 2016 insightfully also plotted $\frac{1}{{T}_{d}}$ vs $\lambda $ and found a collapse of the binned data for all growth conditions onto the y = $\frac{1}{\mathrm{ln}(2)}$ x line. These results are reminiscent of what we previously showed for the relation of $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ and $\u27e8\lambda \u27e9{T}_{d}$.
In the following, we will elucidate why this occurs in this case using an underlying model and predicting the trend based on it. We use simulations of the adder model undergoing exponential growth. The parameters for size added in a cell cycle and mean growth rates are extracted from the experimental data. CV of growth rate is assumed lower in faster growth media as observed by Kennard et al. Using this model, we could obtain the same pattern of flattening at faster growth conditions that is observed in the experiments (Figure 2—figure supplement 2A). The population mean for $\lambda $ and $\frac{1}{{T}_{d}}$ follows the expected y = $\mathrm{ln}(2)$ x equation (shown as black dashed line) as was the case in experiments. Intuitively, such a departure from the expected y = $\mathrm{ln}(2)$ x line for the singlecell data can again be explained by determining the effect of noise on variables plotted on both axes. As previously stated ${T}_{d}$ is affected by both growth rate noise and noise in division timing while growth rate fluctuates independently of other sources of noise. This does not agree with the assumption for binning as noise in division timing affects the xaxis variable rather than the yaxis variable. In such a case, the trend in the binned data might not follow the expected y = $\mathrm{ln}(2)$ x line. However, on interchanging the axes, we would expect the assumptions of binning to be met and the trend to follow the y = $\frac{1}{\mathrm{ln}(2)}$ x line (Figure 2—figure supplement 2B).
Data and simulations
Experimental data
Request a detailed protocolExperimental data obtained by Tanouchi et al., 2017 was used to plot ${L}_{d}$ vs ${L}_{b}$ shown in Figure 1A. E. coli cells were grown at 25 °C in a mother machine device and the length at birth and division were collected for multiple cell cycles. ${L}_{d}$ vs ${L}_{b}$ plot was obtained using these cells and linear regression performed on it provided a best linear fit.
Data from recent mother machine experiments on E. coli was used to make all other plots. Details are provided in the Experimental methods and TiruvadiKrishnan et al., 2021. The experiments were conducted at 28 °C in three different growth conditions  alanine, glycerol, and glucosecas (also see Experimental methods). Cell size trajectories were collected for multiple cell cycles and all of the data collected were considered while making the plots in the paper.
Simulations
Request a detailed protocolMATLAB R2021a was used for simulations. Simulations of the adder model for exponentially growing cells were carried out over a single lineage of 2500 generations (Figure 2C and D, Figure 3—figure supplement 1C). The mean length added between birth and division was set to 1.73 $\mu m$ in line with the experimental results for alanine medium. Growth rate was variable and sampled from a normal distribution at the start of each cell cycle. The mean growth rate was set to $\frac{\mathrm{ln}(2)}{\u27e8{T}_{d}\u27e9}$, where $\u27e8{T}_{d}\u27e9$ = 212 min and coefficient of variation (CV) = $C{V}_{\lambda}$ = 0.15. The noise in division timing was assumed to be time additive with mean 0 and standard deviation $\frac{{\sigma}_{n}}{\u27e8\lambda \u27e9}$, where ${\sigma}_{n}$ = 0.15. The binning data trends and the best linear fits obtained using these simulations could be compared with the analytical results obtained in sections ‘Nonlinearity in binned data’ and ‘Differentiating linear from exponential growth’.
For simulations of linear growth (Figure 3AB, Figure 3—figure supplement 1A and B, Figure 3—figure supplement 3A, B , Figure 4—figure supplement 1D), the mean growth rate was set to $\frac{\u27e8{L}_{d}{L}_{b}\u27e9}{\u27e8{T}_{d}\u27e9}$, with the values of $\u27e8{L}_{d}{L}_{b}\u27e9$ and $\u27e8{T}_{d}\u27e9$ used as mentioned previously. The noise in division timing was size additive with standard deviation = 0.15 $\u27e8{L}_{b}\u27e9$. Noise was also considered to be size additive with the same standard deviation for the simulations of exponentially growing cells shown in Figure 3B, Figure 3—figure supplement Figure 3—figure supplements 2C and 3C. D, and Figure 4—figure supplement 1D.
In the simulations of superexponential growth carried over a single lineage of 2500 generations (Figure 3B), the cells initially grew exponentially but in the later stages of the cell cycle, the growth rate increased as,
where k was fixed to be $\frac{2}{{T}_{d}^{3}}$ and t_{c} was the time from birth at which the growth rate changed from exponential to superexponential growth. t_{c} was fixed to be half of the generation time of the cell or equivalently an age of 0.5. The division size was set by the adder model with a time additive noise with similar parameters as before for exponential growth simulations. The exponential growth rate at the start of each cell cycle was drawn from a normal distribution with mean set to $\frac{\mathrm{ln}(2)}{242}mi{n}^{1}$ and CV = 0.15.
For Figure 3B, Figure 3—figure supplement 3E and F, Figure 4—figure supplement 1D, simulations were carried out over a lineage of 2500 generations for exponentially growing cells following the adder per origin model. In the simulations, the time increment is 0.01 min. The initial condition for the simulations is that cells are born and initiate DNA replication at time t = 0 but the results are independent of initial conditions. The number of origins is also tracked throughout the simulations beginning with an initial value of 2. Cells divide into two daughter cells in a perfectly symmetrical manner (no noise in division ratio), and one of the daughter cells is discarded for the next cell cycle. In simulations, the growth rate was fixed within a cell cycle but varied between different cell cycles. On division, the growth rate for that cell cycle was drawn from a normal distribution with mean $\u27e8\lambda \u27e9$ and coefficient of variation ($C{V}_{\lambda}$) whose values were fixed using the experimental data from alanine medium. The total length at which the next initiation happens is determined by,
where ${\mathrm{\Delta}}_{ii}$ is the length added per origin and O is the number of origins. To determine ${L}_{i}^{tot,next}$, ${\mathrm{\Delta}}_{ii}$ was drawn on reaching initiation length from a normal distribution. The mean and CV of ${\mathrm{\Delta}}_{ii}$ was obtained from experiments done in alanine medium. In the adder per origin model, division happens after a C + D time from initiation. The division length (${L}_{d}$) is obtained to be,
In the simulations, once the initiation length was reached, the corresponding division occurred a time C + D after initiation. C + D timings for each initiation event were again drawn from a normal distribution with the same mean and CV as that of the experiments in alanine medium.
For Figure 3—figure supplement 2A, cells were assumed to grow exponentially in the simulations. The constriction length (${L}_{n}$) was set to be,
The length added (${\mathrm{\Delta}}_{bn}$) was assumed to have a normal distribution with the mean length added between birth and constriction set to 1.18 $\mu m$ and the CV = 0.23, in line with the experimental results for alanine medium. The length at division was set as,
The length added (${\mathrm{\Delta}}_{nd}$) was also assumed to have a normal distribution with the mean length added set to 0.53 µm and the CV = 0.26, again in line with the experimental results for alanine medium.
For Figure 3B, Figure 3—figure supplement 2A–2D, 3 A3F, Figure 4—figure supplement 1D, the cell sizes are recorded within the cell cycle at equal intervals of 4 min, similar to that in the E. coli experiments of TiruvadiKrishnan et al., 2021.
For simulations shown in Figure 4—figure supplement 1D, the cell size trajectories are obtained at intervals of 4 min beyond the current cellcycle. The size after the division event is said to be the sum of the sizes of the daughter cells. It is also further assumed that the daughter cells are equal in size (perfectly symmetric division) and they both grow with the same growth rate (for exponential growth) or elongation speed (for linear growth). The growth rates/elongation speeds for the daughter cells are sampled from a normal distribution with a mean and CV as discussed before. The cell size trajectories are recorded for 80 min after the division event in the current cell cycle.
In Figure 2—figure supplement 2, simulations of the adder model for exponentially growing cells were carried out until a population of 5000 cells was reached. The parameters for size added in a cell cycle and mean growth rates were extracted from the experimental data (Kennard et al., 2016). The value of ${\sigma}_{n}$ used in all growth conditions was 0.17 while $C{V}_{\lambda}$ decreased in faster growth conditions (0.2 in the three slowest growth conditions, 0.12 and 0.07 in the second fastest and fastest growth conditions, respectively).
Appendix 1
Comparing length, surface area and volume growth rate
In the paper, we use cell length to represent cell size. However, other cell size characteristics such as cell surface area and cell volume could also be used to denote cell size. How does the growth rate vary with our choice of cell length, cell surface area, or cell volume to be the cell size?
To study this, we assume a cell morphology as shown in Appendix 1—figure 1. We assume that E. coli cells are cylindrical with hemispherical poles. The total length of the cell is $L$ with a radius $R$. The cell volume ($V$) is then,
The morphology of the cell after constriction is also shown in Appendix 1—figure 1A. The volume in this case is,
If we make the assumption that cell biomass grows exponentially and the total cell surface area is coupled to the biomass (Oldewurtel et al., 2021), then cell surface area grows exponentially with time. Using the morphology in Appendix 1—figure 1A, the total surface area ($S$) before and after constriction is,
Surprisingly, this is independent of h. Since the surface area is proportional to the cell length (Equation (3A1)), the length growth is also exponential with an identical growth rate as surface area growth, assuming the width of the cell is constant. The exponential growth of cell length is shown in Appendix 1—figure 1B using simulations where the cell surface is assumed to grow exponentially. So, for this model of cell growth and morphology, the length and the surface growth rates are found to be identical.
Next, we compare length growth rate to volume growth rate considering the same cell morphology as that in Appendix 1—figure 1A. In this model, the volume growth is assumed to be exponential. The volume before and after the onset of constriction are given by Equations (1A1) and (2A1), respectively.
Before constriction, the volume grows only by an increase in length of the cylindrical part of the cell while the width stays constant. However, after the constriction at midcell starts, the volume grows by an increase in length as well as by adding a septum surface at the midcell. We assume that the septum wall surface grows at a constant rate (c_{1}) (Reshes et al., 2008). We can obtain c_{1} in terms of cell morphology variables to be,
We can solve for $h(t)$ using the following boundary conditions,
where $T}_{n$ is the time from birth at which constriction starts. Using Equation (4A1) and Equation (5A1), we can obtain c_{1} in terms of cell cycle variables $R$, $T}_{n$ and $T}_{d$,
Under these assumptions, for exponential volume growth, we obtain the length growth via simulations. The length growth rate is shown in Appendix 1—figure 1C. The growth rate, the length at birth, the time at constriction from birth and the generation time parameters used in the simulations are obtained from experimental data in alanine growth medium. The width of the cells is assumed to be 0.35 $\mu m$. We find that before constriction, the length growth rate increases to a small extent ($\approx 6\%$). However, after constriction there is a rapid increase in length growth rate. Since the length growth rate increases with age, the length growth is superexponential. Thus, the mode of growth in length and volume are not identical.
Appendix 2
Linear regression on $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8{T}_{d}\u27e9\lambda $ plot and its interchanged axes plot
In section ‘Statistical methods like binning and linear regression should be interpreted based on a model’, we found that binning and linear regression on the plots $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8\lambda \u27e9{T}_{d}$ and its interchanged axes were not a suitable method to identify the underlying mode of growth. In this section, we explore binning and linear regression on similar plots $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8{T}_{d}\u27e9\lambda $ plot and its interchanged axes. We test the usability of these plots to elucidate the mode of growth using the methodology proposed in the paper.
Assuming exponential growth, $\lambda $ for a cell cycle can be calculated as $\frac{1}{{T}_{d}}\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$. On plotting $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8{T}_{d}\u27e9\lambda $ (Appendix 2—figure 1 ) and $\u27e8{T}_{d}\u27e9\lambda $ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ (Appendix 2—figure 1) for the experimental data, we obtain the slope of the best linear fit to be close to zero (values shown in Table 1Appendix 2). Next, using the methodology of the paper, we interpret these results using an underlying model. We consider a model in which cells grow exponentially with the division determined by birth. In the model, growth rate is fixed at the beginning of each cell cycle and is independent of size at birth. The model predicts that $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ will be independent of the growth rate (Equation 19 in main text). Thus, we would expect the slope to be zero for both of the plots $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8{T}_{d}\u27e9\lambda $ and $\u27e8{T}_{d}\u27e9\lambda $ vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$. This is also shown using simulations of the adder model in Appendix 2—figure 1 where the slope of the plots is close to zero. In order to differentiate between exponential growth and linear growth, the best linear fit in case of linear growth for these plots must deviate from y = constant line. However, we find for the simulations of the adder model where cells grow linearly that the slope of the best linear fit for both of the above plots is still zero (Appendix 2—figure 1). Note that $\lambda $ in the case of linear growth is still calculated as $\frac{1}{{T}_{d}}\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$. A slope of zero in case of linear growth can be explained using Equation 62 of the main text. Using the equation, we find that $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ is independent of the underlying growth rate for linear growth. Thus, the best linear fit for both plots have a slope of zero in the case of linear growth. This indicates that binning and linear regression on the $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs $\u27e8{T}_{d}\u27e9\lambda $ and its interchanged axes plots are unsuitable for elucidating the mode of growth.
Data availability
All data generated during this study are deposited in Dataverse (https://doi.org/10.7910/DVN/BNQUDW).

Harvard DataverseDistinguishing different modes of growth using singlecell data.https://doi.org/10.7910/DVN/BNQUDW

figshareData from longterm growth data of Escherichia coli at a singlecell level.https://doi.org/10.6084/m9.figshare.c.3493548.v1
References

Singlecell model of prokaryotic cell cycleJournal of Theoretical Biology 341:78–87.https://doi.org/10.1016/j.jtbi.2013.09.035

Cell Size Regulation in BacteriaPhysical Review Letters 112:208102.https://doi.org/10.1103/PhysRevLett.112.208102

Learning from Noise: How Observing Stochasticity May Aid MicrobiologyTrends in Microbiology 26:376–385.https://doi.org/10.1016/j.tim.2018.02.003

Details Matter: Noise and Model Structure Set the Relationship between Cell Size and Cell Cycle TimingFrontiers in Cell and Developmental Biology 5:92.https://doi.org/10.3389/fcell.2017.00092

The physics of cellsize regulation across timescalesNature Physics 15:993–1004.https://doi.org/10.1038/s415670190629y

Highthroughput measurement of singlecell growth rates using serial microfluidic mass sensor arraysNature Biotechnology 34:1052–1059.https://doi.org/10.1038/nbt.3666

CellSize Homeostasis and the Incremental Rule in a Bacterial PathogenBiophysical Journal 109:521–528.https://doi.org/10.1016/j.bpj.2015.07.002

Controlling cell size through sizer mechanismsCurrent Opinion in Systems Biology 5:86–92.https://doi.org/10.1016/j.coisb.2017.08.010

Using buoyant mass to measure the growth of single cellsNature Methods 7:387–390.https://doi.org/10.1038/nmeth.1452

Simultaneous regulation of cell size and chromosome replication in bacteriaFrontiers in Microbiology 6:662.https://doi.org/10.3389/fmicb.2015.00662

Fundamental principles in bacterial physiologyhistory, recent progress, and the future with focus on cell size control: a reviewReports on Progress in Physics. Physical Society 81:056601.https://doi.org/10.1088/13616633/aaa628

Bilinear cell growth of Escherichia coliJournal of Bacteriology 148:730–733.https://doi.org/10.1128/jb.148.2.730733.1981

Homeostasis of protein and mRNA concentrations in growing cellsNature Communications 9:1–11.https://doi.org/10.1038/s4146701806714z

Concurrent processes set E. coli cell divisionScience Advances 4:eaau3324.https://doi.org/10.1126/sciadv.aau3324

Single cell studies of the cell cycle and some modelsTheoretical Biology & Medical Modelling 2:4.https://doi.org/10.1186/1742468224

Step by Step, Cell by Cell: Quantification of the Bacterial Cell CycleTrends in Microbiology 25:250–256.https://doi.org/10.1016/j.tim.2016.12.005

Cell shape dynamics in Escherichia coliBiophysical Journal 94:251–264.https://doi.org/10.1529/biophysj.107.104398

Mechanistic Origin of CellSize Control and Homeostasis in BacteriaCurrent Biology 29:1760–1770.https://doi.org/10.1016/j.cub.2019.04.062

The Interpretation of Interaction in Contingency TablesJournal of the Royal Statistical Society 13:238–241.https://doi.org/10.1111/j.25176161.1951.tb00088.x

Sampling Bias and the Inspection ParadoxMathematics Magazine 58:96–99.https://doi.org/10.1080/0025570X.1985.11977160

Cellsize control and homeostasis in bacteriaCurrent Biology 25:385–391.https://doi.org/10.1016/j.cub.2014.12.009

E. coli genome manipulation by p1 transductionCurrent Protocols in Molecular Biology 79:1–17.

Robust growth of Escherichia coliCurrent Biology 20:1099–1103.https://doi.org/10.1016/j.cub.2010.04.045

Sizing up the bacterial cell cycleNature Reviews. Microbiology 15:606–620.https://doi.org/10.1038/nrmicro.2017.79

Analysis of Factors Limiting Bacterial Growth in PDMS Mother Machine DevicesFrontiers in Microbiology 9:871.https://doi.org/10.3389/fmicb.2018.00871
Decision letter

Gordon J BermanReviewing Editor; Emory University, United States

Aleksandra M WalczakSenior Editor; École Normale Supérieure, France

Benjamin P BrattonReviewer; Princeton University, United States
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Decision letter after peer review:
Thank you for submitting your article "To bin or not to bin: analyzing singlecell growth data" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Aleksandra Walczak as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Benjamin P Bratton (Reviewer #3).
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
1) The final result (Figure 4) is somewhat disconnected from the majority of the paper that precedes it. Specifically, the authors' procedure that resolves exponential vs. nonexponential growth results in E. coli in alanine being deemed exponential (Figure 2B) only to later be revealed as nonexponential (Figure 4A), albeit weakly. Furthermore, the procedure advertised as distinguishing exponential from linear growth (Figure 3B), when applied to the data, reveals neither (Figure 4). This makes the main point of the paper (the demonstration and resolution of pitfalls) feel disconnected from its application to a particular case, which is more nuanced and likely leaves many questions unanswered.
2) The title ("To bin or not to bin…") implies that binning is the main culprit behind potentially misleading analyses, but one of the reviewers argued that in the end, it is linear regression. Each of the two main pitfalls and their resolution would be unchanged if the data were never binned. Binning affects the apparent curvature of the y vs x relationship, but this reads as a more minor point. Therefore, the title may be a bit misleading in service of its poeticism.
3) The authors look at different choices of binning dimensions but do not sufficiently explore the power of their generative model to perform (un)weighted regression or parameter estimation from the notbinned data. They do explore the unbinned data from an analytical statistics approach in sections 5.4.1 and 5.5 but this is not yet extensively explored in the figures and/or discussion.
4) This manuscript could be improved by including a discussion on binning singlecell length trajectories, those taken on a single cell within one division cycle, either based on the length or on time. Experimentally, my understanding is that most growth measurements are done on a fixed \Δ t basis, and therefore measurements are averaged using length as a function of time. Another option, one that I have not seen used, is to use a \Δ L basis and measure the time to grow either a specific amount or a specific fractional amount. Of course, these are different ways to slice the same cake, and such an extension may be useful to expand the authors' argument that one needs to utilize a generative model to be able to assess their choice in statistical analysis.
5) As the formulation of this approach borrows from and is built heavily on previous data and previous statistical frameworks, it would be helpful to have more discussion on what is new to this particular manuscript. This is always a difficult balance as one needs the context of the previous work to understand the symbology used, but just reiterating what is already in the literature can lead to a dilution of the new/novel/important discoveries of the manuscript.
6) The authors should resolve the use of the term "cell size" to mean both "cell length" and "cell volume" (line 107, line 431/432). The authors briefly touch on a related concept of length growth vs biomass growth in 343357. As a first, naively simple model of cell growth, it is possible to come up with a simple model for linear vs volumetric growth by imagining that a cell grows exponentially in length from birth to death. In this simple model, at the appropriate time, it switches from growing a cylinder to growing hemispherically endcaps, as it needs to perfectly duplicate its shape. In this simple model, one can easily see that the instantaneous, relative volumetric and length growth rates are not identical; vary throughout the cell cycle; and even their average over the whole cell cycle is not 1. While the authors mention that the average diameter is quite robust between growth conditions, only in the limit that the cell is very long and not a spherocylinder does the approximation that shape does not matter come in to play. Granted, this naïve model does not incorporate the wealth of information known about division timing or geometric changes during growth, and underscores the authors' main point that a specific model of the dynamical process should be included.Reviewer #1:
In this manuscript, the authors describe a generative modelbased framework to better analyze stochastic growth data, including bacterial cell growth. They show how this framework can be applied to gain insight into the processes underlying these phenomena. More specifically, they start by showing how binning along different axes in the
This work is wellsupported by simulations and data analysis and will likely be of interest to those trying to understand the processes governing bacterial growth, as well as those studying stochastic growth processes in biology more broadly.
Strengths:
– The choice and execution of the simulations were sensible and welldone, respectively, and they provided clarity as to the overall message of the manuscript.
– The conclusions are wellsupported by the data.
– I found the writing to be clear throughout.
Weaknesses:
– It would be good to have a more extensive discussion about what is specifically new here. This is not my particular field, so having a bit more of an introduction about methods beyond binning (if any) that have emerged to understand these data.Reviewer #2:
The manuscript by Kar et al., uses singlecell experiments, simulations, and theory to investigate a common method of determining whether cell size grows exponentially. Specifically, they show that a relationship that should adhere to y = x for exponential growth on average (where x is the product of the division time and mean growth rate, and y is the log of the ratio of the birth and division size), in fact deviates from y = x with noise, because x contains noise that y does not (growth rate noise). This makes exponential growth seem nonexponential. The resolution, they show, is to plot x vs. y instead. Additionally, they show that when plotting x vs. y for linear growth, the relationship is coincidentally very close to x = y. This makes linear growth seem exponential. The resolution, they show, is to plot the instantaneous growth rate vs. the normalized cell age, which will decrease for linear growth and will be constant for exponential growth. Applying this protocol to E. coli size data, they find that the growth rate actually increases weakly with age, indicating somewhat superexponential growth.
Strengths
The insights in this manuscript are highly important for the field to know. The fact that exponential growth can masquerade as nonexponential growth and vice versa means that much confusion likely exists in a field that is already surprisingly complex given how simple its questions are to state. The fact that the authors offer resolutions to these pitfalls means that this work should add clarity to the field and move it forward in a meaningful way.
The conclusions are well supported by the data. For example, in the case of seemingly nonexponential growth, the problem is presented by the experimental data, reproduced by the simulations, and explained by the theory, while the resolution is inspired by the theory, proven by the simulations, and demonstrated by the experimental data.
The manuscript is well written. While it is rooted in careful analysis, it remains understandable to a largely nonquantitative audience.
Weaknesses
The final result (Figure 4) is somewhat disconnected from the majority of the paper that precedes it. Specifically, the authors' procedure that resolves exponential vs. nonexponential growth results in E. coli in alanine being deemed exponential (Figure 2B) only to later be revealed as nonexponential (Figure 4A), albeit weakly. Furthermore, the procedure advertised as distinguishing exponential from linear growth (Figure 3B), when applied to the data, reveals neither (Figure 4). This makes the main point of the paper (the demonstration and resolution of pitfalls) feel disconnected from its application to a particular case, which is more nuanced and likely leaves many questions unanswered.
The title ("To bin or not to bin…") implies that binning is the main culprit behind potentially misleading analysis, but I would argue that in the end, it is linear regression. Each of the two main pitfalls and their resolution would be unchanged if the data were never binned, I believe. Binning affects the apparent curvature of the y vs x relationship, but this reads as a more minor point. Therefore, the title may be a bit misleading in service of its poeticism.Reviewer #3:
Kar et al., examine an interesting and important question of how to make sense of large sets of observational data, specifically cell length data, which may or may not be consistent with various underlying biological mechanisms. As datasets improve in their technical quality (increasing spatiotemporal resolution, increasing numbers of observations), there is hope that the community will be able to resolve differences between underlying cell biological mechanisms of cell size homeostasis. As the authors point out, these interpretations and analyses require statistical analysis that can accurately perform the model selection or parameter estimation task of interest.
1) The authors succeed in bringing attention to the issue of appropriate binning when analyzing large datasets. The authors focus their figures and discussion on an important, and practical issue, as many researchers perform linear regression on binned data. The title and framing of the manuscript imply that it will provide a comparison with statistical methods that do not involve binning. The authors look at different choices of binning dimensions, but do not sufficiently explore the power of their generative model to perform (un)weighted regression or parameter estimation from the notbinned data. They do explore the unbinned data from an analytical statistics approach in section 5.4.1 and 5.5 but this not yet extensively explored in the figures and/or discussion.
2) The authors succeed in walking the reader through the power of examining a specific mechanism/model and the statistical properties of that model. In this case, the authors do this with both a model of cell length homeostasis that comes from exponential growth or linear growth with homeostatic feedback. This approach rests heavily on previous work from the same group including a reframing of the statistical correlations between different observables that was explored in the included references [13] and [16]. This reframing is complemented by additional experiments and reanalysis of various published experimental datasets.
https://doi.org/10.7554/eLife.72565.sa1Author response
Essential revisions:
1) The final result (Figure 4) is somewhat disconnected from the majority of the paper that precedes it. Specifically, the authors' procedure that resolves exponential vs. nonexponential growth results in E. coli in alanine being deemed exponential (Figure 2B) only to later be revealed as nonexponential (Figure 4A), albeit weakly. Furthermore, the procedure advertised as distinguishing exponential from linear growth (Figure 3B), when applied to the data, reveals neither (Figure 4). This makes the main point of the paper (the demonstration and resolution of pitfalls) feel disconnected from its application to a particular case, which is more nuanced and likely leaves many questions unanswered.
We try to bridge the gap between Figure 4 and the text preceding it by using simulations of superexponential growth as mentioned in response to Comment 1 of Reviewer 2.
2) The title ("To bin or not to bin…") implies that binning is the main culprit behind potentially misleading analyses, but one of the reviewers argued that in the end, it is linear regression. Each of the two main pitfalls and their resolution would be unchanged if the data were never binned. Binning affects the apparent curvature of the y vs x relationship, but this reads as a more minor point. Therefore, the title may be a bit misleading in service of its poeticism.
We have changed the title of the article to a more apt one – ”Distinguishing different modes of growth using singlecell data”.
3) The authors look at different choices of binning dimensions but do not sufficiently explore the power of their generative model to perform (un)weighted regression or parameter estimation from the notbinned data. They do explore the unbinned data from an analytical statistics approach in sections 5.4.1 and 5.5 but this is not yet extensively explored in the figures and/or discussion.
We apologise for the confusion and we try to clarify it in the manuscript as stated previously in response to Comment 1 of Reviewer 3.
4) This manuscript could be improved by including a discussion on binning singlecell length trajectories, those taken on a single cell within one division cycle, either based on the length or on time. Experimentally, my understanding is that most growth measurements are done on a fixed \Δ t basis, and therefore measurements are averaged using length as a function of time. Another option, one that I have not seen used, is to use a \Δ L basis and measure the time to grow either a specific amount or a specific fractional amount. Of course, these are different ways to slice the same cake, and such an extension may be useful to expand the authors' argument that one needs to utilize a generative model to be able to assess their choice in statistical analysis.
In Section 3 of the manuscript, we discuss binning in relation to the growth rate of singlecell length trajectories as a function of age, and time, and the elongation speed as a function of size. In the simulations conducted in Section 3, we consider the cell lengths measurements to be at equal intervals of time within the cell cycle so as to mimic our mother machine experiments. However, using our methodology of analysing based on a model, it is possible to speculate on the outcome of the data analysis methods in the case of mother machine experiments where measurements are done at equal intervals of length (∆L).
We perform simulations in which lengths within the cell cycle are recorded at ∆L = 0.05_{µ}m intervals and the time from birth corresponding to each length increase is also noted. The results from the simulation are shown in Author response image 1. We find that the binned data trend and the best linear fit for the cells growing exponentially following the adder model deviate from the y=x line for ln $\left(\frac{\text{Ld}}{\text{Lb}}\right)vs\text{}\u27e8\lambda \u27e9{T}_{d}$ plot (Author response image 1A). This is expected from Equations. 2, 27 and 28 of the main text. On interchanging the axes, the best linear fit and the binned data trend are expected to follow the y=x line for exponential growth. This is shown in Author response image 1B for simulations of cells following the adder model and undergoing exponential growth. For cells undergoing linear growth and birth controlling the division, the slope and intercept of the best linear fit of ⟨λ⟩T_{d} vs $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ are analytically found to be approximately 1.04 and 0.03 respectively. This is shown in Author response image 1C for simulations of linear growth and cells following the adder model. The binned growth rate trend in growth rate vs age plots are found to be almost constant for cells undergoing exponential growth irrespective of them following the adder model or the adder per origin model (Author response image 1D). In the case of linear growth, the growth rate is expected to decrease with cell age. This is observed for simulations of linearly growing cells following the adder model (Author response image 1D). We also observe an increase in growth rate in agreement with the simulations of cells following the adder model and undergoing faster than exponential growth (Author response image 1D). All of the results presented here are very similar to that of simulations where the cell lengths were measured at equal intervals in time. We also comment about it in p.37 line 685 of the revised manuscript:
“Note that if the measurements were to be carried out at equal length intervals instead of time, the results discussed in the paper would still remain unchanged.”
5) As the formulation of this approach borrows from and is built heavily on previous data and previous statistical frameworks, it would be helpful to have more discussion on what is new to this particular manuscript. This is always a difficult balance as one needs the context of the previous work to understand the symbology used, but just reiterating what is already in the literature can lead to a dilution of the new/novel/important discoveries of the manuscript.
We try to be clearer in the manuscript, as detailed in the response to Comment 1 of Reviewer 1.
6) The authors should resolve the use of the term "cell size" to mean both "cell length" and "cell volume" (line 107, line 431/432). The authors briefly touch on a related concept of length growth vs biomass growth in 343357. As a first, naively simple model of cell growth, it is possible to come up with a simple model for linear vs volumetric growth by imagining that a cell grows exponentially in length from birth to death. In this simple model, at the appropriate time, it switches from growing a cylinder to growing hemispherically endcaps, as it needs to perfectly duplicate its shape. In this simple model, one can easily see that the instantaneous, relative volumetric and length growth rates are not identical; vary throughout the cell cycle; and even their average over the whole cell cycle is not 1. While the authors mention that the average diameter is quite robust between growth conditions, only in the limit that the cell is very long and not a spherocylinder does the approximation that shape does not matter come in to play. Granted, this naïve model does not incorporate the wealth of information known about division timing or geometric changes during growth, and underscores the authors' main point that a specific model of the dynamical process should be included.
Thank you for raising this point. In the revised version, we point out that cell length is used as a measure of cell size. Thus, the growth rate trend shown in the manuscript is the length growth rate.
However, other cell geometry parameters such as surface area and volume can be used to represent cell size and their respective growth rates might differ from the length growth rate. Theoretically, a comparison between the different rates requires choosing a specific model of cell growth as the reviewers point out. In the revised version of the manuscript, we try to compare the different rates assuming a simple cell morphology and using multiple models of cell growth. We add a discussion of these models as Appendix 1 in the revised manuscript:
“In the paper, we use cell length to represent cell size. However, other cell size characteristics such as cell surface area and cell volume could also be used to denote cell size. […] The width of the cells is assumed to be 0.35 µm. We find that before constriction, the length growth rate increases to a small extent (≈ 6%). However, after constriction there is a rapid increase in length growth rate. The mode of growth in length and volume are not identical.”
In conclusion, length growth rate and volume growth rate do vary and we try to state it more clearly in p.6 line 116 of the revised text:
“Note that we consider length to reflect cell size in this paper rather than other cell geometry characteristics such as surface area and volume. The length growth rate that we elucidate in the paper can be different from the cell volume growth rate as shown in Appendix 1 assuming a simple cell morphology and exponential growth. Using the same cell morphology, we also find the length growth rate to be identical to cell surface growth rate.”
We also reiterate this in Section 5, p.20 line 472 of the revised text:
“To reiterate, the length growth is not the same as cell volume growth as shown in Appendix 1.”
Reviewer #1:
In this manuscript, the authors describe a generative modelbased framework to better analyze stochastic growth data, including bacterial cell growth. They show how this framework can be applied to gain insight into the processes underlying these phenomena. More specifically, they start by showing how binning along different axes in the
This work is wellsupported by simulations and data analysis and will likely be of interest to those trying to understand the processes governing bacterial growth, as well as those studying stochastic growth processes in biology more broadly.
We thank Reviewer 1 for appreciating the methods used in the work and its scope.
Strengths:
– The choice and execution of the simulations were sensible and welldone, respectively, and they provided clarity as to the overall message of the manuscript.
– The conclusions are wellsupported by the data.
– I found the writing to be clear throughout.
Weaknesses:
– It would be good to have a more extensive discussion about what is specifically new here. This is not my particular field, so having a bit more of an introduction about methods beyond binning (if any) that have emerged to understand these data.
Binning and linear regression are the most commonly used methods for probing the correlations between variables obtained using single cell data analysis [Ho et al. (2018), Jun et al. (2018)]. In this paper, we try to address the pitfalls associated with applying these simple procedures. In the revised version, we put more emphasis on linear regression throughout the manuscript for example in the ”Introduction” (p.4 line 54 of the revised text):
“While binning may provide a smooth nonlinear relation between variables, linear regression is used to find a linear relationship between the variables. In addition to binning, we use the ordinary least squares regression where the slope and the intercept of the best linear fit line are obtained by minimizing the squared sum of the difference between the dependent variable raw data and the predicted value. Here, the best fit/the best linear fit is obtained using the raw data and not the binned data. Similar to binning, the assumption underlying linear regression is that our knowledge of xaxis variable is precise while the noise is in the yaxis variable.”
Following this comment, we added the best linear fits to the Figures 2C, 2D and 3A.
We also try to clearly point out the sections which are not novel to the paper. For example, the adder model in Figure 1A has been discussed previously. We emphasize that in p.5 line 77 of the revised text:
“This previously discussed example demonstrates and reiterates the use of statistical analysis on singlecell data to understand the underlying cell regulation mechanisms.”
Similarly, the novel results of the paper such as obtaining the best linear fits for the plots $\mathrm{ln}(\frac{{L}_{d}}{{L}_{b}})$ vs ⟨λ⟩T_{d} and its flipped axes are based on a class of models studied by Eun et al. (2018). We try to make it clearer in p.7 line 133 of the revised manuscript:
“For that purpose, we use a previously studied model [Eun et al. (2018)] which considers growth to be exponential with the growth rate distributed normally and independently between cell cycles with mean growth rate ⟨λ⟩ and standard deviation CV_{λ}⟨λ⟩.”
Reviewer #2:
[…]Weaknesses
The final result (Figure 4) is somewhat disconnected from the majority of the paper that precedes it. Specifically, the authors' procedure that resolves exponential vs. nonexponential growth results in E. coli in alanine being deemed exponential (Figure 2B) only to later be revealed as nonexponential (Figure 4A), albeit weakly. Furthermore, the procedure advertised as distinguishing exponential from linear growth (Figure 3B), when applied to the data, reveals neither (Figure 4). This makes the main point of the paper (the demonstration and resolution of pitfalls) feel disconnected from its application to a particular case, which is more nuanced and likely leaves many questions unanswered.
We attempt to resolve the gap between Figure 4 and the text preceding it by showing that growth rate vs age plot can be used to infer the modes of growth, including, but not limited to, exponential and linear growth. To verify this, we simulate the adder model for cells undergoing superexponential growth. The binned growth rate trend as a function of age for superexponential growth is shown in Figure 3B and the following text is added to p.13 line 270 of the revised manuscript:
“Thus, the two growth modes (exponential and linear) could be differentiated using the growth rate vs age plot (for details see Section 5.7). However, the growth rate vs age plots can be used to infer the mode of growth beyond the two discussed above. We show this by using simulations of cells following the adder model and undergoing faster than exponential or superexponential growth (see the Simulations section for details). In such a case, the growth rate is expected to increase. This increase in growth rate is shown in Figure 3B using simulations. The binned data trend (red triangles) again matches the growth rate mode used in the simulations (red dotted line). Thus, the growth rate vs age plots are a consistent method to distinguish linear from exponential and superexponential growths.”
The details about the simulation of superexponential growth has been added to the Simulations section.
We have also added theoretical predictions (dotted lines of same color) to the growth rate vs age curves of different models shown in Figure 3B of the revised manuscript, that agree well with our simulations.
The title ("To bin or not to bin…") implies that binning is the main culprit behind potentially misleading analysis, but I would argue that in the end, it is linear regression. Each of the two main pitfalls and their resolution would be unchanged if the data were never binned, I believe. Binning affects the apparent curvature of the y vs x relationship, but this reads as a more minor point. Therefore, the title may be a bit misleading in service of its poeticism.
We thank the reviewer for the comment and we have changed the title to a more apt one for the manuscript ”Distinguishing different modes of growth using singlecell data”.
Reviewer #3:
Kar et al., examine an interesting and important question of how to make sense of large sets of observational data, specifically cell length data, which may or may not be consistent with various underlying biological mechanisms. As datasets improve in their technical quality (increasing spatiotemporal resolution, increasing numbers of observations), there is hope that the community will be able to resolve differences between underlying cell biological mechanisms of cell size homeostasis. As the authors point out, these interpretations and analyses require statistical analysis that can accurately perform the model selection or parameter estimation task of interest.
We thank Reviewer 3 for the comments. Indeed, our message in the paper is to use underlying biological models to aid the inference of biological mechanisms using various statistical analyses methods.
1) The authors succeed in bringing attention to the issue of appropriate binning when analyzing large datasets. The authors focus their figures and discussion on an important, and practical issue, as many researchers perform linear regression on binned data. The title and framing of the manuscript imply that it will provide a comparison with statistical methods that do not involve binning. The authors look at different choices of binning dimensions, but do not sufficiently explore the power of their generative model to perform (un)weighted regression or parameter estimation from the notbinned data. They do explore the unbinned data from an analytical statistics approach in section 5.4.1 and 5.5 but this not yet extensively explored in the figures and/or discussion.
We apologize for the confusion. Throughout the paper, we perform linear regression on the raw data and not the binned data. We have now added a clarification statement in p.4 line 58 of the revised text:
“Here, the best fit/the best linear fit is obtained using the raw data and not the binned data.”
Further, we have changed the title of the manuscript to ”Distinguishing different modes of growth using singlecell data”. The paper aims to bring forth the issues in both binning and linear regression which arise from similar sources i.e., the intrinsic noise affecting the xaxis variable and the inspection bias. We discuss these issues in relation to mode of growth in single cells as the new title now states.
In addition to section 5.4.1 and 5.5 where we discuss calculating the best fit line for exponential and linear growth, we also try to explicitly mention linear regression on nonbinned data in the figures, and the Results and Discussion section. This is shown by the addition of the best linear fit/best fit in Figures 2C, 2D and 3A.
https://doi.org/10.7554/eLife.72565.sa2Article and author information
Author details
Funding
USIsrael BSF Research Grant (2017004)
 Jaan Männik
National Institutes of Health (R01GM127413)
 Jaan Männik
National Science Foundation (NSF CAREER 1752024)
 Ariel Amir
National Science Foundation (NSF award 1806818)
 Prathitha Kar
National Institutes of Health (R01 AI143611)
 Ariel Amir
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank Ethan Levien, and Jie Lin for useful discussions, Jane Kondev, Xili Liu, and Marco Cosentino Lagomarsino for their useful feedback on the manuscript, Da Yang and Scott Retterer for help in microfluidic chip making, and Rodrigo ReyesLamothe for a kind gift of strain. Authors acknowledge technical assistance and material support from the Center for Environmental Biotechnology at the University of Tennessee. A part of this research was conducted at the Center for Nanophase Materials Sciences, which is sponsored at Oak Ridge National Laboratory by the Scientific User Facilities Division, Office of Basic Energy Sciences, U.S. Department of Energy. This work has been supported by the USIsrael BSF research grant 2017004 (JM), the National Institutes of Health award under R01GM127413 (JM), NSF CAREER 1752024 (AA), NIH grant R01 AI143611 (AA) and NSF award 1806818 (PK).
Senior Editor
 Aleksandra M Walczak, École Normale Supérieure, France
Reviewing Editor
 Gordon J Berman, Emory University, United States
Reviewer
 Benjamin P Bratton, Princeton University, United States
Publication history
 Received: July 28, 2021
 Accepted: November 21, 2021
 Accepted Manuscript published: December 2, 2021 (version 1)
 Accepted Manuscript updated: December 8, 2021 (version 2)
 Version of Record published: January 4, 2022 (version 3)
 Version of Record updated: January 7, 2022 (version 4)
Copyright
© 2021, Kar et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 2,104
 Page views

 272
 Downloads

 3
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading

 Cancer Biology
 Microbiology and Infectious Disease
A small percentage of bladder cancers in the general population have been found to harbor DNA viruses. In contrast, up to 25% of tumors of solid organ transplant recipients, who are at an increased risk of developing bladder cancer and have overall poorer outcome, harbor BK polyomavirus (BKPyV). To better understand the biology of the tumors and the mechanisms of carcinogenesis from potential oncoviruses, we performed whole genome and transcriptome sequencing on bladder cancer specimens from 43 transplant patients. Nearly half of tumors from this patient population contained viral sequences. The most common were from BKPyV (N=9, 21%), JC polyomavirus (N=7, 16%), carcinogenic human papillomaviruses (N=3, 7%), and torque teno viruses (N=5, 12%). Immunohistochemistry revealed variable Large T antigen expression in BKPyVpositive tumors ranging from 100% positive staining of tumor tissue to less than 1%. In most cases of BKPyVpositive tumors, the viral genome appeared to be clonally integrated into the host chromosome consistent with microhomologymediated end joining and coincided with focal amplifications of the tumor genome similar to other virusmediated cancers. Significant changes in host gene expression consistent with the functions of BKPyV Large T antigen were also observed in these tumors. Lastly, we identified four mutation signatures in our cases with those attributable to APOBEC3 and SBS5 being the most abundant. Mutation signatures associated with the antiviral drug, ganciclovir, and aristolochic acid, a nephrotoxic compound found in some herbal medicines, were also observed. The results suggest multiple pathways to carcinogenesis in solid organ transplant recipients with a large fraction being virusassociated.

 Genetics and Genomics
 Microbiology and Infectious Disease
Plasmids enable the dissemination of antimicrobial resistance (AMR) in common Enterobacterales pathogens, representing a major public health challenge. However, the extent of plasmid sharing and evolution between Enterobacterales causing human infections and other niches remains unclear, including the emergence of resistance plasmids. Dense, unselected sampling is highly relevant to developing our understanding of plasmid epidemiology and designing appropriate interventions to limit the emergence and dissemination of plasmidassociated AMR. We established a geographically and temporally restricted collection of human bloodstream infection (BSI)associated, livestockassociated (cattle, pig, poultry, and sheep faeces, farm soils) and wastewater treatment work (WwTW)associated (influent, effluent, waterways upstream/downstream of effluent outlets) Enterobacterales. Isolates were collected between 20082020 from sites <60km apart in Oxfordshire, UK. Pangenome analysis of plasmid clusters revealed shared 'backbones', with phylogenies suggesting an intertwined ecology where wellconserved plasmid backbones carry diverse accessory functions, including AMR genes. Many plasmid 'backbones' were seen across species and niches, raising the possibility that plasmid movement between these followed by rapid accessory gene change could be relatively common. Overall, the signature of identical plasmid sharing is likely to be a highly transient one, implying that plasmid movement might be occurring at greater rates than previously estimated, raising a challenge for future genomic One Health studies.