Skip to main content
Advertisement
  • Loading metrics

Pleiotropic constraints promote the evolution of cooperation in cellular groups

Abstract

The evolution of cooperation in cellular groups is threatened by lineages of cheaters that proliferate at the expense of the group. These cell lineages occur within microbial communities, and multicellular organisms in the form of tumours and cancer. In contrast to an earlier study, here we show how the evolution of pleiotropic genetic architectures—which link the expression of cooperative and private traits—can protect against cheater lineages and allow cooperation to evolve. We develop an age-structured model of cellular groups and show that cooperation breaks down more slowly within groups that tie expression to a private trait than in groups that do not. We then show that this results in group selection for pleiotropy, which strongly promotes cooperation by limiting the emergence of cheater lineages. These results predict that pleiotropy will rapidly evolve, so long as groups persist long enough for cheater lineages to threaten cooperation. Our results hold when pleiotropic links can be undermined by mutations, when pleiotropy is itself costly, and in mixed-genotype groups such as those that occur in microbes. Finally, we consider features of multicellular organisms—a germ line and delayed reproductive maturity—and show that pleiotropy is again predicted to be important for maintaining cooperation. The study of cancer in multicellular organisms provides the best evidence for pleiotropic constraints, where abberant cell proliferation is linked to apoptosis, senescence, and terminal differentiation. Alongside development from a single cell, we propose that the evolution of pleiotropic constraints has been critical for cooperation in many cellular groups.

Introduction

There is widespread cooperation in cellular groups where cells invest in costly traits that benefit all cells in the vicinity, such as bacteria that secrete an extracellular enzyme to digests nutrients or the more complex coordinated phenotypes of multicellular organisms. Cooperative traits can require that cells forego their own reproductive interests in favour of the reproductive interests of the group as a whole [1,2]. This effect, in turn, can lead to the evolution of noncooperative lineages—sometimes known as “cheaters”—that make use of collective benefits without investing in them and threaten cooperative function [3].

The potential for cheater lineages is well documented in microbes. Mutants lacking a range of cooperative traits have been shown to outcompete wild-type cells [410] and occur in the field and clinic [1113]. Such observations beg the question of how cooperation persists over evolutionary time. A key explanation is that many cellular groups, both in microbes and multicellular organsims, are recently derived from a single cell (clonal) [1,2,14,15]. In the terminology of sociobiology, this leads to high relatedness and kin selection, which is a major driver of cooperation across many systems [16,17]. The argument is that, when cell groups are clonal, interactions between cooperative and cheater genotypes are prevented, which allows cooperative genotypes to prosper as cheater genotypes lose the shared benefits of cooperation. While kin selection is undoubtedly important [1,2,14,15], this explanation neglects a key feature of the biology of cellular groups: Even in a group founded from a single cell, cooperation can still break down due to the emergence of mutant noncooperators from within [1820] (Fig 1). Indeed, with nonzero mutation rates, the question is when, not if, these lineages will emerge.

thumbnail
Fig 1. The problem of cooperation and the evolution of pleiotropy.

(A) Mutation–selection dynamics can undermine clonal multicellular groups. Mutation of cooperative cells (blue) can generate noncooperative cells (dark grey) that do not pay the costs of cooperation but reap its rewards. These cheater lineages can spread within multicellular groups because they divide more rapidly than wild-type cells. Invasion of spontaneous mutant lineages can lead to a breakdown in group function (distorted shape). (B) Pleiotropy promotes cooperation in our model because it leads to between-group selection on the rate at which cooperation breaks down within groups as they age. Beginning with groups founded by 1 pleiotropic and 2 nonpleiotropic lineages, we see that groups founded by the nonpleiotropy lineages loose function and are eventually replaced by the pleiotropic lineage. Cells are coloured in 3 parts according to whether they display cooperation (blue), a privately beneficial trait (red), and pleiotropy (yellow), whereby the cooperative and private traits are linked. Mutations can make some traits inactive (grey). When cooperation is lost in a pleiotropic cell lineage (top row), the cells also lose their private trait, which stops them from proliferating. Group function is thereby preserved. When cooperation is lost in a nonpleiotropic lineage (second row), cheater lineages emerge that spread and group function is lost. As a result, pleiotropic groups thrive relative to nonpleiotropic ones and seed more groups, giving rise to between-group selection for pleiotropy.

https://doi.org/10.1371/journal.pbio.3001626.g001

Given the potential for cheater lineages, it has been suggested that genetic architecture can evolve to help stabilise cooperation. When the expression of a cooperative trait is linked to a private trait that helps a cell to survive or divide, mutations that ablate cooperation can also ablate the private trait and, thereby, stop the evolution of cheater lineages. This pleiotropic linkage of cooperative traits and private (personally benefitial) traits has been identified experimentally to be a mechanism that helps to promote cooperation in multiple microbial species [2125]. In addition, an agent-based model of microbes found that links between metabolic and secretion genes helped to stabilise cooperation [26], and a theoretical study suggested that pleiotropy can promote niche construction, which is related to cooperation [27]. These studies raise the possibility that the evolution of pleiotropy might be a general mechanism to promote cooperation in cellular groups [28].

However, a recent theoretical paper argued broadly against the idea that pleiotropy is an explanation for the evolution of cooperation [29]. In particular, the authors argued that pleiotropy only evolves under conditions when kin selection is already operating to stabilise cooperation (Fig 2 in [29]), with, at best, a very minor impact on the evolved level of cooperation (i.e., seen in S14 and S16 Figs but not Fig 2 in [29]). As such, they concluded “Pleiotropy does not help stabilise cooperation over evolutionary time—cooperation is only favoured in the region where Hamilton’s rule is satisfied because of indirect fitness benefits.

As we disuss in detail in the Supporting information (S1 Text), a limitation of this model is that it did not explicitly capture group-level birth and death events or allow groups to develop for long enough to see the importance of pleiotropy for stabilising the evolution of cooperation. The time allowed for groups to develop is important because the problems with cheater lineages only becomes apparent as groups age (Fig 1). Pleiotropy only becomes subject to significant between-group selection, therefore, in longer-lived groups when cheater mutants have time to threaten the group. We show in the Supporting information that increasing the length of time that groups live for increases the levels of cooperation that evolve via pleiotropy in the model of [29]. However, problematic assumptions such as unbounded explonential growth prevented us from exploring this effect further (S1 Text). We, therefore, decided to develop a novel age-structured multilevel selection model for the evolution of cooperation in cellular groups, including both microbes and multicellular organisms. Our model predicts that pleiotropy is a powerful way to promote the evolution of cellular cooperation.

Results

We are interested in understanding how multicellular groups founded by cells with pleiotropic constraints function as compared to groups founded by otherwise similar cells that lack these constraints. We follow the effects of pleiotropic links between cooperative traits (that benefit the whole group) and private traits (that benefit the individual cell that carries them) on multilevel selection dynamics using an age-structured modelling approach (see Methods). A group in our model is intended to capture a group of microbes or a proto-multicellular organism, which lacks the division between germ and soma. Groups start from a single cell and display logistic growth up to a carrying capacity K, which defines the size of the group at maturity. A second parameter, λ, determines the expected life span of a group. This is important because it impacts on the amount of cell turnover that is expected after a group reaches reproductive maturity. Such cell turnover can be major contributor to the number of cell divisions within a multicellular group. For example, high rates of cell turnover occur in bacteria, which commonly live attached to surfaces in structures known as biofilms where dispersing cells are replaced by dividing cells below them [1]. It is also common In multicellular organisms: Members of the genus Hydra (Fig 1) can live for several years, while their epithelial cells are estimated to turnover every few days [30].

The fact that groups start from a single cell in our model ensures high relatedness and strong kin selection, which is consistent with microbes that grow in clonal patches [1] and the biology of multicellular organisms [20]. However, we later reduce this within-group relatedness to study its effects on pleiotropy and cooperation. To study the effects of pleiotropy on cooperation, our modelling has to capture the stochastic effects of mutations. For this reason, the heart of the model is a stochastic simulation that captures populations of cells as they grow, and potentially mutate to other genotypes, within a group. However, as we discuss later, we also need to capture the evolutionary effects of cooperation at the group level, which is done with partial differential equations (PDEs) that allow us to capture a large (infinite) number of competing groups. We hope that this novel approach—stochastic simulations embedded in PDEs—will prove useful to understand a wide range of traits under multilevel selection (Methods).

Pleiotropy slows the breakdown of cooperation within cellular groups

Regulatory networks, and the maps from genotype to phenotype, are often complex [31]. Evolutionary models of cooperation typically overlook this complexity and instead study optimal trait values, an approach known as the phenotypic gambit [32]. Here, we treat genetic architecture as a trait, like any other, that can itself evolve in response to natural selection [26,29]. We do this with a simplified model of pleiotropy. The definition of pleiotropy can vary between disciplines and authors [29,3335], and here we mean the commonly used definition: Pleiotropy is when a single locus affects 2 or more traits [24]. Specifically, our model captures how mutation at a given locus affects 1 cooperative and 1 private trait (Fig 2). While there are a vast range of possible regulatory networks that might influence any 2 traits of interest, the impacts of pleiotropy can be captured by a single value ϕ, which is the probability that a mutation in a network with an active, cooperative, and private trait will give rise to a pleiotropic effect (Fig 2A).

thumbnail
Fig 2. Modelling pleiotropy.

(A) One way to measure the strength of pleiotropy in real-world regulatory networks is to compute the ratio of those mutations in the network that simultaneously impact 2 terminal phenotypic traits to the number of genes. We call this measure ϕ. In principle, this value can be calculated for any network regulating 2 traits. We show illustrative networks that would generate varying levels of pleiotropy across the range of values of ϕ. (B) Genotype–phenotype map for our mathematical model of pleiotropy. To keep things simple, we model 3 traits: a terminal cooperative trait, (blue circles), a terminal private trait, (red circles), and a pleiotropic regulation trait, (yellow circles). All traits in our model can be either active (coloured as above) or inactive (grey circles). This means there are 8 possible genotypes, each labelled gi, where i∈{1,2,…,8}. Although pleiotropy can be active or inactive for a given genotype, we use a parameter ϕ∈[0,1] to tune the strength of its effect. This allows us to use our simple genotype–phenotype map to model biological scenarios when pleiotropy is expected to be weak as well as strong. Furthermore, we assume that pleiotropy only has a functional effect when both other traits are active. Thus, although active pleiotropy is present in genotypes g2, g4, g6, and g8, we assume it only affects genotype g8, where it influences the probability, ϕ, that a g8 cell experiences a pleiotropic mutation given that a loss-of-function mutation has deactivated one of its other traits.

https://doi.org/10.1371/journal.pbio.3001626.g002

There are, therefore, 3 traits in our model: a cooperative trait, a private trait, and a pleiotropy trait, which gives rise to 8 possible genotypes (Fig 2B). The goal of our model is to explore which of these 8 genotypes is favoured by natural selection in the long term due to competition among cells within and between groups. For a given cell, each of the traits is in either an active or inactive state. At the heart of the model is a tension between selection for cooperation within and between groups, i.e., the cooperative trait decreases a cell’s relative division rate within a group but brings benefits to the group function as a whole. This trait might, for example, represent a secreted enzyme that helps nutrient acquisition in a microbial group, or the suppression of cell proliferation in a simple multicellular organism to ensure proper functioning [18]. By contrast, the private trait simply increases a cell’s survival rate within its group. This trait might represent an enzyme involved in central metabolism, for example. With pleiotropy, mutations that ablate the cooperative trait increase the probability that the private trait is also lost and vice versa [29], where ϕ determines this probability (Fig 2A). A mutation matrix specifies the transition probabilities between all 8 genotypes in the model as a function of the strength of pleiotropy (see Methods).

Under a null model in which pleiotropy does nothing, mutations affect each trait independently, where mutations that cause loss-of-function in a trait occur with rate μ, and gain-of-function trait mutations occur with rate νμ, where ν<1. This value reflects the fact that it is typically easier to break trait functionality than to restore or create it, and we typically take ν = 0.01 to capture the strength of this bias. We use μ = 0.0001 per generation for the base mutation rate in most analyses, which describes the probability that a trait is mutated—and function is lost—per cell division. This value is expected to vary widely between systems and traits and is intended only as an illustration. We later perform parameter sweeps of both μ and ν across several orders of magnitude.

We begin by following the evolutionary dynamics within a group. In our first model, each group is founded by a single cell, which gives 8 possible group types corresponding to the 8 cell genotypes (Fig 2B). While groups all start their life with clonal expansion of their founder, mutation–selection processes mean that their genotypic composition may change through time as they age. We can describe this process for each of the 8 group types. The dynamics for groups founded by genotypes 1 to 7 are shown in the Supporting information (S1S7 Figs), and we focus here on genotype 8 groups (Fig 2B), hereafter referred to as “pleiotropic cooperators,” because they capture the effects of pleiotropy on cooperation (Fig 3). Groups with these genotypes initially grow towards their carrying capacity by clonal expansion but, depending on the strength of pleiotropy, have the potential to be invaded by cheater lineages that lack the cooperative trait but express the private trait (genotype 4; see Fig 2A). Importantly, we see that the extent and rate of invasion of the cheater lineage is diminished as the strength of pleiotropy, ϕ, is increased. Cheater lineages make up 25% of the group by approximately day 25 in groups without pleiotropy, by day 40 in groups with intermediate pleiotropy, and never (not before 50 days) in groups with strong pleiotropy (Fig 3B).

thumbnail
Fig 3. The evolution of cooperation within groups depends on the strength of pleiotropy ϕ.

Within-group mutation selection dynamics are shown for a group founded by a cell with genotype g8, which actively expresses a cooperative trait, , a private trait, , and a pleiotropy trait, . Growth of the group as its age, y, increases, is logistic, with a carrying capacity K = 200 at which point cells continue to divide, die, and turnover (see Methods). Dynamics are shown from left to right for 3 different strengths of pleiotropy, ϕ, where pleiotropy is absent/has no effect in the left-hand side column (ϕ = 0). For comparison, the vertical dashed line in (A-C) shows the point at which noncooperative lineages shown in orange make up 25% of the group. Pleiotropy leads to this point being delayed (ϕ = 0.5) or prevented (ϕ = 1). (A) Changes in genotype abundances, nc(y). (B) Changes in genotype relative frequencies, xc(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, . (D) Distribution of fitness effects: Shown is the effect on within-group fitness of the different types of loss-of-function mutations that occur, which are coloured by their phenotypic effect, where blue is loss of cooperative trait, red is loss of the private trait, yellow is loss of pleiotropy, and brown is loss of both cooperation and private trait (due to pleiotropy). When pleiotropy is weak or absent, loss-of-function mutations to cooperative traits increase the fitness of cells within the group, and loss-of-function mutations to private traits decrease the fitness of cells within the group. When the strength of pleiotropy is increased, mutations to either trait tend to have pleiotropic effects, which cancel one another out, meaning mutant cell lineages no longer gain an advantage within the group. Formally, the fitness effect is , where and are the within-group birth rates of the mutant descendant and ancestral wild type, respectively, and and are the death rates of the mutant descendant and ancestral wild type, respectively. Parameters: sc = sg = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.g003

The resistance to invasion by cheater lineages occurs because pleiotropy reduces the frequency with which mutations give rise to a cheater phenotype. As a result, pleiotropy is able to increase the level of cooperation in groups (Fig 3C). The distribution of fitness effects (DFE) of loss-of-function mutations helps to show why pleiotropy is an effective mechanism for limiting cheater cell lineages within a given group (Fig 3D), something also clear from the dos Santos study [29]. When the strength of pleiotropy is relatively weak, mutations to the cooperative trait frequently give rise to mutant descendants that have a competitive advantage over the cooperative cells within the group. By contrast, when the strength of pleiotropy is relatively strong, mutations tend to have have neutral or deleterious effects on cells because a loss of cooperation also comes with a loss of the private trait. In the model, we assume that the effects of expressing the cooperative and private trait on within-group fitness are equal and opposite in magnitude, which is what leads to neutrality when both are lost. Some examples suggest that the loss of a private trait may have a stronger negative effect, such as cell death via apoptosis [3638]. Such examples may lead to a negative change in within-group fitness when both traits are lost. We do not consider this case explicitly here, but it is only expected to strengthen the ability of pleiotropy to remove potential cheater lineages and thereby improve group function.

Pleiotropy evolves to suppress cheater lineages and promote cooperation

Our within-group model supports the established, and intuitive, argument that pleiotropic links between a cooperative and private trait will help to maintain cooperation [2125]. However, this model simply assumes that this pleiotropy exists, rather than explaining how it evolved. One origin of pleiotropy is a result of natural selection on traits unrelated to cooperation [25,39]. Pleiotropy is extremely common in all genotype to phenotype maps, whether or not cooperative traits are involved. As a result, cooperation may become pleiotropically linked to private traits simply through the way that regulatory networks normally evolve. A more intriguing alternative is that pleiotropy evolves to promote cooperation [25,31]. However, as discussed above, a recent theoretical treatment of this idea argued that, despite the abilitiy of pleiotropy to decrease the cheater load within groups, it will not generally evolve to promote cooperation over evolutionary time [29]. Specifically, this earlier study found that “(1) pleiotropy does not stabilise cooperation, unless the cooperative and private traits are linked via a genetic architecture that cannot evolve (mutational constraint); (2) if the genetic architecture is constrained in this way, then pleiotropy favours any type of trait and not especially cooperation; (3) if the genetic architecture can evolve, then pleiotropy does not favour cooperation; and (4) there are several alternative explanations for why traits may be linked, and causality can even be predicted in the opposite direction, with cooperation favouring pleiotropy.”

The authors did find that pleiotropy will often help reduce the prevalence of cheaters within certain groups, in line with experimental evidence in microbes [2123]. However, they also found that between-group selection was not strong enough to notably increase cooperation across the whole population over evolutionary time (Figs S14 and S16 in [29] show a small increase in evolved cooperation with pleiotropy). We discuss this study in detail in the Supporting information (S1 Text), where we conclude that a different type of model will benefit the study of the evolution of cooperation and pleiotropy.

We developed our model, therefore, to study whether pleiotropy will evolve as a mechanism to promote cooperation within multicellular groups over evolutionary time. To do this, we extend the model to capture how groups perform, and compete, across a wider population. In each group, within-group evolution occurs as just discussed (Fig 3), which, in turn, affects a group phenotypic trait, which we call the group’s “function.” For example, this might represent the ability of a bacterial strain to produce a protective biofilm, or the ability of a multicellular organism to coordinate its development. The key is that the group function is assumed to breakdown when within-group selection dynamics take hold.

We capture group function at age y by , where is the average amount of cooperation in the group, and is the average expression of the private trait (see Methods for more details). Thus, we assume that cells expressing both cooperative and private traits contribute fully to group function, whereas cells with either trait missing do not. Specifically, cells lacking only cooperation behave as cheater lineages, as just discussed, while cells lacking the private trait function poorly. As a result, groups founded by these cells (genotypes 1 to 6; Fig 2B) have a greater probability of extinction than higher functioning groups (genotypes 7 and 8; Fig 2B). But within-group evolutionary dynamics also means that group function can degrade in groups founded by cooperative cells (Fig 3).

We study the evolution of these processes with a system of PDEs, which captures an infinitely large age-structured population of groups. We start with a population dominated by individuals lacking all traits (genotype 1) and allow them to evolve by numerically integrating within-group and between-group dynamics over time until a stable age distribution is reached. We then ask if pleiotropy evolves as a function of 3 key parameters: the strength of pleiotropy, ϕ, expected group life span, λ, and group size at maturity, K.

A common approach to model multilevel selection in evolutionary biology is the haystack model [40], which is the approach used in [29]. In the haystack model, group selection is modelled implicitly as an emergent process resulting from fitness differences between individuals within groups; there is no explicit consideration of group-level survival and reproduction events. Instead, all group phenotypes are described in terms of the set of individual phenotypes. In contrast, the PDE approach we use models group-level phenotypes and considers their direct impact on group-level survival and reproduction events. This approach allows us to more explicitly capture microbial groups and multicellular organisms, which have group-level traits that are important for group fitness but are threatened by within group competition.

Our analysis reveals a wide range of parameters where pleiotropy evolves and, in doing so, promotes cooperation over evolutionary time (Fig 4). In particular, when the effects of pleiotropy are absent in the model (ϕ = 0), the levels of cooperation that evolve are often markedly decreased (Fig 4). These results includes cases where cooperation evolves to a lower level without pleiotropy but not to 0, which is also seen to a modest extent in the supplementary figures of dos Santos and colleagues [29]. We also see many cases where cooperation is negligible without pleiotropy, but present with pleiotropy, predicting that pleiotropy can enable the evolution of cooperation in regions where kin selection alone does not maintain it (Fig 4, S8 Fig). Moreover, we do not see the evolution of pleiotropy in a control model of 2 private traits (Methods; S9 Fig). These 2 contrasts—removing pleiotropy and removing cooperation—demonstrate that cooperation can depend on the evolution of pleiotropy and vice versa.

thumbnail
Fig 4. The evolution of pleiotropy promotes cooperation.

Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 group sizes (increasing from top to bottom). Dotted line on each heatmap indicates separation between point where pleiotropy has no effect (control case ϕ = 0) to regions where it has a nonzero influence on the dynamics. Cooperation evolves for a wide range of parameters but is destabilised for longer group life spans, λ, and larger group sizes, K, due to the emergence of cancerous cell lineages. This effect is strongest without pleiotropy (ϕ = 0). When pleiotropy has an effect (ϕ>0), natural selection favours its evolution resulting in higher levels of cooperation, i.e., the blue intensity increases from left to right in the cooperation plot. As expected, pleiotropy is most favoured when it is more effective, i.e., the strength of pleiotropy, ϕ, is higher. Note that the level of expression of the pleiotropic trait can be nonzero even when it has no effect due to stochastic effects. In these situations, however, the expression of the pleiotropic trait does not influence cooperative evolution. We show below that pleiotropy will also evolve with a cost, which greatly reduces this stochasticity (S11 Fig). Parameters: sc = sg = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.g004

Our model predicts that pleiotropy is more important for cooperation as groups become larger and longer lived. This is because, if groups are small or short lived, there is less opportunity for cheater lineages to arise and interfere with group functioning (Fig 3). As expected, the strength of pleiotropy—how protective pleiotropy is against invasion of mutant cheater genotypes—is an important factor in determining when pleiotropy evolves (variation in the x-direction in Fig 4). Nevertheless, we find that pleiotropy evolves and increases cooperation even for low levels of protection (ϕ<0.5, blue plots, Fig 4). In some cases, pleiotropy evolves in the model when it has no phenotypic effect via genetic drift (Fig 4, ϕ = 0). As expected, this effect is strongest when group sizes are small and short lived because this is when natual selection is also generally weakest. Below, we introduce a cost to pleiotropy that largely removes this effect and shows it does not contribute to our findings.

Another assumption of potential importance is how detrimental the invasion of cheater lineages is to overall group function. If the invasion of even a few mutant cheaters can damage group function, then mechanisms that resist that invasion are likely to be more strongly favoured by between-group selection. To be conservative, we do not assume a high detrimental impact where only a few cheater cells are fatal for group functioning but instead consider a linear function where some multicellular function persists in the face of very large numbers of noncooperator cells. Specifically, recall that group function at age y is given by , a function that declines linearly with the invasion of noncooperative mutants (all else being equal). However, even without a high detrimental impact associated with a small numbers of cheaters, we find that pleiotropy is important for cooperation in groups that generate only 5,000 to 10,000 cells across their life span (i.e., 5 × 103 to 104 cell divisions in total). This prediction is borne out in Fig 4 where the total number of cell divisions in the model corresponds roughly to group size multiplied by the number of generations in the life span (i.e., y-axis value multiplied by group size in Fig 4; see Methods). As discussed, these are only rough estimates as they depend on mutation rate and other assumptions that will vary between systems. However, given that many cellular groups undergo many more cell divisions than this estimate, these results predict that the evolution of pleiotropy has the potential for widespread impacts.

Our model also allows us to follow the order in which the traits evolve in the population (Fig 5). Without pleiotropy (ϕ = 0), cooperative genotypes can rapidly evolve, but so too do noncooperative genotypes (red lines in Fig 5A). As a result, these cheater lineages end up making up a substantial proportion of the population such that cooperation and group functioning are limited (Fig 5B, left panel). By contrast, when pleiotropy can influence the distribution of mutational effects, we see that it evolves extremely rapidly after the origin of cooperation itself. Indeed, for the higher strengths of pleiotropy, it evolves alongside the initial origin of cooperation in the population (Fig 5B, middle panel).

thumbnail
Fig 5. Long-term evolutionary dynamics of cooperation and pleiotropy.

Dynamics are shown for the global population of cells over time, t, which encompasses many generations of cell groups. These dynamics encompass both within-group and between-group selection dynamics. We show 3 strengths of pleiotropy, ϕ, to capture 3 scenarios where stronger pleiotropy is associated with more rapid and complete evolution of both pleiotropy and cooperation. (A) Changes in global genotype relative frequencies, xc(t). (B) Changes in the global average levels of cooperation, private trait expression, and pleiotropy, . (C) Average change in traits over a group lifetime, measured as the difference between the average trait values among groups aged y to those expected from their founding cell at birth. (D) Order in which genotypes invade in the global population, where a blue arrow indicates a gain of cooperation (via between-group selection), a grey arrow indicates a loss of cooperation (via within-group selection) and a black arrow indicates no change in cooperation. Parameters: sc = sg = 0.95; K = 200; μ = 0.0001; ν = 0.01; λ = 20. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.g005

Rather than being a secondary adaptation restricted to derived multicellular groups, therefore, our model predicts that pleiotropy can be important at the origin of cooperation within multicellular groups. One can also assess the effects of pleiotropy in terms of trait-based mutational load: the average trait values of a group at age y relative to their trait value at birth (Fig 5C). The mutational load is greatest for the cooperative trait, with groups tending to express reduced cooperation aged y than they do at their origin. However, this mutational load of cooperation is reduced with the evolution of stronger pleiotropy (Fig 5C). In this way, pleiotropy does not just act to reduce the emergence of cheaters in a given group (Fig 3), it acts across the whole population and can greatly increase the evolved level of cooperation (Figs 4 and 5).

Pleiotropy evolves to stabilise cooperation across a wide range of conditions

We have so far assumed that groups are formed from a single cell. While this is realistic for the majority of multicellular organisms, other cellular groups, particularly microbial groups, commonly contain multiple genotypes that meet and mix. If large numbers of different genotypes meet and mix—and relatedness is close to 0—the evolution of pleiotropy and indeed cooperation does not occur in our model. Under these conditions, there is no between-group genetic variation and the outcome of competition is determined solely by within group dynamics (S10 Fig). Here, so long as genotypes that lack the cooperative phenotype can arise at some point, they will take over and pleiotropy serves no function. However, relatedness can often be relatively high in microbial groups due to spatial structure, where a patchwork of groups form, each dominated by a single genotype [1]. We can study the effects of an intermediate level of relatedness in our model by assuming groups are founded by 2 cells (chosen uniformly at random from their parent group), such that there are now up to 32 different group genotypes in the population. This case has an important difference to the single-cell bottleneck case where cheater cells always start a new group alone with little chance of survival. With 2 cells, cheater cells now have the chance of founding groups alongside cooperators that they can exploit, thus greatly improving their prospects. Despite the added complexity, we see again that the evolution of pleiotropy is often favoured and able to promote the evolution of cooperation as it evolves (S11 Fig). While the importance of pleiotropy in our model rests upon some relatedness between cells, therefore, it does not rest upon a single cell origin.

Our conclusions are also robust to changing other assumptions and parameters. One key consideration is that there may be a cost to pleiotropy if, for example, the regulation of 1 trait is compromised by its linkage to another [41]. However, we find that pleiotropy still evolves if it carries such costs to a group’s functioning (S12 Fig), which is further testament to its ultimate importance for improving group function. Another important parameter is the benefit of cooperation (strength of group selection). Reducing the benefit of cooperation in our model reduces the scope for the evolution of cooperation but, importantly, where cooperation can evolve there are broad parameter ranges where pleiotropy evolves to increase cooperation (S13 Fig). Notably, the evolution of pleiotropy is even seen when natural selection for cooperation is very weak, as may have occurred at the inception of multicellular life.

Varying the relative probability of gain-of-function mutations has little impact on outcomes (S14 Fig). However, as expected, the baseline mutation rate is important. Increasing this mutation rate causes the more rapid breakdown of cooperation, which requires stronger pleiotropic effects for cooperation to be maintained. However, so long as strong pleiotropic links are possible, we see that they rapidly evolve and again stabilise cooperation (S15 Fig). For reduced mutation rates, cheater lineages arise less often and so, even in the absence of pleiotropy, cooperation can be maintained more easily. All else being equal, therefore, pleiotropy will now only evolve in larger or longer-lived cellular groups. For example, halving the mutation rate (μ = 0.00005) roughly doubles the number cell divisions where pleiotropy becomes critical for cooperation (compare Fig 4 with S16 Fig). However, even if we lower the mutation rate an order of magnitude (μ = 0.00001), we still observe the widespread evolution of pleiotropy in groups of only 10,000 to 20,000 cell divisions (S16 and S17 Figs show this effect, with and without a cost to pleiotropy, respectively). In summary, we observe that the evolution of pleiotropy promotes cooperation for relatively small multicellular groups across a wide range of parameters.

Pleiotropy is predicted to be important in a simple model of multicellular organisms

Our modelling assumptions are most suited to multicellular groups of microbes and probably some of the ancestral organisms that gave rise to multicellularity in the algae, plants, animals, and fungi. With the evolution of derived multicellularity came many complexities, which our models do not capture. Central among these is the importance of the germ and soma separation for modern multicellular life. We, therefore, next ask whether this biology influences the evolution of pleiotropy and multicellular cooperation.

Here, we assume that multicellular groups can influence which genotype they transmit during reproduction via a germ line that undergoes fewer cell divisions, and mutations, than the soma. Specifically, with probability γ, an ancestral group seeds a new group with a cell of the same genotype to its founding cell, and with probability 1−γ, a cell is chosen uniformly at random from the group to seed a new descendant (as before). As discussed above, we also assume a cost to pleiotropy to be conservative. The model predicts, as expected, that a germ line is generally benefitial for the evolution of cooperation, as it reduces the chance that noncooperative genotypes will start new groups. Nevertheless, with a germ line, we still see the widespread evolution of pleiotropy (Fig 6).

thumbnail
Fig 6. The evolution of pleiotropy in multicellular groups with a germ line.

(A) Cheater lineages that threaten cooperation evolve in multicellular species with a germ line, as well as in simpler multicellular groups. In multicellular organisms, these lineages manifest themselves as tumours and cancers that emerge by mutation. From left to right, bacteria, simple and complex animals. Bacteria: Spontaneous lasR mutant emerge during evolution of the pathogen Pseudomonas aeruginosa. The mutants do not contribute to the production of public goods molecules that their wild-type counterparts do, which enables them to gain a short-term competitive advantage (Image credit: Sheyda Azimi and Steve Diggle). Simple animal: Spontaneous tumour formation in the basal metazoan Hydra, a tumour-bearing Hydra oligactis polyp (right) is shown next to a healthy animal (left). Credit: Alexander Klimovich, Kiel University. Complex animal: Cancer is widespread in long-lived animals, especially zoo animals where other selection pressures are minimised. Nanuq, a 29-year-old polar bear, died from liver cancer in The Columbus Zoo, Ohio, in 2017 (image shows a different polar bear; credit: Alan D. Wilson). (B) To study the impact of pleiotropy on multicellular groups with a germ line, we varied γ, a measure of the strength of the transmission of the germline, modelled as the likelihood that a group propagates the genotype of its founding cell at the age at which it reproduces versus a cell selected at random. Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 germ line strength parameters (increasing the strength of the germline transmission between ancestor and descendent groups from top to bottom). Increasing the strength of the germline has a positive effect on the evolution of cooperation, but pleiotropy continues to be favoured in regions of the parameter space in which cooperation is vulnerable to breakdown. The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = sg = 0.95; K = 200; μ = 0.0001; ν = 0.01; ζ = 0.02. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.g006

A second important characteristic of derived multicellular organisms, like humans, is the need to reach a certain age before reproduction is possible. Up until now, reproduction among our multicellular groups has been age independent. We therefore introduce another parameter, α, which determines the fraction of the expected life span that groups must reach before they can reproduce. For example, if the expected life span is λ = 50 and α = 0.5, then groups can only reproduce after age y = 25. On its own, adding this requirement reduces the levels of cooperation in the population, because it means that organisms reproduce when they have the highest levels of cheater mutations (S18 Fig). However, the negative effects of late life reproduction on cooperation are reduced if we assume that these organisms also have a germ line (Fig 7). Moreover, this effect rests upon the ability of an organism to evolve effective pleiotropic constraints (ϕ > 0). We find, therefore, that organisms with a germ line and delayed reproduction will evolve pleiotropy to protect against cheater lineages and ensure high-level functioning when they reproduce (Fig 7C).

thumbnail
Fig 7. The evolution of pleiotropy in multicellular groups with a germ line and delayed reproductive maturity.

We varied α, the fraction of the expected life span λ groups must have aged to before they can reproduce. We also assumed a germ line, γ = 1, and that the evolution of pleiotropy carries a cost. Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 reproductive maturity parameters (increasing the age at which maturity is reached from top to bottom). Increasing the age of reproductive maturity favours the evolution of pleiotropy. The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = sg = 0.95; K = 200; μ = 0.0001; ν = 0.01; ζ = 0.02; γ = 1. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.g007

Discussion

Our models predict that pleiotropy can greatly increase the levels of cooperation within multicellular groups (Fig 4 and S11S18 Figs). Natural selection favours pleiotropy because it reduces the probability that cheater lineages arise within a cellular group or multicellular organism, which would otherwise damage functioning (Fig 3). We expect this process to be important whenever groups become large or old enough for the emergence and spread of cheater lineages to cause significant harm to multicellular function (Fig 4). These conditions are expected in relatively small and simple multicellular groups. We find pleiotropy can evolve to suppress cheater lineages in groups that undergo as few as 104 cell divisions.

This 104 cell divisions estimate depends on parameters, most notably mutation rate. Mutation rates vary widely, but to give an idea of natural rates, a simple animal like Hydra vulgaris would be expected to have between roughly 1 and 100 mutations per genome per cell division, based upon its approximately 1 Gbp genome and recent estimates of somatic mutation rates taken from mammals (10−9–10−7 mutations per basepair per cell division) [42,43]. If we divide these mutation rates by the number of genes in H. vulgaris (20,000), this gives rates between one-half and 50 times our standard value of μ (0.00005 to 0.005). For a trait affected by a single gene, this range may still overestimate the rate of loss-of-function mutations, because not all mutations in a gene will hit a coding region nor lead to a loss of function if they do. However, the traits we are interested in here—for example, regulated cell proliferation—are often affected by large networks of genes, such that mutations in many genes have the potential to influence the trait, for example, drive unregulated cell proliferation. Our standard value, therefore, appears to be reasonable, but in species with lower mutations rates, a larger number of cell divisions are required before pleiotropy is needed to maintain cooperation (S16S17 Figs). However, even if 104 cell divisions is an underestimate by 1 or even several orders of magnitude, many multicellular groups will remain above the threshold at which pleiotropy is predicted to be important. There are an estimated 105 cells in hydra [44], 105 neurons in the Drosophila brain [45], 1010 cells in a bacterial colony [46], and 1013 cells in humans [47], where these large standing cell numbers can still greatly underestimate total cell divisions because of cell turnover [30]. If relatively few cell divisions are needed for the evolution of pleiotropy, this suggests that pleiotropy may have had an important role in cheater suppression at the origins of multicellularity. Further consistent with this hypothesis, the evolutionary dynamics in the model predict that pleiotropy will evolve rapidly, close to simulatenously, with the first evolution of cooperation (Fig 5).

We have identied a scenario where pleiotropy evolves because of its positive effects on cooperation, but pleiotropy can also arise for other reasons [25,35,39,48]. Pleiotropy is widespread for all types of traits—cooperative and noncooperative—and can readily arise, for example, as a way to coregulate multiple traits [39,48]. Such regulatory structures might help to stabilise a cooperative trait or help cooperation to first arise when there is positive frequency-dependent selection [49]. A challenge for the future is to distinguish between examples where the effects of pleiotropy on cooperation arose independently of cooperation, and cases where cooperation drove the evolution of pleiotropy. Our within-group model underlines that, no matter why pleiotropy first evolves, it has the potential to promote cooperation (Fig 3). Moreover, when individuals are related, this can generate between-group selection for pleiotropy and the stabilisation of cooperation in the long term (Figs 4 and 5).

There are other evolutionary processes than those we have modelled here that may lead to cooperation being associated with pleiotropy. Whenever cooperation is linked to the ability of a species to compete and persist in an ecosystem [50], for example, species-level selection may enrich for species that have evolved pleiotropic links that promote cooperation over those that do not [51]. Another potential route to pleiotropy is when 1 individual evolves to enforce cooperation in another. For enforcement to be effective, the enforcer needs to find a way to constrain the recipient in some way from escaping the enforcement. This need can result, therefore, in a pleiotropic link in a recipient where the cessation of cooperation is linked to a personal cost from enforcement [25,28]. This case is notable in that it does not rest on there also being positive relatedness between interacting individuals (cf., S10 Fig).

Enforcement may be important in microbial examples of pleiotropy, including the case of Dictyostelium discoideum [24] where enforcement of 1 cell type by another appears to be important for cooperation [52]. However, our modelling also suggests that pleiotropy may evolve in microbes to limit the spread of cheater lineages. Mutation rates can be relatively low in microbes [53], and the arrival and mixing of multiple genotypes within a microbial group has the potential to both undermine the evolution of cooperation and the potential for pleiotropy with it [29]. However, cooperative traits are known to be widespread in microbes, particularly in the well-studied bacteria [54]. Moreover, groups of bacteria reach vast numbers and the evolution of noncooperator lineages has been observed [21,55,56]. Whenever these lineages cause significant harm to cooperative function, therefore, there is the potential for pleiotropy to evolve to prevent their emergence.

We predict that the evolution of pleiotropy will be particularly important for multicellular organisms. Development from a single cell (Fig 4 versus S11 Fig), combined with relatively high somatic mutation rates (S15S17 Figs) and the potential for a very large numbers of cell divisions (Fig 4 and S16 and S17 Figs), are all predicted to favour pleiotropy as a means to promote cooperation. Consistent with this prediction, many multicellular organisms suffer from cancer, which has long been viewed as analogous to the evolution of cheater lineages within cellular groups [29,31,5761]. Moreover, the study of cancer has identified powerful mechanisms to supress potential cheater lineages, which function by linking oncogenic traits (loss of cooperation) to secondary traits that limit viability (loss of private trait) [31,38,62,63]. For example, loss of function mutations in the retinoblastoma tumour suppressor gene are linked to the activation of programmed cell death (apoptosis) [36,37], while activation of ras genes—key regulators of cell proliferation—are associated with cellular senescence and cell cycle arrest [64]. Other mechanisms push cells with oncogenic mutations to differentiate from a proliferative cell type (stem cell) to one that does not divide [63]. The links between cooperation and cell viability can be cell intrinsic but they also occur via the actions of immune cells and othe cell types [62,65], which can promote apoptosis in a mutant cell [66,67] or even physically force it out of a proliferative tissue [68]. The biology of multicellularity, therefore, appears to be characterised by a large number of pleiotropic connections between abberant cell proliferation that threatens cooperation and the loss of traits that cells need to survive and divide.

The effects of pleiotropy in our models is testament to the importance of considering molecular mechanisms in the study of cooperation, particularly in cellular systems [31]. In microbes, for example, mechanisms such as prudent regulation of cooperative traits, quorum sensing, and green beard genes can all serve to promote cooperation [6971]. High relatedness is also clearly important for cooperation across a range of cooperative systems, and is particularly important for cooperation in cellular groups. Experiments have demonstrated how relatedness promotes cooperation in microbes, and development from a single cell was likely to have been central to the evolution of complex multicellularity [14,15]. However, the large numbers of cell divisions that occur in many cellular groups means that single-cell ancestry is far from sufficient to maintain cooperation. Indeed, our models predict that even relatively small and simple cellular groups will suffer if cheater lineages are left unchecked. This suggests that, in addition to the single-cell bottleneck, the evolution of pleiotropic constraints may have been important for the origins of multicellularity and the subsequent major transition in evolution that led to the complex multicellular life that we see today.

Methods

A general model of multilevel selection using an age-structured model

Capturing the full dynamics of a multilevel selection process is difficult because the potential for selection within and between groups creates an enormous space of possible outcomes. For these reasons, many early models of multilevel selection, such as the haystack model [40], attempted to simplify the problem by neglecting to model group-level events explicitly. While this simplifies the model, the approach misses the fact that group-level events and individual-level events can happen on different timescales, something that is particularly important for capturing the biology of most multicellular groups.

To capture the 2 levels of selection explicitly, we study the dynamics of natural selection in an age-structured population of multicellular groups and in subpopulations of cells within those groups. Changes in the distribution of groups of different ages and types evolves with time t and age y. Changes in the distribution of cells of different types within each group evolves with age y. Using a similar approach to [72], we use a system of PDEs to model changes in the age-structured population of different types of groups over time. We model changes in abundances of different cell types within the different group types as they age using stochastic simulations. These simulations then allow us to characterise and describe how groups of different types differ in their development as they age. Within-group differences in development further provide the basis for differences in reproductive success between groups founded by different types of cells. All the code used to perform our numerical analyses is open source and can be found in GitHub (https://github.com/euler-mab/pleiotropy) and Zenodo (https://zenodo.org/record/6367788#.YjSBVurP2Uk).

We begin by deriving and describing the structure of our model in general terms, before turning to the specifics of how we use it to study the evolution of cooperation and pleiotropy.

Cell and group types

We assume there is finite set of cell types, i,jI, and a finite set of group types, k,lK. We refer to cells of type i as i-type cells and multicellular groups of type k as k-type groups. To distinguish between other group-level and cell-level variables in our model, we use superscripts g and c, respectively. The type of a group is assumed to be defined by the type(s) of its founding cells. In our standard model, we assume that a group is founded by a single cell. In this case, the set of cell genotypes and group genotypes is the same (i.e., I = K). In a later model, we assume that a group is founded by 2 cells. In this case, the set of group types contains all unordered pairs of cell types, and the index notation becomes more cumbersome, although all the same principles for the model hold true. We therefore focus our description here on the simpler model. Let us now describe how the abundances of cells and groups changes over time.

Between-group dynamics

We begin with the group population. In the limit as the population of groups gets large, we assume that relative density of k-type groups aged y can be modelled as a continuous quantity, . In the absence of births and deaths of whole groups, all groups in the population simply age. What this means is that if there were a population of 10 groups at time t aged y = 1, then at time t+10, there will be 0 groups aged y = 1, and 10 groups aged y = 11.

Consequently, the relative density of k-type groups aged y at time t changes over time according to a system of PDEs satisfying a conservation law of the form (1) where there is no change in the overall density of groups in the population. This law can be derived as follows. First, note that within a particular age range [y1, y2], the total abundance of k-type groups aged y at time t is given by (2)

If we assume that there are no births or deaths of groups within the age range, then the abundance of individuals in the age range [y1, y2] can only change because of a process of ageing. Groups of a younger age may enter this age range at the lower age y1, and those within the age range may get older than the upper age bound y2. If we suppose that groups age at a constant rate v, then the rate of change of k-type groups aged y at time t is just . The rate of change of the total abundance of individuals at time t in the age range [y1, y2] is given by (3) where it is equal to the flux of the ageing population over the boundaries of the age range. Let us now integrate both sides of this equation to get an expression for the abundance of k-type groups in [y1, y2] at time t2>t1 in terms of the abundance of k-type groups aged y at time t1 and the total flux at each boundary during this time period. We have (4)

Assuming is differentiable, then we can use the fundamental theorem of calculus to rewrite the integrands on both sides as (5) and (6) respectively, where we have used Leibniz integral rule to take the derivative operators inside the integrals on the right-hand side. Substituting back in to Eq (4) and rearranging, we have (7)

Since this integral is 0 for any arbitrary age range and any time interval, we must conclude that the integrand itself is exactly 0: (8)

We can always normalise the rate of ageing to a constant v = 1, meaning that this equation simplifies to the conservation law given by Eq (1) we started with above.

These equations are similar to a PDE derived by Burt Simon (see [72]), which was also used to study group selection. An important distinction is that our approach allows one to study group-selection using a system of n PDEs, rather than 1 PDE evolving on an n-dimensional surface. The system of PDEs can simplify the numerical problem of solving group-selection dynamics dramatically but is generally only tractable for biological systems in which groups can be assigned to categories in some simple way. Here, for example, we categorise groups by the genotype of their founding cell. If groups were formed from multiple founding cells, then the number of group types would increase dramatically, making the numerical solution of these equations much more difficult.

To incorporate assumptions about the birth and death of k-type groups in our PDE model, along with their production by other l-type groups due to mutation, we can simply extend Eq (1) to include source and sink terms. The production of k-type groups aged y due to the reproduction of all other l-type groups can be represented by a term where is the rate at which l-type groups aged y′ reproduce at time t, and is the conditional transition probability that a k-type group aged y is produced, given that an l-type group aged y′ just reproduced (in practice, we assume the age of a newly produced group is always y = 0). The death of groups can simply be modelled by a term where is simply the rate at which k-type group aged y die at time t. Incorporating these terms into Eq (1) gives (9)

The extent to which there is between-group competition in our model therefore depends on variation between different types of groups in birth rates, , death rates, , and mutation transition probabilities, . These functions all depend on the within-group dynamics within each type of group.

Within-group dynamics

Now let us describe what happens within groups as they age. We assume all groups are founded by a small number of cells (1 cell in most of our analyses) and undergo logistic growth as they age. The abundance of i-type cells within a k-type group aged y is denoted . We assume that external forces have no impact on the change in abundance of cells of different types within a group. Thus, the rate of change in the abundance of i-type cells within a k-type group aged y depends only on age y, and not on external time t.

Changes in the abundance of different cell types within a particular k-type group are stochastic, but we assume that all groups of type k generally develop in the same way as they age, irrespective of the environment they were born into. In practice, we therefore characterise the expected development of a k-type group by calculating the average behaviour of 10,000 replicates of the stochastic dynamics of cells within the group. We use stochastic simulations rather than deterministic solutions to characterise the within-group dynamics because we are interested in the different rates at which mutant cell lineages invade within a finite subpopulation. If we used ordinary differential equations to characterise the same behaviour, cheater mutants would emerge and spread within each group deterministically very early on in the lifetime of groups due to the assumption of infinite population sizes. In contrast, our averaging approach shows that in finite systems, mutant cell lineages invade at vastly different rates in different types of groups. This is crucial for showing why pleiotropy is important in group selection. Groups founded by pleiotropic cooperators can outcompete groups founded by nonpleiotropic cooperators because pleiotropy slows the rate of invasion of mutant cell cheaters, making groups more competitive.

The stochastic dynamics within groups are characterised by a birth–death process with mutation. The birth rate of i-type cells in a k-type group aged y is given by . The death rate of i-type cells in a k-type group aged y is given by . Finally, the conditional transition probability that an i-type cell is produced, given that a j-type cell just reproduced (akin to a mutation transition probability), is denoted . With these rates defined, we can simulate the stochastic evolutionary dynamics of any population. We generate sample paths using Gillespie’s Direct method [73].

Solving the equations

So far we have described how selection and mutation acts on groups and cells within groups. How are these equations linked? The answer is that the group birth rate, , death rate, , and transition probability, all depend on the outcome of the within-group cell dynamics. Thus, to solve the system of PDEs describing the dynamics of selection and mutation between groups, we first need to resolve the dynamics of selection and mutation between cells within each type of group.

We resolve the within-group dynamics of cells using custom-built stochastic simulations written in Matlab (R2019b). We then numerically solve the between-group dynamics given by Eq (9) using a custom-built finite-volume PDE solver written in Matlab (R2019b). We use a “superbee” flux limiter to control the rate of change of group relative density in regions where the group density function is not smooth [74].

Modelling the evolution of pleiotropy and cooperation

Now that we have described a general model of multilevel selection in an age-structured population, we turn to describing how we use that model to study the evolution of cooperation and pleiotropy. Specifically, we focus on describing the birth, death, and mutation rates for cells and multicellular groups in our model. For reference, all of the parameters used in the model and their default values and ranges are shown in Table 1.

Cells

Cells in our model express 3 phenotypic traits of interest, which can either be in an active or inactive state. There are therefore 8 possible cell genotypes in our model (Fig 2B). We consider a cooperative trait , which is beneficial for the reproductive success of multicellular groups as a whole, but costly for the reproductive success of cells within the lifetime of a group, a private trait, , which is beneficial for the reproductive success of a cell expressing it within a group and beneficial for multicellular groups as a whole, and a “pleiotropy” trait, , which influences the types of cell mutations that can occur. Ordering cell traits as row vectors, [], the genotype–phenotype map for cells can be represented by a matrix (10)

Each of the 8 rows corresponds to a genotype, and each of the 3 columns corresponds to a trait, ordered public (cooperative), private, and pleiotropy. Cells with different phenotypes may vary in their ability to survive and reproduce within groups. The within-group birth and death rates of i-type cells in a k-type group aged y are given by (11) and (12) respectively, where sc is the strength of selection on cell traits, and are the average expression of the cooperative and private traits in a k-type group aged y, respectively, and Nk(y) is the size of a k-type group aged y. The component in the death-rate prevents a group from dying because of stochastic extinction of all its cells, because its death rate is when the size of the group is .

Note that the expression of the cooperative trait places a cost on i-type cell division relative to other cells in the group, the expression of the private trait, , gives i-type cells a survival advantage relative to other cells in the group, and has no impact on cell birth rates or death rates at all. We are assuming here that the cooperative trait positively impacts the birth rate of cells and the loss of the private trait positively impacts the death rate of cells. These assumptions are based on the observations that extracellular growth factors and enzymes in cellular groups (cooperative traits) often promote growth [1], while cells lacking a functional metabolism (a private trait) often die [75].

We generally expect within-group selection to disfavour expression of the cooperative trait , favour the expression of the private noncooperative trait , and be neutral with respect to the pleiotropic trait . Similar to dos Santos and colleagues, we assume that selection is of equal strength with respect to the cooperative and private traits because we did not want to introduce biases by privileging 1 trait over another, but rather focus on the importance of pleiotropy [29]. In a variant of our model, we switch the cooperative trait to another private trait, . To do this, we replace our usual birth rate function with (13) where the expression of the private trait on i-type cell division is positive relative to nonexpressing cells in the group.

Cell mutations can occur during cell division. The conditional transition probability that an i-type ancestral cell produces a j-type descendant cell, given that it has reproduced is given by (14) where ϕ∈[0, 1] is a parameter scaling the strength of pleiotropy, and and are elements of 2 different mutation matrices, the first of which represents a model of pleiotropic mutations, and the second of which represents a null model where pleiotropy has no effect. The parameter ϕ thus scales the likelihood that mutation rates are sampled from a mutation matrix in which pleioptropy has an effect versus a mutation matrix in which pleiotropy is absent. When ϕ is relatively small, mutations to traits are almost always independent events, and pleiotropy is therefore relatively weak or absent, but when ϕ is close to unity, mutations are often nonindependent events, where mutation of 1 trait influences the state of another trait.

The specifics of our asssumptions about the effects of pleiotropy are as follows. Under our null model, we have a mutation matrix , where we simply assume that mutations affect each trait independently. We assume that loss-of-function mutations occur with rate μ and that gain-of-function mutations occur with rate νμ, where ν≤1. The mutation matrix for this model is then given by (15)

Under our model of pleiotropy, we have a mutation matrix , where we assume that a loss-of-function mutation to either or drives a corresponding change in the other trait (i.e., that loss-of-function mutations have pleiotropic effects). However, we assume that pleiotropy is only active in the cell genotype in which all 3 traits are active (i.e., pleiotropic cooperators with phenotype ). The mutation matrix for this model is then given by (16) where the only difference between Pc and Qc is in the final row.

Multicellular groups

Groups in our model express their own phenotypic traits of interest that emerge from interactions of cells within a group. We usually assume that groups are founded by a single cell, in which case there are 8 group genotypes in our model. In a variant of our model, however, we assume that groups are founded by 2 cells. In this scenario, there are 8(8−1)/2+8 = 36 group genotypes. We assume that both the cooperative and private traits of cells are essential components of the ability of the groups to function as a cohesive whole in order to maximise survival and reproduction of the group. The functionality of k-type groups aged y is denoted and defined as (17) where there is a multiplicative effect of the average expression of the cooperative trait, , and average expression of the private trait, , among cells within the group at age y. It is important to note that whereas cell-level traits are all fixed within their lifetime, the functionality trait is expected to change throughout the lifetime of a group as it ages. We usually anticipate that within-group mutation–selection dynamics among cells will degrade group functionality as a function of age. In 1 variant of our model, we add in a cost of pleiotropy on the functioning of groups by defining (18) where the component allows us to explore what happens to the evolution of pleiotropy under the assumption that it carries a cost, ζ.

Whereas the expected growth rate of cells is implicitly defined as 1 cell division per unit time, we assume that groups have an expected life span λ>1, which gives an expected group reproduction rate of ρ = 1/λ. The larger λ is, the more cell divisions a group will experience within its life span. We also consider that group reproduction might depend on reproductive maturity, achieved at a fraction of the expected life span, α∈[0, 1]. For example, if α = 0.5 and λ = 50, then the age of group reproductive maturity will be 25 and groups will only be able to reproduce after that point. The reproductive maturity status of a k-type group aged y is denoted and defined by (19)

The functionality trait and reproductive maturity trait influence the birth and death rates of groups within the group population. Specifically, the between-group birth and death rates of k-type groups aged y at time t are given by (20) and (21) respectively, where sg is the strength of selection on group traits, is the fraction of multicellular groups in the group population that have reached reproductive maturity at time t, and is the average group functionality in the population at time t. Note that groups with a higher level of functionality gain a survival advantage over groups with a lower level of functionality. Thus, between-group selection favours maximal expression of both cooperative and private traits among cells.

An ancestral group will sometimes produce a descendant group with a different set of founding cell genotypes than it started life with. By default we assume that a descendant group is formed by choosing a founding cell at random from the ancestral group at the time at which it reproduced. Under these circumstances, the conditional transition probability that a k-type ancestral group produces an l-type descendant group, given that it has reproduced, is given by (22) where is the relative frequency of l-type cells in a k-type group aged y and is the conditional transition probability that a k-type cell mutates to an l-type cell during the production of the founding cell of a new group. Of course, is nonzero and so the most likely scenario is that a cell does not mutate. Note that and are the absolute abundance of l-type cells in a k-type group and the size of a k-type group aged y, respectively.

In our default scenario, the production of mutant descendant groups tends to become more likely as a group ages, because mutation–selection dynamics among cells degrade its clonality. In a later extension of the model, we consider the alternative possibility that groups might have a germ line. To model this scenario, we assume that with a probability γ, groups are founded by a cell with the same genotype as the ancestral group (except for rare mutations in germ line cells), and with probability 1−γ, they are founded by cells chosen at random with probability proportional to their frequency, as in our default scenario. The conditional transition probability that a k-type ancestral group produces an l-type descendant group, given that it has reproduced, in this scenario is given by (23)

Supporting information

S1 Text. Revisiting the analysis of dos Santos and colleagues [29].

This Supporting information goes through a previous study on the evolution of pleiotropy and concludes that a different modelling approach is needed.

https://doi.org/10.1371/journal.pbio.3001626.s001

(DOCX)

S1 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype i = 7.

Within-group mutation selection dynamics are shown for a group founded by a cell with genotype g7, which actively expresses cooperation, , a private trait, , but no pleiotropy, . Growth of the group as its age, y, increases, is logistic, with a carrying capacity K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, ϕ. The vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 25% of the group. Note that the strength of pleiotropy has no effect on the within-group dynamics. (A) Changes in genotype abundances, nc(y). (B) Changes in genotype relative frequencies, xc(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, . (D) Differential fitness effects of loss-of-function mutations within the group. For example, the blue bar indicates the fitness effect of a loss-of-function mutation in the cooperative trait. Parameters: sc = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s002

(DOCX)

S2 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype i = 6.

Within-group mutation selection dynamics are shown for a group founded by a cell with genotype g6, which actively expresses cooperation, , pleiotropy, , but no private trait, . Growth of the group as its age, y, increases, is logistic, with a carrying capacity K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, ϕ. The vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 25% of the group. Note that the strength of pleiotropy has no effect on the within-group dynamics. (A) Changes in genotype abundances, nc(y). (B) Changes in genotype relative frequencies, xc(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, . (D) Differential fitness effects of loss-of-function mutations within the group. Parameters: sc = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s003

(DOCX)

S3 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype i = 5.

Within-group mutation selection dynamics are shown for a group founded by a cell with genotype g5, which actively expresses cooperation, , but no private trait, , and no pleiotropy, . Growth of the group as its age, y, increases, is logistic, with a carrying capacity K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, ϕ. The vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 25% of the group. Note that the strength of pleiotropy has no effect on the within-group dynamics. (A) Changes in genotype abundances, nc(y). (B) Changes in genotype relative frequencies, xc(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, . (D) Differential fitness effects of loss-of-function mutations within the group. Parameters: sc = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s004

(DOCX)

S4 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype i = 4.

Within-group mutation selection dynamics are shown for a group founded by a cell with genotype g4, which actively expresses a private trait, , and pleiotropy, , but no cooperation, . Growth of the group as its age, y, increases, is logistic, with a carrying capacity K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, ϕ. The vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 25% of the group. Note that the strength of pleiotropy has no effect on the within-group dynamics. (A) Changes in genotype abundances, nc(y). (B) Changes in genotype relative frequencies, xc(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, . (D) Differential fitness effects of loss-of-function mutations within the group. Parameters: sc = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s005

(DOCX)

S5 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype i = 3.

Within-group mutation selection dynamics are shown for a group founded by a cell with genotype g4, which actively expresses a private trait, , but no cooperation, , or pleiotropy, . Growth of the group as its age, y, increases, is logistic, with a carrying capacity K =200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, ϕ. The vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 25% of the group. Note that the strength of pleiotropy has no effect on the within-group dynamics. (A) Changes in genotype abundances, nc(y). (B) Changes in genotype relative frequencies, xc(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, . (D) Differential fitness effects of loss-of-function mutations within the group. Parameters: sc = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s006

(DOCX)

S6 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype i = 2.

Within-group mutation selection dynamics are shown for a group founded by a cell with genotype g2, which actively expresses pleiotropy, , but no cooperation, , a no private trait, . Growth of the group as its age, y, increases, is logistic, with a carrying capacity K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, ϕ. The vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 25% of the group. Note that the strength of pleiotropy has no effect on the within-group dynamics. (A) Changes in genotype abundances, nc(y). (B) Changes in genotype relative frequencies, xc(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, . (D) Differential fitness effects of loss-of-function mutations within the group. Parameters: sc = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s007

(DOCX)

S7 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype i = 1.

Within-group mutation selection dynamics are shown for a group founded by a cell with genotype g1, which does not actively expresses any of cooperation, , a private trait, , or pleiotropy, . Growth of the group as its age, y, increases, is logistic, with a carrying capacity K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, ϕ. The vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 25% of the group. Note that the strength of pleiotropy has no effect on the within-group dynamics. (A) Changes in genotype abundances, nc(y). (B) Changes in genotype relative frequencies, xc(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, . The differential fitness effects of loss-of-function mutations are not shown, because g1 cells have no traits to lose. Parameters: sc = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s008

(DOCX)

S8 Fig. Long-term evolutionary dynamics to illustrate a case where the evolution of cooperation depends entirely on pleiotropy.

Dynamics are shown for the global population of cells over time, t, which encompasses many generations of cell groups. These dynamics, therefore, encompass both within-group and between-group selection dynamics. We show 3 strengths of pleiotropy, ϕ, to capture 3 qualitatively different scenarios. Stronger pleiotropy is associated with more rapid and complete evolution of both pleiotropy and cooperation. (A) Changes in global genotype relative frequencies, xc(t). (B) Changes in the global average levels of cooperation, private trait expression, and pleiotropy, . In the absence of pleiotropy ϕ = 0, cooperation fails to evolve. (C) Average change in traits over a group lifetime, measured as the difference between the average trait values among groups aged y to those expected from their founding cell at birth. Parameters: sc = sg = 0.95; K = 500; μ = 0.0001; ν = 0.01; λ = 25. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s009

(DOCX)

S9 Fig. Private trait evolution does not favour pleiotropy and vice versa.

We explored a model in which the cooperative trait was replaced by another private trait. Heatmaps show average trait values of 2 private traits and pleiotropy among the global population of cells (across all groups) at steady state in our revised model. Results are shown for 3 group sizes (increasing from top to bottom). Both private traits evolve to fixation under all parameter values, but pleiotropy is never favoured. The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s010

(DOCX)

S10 Fig. Within-group evolutionary dynamics of groups founded by 8 cells with all 8 genotypes represented.

Within-group mutation selection dynamics are shown for a group founded by cells with all possible genotypes and phenotypes represented. Growth of the group as its age, y, increases, is logistic, with a carrying capacity K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, ϕ. Note that the strength of pleiotropy has no effect on the within-group dynamics. (A) Changes in genotype abundances, nc(y). (B) Changes in genotype relative frequencies, xc(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, . (D) Differential fitness effects of loss-of-function mutations within the group. Parameters: sc = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s011

(DOCX)

S11 Fig. Pleiotropy stabilises cooperation when groups are founded by 2 cells rather than one.

We explored a model in which 2 cells are selected uniformly at random from the ancestor group to found each descendant group (rather than a single cell, which we assume in the main paper). This assumption lowers the expected relatedness in groups at the point at which they form. Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 group sizes (increasing from top to bottom). As when groups are founded by a single cell, pleiotropy is favoured when the strength of pleiotropy, ϕ, is higher. The overall levels of cooperation are lower when groups are founded by 2 cells rather than one, but the evolution of pleiotropy still promotes the evolution of cooperation. The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = sg = 0.95; K = 200; μ = 0.0001; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s012

(DOCX)

S12 Fig. Pleiotropy evolves to stabilise cooperation even when costly.

We explored a model in which pleiotropy carries a cost on multicellular function. We assumed multicellular group function is reduced by a factor , where ζ is the cost of pleiotropy, and is the average pleiotropy in a group. Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 levels of cost (increasing from top to bottom). Pleiotropy evolves even when it carries a 1%, 2%, or 5% cost to group functionality. When costly, it is especially favoured when groups are longer lived and therefore require mechanisms to limit the spread of noncooperative mutant lineages. The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = sg = 0.95; K = 200; μ = 0.0001; ν = 0.01; K = 200. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s013

(DOCX)

S13 Fig. Pleiotropy evolves to stabilise cooperation even when group selection is weak.

We explored a model in which the strength of group selection, sg, is varied. Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 strengths of group selection (increasing from top to bottom). Stronger group selection favours both cooperation and pleiotropy, but strong pleiotropy can help cooperation evolve even when group selection is weak, (sg = 0.05). The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = 0.95; K = 200; μ = 0.0001; ν = 0.01; K = 200. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s014

(DOCX)

S14 Fig. The rate of gain-of-function mutations has little impact on the dynamics of pleiotropy evolution.

We varied the relative rate of gain-of-function mutations to loss-of-function mutations, ν. Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 gain-of-function rates (increasing from top to bottom). The gain-of-function ratio had a marginal effect on the evolution of pleiotropy, favouring slightly higher rates when ν is higher. The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = sg = 0.95; K = 200; μ = 0.0001; K = 200. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s015

(DOCX)

S15 Fig. Strong pleiotropy can help rescue cooperation even in the face of high mutation rates.

We varied the loss-of-function mutation rate, μ. Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 loss-of-function rates (increasing from top to bottom). Higher mutation rates disfavoured the evolution of cooperation, but pleiotropy still evolved at higher strengths of pleiotropy and for lower group sizes. The evolution of pleiotropy was associated with a stabilisation of cooperation even at higher mutation rates. The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = sg = 0.95; K = 200; ν = 0.01. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s016

(DOCX)

S16 Fig. Lower mutation rates increases cooperation but pleiotropy still evolves.

We varied the loss-of-function mutation rate, μ, using values 10-fold lower than the typical value used in the main text. Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 loss-of-function rates (increasing from top to bottom). With 10-fold lower mutation rates, we found that cooperation breaks down only for longer group life spans, λ. When cooperation is threatened by breakdown, however, we still found that pleiotropy evolves to rescue it. The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = sg = 0.95; K = 200; ν = 0.01; c = 0. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s017

(DOCX)

S17 Fig. Lower mutation rates increase cooperation but pleiotropy still evolves, even when pleiotropy is costly.

We varied the gain-of-function mutation rate, μ, using values 10-fold lower than the typical value used in the main text, and assumed that the evolution of pleiotropy decreases group function by 2%. Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 loss-of-function rates (increasing from top to bottom). With 10-fold lower mutation rates and a cost, we found that pleiotropy evolves most frequently when the mutation rate is higher and cooperation is under the greatest threat. The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = sg = 0.95; K = 200; ν = 0.01; ζ = 0.02. The code used to produce this figure may be found at https://github.com/euler-mab/pleiotropy.

https://doi.org/10.1371/journal.pbio.3001626.s018

(DOCX)

S18 Fig. Reproduction at maturity hinders cooperation and has mixed effects on the evolution of pleiotropy.

We varied the age fraction of the expected life span at which groups reach reproductive maturity, α, and assumed that the evolution of pleiotropy decreases group function by 2%. Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 reproductive maturity parameters (increasing the age at which maturity is reached from top to bottom). Increasing the age of reproductive maturity has a marginal positive effect on the evolution of pleiotropy for shorter life spans, but a marginal negative effect on cooperation. The dotted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: sc = sg = 0.95; K = 200; ν = 0.01; ζ = 0.02. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s019

(DOCX)

S19 Fig. When mutation-driven breakdown of cooperation is a problem, pleiotropy stabilises cooperation.

We modified the individual-based model of dos Santos and colleagues to compare what happens when pleiotropy cannot evolve (left) as compared to when it can evolve (right). We varied the length of the within-group growth phase k during which spontaneous mutants can arise and invade within groups, and the relatedness r at the point at which groups form. We follow the evolution of a private trait, cooperative trait, and pleiotropy trait. (A) Evolutionary dynamics of all 3 traits for a within-group growth phase of k = 30 and r = 1. (B) Steady-state levels of all 3 traits under a when the length of the within-group growth phase is varied (x-axis) for r = 1. (C) Steady-state levels of all 3 traits when relatedness at the point at which groups form is varied for a within-group growth phase of k = 30. Other parameters: b = 0.11, c = 0.1, g = 0.5, mutation rate μ = 0.001, number of groups ng = 1000. All plots are averages of 10 replicates. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s020

(DOCX)

S20 Fig. Pleiotropy evolves to stabilise cooperation even when the loss-of-function mutation rate is 10× greater than the gain-of-function rate.

We performed simulations using the mutation model used by dos Santos and colleagues in the production of their S5 Fig. This again reveals conditions where pleiotropy evolves and increases cooperation as it does. We varied the length of the within-group growth phase k during which spontaneous mutants can arise and invade within groups, and the relatedness r at the point at which groups form. We follow the evolution of a private trait, cooperative trait, and pleiotropy trait. (A) Evolutionary dynamics of all 3 traits for a within-group growth phase of k = 30 and r = 1. (B) Steady-state levels of all 3 traits under a when the length of the within-group growth phase is varied (x-axis) for r = 1. (C) Steady-state levels of all 3 traits when relatedness at the point at which groups form is varied for a within-group growth phase of k = 30. Other parameters: b = 0.11, c = 0.1, g = 0.5, mutation rate μ = 0.001, number of groups ng = 1000. All plots are averages of 10 replicates. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s021

(DOCX)

S21 Fig. Pleiotropy evolves to stabilise cooperation even when the loss-of-function mutation rate is 100× greater than the gain-of-function rate.

We performed simulations using the mutation model used by dos Santos and colleagues in the production of their S6 Fig. This again reveals conditions where pleiotropy evolves and increases cooperation as it does. We varied the length of the within-group growth phase k during which spontaneous mutants can arise and invade within groups, and the relatedness r at the point at which groups form. We follow the evolution of a private trait, cooperative trait, and pleiotropy trait. (A) Evolutionary dynamics of all 3 traits for a within-group growth phase of k = 30 and r = 1. (B) Steady-state levels of all 3 traits under a when the length of the within-group growth phase is varied (x-axis) for r = 1. (C) Steady-state levels of all 3 traits when relatedness at the point at which groups form is varied for a within-group growth phase of k = 30. Other parameters: b = 0.11, c = 0.1, g = 0.5, mutation rate μ = 0.001, number of groups ng = 1000. All plots are averages of 10 replicates. The code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk.

https://doi.org/10.1371/journal.pbio.3001626.s022

(DOCX)

Acknowledgments

With thanks to Dave Queller and Arvid Agren for helpful comments on an earlier version of the main text, to Stu West and Miguel dos Santos for helping us to better understand their own model of pleiotropy and to Andrew Murray for a formative conversation on the evolution of pleiotropy.

References

  1. 1. Nadell CD, Drescher K, Foster KR. Spatial structure, cooperation and competition in biofilms. Nat Rev Microbiol. 2016;14:589–600. pmid:27452230
  2. 2. Fisher RM, Cornwallis CK, West SA. Group formation, relatedness, and the evolution of multicellularity. Curr Biol. 2013;23:1120–5. pmid:23746639
  3. 3. Ghoul M, Griffin AS, West SA. Toward an evolutionary definition of cheating. Evolution. 2014. pmid:24131102
  4. 4. Greig D, Travisano M. The Prisoner’s Dilemma and polymorphism in yeast SUC genes. Proc R Soc B Biol Sci. 2004;271:25–6. pmid:15101409
  5. 5. Diggle SP, Griffin AS, Campbell GS, West S, a. Cooperation and conflict in quorum-sensing bacterial populations. Nature. 2007;450:411–4. pmid:18004383
  6. 6. Harrison F, Browning LE, Vos M, Buckling A. Cooperation and virulence in acute Pseudomonas aeruginosa infections. BMC Biol. 2006;4:1–5.
  7. 7. Griffin AS, West SA, Buckling A. Cooperation and competition in pathogenic bacteria. Nature. 2004;430:1024–7. pmid:15329720
  8. 8. Sandoz KM, Mitzimberg SM, Schuster M. Social cheating in Pseudomonas aeruginosa quorum sensing. Proc Natl Acad Sci U S A. 2007;104:15876–81. pmid:17898171
  9. 9. Velicer GJ, Kroos L, Lenski RE. Developmental cheating in the social bacterium Myxococcus xanthus. Nature. 2000;404:598–601. pmid:10766241
  10. 10. Drescher K, Nadell CD, Stone HA, Wingreen NS, Bassler BL. Solutions to the Public Goods Dilemma in Bacterial Biofilms. Curr Biol. 2014;24:50–5. pmid:24332540
  11. 11. Andersen SB, Marvig RL, Molin S, Johansen HK, Griffin AS. Long-term social dynamics drive loss of function in pathogenic bacteria. Proc Natl Acad Sci U S A. 2015;112:10756–61. pmid:26240352
  12. 12. Raymond B, West S, a, Griffin AS, Bonsall MB. The dynamics of cooperative bacterial virulence in the field. Science. 2012;337:85–8. pmid:22767928
  13. 13. Cordero OX, Ventouras L-A, DeLong EF, Polz MF. Public good dynamics drive evolution of iron acquisition strategies in natural bacterioplankton populations. Proc Natl Acad Sci U S A. 2012;109:20059–64. pmid:23169633
  14. 14. Grosberg RK, Strathmann RR. The Evolution of Multicellularity: A Minor Major Transition? Annu Rev Ecol Evol Syst. 2007;38:621–54.
  15. 15. Grosberg RK, Strathmann RR. One cell, two cell, red cell, blue cell: The persistence of a unicellular stage in multicellular life histories. Trends Ecol Evol. 1998;13:112–6. pmid:21238226
  16. 16. Hamilton WD. The genetical evolution of social behaviour. I & II J Theor Biol. 1964;7:1–52.
  17. 17. Foster KR. A defense of sociobiology. Cold Spring Harb Symp Quant Biol. 2009;74:403–18. pmid:20375319
  18. 18. Michod RE. Evolution of the individual. Am Nat. 1997;150:S5–s21. pmid:18811312
  19. 19. Michod RE, Roze D. Cooperation and conflict in the evolution of multicellularity. Heredity. 2001:1–7. pmid:11298810
  20. 20. Michod RE. Darwinian dynamics: evolutionary transitions in fitness and individuality. Princeton University Press; 2000.
  21. 21. Dandekar AA, Chugani S, Greenberg EP. Bacterial quorum sensing and metabolic incentives to cooperate. Science (80-). 2012;338:264–6. pmid:23066081
  22. 22. Wang M, Schaefer AL, Dandekar AA, Greenberg EP. Quorum sensing and policing of Pseudomonas aeruginosa social cheaters. Proc Natl Acad Sci U S A. 2015;112:2187–91. pmid:25646454
  23. 23. Sathe S, Mathew A, Agnoli K, Eberl L, Kümmerli R. Genetic architecture constrains exploitation of siderophore cooperation in the bacterium Burkholderia cenocepacia. Evol Lett. 2019;3:610–22. pmid:31844554
  24. 24. Foster KR, Shaulsky G, Strassmann JE, Queller DC, Thompson CRL. Pleiotropy as a mechanism to stabilize cooperation. Nature. 2004;431:693–6. pmid:15470429
  25. 25. Foster KR, Parkinson K, Thompson CRL. What can microbial genetics teach sociobiology? Trends Genet. 2007;23:74–80. pmid:17207887
  26. 26. Frénoy A, Taddei F, Misevic D. Genetic architecture promotes the evolution and maintenance of cooperation. PLoS Comput Biol. 2013;9:e1003339. pmid:24278000
  27. 27. Chisholm RH, Connelly BD, Kerr B, Tanaka MM. The role of pleiotropy in the evolutionary maintenance of positive niche construction. Am Nat. 2018;192:35–48. pmid:29897798
  28. 28. Sharp C, Foster KR. 2022. Host control and the evolution of cooperation in host microbiomes. Nature Communications, in press.
  29. 29. dos Santos M, Ghoul M, West SA. Pleiotropy, cooperation, and the social evolution of genetic architecture. PLoS Biol. 2018;16:1–25. pmid:30359363
  30. 30. David CN, Campbell RD. Cell Cycle Kinetics and Development of Hydra Attenuata: I. Epithelial Cells J Cell Sci. 1972;11:557–68. pmid:5076361
  31. 31. Foster KR. The sociobiology of molecular systems. Nat Rev Genet. 2011;12:193–203. pmid:21301472
  32. 32. Grafen A. Natural selection, kin selection and group selection. In: Krebs JR, Davies NB, editors. Behavioural ecology: An evolutionary approach. Oxford: Blackwell; 1984. p. 62–84.
  33. 33. Paaby AB, Rockman MV. The many faces of pleiotropy. Trends Genet. 2013;29:66–73. pmid:23140989
  34. 34. Stearns FW. One hundred years of pleiotropy: a retrospective. Genetics. 2010;186:767–73. pmid:21062962
  35. 35. Guillaume F, Otto SP. Gene functional trade-offs and the evolution of pleiotropy. Genetics. 2012;192:1389–409. pmid:22982578
  36. 36. Classon M, Harlow E. The retinoblastoma tumour suppressor in development and cancer. Nat Rev Cancer. 2002;2:910–7. pmid:12459729
  37. 37. Hickman ES, Moroni MC, Helin K. The role of p53 and pRB in apoptosis and cancer. Curr Opin Genet Dev. 2002;12:60–6. pmid:11790556
  38. 38. Evan G, Lowe S, Cepero E. Intrinsic tumour suppression. Nature. 2004;432:1–9. pmid:15549092
  39. 39. Wagner GP, Zhang J. The pleiotropic structure of the genotype-phenotype map: The evolvability of complex organisms. Nat Rev Genet. 2011;12:204–13. pmid:21331091
  40. 40. Maynard-Smith J. Group selection and kin selection. Nature. 1964;201:1145–7.
  41. 41. Orr HA. Adaptation and the cost of complexity. Evolution (N Y). 2000;54:13–20. pmid:10937178
  42. 42. Milholland B, Dong X, Zhang L, Hao X, Suh Y, Vijg J. Differences between germline and somatic mutation rates in humans and mice. Nat Commun. 2017;8:1–8.
  43. 43. Werner B, Sottoriva A. Variation of mutational burden in healthy human tissues suggests non-random strand segregation and allows measuring somatic mutation rates. PLoS Comput Biol. 2018;14:e1006233. pmid:29879111
  44. 44. Tomczyk S, Fischer K, Austad S, Galliot B. Hydra, a powerful model for aging studies. Invertebr Reprod Dev. 2015;59:11–6. pmid:26120246
  45. 45. Alivisatos AP, Chun M, Church GM, Greenspan RJ, Roukes ML, Yuste R. The Brain Activity Map Project and the Challenge of Functional Connectomics. Neuron. 2012;74:970–4. pmid:22726828
  46. 46. Kim W, Racimo F, Schluter J, Levy SB, Foster KR. Importance of positioning for microbial evolution. Proc Natl Acad Sci U S A. 2014;111.
  47. 47. Bianconi E, Piovesan A, Facchin F, Beraudi A, Casadei R, Frabetti F, et al. An estimation of the number of cells in the human body. Ann Hum Biol. 2013;40:463–71. pmid:23829164
  48. 48. Mitri S, Foster KR. Pleiotropy and the low cost of individual traits promote cooperation. Evolution (N Y). 2016;70:488–94. pmid:26748821
  49. 49. Queller DC. Pleiotropy and synergistic cooperation. PLoS Biol. 2019:3–5. pmid:31226105
  50. 50. Korb J, Foster KR. Ecological competition favours cooperation in termite societies. Ecol Lett. 2010;13:754–60. pmid:20584170
  51. 51. Rankin DJ, López-Sepulcre A, Foster KR, Kokko H. Species-level selection reduces selfishness through competitive exclusion. J Evol Biol. 2007;20:1459–68. pmid:17584239
  52. 52. Ostrowski EA. Enforcing Cooperation in the Social Amoebae. Curr Biol. 2019;29:R474–84. pmid:31163162
  53. 53. Wielgoss S, Schneider D, Barrick JE, Tenaillon O, Cruveiller S, Chane-Woon-Ming B, et al. Mutation rate inferred from synonymous substitutions in a long-term evolution experiment with escherichia coli. G3 Genes, Genomes Genet. 2011;1:183–6. pmid:22207905
  54. 54. Crespi BJ. The evolution of social behavior in microorganisms. Trends Ecol Evol. 2001;16:178–83. pmid:11245940
  55. 55. Harrison F, Buckling A. Hypermutability impedes cooperation in pathogenic bacteria. Curr Biol. 2005;15:1968–71. pmid:16271876
  56. 56. Azimi S, Roberts AEL, Peng S, Weitz JS, McNally A, Brown SP, et al. Allelic polymorphism shapes community function in evolving Pseudomonas aeruginosa populations. ISME J. 2020;14:1929–42. pmid:32341475
  57. 57. Taylor TB, Johnson LJ, Jackson RW, Brockhurst MA, Dash PR. First steps in experimental cancer evolution. Evol Appl. 2013;6:535–48. pmid:23745144
  58. 58. Lean C, Plutynski A. The evolution of failure: explaining cancer as an evolutionary process. Biol Philos. 2016;31:39–57.
  59. 59. Nunney L. Lineage selection and the evolution of multistage carcinogenesis. Proc R Soc B Biol Sci. 1999;266:493–8. pmid:10189713
  60. 60. Nunney L. The real war on cancer: the evolutionary dynamics of cancer suppression. Evol Appl. 2013;6:11–9. pmid:23396311
  61. 61. Singh M, Boomsma JJ. Policing and punishment across the domains of social evolution. Oikos. 2015;124:971–82.
  62. 62. DeGregori J. Adaptive oncogenesis: a new understanding of how cancer evolves inside us. Harvard University Press; 2018.
  63. 63. DeGregori J. Challenging the axiom: does the occurrence of oncogenic mutations truly limit cancer development with age? Oncogene. 2013;32:1869–75. pmid:22751134
  64. 64. Serrano M, Lin AW, McCurrach ME, Beach D, Lowe SW. Oncogenic ras provokes premature cell senescence associated with accumulation of p53 and p16INK4a. Cell. 1997;88:593–602. pmid:9054499
  65. 65. Schreiber RD, Old LJ, Smyth MJ. Cancer immunoediting: integrating immunity’s roles in cancer suppression and promotion. Science (80-). 2011;331:1565–70. pmid:21436444
  66. 66. Alpar L, Bergantiños C, Johnston LA. Spatially restricted regulation of Spätzle/Toll signaling during cell competition. Dev Cell. 2018;46:706–19. pmid:30146479
  67. 67. Russell JH, Ley TJ. Lymphocyte-mediated cytotoxicity. Annu Rev Immunol. 2002;20:323–70. pmid:11861606
  68. 68. Kajita M, Fujita Y. EDAC: Epithelial defence against cancer—cell competition between normal and transformed epithelial cells in mammals. J Biochem. 2015;158:15–23. pmid:25991731
  69. 69. C. QD, Eleonora P, Salvatore B, E. SJ. Single-Gene Greenbeard Effects in the Social Amoeba Dictyostelium discoideum. Science (80-). 2003;299: 105–106. pmid:12511650
  70. 70. Schluter J, Schoech AP, Foster KR, Mitri S. The evolution of quorum sensing as a mechanism to infer kinship. PLoS Comput Biol. 2016;12:e1004848. pmid:27120081
  71. 71. Xavier JB, Kim W, Foster KR. A molecular mechanism that stabilizes cooperative secretions in Pseudomonas aeruginosa. Mol Microbiol. 2011;79:166–79. pmid:21166901
  72. 72. Simon B, Fletcher JA, Doebeli M. Towards a General Theory of Group Selection; 2012.
  73. 73. Gillespie DT. Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977;81:2340–61.
  74. 74. LeVeque RJ. Finite volume methods for hyperbolic problems. Cambridge university press; 2002.
  75. 75. Xavier JC, Patil KR, Rocha I. Metabolic models and gene essentiality data reveal essential and conserved metabolism in prokaryotes. PLoS Comput Biol. 2018;14:e1006556. pmid:30444863