## Abstract

Coevolution of two species is typically thought to favour the evolution of faster evolutionary rates helping a species keep ahead in the Red Queen race, where ‘it takes all the running you can do to stay where you are’. In contrast, if species are in a mutualistic relationship, it was proposed that the Red King effect may act, where it can be beneficial to evolve slower than the mutualistic species. The Red King hypothesis proposes that the species which evolves slower can gain a larger share of the benefits. However, the interactions between the two species may involve multiple individuals. To analyse such a situation, we resort to evolutionary multiplayer games. Even in situations where evolving slower is beneficial in a two-player setting, faster evolution may be favoured in a multiplayer setting. The underlying features of multiplayer games can be crucial for the distribution of benefits. They also suggest a link between the evolution of the rate of evolution and group size.

## 1. Introduction

Mutualistic relationships, interspecific interactions that benefit both species, have been empirically studied for many years [1–7] and also a considerable body of theory has been put forth trying to explain the evolution and maintenance of such relationships [8–15]. Many of these studies use evolutionary game theory for developing the models. The interactions in these models are usually pairwise. A representative of each species is chosen and the outcome of the interactions between these representatives determines the evolutionary dynamics within each of the two species. However, in many cases interactions between species cannot be reduced to such pairwise encounters [16].

For example, in the interaction between ants and aphids or butterfly larvae [17–19] many ants tend to these soft-bodied creatures, providing them with shelter and protection from predation and parasites in exchange for honeydew, a rich source of food for the ants [16,20]. This is not a one-to-one interaction between a larva and an ant, but rather a one-to-many interaction from the perspective of the larva. In this manuscript, we focus on this kind of—possibly—many-to-many interactions between two mutualistic species.

To analyse how the benefits are shared between the two mutualistic species, we make use of evolutionary game theory [21–23]. Following Bergstrom & Lachmann [12], we analyse the interactions between two species with a twist to the standard formulation. The two interacting species can have different evolutionary rates. In coevolutionary arms races, where species are locked in antagonistic relationships such as host–parasite interactions, we observe the Red Queen effect [24]. In these cases, a higher rate of evolution will be beneficial, as it would help the parasite infecting a host or a host escaping from a parasite. However, in mutualistic relationships a slower rate of evolution was predicted to be more favourable [12]. Hence, it is possible that the type of relationship between two species affects the evolution of the rate of evolution. As mentioned in [12,25] the different evolutionary rates could be owing to a multitude of factors ranging from different population sizes to the differing amount of segregating genetic variance. The implications of a difference in evolutionary rates are not limited to mutualism and antagonistic relationships. Epidemiological modelling and data have shown correlations between the rate of epidemic spreading and the evolutionary rate of the spreading pathogen [26,27].

We first recall the mutualistic relationship between two species in a two-player game, as proposed by Bergstrom and Lachmann. Then, we increase the number of players. Note that we do not increase the number of interacting species [28,29], but rather the number of interacting individuals between two species (see [30]). We include asymmetry in evolutionary rates and discuss its effect both in two-player and in multiplayer games. We find that when in a two-player setting it is beneficial to evolve at a slower rate, it can be detrimental in a multiplayer game.

## 2. Model and results

To set the stage, we first recapitulate the two-player model presented in [12] albeit with different notation. They considered two species each with two strategies, ‘Generous’ and ‘Selfish’. Each species is better off being Selfish as long as the other one is Generous. If both are Selfish, then no mutualistic benefit is generated and hence in that case it is better to be Generous.

Under these assumptions, the payoff matrices for the interactions describing the interactions of each type with one member of the other species are

species 2 | species 1 | ||||||

G_{2} | S_{2} | G_{1} | S_{1} | ||||

species 1 | G_{1} | species 2 | G_{2} | ||||

S_{1} | S_{2} |

where, for example, a Generous member of species 1 obtains from an interaction with a Selfish member of species 2, whereas the latter obtains . In our case, we have and for *i*,*j* = 1,2. This ordering of payoffs corresponds to a snowdrift game [31] (see the electronic supplementary material) where there exists a point where the two strategies can coexist. This is because if both the players are Generous, then one can get away with being Selfish, but if both are Selfish, then it actually pays to be Generous and chip in. For a snowdrift game in a single species the coexistence point is stable, i.e. deviations from this point bring the system back to equilibrium, because the deviators would always be disfavoured by selection. However, for a snowdrift game between two species, this coexistence point is unstable as each species would be better off being Selfish, i.e. exploiting the deviation of the other species.

The frequency of players playing strategy Generous (*G*_{1}) in species 1 is given by *x* and in species 2 (*G*_{2}) by *y*. The frequencies of players playing strategy Selfish (*S*_{1} and *S*_{2}) are given by 1 − *x* and 1 − *y* in species 1 and 2, respectively. The fitness of the Generous strategy in species 1, depends on the frequency *y* of Generous players in species 2, . Equivalently, the fitness of the Generous strategy in species 2, depends on the frequency *x* of Generous players in species 1, . The replicator dynamics assumes that the change in frequency of a strategy is proportional to the difference between the fitness of that strategy and the average fitness of the species [23,32,33]. Thus, the time evolution of the frequencies of the Generous players in the two species are
2.1

The parameters *r*_{x} and *r*_{y} are the evolutionary rates of the two species. We first recover the scenarios described in [12]. If the evolutionary rates are equal (*r*_{x} = *r*_{y}) and the evolutionary game is symmetric, then the basins of attraction of (*S*_{1}, *G*_{2}) and (*G*_{1}, *S*_{2}) are of equal size (figure 1*a*). For unequal evolutionary rates, the species which is evolving slower (in our case species 1 with the rate *r*_{x} = *r*_{y}/8) has a larger basin of attraction (figure 1*b*). This asymmetry where most of the initial conditions lead to an outcome favouring the slower evolving species has been termed as the *Red King effect* [12].

We now extend the above approach to multiplayer games. Extending the number of players from 2 to *d* adds a polynomial nonlinearity to the fitness functions of the strategies (see the electronic supplementary material). For a multiplayer game, we no longer have a 2 × 2 payoff matrix, but rather a payoff table. For this, we use the proposal from [34] for a *d*-player snowdrift game, where the costs of being Generous are divided among the Generous players. In addition, only if there are at least *M* Generous players, a benefit is produced. That is, for *k* < *M* Generous players, each one of them exhibits a loss of c/*M* and the Selfish players obtain nothing. If there are at least *M* Generous players, then a benefit is produced. The Generous obtain *b*−*c*/k and the Selfish obtain *b* at no cost. For species 2 we can write down a different payoff setup which could have different values for *b*, *c*, *M*, *d*, etc. thus creating a ‘bi-table’ game. For the time being, we assume that the payoff setups are symmetric for the two species and hence we just elucidate the details for species 1. The exact formulation of the payoffs and the calculations of fitness values are given in the electronic supplementary material.

Note that for *d* = 2, *M* = 1, *b* = 2 and *c* = 1, we recover the matrix used in [12]. Even for these new fitness functions for multiplayer games, the dynamics are still given by the replicator equations [35–37]. Also note that for two-player games with *M* = 1, there are four fixed points in which each species is either Selfish or Generous. In addition, there is an internal fixed point given by . The position and the stability of the fixed points is independent of the evolutionary rates, but as we will see, it depends on the number of players *d*. For a 20-player game (*d* = 20), the basins of attractions are still of the same size, but the dynamics leading to the stable points on the vertices are completely different (figure 1*c*, *r*_{x} = *r*_{y}). The internal equilibrium has now shifted to *x* = *y* = 0.063. As before, we introduce an asymmetry in the evolutionary rates. Interestingly, we find that for a 20-player game (figure 1*d*, *r*_{x} = *r*_{y} / 8) for the same asymmetric values of growth rates as in the two-player case, most of the initial conditions lead to a stable point where species 2 is Selfish and species 1 is Generous (*G*_{1}, *S*_{2}). This reverses the result which we got from *d* = 2. Everything else being the same, in the presence of multiple players, the Red King effect is not observed in this example.

Next, we explore the process owing to which the Red King effect vanishes. The replicator solutions of the two species creates quadrants in the state space (0 ≤ *x*, *y* ≤ 1). Of these quadrants, the top right and the bottom left are of special interest as they contain the curve that separates the two basins of attraction (the blue and red sections in figure 1). Hence all points starting on one side of that curve lead to the same equilibrium. Consider the top-right quadrant. Species 2 is represented by the *y*-axis. A faster evolution by species 2 results in most of the initial conditions leading to the outcome favourable for species 2, i.e. (*G*_{1}, *S*_{2}). An exactly opposite scenario is taking place in the bottom-left quadrant. Hence in this quadrant as species 1 is evolving at a slower pace than species 2, most of the initial conditions here lead to an outcome favouring species 1, i.e. (*S*_{1}, *G*_{2}). As long as the internal equilibrium is on the diagonal, the Red King effect depends on the sizes of these quadrants. Changing the number of players alters the sizes of these two influential quadrants. For example, consider the case of the 20-player game (figure 1*c,d*). The size of the bottom-left quadrant is reduced to such an extent that almost the whole state space leads to the outcome favourable for the faster evolving species.

The bottom-left and the top-right quadrants have equal size when the internal equilibrium is at *x* = *y* = 0.5. For a fixed *b* and *c*, we cannot select any arbitrary number of players *d* to obtain this equilibrium, as *d* is not a continuous variable. If the equilibrium is above *x* = 0.5, then a decrease in the evolutionary rate can be beneficial as demonstrated by the Red King effect. Conversely, if the equilibrium is below *x* = 0.5, then an increase in the evolutionary rate might be favourable. Hence, if the number of players leads to an equilibrium at *x* < 0.5, then the faster evolving species would be favoured.

Asymmetries have been considered in mutualistic species at the level of species or other properties of the system such as interactions [38], interaction lengths [11] or growth rates [12]. Asymmetry in the number of interacting partners has only been recently tackled [30]. Going back to the example of ants and larvae, a single larva is tended to by multiple ants. Thus, while from each ant's point of view this is a two-player game, for the larva this would be a multiplayer game. Where such multiplayer games are feasible, it is also possible that a certain quorum needs to be fulfilled for the game to proceed. Client fish have been shown to choose cleaning stations with two cleaners over solitary cleaners [39]. A certain number of ants are required to save a caterpillar from its predator. It has been shown that the amount of secretions of a lycaenid larva is correlated to the number of attending ants [40]. In the following two paragraphs, we explore these two points, asymmetry in the number of players for the two species and different thresholds in either species to start off the benefits of mutualistic relations.

Instead of the single parameter *d*, now we have *d*_{1} and *d*_{2} as the number of players for the two species 1 and 2. For symmetric evolutionary rates, if the two species play different player games (*d*_{1} ≠ *d*_{2}), then the basins of attraction become asymmetric. Hence, if a species is currently at a disadvantage, a modification of the number of players or the evolutionary rate may put it on equal footing with the other species. Owing to an asymmetric number of players the sizes of the basins of attraction depends not only on the sizes of the quadrants but also on the shape of the curve separating the basins of attraction. Thus, it is possible to counter the Red King effect by changing the number of interacting agents (figure 2).

Until now, we have considered that a single Generous individual can generate the benefit of mutualism (*M* = 1). To begin with the simplest multiplayer case, we consider a symmetric three-player game with different thresholds in either species to start off the benefits of mutualistic relations (say for species 1 the threshold is *M*_{1} and similarly for species 2 is *M*_{2} where in general *M*_{i} can range from 1 to *d*_{i}).

The payoff matrices become asymmetric owing to the different thresholds for the two species. Here, it matters which dynamics we are studying, the usual replicator dynamics or the modified replicator equations (figure 3 and see the electronic supplementary material) as they can result in different sizes of the basins of attraction. The choice of dynamics depends on the details of the model system under consideration. Ultimately, the macroscopic dynamics can be derived from the underlying microscopic process [41,42]. Manipulating the thresholds can also change the nature of the game from coexistence to coordination [34], which implies that different social dilemmas arise in multiplayer games.

## 3. Discussion

Interspecific relationships are exceedingly complex [43]. The development of a game theoretical approach for such multiplayer mutualisms requires an approach beyond that arising typically in multiplayer social dilemmas [44]. In a mutualistic framework, it is best for the two species to cooperate with each other. We do not question how these mutualisms arise. Rather when they do, what is the best strategy to contribute towards the common benefit [7,45]? It would be possible to include the interactions between the individuals of the same species, as has been recently explored experimentally [30]. It has also been shown [40] that the amount of larval secretions is also influenced by the quality of the other larvae in the group. But then, we would be shifting our focus from the problem of interspecific mutualism to intraspecific cooperation [15]. Here, we have focused on the interspecific interactions, where the interacting partners are always picked from the other species [46]. Bergstrom and Lachmann have shown that in such a mutualistic scenario, the species that evolves slower can get away with being Selfish and force the other species to make a Generous contribution. They termed this as the Red King effect. If we include multiple players, then the Red King effect is much more complex.

For simplicity, usually pairwise interactions are assumed in game theoretical arguments. For modelling collective phenomena [47,48], multiplayer games may be necessary. The exact number of players is a matter of choice, though. Group-size distributions give us an idea about the mean group size of a species. Instead of using pairwise interactions or an arbitrary number of individuals to form a group, we could use the mean group sizes as the number of interacting individuals. Group size is known to be of importance in mutualisms [3]. As we have seen here, it can be an influential factor in deciding how the benefits are shared. Countering the Red King or enhancing its effect is possible by altering the group size. Hence, can the group size itself be an evolving strategy? The study of group-size distributions has been tackled theoretically [49–56] and empirically in various species ranging from house sparrows to humans [48,49,57,58]. In our example of ants and butterfly larvae, it has been observed that a larva was most successful in getting more ant attendants in a group of four larvae [18]. It would be interesting to see if the distributions in mutualistic species peak at the group size which is the best response to their symbiont partners choices. This brings forth another of our assumptions also implicit in [12]. The rate at which strategies evolve is assumed to be much faster than the rate at which the evolutionary rates or as just mentioned, the group size evolve. If these traits are genetically determined, then this assumption may no longer hold. The rate of evolution is typically assumed to be constant, but it could well be a variable, subject to evolution. We have seen that the number of players can affect whether evolving slower or faster is favourable. It would then be interesting to determine the interplay between the evolving group size and the evolving evolutionary rate and what effect it has on the dynamics of strategy evolution.

Another method of introducing asymmetry is to have different payoff tables for the two species (i.e. different benefits and costs for the two species). Also we have only considered two strategies per species. An asymmetric number of strategies can induce further asymmetries in the interaction [46]. The intricacies of multiplayer games lend themselves to study such systems, but they also show that mutualistic interactions may be far more complex than often envisioned. Applying multiplayer game theory to mutualism unravels this dynamics between species and can be used to understand the complexity of these non-linear systems.

## Acknowledgments

We thank Noémie Erin, Christian Hilbe, Aniek Ivens, Martin Kalbe, Michael Lachmann, Eric Miller and Istvàn Scheuring for helpful discussions and suggestions. We also thank the referees and the editor for their detailed and constructive input. Financial support from the Emmy-Noether program of the DFG and from the Max Planck Society is gratefully acknowledged.

- Received July 20, 2012.
- Accepted August 22, 2012.

- This journal is © 2012 The Royal Society