Spatial rock-paper-scissors models with inhomogeneous reaction rates

Computing macroscopic reaction rates in reaction-diffusion systems using Monte Carlo simulations

Stochastic reaction-diffusion models are employed to represent many complex physical, biological, societal, and ecological systems. The macroscopic reaction rates describing the large-scale, long-time kinetics in such systems are effective, scale-dependent renormalized parameters that need to be either measured experimentally or computed by means of a microscopic model. In a Monte Carlo simulation of stochastic reaction-diffusion systems, microscopic probabilities for specific events to happen serve as the input control parameters. To match the results of any computer simulation to observations or experiments carried out on the macroscale, a mapping is required between the microscopic probabilities that define the Monte Carlo algorithm and the macroscopic reaction rates that are experimentally measured. Finding the functional dependence of emergent macroscopic rates on the microscopic probabilities (subject to specific rules of interaction) is a very difficult problem, and there is currently no systematic, accurate analytical way to achieve this goal. Therefore, we introduce a straightforward numerical method of using lattice Monte Carlo simulations to evaluate the macroscopic reaction rates by directly obtaining the count statistics of how many events occur per simulation time step. Our technique is first tested on well-understood fundamental examples, namely, restricted birth processes, diffusion-limited two-particle coagulation, and two-species pair annihilation kinetics. Next we utilize the thus gained experience to investigate how the microscopic algorithmic probabilities become coarse-grained into effective macroscopic rates in more complex model systems such as the Lotka-Volterra model for predator-prey competition and coexistence, as well as the rock-paper-scissors or cyclic Lotka-Volterra model and its May-Leonard variant that capture population dynamics with cyclic dominance motifs. Thereby we achieve a more thorough and deeper understanding of coarse graining in spatially extended stochastic reaction-diffusion systems and the nontrivial relationships between the associated microscopic and macroscopic model parameters, with a focus on ecological systems. The proposed technique should generally provide a useful means to better fit Monte Carlo simulation results to experimental or observational data.

Response of a three-species cyclic ecosystem to a short-lived elevation of death rate

A balanced ecosystem with coexisting constituent species is often perturbed by different natural events that persist only for a finite duration of time. What becomes important is whether, in the aftermath, the ecosystem recovers its balance or not. Here we study the fate of an ecosystem by monitoring the dynamics of a particular species that encounters a sudden increase in death rate. For exploration of the fate of the species, we use Monte-Carlo simulation on a three-species cyclic rock-paper-scissor model. The density of the affected (by perturbation) species is found to drop exponentially immediately after the pulse is applied. In spite of showing this exponential decay as a short-time behavior, there exists a region in parameter space where this species surprisingly remains as a single survivor, wiping out the other two which had not been directly affected by the perturbation. Numerical simulations using stochastic differential equations of the species give consistency to our results.

A solution of the Crow-Kimura evolution model on fluctuating fitness landscape

The article discusses the Crow-Kimura model in the context of random transitions between different fitness landscapes. The duration of epochs, during which the fitness landscape is constant over time, is modeled by an exponential distribution. To obtain an exact solution, a system of functional equations is required. However, to approximate the model, we consider the cases of slow or fast transitions and calculate the first-order corrections using either the transition rate or its inverse. Specifically, we focus on the case of slow transitions and find that the average fitness is equal to the average fitness for evolution on static fitness landscapes, but with the addition of a load term. We also investigate the model for a small number of genes and identify the exact transition points to the transient phase.

Eco-evolutionary cyclic dominance among predators, prey, and parasites

Predator-prey interactions are one of ecology's central research themes, but with many interdisciplinary implications across the social and natural sciences. Here we consider an often-overlooked species in these interactions, namely parasites. We first show that a simple predator-prey-parasite model, inspired by the classical Lotka-Volterra equations, fails to produce a stable coexistence of all three species, thus failing to provide a biologically realistic outcome. To improve this, we introduce free space as a relevant eco-evolutionary component in a new mathematical model that uses a game-theoretical payoff matrix to describe a more realistic setup. We then show that the consideration of free space stabilizes the dynamics by means of cyclic dominance that emerges between the three species. We determine the parameter regions of coexistence as well as the types of bifurcations leading to it by means of analytical derivations as well as by means of numerical simulations. We conclude that the consideration of free space as a finite resource reveals the limits of biodiversity in predator-prey-parasite interactions, and it may also help us in the determination of factors that promote a healthy biota.

Universal scaling of extinction time in stochastic evolutionary dynamics

Evolutionary dynamics is well captured by the replicator equations when the population is infinite and well-mixed. However, the extinction dynamics is modified with finite and structured populations. Experiments on the non-transitive ecosystem containing three populations of bacteria found that the ecological stability sensitively depends on the spatial structure of the populations. Based on the Reference-Gamble-Birth algorithm, we use agent-based Monte Carlo simulations to investigate the extinction dynamics in the rock-paper-scissors ecosystem with finite and structured populations. On the fully-connected network, the extinction time in stable and unstable regimes falls into two universal functions when plotted with the rescaled variables. On the two dimensional grid, the spatial structure changes the transition boundary between stable and unstable regimes but doesn't change its extinction trend. The finding of universal scaling in extinction dynamics is unexpected, and may provide a powerful method to classify different evolutionary dynamics into universal classes.

Community formation in wealth-mediated thermodynamic strategy evolution

We study a dynamical system defined by a repeated game on a 1D lattice, in which the players keep track of their gross payoffs over time in a bank. Strategy updates are governed by a Boltzmann distribution, which depends on the neighborhood bank values associated with each strategy, relative to a temperature scale, which defines the random fluctuations. Players with higher bank values are, thus, less likely to change strategy than players with a lower bank value. For a parameterized rock-paper-scissors game, we derive a condition under which communities of a given strategy form with either fixed or drifting boundaries. We show the effect of a temperature increase on the underlying system and identify surprising properties of this model through numerical simulations.

Survival of the weakest in non-transitive asymmetric interactions among strains of E. coli

Hierarchical organization in ecology, whereby interactions are nested in a manner that leads to a dominant species, naturally result in the exclusion of all but the dominant competitor. Alternatively, non-hierarchical competitive dynamics, such as cyclical interactions, can sustain biodiversity. Here, we designed a simple microbial community with three strains of E. coli that cyclically interact through (i) the inhibition of protein production, (ii) the digestion of genomic DNA, and (iii) the disruption of the cell membrane. We find that intrinsic differences in these three major mechanisms of bacterial warfare lead to an unbalanced community that is dominated by the weakest strain. We also use a computational model to describe how the relative toxin strengths, initial fractional occupancies, and spatial patterns affect the maintenance of biodiversity. The engineering of active warfare between microbial species establishes a framework for exploration of the underlying principles that drive complex ecological interactions.

The maintenance of ecological diversity depends on the strength and direction of competitive interactions, but these interactions are difficult to study in microbial communities. Here the authors use engineered E. coli strains to show that competitively weak strains can persist when pairwise interactions are asymmetrical.

Mortality makes coexistence vulnerable in evolutionary game of rock-paper-scissors

Multiple species in the ecosystem are believed to compete cyclically for maintaining balance in nature. The evolutionary dynamics of cyclic interaction crucially depends on different interactions representing different natural habits. Based on a rock-paper-scissors model of cyclic competition, we explore the role of mortality of individual organisms in the collective survival of a species. For this purpose a parameter called "natural death" is introduced. It is meant for bringing about the decease of an individual irrespective of any intra- and interspecific interaction. We perform a Monte Carlo simulation followed by a stability analysis of different fixed points of defined rate equations and observe that the natural death rate is surprisingly one of the most significant factors in deciding whether an ecosystem would come up with a coexistence or a single-species survival.

Population Dynamics in a Changing Environment: Random versus Periodic Switching

Environmental changes greatly influence the evolution of populations. Here, we study the dynamics of a population of two strains, one growing slightly faster than the other, competing for resources in a time-varying binary environment modeled by a carrying capacity switching either randomly or periodically between states of abundance and scarcity. The population dynamics is characterized by demographic noise (birth and death events) coupled to a varying environment. We elucidate the similarities and differences of the evolution subject to a stochastically and periodically varying environment. Importantly, the population size distribution is generally found to be broader under intermediate and fast random switching than under periodic variations, which results in markedly different asymptotic behaviors between the fixation probability of random and periodic switching. We also determine the detailed conditions under which the fixation probability of the slow strain is maximal.

The effect of habitats and fitness on species coexistence in systems with cyclic dominance

Cyclic dominance between species may yield spiral waves that are known to provide a mechanism enabling persistent species coexistence. This observation holds true even in presence of spatial heterogeneity in the form of quenched disorder. In this work we study the effects on spatio-temporal patterns and species coexistence of structured spatial heterogeneity in the form of habitats that locally provide one of the species with an advantage. Performing extensive numerical simulations of systems with three and six species we show that these structured habitats destabilize spiral waves. Analyzing extinction events, we find that species extinction probabilities display a succession of maxima as function of time, that indicate a periodically enhanced probability for species extinction. Analysis of the mean extinction time reveals that as a function of the parameter governing the advantage of one of the species a transition between stable coexistence and unstable coexistence takes place. We also investigate how efficiency as a predator or a prey affects species coexistence.

Emergence of oscillatory coexistence with exponentially decayed waiting times in a coupled cyclic competition system

Interpatch migration between two environments is generally considered as a spatial concept and can affect species biodiversity in each patch by inducing flux of population such as inflow and outflow quantities of species. In this paper, we explore the effect of interpatch migration, which can be generally considered as a spatial concept and may affect species biodiversity between two different patches in the perspective of the macroscopic level by exploiting the coupling of two systems, where each patch is occupied by cyclically competing three species who can stably coexist by exhibiting periodic orbits. For two simple scenarios of interpatch migration either single or all species migration, we found that two systems with independently stable coexisting species in each patch are eventually synchronized, and oscillatory behaviors of species densities in two patches become identical, i.e., the synchronized coexistence emerges. In addition, we find that, whether single or all species interpatch migration occurs, the waiting time for the synchronization is exponentially decreasing as the coupling strength is intensified. Our findings suggest that the synchronized behavior of species as a result of migration between different patches can be easily predicted by the coupling of systems and additional information such as waiting times and sensitivity of initial densities.

Robust coexistence with alternative competition strategy in the spatial cyclic game of five species

Alternative strategy is common in animal populations to promote reproductive fitness by obtaining resources. In spatial dynamics of cyclic competition, reproduction can occur when individuals obtain vacant rooms and, in this regard, empty sites should be resources for reproduction which can be induced by interspecific competition. In this paper, we study the role of alternative competition in the spatial system of cyclically competing five species by utilizing rock-paper-scissors-lizard-spock game. From Monte-Carlo simulations, we found that strong alternative competition can lead to the reemergence of coexistence of five species regardless of mobility, which is never reported in previous works under the symmetric competition structure. By investigating the coexistence probability, we also found that coexistence alternates by passing certain degrees of alternative competition in combination with mobility. In addition, we provided evidences in the opposite scenario by strengthening spontaneous competition, which exhibits the reemergence of coexistence similarly. Our findings may suggest more comprehensive perspectives to interpret mechanisms for biodiversity by alternative strategies in spatially extended systems than previously reported.

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.

Nonlinear dynamics with Hopf bifurcations by targeted mutation in the system of rock-paper-scissors metaphor

The role of mutation, which is an error process in gene evolution, in systems of cyclically competing species has been studied from various perspectives, and it is regarded as one of the key factors for promoting coexistence of all species. In addition to naturally occurring mutations, many experiments in genetic engineering have involved targeted mutation techniques such as recombination between DNA and somatic cell sequences and have studied genetic modifications through loss or augmentation of cell functions. In this paper, we investigate nonlinear dynamics with targeted mutation in cyclically competing species. In different ways to classic approaches of mutation in cyclic games, we assume that mutation may occur in targeted individuals who have been removed from intraspecific competition. By investigating each scenario depending on the number of objects for targeted mutation analytically and numerically, we found that targeted mutation can lead to persistent coexistence of all species. In addition, under the specific condition of targeted mutation, we found that targeted mutation can lead to emergences of bistable states for species survival. Through the linear stability analysis of rate equations, we found that those phenomena are accompanied by Hopf bifurcation which is supercritical. Our findings may provide more global perspectives on understanding underlying mechanisms to control biodiversity in ecological/biological sciences, and evidences with mathematical foundations to resolve social dilemmas such as a turnover of group members by resigning with intragroup conflicts in social sciences.

Three-species competition with non-deterministic outcomes

Theoretical and experimental research studies have shown that ecosystems governed by non-transitive competition networks tend to maintain high levels of biodiversity. The theoretical body of work, however, has mainly focused on competition networks in which the outcomes of competition events are predetermined and hence deterministic, and where all species are identical up to their competitive relationships, an assumption that may limit the applicability of theoretical results to real-life situations. In this paper, we aim to probe the robustness of the link between biodiversity and non-transitive competition by introducing a three-dimensional winning probability parameter space, making the outcomes of competition events in a three-species in silico ecosystem uncertain. While two degenerate points in this parameter space have been the subject of previous studies, we investigate the remaining settings, which equip the species with distinct competitive abilities. We find that the impact of this modification depends on the spatial dimension of the system. When the system is well mixed, it collapses to monoculture, as is also the case in the non-transitive deterministic setting. In one dimension, chaotic patterns emerge, which tend to maintain biodiversity, and a power law relates the time that species manage to coexist to the degree of uncertainty regarding competition event outcomes. In two dimensions, the formation of spiral wave patterns ensures that biodiversity is maintained for moderate degrees of uncertainty, while considerable deviations from the non-transitive deterministic setting have strong negative effects on species coexistence. It can hence be concluded that non-transitive competition can still produce coexistence when the assumption of deterministic competition is abandoned. When the system collapses to monoculture, one observes a "survival of the strongest" law, as the species that has the highest probability of defeating its competitors has the best odds to become the sole survivor.

Eco-evolutionary dynamics of a population with randomly switching carrying capacity

Environmental variability greatly influences the eco-evolutionary dynamics of a population, i.e. it affects how its size and composition evolve. Here, we study a well-mixed population of finite and fluctuating size whose growth is limited by a randomly switching carrying capacity. This models the environmental fluctuations between states of resources abundance and scarcity. The population consists of two strains, one growing slightly faster than the other, competing under two scenarios: one in which competition is solely for resources, and one in which the slow (cooperating) strain produces a public good (PG) that benefits also the fast (free-riding) strain. We investigate how the coupling of demographic and environmental (external) noise affects the population's eco-evolutionary dynamics. By analytical and computational means, we study the correlations between the population size and its composition, and discuss the social-dilemma-like 'eco-evolutionary game' characterizing the PG production. We determine in what conditions it is best to produce a PG; when cooperating is beneficial but outcompeted by free riding, and when the PG production is detrimental for cooperators. Within a linear noise approximation to populations of varying size, we also accurately analyse the coupled effects of demographic and environmental noise on the size distribution.

How directional mobility affects coexistence in rock-paper-scissors models

This work deals with a system of three distinct species that changes in time under the presence of mobility, selection, and reproduction, as in the popular rock-paper-scissors game. The novelty of the current study is the modification of the mobility rule to the case of directional mobility, in which the species move following the direction associated to a larger (averaged) number density of selection targets in the surrounding neighborhood. Directional mobility can be used to simulate eyes that see or a nose that smells, and we show how it may contribute to reduce the probability of coexistence.

Survival behavior in the cyclic Lotka-Volterra model with a randomly switching reaction rate

We study the influence of a randomly switching reproduction-predation rate on the survival behavior of the nonspatial cyclic Lotka-Volterra model, also known as the zero-sum rock-paper-scissors game, used to metaphorically describe the cyclic competition between three species. In large and finite populations, demographic fluctuations (internal noise) drive two species to extinction in a finite time, while the species with the smallest reproduction-predation rate is the most likely to be the surviving one (law of the weakest). Here we model environmental (external) noise by assuming that the reproduction-predation rate of the strongest species (the fastest to reproduce and predate) in a given static environment randomly switches between two values corresponding to more and less favorable external conditions. We study the joint effect of environmental and demographic noise on the species survival probabilities and on the mean extinction time. In particular, we investigate whether the survival probabilities follow the law of the weakest and analyze their dependence on the external noise intensity and switching rate. Remarkably, when, on average, there is a finite number of switches prior to extinction, the survival probability of the predator of the species whose reaction rate switches typically varies nonmonotonically with the external noise intensity (with optimal survival about a critical noise strength). We also outline the relationship with the case where all reaction rates switch on markedly different time scales.

Evolution of a Fluctuating Population in a Randomly Switching Environment

Environment plays a fundamental role in the competition for resources, and hence in the evolution of populations. Here, we study a well-mixed, finite population consisting of two strains competing for the limited resources provided by an environment that randomly switches between states of abundance and scarcity. Assuming that one strain grows slightly faster than the other, we consider two scenarios-one of pure resource competition, and one in which one strain provides a public good-and investigate how environmental randomness (external noise) coupled to demographic (internal) noise determines the population's fixation properties and size distribution. By analytical means and simulations, we show that these coupled sources of noise can significantly enhance the fixation probability of the slower-growing species. We also show that the population size distribution can be unimodal, bimodal, or multimodal and undergoes noise-induced transitions between these regimes when the rate of switching matches the population's growth rate.

Basins of distinct asymptotic states in the cyclically competing mobile five species game

We study the dynamics of cyclic competing mobile five species on spatially extended systems originated from asymmetric initial populations and investigate the basins for the three possible asymptotic states, coexistence of all species, existences of only two independent species, and the extinction. Through extensive numerical simulations, we find a prosperous dependence on initial conditions for species biodiversity. In particular, for fixed given equal densities of two relevant species, we find that only five basins for the existence of two independent species exist and they are spirally entangled for high mobility. A basin of coexistence is outbreaking when the mobility parameter is decreased through a critical value and surrounded by the other five basins. For fixed given equal densities of two independent species, however, we find that basin structures are not spirally entangled. Further, final states of two independent species are totally different. For all possible considerations, the extinction state is not witnessed which is verified by the survival probability. To provide the validity of basin structures from lattice simulations, we analyze the system in mean-field manners. Consequently, results on macroscopic levels are matched to direct lattice simulations for high mobility regimes. These findings provide a good insight into the fundamental issue of the biodiversity among many species than previous cases.

Strategic tradeoffs in competitor dynamics on adaptive networks

Recent empirical work highlights the heterogeneity of social competitions such as political campaigns: proponents of some ideologies seek debate and conversation, others create echo chambers. While symmetric and static network structure is typically used as a substrate to study such competitor dynamics, network structure can instead be interpreted as a signature of the competitor strategies, yielding competition dynamics on adaptive networks. Here we demonstrate that tradeoffs between aggressiveness and defensiveness (i.e., targeting adversaries vs. targeting like-minded individuals) creates paradoxical behaviour such as non-transitive dynamics. And while there is an optimal strategy in a two competitor system, three competitor systems have no such solution; the introduction of extreme strategies can easily affect the outcome of a competition, even if the extreme strategies have no chance of winning. Not only are these results reminiscent of classic paradoxical results from evolutionary game theory, but the structure of social networks created by our model can be mapped to particular forms of payoff matrices. Consequently, social structure can act as a measurable metric for social games which in turn allows us to provide a game theoretical perspective on online political debates.

Emergence of unusual coexistence states in cyclic game systems

Evolutionary games of cyclic competitions have been extensively studied to gain insights into one of the most fundamental phenomena in nature: biodiversity that seems to be excluded by the principle of natural selection. The Rock-Paper-Scissors (RPS) game of three species and its extensions [e.g., the Rock-Paper-Scissors-Lizard-Spock (RPSLS) game] are paradigmatic models in this field. In all previous studies, the intrinsic symmetry associated with cyclic competitions imposes a limitation on the resulting coexistence states, leading to only selective types of such states. We investigate the effect of nonuniform intraspecific competitions on coexistence and find that a wider spectrum of coexistence states can emerge and persist. This surprising finding is substantiated using three classes of cyclic game models through stability analysis, Monte Carlo simulations and continuous spatiotemporal dynamical evolution from partial differential equations. Our finding indicates that intraspecific competitions or alternative symmetry-breaking mechanisms can promote biodiversity to a broader extent than previously thought.

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.

Zealotry promotes coexistence in the rock-paper-scissors model of cyclic dominance

Cyclic dominance models, such as the classic rock-paper-scissors (RPS) game, have found real-world applications in biology, ecology, and sociology. A key quantity of interest in such models is the coexistence time, i.e., the time until at least one population type goes extinct. Much recent research has considered conditions that lengthen coexistence times in an RPS model. A general finding is that coexistence is promoted by localized spatial interactions (low mobility), while extinction is fostered by global interactions (high mobility). That is, there exists a mobility threshold which separates a regime of long coexistence from a regime of rapid collapse of coexistence. The key finding of our paper is that if zealots (i.e., nodes able to defeat others while themselves being immune to defeat) of even a single type exist, then system coexistence time can be significantly prolonged, even in the presence of global interactions. This work thus highlights a crucial determinant of system survival time in cyclic dominance models.

Evolutionary game theory: cells as players

In two papers we review game theory applications in biology below the level of cognitive living beings. It can be seen that evolution and natural selection replace the rationality of the actors appropriately. Even in these micro worlds, competing situations and cooperative relationships can be found and modeled by evolutionary game theory. Also those units of the lowest levels of life show different strategies for different environmental situations or different partners. We give a wide overview of evolutionary game theory applications to microscopic units. In this first review situations on the cellular level are tackled. In particular metabolic problems are discussed, such as ATP-producing pathways, secretion of public goods and cross-feeding. Further topics are cyclic competition among more than two partners, intra- and inter-cellular signalling, the struggle between pathogens and the immune system, and the interactions of cancer cells. Moreover, we introduce the theoretical basics to encourage scientists to investigate problems in cell biology and molecular biology by evolutionary game theory.

Synchronization and extinction in cyclic games with mixed strategies

We consider cyclic Lotka-Volterra models with three and four strategies where at every interaction agents play a strategy using a time-dependent probability distribution. Agents learn from a loss by reducing the probability to play a losing strategy at the next interaction. For that, an agent is described as an urn containing β balls of three and four types, respectively, where after a loss one of the balls corresponding to the losing strategy is replaced by a ball representing the winning strategy. Using both mean-field rate equations and numerical simulations, we investigate a range of quantities that allows us to characterize the properties of these cyclic models with time-dependent probability distributions. For the three-strategy case in a spatial setting we observe a transition from neutrally stable to stable when changing the level of discretization of the probability distribution. For large values of β, yielding a good approximation to a continuous distribution, spatially synchronized temporal oscillations dominate the system. For the four-strategy game the system is always neutrally stable, but different regimes emerge, depending on the size of the system and the level of discretization.

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.

Mesoscopic interactions and species coexistence in evolutionary game dynamics of cyclic competitions

Evolutionary dynamical models for cyclic competitions of three species (e.g., rock, paper, and scissors, or RPS) provide a paradigm, at the microscopic level of individual interactions, to address many issues in coexistence and biodiversity. Real ecosystems often involve competitions among more than three species. By extending the RPS game model to five (rock-paper-scissors-lizard-Spock, or RPSLS) mobile species, we uncover a fundamental type of mesoscopic interactions among subgroups of species. In particular, competitions at the microscopic level lead to the emergence of various local groups in different regions of the space, each involving three species. It is the interactions among the groups that fundamentally determine how many species can coexist. In fact, as the mobility is increased from zero, two transitions can occur: one from a five- to a three-species coexistence state and another from the latter to a uniform, single-species state. We develop a mean-field theory to show that, in order to understand the first transition, group interactions at the mesoscopic scale must be taken into account. Our findings suggest, more broadly, the importance of mesoscopic interactions in coexistence of great many species.

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.

Self-organizing patterns in an evolutionary rock-paper-scissors game for stochastic synchronized strategy updates

We study a spatial evolutionary rock-paper-scissors game with synchronized strategy updating. Players gain their payoff from games with their four neighbors on a square lattice and can update their strategies simultaneously according to the logit rule, which is the noisy version of the best-response dynamics. For the synchronized strategy update two types of global oscillations (with an ordered strategy arrangement and periods of three and six generations) can occur in this system in the zero noise limit. At low noise values, all nine oscillating phases are present in the system by forming a self-organizing spatial pattern due to the comprising invasion and speciation processes along the interfaces separating the different domains.

Characterization of spiraling patterns in spatial rock-paper-scissors games

The spatiotemporal arrangement of interacting populations often influences the maintenance of species diversity and is a subject of intense research. Here, we study the spatiotemporal patterns arising from the cyclic competition between three species in two dimensions. Inspired by recent experiments, we consider a generic metapopulation model comprising "rock-paper-scissors" interactions via dominance removal and replacement, reproduction, mutations, pair exchange, and hopping of individuals. By combining analytical and numerical methods, we obtain the model's phase diagram near its Hopf bifurcation and quantitatively characterize the properties of the spiraling patterns arising in each phase. The phases characterizing the cyclic competition away from the Hopf bifurcation (at low mutation rate) are also investigated. Our analytical approach relies on the careful analysis of the properties of the complex Ginzburg-Landau equation derived through a controlled (perturbative) multiscale expansion around the model's Hopf bifurcation. Our results allow us to clarify when spatial "rock-paper-scissors" competition leads to stable spiral waves and under which circumstances they are influenced by nonlinear mobility.

Nutrient-responsive regulation determines biodiversity in a colicin-mediated bacterial community

BACKGROUND

Antagonistic interactions mediated by antibiotics are strong drivers of bacterial community dynamics which shape biodiversity. Colicin production by Escherichia coli is such an interaction that governs intraspecific competition and is involved in promoting biodiversity. It is unknown how environmental cues affect regulation of the colicin operon and thus influence antibiotic-mediated community dynamics.

RESULTS

Here, we investigate the community dynamics of colicin-producing, -sensitive, and -resistant/non-producer E. coli strains that colonize a microfabricated spatially-structured habitat. Nutrients are found to strongly influence community dynamics: when growing on amino acids and peptides, colicin-mediated competition is intense and the three strains do not coexist unless spatially separated at large scales (millimeters). Surprisingly, when growing on sugars, colicin-mediated competition is minimal and the three strains coexist at the micrometer scale. Carbon storage regulator A (CsrA) is found to play a key role in translating the type of nutrients into the observed community dynamics by controlling colicin release. We demonstrate that by mitigating lysis, CsrA shapes the community dynamics and determines whether the three strains coexist. Indeed, a mutant producer that is unable to suppress colicin release, causes the collapse of biodiversity in media that would otherwise support co-localized growth of the three strains.

CONCLUSIONS

Our results show how the environmental regulation of an antagonistic trait shapes community dynamics. We demonstrate that nutrient-responsive regulation of colicin release by CsrA, determines whether colicin producer, resistant non-producer, and sensitive strains coexist at small spatial scales, or whether the sensitive strain is eradicated. This study highlights how molecular-level regulatory mechanisms that govern interference competition give rise to community-level biodiversity patterns.

From pairwise to group interactions in games of cyclic dominance

We study the rock-paper-scissors game in structured populations, where the invasion rates determine individual payoffs that govern the process of strategy change. The traditional version of the game is recovered if the payoffs for each potential invasion stem from a single pairwise interaction. However, the transformation of invasion rates to payoffs also allows the usage of larger interaction ranges. In addition to the traditional pairwise interaction, we therefore consider simultaneous interactions with all nearest neighbors, as well as with all nearest and next-nearest neighbors, thus effectively going from single pair to group interactions in games of cyclic dominance. We show that differences in the interaction range affect not only the stationary fractions of strategies but also their relations of dominance. The transition from pairwise to group interactions can thus decelerate and even revert the direction of the invasion between the competing strategies. Like in evolutionary social dilemmas, in games of cyclic dominance, too, the indirect multipoint interactions that are due to group interactions hence play a pivotal role. Our results indicate that, in addition to the invasion rates, the interaction range is at least as important for the maintenance of biodiversity among cyclically competing strategies.

Mobility-dependent selection of competing strategy associations

Standard models of population dynamics focus on the interaction, survival, and extinction of the competing species individually. Real ecological systems, however, are characterized by an abundance of species (or strategies, in the terminology of evolutionary-game theory) that form intricate, complex interaction networks. The description of the ensuing dynamics may be aided by studying associations of certain strategies rather than individual ones. Here we show how such a higher-level description can bear fruitful insight. Motivated from different strains of colicinogenic Escherichia coli bacteria, we investigate a four-strategy system which contains a three-strategy cycle and a neutral alliance of two strategies. We find that the stochastic, spatial model exhibits a mobility-dependent selection of either the three-strategy cycle or of the neutral pair. We analyze this intriguing phenomenon numerically and analytically.

Global attractors and extinction dynamics of cyclically competing species

Transitions to absorbing states are of fundamental importance in nonequilibrium physics as well as ecology. In ecology, absorbing states correspond to the extinction of species. We here study the spatial population dynamics of three cyclically interacting species. The interaction scheme comprises both direct competition between species as in the cyclic Lotka-Volterra model, and separated selection and reproduction processes as in the May-Leonard model. We show that the dynamic processes leading to the transient maintenance of biodiversity are closely linked to attractors of the nonlinear dynamics for the overall species' concentrations. The characteristics of these global attractors change qualitatively at certain threshold values of the mobility and depend on the relative strength of the different types of competition between species. They give information about the scaling of extinction times with the system size and thereby the stability of biodiversity. We define an effective free energy as the negative logarithm of the probability to find the system in a specific global state before reaching one of the absorbing states. The global attractors then correspond to minima of this effective energy landscape and determine the most probable values for the species' global concentrations. As in equilibrium thermodynamics, qualitative changes in the effective free energy landscape indicate and characterize the underlying nonequilibrium phase transitions. We provide the complete phase diagrams for the population dynamics and give a comprehensive analysis of the spatio-temporal dynamics and routes to extinction in the respective phases.

Labyrinthine clustering in a spatial rock-paper-scissors ecosystem

The spatial rock-paper-scissors ecosystem, where three species interact cyclically, is a model example of how spatial structure can maintain biodiversity. We here consider such a system for a broad range of interaction rates. When one species grows very slowly, this species and its prey dominate the system by self-organizing into a labyrinthine configuration in which the third species propagates. The cluster size distributions of the two dominating species have heavy tails and the configuration is stabilized through a complex spatial feedback loop. We introduce a statistical measure that quantifies the amount of clustering in the spatial system by comparison with its mean-field approximation. Hereby, we are able to quantitatively explain how the labyrinthine configuration slows down the dynamics and stabilizes the system.

Discriminating the effects of spatial extent and population size in cyclic competition among species

We introduce a population model for species under cyclic competition. This model allows individuals to coexist and interact on single cells while migration takes place between adjacent cells. In contrast to the model introduced by Reichenbach, Mobilia, and Frey [Reichenbach, Mobilia, and Frey, Nature (London) 448, 1046 (2007)], we find that the emergence of spirals results in an ambiguous behavior regarding the stability of coexistence. The typical time until extinction exhibits, however, a qualitatively opposite dependence on the newly introduced nonunit carrying capacity in the spiraling and the nonspiraling regimes. This allows us to determine a critical mobility that marks the onset of this spiraling state sharply. In contrast, we demonstrate that the conventional finite size stability analysis with respect to spatial size is of limited use for identifying the onset of the spiraling regime.

Extinction in neutrally stable stochastic Lotka-Volterra models

Populations of competing biological species exhibit a fascinating interplay between the nonlinear dynamics of evolutionary selection forces and random fluctuations arising from the stochastic nature of the interactions. The processes leading to extinction of species, whose understanding is a key component in the study of evolution and biodiversity, are influenced by both of these factors. Here, we investigate a class of stochastic population dynamics models based on generalized Lotka-Volterra systems. In the case of neutral stability of the underlying deterministic model, the impact of intrinsic noise on the survival of species is dramatic: It destroys coexistence of interacting species on a time scale proportional to the population size. We introduce a new method based on stochastic averaging which allows one to understand this extinction process quantitatively by reduction to a lower-dimensional effective dynamics. This is performed analytically for two highly symmetrical models and can be generalized numerically to more complex situations. The extinction probability distributions and other quantities of interest we obtain show excellent agreement with simulations.

Saddles, arrows, and spirals: deterministic trajectories in cyclic competition of four species

Population dynamics in systems composed of cyclically competing species has been of increasing interest recently. Here we investigate a system with four or more species. Using mean field theory, we study in detail the trajectories in configuration space of the population fractions. We discover a variety of orbits, shaped like saddles, spirals, and straight lines. Many of their properties are found explicitly. Most remarkably, we identify a collective variable that evolves simply as an exponential: Q ∝ e(λt), where λ is a function of the reaction rates. It provides information on the state of the system for late times (as well as for t→-∞). We discuss implications of these results for the evolution of a finite, stochastic system. A generalization to an arbitrary number of cyclically competing species yields valuable insights into universal properties of such systems.