Journal of Theoretical Biology 457 (2018) 170–179
Contents lists available at ScienceDirect
Journal of Theoretical Biology
journal homepage: www.elsevier.com/locate/jtb
Selection for synchronized cell division in simple multicellular organisms
∗ Jason Olejarz a, , Kamran Kaveh a, Carl Veller a,b, Martin A. Nowak a,b,c a Program for Evolutionary Dynamics, Harvard University, Cambridge, MA 02138, USA b Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA c Department of Mathematics, Harvard University, Cambridge, MA 02138, USA
a r t i c l e i n f o a b s t r a c t
Article history: The evolution of multicellularity was a major transition in the history of life on earth. Conditions un- Received 15 March 2018 der which multicellularity is favored have been studied theoretically and experimentally. But since the
Revised 30 July 2018 construction of a multicellular organism requires multiple rounds of cell division, a natural question is Accepted 29 August 2018 whether these cell divisions should be synchronous or not. We study a population model in which there Available online 30 August 2018 compete simple multicellular organisms that grow by either synchronous or asynchronous cell divisions. Keywords: We demonstrate that natural selection can act differently on synchronous and asynchronous cell division, Evolutionary dynamics and we offer intuition for why these phenotypes are generally not neutral variants of each other. Multicellularity ©2018 The Authors. Published by Elsevier Ltd. Synchronization
Cell division This is an open access article under the CC BY license. (http://creativecommons.org/licenses/by/4.0/)
1. Introduction the organism is large enough, it starts to produce progenitor cells, which disperse to seed the growth of new multicellular organisms. The evolution of multicellular organisms from unicellular an- In simplest terms, the rate of growth (and ultimately the pro- cestors was one of the major transitions in the evolutionary his- ductivity) of a multicellular organism depends on the rate of di- tory of life, and occurred at least 25 times independently, begin- vision of its cells at various stages of its growth. For example, if ning as far back as 3–3.5 billion years ago ( Bonner, 1988, 1998; a k -complex produces new cells k times faster than a unicellular Grosberg and Strathmann, 2007; Knoll, 2003; Niklas and Newman, organism does, then the per-cell division rate of the multicellular 2013; Rokas, 2008; Maynard Smith and Szathmary, 1995 ). Progress organism is always equal to that of the unicellular organism, and has been made in elucidating conditions that select for simple, multicellularity is no more productive than unicellularity ( Bonner, undifferentiated multicellularity over unicellularity, both in theory 1998; Tarnita et al., 2013; Willensdorfer, 2009 ). But this represents ( Driscoll and Travisano, 2017; Ghang and Nowak, 2014; Michod, only a single possibility. More generally, it is natural to consider 2007; Niklas, 2014; Olejarz and Nowak, 2014; Pfeiffer and Bon- cases where selection acts differently on complexes of different hoeffer, 2003; van Gestel and Tarnita, 2017 ) and in experimental sizes ( Tarnita et al., 2013; Willensdorfer, 20 08, 20 09 ). For exam- work ( Herron et al., 2018; Ratcliff et al., 2012, 2013a; Solari et al., ple, if each k -complex produces new cells at a rate more than k 2006a,b; Tarnita et al., 2015 ). times faster than a unicellular organism, then the ST phenotype Staying together and coming together are mechanisms for the outcompetes the solitary phenotype, and multicellularity evolves. formation of biological complexity ( Tarnita et al., 2013 ). Many sim- Natural selection may also act in non-linear, non-monotonic, or ple multicellular organisms grow by the staying together (ST) of frequency-dependent ways on complexes of different sizes ( Celiker dividing cells, starting from a single progenitor cell ( Grosberg and and Gore, 2013; Julou et al., 2013; Koschwanez et al., 2013; Lavren- Strathmann, 1998 ). Thus, a progenitor cell divides, and the daugh- tovich et al., 2013; Ratcliff et al., 2013b; Tarnita, 2017 ), and for ter cell stays attached to the parent cell to form a complex of two many interesting cases, the population dynamics of ST are well cells (a ‘2-complex’). Further cell divisions produce new cells that characterized ( Allen et al., 2013; Ghang and Nowak, 2014; Kaveh stay with the group, and the organism grows. Eventually, when et al., 2016; Maliet et al., 2015; Michod, 2005; Michod et al., 2006; Momeni et al., 2013; Olejarz and Nowak, 2014; van Gestel and Nowak, 2016 ). ∗ Corresponding author. Against the background of this rich set of possibilities for the E-mail addresses: [email protected] (J. Olejarz), [email protected] (K. Kaveh), [email protected] (C. Veller), fitness effects of multicellularity, a question that has been ignored [email protected] (M.A. Nowak). (to our knowledge) concerns the timing of cell divisions in the con- https://doi.org/10.1016/j.jtbi.2018.08.038 0022-5193/© 2018 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license. ( http://creativecommons.org/licenses/by/4.0/ ) J. Olejarz et al. / Journal of Theoretical Biology 457 (2018) 170–179 171 struction of a multicellular organism. Specifically, should their tim- φy ( α ; y ) such that y 1 + 2 y 2 + 3 y 3 + 4 y 4 = 1 at all times. We obtain ing be independent or temporally correlated? That is, can there be selection for synchrony in cell division? Here, we study a model of 4 φ ( α ; ) = + (α − ) . simple multicellularity to determine the conditions under which y y 1 kyk k 1 (2) = synchronized cell division is favored or disfavored. k 1 The dynamics of synchronous cell division and staying together = 2. Model for n 4 are described by the following system of differential equations: We suppose that new cells remain attached to their parent cell x˙ 1 −α1 0 4 α4 x 1 x 1 after cell division. This process continues until a complex reaches x˙ 2 = α1 −α2 0 x 2 − φx ( α ; x ) x 2 . (3) its maximum size, n . A complex of size n then produces new soli- x˙ 4 0 α2 0 x 4 x 4 tary cells. Here, the variables x1 , x2 , and x4 denote the frequencies of com-
First, consider a population of asynchronously dividing cells. For plexes with 1, 2, and 4 synchronously dividing cells, respectively. asynchronous cell division, the reproduction of each individual cell α ≤ ≤ (The parameters k for 1 k 4 do not depend on synchroniza- is a Poisson process, and cells divide independently. For illustra- tion or asynchronization in cell division. Therefore, α = { α } is k tion, consider the case of neutrality. The distribution of time inter- defined exactly as for the case of asynchronous cell division, as vals between production of new cells is exponential, with an aver- α described above, although in the case of synchronization, the 3 age rate of a single cell division in one time unit. In one time unit, value in α is irrelevant.) In Eq. (3) , we choose φx ( α ; x ) such that on average, a single cell reproduces to form a complex containing x 1 + 2 x 2 + 4 x 4 = 1 at all times. We obtain two cells (the parent and the offspring). With asynchronous cell 2 → → division, it takes only another 1/2 time unit, on average, for either p φ α; x = 1 + kx ( α − 1 ), where k = 2 . (4) of the cells of the 2-complex to reproduce and form a 3-complex. x k k p=0 Once the 3-complex appears, in another 1/3 time unit, on average, one of the three cells of the 3-complex will reproduce to form a 4- In what follows, we use an asterisk to denote quantities that { ∗} complex. If n = 4 , then each 4-complex produces new solitary cells are in steady-state. For asynchronously dividing cells, yk denote = at a rate of 4 cells per time unit, and the cell division and staying the frequencies of k-complexes when y˙ k 0 for all k. Also for φ∗( α ) together process starting from each new solitary cell is repeated. asynchronously dividing cells, y denotes the population fitness = φ∗( α ) (For a more detailed explanation, see Appendix A .) when y˙ k 0 for all k. (Notice that y is equal to the largest real Next, consider a population of synchronously dividing cells. For eigenvalue of the matrix on the right-hand side of Eq. (1) , and this synchronous cell division, all cells in a k -complex divide simul- quantity represents the growth rate of the population (if we ne- taneously, and simultaneous division of an entire cluster of cells glect death of cells) when that matrix multiplies the vector of com- is a Poisson process. The growth process starting from a single plex frequencies. A higher growth rate then requires a larger com- φ∗( α ) cell is subtly different if cells divide synchronously. For illustration, pensating value of y in order to keep the population size con- φ∗( α ) again consider the case of neutrality. The distribution of time inter- stant. As such, y can be viewed as an overall death rate due to { ∗} vals between doubling of an entire cluster of cells is exponential, overcrowding.) Similarly, for synchronously dividing cells, xk and φ∗( α ) with an average rate of one doubling of a cluster’s size in a single x denote the frequencies of k-complexes and the population = time unit. In one time unit, on average, a single cell reproduces to fitness, respectively, when x˙ k 0 for all k. The processes of staying form a 2-complex. In one time unit, on average, the two cells in together with synchronous and asynchronous cell division for the the 2-complex simultaneously divide, the result being a new com- case n = 4 are shown schematically in Fig. 1 . plex with four cells—the two parent cells and the two offspring. Suppose that we have a mixed population of cells with the syn- (Notice that 3-complexes do not form if cell division is perfectly chronously and asynchronously dividing phenotypes. Notice that φ∗ < φ∗, synchronous.) If n = 4 , then each 4-complex produces new solitary there is competitive exclusion. If x y then asynchronously di- φ∗ > φ∗, cells at a rate of 4 cells per time unit, and each new solitary cell viding cells outcompete synchronously dividing cells. If x y repeats the cell division and staying together process. then synchronously dividing cells outcompete asynchronously di- φ∗ = φ∗, viding cells. If x y then synchronously and asynchronously di-
3. Results viding cells can coexist. From Eqs. (2) and (4) , for the partic- α = α = α = α = , 3.1. n = 4 cells ular case 1 2 3 4 1 we have φ ( , , , ; ) = φ ( , , , ; ) = φ∗( , , , ) = φ∗( , , , ) = , x 1 1 1 1 x y 1 1 1 1 y x 1 1 1 1 y 1 1 1 1 1
We begin by studying the evolutionary dynamics for n = 4 . The and there is neutrality between the synchronously and asyn- α = α = α = α = dynamics of asynchronous cell division and staying together for chronously dividing phenotypes. But the case 1 2 3 4 α = n = 4 are described by the following system of differential equa- 1 is nongeneric. What happens if k 1 for some k? α tions: To understand the effect of k on the evolutionary dynam-
⎛ ⎞ ⎛ ⎞ ⎛ ⎞ ⎛ ⎞ ics, we consider a couple of simple cases. First, we consider (α , α , α , α ) = ( , , , α) y˙ 1 −α1 0 0 4 α4 y 1 y 1 the case 1 2 3 4 1 1 1 . The difference between the ⎜ y˙ ⎟ ⎜ α −2 α 0 0 ⎟ ⎜ y ⎟ ⎜y ⎟ steady-state growth rates of the synchronously and asynchronously 2 = 1 2 2 − φ ( α ; ) 2 . ⎝ ⎠ ⎝ ⎠ ⎝ ⎠ y y ⎝ ⎠ φ∗( , , , α) φ∗( , , , α) , y˙ 3 0 2 α2 −3 α3 0 y 3 y 3 dividing subpopulations, x 1 1 1 and y 1 1 1 respec- α = , φ∗( , , , α) y˙ 4 0 0 3 α3 0 y 4 y 4 tively, is plotted in Fig. 2A. If 1 then x 1 1 1 and φ∗(1 , 1 , 1 , α) are exactly equal, as already noted. If α < 1, then (1) y φ∗( , , , α) < φ∗( , , , α) , x 1 1 1 y 1 1 1 and asynchronously dividing cells α < The notation y˙ i indicates the time derivative. Here, the variables outcompete synchronously dividing cells. But if 1, then a soli- y1 , y2 , y3 , and y4 denote the frequencies of complexes with 1, 2, tary phenotype would outcompete the staying-together phenotype, 3, and 4 asynchronously dividing cells, respectively. The parame- and there would be no formation of clusters in the first place. If α ≤ ≤ α > φ∗( , , , α) > φ∗( , , , α) , ters k for 1 k 4 represent the consequences of staying together 1, then x 1 1 1 y 1 1 1 and synchronously divid- on the fitness of cells in k -complexes. We use the shorthand nota- ing cells outcompete asynchronously dividing cells. Thus, for α > 1, α = { α } α tion k to denote the set of k values. In Eq. (1), we choose multicellularity with synchronized cell division evolves. 172 J. Olejarz et al. / Journal of Theoretical Biology 457 (2018) 170–179
Fig. 1. Growth of multicellular organisms by synchronous and asynchronous cell division, when maximum size is n = 4 cells. (A) In synchronous cell division, a solitary cell divides to form a 2-complex, and then both cells of the 2-complex divide simultaneously to form a 4-complex. Further cell divisions in the 4-complex give rise to dispersing single cells. When all cells are equally productive in terms of their own division rates, no matter what size complex they are in (the neutral case), a steady state is reached where 1/2 of organisms are single cells (accounting for 1/4 of all cells), 1/4 are 2-complexes (1/4 of all cells), and 1/4 are 4-complexes (1/2 of all cells). 3-complexes are never produced. (B) In asynchronous cell division, a solitary cell divides to form a 2-complex, after which one of the cells in the 2-complex divides to form a 3-complex, after which one of the cells in the 3-complex divides to form a 4-complex. Further divisions lead to dispersing single cells. In the neutral case, an asynchronous population reaches a steady state where 1/2 of complexes are single cells (6/25 of all cells), 1/6 are 2-complexes (4/25 of all cells), 1/12 are 3-complexes (3/25 of all cells), and 1/4 are 4-complexes (12/25 of all cells).
Fig. 2. The difference in the steady-state fitnesses of the synchronously and asynchronously dividing subpopulations is shown for two choices of α for n = 4 . When multi- cellularity is selected for ( α > 1), synchronization of cell divisions is the more successful phenotype in some cases (A), and asynchronization in others (B).
What is the intuition for this result? To develop insight, we re- chronously dividing cells produce new cells at the same reduced turn to the simple case of α = 1 . In steady state, what fraction of rate 1 − . In this case, asynchronization is the more successful the total number of cells in the synchronously and asynchronously phenotype. α dividing subpopulations belong to 4-complexes? The steady-state Our intuition further suggests that, for different values of ( 1 , α α α distributions of cluster sizes are shown in Fig. 1. For the syn- 2 , 3 , 4 ), the asynchronous phenotype can outcompete the syn- chronous phenotype, exactly 1/2 of all cells belong to 4-complexes chronous phenotype, under conditions in which multicellularity in steady state. For the asynchronous phenotype, exactly 12/25 will evolve. Consider again the steady-state distributions of k - ( < 1/2) of all cells belong to 4-complexes in steady state. The frac- complexes for the case of α = 1 , as shown in Fig. 1 . What frac- tion of cells that belong to 4-complexes—and are therefore affected tion of all cells in the synchronous and asynchronous subpopula- by small changes in α—is larger for the synchronous phenotype. tions belong to k -complexes of size k ≥ 2? In steady state, 3/4 of This suggests that, if α differs from 1 by a small amount, then the all synchronously dividing cells belong to complexes with at least corresponding effect on the population’s fitness—either positive or 2 cells, while 19/25 ( > 3/4) of all asynchronously dividing cells be- negative—is amplified for synchronously dividing cells. long to complexes with at least 2 cells. Therefore, we anticipate For example, if α = 1 + with 0 < 1, then approximately that, for the fitness values (α1 , α2 , α3 , α4 ) = (1 , α, α, α) , and for 1/2 of all synchronously dividing cells produce new cells at an α = 1 + with 0 < 1, the asynchronous phenotype is more suc- enhanced rate 1 + , while only approximately 12/25 of all asyn- cessful, and multicellularity is also evolutionarily preferred over chronously dividing cells produce new cells at the same enhanced the solitary phenotype. rate 1 + . In this case, synchronization is the more successful phe- The difference between the steady-state growth rates of notype. If, instead, α = 1 − with 0 < 1, then approximately the synchronously and asynchronously dividing subpopulations, φ∗( , α, α, α) φ∗( , α, α, α) , 1/2 of all synchronously dividing cells produce new cells at a x 1 and y 1 respectively, is plotted in Fig. 2B. reduced rate 1 − , while only approximately 12/25 of all asyn- Our expectation is correct: If α < 1, then the synchronous phe- J. Olejarz et al. / Journal of Theoretical Biology 457 (2018) 170–179 173 notype outcompetes the asynchronous phenotype, but a solitary then the synchronous phenotype would be more successful, but phenotype would also outcompete the staying-together phenotype, multicellularity is not selected. As shown in Fig. 4 B, our expecta- and there would be no formation of clusters. If α > 1, then the tion is again correct. asynchronous phenotype outcompetes the synchronous phenotype, and evolutionary construction develops. 3.3. Any number of cells 3.2. n = 8 cells We have considered the cases n = 4 and n = 8 , but more gen-
We can also consider the case n = 8 , for which a complex con- erally, we can describe the evolutionary dynamics for any value of tains a maximum of 8 cells. The dynamics of asynchronous cell the maximum complex size, n. The dynamics of asynchronous cell division and staying together for n = 8 are described by the fol- lowing equations: ⎛ ⎞ ⎛ ⎞ ⎛ ⎞ ⎛ ⎞ y˙ 1 −α1 0 0 0 0 0 0 8 α8 y 1 y 1 ⎜ y˙ 2 ⎟ ⎜ α1 −2 α2 0 0 0 0 0 0 ⎟ ⎜ y 2 ⎟ ⎜ y 2 ⎟ ⎜ ⎟ ⎜ α − α ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ y˙ 3 ⎟ ⎜ 0 2 2 3 3 0 0 0 0 0 ⎟ ⎜ y3 ⎟ ⎜ y3 ⎟ ⎜ y˙ 4 ⎟ ⎜ 0 0 3 α3 −4 α4 0 0 0 0 ⎟ ⎜ y 4 ⎟ ⎜ y 4 ⎟ ⎜ ⎟ = ⎜ ⎟ ⎜ ⎟ − φy ( α ; y ) ⎜ ⎟ . (5) ⎜ y˙ 5 ⎟ ⎜ 0 0 0 4 α4 −5 α5 0 0 0 ⎟ ⎜ y 5 ⎟ ⎜ y 5 ⎟ ⎜ ⎟ ⎜ α − α ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ y˙ 6 ⎠ ⎝ 0 0 0 0 5 5 6 6 0 0 ⎠ ⎝ y6 ⎠ ⎝ y6 ⎠ y˙ 7 0 0 0 0 0 6 α6 −7 α7 0 y 7 y 7 y˙ 8 0 0 0 0 0 0 7 α7 0 y 8 y 8
In Eq. (5) , we choose φy ( α ; y ) such that the total number of cells equals one at all times. We obtain division and staying together are described by the following equa- tions: 8 φ ( α ; y ) = 1 + ky (α − 1) . (6) y k k y˙ 1 = −α1 y 1 + nαn y n − φy ( α ; y ) y 1 ; k =1 = ( − ) α − α − φ ( α ; ) ∀ < < ; y˙ k k 1 k −1 yk −1 k k yk y y yk 1 k n The dynamics of synchronous cell division and staying together for y˙ = (n − 1) α − y − − φ ( α ; y ) y . n = 8 are described by the following equations: n n 1 n 1 y n ⎛ ⎞ ⎛ ⎞ ⎛ ⎞ ⎛ ⎞ We choose φy ( α ; y ) such that the total number of cells equals one x˙ −α 0 0 8 α x x 1 1 8 1 1 at all times. We obtain ⎜ x˙ ⎟ ⎜ α −α 0 0 ⎟ ⎜ x ⎟ ⎜x ⎟ 2 = 1 2 2 − φ ( α ; ) 2 . ⎝ ⎠ ⎝ ⎠ ⎝ ⎠ x x ⎝ ⎠ n x˙ 4 0 α2 −α4 0 x 4 x 4 φy ( α ; y ) = 1 + ky (α − 1) . (9) x˙ 8 0 0 α4 0 x 8 x 8 k k k =1 (7) The dynamics of synchronous cell division and staying together are In Eq. (7) , we choose φ ( α ; x ) such that the total number of cells x described by the following equations: equals one at all times. We obtain = −α + α − φ ( α ; ) ; 3 x˙ 1 1 x1 n n xn x x x1 → → p φ α; = + ( α − ), = . = α − α − φ ( α ; ) ∀ = p, < < ; x x 1 kxk k 1 where k 2 (8) x˙ k k/ 2 xk/ 2 k xk x x xk k 2 0 p log2 n p=0 x˙ n = αn/ 2 x n/ 2 − φx ( α ; x ) x n . α = Notice that, from Eqs. (6) and (8), if k 1 for all k, then φ ( α ; ) = φ ( α ; ) = φ∗( α ) = φ∗( α ) = , We choose φx ( α ; x ) such that the total number of cells equals one y y x x y x 1 and there is neutrality. The processes of staying together with synchronous and asyn- at all times. We obtain = chronous cell division for the case n 8 are shown schematically log 2 n → → in Fig. 3 . φ α; = + α − , = p. x x 1 kxk ( k 1) where k 2 (10)
To see what this means for the evolutionary dynamics, con- p=0 sider the steady-state fraction of total cells in the synchronously α = and asynchronously dividing subpopulations that belong to 8- From Eqs. (9) and (10), it follows that, if k 1 for all k, then φ ( α ; ) = φ ( α ; ) = φ∗( α ) = φ∗( α ) = , complexes if there is neutrality. The steady-state distributions of y y x x y x 1 and there is neutrality. cluster sizes are shown in Fig. 3 . For the synchronous case, that But the same intuition that applies for the cases n = 4 and n = 8 α = fraction is 2/5, while for the asynchronous case, that fraction also applies for larger values of n: If k 1 for some value of is 280/761 ≈ 0.368 ( < 2/5). This suggests that, for the case α = k , then, in general, natural selection will act differently on syn- (1 , 1 , 1 , 1 , 1 , 1 , 1 , α) , and for α = 1 + with 0 < 1, multicellu- chronous and asynchronous phenotypes. larity is selected and the synchronous phenotype is more success- ful. If instead α = 1 − with 0 < 1, then the asynchronous phe- 3.4. Different fitnesses for cells of asynchronous and synchronous notype would be more successful, but multicellularity is not se- phenotypes lected. As shown in Fig. 4 A, our expectation is correct. Consider also the steady-state fraction of total cells in the syn- The consideration of the same fitness, α , of cells in a k - chronously and asynchronously dividing subpopulations that be- k complex between the asynchronous and synchronous phenotypes long to complexes containing at least 2 cells if there is neutrality. thus reveals that natural selection acts differently on the two For the synchronous case, that fraction is 4/5, while for the asyn- modes of cellular division. More generally, one can consider chronous case, that fraction is 621/761 ≈ 0.816 ( > 4/5). This sug- that cells in k -complexes that divide asynchronously versus syn- gests that, for the case α = (1 , α, α, α, α, α, α, α) , and for α = 1 + chronously do not necessarily have the same fitness. Denote by β with 0 < 1, multicellularity is selected and the asynchronous k the fitness of cells in an asynchronously dividing k -complex. We phenotype is more successful. If instead α = 1 − with 0 < 1, 174 J. Olejarz et al. / Journal of Theoretical Biology 457 (2018) 170–179
Fig. 3. Growth of multicellular organisms by synchronous and asynchronous cell division, when maximum size is n = 8 cells.
Fig. 4. The difference in the steady-state fitnesses of the synchronously and asynchronously dividing subpopulations is shown for two choices of α for n = 8 . When multi- cellularity is selected for ( α > 1), synchronization of cell divisions is the more successful phenotype in some cases (A), and asynchronization in others (B).
β = γ have the following dynamics: Here, in general, k k . What is the condition for the syn- chronous phenotype to be favored over the asynchronous pheno- ˙ = −β + β − ( β ; ) ; Y1 1 Y1 n n Yn Y Y Y1 type, or vice versa? = ( − ) β − β − ( β ; ) ∀ < < ; ˙ = Y˙k k 1 k −1 Yk −1 k k Yk Y Y Yk 1 k n For the asynchronously dividing phenotype, setting Yk 0 for all k , we have the following homogeneous linear system: Y˙ n = (n − 1) βn −1 Y n −1 − Y ( β;Y )Y n . n ( β ; ) We choose Y Y such that the total number of cells equals one ∗ ∗ 0 = B ( β, ( β ))Y ∀ 1 ≤ i ≤ n. at all times. We obtain ij Y j j=1 n ( β ; ) = + (β − ) . Y Y 1 kYk k 1 For the synchronously dividing phenotype, setting X˙ = 0 for all k , k =1 k we have the following homogeneous linear system: γ Denote by k the fitness of cells in a synchronously dividing k- + complex. We have the following dynamics: 1 log2 n ∗ ∗ − = ( γ , ( γ )) = j 1 0 Cij X Xk where k 2 X˙ 1 = −γ1 X 1 + nγn X n − X ( γ ; X ) X 1 ; j=1
˙ = γ − γ − ( γ ; ) ∀ = p, < < ; Xk k/ 2 Xk/ 2 k Xk X X Xk k 2 0 p log2 n ∀ ≤ ≤ + . 1 i 1 log2 n X˙ n = γn/ 2 X n/ 2 − X ( γ ; X ) X n . For the asynchronously dividing phenotype, using a cofactor ex- We choose ( γ ; X ) such that the total number of cells equals one ∗ X pansion, we can solve implicitly for : at all times. We obtain Y