Emergent multicellular life cycles in filamentous bacteria owing to density-dependent population dynamics

Filamentous bacteria are the oldest and simplest known multicellular life forms. By using computer simulations and experiments that address cell division in a filamentous context, we investigate some of the ecological factors that can lead to the emergence of a multicellular life cycle in filamentous life forms. The model predicts that if cell division and death rates are dependent on the density of cells in a population, a predictable cycle between short and long filament lengths is produced. During exponential growth, there will be a predominance of multicellular filaments, while at carrying capacity, the population converges to a predominance of short filaments and single cells. Model predictions are experimentally tested and confirmed in cultures of heterotrophic and phototrophic bacterial species. Furthermore, by developing a formulation of generation time in bacterial populations, it is shown that changes in generation time can alter length distributions. The theory predicts that given the same population growth curve and fitness, species with longer generation times have longer filaments during comparable population growth phases. Characterization of the environmental dependence of morphological properties such as length, and the number of cells per filament, helps in understanding the pre-existing conditions for the evolution of developmental cycles in simple multicellular organisms. Moreover, the theoretical prediction that strains with the same fitness can exhibit different lengths at comparable growth phases has important implications. It demonstrates that differences in fitness attributed to morphology are not the sole explanation for the evolution of life cycles dominated by multicellularity.


INTRODUCTION
Multicellularity is an organizational characteristic present in the majority of organisms whose size surpasses microbial scales. Phylogenetic inference suggests that evolutionary transitions to multicellularity have occurred several times during the history of life [1 -6]. Nonetheless, despite its fundamental importance, it is difficult to empirically study the evolutionary and ecological forces that may lead to a transition from single-celled to multicellular organization. One set of theoretical explanations are based on examining the consequences of a shift between units of selection. These approaches consider the changing boundaries of an individual after a transition to multicellularity, and how the different activities of component cells can potentially lead to synergies that increase the fitness of the multicellular assemblage [1,5,7,8]. Experimental results concerning social evolution in microbes are also largely compatible with this perspective [9 -14].
Multicellular bacteria with a filamentous form are ubiquitous in many environments, including freshwater, oceans, soil, extreme habitats and the human body [30][31][32][33]. Extensive empirical work has been done to examine and monitor filamentous bacteria that can be toxic or problematic in the environment [34][35][36][37]. Some filament-forming cyanobacteria also develop specialized terminally differentiated cells, named heterocysts, that fix nitrogen and allow for the division of labour [38,39]. Although differentiation can represent a clear evolutionary advantage, theoretical and phylogenetic evidence suggests that in the cyanobacterial case, undifferentiated multicellularity evolved prior to differentiated multicellularity [40]. In aquatic bacteria, the most common hypothesis for the advantage of filamentation-and hence undifferentiated multicellularity-is an increase in size and the concomitant defence against predation by grazers [25][26][27][41][42][43]. However, experiments indicate that the avoidance of predatory grazers is not the only factor causing an increased frequency of filamentous bacteria in aquatic environments [44]. Notably, in some bacterial species, filament formation appears to be dependent on the growth state of the population, whereby an increase in the dilution rate of chemostat cultures leads to longer filament lengths [27,45].
No theoretical studies address the distribution of filament lengths and the population dynamics leading to shorter or longer filaments, although differences in length can reflect the extent to which a species is able to maintain multicellularity. Environmental conditions such as temperature, solar irradiation and nutrient concentrations have been detected as factors in determining the mean size (filament length) of different cyanobacterial species [46 -48]. Several of these factors also contribute to competition between species and adaptation to different niches [49,50]. Filament breakage can occur because of external mechanical stress, lytic processes initiated by pathogens [51][52][53], or programmed cell death [54][55][56][57][58].
In the present study, we deliberately avoid a modelling framework in which one assumes an a priori fitness advantage to multicellularity. We also do not aim to provide mechanistic explanations for the 'origin' of multicellularity. Rather, by considering cellular growth and division in a filamentous context, we investigate the constraints that would affect filament formation and the maintenance of multicellularity. This gives an indication of the capability of bacterial organisms to evolve multicellular life cycles according to their life-history traits. Our model is built on the basic assumption that the growth of a population of filamentous bacteria and the changes in the length of the filaments can be set in an ecological framework. In classical population dynamics, the change in population size is governed by the processes of birth and death (and sometimes migration), combined by the well-known logistic equation owing to Verhulst. In this model, birth and death rates are usually assumed to be decreasing and increasing functions of the population density, respectively. When the birth and death rates are the same, the population density is at a steady state and is said to be at its carrying capacity. The value at which the two rates are equal will be here referred to as the turnover rate u. The same carrying capacity can be achieved with different birth and death rate functions, and therefore at different corresponding turnover rates (figure 1a).
Our approach to study distributions of bacterial filament length places the described density-dependent concepts in the context of a population of cells that can form filamentous individuals. Although prokaryotes had been commonly thought to be in principle immortal [59,60], there is now evidence for programmed cell death and lysis [61], and in addition ageing and senescence [62]. Taking this into account, a multicellular filament is made up of several cells, and a population of filaments can also be considered as a population of cells. In this case, the population size is given by the total number of cells, while the birth and death rates coincide with the birth and death rates of each cell of the filaments. In such a population, one can vary the turnover at which a given fixed carrying capacity is achieved and ask the question if different turnovers correspond to different growth strategies; namely whether during its growth cycle the population is composed mostly of a few long filaments, or of many short filaments. Such a line of questioning is possible because as it will be shown, the rate at which a population reaches its carrying capacity can be different from the rate at which it reaches its stationary filament length distribution. The present work is an attempt to model and experimentally test the dynamics of filamentous bacteria by only considering basic life-history traits of the population. The distribution of filament lengths in the different growth phases can provide insights into the developmental strategies of different filament-forming species, and the population conditions influencing multicellular development.
In the following sections, we first state the model assumptions and we specify the experimental set-up for the culture of five filamentous bacterial species. Since turnover is a key factor in our model, we chose species that also exhibit different turnovers. We hence cultured heterotrophic species with fast generation time and high turnover, and photoautotrophic species characterized by a slower generation time and a low turnover. Heterotrophs differ from photoautotrophs in the way they obtain organic carbon for growth. While the former absorb organic compounds from the environment, the latter autonomously fix carbon using solar light as an energy source. We present the simulation results and test the model through a comparison with our experimental results.

Theoretical model
We consider a population of multicellular, undifferentiated, non-branching filamentous bacteria, whose growth is regulated only by the density-dependent birth and death rates. In this context, the total number of cells of the population is used as the population size. The birth and death rates are represented by birth and death probabilities per iteration per-cell. At every iteration, each cell in a filament has a certain probability of dividing into two daughter cells or of dying. Based on these probabilities, the number of cell divisions and cell lyses are calculated in each filament. Accordingly, the filaments are elongated and then broken into smaller pieces. We simulate the change in time of the population until it has stabilized at its carrying capacity. At each iteration, the following steps are where N c is the total number of cells and N c * the carrying capacity. The latter is defined as the value of the population size, such that b (N c  . We refer to the value of the rates at the carrying capacity as the turnover rate u. The plot illustrates that the same carrying capacity can be achieved by species that have a different turnover (solid lines, species 1; dashed lines, species 2). (b) For all birth and death functions following the scheme in (a) N c reaches the steady state following the same growth curve despite the differing turnovers. The dashed line indicates the time at which the carrying capacity is reached (see mathematical definition of N c * in §3.1.1). Error bars represent standard deviation of the mean. (c) When the population size reaches the carrying capacity (dashed line), the higher the turnover, the higher the total number of filaments N f . Nonetheless, after many more iterations, N f reaches a steady state of about 2500 filaments for all sufficiently small turnovers (plot not shown). Error bars represent standard deviation of number of filaments. (d) Box plots of filament lengths when N c reaches the carrying capacity. Populations with low turnover have higher median length values than populations with high turnover. The corresponding pie charts on the top show the proportion of single-celled, double-celled and longer filaments in the population.
As b and d are the per-cell probabilities in a filament of length i, the total number of events b or d can be viewed as the number of successes in i trials, each one with probability of success b or d. By picking from a binomial distribution B, we determine b and d for each filament of length i: b ¼ random pick from Bði; bÞ and d ¼ random pick from Bði; dÞ: The filament is elongated by b cells. Then, a number of d cells computed on the basis of the original filament length i are randomly selected for lysis in the elongated filament. The filaments resulting from the breakages are stored for the next generation. After completing the scan of all the filaments, the new population is set for the next iteration. The choice of initial conditions and the order at which birth and death are applied to filaments do not significantly affect the results of the simulations. The default initial condition for the simulations is a single cell. Alternative algorithms implementing a reverse or a random order of birth and death steps have been tested and produced similar results (figures not shown).
As stated in point (ii), this model implements linear birth and death rate functions, with fixed carrying capacity. Alternatively, nonlinear functions have been implemented and provided qualitatively comparable results (see the electronic supplementary material, §S2). Considering the difference between those rates, given by b 2 d ¼ 2(c 1 þ c 2 )N c þ 1, one can observe that the net growth rate (c 1 þ c 2 ) is independent of the turnover u, as c 1 þ c 2 ¼ 1/N c * . This implies that regardless of the turnover, the population size is expected to grow with the same curve for any turnover.

Relation between turnover and generation time
Before presenting the results, we provide a formula that helps to understand the connection between the model parameters and the experimental data. In the theoretical model, the turnover is the characterizing property of a strain. However, in natural bacterial populations, turnover is not an easily measurable quantity. For bacteria, doubling time during the exponential growth phase is usually calculated in lieu of generation time. Two populations with the same growth curve also have the same doubling time. Nonetheless, the reality is that the internal dynamics leading to this growth curve can vary. In the same time period, the same concentration of bacteria can be achieved by undertaking many births and deaths, or by conversely having less births and also less deaths. This leads to different turnovers and generation times. Moreover, at carrying capacity, the population density is stable, hence there is no indication of generation time at that stage by measuring doubling time. Here, we calculate the generation time of bacterial populations at carrying capacity in relation to turnover rate.
According to a common approach in ecology to describe the dynamics of an age-structured population [63], the mean age of parents at childbirth in a population with stable age distribution is calculated as: where x is the age, l x is the probability of survival from conception to age x and m x is the mean number of offspring at age x. G can also be interpreted as a measure of generation time without making any assumption about the growth rate of the population [64]. In our case, age is measured in a number of iterations. We consider bacterial cells as individuals that give birth by cell division and die by a lytic process. Suppose that a cell survives a given iteration with probability 1 2 d. Then the probability that a cell is living after x iterations if d remains constant is given by If also the probability of cell division is b at every iteration, the average number of offspring at a given age x coincides with the probability of birth b, whereby m x ¼ b. The generation time of a cell can then be expressed as: where the last equivalence was obtained by using geometric series. The standard formula in equation (2.3) applies when the growth rate is constant. Given the logistic nature of our model however, the functions b and d are actually dependent on the population size N c . Only when the population approaches the carrying capacity, namely when N c ! N c * , we have that b, d ! u, hence the growth rate is constant (equal to 0). We then present a formula for the generation time at the carrying capacity. Substituting b ¼ d ¼ u in equation (2.4), the generation time is given by which is a decreasing function of the parameter u (graphical illustrations in the electronic supplementary material, figure S2). Hence, a population with low turnover has a long generation time, whereby cell birth and deaths occur at a low rate. On the other hand, populations with high turnover have a short generation time, whereby births and deaths occur rapidly.  A and B).
To extrapolate the cell number from the filament length, an estimated size of cells forming a filament was measured from several light microscopy pictures.
From at least 60 observations, the estimated mean cell lengths are 6.0, 6.7 and 5.9 mm for R. lutea, F. aestuarina and F. limi, respectively. These values are higher than those measured from the inoculum of the same bacteria, indicating that cell size of heterotrophs may change over time. For cyanobacteria, the mean cell length is less variable, and an approximate cell length of 4.0 mm was used for the conversion.

Simulation results
3.1.1. Turnover and fitness. For the rest of this article, we will refer to the time at which the population reaches its carrying capacity as the first time t, such that N c (t þ 1)/N c (t) 1 þ e, where N c (t) is the population size at time t and e ¼ 1%. The growth curve of the total number of cells gives information on the fitness of the strains. Figure 1b shows that for any turnover, the total number of cells reaches the carrying capacity after less than 20 generations with the same growth curve (average and standard deviation of 1000 runs is shown). In our case, fitness is represented by net growth rate and does not depend on the turnover ( §2.1). Hence, populations with different turnovers have the same fitness. On the other hand, when the total number of cells reaches the carrying capacity, the total number of filaments (N f ) differs significantly according to the turnover (figure 1c). At this point, N f ranges from values close to the carrying capacity, corresponding to almost unicellular filaments (u ¼ 0.9), to values under 20 units (u ¼ 0.001). Figure 1b,c shows that the steady states of number of cells and number of filaments are not reached at the same time. This is because the stationary filament length distribution is reached much later than the steady state population density (carrying capacity). In the long run, the number of cells per filament reaches a stable equilibrium of about two cells per filament for sufficiently small turnovers, typically less than 0.1 (figure not shown). An analytical support for this result has been obtained by deriving an analogous model in continuous time, as presented in the electronic supplementary material, §S9. For higher turnovers, the effect of discretization leads to a deviation from this equilibrium. Hence, filament length distributions can be affected by turnover rates in temporary phases of the population dynamics, before a stationary distribution is reached. An example is given by the time when the carrying capacity of the population is reached (17 iterations in this case).
3.1.2. Length distribution at the time when carrying capacity is reached. Figure 1d shows the box plot of filament lengths when the total number of cells in the population (N c ) has reached the carrying capacity N c * ¼ 5000 cells. This time point however represents an intermediate stage before a stationary distribution of filament lengths is reached. The corresponding piecharts indicate the fraction of single cells, doublecelled filaments and filaments with more than two cells in the population. When turnover is close to 1 (u ¼ 0.9), the population consists mainly of unicellular and double-celled filaments. Note that u ¼ 0.9 means that at carrying capacity, each cell has a 90 per cent chance of dividing as well as a 90 per cent chance of subsequently dying. When u ¼ 0.1, filaments with one or two cells represent about half of the population. For lower turnovers (u ¼ 0.01, 0.001), the median values of the length are significantly higher. In the latter cases, filaments with more than two cells are more than 90 per cent. Figure 1d shows that, as the turnover increases, the fraction of short filaments when the population reaches its carrying capacity also increases. The distribution means at this growth stage were statistically compared by means of a one-way ANOVA and a multiple comparison (see the electronic supplementary material, §S6). The means obtained for u ¼ 0.9 and u ¼ 0.1 were not significantly different. All the other pairwise comparisons were significant. Figure 2 shows two examples where the mean length of the filaments was calculated every iteration. The plot shows the average of 1000 runs, with error bars indicating the 2.5th and 97.5th percentiles for the mean. The dotted line indicates the iteration at which the total number of cells reaches the carrying capacity. An increase to a peak of the mean length followed by a decay towards a steady state is observed in all cases. At the end of the decreasing phase, the mean length is comparable to that in the beginning. The same pattern was also observed for u ¼ 0.01, 0.001 (see the electronic supplementary material, §S3). An analytical characterization of the filament length dynamics has been carried out in the electronic supplementary material ( §S8), where an approximation for the iteration of first filament breakage has been derived. In order to statistically analyse the average length trend, we considered three reference time points and we compared the distribution of the means at those points using a non-parametric hypothesis test. As reference points we chose the starting point, the peak and the carrying capacity. The outcome of the tests provides statistical evidence of a significant increase and subsequent decrease of the average length. Details and results are given in the electronic supplementary material, §S6.

Length distribution during the transient phase.
Significant differences are found in the maximum average length achieved. Long filaments, with a mean length up to 250 cells, are observed at a turnover of u ¼ 0.1. Longer filaments can be achieved with lower turnovers (electronic supplementary material, figure S4). For a higher turnover (u ¼ 0.9), the average length of filaments at its peak does not exceed 100 cells. The observed trend indicates the prevalence of a cycle in the mean length induced by the birth and death rates, whereby filaments are short in the beginning, longer in the transient phase and again short at the carrying capacity. Moreover, in the transient phase, filaments with low turnover can elongate more than filaments with high turnover.

Experimental results and model validation
In the theoretical model, the population stays at the carrying capacity once this has been reached. However, as bacterial strains were not cultured in a chemostat, carrying capacities could not be maintained indefinitely. In order to assure that the experimental dataset and the simulation data are comparable, we indicate with a dotted line the time point where the bacterial populations reach their carrying capacity on each plot of the experimental results. This point was determined from the OD curves (not shown). Subsequent to reaching the carrying capacity, the process of filament breakage is eventually accelerated in the empirical case owing to a progressive decrease in the available nutrients.
In order to prove the significant increase and decrease of bacterial filament length, we performed an ANOVA test associated with a multiple comparison on the data based on Tukey's HSD test [69]. Details of the tests for each species and within species are provided in the electronic supplementary material, §S6.
In the following sections, the time point 0 is represented by the inoculum. The starting culture of  figure S6). The mean length of cyanobacteria, tracked for two months, shows a qualitative trend similar to that of heterotrophs, albeit much slower. Both the A. variabilis and N. muscorum species have a mean-length peak in the transient phase. However, because of the low turnover, they need more time than heterotrophs to reach the stationary distribution, as predicted by the model. When they reach the carrying capacity, they are still long. The pattern observed in heterotrophs is comparable with the simulation results shown in figure 2. The trend of the cyanobacteria, although not so markedly conclusive, gives an indication that the model predictions also hold in autotrophs. An additional illustration of the mean length trend of the cultured species is provided in the electronic supplementary material, figure S5.

Length in the transient phase.
Generation times in heterotrophs and cyanobacteria are very different. While the former are observed on the scale of hours, the latter grow on a scale of days/weeks. The difference could be, in principle, owing to the different trophic state of the bacteria. This hypothesis could be tested by means of a comparison between different generation times intraspecifically in the autotrophs or in the heterotrophs. However, this goes beyond the scope of the present work and we proceed relying on the mathematical derivation of §2.2. According to the relationship between generation time and turnover established in §2.2, we assume that the heterotrophs and cyanobacteria have high and low turnovers, respectively. The simulation results presented in §3.1.3 show that the peak value of the mean length in the transient phase increases with decreasing turnover. Figure 3 shows that heterotrophic filaments will rarely be longer than 60 cells. The length of cyanobacteria can be significantly higher than that of heterotrophs. Nostoc muscorum and A. variabilis filaments have up to around 350 and 600 cells, respectively. Hence, cultured bacterial populations  Life cycles in filamentous bacteria V. Rossetti et al. 1779 support the model predictions, indicating that bacteria with a putatively higher turnover (in this case heterotrophic species) are able to elongate less than those with lower turnover (cyanobacteria).

Distribution of lengths at the time when carrying
capacity is reached. In figure 3f, the box plots and the pie charts show the proportion of short and long filaments in the cultured bacterial populations when they reach their carrying capacity. As for the simulation results, the pie charts illustrate three length classes, namely single cells, double-celled filaments and filaments with more than two cells. The heterotrophs show a considerable amount of short filaments, with a proportion of single cells and double-celled filaments varying across species in the range 8-16% and 29-42%, respectively. In cyanobacteria, filaments with one or two cells are absent. Figure 3f can be compared with figure 1d. Again, the experimental results confirm the theoretical predictions, according to which short filaments at the carrying capacity are more abundant in populations with high turnover.

Serial transfer experiment and simulations.
If filament elongation is governed by density-dependent processes as hypothesized in the model in §2.1, then interventions on such underlying processes should produce testable hypotheses. Experimental manipulation of cell densities is one simple possibility. If the population drops significantly below its carrying capacity, the filaments that survive in this new condition should be able to elongate or at least maintain their original length until the carrying capacity is reached again. We simulated and experimentally carried out a series of successive transfers of bacterial populations, whereby a fraction of the population was transferred to a fresh medium whenever the average length of the original population was close to its peak. Filaments developing in the fresh medium were expected to keep or increase the length reached at the previous stage. Details of the computational and experimental settings are provided in the electronic supplementary material, §S7. Figure 5a shows the average length of filaments obtained by successive transfer simulations. Figure 5b shows the average length of F. aestuarina tracked along successive transfers. In this case, filament size When the average filament length is at its peak, the population density is below the carrying capacity. Whenever a culture reaches this peak, the corresponding average length is used as a starting condition of a new simulation at lower density (transfers [1][2][3][4]. This procedure simulates the transfer of an aliquot of a higher density population to a fresh medium. Edges of the error bars correspond to the 2.5th and 97.5th percentiles of the distribution of the mean filament lengths in the 1000 runs (black circles, initial curve; grey circles, transfer 1; black squares, transfer 2; grey squares, transfer 3; black triangles, transfer 4). (b) Average length of F. aestuarina during serial transfers to fresh medium. The first transfer was done after 17 h. For the first five transferred cultures, the mean length was measured and plotted at several successive time points. For subsequent cultures, only one measurement was taken. Edges of the error bars correspond to the 2.5th and 97.5th percentiles of the measured lengths. slightly increases or remains almost constant in most of the cases. The fact that filaments do not elongate so strongly as in the simulations could be owing to physical handling (breakage of long filaments while pipetting) or to other culture conditions influencing the elongation that were not considered in the model. The model results indicate that the filaments transferred to the fresh medium maintain or increase their length, and eventually, start shortening when they reach their carrying capacity again. The experimental results, although not as strong as the simulations, show that the transferred bacteria are generally able to maintain their length.

DISCUSSION
Filamentous bacteria can be viewed as a simple example of multicellular individuals. Some of them, such as filamentous cyanobacteria, have been present in the environment since the early evolution of life, as indicated by fossil records and phylogenetic analyses [70][71][72][73]. The theoretical model shows that at various life stages, two populations with the same fitness and carrying capacity can differ in the composition of filament length classes according to the life-history traits regulating their growth. Although the long-term stationary distribution of filament length is the same for any sufficiently small turnover rate, the differences owing to turnover observed at other growth stages (e.g. transient phase, when reaching carrying capacity) are quantitatively significant and biologically relevant. Filaments belonging to a population with high turnover (and hence short generation time) have moderate lengths at their maximum growth phase. Furthermore, when the population reaches carrying capacity, they are chiefly unicellular. Individuals of a population with low turnover (and hence long generation time) can instead become very long in the transient phase, and still be multicellular when reaching the carrying capacity. The species that we cultured show evidence that the predictions of the model hold. The relationship between bacterial generation time and turnover derived in §2.2 helps us to understand the relationship between simulations and experiments. The heterotrophs have a faster generation time than the cyanobacteria (and hence a higher turnover). As a consequence, while the heterotrophic filaments are rarely longer than 60 cells, the cyanobacteria can reach lengths of up to 600 cells. Both theory and experiments show that according to different life-history traits, bacteria can cover a wide spectrum of lengths, indicating that multicellularity can be achieved to different degrees. Moreover, a key feature of this model is that the predominance of multicellularity during the life cycle is not necessarily determined by morphology-dependent differences in fitness. Here, the factors affecting multicellularity are birth and death rates, and hence turnover rate and generation time.
The model indicates a common temporal behaviour for filament length distributions. The mean length increases and peaks during the transient growth phase, and then decreases until it reaches a value approaching the initial conditions. The results from cultured bacterial populations of heterotrophic and photoautotrophic species show that the predicted trend in filament length distribution is robust across species, and is representative of what occurs in nature. The fact that individuals eventually become shorter (and of average length 2 at the stationary distribution) regardless of the turnover is because of the effect of cell lysis on a filament. Consider a filament of length L. A cell division (birth) increases its length by one unit, namely the filament length becomes L þ 1. By contrast, after a cell lysis, the expected filament length will be on average L/2. The average effect of cell death on length is hence much larger than that of cell birth. This is consistent with the expectation that the average filament length in a population can keep increasing if the death rate remains sufficiently low, as shown in the estimation of the time of first filament split (electronic supplementary material, §8). When the average filament length in the population is at its stationary distribution, cell birth and death have the same effect. This is achievable only if the filaments are of average length 2. This argument is supported by the analysis of the analogous continuous time model derived in the electronic supplementary material. In this context, a way to avoid the breakage into small filaments is to decrease the effect of cell lysis on organism size. A solution in this direction is the evolution of a second or third dimension, as in spherical multicellular bacteria [74] or microbial biofilms. In this case, if we think of a square or a ball of N cells, a cell division will still lead to an individual of N þ 1 cells, but cell death will also only decrease size by one unit, resulting in a size of N 2 1 (as opposed to N/2).
There are, to our knowledge, no previous theoretical and experimental results that are directly comparable to the results presented here. However, the extensive work of Hahn & Hoefle [27] and Hahn et al. [45] set important precedents and reference points. These works focus on the effects of predators on selection for longer filament lengths. Their results are not directly comparable to those presented here because the death rates owing to predators preferentially afflict smaller filaments and single cells. However, they also provide revealing information on the dynamics of filament formation in the absence of predators. Hahn et al. [45] show that increasing the nutrient flow and the dilution rate of a chemostat in the absence of predators leads to longer filament lengths. This corresponds to our prediction that in transient growth conditions with lower population densities, filament lengths would get longer. Furthermore, although they do not present extensive data on this, [45] they mention that in batch culture, the percentage of multicellular filaments is higher in the exponential phase in comparison with the stationary phase. This situation is again consistent with the predictions presented here.
The model studied here is a general model of filamentous growth, and is not restricted to prokaryotes. Many simple forms of eukaryotic algae grow in a multicellular filamentous form, and the theoretical results presented here are in principle also valid for such eukaryotes. However, it is clear that eukaryotic algae have surpassed their prokaryotic counterparts in their Life cycles in filamentous bacteria V. Rossetti et al. 1781 morphological diversity. There are at least six major algal lineages in which filamentous multicellularity has evolved: the Phaeophyta, Rhodophyta, Chrysophyta, Chlorophyta, Charophyta and Embryophyta [75]. However, in each of the latter cases, higher dimensional variants have also evolved, whereupon cell division can occur in more than one planar dimension. This diversity can for example range from the simple globular multicellularity of some of the Volvocales (within the Cholorphyta), to the highly differentiated and developmentally complex vascular land plants (within the Embryophyta). Furthermore, plants can differentiate into diverse developmental modules such as branching nodes, leaves, inflorescences and root structures. By combining such modules in different combinations, they can access a wider variety of morphologies during their growth process [76,77]. This flexibility in morphology can translate into both phenotypic plasticity and evolvability.
Developmental life cycles are not unique to multicellular eukaryotes, and an understanding of how they evolved is lacking. Theory and experiments in this article indicate that the size of the individual, in terms of number of cells, can follow a cycle: from unicellular (or short length) to multicellular and back to unicellular (or short length) again. Each time, the initiation of the cycle can be hypothetically achieved by moving cell density below the carrying capacity (e.g. increases in fluid or nutrient availability). The transfer experiment and simulations do not contradict this hypothesis.
In our model, only the birth and death rates play a role in the change in size during the bacterial life cycle. These life-history traits are intrinsic properties of every living organism. This pattern, which automatically arises from the interplay between ecology and the filamentous nature of the bacteria, is an emergent property. As such, it can serve as the basis for a primitive life cycle, upon which a more complex developmental programme can be subsequently built. As a case in point, some species of Nostoc can differentiate into hormogonia during low growth phases characterized by high stress or scarce resources. These hormogonia are short filaments specialized for survival in harsh conditions, which subsequently migrate to a new area before growing into full-grown filaments again [78]. Our results suggest that the cycling of filament lengths between growth and static phases of the population is a pre-existing context within which multicellular developmental cycles and differentiation can evolve.