## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Noise-driven growth rate gain in clonal cellular populations

Edited by Daniel L. Hartl, Harvard University, Cambridge, MA, and approved February 5, 2016 (received for review October 19, 2015)

## Significance

Differences between individuals exist even in the absence of genetic differences, e.g., in identical twins. Over the last decade, experiments have shown that even genetically identical microbes exhibit large cell-to-cell differences. In particular, the timing of cell division events is highly variable between single bacterial cells. The effect of this variability on long-term growth and survival of bacteria, however, remains elusive. Here, we present a striking finding showing that a bacterial population grows faster on average than its constituent cells. To explain this counterintuitive result, we present a mathematical model that precisely predicts our measurements. Furthermore, we show an empirical growth law that constrains the maximal growth rate of *Escherichia coli*.

## Abstract

Cellular populations in both nature and the laboratory are composed of phenotypically heterogeneous individuals that compete with each other resulting in complex population dynamics. Predicting population growth characteristics based on knowledge of heterogeneous single-cell dynamics remains challenging. By observing groups of cells for hundreds of generations at single-cell resolution, we reveal that growth noise causes clonal populations of *Escherichia coli* to double faster than the mean doubling time of their constituent single cells across a broad set of balanced-growth conditions. We show that the population-level growth rate gain as well as age structures of populations and of cell lineages in competition are predictable. Furthermore, we theoretically reveal that the growth rate gain can be linked with the relative entropy of lineage generation time distributions. Unexpectedly, we find an empirical linear relation between the means and the variances of generation times across conditions, which provides a general constraint on maximal growth rates. Together, these results demonstrate a fundamental benefit of noise for population growth, and identify a growth law that sets a “speed limit” for proliferation.

Cell growth is an important physiological process that underlies the fitness of organisms. In exponentially growing cell populations, proliferation is usually quantified using the bulk population growth rate, which is assumed to represent the average growth rate of single cells within a population. In addition, basic growth laws exist that relate ribosome function and metabolic efficiency, macromolecular composition, and cell size of the culture as a whole to the bulk population growth rate (1⇓–3). Population growth rate is therefore a quantity of primary importance that reports cellular physiological states and fitness.

However, at the single-cell level, growth-related parameters such as the division time interval and division cell size are heterogeneous even in a clonal population growing at a constant rate (4⇓⇓⇓⇓–9). Such “growth noise” causes concurrently living cells to compete within the population for representation among its future descendants. For example, if two sibling cells born from the same mother cell had different division intervals, the faster dividing sibling is likely to have more descendants in the future population compared with its slower dividing sister, despite the fact that progenies of both siblings may proliferate equally well (Fig. 1). Intrapopulation competition complicates single-cell analysis because any growth-correlated quantities measured over the population deviate from intrinsic single-cell properties (10⇓–12). In the case of the toy model described in Fig. 1, cells are assumed to determine their generation times (division interval) randomly by roll of a dice. The mean of intrinsic cellular generation time is thus

A classical study of theoretical and experimental microbiology attempted to reveal the discrepancies between mean cellular generation times and population doubling times in real bacterial populations (13). Experimental methods and techniques available at that time, however, hampered reliable tests. Recently, the techniques of single-cell time-lapse microscopy have advanced to a great extent, revealing the heterogeneous and stochastic nature of single-cell dynamics quantitatively (14). With the aid of microfluidic platforms, tracking single cells over many generations in controlled constant or changing environments has also become feasible, providing insights into the mechanisms of cell size homeostasis and stress responses (3, 6, 7, 9, 15⇓⇓⇓⇓⇓⇓–22).

In this study, through microfluidics time-lapse microscopy and single-cell analysis on large-scale, single-cell lineage trees, we reveal that clonal populations of *Escherichia coli* indeed grow with a doubling time that is smaller than the mean doubling time of their constituent cells under broad, balanced-growth conditions. We show that the observed growth rate gains and population age structures are predictable from cellular generation time distributions based on a simple age-structured population model. Furthermore, we reveal unique features of long single-cell lineages within populations in competition, and provide a history-based formulation that connects growth rate gain with a measure of statistical deviation between isolated and competing lineages. Finally, we demonstrate a linear relation between the means and the variances of generation time across conditions, which constrains the maximum growth rate of this organism.

## Results

### Microfluidics Time-Lapse Microscopy Reveals Long-Term Single-Cell Lineage Tree Structures.

To investigate the role and the consequences of intrapopulation growth noise, we developed a custom microfluidics device for single-cell analysis, which we call the “dynamics cytometer” (Fig. 2 *A–C*, *SI Appendix*, Fig. S1, and Movie S1). The dynamics cytometer is similar to the previously reported polydimethylsiloxane (PDMS)-based microfluidics devices designed for long-term single-cell observation (7, 9, 15, 18⇓⇓–21). The important differences are (*i*) the microchannels are directly created on a glass coverslip, not on PDMS; and (*ii*) the channel region is covered by a semipermeable membrane (*SI Appendix*, Fig. S1) via biotin–streptavidin bonding by chemically decorating the surface of a microfabricated glass coverslip with biotin and membrane with streptavidin (23). The narrow and shallow growth channels can harbor 25 ∼ 40 cells at a time. A small fraction of the cells within growth channels are constantly removed to the wide and deep flow channels, which maintains the number of cells in the growth channels nearly constant. In contrast to the “mother machine” (7), the width of the growth channels is wider than a single-cell width, and both ends are open to flow channels. This configuration introduces competition among cells to remain within the growth channels, which is an important feature of this device for evaluating growth properties of both single cells and the population. Fresh medium is constantly supplied to the cells both from the flows in the flow channels and above the membrane, allowing efficient environmental control. Cellular growth rates in this device were indeed stable during the entire period of observation (*SI Appendix*, Fig. S2). Furthermore, the statistics of cellular growth were indistinguishable across different locations within the growth channels (*SI Appendix*, Fig. S3), which confirms efficient control of environmental conditions around the cells. Using this device, we observed clonal proliferation processes of *E. coli* in constant environments by time-lapse microscopy, obtaining large-scale single-cell lineage tree structures (Fig. 2*D* and *SI Appendix*, Fig. S4) as well as dynamics of cell size along long, single-cell lineages (Fig. 2*E* and *SI Appendix*, Fig. S5).

### Analysis of Single-Cell Lineage Trees Reveals That Cell Populations Grow Faster than Their Constituent Single Cells.

Large-scale single-cell lineage trees obtained by the dynamics cytometer contain sufficient information for determining population growth rates and cellular generation time distributions, to which we applied a careful measurement analysis described in *SI Appendix*. For population growth rate, we first evaluated the instantaneous division rate as *t* is the time-lapse interval; *N*(*t*) is the number of cells in the growth channels at time *t*; and *D*(*t*) is the number of cells that divided between *t* and *t* + Δ*t*. Note that *N*(*t*) + *D*(*t*) is the number of cells at *t* + Δ*t* when no cell is removed from the population. Population growth rate Λ_{p} is then measured as the average of Λ(*t*), i.e., *n* is the number of time points in the measurements, yielding the population doubling time

A cell’s generation time is the age *τ* at which a newborn cell will divide next (here “age” means time elapsed since the previous division), which is a random variable with probability distribution function *g*(*τ*). The generation time distribution is related to the age-specific division rate *b*(*τ*) by *SI Appendix*). Thus, knowing *b*(*τ*) determines *g*(*τ*) as well. We measured the age-specific division rate from the single-cell data using *N*_{a}(*τ*) is the number of cells that reached age *τ*, and *D*_{a}(*τ*) is the number of cells that divided between age *τ* and *τ* + Δ*t*. We confirmed by simulation that these statistical estimators provide precise measures of population growth rate and generation time distributions both with and without cell removal (*SI Appendix*, Figs. S6–S10). Note that all of the information above is obtained directly from lineage tree structures without any fitting to elongation curves. Indeed, our precise measurement of Λ_{p} demonstrates that comparison between population growth rate and simple mean of individual cells’ elongation rate within a population requires attention because they often report significantly different values (*SI Appendix*, Table S1).

Using the validated growth parameters above, we evaluated generation time distributions (Fig. 3*A* and *SI Appendix*, Fig. S11) and population growth rates of four *E. coli* strains listed in *SI Appendix*, Table S2, under several culture conditions, in which temperature and nutrient supplements were altered. Under those conditions, the mean of unbiased cellular generation time ranged from 26 to 170 min. The measured generation time distributions were asymmetric and similar to gamma distributions [=*k* and *θ* are shape and scale parameters], although significant deviations existed especially in the tails (*SI Appendix*, Fig. S11). Comparison showed that the mean generation time *T*_{d} in most cases (Fig. 3*B*), which means that clonal populations grow faster than the constituent single cells on average under a broad range of constant environmental conditions. The relative gain of population growth rate was positively correlated with the noise [coefficient of variation (CV)] of the generation time distributions (Fig. 3*C*), indicating that generation time noise increases population growth rate.

### Growth Rate Gain and Population Age Structures Can Be Predicted from Generation Time Distributions.

To gain further insights into the observed growth rate gain, we compared the experimental results with the prediction of an age-structured population model (Fig. 3*D*) (12, 13, 24) in which cells divide according to an age-dependent division probability. An important prediction of this model is that the generation time distribution *g*(*τ*) and population growth rate Λ_{p} must satisfy the so-called Euler–Lotka equation:*SI Appendix* for derivation). This equation allows us to calculate Λ_{p} using knowledge of *g*(*τ*). We thus calculated Λ_{p}^{(theory)} [=ln(2)/*T*_{d}^{(theory)}] from experimentally obtained *g*(*τ*) according to the computational method described in *SI Appendix*, finding good agreement with the directly measured value (Fig. 3*E*). When generation time follows a gamma distribution, i.e., *k* and *θ* are shape and scale parameters, Eq. **1** is analytically solvable, and growth rate gain is determined solely by the CV. Despite the fact that the generation time distributions were not perfectly gamma, the result in Fig. 3*C* agrees in trend with this relation.

In addition to population growth rate, this simple model predicts the steady-state age distribution of the population as follows:*SI Appendix*). This equation shows that the generation time distribution of single cells determines the age distribution of the population because Λ_{p} is itself determined by *g*(*τ*) using [**1**]. We compared experimentally measured age distributions of the population with those predicted from *g*(*τ*) using [**2**] according to the method described in *SI Appendix*, finding excellent agreement in all of the conditions without any adjustable parameters (Fig. 4).

### Age Distributions Along Long Single-Cell Lineages in Populations Exhibit Optimal Lineage Statistics.

In addition to population growth rates and age distributions, the model makes testable predictions regarding the histories of individual cells within a population. Theoretical analysis of the cell proliferation process from the viewpoint of cellular lineages shows that the distribution **1**], is in fact the typical generation time distribution along a lineage extracted by randomly choosing a descendant cell and tracking back to its ancestor over many generations (*SI Appendix*). We previously showed that the distribution *g**(*τ*) corresponds to a set of lineages that optimize the population growth rate (12); hence we call such lineages “optimal lineages.” Importantly, optimal lineages possess a generation time distribution that is distinct from that of single cells in isolation. Focusing on the experimental datasets in which we were able to obtain lineages longer than 100 generations, we confirmed that their age distributions correspond to that of optimal lineages, rather than that of single cells in isolation (Fig. 5).

Analyzing the cell proliferation processes from the viewpoint of cell lineages offers an illuminating mathematical interpretation of these differences, as we prove within the context of the simple age-structured model that*SI Appendix*). Noting that

### The Means and the Variances of Generation Time Across Conditions Are Linearly Related.

The relations we examined thus far do not impose any constraints on generation time distributions. However, we noticed that the means and variances of the three *E. coli* strains derived from W3110 (F3 series) were approximately linearly related (Fig. 6). In general, such a relation suggests that the dominant contribution to stochasticity in the cell division time may arise from a number of discrete steps within the cell division cycle, each contributing comparable and largely independent portions of the variance in timing. If the number of steps increases as growth conditions deteriorate, a linear mean–variance scaling relation would arise. It is noteworthy that the *X* intercept of the linear scaling relation is positive (24.4 ± 0.8 min) and close to the minimum generation time of *E. coli* in rich media (∼20 min). Thus, this relation may be relevant for establishing the upper bound of growth rate of these strains in any constant environments. Because the *X* intercept is positive, the CV of generation time is expected to take the maximum value at an intermediate position, where the mean generation time is twice the minimum (*X*-intercept value). This is indeed seen in Fig. 3*C* (*Inset*). The variances of BrF2 strain were significantly smaller than this trend, suggesting that the intrinsic level of generation time heterogeneity is strain specific.

## Discussion

The growth of cells constitutes an important physiological process that contributes critically to the fitness of an organism, while also imposing global constraints on intracellular gene expression levels (29), metabolic efficiency (1), and macromolecular composition (2). Most studies on cellular growth have referred to bulk population growth rates in measurements, but recent research on the heterogeneous and stochastic nature of growth at the single-cell level poses new, fundamental questions regarding (*i*) how precisely the insights gained through population measurements apply to single-cell properties (3), and (*ii*) how single cells as a whole bring about the observed population properties (12, 13). Our study tackled the latter issue both experimentally and theoretically, showing that a simple age-structured population model explains the intricate relationships between the growth rates of the cells and the population.

The result that the growth rate gain increases with the CV of generation time (Fig. 3*C*) clearly demonstrates the importance of growth noise in determining population growth rate. This suggests that there are two routes for an organism to evolve to a higher population growth rate in a constant environment: one is simply to decrease the mean generation time; and the other is to increase the heterogeneity of generation times in a population. In reality, a single mutation may affect both mean and variance of generation times simultaneously, resulting in a synergetic or antagonistic effect on the population growth rate. Discriminative measurements on these two contributions have never been performed to our knowledge. In experimental evolution, an increase of population growth rate has usually been attributed exclusively to the change of mean growth properties; measuring single-cell growth heterogeneity over the course of an evolutionary process may reveal the dominant contribution to an increase of population growth rate at each step of the evolutionary trajectories.

The measurements in this study are restricted to the stress-free conditions where cells can grow and divide stably for hundreds of generations. Generally, cells placed in stressful environments show increased levels of phenotypic heterogeneity (16, 30). One apparent role of large phenotypic heterogeneity is to increase the chance of producing fitter individuals to be selectively amplified within a population. Additionally, as our results suggest, heterogeneity increases the population growth rate achievable for a fixed, mean generation time, a result that may clarify the role of stress-induced phenotypic heterogeneity observed in previous experiments.

It is remarkable that a simple population model, using only an age-dependent cell division process, can explain the relationships between the growth rates of cells and the population across such broad culture conditions. We attribute this to having focused only on the age-related growth parameters and cell division statistics, while not modeling cell-size related parameters. Indeed, cell size at birth and generation time are weakly correlated negatively (Pearson correlation coefficient, −0.48 ∼ −0.22; *SI Appendix*, Fig. S12), and this effect must be considered to account for cell size stability (3, 31⇓⇓⇓–35). Nevertheless, our results suggest that, as far as age-related parameters and population growth rates are concerned, cell size information does not play a predominant role in *E. coli* over a broad set of culture conditions. It is also important to note that this model assumes no correlation of generation time between mother and daughter cells. Generally, positive mother–daughter correlation increases population growth rate, and negative correlation decreases it (13). However, our experimental results show that mother–daughter correlations of generation time are weak: a positive correlation of 0.2 would increase population growth rate by ∼2% (13), and most of the observed correlations are below this level (*SI Appendix*, Table S5).

With regard to cell size stability, our study highlights the importance of lineage differences. As recognized by previous studies (6, 7), a small fraction of *E. coli* cells in clonal populations often become filamentous even in constant environments. We also observed such filamentous cells in our measurements across the population (*SI Appendix*, Fig. S5 *A* and *B*), but importantly these cells are rarely observed along the optimal lineages (Fig. 2*E* and *SI Appendix*, Fig. S5). This is likely due to the tail of *g**(*τ*), which decays faster than *g*(*τ*) by the factor

The linear relation between the means and the variances of generation time in Fig. 6 suggests that we can predict the noise levels of generation time only from the means without knowing the details about the environments. A similar idea has been proposed in several recent studies. For example, Iyer-Biswas et al. (31, 36) proposed that the distributions of generation time normalized by their means collapses onto the universal curve based on their experiments with *Caulobacter crescentus* cultured under various temperature conditions, and that an autocatalytic cycle model for cell cycle control might account for this scaling property. Similarly, Taheri-Araghi et al. (3) proposed the scaling of generation time distribution based on their single-cell analysis on *E. coli* growth under different nutrient conditions at a fixed temperature, and presented the adder model, which assumes that cells achieve cell size homeostasis by adding a constant size between birth and division irrespective of birth size. For this simple scaling rule to be valid, the CVs of generation time must be constant irrespective of the environmental conditions, and the SDs become proportional to the means. However, our results in Fig. 3*C* (*Inset*) and Fig. 6 do not follow this rule. Thus, our results suggest that the scaling rule breaks down in some culture conditions and/or organisms not explored in the previous studies. It should be noted that the mean–variance linear relation is tested only for the W3110-derived strains in our study; characterizing which rule is obeyed by different strains and organisms might be a fruitful subject for future evolutionary study.

In summary, our measurements have established a predictive framework that bridges from single-cell properties to population-level dynamics in conditions where cells can grow stably. The same approach can be used to discover rules that govern growth in severe stress environments, including antibiotics and chemotherapeutics, and to investigate the generality and prevalence of single-cell and population growth laws across strains, species, cell types, and environments.

## Materials and Methods

### Cell Strains and Culture Conditions.

We used W3110 and B/r derivative *E. coli* strains in this study (*SI Appendix*, Table S2). Cells were grown in M9 minimal medium. The culture conditions were changed by different nutrient supplements and temperature as listed in *SI Appendix*, Table S3. See *SI Appendix* for details.

### Single-Cell Time-Lapse Observation.

We implemented time-lapse observation of *E. coli* proliferating in the dynamics cytometer using Nikon Ti-E microscope controlled by micromanager (https://micro-manager.org/). The time-lapse images were analyzed by a custom macro of ImageJ (imagej.nih.gov/ij/), and the exported data were further analyzed using a custom C-program. See *SI Appendix* for details.

## Acknowledgments

We thank Ippei Inoue (Ajinomoto Company) for technical advice, Atsushi Miyawaki (RIKEN) for providing Venus/pCS2 plasmid, Hironori Niki (National Institute of Genetics) for providing pKP2375 plasmid, CGSC at Yale for providing *E. coli* B/r strain, and Kenji Yasuda (Tokyo Medical and Dental University) for the microfabrication facilities. This work was supported by Japan Society for the Promotion of Science KAKENHI Grants 25711008 (to Y.W.), 15H05746 (to K.K. and Y.W.), 13J09314 (to M.H.), and 14J01376 (to T.N.); Japan Science and Technology Agency Precursory Research for Embryonic Science and Technology Program (Y.W.); NIH Grant R01-GM-097356 (to E.K.); and Platform for Dynamic Approaches to Living System from Ministry of Education, Culture, Sports, Science and Technology, Japan and Japan Agency for Medical Research and Development, (K.K. and Y.W.).

## Footnotes

- ↵
^{1}To whom correspondence should be addressed. Email: cwaka{at}mail.ecc.u-tokyo.ac.jp.

Author contributions: M.H. and Y.W. designed research; M.H., T.N., H.N., R.O., K.K., E.K., and Y.W. performed research; Y.W. contributed new reagents/analytic tools; M.H., R.O., S.A., and Y.W. analyzed data; and M.H., E.K., and Y.W. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1519412113/-/DCSupplemental.

## References

- ↵.
- Scott M,
- Gunderson CW,
- Mateescu EM,
- Zhang Z,
- Hwa T

- ↵
- ↵
- ↵.
- Powell EO

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Wakamoto Y, et al.

- ↵.
- Balaban NQ,
- Merrin J,
- Chait R,
- Kowalik L,
- Leibler S

- ↵.
- Aldridge BB, et al.

- ↵
- ↵
- ↵
- ↵.
- Lambert G,
- Kussell E

- ↵
- ↵.
- Bellman R,
- Harris TE

- ↵.
- Leibler S,
- Kussell E

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Iyer-Biswas S, et al.

- ↵
- ↵.
- Osella M,
- Nugent E,
- Cosentino Lagomarsino M

- ↵
- ↵
- ↵

## Citation Manager Formats

### More Articles of This Classification

### Biological Sciences

### Biophysics and Computational Biology

### Related Content

- No related articles found.