Globally synchronized oscillations in complex cyclic games

Nutrients and flow shape the cyclic dominance games between Escherichia coli strains

Evolutionary game theory has provided various models to explain the coexistence of competing strategies, one of which is the rock-paper-scissors (RPS) game. A system of three Escherichia coli strains-a toxin-producer, a resistant and a sensitive-has become a classic experimental model for studying RPS games. Previous experimental and theoretical studies, however, often ignored the influence of ecological factors such as nutrients and toxin dynamics on the evolutionary game dynamics. In this work, we combine experiments and modelling to study how these factors affect competition dynamics. Using three-dimensional printed mini-bioreactors, we tracked the frequency of the three strains in different culturing media and under different flow regimes. Although our experimental system fulfilled the requirements of cyclic dominance, we did not observe clear cycles or long-term coexistence between strains. We found that both nutrients and flow rates strongly impacted population dynamics. In our simulations, we explicitly modelled the release, removal and diffusion of toxin. We showed that the amount of toxin that is retained in the system is a simple indicator that can predict competition outcomes across broad parameter space. Moreover, our simulation results suggest that high rates of toxin diffusion might have prevented cyclic patterns from emerging in our experimental system. This article is part of the theme issue 'Half a century of evolutionary games: a synthesis of theory, application and future directions'.

Lotka-Volterra versus May-Leonard formulations of the spatial stochastic rock-paper-scissors model: The missing link

The rock-paper-scissors (RPS) model successfully reproduces some of the main features of simple cyclic predator-prey systems with interspecific competition observed in nature. Still, lattice-based simulations of the spatial stochastic RPS model are known to give rise to significantly different results, depending on whether the three-state Lotka-Volterra or the four-state May-Leonard formulation is employed. This is true independently of the values of the model parameters and of the use of either a von Neumann or a Moore neighborhood. In this paper, we introduce a simple modification to the standard spatial stochastic RPS model in which the range of the search of the nearest neighbor may be extended up to a maximum Euclidean radius R. We show that, with this adjustment, the Lotka-Volterra and May-Leonard formulations can be designed to produce similar results, both in terms of dynamical properties and spatial features, by means of an appropriate parameter choice. In particular, we show that this modified spatial stochastic RPS model naturally leads to the emergence of spiral patterns in both its three- and four-state formulations.

Environment driven oscillation in an off-lattice May-Leonard model

Cyclic dominance of competing species is an intensively used working hypothesis to explain biodiversity in certain living systems, where the evolutionary selection principle would dictate a single victor otherwise. Technically the May–Leonard models offer a mathematical framework to describe the mentioned non-transitive interaction of competing species when individual movement is also considered in a spatial system. Emerging rotating spirals composed by the competing species are frequently observed character of the resulting patterns. But how do these spiraling patterns change when we vary the external environment which affects the general vitality of individuals? Motivated by this question we suggest an off-lattice version of the tradition May–Leonard model which allows us to change the actual state of the environment gradually. This can be done by introducing a local carrying capacity parameter which value can be varied gently in an off-lattice environment. Our results support a previous analysis obtained in a more intricate metapopulation model and we show that the well-known rotating spirals become evident in a benign environment when the general density of the population is high. The accompanying time-dependent oscillation of competing species can also be detected where the amplitude and the frequency show a scaling law of the parameter that characterizes the state of the environment. These observations highlight that the assumed non-transitive interaction alone is insufficient condition to maintain biodiversity safely, but the actual state of the environment, which characterizes the general living conditions, also plays a decisive role on the evolution of related systems.

Exit rights open complex pathways to cooperation

We study the evolutionary dynamics of the Prisoner's Dilemma game in which cooperators and defectors interact with another actor type called exiters. Rather than being exploited by defectors, exiters exit the game in favour of a small pay-off. We find that this simple extension of the game allows cooperation to flourish in well-mixed populations when iterations or reputation are added. In networked populations, however, the exit option is less conducive to cooperation. Instead, it enables the coexistence of cooperators, defectors, and exiters through cyclic dominance. Other outcomes are also possible as the exit pay-off increases or the network structure changes, including network-wide oscillations in actor abundances that may cause the extinction of exiters and the domination of defectors, although game parameters should favour exiting. The complex dynamics that emerges in the wake of a simple option to exit the game implies that nuances matter even if our analyses are restricted to incentives for rational behaviour.

Performance of weak species in the simplest generalization of the rock-paper-scissors model to four species

We investigate the problem of the predominance and survival of "weak" species in the context of the simplest generalization of the spatial stochastic rock-paper-scissors model to four species by considering models in which one, two, or three species have a reduced predation probability. We show, using lattice based spatial stochastic simulations with random initial conditions, that if only one of the four species has its probability reduced, then the most abundant species is the prey of the "weakest" (assuming that the simulations are large enough for coexistence to prevail). Also, among the remaining cases, we present examples in which "weak" and "strong" species have similar average abundances and others in which either of them dominates-the most abundant species being always a prey of a weak species with which it maintains a unidirectional predator-prey interaction. However, in contrast to the three-species model, we find no systematic difference in the global performance of weak and strong species, and we conjecture that a similar result will hold if the number of species is further increased. We also determine the probability of single species survival and coexistence as a function of the lattice size, discussing its dependence on initial conditions and on the change to the dynamics of the model which results from the extinction of one of the species.

A novel route to cyclic dominance in voluntary social dilemmas

Cooperation is the backbone of modern human societies, making it a priority to understand how successful cooperation-sustaining mechanisms operate. Cyclic dominance, a non-transitive set-up comprising at least three strategies wherein the first strategy overrules the second, which overrules the third, which, in turn, overrules the first strategy, is known to maintain biodiversity, drive competition between bacterial strains, and preserve cooperation in social dilemmas. Here, we present a novel route to cyclic dominance in voluntary social dilemmas by adding to the traditional mix of cooperators, defectors and loners, a fourth player type, risk-averse hedgers, who enact tit-for-tat upon paying a hedging cost to avoid being exploited. When this cost is sufficiently small, cooperators, defectors and hedgers enter a loop of cyclic dominance that preserves cooperation even under the most adverse conditions. By contrast, when the hedging cost is large, hedgers disappear, consequently reverting to the traditional interplay of cooperators, defectors, and loners. In the interim region of hedging costs, complex evolutionary dynamics ensues, prompting transitions between states with two, three or four competing strategies. Our results thus reveal that voluntary participation is but one pathway to sustained cooperation via cyclic dominance.

Predominance of the weakest species in Lotka-Volterra and May-Leonard formulations of the rock-paper-scissors model

We revisit the problem of the predominance of the "weakest" species in the context of Lotka-Volterra and May-Leonard formulations of a spatial stochastic rock-paper-scissors model in which one of the species has its predation probability reduced by 0<P_{w}<1. We show that, despite the different population dynamics and spatial patterns, these two formulations lead to qualitatively similar results for the late time values of the relative abundances of the three species (as a function of P_{w}), as long as the simulation lattices are sufficiently large for coexistence to prevail-the "weakest" species generally having an advantage over the others (specially over its predator). However, for smaller simulation lattices, we find that the relatively large oscillations at the initial stages of simulations with random initial conditions may result in a significant dependence of the probability of species survival on the lattice size.

Invasion-controlled pattern formation in a generalized multispecies predator-prey system

Rock-scissors-paper game, as the simplest model of intransitive relation between competing agents, is a frequently quoted model to explain the stable diversity of competitors in the race of surviving. When increasing the number of competitors we may face a novel situation because beside the mentioned unidirectional predator-prey-like dominance a balanced or peer relation can emerge between some competitors. By utilizing this possibility in the present work we generalize a four-state predator-prey-type model where we establish two groups of species labeled by even and odd numbers. In particular, we introduce different invasion probabilities between and within these groups, which results in a tunable intensity of bidirectional invasion among peer species. Our study reveals an exceptional richness of pattern formations where five quantitatively different phases are observed by varying solely the strength of the mentioned inner invasion. The related transition points can be identified with the help of appropriate order parameters based on the spatial autocorrelation decay, on the fraction of empty sites, and on the variance of the species density. Furthermore, the application of diverse, alliance-specific inner invasion rates for different groups may result in the extinction of the pair of species where this inner invasion is moderate. These observations highlight that beyond the well-known and intensively studied cyclic dominance there is an additional source of complexity of pattern formation that has not been explored earlier.

Biodiversity in models of cyclic dominance is preserved by heterogeneity in site-specific invasion rates

Global, population-wide oscillations in models of cyclic dominance may result in the collapse of biodiversity due to the accidental extinction of one species in the loop. Previous research has shown that such oscillations can emerge if the interaction network has small-world properties, and more generally, because of long-range interactions among individuals or because of mobility. But although these features are all common in nature, global oscillations are rarely observed in actual biological systems. This begets the question what is the missing ingredient that would prevent local oscillations to synchronize across the population to form global oscillations. Here we show that, although heterogeneous species-specific invasion rates fail to have a noticeable impact on species coexistence, randomness in site-specific invasion rates successfully hinders the emergence of global oscillations and thus preserves biodiversity. Our model takes into account that the environment is often not uniform but rather spatially heterogeneous, which may influence the success of microscopic dynamics locally. This prevents the synchronization of locally emerging oscillations, and ultimately results in a phenomenon where one type of randomness is used to mitigate the adverse effects of other types of randomness in the system.

Zealots tame oscillations in the spatial rock-paper-scissors game

The rock-paper-scissors game is a paradigmatic model for biodiversity, with applications ranging from microbial populations to human societies. Research has shown, however, that mobility jeopardizes biodiversity by promoting the formation of spiral waves, especially if there is no conservation law in place for the total number of competing players. First, we show that even if such a conservation law applies, mobility still jeopardizes biodiversity in the spatial rock-paper-scissors game if only a small fraction of links of the square lattice is randomly rewired. Secondly, we show that zealots are very effective in taming the amplitude of oscillations that emerge due to mobility and/or interaction randomness, and this regardless of whether the later is quenched or annealed. While even a tiny fraction of zealots brings significant benefits, at 5% occupancy zealots practically destroy all oscillations regardless of the intensity of mobility, and regardless of the type and strength of randomness in the interaction structure. Interestingly, by annealed randomness the impact of zealots is qualitatively the same as by mobility, which highlights that fast diffusion does not necessarily destroy the coexistence of species, and that zealotry thus helps to recover the stable mean-field solution. Our results strengthen the important role of zealots in models of cyclic dominance, and they reveal fascinating evolutionary outcomes in structured populations that are a unique consequence of such uncompromising behavior.

How turbulence regulates biodiversity in systems with cyclic competition

Cyclic, nonhierarchical interactions among biological species represent a general mechanism by which ecosystems are able to maintain high levels of biodiversity. However, species coexistence is often possible only in spatially extended systems with a limited range of dispersal, whereas in well-mixed environments models for cyclic competition often lead to a loss of biodiversity. Here we consider the dispersal of biological species in a fluid environment, where mixing is achieved by a combination of advection and diffusion. In particular, we perform a detailed numerical analysis of a model composed of turbulent advection, diffusive transport, and cyclic interactions among biological species in two spatial dimensions and discuss the circumstances under which biodiversity is maintained when external environmental conditions, such as resource supply, are uniform in space. Cyclic interactions are represented by a model with three competitors, resembling the children's game of rock-paper-scissors, whereas the flow field is obtained from a direct numerical simulation of two-dimensional turbulence with hyperviscosity. It is shown that the space-averaged dynamics undergoes bifurcations as the relative strengths of advection and diffusion compared to biological interactions are varied.

Cyclic dominance in evolutionary games: a review

Rock is wrapped by paper, paper is cut by scissors and scissors are crushed by rock. This simple game is popular among children and adults to decide on trivial disputes that have no obvious winner, but cyclic dominance is also at the heart of predator-prey interactions, the mating strategy of side-blotched lizards, the overgrowth of marine sessile organisms and competition in microbial populations. Cyclical interactions also emerge spontaneously in evolutionary games entailing volunteering, reward, punishment, and in fact are common when the competing strategies are three or more, regardless of the particularities of the game. Here, we review recent advances on the rock-paper-scissors (RPS) and related evolutionary games, focusing, in particular, on pattern formation, the impact of mobility and the spontaneous emergence of cyclic dominance. We also review mean-field and zero-dimensional RPS models and the application of the complex Ginzburg-Landau equation, and we highlight the importance and usefulness of statistical physics for the successful study of large-scale ecological systems. Directions for future research, related, for example, to dynamical effects of coevolutionary rules and invasion reversals owing to multi-point interactions, are also outlined.