The evolution of division of labour in structured and unstructured groups
Abstract
Recent theory has overturned the assumption that accelerating returns from individual specialisation are required to favour the evolution of division of labour. Yanni et al., 2020, showed that topologically constrained groups, where cells cooperate with only direct neighbours such as for filaments or branching growths, can evolve a reproductive division of labour even with diminishing returns from individual specialisation. We develop a conceptual framework and specific models to investigate the factors that can favour the initial evolution of reproductive division of labour. We find that selection for division of labour in topologically constrained groups: (1) is not a single mechanism to favour division of labour—depending upon details of the group structure, division of labour can be favoured for different reasons; (2) always involves an efficiency benefit at the level of group fitness; and (3) requires a mechanism of coordination to determine which individuals perform which tasks. Given that such coordination must evolve prior to or concurrently with division of labour, this could limit the extent to which topological constraints favoured the initial evolution of division of labour. We conclude by suggesting experimental designs that could determine why division of labour is favoured in the natural world.
Editor's evaluation
This manuscript presents a theoretical study of the evolution of division of labor, exploring the impact of topology, the convexity and concavity of fitness returns on investment, and different biological modes through which division of labor may arise. This is a difficult topic to study as division of labor evolved long ago, and many theoretical predictions have proven difficult to directly test. The results presented here may provide the next step necessary to produce truly testable hypotheses on how division of labor evolves, and will be of interest to evolutionary biologists, mathematical biologists, and biophysicists.
https://doi.org/10.7554/eLife.71968.sa0Introduction
Division of labour, where cooperating individuals specialise to carry out distinct tasks, plays a key role at all levels of biology (Bourke, 2011; Queller, 1997; Maynard Smith and Szathmáry, 1995; West et al., 2015). Cells are built by genes carrying out different functions (Bourke, 2011; Levin and West, 2017). In clonal groups of bacteria, cells specialise to produce and secrete different factors that facilitate growth (Dragoš et al., 2018a; Veening et al., 2008; West and Cooper, 2016). Pathogens rely on division of labour for protection from the host immune response and competitors (Ackermann et al., 2008; Diard et al., 2013). Multicellular organisms are composed of reproductive germ cells and sterile somatic cells that are not passed to the next generation (Bourke, 2011; Maynard Smith and Szathmáry, 1995). The ecological dominance of the social insects arises from division of labour between queens and the different types of workers (castes) (Hölldobler and Wilson, 1990; Oster and Wilson, 1978).
It has long been established that the evolution of division of labour requires an efficiency benefit from individual specialisation (Figure 1A and B for reproductive division of labour) (Bourke, 2011; Cooper and West, 2018; Ispolatov et al., 2012; Michod, 2006; Oster and Wilson, 1978; Schiessl et al., 2019; Biggart, 1776; Maynard Smith and Szathmáry, 1995; Solari et al., 2013). In particular, that there is an accelerating (convex) return when individuals commit more effort to a particular task, such that twice the investment more than doubles the return (Bourke, 2011; Cooper and West, 2018; Ispolatov et al., 2012; Michod, 2006; Solari et al., 2013). An accelerating return from individual investment can exist for several reasons. A task could become more effective as more effort is put into it, or it could be carried out with diminishing costs. This could occur if there are large upfront costs from performing a task. For instance, any reproduction by a cell in Volvocine groups first requires individual growth to the size of a daughter colony (Michod, 2006). Alternatively, there could be a disruptive cost to carrying out multiple tasks at the same time if the tasks do not mix well. For instance, in cyanobacteria the enzymes that fix environmental nitrogen are degraded by oxygen, a biproduct of photosynthesis (Flores and Herrero, 2010).
In contrast, Yanni et al. found that division of labour between helpers and reproductives can sometimes be favoured even when there are diminishing (concave) returns from individual specialisation (Yanni et al., 2020). Specifically, reproductive division of labour could arise in topologically constrained groups—where each cell in a spatially structured group shares cooperative benefits with only their direct neighbours (Staps and Tarnita, 2020; Yanni et al., 2020). Their analyses suggested that this is particularly likely to occur in sparsely structured groups, where cells have a small number of neighbours (Yanni et al., 2020). This is a novel result. Diminishing returns means that specialised individuals are inefficient, and earlier work suggested that division of labour could not be favoured in this situation (Figure 1A, Cooper and West, 2018; Michod, 2006; Schiessl et al., 2019). Consequently, this result has the potential to overturn our understanding of the factors that favour the evolution of division of labour.
However, there are several issues that still need to be resolved with how topological constraints can favour division of labour. Why exactly do the predictions of this new theory differ from previous theory? Are special group structures the only way to alter the predictions of the previous theory, or is this an example of a more general phenomenon (Rueffler et al., 2012)? Do these findings rely upon implicit assumptions, which may not be reasonable during the initial evolution of division of labour? Answering these questions is not only of theoretical importance: it is also key for planning future empirical studies. Quantifying the shape of the returns from individual specialisation has been assumed to be a fundamental step in determining why division of labour was favoured in some species, but not others (Diard et al., 2013; Dragoš et al., 2018a; Flores and Herrero, 2010; Koufopanou, 1994; Mridha and Kummerli, 2021; Strassmann et al., 2000; Veening et al., 2008).
We first use the methodology developed by Rueffler et al., 2012, to derive the general conditions that favour the initial evolution of reproductive division of labour between helpers and reproductives. We then use this framework to examine when and why topological constraints can favour division of labour. More specifically, we determine the ultimate cause of division of labour in specific topologically constrained groups, such as filaments and branching growths, as well as in a general analysis of arbitrary group structures. We then ask whether division of labour without an accelerating return from individual specialisation could arise in groups without topological constraints. To test our hypothesis that betweenindividual coordination is required for division of labour in these cases, we reexamine our models while assuming that cells adopt helper and reproductive roles randomly (no coordination). We finish by suggesting experimental designs for determining why division of labour has evolved in specific species.
Results and discussion
General invasion analysis
We follow previous studies by assuming that individual fitness is the product of individual viability, which is the chance of surviving to maturity, and individual fecundity, which is proportional to the number of offspring if the individual reaches maturity (Cooper and West, 2018; Michod, 2006; Yanni et al., 2020). We examine the specific case of reproductive division of labour between helpers and reproductives, where helpers are more cooperative, contributing to a higher viability for group members, and reproductives are less cooperative, contributing to higher individual fecundity.
We consider an initial population of clonal groups each containing $n$ individuals, in which all individuals cooperate at the evolutionarily stable (ES) level ($z}^{\ast$), which is the level that cannot be outcompeted by a mutant strain that uses a different level of uniform cooperation across the group (Maynard Smith, 1982). We then ask when this population of uniform cooperators can be invaded by a mutant strain that employs a reproductive division of labour.
Without loss of generality, we assume that the mutant strain is composed of $n}_{h$ helpers that invest ${z}_{h}\ge 0$ into cooperation and $n}_{r$ reproductives that invest ${z}_{r}\ge 0$ into cooperation (where ${n}_{h}+{n}_{r}=n>2$ and $z}_{h}\text{}{z}_{r$). We set individual fitness as the product of individual fecundity, $F>0$, and individual viability, where helpers and reproductives may in principle have different viability functions, ${V}_{h}>0$ and ${V}_{r}>0$ (but see Appendix C.2) (Michod, 2006; Yanni et al., 2020). The fitness of the clonal group is given by the sum of individual fitness:
where the first term on the righthand side is the total fitness of the prospective helpers and the second term on the righthand side is the total fitness of the prospective reproductives.
Fecundity is determined by an individual’s investment in cooperation ($F=F\left(z\right),$ where $z$ is the focal individual’s level of cooperation) and viability is determined by the level of cooperation at the level of the group (${V}_{h}={V}_{h}\left({z}_{h},{z}_{r}\right)$ and ${V}_{r}={V}_{r}\left({z}_{h},{z}_{r}\right)$). We assume that there is a tradeoff between fecundity and viability such that higher individual cooperation leads to lower individual fecundity (${F}^{\prime}\left(z\right)<0$), but that more cooperation leads to a higher viability for all individuals (i.e. ${V}_{h}^{{z}_{h}},{V}_{h}^{{z}_{r}},\text{}{V}_{r}^{{z}_{h}},\text{}{V}_{r}^{{z}_{r}}0$, where superscripts denote partial derivatives). We assume that viability selection occurs just prior to reproduction. This is consistent with previous models and ensures that there there is no feedback between a cell’s viability and its ability to produce cooperative benefits for the group (Cooper and West, 2018; Michod, 2006; Yanni et al., 2020).
We determine the invadability conditions that favour reproductive division of labour by applying the general approach of Rueffler et al., 2012. The key step is to approximate the relative fitness of a reproductive division of labour mutant by taking a secondorder Taylor expansion of fitness, centred on the resident strategy of uniform cooperation, $z}^{\ast$:
where $\mathrm{\Delta}{z}_{h}>\mathrm{\Delta}{z}_{r}$ captures the change in the level of cooperation for mutant helpers and reproductives, respectively, which we assume are small in magnitude. The superscripts represent first and secondorder partial derivatives, where all partial derivatives are evaluated at the resident strategy of uniform cooperation ($z}_{h}={z}_{r}={z}^{\ast$). If a mutant strain exists such that Equation 2 is positive ($W\left({z}^{\ast}+\mathrm{\Delta}{z}_{h},\text{}{z}^{\ast}+\mathrm{\Delta}{z}_{r}\right)W\left({z}^{\ast},\text{}{z}^{\ast}\right)$), then division of labour between helpers and reproductives is favoured to evolve. Conversely, if for all possible mutant strains, Equation 2 is negative ($W\left({z}^{\ast},\text{}{z}^{\ast}\right)W\left({z}^{\ast}+\mathrm{\Delta}{z}_{h},\text{}{z}^{\ast}+\mathrm{\Delta}{z}_{r}\right)$), then uniform cooperation is evolutionarily stable.
The three pathways to division of labour
We found that reproductive division of labour could be favoured for three distinct reasons, corresponding to different subsets of terms on the righthand side of Equation 2. Our results for reproductive division of labour, where fitness is partitioned as the product of fecundity and viability, align with those found by Rueffler et al., 2012, for division of labour more generally. We now go through these three distinct scenarios.
Scenario 1: Accelerating returns from individual specialisation
The first and most studied scenario that can favour division of labour is when there are accelerating returns from individual specialisation. This occurs if there is an accelerating fitness return from either helper specialisation in cooperation or reproductive specialisation in fecundity (Figure 1; Cooper and West, 2018; Michod, 2006; Oster and Wilson, 1978).
Mathematically, this scenario is a consequence of the third and fourth terms of the Taylor expansion ($2c$ and $2d$), which capture the secondorder fitness effect of a small, unilateral change in cooperation by either prospective helpers or reproductives, respectively. Division of labour is favoured to evolve whenever at least one of $2c$ and $2d$ is greater than zero (${W}^{{z}_{h}{z}_{h}}>0$ or ${W}^{{z}_{r}{z}_{r}}\text{}0$), and where we have assumed that the first two terms are both zero (${W}^{{z}_{h}}=0$ and ${W}^{{z}_{r}}=0$; see betweenindividual differences below; Figure 2A). In either scenario, an efficiency benefit to group fitness arises from individual specialisation because the more effort that each individual puts into a task, the better they can perform that task. Rueffler et al. termed these kinds of scenario as ‘accelerating performance functions’ (Rueffler et al., 2012).
Scenario 2: Betweenindividual differences
The second scenario that can favour reproductive division of labour is when there are preexisting differences between individuals in the group, such that some individuals are predisposed to one task or the other. For example, if some individuals can secure larger viability benefits for the group at the same fecundity cost as others (Figure 2B).
This scenario is captured by the first two terms of the Taylor expansion $\left(2a\text{}\text{and}\text{}2b\right)$ , which are the firstorder fitness effects from a small, unilateral change in the level of cooperation by prospective helpers or reproductives. If the direct fitness effects are nonzero (positive or negative) at the resident strategy of uniform cooperation (${W}^{{z}_{h}}\ne 0$ or ${W}^{{z}_{r}}\ne 0$), then division of labour can invade independently of any higherorder effects (the remaining terms in Equation 2).
We term this scenario ‘betweenindividual differences’ because it requires that there is preexisting phenotypic or environmental variation between individuals in the group. For the withinspecies case, ancestral groups are usually composed of clonal or highly related individuals, who will be phenotypically similar or identical. Consequently, this mechanism could be less important for the division of labour except when there are consistent differences in the microenvironment experienced by different individuals (Tverskoi et al., 2018; Tverskoi and Gavrilets, 2021). In contrast, this scenario is likely to be widespread in the evolution of nonreproductive division of labour between species, such as for mutualisms or symbioses (Kiers et al., 2011; Rueffler et al., 2012; Wyatt et al., 2014). Individuals of different species often differ in their abilities to perform certain tasks (Kiers et al., 2011). Rueffler et al. termed this scenario ‘positional’, but we avoid that term to prevent confusion with topological position (Rueffler et al., 2012).
Betweenindividual differences provide a firstorder fitness benefit to dividing labour, and so it does not matter whether the subsequent benefits of increased cooperation or fecundity are accelerating or diminishing (Figure 1), so long as these benefits are different for different individuals (Figure 2). When some individuals are predisposed to being either helpers or reproductives, then individual specialisation provides an efficiency benefit to group fitness by capitalising on these inherent differences.
Scenario 3: Reciprocal specialisation
The final scenario that can favour division of labour is when reciprocal specialisation by both helpers and reproductives provides a fitness benefit to the group (Figure 3). This scenario requires two key conditions. First, simultaneous specialisation, where some individuals invest more in cooperation (more viability benefits for the group), and others invest less in cooperation (greater individual fecundity; but see below). Second, this reciprocal specialisation must provide a grouplevel fitness benefit, because the increased benefits of cooperation are preferentially directed towards reproductives.
Mathematically, this scenario involves the last term of the Taylor expansion ($2e;\text{}{W}^{{z}_{h}{z}_{r}}\mathrm{\Delta}{z}_{h}\mathrm{\Delta}{z}_{r}$). This term is generated by a betweenindividual, secondorder fitness effect, capturing how increased investment in viability by some individuals affects the returns from increased investment in fecundity by others, and vice versa. Rueffler et al. referred to this as a ‘synergistic benefit’ to division of labour (Queller, 1985; Queller, 2011; Rueffler et al., 2012).
Critically, this scenario still involves an efficiency benefit to specialisation, but at the level of group fitness rather than in each fitness component separately (Appendix C.1). By this we mean that there is an accelerating fitness benefit to the group when helpers and reproductives reciprocally specialise, leading to a higher group fitness than in groups with uniform cooperation (generalists). This occurs if the increased help given to reproductives is sufficiently amplified by the increased fecundity of reproductives (Yanni et al., 2020). This synergistic efficiency benefit can favour division of labour even if there are diminishing returns from individual specialisation.
Division of labour by reciprocal specialisation can also evolve without a joint mutation in the level of cooperation of both helpers and reproductives (no simultaneous specialisation). In this case, the chance invasion (to fixation) of a slightly deleterious mutant that specialises in only one phenotype can destabilise uniform cooperation, creating a selection pressure for the other phenotype to also specialise that is greater than the selection pressure to purge the initial mutant. In this twostep scenario, it is nevertheless the synergistic benefit from reciprocal specialisation that makes division of labour more efficient.
Group structure in the general framework
Our above analysis has shown that reproductive division of labour can be favoured for three reasons: (1) accelerating returns make individual specialisation more efficient; (2) betweenindividual differences make individual specialisation more efficient; and (3) there is a synergistic efficiency benefit from reciprocal specialisation. These results agree with previous analyses by Rueffler et al., 2012.
We now use this framework to examine how and why topological constraints can favour division of labour in the absence of an accelerating return from individual specialisation (i.e. when scenario 1 does not hold). We ask three questions. First, can topological constraints favour division of labour by betweenindividual differences (scenario 2), and/or by reciprocal specialisation (scenario 3)? Second, are topological constraints the only way to evolve a division of labour without an accelerating return from individual specialisation? Third, does the evolution of division of labour by betweenindividual differences (scenario 2) and reciprocal specialisation (scenario 3) require coordination between individuals to determine which cells become helpers or reproductives?
Question 1: How do topological constraints favour division of labour?
We consider two spatial models, based on the group structures proposed by Yanni et al., to examine whether topologically constrained groups favour division of labour by: (a) betweenindividual differences and/or (b) reciprocal specialisation (Yanni et al., 2020).
Can topological constrains lead to division of labour by betweenindividual differences?
Consider a group in which cells alternately have either two or three neighbours, in a branching structure (Figure 4). Such a group structure might have occurred for some early forms of multicellular life (Yanni et al., 2020). We term cells with three neighbours ‘node’ cells and cells with two neighbours ‘edge’ cells. We assume that cells investing an amount $z\ge 0$ into cooperation produce an amount $H\left(z\right)$ of a public good. We assume nonaccelerating returns from individual specialisation (i.e. ${H}^{\prime \prime}\left(z\right)\le 0$ or $F\prime \prime \text{}\left(z\right)\le 0$). The cell keeps a fraction $1\lambda$ of the public good that it produces, and the remaining fraction $\lambda$ is shared equally between its direct neighbours (the ‘shareability’ of cooperation: $0<\lambda \le 1$). We assume that the viability of a cell is equal to the sum of the public good that it absorbs.
For this model, we find that for all social traits ($\lambda >0$), reproductive division of labour by betweenindividual differences can evolve (Figure 4A; Appendix A.1). This occurs because different cells have different viabilityfecundity tradeoffs depending on their position in the group. Edge cells receive relatively less public good from their (fewer) neighbours, and so pay a smaller opportunity cost from decreased fecundity (increased cooperation). In contrast, node cells receive relatively more public good from their (more numerous) neighbours, and so pay a larger opportunity cost from decreased fecundity (increased cooperation). Consequently, this betweencell difference favours node cells to specialise in fecundity (reproductives) and edge cells to specialise in increased cooperation (helpers). Importantly, because this pathway to division of labour is driven entirely by a firstorder effect (2a and 2b), it does not require a secondorder efficiency benefit from specialisation (2c, 2d, or 2e).
More generally, a formal analysis of arbitrary group structures reveals that division of labour by betweenindividual differences can always evolve whenever the number of neighbours varies for different cells in the group (Appendix A.2).
Can topological constraints lead to division of labour by reciprocal specialisation?
Consider a onedimensional chain of cells, as examined by Yanni et al. Chains are found in species like cyanobacteria that form filaments of cells, and such structures might have been important at the onset of the evolution of multicellularity (Figure 4, Yanni et al., 2020). We assume arbitrarily that ‘odd’ cells along the filament are putative helpers and ‘even’ cells are putative reproductives. We otherwise make the same assumptions as for the branching structure model: there is a nonaccelerating return from individual specialisation (i.e. ${H}^{\prime \prime}\left(z\right)\le 0$ or $F\prime \prime \text{}\left(z\right)\le 0)$, and the cell keeps a fraction $1\lambda$ of the public good that it produces, with the remaining fraction $\lambda$ being shared equally by its direct neighbours.
If the amount of public good shared with neighbours is sufficiently large (high $\lambda$), then we find that division of labour via reciprocal specialisation can evolve (Figure 4B; Appendix A. 3). For instance, in the case of linear fecundity and public good returns (${H}^{\prime \prime}\left(z\right)={F}^{\prime \prime}\left(z\right)=0$), division of labour by reciprocal specialisation can evolve if helpers share more of the public good that they produce with their neighbours than they keep for themselves ($\lambda >\frac{1}{2}$). If there are diminishing returns from specialisation (${H}^{\prime \prime}\left(z\right)<0$ or ${F}^{\prime \prime}\left(z\right)<0$), then division of labour can still be favoured but then the amount of the public good preferentially shared with neighbours must be even greater still (higher $\lambda ;$ Figure 4E).
For an arbitrary group structure, our analysis in the previous section implies that division of labour can evolve by betweenindividual differences, unless every cell in the group has the same number of neighbours. Consider a group in which every cell has exactly d neighbours. In this case, we show (Appendix A.5) that division of labour can still be favoured due to reciprocal specialisation if:
λ is the shareability of cooperation as defined previously and μ captures how easily the group can be ‘bipartitioned’. That is, μ is a measure of the extent that the group can be divided into two classes of cells such that cells are neighbours with many cells of the opposing class and few neighbours of their own class. Thus, reciprocal specialisation can favour division of labour if: (1) groups are more sparse (low d); (2) groups are structured such that helpers can be neighbours with reproductives more than with other helpers, and vice versa (high μ); and/or (3) when the benefits of cooperation are preferentially shared with neighbours (high $\lambda$). In combination, these three factors amplify the synergistic benefits of reciprocal helper and reproductive specialisation, which can produce an accelerating fitness return for the group, even when there are nonaccelerating returns from individual specialisation (Appendix C.1).
If one or two of these factors are particularly favourable for reciprocal specialisation, then the condition(s) on the remaining factor(s) can be relaxed. For instance, cells in a filament have only two neighbours ($d=2$), and the potential alternation of helpers and reproductives in the filament means that helpers can share their cooperative public goods with reproductives exclusively (maximal $\mu ).$ Consequently, reciprocal specialisation is possible even when the shareability of cooperation is reasonably low (e.g. $\lambda >\frac{1}{2}$ for linear benefits).
To conclude, topological constraints are not a single explanation for division of labour, in that they can favour division of labour for two different biological reasons. Different group structures can lead to either betweenindividual differences favouring division of labour (scenario 2) or reciprocal specialisation favouring division of labour (scenario 3). In all cases, there is an efficiency benefit from specialisation at the level of group fitness even if the returns from individual specialisation are nonaccelerating.
Question 2: Are topological constraints required for division of labour without accelerating returns from individual specialisation?
We considered a wellmixed social group of $n$ cells, where all cells share the benefits of cooperation with one another, and so there are no topological constraints (Figure 4). We then examined whether division of labour could be favoured by: (a) betweenindividual differences; and/or (b) reciprocal specialisation. In both cases, we assume that when a cell invests $z$ into cooperation, it produces an amount $H\left(z\right)$ of a public good. A cell keeps a fraction $1\lambda$ of the public good that it produces and the remaining fraction $\lambda$ is shared by the rest of the social group members equally. We again consider the case where there is a nonaccelerating return from individual specialisation ($H{\left(z\right)}^{\prime \prime}\le 0$; $F\prime \prime \text{}\left(z\right)\le 0)$.
Can betweenindividual differences favour division of labour without a topological constraint?
In the wellmixed group of identical cells, we find that division of labour cannot arise by betweenindividual differences (Appendix A.4). This is because all cells have the same number of neighbours, which we have shown more generally can never produce betweenindividual differences. This prediction could be violated if one of our assumptions do not hold: for instance, if there are consistent differences in the microenvironment that predispose some cells to one task or the other (Tverskoi et al., 2018; Tverskoi and Gavrilets, 2021; Yanni et al., 2020).
Can reciprocal specialisation favour division of labour without a topological constraint?
In the wellmixed group of identical cells, if the amount of public good shared with neighbours is sufficiently large (high $\lambda$), then we find that division of labour via reciprocal specialisation can evolve (Figure 4F; Appendix A.4). If there are linear returns from increased specialisation (${H}^{\prime \prime}\left(z\right)={F}^{\prime \prime}\left(z\right)=0$), then division of labour can evolve when the public good produced by an individual benefits an average group member more than the producer ($\lambda >\frac{n1}{n}$; Figure 4C). These results are like those found for a filament of cells (Figure 4B). In both cases, more generous sharing (higher $\lambda$) means that the synergistic benefits of reciprocal specialisation can be great enough to compensate for the nonaccelerating returns from individual specialisation. In wellmixed groups, very generous sharing ($\lambda \approx 1$) also compensates for the fact that helpers are neighbours with all other helpers (no sparsity and minimally ‘bipartionable’).
To conclude, the wellmixed model shows that a topological constraint is not required for the evolution of division of labour with nonaccelerating returns from individual specialisation. This result is in direct contradiction to that of Yanni et al., 2020. This difference arises because helpers in their model always benefit at least as much as any of its neighbours from its own public good production ($\lambda \le \frac{n1}{n}$) (Yanni et al., 2020). Our model allows for biological scenarios where the public good benefits an average neighbour more than the producer ($\lambda >\frac{n1}{n}$). For instance, reproductive cells in cyanobacteria may absorb more of the fixed nitrogen produced by helpers than helpers do to meet the large energetic requirements of cell duplication and division (Flores and Herrero, 2010; Herrero et al., 2016; Meeks and Elhai, 2002). At the extreme, the public good can be an ‘othersonly’ trait that benefits neighbours but not the producer at all ($\lambda =1$) (Pepper, 2000). An example of this are the dispersal benefits provided by stalk cells in Dyctiostelium discodeum fruiting bodies or the selfsacrificing behaviour of helper cells in Salmonella enterica infections (Ackermann et al., 2008; Strmecki et al., 2005). Consequently, our model allows for a wider spectrum of biologically realistic scenarios. Critically, division of labour can be favoured in a group of wellmixed cells because it provides an efficiency benefit at the group level, via reciprocal specialisation (scenario 3), in an analogous way to our model with a onedimensional chain of cells (question 1b).
Question 3: Is coordination required to favour division of labour without accelerating returns from individual specialisation?
We hypothesised that the benefits of betweenindividual differences (scenario 2) or reciprocal specialisation (scenario 3) rely on the implicit assumption that cells are coordinating which individuals specialise to become reproductive and helpers. This matters because mechanisms for coordinating division of labour, such as between cell signalling, might not be expected to exist before division of labour has evolved (Cooper et al., 2022; Liu et al., 2021). Consequently, if coordination was required, then this could limit the extent to which topological constrains favour the initial evolution of division of labour.
We investigated this hypothesis by repeating our above analyses, while assuming that cells do not have access to information that allows them to coordinate their phenotypes. Specifically, cells do not know if they are ‘odd’ or ‘even’, or if they are ‘edge’ or ‘node’. We assumed instead that a reproductive division of labour mutant induces each cell in the group to adopt the role of a helper or reproductive with a uniform probability (random specialisation). Random specialisation has been observed in a number of microbes (Ackermann et al., 2008; Diard et al., 2013; Veening et al., 2008). For filaments, branching group structures, and wellmixed groups, we found that division of labour can no longer evolve with nonaccelerating returns from individual specialisation (Figure 4D–E; Appendices B.1 and B.2). In Appendix B.3, we have shown that this result holds for any group structure.
Consequently, for division of labour to evolve with nonaccelerating returns from individual specialisation, there must exist some mechanism to coordinate which cells specialise to perform which tasks. It is possible that the mechanism need not produce a perfect allocation of labour across the group as analysed in our models (Liu et al., 2021). However, because division of labour cannot be favoured to evolve if role allocation is fully random, at least some degree of even imperfect betweencell coordination will be required.
A clear example of coordinated division of labour in topologically constrained groups is the use of betweencell signalling in some cyanobacteria filaments to determine which cells become sterile nitrogen fixing heterocysts and which cells become reproductive photosynthesisers (Flores and Herrero, 2010; Meeks and Elhai, 2002). However, a signal to coordinate distinct phenotypes must exist prior to or concurrently with the emergence of division of labour, and so a topological constraint is less likely to have favoured the initial evolution of division of labour in cyanobacteria. Alternatively, division of labour could have been favoured by an accelerating return from individual specialisation (scenario 1), with coordination only being favoured to evolve subsequently. Empirically, an accelerating return seems likely, as the key tasks performed by reproductives and helpers do not mix well (photosynthesis and nitrogen fixation) (Flores and Herrero, 2010; Meeks and Elhai, 2002).
These analyses do not suggest that topological constraints could never favour the initial evolution of division of labour. For instance, a preexisting cue could allow division of labour to initially evolve with a metabolically cheaper form of coordination. More specifically, phenotype could be determined in response to the number of neighbours or the local concentration of some resource. Further, if there are preexisting differences between individuals due to a preexisting mechanism of coordination, then this mechanism can be coopted to coordinate division of labour. However, the biological plausibility of any preexisting mechanism would need to be explicitly justified and modelled on a casebycase basis (Duarte et al., 2011). This would include modelling the metabolic cost, benefits and effectiveness of the mechanism (Cooper et al., 2022; Duarte et al., 2012; Liu et al., 2021). Empirically, while betweencell coordination has evolved in several labourdividing microbial species, further studies—such as ancestralstate reconstructions—are needed to show whether coordination evolved prior to, concurrently with, or subsequent to division of labour in individual species.
In contrast, an accelerating return from individual specialisation depends on nonadaptive factors such as the physics, chemistry, or external constraints associated with the public good and its production. For instance, an accelerating return can arise if some intermediate products associated with cooperation and fecundity do not ‘mixwell’ on a chemical level. Consequently, no additional adaptive or preadaptation argument is needed to explain this pathway to division of labour.
Distinguishing the ultimate causes of division of labour in the wild
How can we distinguish empirically which of the different scenarios favoured realworld examples of division of labour (Figure 5)? We suggest experimental designs for microbial systems, where mixtures of helper and reproductive cells are grown together, and which make use of methods to genetically manipulate and measure the relative levels of cooperation and reproduction of each phenotype (Ackermann et al., 2008; Diard et al., 2013; Dragoš et al., 2018a; Dragoš et al., 2018b; Mavridou et al., 2018; Mridha and Kummerli, 2021; van Gestel et al., 2015). These are rough suggestions for the kind of experiments required, as details and possibilities will vary system from system, depending upon factors such as the degree of specialisation, the mechanism by which labour is divided, and what manipulations are possible. In addition, these experiments would need to follow from key first steps, such as demonstrating division of labour and a tradeoff between reproduction and cooperation (Diard et al., 2013; Dragoš et al., 2018a; Dragoš et al., 2018b; Veening et al., 2008; Zhang et al., 2020; Figure 5).
Testing for accelerating returns from individual specialisation (scenario 1): In at least three treatments, vary the level of cooperation performed by the helpers (Figure 5A topleft), to test whether the benefits of increased cooperation are accelerating (Figure 5A right). Across at least three other treatments, vary the level of reproduction by the reproductives (Figure 5A bottomleft), to test whether the benefits of increased fecundity are accelerating (Figure 5A right). At least three treatments are required to be able to test for nonlinear (accelerating) benefits.
Testing for betweenindividual differences (scenario 2): In some treatments, vary the level of cooperation of helpers (Figure 5B topleft). In other treatments vary the level of cooperation of reproductives (Figure 5B bottomleft). If division of labour evolves in this system by betweenindividual differences, then we should observe that group fitness varies differently depending on whether it is helpers or reproductives that are cooperating at a higher rate (Figure 5B right).
Testing for reciprocal specialisation (scenario 3): Use a classic $N\times N$ factorial experiment where the level of cooperation performed by helpers and the level of reproduction performed by reproductives are both varied (Figure 5C left). Division of labour is favoured by reciprocal specialisation if there is a significant interaction term between these two factors, with at least one treatment that produces a group fitness larger than that for uniformly cooperating cells (Figure 5C right).
Conclusion
Division of labour can be favoured to evolve without accelerating returns from individual specialisation. Nevertheless, for this to occur requires: (a) betweenindividual differences in taskefficiency or synergistic benefits from reciprocal specialisation and (b) a mechanism to coordinate which individuals perform which tasks. In contrast, accelerating returns can favour division of labour without a mechanism to coordinate task allocation, possibly making it more likely to favour the initial evolution of division of labour. Ultimately, determining the relative importance of these different pathways to division of labour is an empirical question, requiring experimental studies of the type we have outlined above.
Materials and methods
Resident strategy of uniform cooperation
We start by solving for the ESS strategy where both types of individuals invest the same amount in cooperation (${z}_{h}={z}_{r}=z$; uniform cooperation). This is the level of uniform investment in cooperation, ${z}^{*}$ , for which there is no selection for a uniform change in the amount of cooperation by all individuals in the group:
More explicitly, we can write this as:
where we have suppressed the functional dependencies for ease of presentation. The first term gives the group fitness change due to a marginal increase in ‘helper’ cooperation, and the second term gives the group fitness change due to a marginal increase in ‘reproductive’ cooperation. So, at the uniform strategy, ${z}^{*}$ , any increase in the fitness caused by the increased cooperation of one subgroup of the population is balanced by a commensurate decrease in fitness caused by the same increase in cooperation for the other subgroup $\frac{\mathrm{\partial}W\left({z}_{h},{z}_{r}\right)}{\mathrm{\partial}\text{}{z}_{h}}{}_{{z}_{h}={z}_{r}={z}^{\ast}}=\frac{\mathrm{\partial}W\left({z}_{h},{z}_{r}\right)}{\mathrm{\partial}\text{}{z}_{r}}{}_{{z}_{h}={z}_{r}={z}^{\ast}}$.
The constrained optimum, $\left({z}_{h},{z}_{r}\right)=\left({z}^{*},{z}^{*}\right),$ computed using Equation 4 does not necessarily correspond to a critical point of $W$, that is, the first derivates $\frac{\partial W\left({z}_{h},{z}_{r}\right)}{\partial {z}_{h}}$ and $\frac{\partial W\left({z}_{h},{z}_{r}\right)}{\partial {z}_{r}}$ might not vanish at this point. However, if the functions ${V}_{h}$ and ${V}_{r}$ are symmetrical, in the sense that
then the fitness function satisfies $W({z}_{h},{z}_{r})=W({z}_{r},{z}_{h})$, and it can be seen that this implies that the point $({z}_{h},{z}_{r})=({z}^{*},{z}^{*})$ is actually a critical point of $W({z}_{h},{z}_{r})$.
Division of labour by betweenindividual differences
Request a detailed protocolDivision of labour by betweenindividual differences occurs if either of the first two terms of the Taylor expansion (Equation 2) is nonzero (${W}^{{z}_{h}}\ne 0$ or ${W}^{{z}_{r}}\ne 0$). In this case, directional selection will increase or decrease the level of cooperation of one of the individual types. We give here the associated partial differentials of fitness (Equation 1):
These expressions capture the fitness consequences of a marginal increase in cooperation by helpers and reproductives, respectively. The first term of each captures the fecundity cost to own type of producing more public good, whereas the second term and third term are the viability benefits that accrue to both types from this increased cooperation. If directional selection in both traits is zero (${\frac{\partial W\left({z}_{h},{z}_{r}\right)}{\partial {z}_{h}}}_{{z}_{h}={z}_{r}={z}^{*}}=0$ and ${\frac{\partial W\left({z}_{h},{z}_{r}\right)}{\partial {z}_{r}}}_{{z}_{h}={z}_{r}={z}^{*}}=0$), then ${z}_{h}={z}_{r}={z}^{*}$ is a critical point, and Equations 6 and 7 imply that
These equations mean that, if ${z}_{h}={z}_{r}={z}^{*}$ is a critical point, then any marginal viability benefit to the group of increased cooperation by one subgroup is cancelled by the fecundity cost to that same subgroup. Moreover, Equations 8 and 9 together imply that
If this equation does not hold, then ${z}_{h}={z}_{r}={z}^{*}$ is not a critical point, that is, there is a difference in the viabilityfecundity tradeoffs between subgroups such that some individuals (without loss of generality, helpers) can secure larger benefits for the group at the same fecundity cost as others (reproductives). This gives our first condition for division of labour being able to evolve:
The betweenindividual differences condition for division of labour
If individuals are indistinguishable when both types invest equally in cooperation (${z}_{h}={z}_{r}=z$), then the viability functions satisfy ${V}_{h}(z,z)={V}_{r}(z,z)$. In this case, Condition 11 can be restated as:
This says that the contribution to total viability from the increased specialisation of helper individuals is strictly larger than the contribution to total viability from the increased specialisation of reproductives. As a result, helpers are predisposed to become more helperlike as they can gain larger viability gains for the group than the other type of individual.
Division of labour by an accelerating return from individual specialisation
Request a detailed protocolDivision of labour by an accelerating return from individual specialisation can occur if either of the third or fourth terms in the Taylor expansion (Equation 2) are positive in value (${W}^{{z}_{h}{z}_{h}}>0$ or ${W}^{{z}_{r}{z}_{r}}>0$). Taking the second derivative of fitness (Equation 24) with respect to each trait and evaluating at the critical point of uniform cooperation ${z}^{*}$ gives:
The terms of Equations 13 and 14 capture the secondorder effects of increased investment in cooperation. The first term of each captures the decline in the fitness benefit of increased cooperation due to the crossinteraction between fecundity and viability. For instance, as a helper invests more in cooperation (higher $x$), it increases its own viability (higher ${V}_{h}$), but its fecundity declines as well (lower $F$) and so the relative benefit of this increased viability is lessened (the cross term ${F}^{`}{V}_{h}^{{z}_{h}}$ is negative). This represents a kind of decelerating return from cooperation. The second term of each captures the secondorder effect of decreased investment in fecundity. If this term is positive, then this means that there is a diminishing fecundity cost to increased investment in cooperation, which can favour division of labour. The third and fourth terms capture the secondorder effect of increased investment in viability, that is, does each successive investment in the public good lead to a larger or smaller increase in viability than the previous investment of the same size? The return on investment (ROI) in viability is accelerating if ${V}_{h}^{{z}_{h}{z}_{h}}>0$, ${V}_{r}^{{z}_{h}{z}_{h}}>0$, ${V}_{h}^{{z}_{r}{z}_{r}}>0$, and ${V}_{r}^{{z}_{r}{z}_{r}}>0$. The ROI is diminishing if these second derivates are negative: ${V}_{h}^{{z}_{h}{z}_{h}}<0$, ${V}_{r}^{{z}_{h}{z}_{h}}<0$, ${V}_{h}^{{z}_{r}{z}_{r}}<0$, and ${V}_{r}^{{z}_{r}{z}_{r}}<0$.
Thus if either Equation 13 or Equation 14 is positive, then division of labour is favoured to evolve. This gives the second condition for division of labour.
The accelerating returns from individual specialisation condition for division of labour
Fixing our attention on just Equation 15 (or equivalently on 16), this condition states that the sum of the secondorder viability effect from increased cooperation (first term on lefthand side) and the secondorder fecundity effect of increased cooperation (second term on lefthand side) must be larger than the marginal fecundity cost of increased investment in viability (third term on lefthand side). Note that the third term on lefthand side is always negative because increased investment in viability decreases the value of increased investment in fecundity. Therefore, division of labour by a singletrait mutation can only happen if there is an accelerating ROI in at least one of fecundity, $F\left(z\right)$; helper viability, ${V}_{h}({z}_{h},{z}_{r})$; or reproductive viability, ${V}_{r}({z}_{h},{z}_{r})$.
Division of labour by reciprocal specialisation
Request a detailed protocolThe last remaining scenario for division of labour is that the resident strategy of uniform cooperation is unstable to mutations in both traits, which can occur depending on the value of the last term of the Taylor expansion (Equation 2; $W}^{{z}_{h}{z}_{r}}\mathrm{\Delta}{z}_{h}\mathrm{\Delta}{z}_{r$). This kind of instability can arise if there is a joint mutation that affects the level of cooperation of both helpers and reproductives at the same time ($\mathrm{\Delta}{z}_{h}\ne \text{}0$ and $\mathrm{\Delta}{z}_{r}\ne 0$). However, it could also occur if a slightly deleterious mutation in one trait invades by drift and destabilises the other trait so much that the population evolves away from the critical point. In either case, it will be found that same condition must be satisfied in order for division of labour to evolve. In the rest of this section, we give the general analysis of whether $({z}_{h},{z}_{r})=({z}^{*},{z}^{*})$ is unstable to twotrait mutations, and then consider a simplifying special case to clarify the biological interpretation of this analysis.
Suppose that the resident strategy of uniform cooperation ($({z}_{h},{z}_{r})=({z}^{*},{z}^{*})$) is stable against singletrait mutations (i.e. Conditions 15 and 16 not satisfied). Then the resident strategy is unstable to twotrait mutations if and only if the determinant of the Hessian is negative:
This condition is satisfied if the strength of directional selection pushing the population back to the critical point along either of the traitvalue directions is less than the strength of directional selection on a trait when moved off of the critical point along the other trait direction. To evaluate the Hessian condition, we first compute the secondorder cross derivatives:
Here, we have used that ${V}_{i}^{{z}_{h}{z}_{r}}={V}_{i}^{{z}_{r}{z}_{h}}$ . Substituting this, and Equations 13 and 14, into the Hessian condition gives:
The reciprocal specialisation condition for division of labour
Assume that neither Condition 15 nor 16 is satisfied, that is, the individual ROI is nonaccelerating. Then the lefthand side of the inequality is strictly positive, which means that Condition 19 is nontrivial. We will see in examples that Condition 19 can be satisfied, which means that division of labour can evolve by reciprocal specialisation even when the individual ROI is diminishing.
To clarify further the biological meaning of Condition 19, consider a simple family of models in which viability and fecundity are linear functions. In this case, all second derivates are zero, and so we get the simplified condition:
Simplified reciprocal specialisation condition for division of labour
In the case that the viability functions are (Equation 5) this condition further simplifies to:
This inequality is satisfied if the viability of reproductives increases faster with increased cooperation from helpers, than it does from increased cooperation from reproductives. This makes clear that reciprocal specialisation can evolve if reproductives stand to gain more from help from helpers than they would gain by helping themselves.
Appendix 1
A Group structure models
We now consider models that have specified spatial (or ‘network’) structures. We apply the division of labour scenarios outlined in the Materials and methods to study division of labour in these graph models. We find simple, general conditions on the graph for when division of labour can invade.
A.1 Division of labour by betweenindividual differences: stars and branching groups
Consider a star of $n$ cells, with 1 focal cell, and $n1$ peripheral cells (Appendix 1—figure 1). Assume that some function, $H(z)$, determines the amount of a shared resource produced by a cell with level of cooperation $z$. And let $1\lambda $ be the proportion of benefits produced by an individual that it keeps for itself, so that $\lambda $ is the proportion of benefits that is shared equally by a cell’s neighbours. If the peripheral cells in the star have phenotype z_{h} and the focal cell has phenotype z_{r}, then
The viability of a peripheral cell is ${V}_{h}({z}_{h},{z}_{r})=(1\lambda )H({z}_{h})+\lambda H({z}_{r})/(n1)$, and the viability of the focal cell is ${V}_{r}({z}_{h},{z}_{r})=(1\lambda )H({z}_{r})+(n1)\lambda H({z}_{h})$. Equation 12 then gives the following condition for division of labour by betweenindividual differences:
which holds for every positive $\lambda >0$.
In a similar vein, consider a group with a branching structure, where alternating cells have either $a$ or $b$ neighbours, with $a<b$. We term the cells with $a$ neighbours as edge cells and cells with $b$ neighbours as node cells. (In order to avoid edge effects, we assume that the group eventually wraps around on itself, creating the shape of a regular polyhedron: in the case of $a=2$, $b=3$, this would be the dodecahedron, with ${n}_{h}=30$ edge cells and ${n}_{r}=20$ node cells.) This gives the following expected fitness:
where n_{h} is the number of edge cells and n_{r} is the number of node cells. Note that for any graph in which n_{h} cells have degree $a$, and n_{r} cells have degree $b$, we have that ${n}_{h}/{n}_{r}=b/a$. The viability of an edge cell is ${V}_{h}({z}_{h},{z}_{r})=(1\lambda )H({z}_{h})+(a/b)\lambda H({z}_{r})$ and the viability of a node cell is ${V}_{r}({z}_{h},{z}_{r})=(1\lambda )H({z}_{r})+(b/a)\lambda H({z}_{h})$. Equation 12 then gives the following condition for division of labour by betweenindividual differences:
So, since $b>a$, we find that division of labour by betweenindividual differences occurs for all values of $\lambda >0$ (any social trait).
In order to calculate the benefit of division of labour in Figure 4 of the main text, we evaluate $\text{max}(0,W({z}^{*}+\mathrm{\Delta}z,{z}^{*})W({z}^{*},{z}^{*}),W({z}^{*},{z}^{*}\mathrm{\Delta}z)W({z}^{*},{z}^{*}))$, where we approximate fitness using a firstorder Taylor expansion and setting $\mathrm{\Delta}z=0.01$.
A.2 General graph analysis: betweenindividual differences
Consider a graph with vertices $i=1,2,\mathrm{\dots},n$, each vertex with degree d_{i}. Let ${A}_{ij}$ be the matrix of adjacencies: ${A}_{ij}=0$ if there is no edge $ij$, and for every edge $ij$ we have ${A}_{ij}=1/{d}_{j}$. Moreover, ${A}_{ii}=1$, for each $i$. Notice that the column sums of the matrix $A$ are zero:
The fitness function of the group associated with this graph is
where z_{i} is the level of cooperation of the individual at vertex $i$, and ${\delta}_{ij}=1$ if $i=j$ and 0 otherwise. The uniform cooperation strategy, ${z}_{i}=z\forall i$, has fitness
which is maximised for some $z={z}^{*}$ between 0 and 1, by our assumption that F is decreasing in $z$ and H is increasing in $z$.
We now restrict to the ‘marginal’ case that $F(z)$ and $H(z)$ are both linear functions. In this case we can, without loss of generality, take $F(z)=1z$ and $H(z)=z$. The uniform cooperation strategy is given by ${z}^{*}=1/2$. The first derivatives of $W$ at the uniform strategy are given by:
Notice that, for fixed i,
If all the vertices of the graph have the same degree, ${d}_{1}={d}_{2}=\mathrm{\dots}=d$, then ${\sum}_{j}{A}_{ij}=0$, and so the first derivatives all vanish. If one or more of the vertices, $j$, connected to $i$, does not have the degree d_{i}, then the sum, Equation 30, does not vanish. In this case, the uniform strategy is not stable, and division of labour can invade.
For example, consider the branching tree with seven cells, shown in Appendix 1—figure 2. The $A$ matrix of this tree is:
Using Equation 29, it follows that
So vertices 2 and 3 are ‘predisposed’ towards helping less, whereas the more peripheral vertices (i = 1, 4, 5, 6, 7) are predisposed towards helping more. The vertices with the fewest neighbours (i = 1, 4, 5, 6, 7) are also the most predisposed towards helping more.
A.3 A ring of cells
Let us assume that the group is composed of $n$ cells that form a onedimensional filament wrapped into a ring, where $n$ is an even number (Appendix 1—figure 3). Label the cells 1 through $n$, and assume that 1 is neighbours with $n$. We assume that each cell in the filament shares social benefits with only its direct neighbours. Suppose that ‘odd’ cells are putative helper cells, and ‘even’ cells are reproductives. This gives the following fitness for the group:
The viability of helpers is ${V}_{h}({z}_{h},{z}_{r})=(1\lambda )H({z}_{h})+\lambda H({z}_{r})$ and the viability of reproductives is ${V}_{r}({z}_{h},{z}_{r})=(1\lambda )H({z}_{r})+\lambda H({z}_{h})$. We now consider each of the possible routes by which division of labour can evolve.
Using Equations 6 and 7, we find that ${F}^{\prime}/F={H}^{\prime}/H$ at the resident strategy of uniform cooperation for both helpers and reproductives and therefore that the condition for division of labour by betweenindividual differences (Equation 11) can never hold.
Assume that $F$ is a linear function. If the returns from cooperation are nonincreasing (${H}^{\prime \prime}\le 0$), then we use Condition 17 (Figure 5a), to find that division of labour is favoured if and only if
In the case that ${H}^{\prime \prime}=0$, division of labour is favoured only if $\lambda >1/2$. If ${H}^{\prime \prime}$ is less than zero, then the ROIs are decelerating, and $\lambda $ must be larger than $1/2$ in order for division of labour to be favoured.
In order to calculate the fitness benefit of division of labour in Figure 4 of the main text, we evaluate ${max}_{\pi \le \theta \le \pi}W({z}^{\ast}+\mathrm{\Delta}z\mathrm{cos}(\theta ),{z}^{\ast}+\mathrm{\Delta}z\mathrm{sin}(\theta ))W({z}^{\ast},{z}^{\ast})$, where we approximate fitness using a secondorder Taylor expansion and setting $\mathrm{\Delta}z=0.01$. In this case, we can show analytically that the direction of mutation giving the maximal increase in fitness satisfies $\theta =\{\pi /4,3\pi /4\}$.
A.4 A wellmixed group of cells
We now consider a group that is ‘well mixed’ such that cells share the benefits of cooperation with all other group members. That is, we let $1\lambda $ be the amount of benefits produced by an individual that it keeps for itself and $\lambda $ is the amount of benefits that is shared equally by all other cells in the group. This produces the following expected fitness function:
where $n={n}_{h}+{n}_{r}$ and we have that the viability of a helper is ${V}_{h}({z}_{h},{z}_{r})=(1\lambda )H({z}_{h})+\lambda \frac{({n}_{h}1)H({z}_{h})+{n}_{r}H({z}_{r})}{n1}$ and the viability of a reproductive is ${V}_{r}({z}_{h},{z}_{r})=(1\lambda )H({z}_{r})+\lambda \frac{({n}_{r}1)H({z}_{r})+{n}_{h}H({z}_{h})}{n1}$. Plugging these into Equation 11, we find that division of labour cannot evolve by betweenindividual differences as both sides are equal to ${H}^{\prime}/H$, and thus helpers and reproductives face the same viabilityfecundity tradeoff.
If the returns from cooperation and fecundity are linear (${F}^{\prime \prime}={H}^{\prime \prime}=0$), then we can evaluate condition 20 directly to find that division of labour by reciprocal specialisation can evolve if (Figure 5c)
Otherwise, if the returns from cooperation are diminishing (say, ${H}^{\mathrm{\prime}\mathrm{\prime}}<0$, with $F$ linear), then we can use the general condition 17 to find,
Notice that, for diminishing returns (${H}^{\mathrm{\prime}\mathrm{\prime}}<0$), we have that $x>0$, and in this case it follows that
In other words, if ${H}^{\mathrm{\prime}\mathrm{\prime}}<0$, then $\lambda $ must be even higher for division of labour to be favoured, than when ${H}^{\prime \prime}=0$.
In order to calculate the fitness benefit of division of labour in Figure 4 of the main text, we evaluate ${\text{max}}_{\pi \le \theta \le \pi}W({z}^{*}+\mathrm{\Delta}z\mathrm{cos}(\theta ),{z}^{*}+\mathrm{\Delta}z\mathrm{sin}(\theta ))W({z}^{*},{z}^{*})$, where we approximate fitness using a secondorder Taylor expansion and setting $\mathrm{\Delta}z=0.01$.
A.5 General graph analysis: reciprocal specialisation
We now return to the case of a general graph model, and consider when division of labour is possible via reciprocal specialisation. Once again, fix some graph $G$. Assume that the first derivatives, ${\partial}_{i}W$, vanish at the strategy of uniform cooperation (${z}_{i}={z}^{*}$). Recall (Equation 29) that this is equivalent to assuming that all the vertices of the graph $G$ have the same degree, $d$. The matrix of second derivatives of $W$ is
where $M$ is a symmetrised version of the $A$ matrix:
All vertices have the same degree, say $d$, so that the matrix $M$ is given by:
where $L$ is the Laplacian matrix:
The uniform strategy ($z={z}^{*}$) is unstable if and only if the matrix Hess has one or more positive eigenvalues. When $\lambda =0$, $\text{Hess}=2I$, where $I$ is the identity. So all eigenvalues of Hess for $\lambda =0$ are negative, and the uniform strategy is stable. This agrees with the biological intuition: when $\lambda =0$, the group is not interacting socially, and so division of labour cannot evolve. We would like to find ${\lambda}^{*}$ (if it exists), such that division of labour is possible for $\lambda >{\lambda}^{\ast}$. This corresponds to the smallest value of $\lambda $ for which $det(\text{Hess})$ vanishes. We can explicitly evaluate the determinant as:
where the characteristic polynomial $\mathcal{P}(x)$ is a degree $n$ polynomial defined by:
The roots of $\mathcal{P}(x)$ are all nonnegative (this is because $L$ is a symmetric matrix). Moreover, $x=0$ is a root of $\mathcal{P}(x)$. Let μ be the largest root of $\mathcal{P}(x)$, that is, μ is the largest eigenvalue of $L$. Then the smallest value of $\lambda $ for which $det(\text{Hess})(\lambda )$ vanishes is ${\lambda}^{*}=d/\mu $. It follows that division of labour can evolve if and only if
Recall that $d$ is the number of neighbours that an individual cell has in the graph. The eigenvalue μ can be thought of (see below) as a measure of how ‘bipartitionable’ the graph is. The inequality can then be interpreted loosely as follows: if the graph is made ‘more bipartite’, $\lambda $ can decrease. If the number of neighbours, $d$, is increased, $\lambda $ must increase.
To further understand μ, a useful property of $L$ is that, for any vector of numbers $x=({x}_{1},\mathrm{\dots},{x}_{N})$:
The largest eigenvalue of $L$ is
On the other hand,
It follows that
Moreover, if the graph is bipartite, choose any bipartite colouring of its vertices. Then, assigning $x=+1$ to vertices of one colour, and $x=1$ to vertices of the other colour, we find that
This means that μ achieves its maximum value, $\mu =2d$, for bipartite graphs. It can be shown that $\mu =2d$ if and only if the graph is bipartite. μ can be regarded, therefore, as a measure of ‘how bipartite’ a given graph is. Some representative examples for small ($N\le 8$) graphs are given in Appendix 1—figure 4, Appendix 1—figure 5, Appendix 1—figure 6, Appendix 1—figure 7. The rest of this section applies the above analysis to several basic families of examples.
Filaments of cells
A ring of $n$ cells has degree $d=2$, and the largest eigenvalue of the Laplacian matrix depends on whether $N$ is even or odd:
For $n$ even, division of labour is possible if $\lambda >1/2$. For large odd $n$, μ is approximately 4, and so division of labour is again possible for $\lambda >1/2$. However, for small odd $n$, μ is less than 4. This means that, for small $n$, division of labour is ‘easier’ for even numbers of cells than for odd numbers of cells. For some examples, see Appendix 1—figure 5.
Complete graphs
The complete graph with $n$ cells has degree is $d=n1$ and the largest eigenvalue is (In terms of Equation 47, a possible vector that achieves the maximum is ${x}_{1}=+1$, ${x}_{2}=1$, and ${x}_{3},{x}_{4},\mathrm{\dots}=0$. For this vector, ${x}^{T}Lx=4+2(N2)=2N$ and ${x}^{T}x=2$.)
So reciprocal specialisation is possible only if $\lambda >(n1)/n$.
Tessellations
We can carry out similar computations for (infinitely extended) tessellation patterns. There are many bipartite tessellations, in both two and three dimensions. All of these have the maximum possible μ, that is, $\mu =2d$. Examples of this include the hexagonal lattice, the square lattice, and the 3D cube lattice: Appendix 1—figure 4(a), (b), and (e). Reciprocal specialisation is also possible in many tessellations which are not bipartite. An example with $d=6$ is shown in Appendix 1—figure 4(c), which has $\mu =8$. Finally, there are tessellations which are so far from being bipartite that reciprocal specialisation is not possible. An example with $d=8$ is shown in Appendix 1—figure 4(d), which has $\mu =8$. The condition $\lambda \mu >d$ cannot be satisfied with $\mu =d$, since $\lambda $ is at most 1.
B Random role allocation
We repeat the above analyses where we now assume that individuals adopt their phenotypes with a particular probability, independent of other individuals in the group. We find that, with random role allocation, division of labour is not possible (in the absence of accelerating individual ROIs).
B.1 A ring of cells
Let p be the probability of adopting the role of a helper and investing $x$ in cooperation and $1p$ as the probability of adopting the role of a reproductive and investing $y$ in cooperation. In the case of a ring of cells, the number of helpers among a given cell’s two neighbours is a binomial random variable, $k\in \{0,1,2\}$. This gives the expected fitness:
The first term is the expected fitness if the focal cell adopts the role of a helper and the second term is the expected fitness if the focal cell adopts the role of a reproductive. We can simplify the above fitness equation to:
If we multiply the fitness equation by ${n}^{2}$ and approximate ${n}_{h}=pn$ and ${n}_{r}=(1p)n$, then we arrive at the final fitness equation:
where we now have that the expected viability of a helper is ${V}_{h}({z}_{h},{z}_{r})=(1\lambda )H({z}_{h})n+\lambda ({n}_{h}H({z}_{h})+{n}_{r}H({z}_{r}))$ and the expected viability of a reproductive is ${V}_{r}({z}_{h},{z}_{r})=(1\lambda )H({z}_{r})n+\lambda ({n}_{h}H({z}_{h})+{n}_{r}H({z}_{r}))$.
The function, Equation 53, can be analysed using the results of Section B. We find (as expected) that division of labour by betweenindividual differences is not possible for this function:
Moreover, we can evaluate Condition 20 directly to find that division of labour by reciprocal specialisation can only evolve if (for linear functions $F$ and $H$)
It is impossible to satisfy this inequality for $\lambda \le 1$. Thus, assuming linear or decelerating ROIs, we find that division of labour cannot evolve by reciprocal specialisation.
In order to calculate the fitness benefit of division of labour in Figure 4 of the main text, we evaluate ${\text{max}}_{\pi \le \theta \le \pi ,0\le p\le 1}W({z}^{*}+\mathrm{\Delta}z\mathrm{cos}(\theta ),{z}^{*}+\mathrm{\Delta}z\mathrm{sin}(\theta ))W({z}^{*},{z}^{*})$, where we approximate fitness using a secondorder Taylor expansion and setting $\mathrm{\Delta}z=0.01$. In this case, the above expression must now be maximised over two variables. This is done numerically to produce the main text figure.
B.2 A branching group of cells
We approximate the fitness of the group as the expected fitness of a randomly chosen cell in the branching group. With probability ${n}_{h}/n$ we choose an ‘edge cell’ with $a$ neighbours and with probability ${n}_{r}/n$, we choose a ‘node’ cell with $b$ neighbours (here ${n}_{h}+{n}_{r}=n$, $a<b$ and ${n}_{h}/{n}_{r}=b/a$). If the cell is an edge cell, the number of neighbouring cells that are helpers is a binomial random variable $k\in \{0,\mathrm{\dots},a\}$. If the cell is a node cell, the number of neighbouring cells that are helpers is a binomial random variable $k\in \{0,\mathrm{\dots},b\}$. Again, we let p be the probability that the cell adopts the role of a helper. Similar to the process above, we can derive a fitness equation of the form:
where we have approximated ${n}_{h}=pn$ and ${n}_{r}=(1p)n$. We thus have that the viability of a helper is ${V}_{h}({z}_{h},{z}_{r})=(1\lambda )H({z}_{h})n$ $+(\frac{{n}_{h}}{n}\frac{a\lambda}{b}+\frac{{n}_{r}}{n}\frac{b\lambda}{a})$ $({n}_{h}H({z}_{h})+{n}_{r}H({z}_{r}))$ and the viability of a reproductive is ${V}_{r}({z}_{h},{z}_{r})=(1\lambda )H({z}_{r})n$ $+(\frac{{n}_{h}}{n}\frac{a\lambda}{b}+\frac{{n}_{r}}{n}\frac{b\lambda}{a})$ $({n}_{h}H({z}_{h})+{n}_{r}H({z}_{r}))$. For the particular case of $a=2$ and $b=3$ (a branching filament), we can use Equation 11 to show that division of labour by betweenindividual differences can never evolve.
Moreover, if the returns from cooperation and fecundity are linear (${F}^{\prime \prime}={H}^{\prime \prime}=0$), then we can evaluate Condition 20 directly to find that division of labour by reciprocal specialisation can only evolve if $\lambda >1$, which is outside the physically permissible range of the parameter. Thus, division of labour can only evolve in this scenario if there are individual efficiency benefits from specialisation.
In order to calculate the fitness benefit of division of labour in Figure 4 of the main text, we evaluate ${\text{max}}_{\pi \le \theta \le \pi ,0\le p\le 1}W({z}^{*}+\mathrm{\Delta}z\mathrm{cos}(\theta ),{z}^{*}+\mathrm{\Delta}z\mathrm{sin}(\theta ))W({z}^{*},{z}^{*})$, where we approximate fitness using a secondorder Taylor expansion and setting $\mathrm{\Delta}z=0.01$. In this case, the above expression must now be maximised over two variables. This is done numerically to produce the main text figure.
B.3 General graph analysis: random allocation
In general, suppose that the assignment of phenotype is done randomly, according to a probability density function $\rho (z)$, such that $\int dz\rho (z)=1$. Write ${z}_{0}=\u27e8z\u27e9$ for the average, and $\sigma}^{2}=\u27e8{z}^{2}\u27e9{\u27e8z\u27e9}^{2$ for the variance of the distribution. The expected fitness for an arbitrary graph model (Equation 27) is
Putting $F(z)=1z$ and $H(z)=z$, this becomes
where we have used that
This calculation shows that the expected fitness of the group, $\u27e8W\u27e9$, is independent of the group’s graph structure. Moreover, notice that $\u27e8W\u27e9$ is largest (for fixed z_{0}) if the distribution is chosen so that ${\sigma}^{2}=0$, that is, the state of uniform cooperation. Division of labour would correspond to ${\sigma}^{2}>0$. So division of labour cannot be favoured (absent accelerating returns) if phenotypes are assigned randomly.
To see this in greater detail, consider the distribution
Under this distribution, cells become a helper (phenotype z_{h}) with probability $1p$, and a reproductive (phenotype z_{r}), with probability $p$. Then ${z}_{0}=(1p){z}_{h}+p{z}_{r}$, and
The total expected fitness is then
If $\lambda \le 1/2$, the function has a critical point at ${z}_{h}={z}_{r}=z=1/(2(1\lambda ))$, where its first derivatives vanish. We compute that the determinant of the Hessian of $\u27e8W\u27e9$ at the critical point is
which is positive. So
is a maximum of $\u27e8W\u27e9$, for all choices of p, and $\lambda <1/2$. If $\lambda \ge 1/2$, then $\u27e8W\u27e9$ is maximised for ${z}_{h}={z}_{r}=z=1$. In either case, division of labour is not favoured.
C Additional questions
C.1 Does division of labour by reciprocal specialisation produce an accelerating fitness benefit?
We consider a population of uniform cooperators with ESS level of cooperation ${z}^{*}$. We consider a mutation where both helpers and reproductives specialise in their respective functions such that $\mathrm{\Delta}{z}_{h}=\beta \mathrm{\Delta}{z}_{r}=\mathrm{\Delta}z>0$, where $\beta >0$ is the degree to which reproductives specialise more than helpers ($\beta >1$) or less than helpers ($\beta <1$). This gives the fitness of a mutant with reciprocal specialisation as:
where we have assumed that there are no betweenindividual differences and so ${W}^{{z}_{h}}={W}^{{z}_{r}}=0$ at the ESS strategy of uniform cooperation, ${z}^{*}$. What is the shape of the return from increased specialisation in this case? We solve for the second derivative of mutant fitness with respect to the degree of specialisation $\mathrm{\Delta}z$:
If this expression is positive, then the group fitness returns from reciprocal specialisation are accelerating. In contrast, if this expression is negative, then the group fitness returns from reciprocal specialisation are diminishing.
Let us now assume that the mutant has a higher fitness than the resident strategy of uniform cooperation, such that division of labour by reciprocal specialisation evolves ($W({z}^{\ast}+\mathrm{\Delta}z,{z}^{\ast}\beta \mathrm{\Delta}z)>W({z}^{\ast},{z}^{\ast})$). This means that the second term on the righthand side of Equation 66 is strictly positive, which occurs when $\frac{1}{2}{W}^{{z}_{h}{z}_{h}}+\frac{1}{2}{\beta}^{2}{W}^{{z}_{r}{z}_{r}}\beta {W}^{{z}_{h}{z}_{r}}>0$. Importantly, this is the exact same condition as for an accelerating fitness return from reciprocal specialisation. Consequently, division of labour by reciprocal specialisation evolves if and only if the group fitness return from reciprocal specialisation is accelerating.
C.2 Does division of labour without accelerating returns from individual specialisation require that individuals have different viabilities?
In many cases, it makes more sense to consider the viability of the group as a whole rather than the viabilities of the individual cells in the group. For instance, in Volvocine algae, sterile flagella beaters provide a benefit to the group by keeping it afloat at the appropriate height in the water column for light absorption. In this instance, individuals literally sink or swim as a whole group and so all individuals have the same viability. Here, we go through the general framework where we now assume that ${V}_{h}={V}_{r}=V$.
Considering Equation 11, we find that division of labour by betweenindividual differences occurs when:
This specifies that the average contribution by helpers to group viability by a small increase in cooperation is greater than the average contribution by reproductives to group viability by a small increase in cooperation. Thus, putative helpers are predisposed to cooperation.
Let us now assume that there are no betweenindividual differences (${V}^{{z}_{h}}/{n}_{h}={V}^{{z}_{r}}/{n}_{r}$). Considering Equations 15 and 16, we find that division of labour by an accelerating return from individual specialisation can arise if:
We thus find that division of labour can evolve if ${V}^{{z}_{h}{z}_{h}}>0,{V}^{{z}_{r}{z}_{r}}>0$ or ${F}^{\mathrm{\prime}\mathrm{\prime}}>0$. We now assume that there is a nonaccelerating return from individual specialisation ${V}^{{z}_{h}{z}_{h}}=0,{V}^{{z}_{r}{z}_{r}}=0$ and ${F}^{\prime \prime}=0$. Considering Equation 20, we find that the condition for division of labour by reciprocal specialisation is: $V}^{{z}_{h}}/{n}_{h}>{V}^{{z}_{r}}/{n}_{r$. However, this contradicts our assumption that there are no betweenindividual differences (${V}^{{z}_{h}}/{n}_{h}={V}^{{z}_{r}}/{n}_{r}$), and so reciprocal specialisation is not a pathway to division of labour. Consequently, if we assume that all individuals in the group have the same viability (${V}_{h}={V}_{r}=V$), then division of labour evolves by either: (1) betweenindividual differences or (2) an accelerating return from individual specialisation.
In our main analysis, where different individuals may have different viabilities, the consequence to the fitness of the group when helpers have a lower viability is not modelled. Indeed, a helper produces the same amount of public good for the group regardless of its expected viability. This approximation is also made in the Yanni et al., 2020, models and in earlier models by Michod, 2006. More realistically, helpers with a lower viability may be expected to produce a smaller total amount of public good over the course of the group life cycle. This in principle could limit the efficiency benefits of betweenindividual differences and of synergistic reciprocal specialisation. Capturing this feedback between cooperation and viability with a formal model would require additional assumptions, which we leave to future work.
C.3 What happens if the fitness costs and benefits of cooperation are additive?
In the main analysis, we have assumed that cooperation provides a viability benefit at a personal cost to fecundity, where individual fitness is the product of viability and fecundity. Here, we assume instead that the benefits and costs of cooperation are additive, such that increased cooperation leads to an increase in group fecundity at a cost to personal fecundity. For ease of comparison, we maintain the notation ${V}_{h}$ and ${V}_{r}$, which now correspond to the fecundity benefit provided to helpers and reproductives by the cooperation of group members.
This gives the following fitness of the group:
We can employ the Taylor expansion of fitness (Equation 2) to determine the three conditions leading to division of labour. First, we find that division of labour by betweenindividual differences can arise if:
This shows that division of labour can evolve if some individuals (helpers) can provide larger fecundity benefits for the group at the same fecundity cost as others (reproductives). Division of labour by an accelerating return from individual specialisation can arise if:
This shows that division of labour by an accelerating return from individual specialisation can arise if there is an accelerating return from decreased cooperation due to increased personal fecundity ${F}^{\mathrm{\prime}\mathrm{\prime}}>0$, or an accelerating return from increased cooperation due to increased viability of the group (${V}_{h}^{{z}_{h}{z}_{h}},{V}_{r}^{{z}_{h}{z}_{h}},{V}_{h}^{{z}_{r}{z}_{r}},{V}_{r}^{{z}_{r}{z}_{r}}>0$). If none of the above hold, we find that reciprocal specialisation can evolve if:
where we have assumed that $\mathrm{\Delta}{z}_{h}=\mathrm{\Delta}{z}_{r}$. This states that there is interference between the cooperative efforts of helpers and reproductives such that increased cooperation by helpers leads to larger relative fecundity benefits for the group when reproductives cooperate less (${V}_{h}^{{z}_{h}{z}_{r}},{V}_{r}^{{z}_{h}{z}_{r}}<0$). There is an analogous effect when fitness is the product of fecundity and viability (first term in small brackets on righthand side of Condition 17).
Data availability
The Matlab (R2020b) source code used to generate Figures 4DF is available at https://osf.io/nw8gz/.

Open Science FrameworkDoes the evolution of division of labour require accelerating returns from individual specialisation?.https://doi.org/10.17605/OSF.IO/NW8GZ
References

BookPrinciples of Social EvolutionOxford: Oxford University Press.https://doi.org/10.1093/acprof:oso/9780199231157.001.0001

Division of labour and the evolution of extreme specializationNature Ecology & Evolution 2:1161–1167.https://doi.org/10.1038/s4155901805649

The evolution of mechanisms to produce phenotypic heterogeneity in microorganismsNature Communications In press.

Division of Labor during Biofilm Matrix ProductionCurrent Biology 28:1903–1913.https://doi.org/10.1016/j.cub.2018.04.046

An Evolutionary Perspective on SelfOrganized Division of Labor in Social InsectsAnnual Review of Ecology, Evolution, and Systematics 42:91–110.https://doi.org/10.1146/annurevecolsys102710145017

Evolution of selforganized division of labor in a response threshold modelBehavioral Ecology and Sociobiology 66:947–957.https://doi.org/10.1007/s0026501213432

Compartmentalized function through cell differentiation in filamentous cyanobacteriaNature Reviews. Microbiology 8:39–50.https://doi.org/10.1038/nrmicro2242

The multicellular nature of filamentous heterocystforming cyanobacteriaFEMS Microbiology Reviews 40:831–854.https://doi.org/10.1093/femsre/fuw029

Division of labour and the evolution of multicellularityProceedings. Biological Sciences 279:1768–1776.https://doi.org/10.1098/rspb.2011.1999

The Evolution of Soma in the VolvocalesThe American Naturalist 143:907–931.https://doi.org/10.1086/285639

The evolution of cooperation in simple molecular replicatorsProceedings of the Royal Society B 284:20171967.https://doi.org/10.1098/rspb.2017.1967

The evolution of mechanisms to divide labour in nonclonal populationsEcology and Evolution 11:14475–14489.https://doi.org/10.1002/ece3.8067

BookEvolution and the Theory of GamesCambridge: Cambridge University Press.https://doi.org/10.1017/CBO9780511806292

Regulation of Cellular Differentiation in Filamentous Cyanobacteria in FreeLiving and PlantAssociated Symbiotic Growth StatesMicrobiology and Molecular Biology Reviews 66:94–121.https://doi.org/10.1128/MMBR.66.1.94121.2002

BookCaste and Ecology in the Social InsectsUnited States: Princeton University Press.

Relatedness in trait group models of social evolutionJournal of Theoretical Biology 206:355–368.https://doi.org/10.1006/jtbi.2000.2132

Cooperators since life beganThe Quarterly Review of Biology 72:184–188.https://doi.org/10.1086/419766

Expanded social fitness and Hamilton’s rule for kin, kith, and kindPNAS 108 Suppl 2:10792–10799.https://doi.org/10.1073/pnas.1100298108

Individual versus groupoptimality in the production of secreted bacterial compoundsEvolution; International Journal of Organic Evolution 73:675–688.https://doi.org/10.1111/evo.13701

A General Allometric and LifeHistory Model for Cellular Differentiation in the Transition to MulticellularityThe American Naturalist 181:369–380.https://doi.org/10.1086/669151

Developmental decisions in Dictyostelium discoideumDevelopmental Biology 284:25–36.https://doi.org/10.1016/j.ydbio.2005.05.011

Transient heterogeneity in extracellular protease production by Bacillus subtilisMolecular Systems Biology 4:184.https://doi.org/10.1038/msb.2008.18

Division of labour in microorganisms: An evolutionary perspectiveNature Reviews. Microbiology 14:716–723.https://doi.org/10.1038/nrmicro.2016.111

A biological market analysis of the plantmycorrhizal symbiosisEvolution; International Journal of Organic Evolution 68:2603–2618.https://doi.org/10.1111/evo.12466
Decision letter

Raymond E GoldsteinReviewing Editor; University of Cambridge, United Kingdom

Detlef WeigelSenior Editor; Max Planck Institute for Developmental Biology, Germany

Peter J YunkerReviewer; Georgia Institute of Technology, United States
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Decision letter after peer review:
Thank you for submitting your article "Does the evolution of division of labour require accelerating returns from individual specialisation?" for consideration by eLife. Your article has been reviewed by 2 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Detlef Weigel as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Peter J Yunker (Reviewer #2).
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
1) I don't think that the title adequately reflect the novelty of the paper. The work by Rueffler et al. 2012 and Yanni et al. 2020 already provide an answer to this question. We know that accelerating returns are not required. The indepth analysis on the role of topological constraints and coordination are the key novelties and this should be reflected in the title.
2) Figure 4 and the conclusion on the efficiency benefits at the group level (lines 355357) are central to the argumentation of the paper. However, the actual group benefits as a function of λ and individual return functions are nowhere shown in the paper. I think it would be very beneficial if the authors add additional panels to Figure 4, perhaps in the form of heat maps, showing how the group benefit is affected by the two main parameters altered.
3) The contradiction with Yanni et al. 2020 boils down to whether λ > (n1)/n is a realistic scenario or not. Can the authors provide biological examples where this condition could be fulfilled?
4) The findings on lines 431434 are very important, but the data is not presented in a visual form. I recommend to add an additional figure displaying the actual results, again together with group productivity effects as mentioned in my comment 3. above.
5) It is shown that division of labor can evolve with concave returns on investment if an optimal pattern is produced. Further, it is shown that if specialists are randomly distributed on a graph, specialization can only evolve if returns on investment are convex. This is a very nice result. However, while I really like this result, the discussion of it overstates what has been learned. In particular, perfect realizations of patterning and completely random distributions of specialists represent the two extreme cases. While the authors are likely correct that mechanisms that perfectly reproduce patterns are unlikely to evolve prior to division of labor, mechanisms that produce patterns with random errors very well are different. Single celled organisms have many mechanisms through which they modify behavior based on interactions with their environment and each other. Further, alternating patterns are common in biology and physics, and can emerge from simple lowlevel rules. So, the possibility of producing a pattern of alternating specialists – with errors – is plausible. The question then becomes, for a given topological structure, how precise must the pattern be reproduced? If perfect (or nearly perfect) patterning is required, then concave specialization via topology seems unlikely. However, if nearly random patterns favor concave specialization, this barrier would be significantly reduced. The 'error tolerance' of a given topological structure seems likely to depend on its fitness advantage over generalists. While a thorough analysis of this question would be fantastic to see, I do not think it is necessary for publication. Instead, I think it would be sufficient to expand the discussion of this topic to reflect these uncertainties.
In a related vein, "between individual differences" share some similarities with these imperfect patterns. It is unlikely that between individual differences would arise via sophisticated intercellular communication/interactions. But simple mechanisms, such as stochastic switches or persistent environmental interactions can lead to phenotypic differences. If these phenotypic differences arise via spatial interactions, then a pattern emerges as well.
6) I wonder if the way sharing is modeled here is too generous. If a cell shares all of its viability, shouldn't it be dead? Either (1) the viability that is being shared is not all of that individual's viability, in which case the fitness written is not the actual fitness of the organism, (2) this sharing is a onetime act after which the altruistic cell is dead, (3) or there is some other reason why this scenario is reasonable which needs to be further explained.
7) Line 408 These results do not seem contradictory to me. Different assumptions produce different models, which, in turn, produce different results. For example, if I construct a model with no sharing, division of labor will not evolve, even for accelerating returns. This does not contradict models that include sharing, but is the result of a different assumption / scenario. It can, of course, be discussed which approach is applicable in different scenarios, and which may be more broadly applicable.
8) Line 453 "However, the biological plausibility of any mechanism based on preexisting cues would need to be explicitly justified and modelled…"
Line 458 "further studiessuch as ancestralstate reconstructionsare needed to show whether coordination preceded division of labour in individual species"
These are true statements. However, similarly true statements can be made about nearly any mechanism that facilitates the evolution of division of labor. For example, accelerating returns. Placing all of these mechanisms within this context, or presenting the evidence that exists for some but not others, would be helpful.
9) Equation 3 was another very nice result. However, I think the discussion here could also be clarified. This model appears likely to hold as long as the standard deviation of d is small compared to the mean d. Consider an organism that has 1000 somatic cells that each only connect to one central germ cell; consider also a linear return on investment. Equation 3 suggests that λ must be greater than 0.5. In this absurd case, specialists should produce a higher fitness than generalists even if λ is less than 0.5. Again, while additional work determining the exact limits in which Equation 3 applies would be nice, I do not think they are necessary. Instead, appropriate caveats could be added to this effect.
Reviewer #1:
The paper consists of two parts.
The first part deals with "the three pathways to division of labour". It builds on a mathematical model showing that division of labour can evolve when (1) there are accelerating returns for individuals from specialisation, (2) there are preexisting differences between individuals such that some individuals are predisposed for one task or the other, and (3) there is reciprocal specialisation leading to synergistic efficiency benefits at the group level. This model recovers the findings by Rueffler et al. 2012. This part is very well written and reads more like a review, bringing specialised and nonspecialised readers onto the same page.
The second part focusses on pathways (2) and (3), where individual returns can be diminishing, and thus the evolution of division of labour is more difficult to explain. This part is motivated by a recent paper by Yanni et al. 2020 in eLife showing that topological constraints (e.g. small network of individuals, limited number of neighbours) are essential to favour division of labour with diminishing returns. The authors challenge this view and provide an indepth analysis on topological constraints. They show that such constraints in combination with preexisting differences between individuals and reciprocal specialisation are indeed conducive for the evolution of division of labour, but not essential. They show that when the benefit of cooperation is larger for neighbours than for the cooperator than division of labour can evolve with diminishing returns even in the absence of topological constraints. This is a key new insight. But perhaps even more important, the authors highlight that pathways (2) and (3) rely on the assumption that individuals have access to information from neighbours to coordinate their actions at the group level. The authors show that division of labour cannot evolve with diminishing returns when such information is absent. And because mechanisms of information collection and coordination are likely to only evolve once division of labour is already in place, the authors argue that topological constraints might play a minor role in driving the initial evolutionary steps towards division of labour.
In brief, this is a very insightful paper and significantly advances our theoretical and conceptual understanding of division of labour. It will spur future theoretical and empirical work in the field, and for the latter, the authors present guidelines of how to test the theory.
Reviewer #2:
This is a very interesting paper on the evolution of division of labor. In particular, the authors explore the impact of topology, the convexity and concavity of fitness returns on investment, and different biological 'modalities' through which division of labor may arise. This is a difficult topic to study as, in most lineages, division of labor evolved long ago, and thus cannot be directly studied in the lab. Further, many theoretical predictions have proven difficult to directly test. This manuscript furthers our understanding of the underlying theory of the evolution of division of labor, and presents a means to test which modality is responsible for the emergence of division of labor in different cases.
However, there are a few caveats worth mentioning. Comparisons to previous works are not always clear. Different models built with different assumptions can produce different predictions; however, that does not mean they disagree, only that they describe different scenarios. Further, the model used here allows an entity to give away all of its 'viability,' making it unclear how it continues to live and function. Finally, the order in which division of labor and 'patterning' evolve is presented as definitive, when it is ultimately a postulate.
https://doi.org/10.7554/eLife.71968.sa1Author response
Essential revisions:
1) I don't think that the title adequately reflect the novelty of the paper. The work by Rueffler et al. 2012 and Yanni et al. 2020 already provide an answer to this question. We know that accelerating returns are not required. The indepth analysis on the role of topological constraints and coordination are the key novelties and this should be reflected in the title.
We have changed the title. There are several novel results in the analysis and so we have chosen a title that we think captures its scope, without focus on any one particular result or use of any jargon.
2) Figure 4 and the conclusion on the efficiency benefits at the group level (lines 355357) are central to the argumentation of the paper. However, the actual group benefits as a function of λ and individual return functions are nowhere shown in the paper. I think it would be very beneficial if the authors add additional panels to Figure 4, perhaps in the form of heat maps, showing how the group benefit is affected by the two main parameters altered.
This is a really nice idea. We have replaced the relevant panels in Figure 4 with heat maps as suggested. This now shows both when division of labour is favoured to invade (nonwhite shading) and the relative fitness benefit of division of labour when it does (relative darkness of shading).
3) The contradiction with Yanni et al. 2020 boils down to whether λ > (n1)/n is a realistic scenario or not. Can the authors provide biological examples where this condition could be fulfilled?
We have expanded discussion of this point, highlighting why there are differences in the predictions of each model, and giving biological examples of cooperation that benefits neighbours more than producers. Lines 420438.
4) The findings on lines 431434 are very important, but the data is not presented in a visual form. I recommend to add an additional figure displaying the actual results, again together with group productivity effects as mentioned in my comment 3. above.
Thank you for the suggestion, we have expanded Figure 4 to include new panels that show when division of labour is favoured to evolve under fully random specialisation. These illustrate that an accelerating return form specialisation is required in all three models.
5) It is shown that division of labor can evolve with concave returns on investment if an optimal pattern is produced. Further, it is shown that if specialists are randomly distributed on a graph, specialization can only evolve if returns on investment are convex. This is a very nice result. However, while I really like this result, the discussion of it overstates what has been learned. In particular, perfect realizations of patterning and completely random distributions of specialists represent the two extreme cases. While the authors are likely correct that mechanisms that perfectly reproduce patterns are unlikely to evolve prior to division of labor, mechanisms that produce patterns with random errors very well are different. Single celled organisms have many mechanisms through which they modify behavior based on interactions with their environment and each other. Further, alternating patterns are common in biology and physics, and can emerge from simple lowlevel rules. So, the possibility of producing a pattern of alternating specialists – with errors – is plausible. The question then becomes, for a given topological structure, how precise must the pattern be reproduced? If perfect (or nearly perfect) patterning is required, then concave specialization via topology seems unlikely. However, if nearly random patterns favor concave specialization, this barrier would be significantly reduced. The 'error tolerance' of a given topological structure seems likely to depend on its fitness advantage over generalists. While a thorough analysis of this question would be fantastic to see, I do not think it is necessary for publication. Instead, I think it would be sufficient to expand the discussion of this topic to reflect these uncertainties.
We have expanded the discussion in the section to highlight that a perfect allocation of labour may not in principle be required, but simply that cell specialisation cannot be fully random (Lines 461466).
In a related vein, "between individual differences" share some similarities with these imperfect patterns. It is unlikely that between individual differences would arise via sophisticated intercellular communication/interactions. But simple mechanisms, such as stochastic switches or persistent environmental interactions can lead to phenotypic differences. If these phenotypic differences arise via spatial interactions, then a pattern emerges as well.
We have expanded the discussion of this point to include the possibility that the mechanism that produces between individual differences may be coopted to coordinate division of labour (Lines 484486).
6) I wonder if the way sharing is modeled here is too generous. If a cell shares all of its viability, shouldn't it be dead? Either (1) the viability that is being shared is not all of that individual's viability, in which case the fitness written is not the actual fitness of the organism, (2) this sharing is a onetime act after which the altruistic cell is dead, (3) or there is some other reason why this scenario is reasonable which needs to be further explained.
This is a good point, which we had discussed in supplementary section C.2. Our assumptions align with scenario (2) above and are consistent with the modelling assumptions of the Yanni et al. (2020) paper and Michod (2006) paper. We have added a description of this assumption in the main text (Lines 135138) and expanded the discussion in the supplementary information (section C.2).
7) Line 408 These results do not seem contradictory to me. Different assumptions produce different models, which, in turn, produce different results. For example, if I construct a model with no sharing, division of labor will not evolve, even for accelerating returns. This does not contradict models that include sharing, but is the result of a different assumption / scenario. It can, of course, be discussed which approach is applicable in different scenarios, and which may be more broadly applicable.
We have expanded discussion of the assumptions made in the different models. Please see our response to revision point (3) above.
8) Line 453 "However, the biological plausibility of any mechanism based on preexisting cues would need to be explicitly justified and modelled…"
Line 458 "further studiessuch as ancestralstate reconstructionsare needed to show whether coordination preceded division of labour in individual species"
These are true statements. However, similarly true statements can be made about nearly any mechanism that facilitates the evolution of division of labor. For example, accelerating returns. Placing all of these mechanisms within this context, or presenting the evidence that exists for some but not others, would be helpful.
We have added a paragraph to highlight why we think that a larger burden of evidence is needed for betweencell coordination then for accelerating returns from individual specialisation (Lines 495500). In short, coordination requires an adaptive argument for why such betweencell behaviour will have evolved prior to or concurrently with the emergence of division of labour. In contrast, the shape of the return from more investment is simply an inherent aspect of the biotic or abiotic environment and so does not require an adaptive justification.
9) Equation 3 was another very nice result. However, I think the discussion here could also be clarified. This model appears likely to hold as long as the standard deviation of d is small compared to the mean d. Consider an organism that has 1000 somatic cells that each only connect to one central germ cell; consider also a linear return on investment. Equation 3 suggests that λ must be greater than 0.5. In this absurd case, specialists should produce a higher fitness than generalists even if λ is less than 0.5. Again, while additional work determining the exact limits in which Equation 3 applies would be nice, I do not think they are necessary. Instead, appropriate caveats could be added to this effect.
We have added text to clarify this point (Lines 332336). Equation 3 is arrived at by assuming that all cells have the same number of neighbours (standard deviation of d is zero). If this is not the case, then different cells have different numbers of neighbours and therefore division of labour by between individual differences will evolve (as established in the section just prior to this one). If that is the case, the second order analysis that goes into deriving Equation 3 is inconsequential to whether division of labour can evolve.
https://doi.org/10.7554/eLife.71968.sa2Article and author information
Author details
Funding
H2020 European Research Council (834164)
 Stuart Andrew West
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Thomas Scott and Zheren Zhang for their helpful comments and suggestions; St John’s College, Oxford (GAC), and the ERC (Horizon 2020 Advanced Grant 834164; SAW) for funding.
Senior Editor
 Detlef Weigel, Max Planck Institute for Developmental Biology, Germany
Reviewing Editor
 Raymond E Goldstein, University of Cambridge, United Kingdom
Reviewer
 Peter J Yunker, Georgia Institute of Technology, United States
Publication history
 Received: July 6, 2021
 Accepted: October 26, 2021
 Accepted Manuscript published: October 29, 2021 (version 1)
 Version of Record published: January 25, 2022 (version 2)
Copyright
© 2021, Cooper et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 855
 Page views

 218
 Downloads

 1
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading

 Evolutionary Biology
 Genetics and Genomics
The exceptionally rich fossil record available for the equid family has provided textbook examples of macroevolutionary changes. Horses, asses, and zebras represent three extant subgenera of Equus lineage, while the Sussemionus subgenus is another remarkable Equus lineage ranging from North America to Ethiopia in the Pleistocene. We sequenced 26 archaeological specimens from Northern China in the Holocene that could be assigned morphologically and genetically to Equus ovodovi, a species representative of Sussemionus. We present the first highquality complete genome of the Sussemionus lineage, which was sequenced to 13.4× depth of coverage. Radiocarbon dating demonstrates that this lineage survived until ~3500 years ago, despite continued demographic collapse during the Last Glacial Maximum and the great human expansion in East Asia. We also confirmed the Equus phylogenetic tree and found that Sussemionus diverged from the ancestor of noncaballine equids ~2.3–2.7 million years ago and possibly remained affected by secondary gene flow postdivergence. We found that the small genetic diversity, rather than enhanced inbreeding, limited the species’ chances of survival. Our work adds to the growing literature illustrating how ancient DNA can inform on extinction dynamics and the longterm resilience of species surviving in cryptic population pockets.

 Evolutionary Biology
 Genetics and Genomics
Estimating the complex relationship between fitness and genotype or phenotype (i.e. the adaptive landscape) is one of the central goals of evolutionary biology. However, adaptive walks connecting genotypes to organismal fitness, speciation, and novel ecological niches are still poorly understood and processes for surmounting fitness valleys remain controversial. One outstanding system for addressing these connections is a recent adaptive radiation of ecologically and morphologically novel pupfishes (a generalist, molluscivore, and scaleeater) endemic to San Salvador Island, Bahamas. We leveraged wholegenome sequencing of 139 hybrids from two independent field fitness experiments to identify the genomic basis of fitness, estimate genotypic fitness networks, and measure the accessibility of adaptive walks on the fitness landscape. We identified 132 single nucleotide polymorphisms (SNPs) that were significantly associated with fitness in field enclosures. Six out of the 13 regions most strongly associated with fitness contained differentially expressed genes and fixed SNPs between trophic specialists; one gene (mettl21e) was also misexpressed in labreared hybrids, suggesting a potential intrinsic genetic incompatibility. We then constructed genotypic fitness networks from adaptive alleles and show that scaleeating specialists are the most isolated of the three species on these networks. Intriguingly, introgressed and de novo variants reduced fitness landscape ruggedness as compared to standing variation, increasing the accessibility of genotypic fitness paths from generalist to specialists. Our results suggest that adaptive introgression and de novo mutations alter the shape of the fitness landscape, providing key connections in adaptive walks circumventing fitness valleys and triggering the evolution of novelty during adaptive radiation.