Skip to main content
Advertisement
  • Loading metrics

Characterizing non-exponential growth and bimodal cell size distributions in fission yeast: An analytical approach

Abstract

Unlike many single-celled organisms, the growth of fission yeast cells within a cell cycle is not exponential. It is rather characterized by three distinct phases (elongation, septation, and reshaping), each with a different growth rate. Experiments also showed that the distribution of cell size in a lineage can be bimodal, unlike the unimodal distributions measured for the bacterium Escherichia coli. Here we construct a detailed stochastic model of cell size dynamics in fission yeast. The theory leads to analytic expressions for the cell size and the birth size distributions, and explains the origin of bimodality seen in experiments. In particular, our theory shows that the left peak in the bimodal distribution is associated with cells in the elongation phase, while the right peak is due to cells in the septation and reshaping phases. We show that the size control strategy, the variability in the added size during a cell cycle, and the fraction of time spent in each of the three cell growth phases have a strong bearing on the shape of the cell size distribution. Furthermore, we infer all the parameters of our model by matching the theoretical cell size and birth size distributions to those from experimental single-cell time-course data for seven different growth conditions. Our method provides a much more accurate means of determining the size control strategy (timer, adder or sizer) than the standard method based on the slope of the best linear fit between the birth and division sizes. We also show that the variability in added size and the strength of size control in fission yeast depend weakly on the temperature but strongly on the culture medium. More importantly, we find that stronger size homeostasis and larger added size variability are required for fission yeast to adapt to unfavorable environmental conditions.

Author summary

Advances in microscopy enable us to follow single cells over long timescales from which we can understand how their size varies with time and the nature of innate strategies developed to control cell size. These data show that in many cell types, growth is exponential and the distribution of cell size has one peak, namely there is a single characteristic cell size. However data for fission yeast show remarkable differences: growth is non-exponential and the distribution of cell sizes has two peaks, corresponding to different growth phases. Here we construct a detailed stochastic mathematical model of this organism; by solving the model analytically, we show that it is able to predict the two peaked distributions of cell size seen in data and provide an explanation for each peak in terms of various growth phases of the single-celled organism. Furthermore, by fitting the model to the data, we infer values for the rates of all microscopic processes in our model. This method is shown to provide a much more reliable inference than current methods and shed light on how the strategy used by fission yeast cells to control their size varies with external conditions.

Introduction

The fission yeast Schizosaccharomyces pombe is a single-cell eukaryote whose shape is well approximated by a cylinder with hemispherical ends [13]. The length of the rod-shaped cell increases during the G2 phase of the cell cycle, while its width (diameter) remains almost constant. In experiments, length, area, and volume have all been used to characterize cell size. The recent advent of microfluidic techniques allows the tracking of thousands of individual cells over hundreds of cell cycles which potentially enables a detailed investigation of cell growth and size control strategies [4].

It has been reported that cell size grows exponentially in many cell types such as bacteria, cyanobacteria, archaea, budding yeast, and mammalian cells [517]. However, fission yeast undergoes a complex non-exponential growth pattern in each cell cycle, as illustrated by the time-course data of cell size along a typical cell lineage (Fig 1A). At the beginning of the cell cycle, the rod-shaped cell starts to grow by extension at its old cell end (the end that existed before the last division). Later in mid G2 phase, the cell exhibits a transition in cell polarization, and growth is also initiated at the new cell end (the end created during the last division), in a process called new end take-off (NETO) [1, 3]. Cell length increases during the first ∼75% of the cell cycle. Cell elongation stops during the remaining ∼25% of the cell cycle, when mitosis and cytokinesis occur, and the cell subsequently divides into two almost identical progenies [1, 3].

thumbnail
Fig 1. Cell size dynamics in fission yeast.

A: Single-cell time-course data of cell size along a typical cell lineage cultured in the yeast extract medium at 34°C. Here the size of a cell is characterized by its length. The data shown are published in [4]. The green dots show cell sizes at birth. B: Histogram of cell sizes along all cell lineages. The cell size distribution of lineage measurements has a bimodal shape. C: Scatter plot of the birth size versus the division size and the associated regression line. When plotting B, C, we use the data of all 1500 cell lineages cultured in the yeast extract medium at 34°C, each of which is recorded every 3 minutes and is typically composed of 50 − 70 generations. The generation time is 114 ± 23 minutes.

https://doi.org/10.1371/journal.pcbi.1009793.g001

There has been a long-standing controversy about the growth pattern of fission yeast before mitosis [1820]. In earlier studies, exponential [19, 2125], linear [26, 27], and bilinear [1, 3, 2830] growth models have been proposed. The bilinear model consists of two linear growth regimes with different growth rates separated by a rate change point at ∼34% of the cell cycle in wild-type cells, which coincides with NETO [1]. While the predominant viewpoint is that the growth before mitosis is bilinear, more recent data has confirmed exponential growth of mass with some changes in density through the cell cycle [25]. In practice, however, it is very difficult to distinguish the exponential and bilinear growth patterns due to the stochasticity of growth dynamics and the relatively low temporal and spatial resolution of the data.

A remarkable feature of lineage measurements is the bimodal shape of the cell size distribution computed over many cell cycles (Fig 1B). Note that this was not reported in the original paper [4] but stems from an analysis of their published data. Bimodal lineage distributions of fission yeast have not been previously reported in the literature, possibly because such high throughput data have become available only recently. As well, previous studies focused more on the distributions of birth and division sizes [3133], instead of the distribution of cell size over the whole cell lineage. Actually, the latter contains much more information than the former, since it reflects the full cell cycle dynamics. Recent studies have shown that if cell size grows exponentially in each generation, then the distribution of cell size must be unimodal [34, 35]. The main aim of the present paper is to propose a detailed model of cell size dynamics in fission yeast that can characterize its non-exponential growth, cell division, and size homeostasis, as well as develop an analytical theory that can account for the bimodal shape of the cell size distribution.

In the study of cell size dynamics, a core issue is to understand the size homeostasis strategies in various cell types, especially in fission yeast [3651]. There are three popular phenomenological models of cell size control leading to size homeostasis [52]: (i) the timer strategy which implies a constant time between successive divisions regardless of initial size; (ii) the sizer strategy which implies cell division upon attainment of a critical size, and (iii) the adder strategy which implies a constant size addition between consecutive generations. The adder or near-adder behavior has been observed in bacteria, budding yeast, and mammalian cells [10, 11, 16]. However, fission yeast exhibits a sizer-like behavior [36, 53], where cells grow during interphase to a target size of ∼14μm in length before entering mitosis and dividing medially, with the standard deviation of the division size being only 7% of the mean [54, 55]. Recently, significant progress has been made to decipher the molecular mechanism responsible for size control in fission yeast and some key proteins have been found to sense cell size and promote mitotic entry [5458].

A conventional method of inferring the size control strategy is to use the information of cell sizes at birth and at division [36, 59]. This method assumes that the birth size Vb and the division size Vd in each generation are related linearly by (1) where 0 ≤ β ≤ 2 and γ ≥ 0 are two constants and ϵ is a noise term independent of the birth size. Here β characterizes the strength of size control with β = 0, β = 1, and β = 2 corresponding to the sizer, adder, and timer strategies, respectively. Using the data of birth and division sizes across generations, the parameter β can be determined as the slope of the regression line of the division size on the birth size. The linear relationship with a slope less than 1 between the birth and division sizes in Fig 1C is suggestive of a sizer-like mechanism. However it is also clear that the relationship is very weak with numerous outliers and an exceptionally low R2 around 0.1. This makes the inference of the parameter β highly unreliable. Hence another aim of the present paper is to develop a more reliable technique that can be used to accurately infer the size control strategy in fission yeast using a stochastic dynamic approach.

Results

Model specification

Here we consider a detailed model of cell size dynamics in fission yeast across many generations, including a complex three-stage growth pattern, asymmetric and stochastic cell division, and size homeostasis (see Fig 2B for an illustration). In our model, cell size can be interpreted as either cell length, area, or volume. Because the cell width is approximately a constant and the cell is rod-like in shape, these three size measures are proportional to each other. The model is based on a number of assumptions that are closely related to lineage data obtained using microfluidic devices. The assumptions are as follows and the specific meaning of all model parameters is listed in Table 1.

thumbnail
Fig 2. A detailed model of cell size dynamics in fission yeast.

A: Three-stage growth pattern of fission yeast: an elongation (G2) phase where cell size grows exponentially with rate g0, followed by a septation phase during which the septum is formed and cell size remains constant, and then followed by a reshaping phase where cell length increases abruptly with a higher exponential growth rate g1 > g0. Here Vb is the size at birth, Vs is the size at septation, and Vd is the size at division. B: Schematic illustrating a detailed model of cell size dynamics describing three cell growth phases, size homeostasis, and symmetric or asymmetric partitioning at cell division (see inset graph). Each cell can exist in N effective cell cycle stages. Cell elongation occurs during the first N0 stages, septation occurs during the intermediate NN0N1 stages, and reshaping of the new cell end occurs during the last N1 stages. The transition rate from one stage to the next at time t is proportional to the αth power of the cell size V(t) with α > 0 being the strength of size control and a > 0 being the proportionality constant. This guarantees that larger cells at birth divide faster than smaller ones to achieve size homeostasis. At stage N, a mother cell divides into two daughters that are typically different in size via asymmetric cell division. Symmetric division is the special case where daughters are equisized.

https://doi.org/10.1371/journal.pcbi.1009793.g002

1) The growth of cell size in fission yeast is very different from the exponential growth observed in many other cell types [6]. Actually, fission yeast undergoes a non-exponential three-stage growth pattern: an elongation phase followed by a septation phase and a reshaping phase (Fig 2A) [26]. During the elongation (G2) phase, we assume that the size of each cell grows exponentially with rate g0, which is supported by experiments [19, 2125]. In some previous papers, the growth before mitosis is assumed to be bilinear with a change in the slope in the mid G2 phase [1, 3, 2830]. However, bilinear growth is very close to exponential growth and it is very difficult to distinguish them experimentally due to the stochasticity of growth dynamics and the relatively low temporal and spatial resolution of the data.

After the elongation phase, the growth ceases for a period during which the septum is formed [26, 54]. From the lineage data in Fig 1A, it seems also reasonable to assume that there is a small non-zero growth rate in the septation phase. However, according to the principle of parsimony, we choose not to introduce an extra parameter and assume zero growth rate during septation.

At the end of the septation phase, there is a sharp increase in cell length for a short period (a few minutes), during which the new cell end pops out and forms a hemisphere due to turgor pressure [60]. This period is referred to as the reshaping phase since it corresponds to the rounding off of the new end from the septum. Hence intuitively, the mean added size in this phase should roughly equal the size of the two hemispherical end caps. During the reshaping phase, we assume that cell size grows exponentially with a higher rate g1 > g0. In general, the duration of this phase is very variable and thus this part is often omitted from cell growth studies. Here we choose to include this part into our modelling since it naturally appears in the lineage data (Fig 1A). Since the period of the reshaping phase is very short, the specific growth pattern in this phase does not have much effect on the overall dynamics. The choice of exponential growth during this phase is convenient for an analytical treatment.

2) Each cell can exist in N effective cell cycle stages, denoted by 1, 2, …, N. Note that the effective cell cycle stages introduced here do not directly correspond to the four biological cell cycle phases (G1, S, G2, and M) or the three growth phases (elongation, septation, and reshaping). Rather, a cell cycle phase or a growth phase corresponds to multiple effective cell cycle stages (Fig 2). Similar assumptions have been successfully used to reproduce the measured variability in cell cycle phase durations in other cell types [61]. We assume that the cell stays in the elongation phase in the first N0 stages, in the reshaping phase in the last N1 stages, and in the septation phase in the intermediate NN0N1 stages (Fig 2). The transition rate from one stage to the next at a particular time is proportional to the αth power of cell size at that time [35, 62]. In other words, the transition rate between stages at time t is equal to aV(t)α, where V(t) is the cell size at that time, α > 0 is the strength of size control, and a > 0 is a proportionality constant. Under this assumption, larger cells at birth, on average, have shorter cell cycle duration and lesser volume change than smaller ones; in this way size homeostasis is achieved. Interestingly, under symmetric division and small cell size variability, the size control strength α in our model and the size control strength β in the conventional model given in Eq (1) are related by β = 21−α (see Section A in S1 Appendix for the proof).

The entry to mitosis is controlled by a complex gene regulatory network. The cyclin dependent kinase Cdk1, the central mitotic regulator, is regulated by many proteins such as the peripheral membrane kinase Cdr2, the kinase Wee1, the phosphatase Cdc25, and the cyclin B Cdc13 [54]. The prevailing consensus is that the accumulation of regulators upstream of Cdk1, such as Cdr2, Cdc25, and Cdc13, to a critical threshold is required to trigger mitotic entry and cell division, a strategy known as the activator accumulation mechanism [5458]. Biophysically, the N effective cell cycle stages in our model can be understood as different levels of the key protein that triggers cell division. The power law form for the rate of cell cycle progression may come from cooperation of the key protein, as explained in detail in [35, 62]. This power law not only coincides with certain biophysical mechanisms, but also results in a natural scaling transformation among the timer, adder, and sizer, as will be explained later. Besides, we point out that another strategy called the inhibitor dilution mechanism may also be used for size control. This strategy has been observed in budding yeast [54, 63, 64], where the concentration of Whi5 decreases during the G1 phase due to dilution caused by an increase in cell size. When cell volume is sufficiently large, the Whi5 concentration drops below a given threshold to trigger the G1/S transition allowing subsequent DNA replication and budding.

Let Vb and Vd denote the cell sizes at birth and at division in a particular generation, respectively, and let Vs denote the cell size in the septation phase, which is assumed to be a constant. Then the increment in the αth power of cell size, which is referred to as generalized added size, in the elongation phase, , has an Erlang distribution with shape parameter N0 and mean M0 = N0 g0 α/a (see Section A in S1 Appendix for the proof). Similarly, the generalized added size in the reshaping phase, , also has an Erlang distribution with shape parameter N1 and mean M1 = N1 g1 α/a. Therefore, the total generalized added size across the cell cycle, , is the sum of two independent Erlang distributed random variables and has a hypoexponential distribution (also called generalized Erlang distribution) whose Laplace transform is given by (2)

Note that as N → ∞. This means that the generalized added size Δ = M0 + M1 becomes deterministic when N is very large. In other words, a higher threshold for the division protein (the key regulator triggering cell division) level results in a less noisy added size and thus a less noisy cell cycle duration. When N is small, the variability in Δ is much larger. Hence, our model allows the investigation of the influence of added size variability on cell size dynamics. Since the strength of size control and the variability in added size may strongly depend on the culturing condition (strain, medium, temperature) applied, the specific values of α and N may vary a lot in different culturing conditions.

Three special cases deserve special attention. When α → 0, the transition rate between stages is a constant and thus the doubling time has an Erlang distribution that is independent of the birth size; this corresponds to the timer strategy. When α = 1, the added size VdVb has an hypoexponential distribution that is independent of the birth size; this corresponds to the adder strategy. When α → ∞, the αth power of the division size, , has a hypoexponential distribution that is independent of the birth size; this corresponds to the sizer strategy. Intermediate strategies are naturally obtained for intermediate values of α; timer-like control is obtained when 0 < α < 1 and sizer-like control is obtained when 1 < α < ∞ [62].

3) Cell division occurs when the cell transitions from stage N to the next stage 1. At division, most previous papers assume that the mother cell divides into two daughters that are exactly the same in size via symmetric partitioning [31, 6567]. Experimentally, fission yeast in general do not divide perfectly in half. There has been some evidence indicating that there is a small asymmetry in the position of the septum that is slightly nearer the new end [68, 69]. Here we follow the methodology that we devised in [35, 70] and extend previous models by considering asymmetric partitioning at division: the mother cell divides into two daughters with different sizes.

In this paper, division is assumed to occur after the reshaping phase, i.e. after the rounding off of the two new ends. This is consistent with the lineage data in [4, 26] (see also Fig 1A), as well as Fig 2A in [54]. However, in previous papers [1, 3, 36], division is often assumed to occur after the septation phase—when the septum starts to be degraded, the mother cell has divided into two progenies. Under the latter definition, the division size should be the septation size Vs and the birth size should be the size of a progeny from the septum to the old end. However, since our model is based on the lineage data shown in Fig 1A, where the size of a daughter cell from the septum to the old end is not explicitly given, we choose not to use this definition.

If the partitioning of cell size is symmetric, we track one of the two daughters randomly after division [71, 72]; if the partitioning is asymmetric, we either track the smaller or the larger daughter after division [73, 74]. Let Vd and denote the cell sizes at division and just after division, respectively. If the partitioning is deterministic, then we have , where 0 < p < 1 is a constant with p = 0.5 corresponding to symmetric division, p < 0.5 corresponding to smaller daughter tracking, and p > 0.5 corresponding to larger daughter tracking. The value of p can be inferred from experiments. However, in fission yeast, the partitioning of cell size is appreciably stochastic. In this case, we assume that the partition ratio has a beta distribution with mean p [75], whose probability density function is given by (3) where B is the beta function, q = 1 − p, and ν > 0 is referred to as the sample size parameter. When ν → ∞, the variance of the beta distribution tends to zero and thus stochastic partitioning reduces to deterministic partitioning, i.e. f(r) = δ(rp).

We next describe our stochastic model of cell size dynamics. The microstate of the cell can be represented by an ordered pair (k, x), where k is the effective cell cycle stage which is a discrete variable and x is the cell size which is a continuous variable. Note that the cell undergoes deterministic growth in each stage (exponential growth in the first N0 and the last N1 stages and no growth in the remaining NN0N1 stages), and the system can hop between successive stages stochastically. Let pk(x) denote the probability density function of cell size when the cell is in stage k. Then the evolution of cell size dynamics in fission yeast can be described by a piecewise deterministic Markov process whose master equation is given by (4) where h(r) is the function given in Eq (3). In the first, second, and fourth equations, the first term on the right-hand side describe cell growth and the remaining two terms describe transitions between cell cycle stages. In the third equation, the two terms on the right-hand side describe cell cycle stage transitions. In the first equation, the middle term on the right-hand side describes the partitioning of cell size at division.

Analytical distribution of cell size for lineage measurements

Let denote the probability density function of cell size V (here we use V to represent a random variable and use x to represent a realization of V). In our model, we assume that the rate of cell cycle progression has a power law dependence on cell size. This assumption implies an important scaling property of our model: if the dynamics for cell size V has a control strength α (with α < 1 corresponding to timer-like and α > 1 corresponding to sizer-like strategies), then the dynamics for the αth power of cell size, Vα, has an adder strategy. This scaling property serves as the key to our analytical theory.

Recall that the probability distribution of any random variable with nonnegative values is fully determined by its Laplace transform. To obtain the analytical distribution of cell size along a cell lineage, we introduce , which is nothing but the Laplace transform for the αth power of cell size. For simplicity, we first focus on the case of deterministic partitioning. Despite the biological complexity described by our model, the Laplace transform can still be solved exactly in steady-state conditions as (see Section B in S1 Appendix for the proof) (5) where b(λ) is the function given in Eq (2), (6) is another function with A0 = M0/N0 and A1 = M1/N1, and is a normalization constant. From the definition of f(λ) in Eq (6), it is clear that f(λ) tends to infinity as λ → ∞. However, from the definition of b(λ) in Eq (2), the infinite product decays to zero as λ → ∞ at a faster exponential speed. Hence the integral in Eq (5) is always well defined.

In principle, taking the inverse Laplace transform gives the probability density function of Vα, from which the distribution of cell size V can be obtained. Next we introduce how to compute the cell size distribution more effectively using our analytical results. Taking the derivative with respect to λ on both sides of Eq (5), using the change of variables formula, and finally replacing λ by yield (see Section B in S1 Appendix for the proof) (7) where is the probability density function of Vα. This shows that the Fourier transform of is exactly G(λ). Since the Fourier transform and inverse Fourier transform are inverses of each other, we only need to take the inverse Fourier transform of G(λ) so that we can obtain . Finally, the cell size distribution p(x) can be recovered from as (8)

In general, the cell size distribution along a cell lineage can also be numerically computed by carrying out stochastic simulations of the piecewise deterministic Markovian model. However, under the complex three-stage growth pattern of fission yeast, according to our simulations, over 107 stochastic trajectories must be generated in order to obtain an accurate computation of the size distribution (S1 Fig)—this turns out to be very slow. The analytical solution is thus important since it allows a fast exploration of large swathes of parameter space without performing stochastic simulations.

To gain deeper insights into the cell size distribution, we next consider two important special cases. For the case of exponential growth of cell size, there is only the elongation phase and the remaining two phases vanish. In this case, we have N1 = 0 and N = N0; the cell size distribution is still determined by Eq (5) with the functions b(λ) and f(λ) being simplified greatly as (9)

In fact, the analytical cell size distribution for exponentially growing cells has been studied previously in [35], where the distribution of the logarithmic cell size, instead of the original cell size, is obtained. We emphasize that the analytical expression given here is not only much simpler, but also numerically more accurate than the one given in that paper, which includes the integral of an infinite product term which is difficult to compute accurately.

The second case occurs when N → ∞, while keeping r0 = N0/N and r1 = N1/N as constant, where r0 and r1 represent the proportions of cell cycle stages in the elongation and reshaping phases, respectively. In this case, the generalized added size Δ becomes deterministic and the system does not involve any stochasticity. As N → ∞, the Laplace transform given in Eq (5) can be simplified to a large extent as (see Section B in S1 Appendix for the proof) (10) where is the exponential integral, are the birth size, septation size, and division size, respectively, and K = (T0 + Ts + T1)−1 is a normalization constant with being the durations of the three phases, respectively. Note that in [26], the septation size Vs (size after septation) is called the division size and the division size Vd (size after reshaping) is called the fission size. Here the terminology is slightly different. Taking the inverse Laplace transform finally gives the cell size distribution: (11) where IA(x) is the indicator function which takes the value of 1 when xA and takes the value of 0 otherwise, δ(x) is Dirac’s delta function, and (12) are the proportions of subpopulations in the three phases, respectively. This indicates that when added size variability is small, cell size has a distribution that is concentrated on a finite interval between vb and vd.

To validate our theory, we compare the analytical cell size distribution with the one obtained from stochastic simulations under different choices of N (Fig 3A). Clearly, they coincide perfectly with each other. It can be seen that as added size variability become smaller (N increases), the analytical distribution given in Eq (8) converges to the limit distribution given in Eq (11). When N is small, the size distribution is unimodal. As N increases, the size distribution becomes bimodal with the right peak becoming higher and narrower. The bimodality of the size distribution can be attributed to cells in different phases: the left peak corresponds to cells in the elongation phase and the right peak corresponds to cells in the septation and reshaping phases. When N is very large, the size distribution is the superposition of three terms, corresponding to the three phases of cell growth (see the rightmost panel of Fig 3A).

thumbnail
Fig 3. Influence of model parameters on the cell size distribution.

A: Cell size distribution as N varies. The blue curve shows the analytical distribution obtained by taking the inverse Laplace transform of Eq (5) (e.g. using the technique described by Eqs (7) and (8)) and the red circles show the distribution obtained from stochastic simulations. The parameters are chosen as r0 = 0.6, r1 = 0.1, g1 = 2g0, α = 2. B: Cell size distribution as α varies. The parameters are chosen as N = 30, r0 = 0.6, r1 = 0.1, g1 = 2g0. C: Cell size distribution as r0 varies. The parameters are chosen as N = 30, r1 = 0.1, g1 = 2g0, α = 2. D: Cell size distribution as r1 varies. The parameters are chosen as N = 30, r0 = 0.6, g1 = 2g0, α = 2. E: Cell size distribution as g1/g0 varies. The parameters are chosen as N = 30, r0 = 0.6, r1 = 0.1, α = 2. In A-E, the parameters g0 and p are chosen as g0 = 0.01, p = 0.5 and the parameters a, M0, M1 are chosen so that the mean cell size 〈V〉 = 3.

https://doi.org/10.1371/journal.pcbi.1009793.g003

To gain a deeper insight, we illustrate the cell size distribution as a function of the parameters α, r0, r1, and g1 when N is relatively large (Fig 3B–3E). It can be seen that as size control becomes stronger (α increases), the size distribution changes from the unimodal to the bimodal shape (Fig 3B). The size distribution is generally unimodal for timer-like strategies and bimodal for sizer-like strategies. The dependence of the size distribution on r0 is expected—a small r0 results in a small fraction of cells in the elongation phase and thus the left peak is much lower than the right peak, while a large r0 gives rise to the opposite effect (Fig 3C). Bimodality is the most apparent when r0 is neither too large nor too small.

The influence of r1 on the cell size distribution is more complicated. Recall that a larger r1 means a larger fraction of cells in the reshaping phase and a smaller fraction of cells in the septation phase. Here since we fix r0 to be a constant and tune r1, there is little change in the fraction of cells in the elongation phase. As the septation phase becomes shorter (r1 increases), the size distribution changes from being bimodal to being unimodal and becomes more concentrated (Fig 3D). In particular, bimodality is apparent when the septation phase is relatively long, while a very short septation phase may even destroy bimodality.

Finally, we examine how the cell size distribution depends on the ratio of the growth rate in the reshaping phase to the one in the elongation phase, g1/g0, which characterizes the sharpness of the size increase in the reshaping phase. As the size addition in the reshaping phase becomes sharper (g1/g0 increases), the size distribution changes from being bimodal to being unimodal and becomes more concentrated (Fig 3E). Here we fix the mean cell size to be a constant by tuning the parameter a and thus the increase in g1 does not make the right peak shift more to the right. To our surprise, we find that bimodality is the most apparent when the growth rates in the two phases are close to each other, while a very abrupt size addition in the reshaping phase may even destroy bimodality.

To summarize, we find that small added size variability, strong size control, moderate length in the elongation phase, long septation phase, short reshaping phase, and mild size addition in the reshaping phase are capable of producing more apparent bimodality.

Analytical distribution of the birth size

In our model, the distribution of the birth size Vb can also be derived analytically in steady-state conditions. In fact, the Laplace transform for the αth power of the birth size, , is given by (see Section C in S1 Appendix for the proof) (13)

Taking the inverse Laplace transform gives the probability density function of , from which the distribution of Vb can be obtained. A special case takes place when α is large (strong size control) or when p is small (smaller daughter tracking). Under the large α or small p approximation, the term pαn is negligible for n ≥ 2 and it suffices to keep only the first term in the infinite product given in Eq (13). In this case, the laplace transform of reduces to

Taking the inverse Laplace transform gives the birth size distribution (14) where 1F1 is the confluent hypergeometric function, β0 = N0/M0 pα, and β0 = N1/M1 pα.

Actually, the birth size distribution has also been computed analytically in some simpler models. It has been shown that the birth size in those models approximately has a log-normal distribution [31] or a gamma distribution [33]. Therefore it is natural to ask whether the birth size in our model shares the same property. To see this, we illustrate the birth size distribution and its approximation by the log-normal and gamma distributions as N and α vary (Fig 4A and 4B). We find that under a wide range of model parameters, the true distribution is in excellent agreement with its log-normal approximation. However, when N and α are both small, the true distribution is severely right-skewed and deviates significantly from its gamma approximation. When N and α are both large, the true distribution becomes more symmetric and the three distributions become almost indistinguishable.

thumbnail
Fig 4. Further properties of the birth size and cell size distributions.

A: Comparison of the birth size distribution (blue curve and red circles) with its log-normal (solid grey region) and gamma (dashed green curve) approximations when N and α are small. The blue curve shows the analytical distribution obtained by taking the inverse Laplace transform of Eq (13) and the red circles show the distribution obtained from stochastic simulations. B: Same as A but when N and α are large. In A, B, the parameters are chosen as r0 = 0.6, r1 = 0.1, g0 = 0.01, g1 = 4g0, p = 0.5. The parameters N and α are chosen as N = 10, α = 0.5 in A and N = 30, α = 2 in B. C: Comparison between the cell size distributions for the model with deterministic partitioning (solid grey region) and the model with stochastic partitioning (blue curve and red circles). The blue curve shows the analytical distribution obtained by taking the inverse Laplace transform of Eq (15) and the red circles show the simulated distribution. The parameters are chosen as N = 30, r0 = 0.6, r1 = 0.1, g0 = 0.01, g1 = 4g0, α = 2, p = 0.5. For the model with stochastic partitioning, the parameter ν is chosen as ν = 200. In A-C, the parameters a, M0, M1 are chosen so that the mean cell size 〈V〉 = 3 for the model with deterministic partitioning.

https://doi.org/10.1371/journal.pcbi.1009793.g004

Influence of stochastic partitioning on the cell size distribution

Thus far, the analytical distribution of cell size is derived when the partitioning at division is deterministic. In the presence of noise in partitioning, we can also obtain an explicit expression for the cell size distribution, whose Laplace transform is given by (see Section B in S1 Appendix for the proof) (15) where f(λ) is the function given in Eq (6), is a normalization constant, and an is a sequence that can be determined by the following recursive relations: (16)

Here bn and cn are two other sequences that are defined by with (x)m = x(x + 1)…(x + m − 1) being the Pochhammer symbol. For the special case of exponential growth of cell size, there is only the elongation phase and the remaining two phases vanish. In this case, we have N1 = 0 and N = N0; the cell size distribution is still determined by Eq (15) with the sequence bn and the function f(λ) being greatly simplified as

Clearly, fluctuations in partitioning at division lead to a much more complicated analytical expression of the cell size distribution. Actually, when partitioning is stochastic, the size distribution for exponentially growing cells has been derived approximately in [35] under the assumption that noise in partitioning is very small. Here we have removed this assumption and obtained a closed-form solution of the size distribution for general non-exponentially growing cells even if noise in partitioning is very large. Recent cell lineage data suggest that the coefficient of variation of the partition ratio in fission yeast is 6%—8% under different growth conditions [4].

To see the effect of stochastic partitioning, we illustrate the cell size distributions under deterministic and stochastic partitioning in Fig 4C with the standard deviation of the partition ratio R being 7% of the mean for the latter. Clearly, the analytical solution given in Eq (15) matches the simulation results very well. In addition, it can be seen that noise in partitioning gives rise to larger fluctuations in cell size, characterized by a smaller slope of the left shoulder, an apparent decrease in the height of the left peak, and a slight decrease in the height of the right peak. The valley between the two peaks and the right shoulder are almost the same for the two models.

Correlation between cell sizes at birth and at division

In [31], it has been shown that the correlation between cell sizes at birth and at division can be used to infer the size control strategy. For the case of deterministic partitioning, since the generalized added size is hypoexponentially distributed, it is easy to obtain (see Section D in S1 Appendix for the proof) (17) where ρ(X, Y) denotes the correlation coefficient between X and Y. This characterizes the correlation between birth and division sizes, which only depends on the asymmetry of partitioning (p) and the strength of size control (α). In particular, we find that if partitioning is deterministic, the correlation is independent of the growth pattern of the cell—both exponentially and non-exponentially growing cells share the same correlation coefficient whenever they have the same p and α. Note that in Eq (17), the correlation between and , instead of the correlation between Vb and Vd, is computed. This is because of the scaling property of our model: only the generalized added size has good analytical properties, instead of the real added size VdVb.

In the presence of noise in partitioning, the formula for the correlation coefficient should be modified as (see Section D in S1 Appendix for the proof) (18) where (19)

In this case, is generally lower than pα due to partitioning noise. Interestingly, if partitioning is stochastic, the correlation between birth and division sizes not only depends on p and α, but also depends on the parameters N0, M0, N1, M1, which describe the growth pattern of fission yeast. This is very different from the case of deterministic partitioning.

Experimental validation of the theory

To test our theory, we apply it to the lineage data of cell size in haploid fission yeast that are published in [4]. This data set contains high throughput data of the whole time series of thousands of individual cells over many cell cycles, instead of data at some particular time points [26]. The monitoring of the whole time series allows an accurate inference of all model parameters as well as a deeper understanding of the full cell cycle dynamics.

In this data set, the single-cell time traces of cell area (with unit μm2) were recorded every three minutes using microfluidic devices. The experiments were performed under seven growth conditions with different media (Edinburgh minimal medium (EMM) and yeast extract medium (YE)) and different temperatures. For EMM, cells were cultured at four different temperatures (28°C, 30°C, 32°C, and 34°C), while for YE, three different temperatures (28°C, 30°C, and 34°C) were applied. For each growth condition, 1500 cell lineages were tracked and each lineage is typically composed of 50 − 70 generations. Note that for a particular cell lineage, it may occur that the cell was dead or disappeared from the channel during the measurement [4, 76]. Such lineages are removed from the data set and thus the number of lineages used for data analysis for each growth condition is actually less than 1500. In addition, we emphasize that in this data set, the size of a cell is characterized by its area, which has rarely been measured in previous experiments; more commonly used quantities are cell length and cell volume. To make our results more easily comparable to those in the literature, we convert the cell area data to cell length data (with unit μm) by using the information of mean cell diameter for each growth condition, which can be estimated from the fluorescence images provided to us by the authors of [4].

Based on the cell length data, it is possible to estimate all the parameters involved in our model for the seven growth conditions. Parameter inference is crucial since it provides insights into the size control strategy, added size variability, and complex growth pattern in fission yeast. We perform parameter inference by fitting the noisy data to two models: the model with deterministic partitioning (model I) and the model with stochastic partitioning (model II). The basic statistics of some important quantities including the birth size Vb, the septation size Vs, the division size Vd, the cell size V, the cell diameter D, and the cell cycle duration T, as well as the estimated values of all parameters for the two models are listed in Table 2. In the following, we briefly describe our parameter estimation method.

  1. 1) Estimation of p and ν. Note that the data of cell sizes just before division and just after division, Vd and , across different generations can be easily extracted from the lineage data and thus for model I, the parameter p can be simply estimated as the mean partition ratio . For model II, the parameters p and ν can be inferred by fitting the partition ratio data to a beta distribution.
  2. 2) Estimation of α. Note that the data of cell sizes at birth and at division, Vb and Vd, across different generations can be easily extracted from the lineage data. For model I, since the parameter p has been determined, the strength α of cell size control can be estimated by finding the unique value of α satisfying the equality . The inference of α for model II is much more complicated. Note that once α is determined, both K1 and K2 can be computed via Eq (19). For model II, the mean and variance for the αth power of the birth size are given by (see Section D in S1 Appendix for the proof)
    Since K1 and K2 have been determined (assuming α is known), it is possible to estimate both M0 + M1 and using the data of birth sizes. Finally, the control strength α can be estimated by finding the unique value of α satisfying Eq (18).
  3. 3) Estimation of g0/a and g1/a. For model I, the mean and variance for the αth power of the birth size are given by (see Section D in S1 Appendix for the proof)
    Since the parameters p and α have been determined, using the data of birth sizes, we are able to estimate the following two quantities: where and . Once N0 and N1 are known, both and can be solved from the above two equations and thus can be inferred. For model II, we have shown how to estimate M0 + M1 and in step 2).
  4. 4) Estimation of a, g0, and g1. For each generation, say, the kth generation, we fit the time-course data of cell size to a three-stage growth model: an exponential growth in the elongation phase, followed by a constant size in the septation phase and another round of exponential growth in the reshaping phase: where Tk and Tk+1 are two successive division times, g0(k) and g1(k) are the growth rates in the elongation and reshaping phases for the kth generation, respectively, and t0 and t1 are the initial and end times of the septation phase, respectively. By carrying out least-squares optimal fitting, we can estimate the growth rate g0(k) in the elongation phase and the growth rate g1(k) in the reshaping phase for the kth generation. Fig 5A illustrates the fitting of the time-course data to the three-stage growth model for three typical cell lineages, from which we can see that the model matches the data reasonably well. Then the parameter g0 can be determined as the mean of g0(k) across different generations. Since the time that the cell stays in the reshaping phase is very short, the estimate of g1(k) in general is not accurate. Therefore, we do not adopt this method to estimate the parameter g1. Since both g0 and have been determined, the parameter a can also be inferred. Since both a and have been estimated, the parameter g1 can be determined.
  5. 5) Estimation of N, N0, and N1. Note that once the parameters N, N0, and N1 are known, all other parameters can be inferred by carrying out steps 1)—4). Finally, we determine these three parameters by solving the following optimization problem: (20) where p(x) is the theoretical cell size distribution, is the sample cell size distribution obtained from lineage data, xi are some reference points, and M is the number of bins chosen. In other words, we estimate the three parameters by matching the theoretical and experimental cell size distributions. For model I, the theoretical distribution is determined using Eq (5), while for model II, the theoretical distribution is determined using Eq (15). Thus far, all model parameters have been determined.
thumbnail
Fig 5. Fitting experimental data to theory based on the model with deterministic partitioning (model I).

A: Fitting the time-course data of cell size (grey curve) to a three-stage growth model (red curve) for three typical cell lineages cultured in YE at 34°C. B: Experimental cell size distributions (blue bars) and their optimal fitting to model I (red curve) for seven growth conditions. Here the theoretical distributions are computed using Eq (5). C: Same as B but for the birth size distributions. Here the theoretical distributions are computed using Eq (13). D: Fitting the experimental cell size distributions (blue bars) to the one-stage model with only the elongation phase (orange curve) and the two-stage model with only the elongation and septation phases (red curve). Here the theoretical distributions are computed using Eq (5) by taking r0 = 1 and r1 = 0 for the one-stage model and by taking r1 = 0 for the two-stage model. Both simpler models fail to capture the shape of the cell size distribution. E: Scatter plot of the birth length versus the added length and the associated regression line. The slope of the regression line is significantly affected by outliers. To avoid this influence, we removed the outliers identified using Hadi’s potential-residual plot [77] and grouped the data into seven bins.

https://doi.org/10.1371/journal.pcbi.1009793.g005

thumbnail
Table 2. Parameters estimated using lineage data of cell length under seven growth conditions.

The mean and standard deviation (upper-right corner) of six variables are given: the birth size Vb, the septation size Vs, the division size Vd, the cell size V, the cell diameter D, and the doubling time T. We perform parameter inference for both model I and model II. The estimation error for each parameter was computed using bootstrap. Specifically, we performed parameter inference 50 times; for each estimation, the theoretical model was fitted to the data of 50 randomly selected cell lineages. The estimation error was then calculated as the standard deviation over the 50 repeated samplings. The estimates of p and g0 are the same for both models and thus we only list their values for model I.

https://doi.org/10.1371/journal.pcbi.1009793.t002

To test our inference method, we compare the experimental cell size and birth size distributions obtained from the lineage data (blue bars) with the theoretical ones based on the estimated parameters (red curves) under the seven growth conditions for both model I (Fig 5B and 5C) and model II (S2(A) and S2(B) Fig). It can be seen that the cell size distributions of lineage measurements for the seven growth conditions are all bimodal, while the birth size distributions are all unimodal. For the latter model, we also compare the distribution of the partition ratio with its approximation using the beta distribution (S2(C) Fig). Clearly, the theory reproduces the experimental data of fission yeast excellently. Interestingly, while our inference method only involves the matching of the theoretical and experimental cell size distributions, the theoretical birth size distribution also matches the experimental one reasonably well.

To further evaluate the performance of our model, we examine the correlation between the birth size Vb and the division size Vd, as well as the correlation between the birth size Vb and the added size VdVb. Based on the lineage data, the correlation coefficients for the seven growth conditions are listed in Table 3. The theoretical predictions of the correlation coefficients based on stochastic simulations of model I and model II with the estimated parameters are also listed in Table 3. Clearly, both models capture the birth and division size correlations perfectly. Model I slightly underestimates the birth and added size correlations, while model II slightly overestimates these correlations.

thumbnail
Table 3. Correlation coefficients between the birth and division (added) sizes for the seven growth conditions.

The experimental correlation coefficients are computed using the lineage data, while the theoretical correlation coefficients are computed using stochastic simulations based on model I and model II.

https://doi.org/10.1371/journal.pcbi.1009793.t003

A natural question is whether the lineage data used here can be described by simpler models, such as the one-stage model with only the elongation phase [35, 62] or the two-stage model with only the elongation and septation phases. Here the former only describes the exponential growth in the G2 phase, while the latter ignores the abrupt increase in cell length due to the rounding off of the new ends. To see this, we also fit the cell size distribution to the one-stage and two-stage models by using the inference method introduced above (Fig 5D). Both simpler models fail to capture the unusual shape of the cell size distribution. The one-stage model always predicts a unimodal distribution; the two-stage model predicts a unimodal distribution for EMM and a bimodal distribution for YE. While the two-stage model excellently captures the left peak, it fails to reproduce the right peak due to neglection of the reshaping phase. This suggests that our model is the simplest model that can describe the lineage data of fission yeast. The perfect match between experiments and the three-stage model and the poor match between experiments and simpler models support the main assumptions of the three-stage growth model and the choice of the rate of moving from one cell cycle stage to the next to be a power law of cell size.

The parameters estimated in Table 2 also provide some useful insights of biological interest. Compared with model II, model I gives rise to a lower estimate of α and r1, as well as a higher estimate of N, a, and g1; the estimates of other parameters are very similar for the two models. While both model I and model II capture lineage data very well, we next base our discussion on the parameters of model II since this is biologically more realistic and thus its parameter estimates are more reliable.

First, our data analysis reveals some significant differences between the two media applied. From Table 2, it can be seen that cells cultured in EMM have a relatively strong size control (large α) and a relatively large added size variability (small N), while cells cultured in YE have a relatively weak size control (small α) and a relatively small added size variability (large N). Furthermore, we find that the size control strategy in fission yeast is sizer-like for all the seven growth conditions: for model II, the strength α of size control is typically 2.0 for EMM and is typically 1.5 for YE. This is in sharp contrast to the adder strategy found in E. coli, where α is estimated to be 0.8 − 1.2 for different growth conditions [35]. Our result confirms the previous finding that fission yeast uses a sizer-like strategy to achieve size homeostasis [36, 53]. The sizer-like strategy also agrees with our theoretical result that bimodal size distributions are more likely to occur in a sizer than in a timer. Based on the lineage data, the mean septation length 〈Vs〉, i.e. the mean cell length before mitotic entry, is estimated to be 13.5 − 14.5 μm for all growth conditions, which is very stable. This is consistent with the experimental value measured in many previous papers [28, 36, 42, 78, 79]. Recall that the septation length defined here is the division size defined in those papers since the reshaping phase is assumed to belong to the previous cell cycle in the present paper.

Based on the estimated parameters for the seven growth conditions, we observed a strong linear relationship between g0 and N, as well as a weaker linear relationship between g0 and α (Fig 6A and 6B). A higher growth rate in the elongation phase gives rise to a larger number of cell cycle stages and a smaller strength of size control. The positive correlation between g0 and N implies that under unfavorable environmental conditions, a lower threshold level of the division protein is used to promote mitotic entry and trigger cell division, which causes larger added size variability. The negative correlation between g0 and α suggests that a stronger size control is required for fission yeast to adapt to unfavorable conditions.

thumbnail
Fig 6. Connection between the estimated model parameters of model II.

A: Scatter plot of the estimated N versus the estimated g0 under seven growth conditions. B: Scatter plot of the estimated α versus the estimated g0. C: Scatter plot of the estimated w0 versus the estimated g0. D: Scatter plot of the estimated g1 versus the estimated g0. The red lines shown in A-D are regression lines.

https://doi.org/10.1371/journal.pcbi.1009793.g006

To further validate our theory, we illustrate the scatter plots of the birth length versus the added length and the corresponding regression lines in Fig 5E for all growth conditions. A slope of −1 for this regression line is indicative of the sizer strategy, while a slope of 0 implies the adder strategy. Such plots have been reported many times in the literature for wild-type haploid cells and generally show a clear negative correlation with a slope between −0.7 and −0.9 [28, 36, 51, 53], while there is also a study showing that the slope can be as low as −0.37 for wild-type diploids [28]. Based on the lineage data, we find that the slope is typically −0.7 for EMM, which is consistent with the values reported in previous studies, while the slope is typically −0.5 for YE. The reason why YE has a smaller slope is probably due to the fact that cells cultured in YE have faster growth rates, which give rise to weaker size control (Fig 6B). Note that similar phenomenon has also been observed in E. coli, where an adder-like behavior was found in fast growth conditions, while a sizer-like behavior was found at low growth rates [80].

In addition, our data analysis also provides the estimation of the fractions of cell cycle in the three phases. Note that the fraction of the elongation (reshaping) phase is not simply given by r0 (r1). This is because the transition rate between cell cycle stages is an increasing function of cell size, which means that earlier (later) stages have longer (shorter) durations. The real fraction of each growth phase can be estimated approximately via Eq (12). According to our estimation, the proportion w0 of the elongation phase is about 72% − 80% of the cell cycle and the proportion w1 of the reshaping phase is only 3% − 5% for all growth conditions (Table 2). This is consistent with the previous result that cells elongate during the first ∼75% of the cell cycle [1]. In addition, the mean duration w1T〉 in the reshaping phase, i.e. the time needed for cells to form hemispherical new ends from the septum, is estimated to be 5 − 8 min for all growth conditions (except EMM at 28°C) and the mean length increase 〈Vd〉 − 〈Vs〉 in the reshaping phase, i.e. the total length of the two hemispherical end caps, is estimated to be 1.2 − 2 μm for all growth conditions. This coincides with the previous observation that the length of each daughter cell grows by ∼1 μm within 5 min after septation [3].

Interestingly, from the estimated parameters, we also observed a negative correlation between g0 and w0 and a positive correlation between g0 and g1 (Fig 6C and 6D). This means that a higher growth rate in the elongation phase is associated with a smaller proportion of the elongation phase and a higher growth rate in the reshaping phase. The negative correlation between g0 and w0 can be explained as follows: a higher growth rate in the G2 phase gives rises to a shorter time to reach the target length of ∼14 μm before mitosis and thus results in the smaller proportion of that phase. In addition, the slope of Fig 6D is estimated to be 1.88, implying that the growth rate of cell length in the rephrasing phase is about twice as large as that in the elongation phase.

An interesting feature implied by the fission yeast data is that the mean partition ratio p is 0.46 − 0.48 for all growth conditions, implying that cell division is possibly asymmetric (Table 2). Although p only slightly deviates from 0.5, this difference is significant with a p-value less than 0.001 for each growth condition according to the sensitive Z-test. This deviation may result from the asymmetry in the position of the septum which is slightly nearer the new end [68, 69]. However, we cannot exclude the possibility that the partitioning is actually symmetric and the deviation of p from 0.5 is an artifact due to the segmentation algorithm used in [4], where the old-pole tips tend to be cut (S3 Fig).

Conclusions and discussion

In this work, we proposed two detailed models of cell size dynamics in fission yeast across many generations and analytically derived the cell size and birth size distributions of measurements obtained from a cell lineage. The main feature of cell size dynamics in fission yeast is its three-stage non-exponential growth pattern: a slow growth in the elongation phase, an arrest of growth in the septation phase, and a rapid elongation in the reshaping phase. The first model assumes that (i) the cell undergoes deterministic exponential growth in the elongation and reshaping phases with the growth rate in the latter phase being greater than that in the former phase; (ii) the size remains constant in the septation phase; (iii) the size just after division is a fixed fraction of the one just before division; (iv) the cell cycle is divided into multiple effective cell cycle stages which correspond to different levels of the division protein which triggers cell division; (v) the rate of moving from one stage to the next has a power law dependence on cell size. A second model was also solved which relaxes assumption (iii) by allowing the size just after division to be a stochastic fraction of the one just before division with the fraction being sampled according to a beta distribution. Under assumptions (iv) and (v), the three typical strategies of size homeostasis (timer, adder, and sizer) are unified.

Experimentally, the cell size distribution of lineage data in fission yeast is typically bimodal under various growth conditions. This is very different from the unimodal size distribution observed in many other cell types [35]. Interestingly, the bimodal cell size distribution of fission yeast can be excellently reproduced by the analytical solutions of both models. The origin of bimodality is further investigated and clarified in detail; we find that bimodality becomes apparent when (i) the variability in added size is not too large, (ii) the strength of size control is not too weak, which implies that adder or sizer-like strategies enforce size homeostasis, (iii) the proportion of the elongation phase is neither too large nor too small, (iv) the proportion of the septation phase is large, (v) the proportion of the reshaping phase is small, and (vi) the size addition in the reshaping phase is not too sharp. We also find that fluctuations in partitioning at division has a considerable influence on the shape of the cell size distribution by declining the slope of the left shoulder, as well as lowering the heights of the two peaks.

Furthermore, we have developed an effective method of inferring all the parameters involved in both models using single-cell lineage measurements of fission yeast based on the information of (i) the partition ratio, namely, the ratio of the size just after division to the size just before division, across different generations, (ii) the mean and variance of the birth size across different generations, (iii) the correlation of cell sizes at birth and at division, and (iv) the cell size distribution. Specifically, we infer the parameters except the numbers of cell cycle stages in different phases using the information (i)-(iii) and then determine the remaining parameters by matching the theoretical and experimental cell size distributions.

We have shown that the theoretical cell size and birth size distributions provide an excellent fit to the experimental ones of fission yeast reported in [4] under seven different growth conditions. This match provides support for two implicit important assumptions of our model: (i) the cell undergoes a complex three-stage growth pattern and (ii) the speed of the cell cycle progression (the transition rate between cell cycle stages) depends on cell size in a power law form. Finally, based on matching the experimental to the theoretical cell size distributions, we have estimated all model parameters from lineage data. Simulations with the inferred parameters using distribution matching also captured the correlation between birth and division sizes, and between birth and added sizes—this provides further evidence of the accuracy of our detailed model.

Based on the estimated parameters, we confirmed the previous result that fission yeast has a sizer-like control strategy. Cells cultured in EMM have a larger added size variability and a stronger size control than those cultured in YE. The estimated values of the mean septation length, the proportion of the elongation phase, as well as the duration and length increase in the reshaping phase are all consistent with the literature. Moreover, we also observed a negative correlation between the growth rate and (i) the added size variability, (ii) the size control strength, and (iii) the proportion of the elongation phase. This reveals that stronger size homeostasis and larger cell cycle duration variability are required in slow growth conditions. The growth rate in the reshaping phase was found to be twice that in the elongation phase.

Further research aims to develop more realistic models which coordinate cell size dynamics with gene expression dynamics in fission yeast and investigate the corresponding concentration homeostasis of mRNAs and proteins [81].

Methods

All methods can be found in the main text and in S1 Appendix.

Supporting information

S1 Appendix. Mathematical details.

This file contains the mathematical details of the stochastic cell size model, as well as the detailed derivations of the cell size distribution, the birth size distribution, and the correlation between birth and division sizes.

https://doi.org/10.1371/journal.pcbi.1009793.s001

(PDF)

S1 Fig. Cell size distributions obtained using stochastic simulations.

The three distributions are obtained by generating 105, 106, and 107 stochastic trajectories, respectively. The parameters are chosen as N = 50, r0 = 0.6, r1 = 0.1, g0 = 0.01, g1 = 2g0, α = 2, p = 0.5. The parameters a, M0, M1 are chosen so that the mean cell size 〈V〉 = 3.

https://doi.org/10.1371/journal.pcbi.1009793.s002

(EPS)

S2 Fig. Fitting experimental data to theory based on the model with stochastic partitioning (model II).

A: Experimental cell size distributions (blue bars) and their optimal fitting to model II (red curve) for seven growth conditions. Here the theoretical distributions are computed using Eq (15). B: Same as A but for the birth size distributions. Here the theoretical distributions are computed using stochastic simulations. C: Same as A but for the partition ratio distributions. Here the theoretical distributions are computed using Eq (3).

https://doi.org/10.1371/journal.pcbi.1009793.s003

(EPS)

S3 Fig. Fluorescence image of fission yeast cells and the segmentation algorithm used in [4] to identify the outline of a new born cell.

At division, the segmentation algorithm tends to cut old-pole tips.

https://doi.org/10.1371/journal.pcbi.1009793.s004

(TIFF)

Acknowledgments

We are grateful to Professor Nakoaka and Professor Wakamoto for sending us the fluorescence images of fission yeast that we used to estimate the cell diameter.

References

  1. 1. Mitchison JM, Nurse P. Growth in cell length in the fission yeast Schizosaccharomyces pombe. J Cell Sci. 1985;75(1):357–376. pmid:4044680
  2. 2. Mitchison J. Growth during the cell cycle. Int Rev Cytol. 2003; p. 166–258. pmid:12921238
  3. 3. Baumgärtner S, Tolić-Nørrelykke IM. Growth pattern of single fission yeast cells is bilinear and depends on temperature and DNA synthesis. Biophys J. 2009;96(10):4336–4347. pmid:19450504
  4. 4. Nakaoka H, Wakamoto Y. Aging, mortality, and the fast growth trade-off of Schizosaccharomyces pombe. PLoS Biol. 2017;15(6):e2001109. pmid:28632741
  5. 5. Tzur A, Kafri R, LeBleu VS, Lahav G, Kirschner MW. Cell growth and size homeostasis in proliferating animal cells. Science. 2009;325(5937):167–171. pmid:19589995
  6. 6. Godin M, Delgado FF, Son S, Grover WH, Bryan AK, Tzur A, et al. Using buoyant mass to measure the growth of single cells. Nat Methods. 2010;7(5):387–390. pmid:20383132
  7. 7. Wang P, Robert L, Pelletier J, Dang WL, Taddei F, Wright A, et al. Robust growth of Escherichia coli. Curr Biol. 2010;20(12):1099–1103. pmid:20537537
  8. 8. Mir M, Wang Z, Shen Z, Bednarz M, Bashir R, Golding I, et al. Optical measurement of cycle-dependent cell growth. Proc Natl Acad Sci USA. 2011;108(32):13124–13129. pmid:21788503
  9. 9. Campos M, Surovtsev IV, Kato S, Paintdakhi A, Beltran B, Ebmeier SE, et al. A constant size extension drives bacterial cell size homeostasis. Cell. 2014;159(6):1433–1446. pmid:25480302
  10. 10. Taheri-Araghi S, Bradde S, Sauls JT, Hill NS, Levin PA, Paulsson J, et al. Cell-size control and homeostasis in bacteria. Curr Biol. 2015;25(3):385–391. pmid:25544609
  11. 11. Soifer I, Robert L, Amir A. Single-cell analysis of growth in budding yeast and bacteria reveals a common size regulation strategy. Curr Biol. 2016;26(3):356–361. pmid:26776734
  12. 12. Cermak N, Olcum S, Delgado FF, Wasserman SC, Payer KR, Murakami MA, et al. High-throughput measurement of single-cell growth rates using serial microfluidic mass sensor arrays. Nat Biotechnol. 2016;34(10):1052–1059. pmid:27598230
  13. 13. Priestman M, Thomas P, Robertson BD, Shahrezaei V. Mycobacteria modify their cell size control under sub-optimal carbon sources. Front Cell Dev Biol. 2017;5:64. pmid:28748182
  14. 14. Yu FB, Willis L, Chau RMW, Zambon A, Horowitz M, Bhaya D, et al. Long-term microfluidic tracking of coccoid cyanobacterial cells reveals robust control of division timing. BMC Biol. 2017;15(1):1–14. pmid:28196492
  15. 15. Eun YJ, Ho PY, Kim M, LaRussa S, Robert L, Renner LD, et al. Archaeal cells share common size control with bacteria despite noisier growth and division. Nat Microbiol. 2018;3(2):148–154. pmid:29255255
  16. 16. Cadart C, Monnier S, Grilli J, Sáez PJ, Srivastava N, Attia R, et al. Size control in mammalian cells involves modulation of both growth rate and cell cycle duration. Nat Commun. 2018;9(1):1–15. pmid:30115907
  17. 17. Cao Z, Grima R. Analytical distributions for detailed models of stochastic gene expression in eukaryotic cells. Proc Natl Acad Sci USA. 2020;117(9):4682–4692. pmid:32071224
  18. 18. Mitchison J, Sveiczer A, Novak B. Length growth in fission yeast: is growth exponential?-No. Microbiology. 1998;144(2):265–266. pmid:9556367
  19. 19. Cooper S. Length extension in growing yeast: is growth exponential?-Yes. Microbiology. 1998;144(2):263–265. pmid:9493363
  20. 20. Sveiczer Á, Horváth A, Buchwald P. Is there a universal rule for cellular growth?–Problems in studying and interpreting this phenomenon. FEMS Yeast Res. 2014;14(5):679–682. pmid:24866294
  21. 21. Cooper S. Distinguishing between linear and exponential cell growth during the division cycle: single-cell studies, cell-culture studies, and the object of cell-cycle research. Theor Biol Med Model. 2006;3(1):1–15. pmid:16504098
  22. 22. Cooper S. Schizosaccharomyces pombe grows exponentially during the division cycle with no rate change points. FEMS Yeast Res. 2013;13(7):650–658. pmid:23981297
  23. 23. Pickering M, Hollis LN, D’Souza E, Rhind N. Fission yeast cells grow approximately exponentially. Cell Cycle. 2019;18(8):869–879. pmid:30957637
  24. 24. Knapp BD, Odermatt P, Rojas ER, Cheng W, He X, Huang KC, et al. Decoupling of rates of protein synthesis from cell expansion leads to supergrowth. Cell Syst. 2019;9(5):434–445. pmid:31706948
  25. 25. Odermatt PD, Miettinen TP, Lemiere J, Kang JH, Bostan E, Manalis SR, et al. Variations of intracellular density during the cell cycle arise from tip-growth regulation in fission yeast. Elife. 2021;10:e64901. pmid:34100714
  26. 26. Nobs JB, Maerkl SJ. Long-term single cell analysis of S. pombe on a microfluidic microchemostat array. PloS one. 2014;9(4):e93466. pmid:24710337
  27. 27. Nagy Z, Medgyes-Horváth A, Vörös E, Sveiczer Á. Strongly oversized fission yeast cells lack any size control and tend to grow linearly rather than bilinearly. Yeast. 2021;38(3):206–221. pmid:33244789
  28. 28. Sveiczer A, Novak B, Mitchison J. The size control of fission yeast revisited. J Cell Sci. 1996;109(12):2947–2957. pmid:9013342
  29. 29. Buchwald P, Sveiczer A. The time-profile of cell growth in fission yeast: model selection criteria favoring bilinear models over exponential ones. Theor Biol Med Model. 2006;3(1):1–10. pmid:16566825
  30. 30. Horváth A, Rácz-Mónus A, Buchwald P, Sveiczer Á. Cell length growth in fission yeast: an analysis of its bilinear character and the nature of its rate change transition. FEMS Yeast Res. 2013;13(7):635–649. pmid:23848460
  31. 31. Amir A. Cell size regulation in bacteria. Phys Rev Lett. 2014;112(20):208102.
  32. 32. Iyer-Biswas S, Wright CS, Henry JT, Lo K, Burov S, Lin Y, et al. Scaling laws governing stochastic growth and division of single bacterial cells. Proc Natl Acad Sci USA. 2014;111(45):15912–15917. pmid:25349411
  33. 33. Iyer-Biswas S, Crooks GE, Scherer NF, Dinner AR. Universality in stochastic exponential growth. Phys Rev Lett. 2014;113(2):028101. pmid:25062238
  34. 34. Thomas P. Analysis of cell size homeostasis at the single-cell and population level. Front Phys (Lausanne). 2018;6:64.
  35. 35. Jia C, Singh A, Grima R. Cell size distribution of lineage data: analytic results and parameter inference. iScience. 2021;24(3):102220. pmid:33748708
  36. 36. Fantes P, Nurse P. Control of cell size at division in fission yeast by a growth-modulated size control over nuclear division. Exp Cell Res. 1977;107(2):377–386. pmid:872891
  37. 37. Nurse P, Fantes P. Cell cycle controls in fission yeast: a genetic analysis. The cell cycle. 1981; p. 85–98.
  38. 38. Moreno S, Nurse P. Regulation of progression through the Gl phase of the cell cycle by the rum1+ gene. Nature. 1994;367(6460):236–242. pmid:8121488
  39. 39. Sveiczer A, Novak B, Mitchison J. Mitotic control in the absence of cdc25 mitotic inducer in fission yeast. J Cell Sci. 1999;112(7):1085–1092. pmid:10198290
  40. 40. Neumann FR, Nurse P. Nuclear size control in fission yeast. J Cell Biol. 2007;179(4):593–600. pmid:17998401
  41. 41. Martin SG, Berthelot-Grosjean M. Polar gradients of the DYRK-family kinase Pom1 couple cell length with the cell cycle. Nature. 2009;459(7248):852–856. pmid:19474792
  42. 42. Moseley JB, Mayeux A, Paoletti A, Nurse P. A spatial gradient coordinates cell size and mitotic entry in fission yeast. Nature. 2009;459(7248):857–860. pmid:19474789
  43. 43. Coudreuse D, Nurse P. Driving the cell cycle with a minimal CDK control network. Nature. 2010;468(7327):1074–1079. pmid:21179163
  44. 44. Navarro FJ, Weston L, Nurse P. Global control of cell growth in fission yeast and its coordination with the cell cycle. Curr Opin Cell Biol. 2012;24(6):833–837. pmid:23182517
  45. 45. Horváth A, Rácz-Mónus A, Buchwald P, Sveiczer Á. Cell length growth patterns in fission yeast reveal a novel size control mechanism operating in late G2 phase. Biol Cell. 2016;108(9):259–277. pmid:27063594
  46. 46. Schmoller KM. The phenomenology of cell size control. Curr Opin Cell Biol. 2017;49:53–58. pmid:29232627
  47. 47. Allard CA, Opalko HE, Liu KW, Medoh U, Moseley JB. Cell size–dependent regulation of Wee1 localization by Cdr2 cortical nodes. J Cell Biol. 2018;217(5):1589–1599. pmid:29514920
  48. 48. Opalko HE, Nasa I, Kettenbach AN, Moseley JB. A mechanism for how Cdr1/Nim1 kinase promotes mitotic entry by inhibiting Wee1. Mol Biol Cell. 2019;30(25):3015–3023. pmid:31644361
  49. 49. Saint M, Bertaux F, Tang W, Sun XM, Game L, Köferle A, et al. Single-cell imaging and RNA sequencing reveal patterns of gene expression heterogeneity during fission yeast growth and adaptation. Nat Microbiol. 2019;4(3):480–491. pmid:30718845
  50. 50. Gerganova V, Bhatia P, Vincenzetti V, Martin SG. Direct and indirect regulation of Pom1 cell size pathway by the protein phosphatase 2C Ptc1. Mol Biol Cell. 2021;32(8):703–711. pmid:33625871
  51. 51. Scotchman E, Kume K, Navarro FJ, Nurse P. Identification of mutants with increased variation in cell size at onset of mitosis in fission yeast. J Cell Sci. 2021;134(3):jcs251769. pmid:33419777
  52. 52. Vargas-Garcia CA, Ghusinga KR, Singh A. Cell size control and gene expression homeostasis in single-cells. Curr Opin Syst Biol. 2018;8:109–116. pmid:29862376
  53. 53. Wood E, Nurse P. Pom1 and cell size homeostasis in fission yeast. Cell cycle. 2013;12(19):3417–3425. pmid:24047646
  54. 54. Facchetti G, Chang F, Howard M. Controlling cell size through sizer mechanisms. Curr Opin Syst Biol. 2017;5:86–92. pmid:32984663
  55. 55. Facchetti G, Knapp B, Flor-Parra I, Chang F, Howard M. Reprogramming Cdr2-dependent geometry-based cell size control in fission yeast. Curr Biol. 2019;29(2):350–358. pmid:30639107
  56. 56. Pan KZ, Saunders TE, Flor-Parra I, Howard M, Chang F. Cortical regulation of cell size by a sizer cdr2p. Elife. 2014;3:e02040. pmid:24642412
  57. 57. Keifenheim D, Sun XM, D¡¯Souza E, Ohira MJ, Magner M, Mayhew MB, et al. Size-dependent expression of the mitotic activator Cdc25 suggests a mechanism of size control in fission yeast. Curr Biol. 2017;27(10):1491–1497. pmid:28479325
  58. 58. Patterson JO, Rees P, Nurse P. Noisy cell-size-correlated expression of cyclin b drives probabilistic cell-size homeostasis in fission yeast. Curr Biol. 2019;29(8):1379–1386. pmid:30955932
  59. 59. Tanouchi Y, Pai A, Park H, Huang S, Stamatov R, Buchler NE, et al. A noisy linear map underlies oscillations in cell size and gene expression in bacteria. Nature. 2015;523(7560):357–360. pmid:26040722
  60. 60. Atilgan E, Magidson V, Khodjakov A, Chang F. Morphogenesis of the fission yeast cell through cell wall expansion. Curr Biol. 2015;25(16):2150–2157. pmid:26212881
  61. 61. Chao HX, Fakhreddin RI, Shimerov HK, Kedziora KM, Kumar RJ, Perez J, et al. Evidence that the human cell cycle is a series of uncoupled, memoryless phases. Mol Syst Biol. 2019;15(3). pmid:30886052
  62. 62. Nieto C, Arias-Castro J, Sánchez C, Vargas-García C, Pedraza JM. Unification of cell division control strategies through continuous rate models. Phys Rev E. 2020;101(2):022401. pmid:32168656
  63. 63. Schmoller KM, Turner J, Kõivomägi M, Skotheim JM. Dilution of the cell cycle inhibitor Whi5 controls budding-yeast cell size. Nature. 2015;526(7572):268–272. pmid:26390151
  64. 64. Barber F, Amir A, Murray AW. Cell-size regulation in budding yeast does not depend on linear accumulation of Whi5. Proc Natl Acad Sci USA. 2020;117(25):14243–14250. pmid:32518113
  65. 65. Vargas-García CA, Singh A. Elucidating cell size control mechanisms with stochastic hybrid systems. In: 2018 IEEE Conference on Decision and Control (CDC). IEEE; 2018. p. 4366–4371.
  66. 66. Nieto-Acuna CA, Vargas-Garcia CA, Singh A, Pedraza JM. Efficient computation of stochastic cell-size transient dynamics. BMC Bioinformatics. 2019;20(23):1–6. pmid:31881826
  67. 67. Totis N, Nieto C, Küper A, Vargas-García C, Singh A, Waldherr S. A population-based approach to study the effects of growth and division rates on the dynamics of cell size statistics. IEEE Control Syst Lett. 2020;5(2):725–730.
  68. 68. Johnson BF, Calleja G, Boisclair I, Yoo BY. Cell division in yeasts: III. The biased, asymmetric location of the septum in the fission yeast cell, Schizosaccharomyces pombe. Exp Cell Res. 1979;123(2):253–259. pmid:499357
  69. 69. May J, Mitchison J. Pattern of polar extension of the cell wall in the fission yeast Schizosaccharomyces pombe. Can J Microbiol. 1995;41(3):273–277. pmid:7736356
  70. 70. Jia C, Grima R. Frequency domain analysis of fluctuations of mRNA and protein copy numbers within a cell lineage: theory and experimental validation. Phys Rev X. 2021;11:021032.
  71. 71. Brenner N, Braun E, Yoney A, Susman L, Rotella J, Salman H. Single-cell protein dynamics reproduce universal fluctuations in cell populations. Eur Phys J E. 2015;38(9):102. pmid:26410847
  72. 72. Robert L, Ollion J, Robert J, Song X, Matic I, Elez M. Mutation dynamics and fitness effects followed in single cells. Science. 2018;359(6381):1283–1286. pmid:29590079
  73. 73. Zopf C, Quinn K, Zeidman J, Maheshri N. Cell-cycle dependence of transcription dominates noise in gene expression. PLoS Comput Biol. 2013;9(7):e1003161. pmid:23935476
  74. 74. Crane MM, Clark IB, Bakker E, Smith S, Swain PS. A microfluidic system for studying ageing and dynamic single-cell responses in budding yeast. PloS one. 2014;9(6):e100042. pmid:24950344
  75. 75. Nieto-Acuña C, Arias-Castro JC, Vargas-García C, Sánchez C, Pedraza JM. Correlation between protein concentration and bacterial cell size can reveal mechanisms of gene expression. Phys Biol. 2020;17(4):045002. pmid:32289764
  76. 76. Kamimura A, Kobayashi TJ. Representation and inference of size control laws by neural-network-aided point processes. Phys Rev Research. 2021;3:033032.
  77. 77. Chatterjee S, Hadi AS. Regression analysis by example. John Wiley & Sons; 2015.
  78. 78. Kamasaki T, Arai R, Osumi M, Mabuchi I. Directionality of F-actin cables changes during the fission yeast cell cycle. Nat Cell Biol. 2005;7(9):916–917. pmid:16136186
  79. 79. Das M, Wiley DJ, Medina S, Vincent HA, Larrea M, Oriolo A, et al. Regulation of cell diameter, For3p localization, and cell symmetry by fission yeast Rho-GAP Rga4p. Mol Biol Cell. 2007;18(6):2090–2101. pmid:17377067
  80. 80. Wallden M, Fange D, Lundius EG, Baltekin Ö, Elf J. The synchronization of replication and division cycles in individual E. coli cells. Cell. 2016;166(3):729–739. pmid:27471967
  81. 81. Sun XM, Bowman A, Priestman M, Bertaux F, Martinez-Segura A, Tang W, et al. Size-dependent increase in RNA Polymerase II initiation rates mediates gene expression scaling with cell size. Curr Biol. 2020;30(7):1217–1230. pmid:32059768