Abstract
We consider a size-structured bacterial population model in which the rate of cell growth is both size- and time-dependent and the average per capita reproduction rate is specified as a model parameter. It is shown that the model admits classical solutions. The population-level and distribution-level behaviours of these solutions are then determined in terms of the model parameters. The distribution-level behaviour is found to be different from that found in similar models of bacterial population dynamics. Rather than convergence to a stable size distribution, we find that size distributions repeat in cycles. This phenomenon is observed in similar models only under special assumptions on the functional form of the size-dependent growth rate factor. Our main results are illustrated with examples, and we also provide an introductory study of the bacterial growth in a chemostat within the framework of our model.
1. Introduction
When a bacterial population grows, there are two fundamental types of ‘growth’ that take place: growth in biomass and growth in cell number. The former type of growth results from the increase in mass of individual cells and the latter type results from binary fission. In his pioneering work on bacterial growth, Monod Citation19 pointed out that
When the average size of cells does not change in the time interval considered, the increase in bacterial density is proportional to the increase in cell concentration. Whether growth is estimated in terms of one variable of the other, the growth rate is the same. However, as established particularly by the classical studies of Henrici, the average size of the cells may vary considerably from one phase to another of a growth cycle. It follows that the two variables, cell concentration and bacterial density, are not equivalent. Much confusion has been created because this important distinction has been frequently overlooked.
The benefit gained by ‘overlooking’ the distinction between biomass and cell number is that this simplification allows for the formulation of tractable mathematical models in which the state of a population can be described by a single time-dependent variable. Such models essentially view the population as a single ‘blob’ whose mass (or some other particular characteristic) is continuously changing. Models that do distinguish between growth in biomass and growth in cell number must include separate mechanisms for describing cell growth and fission. Such size-structured models were first developed by Sinko and Streifer Citation22 Citation23, Bell Citation3 and Bell and Anderson Citation4 and were later modified and investigated extensively by Diekmann et al. Citation7 Citation8 and Metz and Diekmann Citation17.
The model that we develop and analyse in this paper is a size-structured model that is similar to the model studied in Citation7
Citation8
Citation17. The major difference is that our model contains a balance law (a linear partial differential equation) that describes only the cell growth process. The reproductive process is described by a parameter, q, which we call the per capita reproductive quota, and which is introduced as a boundary condition to accompany the balance law. In contrast, both the growth and reproductive processes are described by a single balance law in the model studied in Citation7
Citation8
Citation17. The different modelling approach that we use is found to result in different dynamics. We find that the cell size distributions that occur in all solutions of our model repeat in cycles, whereas the generic behaviour of solutions of the model in Citation7
Citation8
Citation17 is the convergence of all solutions to a stable size distribution (with cyclic distributions occurring only under certain assumptions on the cell growth rate). The different dynamics that we find at the distribution level also results in different dynamics at the overall population level. In particular, we find that, in accordance with the observations of Monod given above, our model admits solutions for which the average cell size does not remain constant (or approach a constant value as ).
The remainder of this paper is organized as follows. In Section 2, we formulate our model and make some preliminary observations. The general solution of the model is constructed in Section 3, and the behaviour of the general solution is described in terms of the model parameters in Section 4. Two examples that illustrate our results are provided in Section 5. In Section 6, we adapt our model to a chemostat setting and provide a preliminary analysis of this chemostat model. A summary and discussion that compares our model to the model studied in Citation7 Citation8 Citation17 is given in Section 7.
2. Formulation of the model
To formulate our model of bacterial growth and reproduction, we begin by assuming that all cells in the population are born with size m>0 and grow to size 2m, at which point they undergo fission. We also assume that cell growth takes place according to a law of the form
(H m ) The minimum cell size, m, satisfies m>0. | |||||
(H
g
) The size-dependent growth factor, | |||||
(H
f
) The time-dependent growth factor, | |||||
(H q ) The per capita reproductive quota, q, satisfies 0<q≤2. | |||||
(H D ) The per capita removal rate, D, satisfies D≥0. |
Throughout our work, we will use the term ‘size’ to mean cell biomass. Thus m and x have units of mass. In addition, we will assume that g has units of mass and that f has units of time−1. In H
q
, the term ‘per capita’ refers to the average number of daughters produced per mother and hence q has units of number/number (or no units). In H
D
, the term ‘per capita’ refers to the percentage of the population that is removed per unit of time and thus D has units of .
EquationEquation (1) together with hypotheses H
g
and H
f
implies that a cell that is born with size m at a given time t will (if it is viable and is not removed from the population) reach size 2m and undergo fission at time t* where
As is customary, we define ρ(t, x) to be the density with respect to the cell size, x, of the population at time t. Thus ρ has units of number/mass. We take the domain of ρ to be the closed set . At any time t≥0, the total number of cells in the population, P(t), and the total biomass of the population, B(t), are
By incorporating H m –H D , we obtain the model
Equation Equation(2a), which gives a balance law for cell growth and removal, can be derived using the approach of Diekmann et al. [Citation8, p. 247]. These authors give the derivation in the case that the cell growth rate depends only on cell size, x, and not explicitly on t, but their derivation is easily modified to yield Equation Equation(2a)
in the non-autonomous case. The balance law derived in Citation8 also includes terms that account for removal of larger cells and creation of smaller cells due to fission, but we do not include these terms because the reproductive process in our model is described completely by the boundary condition (2b), for which we provide a derivation in the appendix. The initial size distribution, φ, is assumed to be a function of class C
1 with
for all x∈[m, 2m]. Throughout this paper, the term ‘distribution’ will always refer to a function of this class.
Our objective is to describe the behaviour of all solutions of model (2). We will investigate the behaviour of solutions at both the population level and the distribution level. By behaviour at the population level, we mean the evolution over time of the total population, P, and the total biomass, B. By behaviour at the distribution level, we mean the evolution over time of the size distributions, , where for each t≥0,
is the distribution that has value ρ(t, x) for each x∈[m, 2m].
2.1 Some basic observations
If ρ is a positive-valued solution of model (2), then
In addition, the biomass satisfies
3. The general solution
To set the stage for constructing the general solution of model (2), we define
3.1 Change of variables
Model (2) can be simplified by introducing new variables, τ and ξ, defined as
3.2 Invariant size distributions
We say that a distribution ψ≥0 is an invariant size distribution (or, for brevity, an invariant distribution) on [0, G] if system (3) admits a solution of the form . Clearly ψ≡0 is an invariant distribution. By separation of variables, we find that any non-zero invariant distribution must satisfy
In terms of the original variables from model (2), solutions corresponding to invariant distributions can be expressed as
At the distribution level, Equation Equation(5a) implies that φ can have an extremum at any value of x, where
. It is thus clearly possible to construct functions, g, for which φ has one or more extrema in [m, 2m]. This observation is relevant to solving the ‘inverse problem’ of finding appropriate g (and other model parameters) that will fit a given empirically observed distribution. Solving the inverse problem has been found to be problematic in some other size-structured bacterial population models. For example, the model of Cushing Citation6 discussed in Smith and Waltman [Citation26, Chapter 9] was shown to admit only monotone-decreasing invariant distributions and this precluded being able to fit the model to the unimodal invariant distributions observed in the experiments of Williams Citation28. In the present model, the size-dependent growth factor, g(x), can be reconstructed directly from the invariant distribution, φ(x), by solving Equation Equation(5a)
:
3.3 The general solution
Any solution, , of Equation Equation(3a)
remains constant along the characteristic curves
. If
also satisfies Equation (3b), then for all τ≥0, we have
Now we can express the general solution in terms of the original variables. Since and
, the compatibility conditions (6) correspond to the following compatibility conditions for φ:
At this point, we introduce two auxiliary extensions. We let g
m
(x) be the m-periodic extension of g from [m, 2m) onto R. We also extend the initial distribution, φ, by defining for all
and all n∈Z. The (global) characteristics now take the form
Let X(t, x) be the point of intersection of the characteristic passing through (t, x) with the x-axis. Due to EquationEquation (8), X(t, x) must satisfy
4. Behaviour of solutions
We are now prepared to study the population and distribution-level behaviour of all solutions of Equation (2). Throughout the remainder, we assume that the initial distribution, φ, satisfies the compatibility conditions (7) and hence that ρ, as defined by EquationEquation (9), satisfies Equation (2) in the classical sense throughout Ω.
4.1 Population-level dynamics
In Section 3.2, it was seen that solutions of Equation (2) that correspond to invariant distributions have the form where
The total population is given by
Theorem 1
Let ρ be the solution of Equation (2) corresponding to the initial distribution φ and
(i) For all t≥0, | |||||
(ii) If q¬=1, the ratio P(t)/π(t) remains constant if and only if the solution ρ corresponds to an invariant distribution. If q=1, the ratio P(t)/π(t) remains constant for all solutions of Equation (2). | |||||
(iii) For all t≥0, |
Proof
First we will prove claim (i) in the case q≥1: let , where
and n is a non-negative integer. Since
for all ξ∈R, we observe that
Theorem 1 tells us that the population-level behaviour of all solutions of Equation (2) is determined by the time evolution of π(t), which is in turn determined by the parameters f, g, q, and D. If these parameters are such that (∞) as
, then also
(∞) as
. If π(t) is bounded, then so is P(t), and we obtain an estimate on its bounds in terms of m, g, q, and φ. Since π(t) does not depend on φ, the population-level behaviour of all solutions is qualitatively the same (not dependent on the initial size distribution). In many models of structured population dynamics, the distribution-level behaviour of all solutions is also found to be the same in that all solutions are seen to converge to an invariant distribution as
. In what follows, we will see that this behaviour does not occur in model (2).
4.2 Distribution-level dynamics
It will now be shown that the size distributions of all solutions of Equation (2) occur in cycles. Since for all ξ∈R, we have
Theorem 2
If ρ is a solution of Equation (2) and
are such that
Theorem 2 tells us that, for any given compatible φ, all possible size distributions that are achieved by the corresponding solution of Equation (2) occur in the interval [0, t(G)). These same distributions then occur again (in the sense of equivalence) in each of the intervals . In particular, no solutions converge to invariant distributions (other than those solutions that actually correspond to invariant distributions). Furthermore, for a given φ, the actual size distributions that occur do not depend on f but the frequencies with which these distributions repeat themselves do depend on f. If f
1 and f
2 are two different functions (both satisfying H
f
) and ρ1 and ρ2 are the corresponding solutions of Equation (2) with
, then for any t≥0, we will have
for infinitely many t*≥0. However, the population-level behaviour (growth, decay, or boundedness) of ρ1 and ρ2 might differ greatly depending on the behaviour of f
1 and f
2.
5. Examples
To illustrate the results given in Sections 3 and 4, we now consider two examples. In each of these examples, we assume that f(t) is constant for all t≥0. Each example could easily be modified to allow other choices of f(t), but we do not do so here for the sake of brevity.
5.1 Size-proportional growth rate
Assume that cells grow according to , where k>0 is constant. In this case, we can take f(t)=k and g(x)=x to obtain
The compatibility conditions (7) from which classical solutions of Equation (2) can be constructed are
(1) If | |||||
If | |||||
If |
To investigate statement 3 a bit further, we will show that the critical case does not imply that either the population, P(t), or the biomass, B(t), remains constant. In fact, we will show that both P(t) and B(t) are periodic functions of period ln(2)/k which are, in general, not constant.
As we did previously, it is more convenient to work with P and B as functions of τ:
Assuming that q≠1 (and hence D≠0), we obtain from part (ii) of Theorem 1 that P(t) does not remain constant unless φ is an invariant distribution. A similar argument, which requires the assumption that k≠D and q≠2, shows that B(t) does not remain constant unless φ is an invariant distribution. The two subcases of the critical case which have to be checked separately are Equation(1)
q=1, D=0 and (2) q=2, k=D. In the former subcase, P(t) remains constant and B(t) is periodic (and not constant unless φ is an invariant distribution). In the latter subcase, B(t) remains constant and P(t) is periodic (and not constant unless φ is an invariant distribution).
5.2 Constant growth rate
If cells grow according to (where a>0 is constant), then we can take f(t)=1 and g(x)=a. In this case,
In addition, the compatibility conditions are
(1) If aln(q)<Dm, then | |||||
If aln(q)>Dm, then | |||||
If aln(q)=Dm (implying that q≥1), then |
A noteworthy feature of this example is that the minimum cell size, m, is found to have a role in determining the population-level behaviour. For given D>0, a>0, and q>1, we will have aln(q)<Dm (and ) if m is sufficiently large and aln(q)>Dm (and
) if m is sufficiently small. The intuitive reason for this is that, with a constant cell growth rate, smaller cells take less time (than do larger cells) to reach maturity and are thus able to reproduce at a rate that exceeds the removal rate.
6. Application to the chemostat
The standard unstructured model for microbial growth in a chemostat is
To model a bacterial population in a chemostat within the framework of the size-structured model (2), we begin with the assumption that cell size increases according to a law of the form where μ is a Monod function. The equations that describe growth and reproduction are thus
A general study of model Equation(15)–(19) will be the subject of a future investigation. In order to pave the way for this future work, we will examine the chemostat dynamics that arise from the model in the important case of perfect reproductive efficiency (q=2) and size-proportional growth rate (g(x)=x). In this case, EquationEquations (15)–(17) take the form
We remark that Equations Equation(20a)–Equation(20d)
were derived from the first principles by Smith in Citation24 by assuming the size-proportional cell growth rate (g(x)=x) and perfect reproductive efficiency (q=2) and by then passing from a discrete version of the model to Equations Equation(20a)
–(20d). As was shown in the analysis provided in [Citation24, p. 753], since Equations Equation(20c)
and Equation(20d)
are identical to the classical chemostat model (14), the biomass–substrate dynamics are the same in the structured model as in the unstructured model: if μ(S
f
)>D, then
where
. However, we should point out that the corresponding size distribution will oscillate, unless the initial distribution is invariant. Thus, the attractivity of the steady state claimed by Smith is valid only for solutions corresponding to invariant distributions, e.g. precisely those considered in Citation24. However, the structured model provides a richer description of the population dynamics as we now illustrate by examining solutions for which the biomass–substrate steady state has already been achieved. To construct such solutions, we first set
in Equation Equation(20a)
and choose any compatible distribution φ, fixed but arbitrary, such that
7. Discussion
Early models of bacterial population dynamics attempted to describe the behaviour of such populations in terms of a single observable feature, usually the cell number or biomass. These models were improved upon by taking the population structure into account. The development of structured models has proceeded according to two fundamental approaches which we will refer to very loosely as state structuring and age/size structuring. In the state structure approach, the cell cycle is viewed as being comprised of a discrete set of states and the process by which cells transfer between these states is modelled. The most well-known model of this type is the one introduced by Smith and Martin Citation25 in which the cell cycle is viewed as consisting of two states, A and B, with the probabilistic time of residence in the A state and deterministic time of residence in the B state. The Smith–Martin idea has been modified and expanded upon in the work of Tyson and Hannsgen Citation27, Grasman Citation9, Pilyugin et al. Citation21, and others. The age/size structure approach differs from the state structure approach in that it classifies cells according to age and/or body size and/or some other identifiable physical characteristics. Some structured models, such as those developed by Lasota and Mackey Citation16 and Basse et al. Citation2, incorporate elements of both state structure and age/size structure. A general survey of the structured cell population models is given in Arino Citation1.
The model (2) that we have investigated here is a size-structured model where ‘size’ is interpreted to mean cell biomass. Under very general assumptions on the model parameters, we show that the model admits classical solutions and determine the behaviour of these solutions in terms of the model parameters. This behaviour is found to be fundamentally different than that typically found in similar models. Instead of convergence of all solutions to a stable size distribution, we obtain cyclic size distributions that correspond to populations in which the average cell size fluctuates in cycles that are possibly of variable length.
A general concern in the process of mathematical modelling is the trade-off that must be made between realism and the feasibility of mathematical analysis. Model (2) incorporates the processes of bacterial growth and reproduction in a fundamental way that allows for the construction of classical solutions whose behaviour can be clearly elucidated. However, we sacrifice some realism by hypothesizing that all cells undergo fission at the same critical size (2m). A similar simplifying assumption was made in the model studied by Cushing Citation6 and Smith and Waltman Citation26. This model allows for variation of size at fission but requires that all cells have the same birth size. In reality, empirical studies performed on cultures of Escherichia coli Citation10 Citation11 have shown that the coefficient of variation for size at fission is generally around 10% to 16% and that this does not depend significantly on the medium in which the cells are being cultured. The coefficient of variation of age at fission is much larger. By assuming that all cells divide at the same size, we are able to introduce the reproductive parameter, q, which keeps track of the reproductive process and enables us to determine how the interplay between growth and reproduction governs the population dynamics.
A model that bears close resemblance to model (2) is model (21) which was studied by Diekmann et al. in Citation7 Citation8 and by Metz and Diekmann in Citation17.
Model (21) contains a balance law that describes both the growth and reproductive processes and hence allows for the fact that fission can occur over a range of sizes. Specifically, the maximum possible cell size is normalized to be 1, it is assumed that cells divide into two equal parts, and the minimum size at which cells can undergo fission is assumed to be a<1. This implies that all newborn cells lie in the size range [a/2, 1/2] and that fission takes place only in cells that are in the size range [a, 1]. Cells of size x are assumed to divide with probability b(x). The detailed analysis of this model given in Citation7 Citation8 Citation17 proceeds by first transforming the model into an abstract evolution equation in an appropriate Banach space and by then applying the theory of semigroups of linear operators. However, solutions of the abstract evolution equation cannot typically be transformed back into classical solutions of Equation (21) and are thus interpreted as satisfying Equation (21) only in a certain weak sense. The generic dynamical behaviour of model (21) is the convergence of all solutions to an asymptotically stable invariant size distribution in which the average cell size is constant. However, if the typical assumption of the size-proportional cell growth rate (g(x)=x) is made, then model (21) does not have a stable size distribution and all solutions have cyclic size distributions just as in model (2).
We conclude with a discussion of the modelling role of the function g(x). In most studies, it is assumed that cell size increases at a size-proportional rate (g(x)=x). In this case, models (2) and (21) are in agreement on their predictions of the population dynamics. However, the question of whether the cell growth rate can be generally expected to be size-proportional appears to have still not been settled. Based on an extensive set of experiments performed with Escherichia coli and other bacteria, Kubitschek Citation12 Citation13 Citation14 Citation15 argues that the rate of the increase in size of these bacteria is constant throughout most of the cell cycle. On the other hand, Cooper Citation5 strongly takes the position that the rate of cell growth is indeed size proportional. Empirical results obtained by Meyer et al. Citation18 indicate that the growth rate is probably neither constant nor size proportional but that it varies across different stages of the cell cycle and for different species. The ‘inverse problem’ of finding g(x) based on empirical observations of cell size distributions is thus clearly an issue that is of interest to microbiologists. We are hopeful that size-structured models such as the ones we have introduced here can be helpful in furthering such investigations.
Acknowledgements
The research by Pilyugin was partially supported by NSF grant NSF DMS–0818050.
References
- Arino , O. 1995 . A survey of structured cell population dynamics . Acta Biotheor. , 43 : 3 – 25 .
- Basse , B. , Baguley , B. C. , Marshall , E. S. , Joseph , W. R. , van Brunt , B. , Wake , G. and Wall , D. J.N. 2003 . A mathematical model for analysis of the cell cycle in cell lines derived from human tumors . J. Math. Biol. , 47 : 295 – 312 .
- Bell , G. I. 1968 . Cell growth and division III. Conditions for balanced exponential growth in a mathematical model . Biophysics , 8 : 431 – 444 .
- Bell , G. I. and Anderson , E. C. 1967 . Cell growth and division I. A mathematical model with applications to cell volume distributions in mammalian suspension cultures . Biophysics , 7 : 329 – 351 .
- Cooper , S. 1988 . What is the bacterial growth law during the division cycle? . J. Bacteriol. , 170 ( 11 ) : 5001 – 5005 .
- Cushing , J. M. 1989 . A competition model for size-structured species . SIAM J. Appl. Math. , 49 : 838 – 858 .
- Diekmann , O. , Lauwerier , H. A. , Aldenberg , T. and Metz , J. A.J. 1983 . Growth, fission and the stable size distribution . J. Math. Biol. , 18 : 135 – 148 .
- Diekmann , O. , Heijmans , H. J.A.M. and Thieme , H. R. 1984 . On the stability of the cell size distribution . J. Math. Biol. , 19 : 227 – 248 .
- Grasman , J. 1990 . A deterministic model of the cell cycle . Bull. Math. Biol. , 56 : 535 – 547 .
- Koch , A. L. and Schaechter , M. 1962 . A model for statistics of the cell division process . J. Gen. Microbiol. , 29 : 435 – 454 .
- Koppes , L. J.H. , Woldringh , C. L. and Nanninga , N. 1978 . Size variations and correlation of different cell cycle events in slow-growing Escherichia coli . J. Bacteriol. , 134 ( 2 ) : 423 – 433 .
- Kubitschek , H. E. 1968 . Linear growth in Escherichia coli . Biophys. J. , 8 ( 7 ) : 792 – 804 .
- Kubitschek , H. E. 1968 . Constancy of uptake during the cell cycle in Escherichia coli . Biophys. J. , 8 ( 12 ) : 1401 – 1412 .
- Kubitschek , H. E. 1969 . Growth during the bacterial cycle: Analysis of cell size distribution . Biophys. J. , 9 ( 6 ) : 792 – 809 .
- Kubitschek , H. E. 1986 . Increase in cell mass during the division cycle of Escherichia coli B/rA . J. Bacteriol. , 168 ( 2 ) : 613 – 618 .
- Lasota , A. and Mackey , M. C. 1999 . Cell division and the stability of cellular populations . J. Math. Biol. , 38 : 241 – 261 .
- Metz , J. A.J. and Diekmann , O. 1986 . The Dynamics of Physiologically Structured Populations , Berlin : Springer . Lecture Notes in Biomathematics 68
- Meyer , M. , De Jong , M. A. , Demets , R. and Nanninga , N. 1979 . Length growth of two Escherichia coli B/r strains . J. Bacteriol. , 138 ( 1 ) : 17 – 23 .
- Monod , J. 1949 . The growth of bacterial cultures . Ann. Rev. Microbiol. , 3 : 371 – 394 .
- Panikov , N. S. 1995 . Microbial Growth Kinetics , London : Chapman and Hall .
- Pilyugin , S. S. , Ganuov , V. V. , Murali-Krishna , K. , Ahmed , R. and Antia , R. 2003 . The rescaling method for quantifying the turnover of cell populations . J. Theor. Biol. , 225 : 275 – 283 .
- Sinko , J. W. and Streifer , W. 1967 . A new model for age-size structure of a population . Ecology , 48 : 910 – 918 .
- Sinko , J. W. and Streifer , W. 1971 . A model for populations reproducing by fission . Ecology , 52 : 330 – 335 .
- Smith , H. L. 1996 . A discrete size-structured model of microbial growth in a chemostat . J. Math. Biol. , 34 : 734 – 754 .
- Smith , J. A. and Martin , L. 1973 . Do cells cycle? . Proc. Natl Acad. Sci. USA , 70 : 1263 – 1267 .
- Smith , H. L. and Waltman , P. 1995 . The Theory of the Chemostat , New York : Cambridge University Press .
- Tyson , J. J. and Hannsgen , K. B. 1986 . Cell growth and division: A deterministic/probabilistic model of the cell cycle . J. Math. Biol. , 23 : 231 – 246 .
- Williams , F. M. 1971 . Dynamics of microbial populations, in Systems Analysis and Simulations in Ecology , Edited by: Patten , B. C. 197 – 267 . New York : Academic Press .
Appendix. Derivation of boundary condition (2b)
For fixed t≥0 and small h>0, the time that it takes for a cell that has size 2m−h at time t to grow to size 2m is τ(h, t), where
During the time interval , a cell that has size m will grow to size
where x is the solution of
, x(t)=m.
Our assumption that each mother produces q daughters is thus modelled by the requirement that
The right-hand side of the above equation is equal to qρ(2m, t). Also, by L'Hopital's rule, we have
We thus require that