Competing associations in bacterial warfare with two toxins

Influence and influenceability: global directionality in directed complex networks

Knowing which nodes are influential in a complex network and whether the network can be influenced by a small subset of nodes is a key part of network analysis. However, many traditional measures of importance focus on node level information without considering the global network architecture. We use the method of trophic analysis to study directed networks and show that both 'influence' and 'influenceability' in directed networks depend on the hierarchical structure and the global directionality, as measured by the trophic levels and trophic coherence, respectively. We show that in directed networks trophic hierarchy can explain: the nodes that can reach the most others; where the eigenvector centrality localizes; which nodes shape the behaviour in opinion or oscillator dynamics; and which strategies will be successful in generalized rock-paper-scissors games. We show, moreover, that these phenomena are mediated by the global directionality. We also highlight other structural properties of real networks related to influenceability, such as the pseudospectra, which depend on trophic coherence. These results apply to any directed network and the principles highlighted-that node hierarchy is essential for understanding network influence, mediated by global directionality-are applicable to many real-world dynamics.

Landscapes of Biochemical Warfare: Spatial Self-Organization Woven from Allelopathic Interactions

Evidence shows that diversity and spatial distributions of biological communities are largely driven by the race of living organisms in their adaptation to chemicals synthesized by their neighbors. In this report, the emergence of mathematical models on pure spatial self-organization induced by biochemical suppression (allelopathy) and competition between species were investigated through numerical analysis. For both random and patched initial spatial distributions of species, we demonstrate that warfare survivors are self-organized on the landscape in Turing-like patterns driven by diffusive instabilities of allelochemicals. These patterns are simple; either all species coexist at low diffusion rates or are massively extinct, except for a few at high diffusivities, but they are complex and biodiversity-sustained at intermediate diffusion rates. "Defensive alliances" and ecotones seem to be basic mechanisms that sustain great biodiversity in our hybrid cellular automata model. Moreover, species coexistence and extinction exhibit multi-stationarity.

Emerging spatiotemporal patterns in cyclic predator-prey systems with habitats

Three-species cyclic predator-prey systems are known to establish spiral waves that allow species to coexist. In this study, we analyze a structured heterogeneous system which gives one species an advantage to escape predation in an area that we refer to as a habitat and study the effect on species coexistence and emerging spatiotemporal patterns. Counterintuitively, the predator of the advantaged species emerges as dominant species with the highest average density inside the habitat. The species given the advantage in the form of an escape rate has the lowest average density until some threshold value for the escape rate is exceeded, after which the density of the species with the advantage overtakes that of its prey. Numerical analysis of the spatial density of each species as well as of the spatial two-point correlation function for both inside and outside the habitats allow a detailed quantitative discussion. Our analysis is extended to a six-species game that exhibits spontaneous spiral waves, which displays similar but more complicated results.

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.

Dynamically generated hierarchies in games of competition

Spatial many-species predator-prey systems have been shown to yield very rich space-time patterns. This observation begs the question whether there exist universal mechanisms for generating this type of emerging complex patterns in nonequilibrium systems. In this work we investigate the possibility of dynamically generated hierarchies in predator-prey systems. We analyze a nine-species model with competing interactions and show that the studied situation results in the spontaneous formation of spirals within spirals. The parameter dependence of these intriguing nested spirals is elucidated. This is achieved through the numerical investigation of various quantities (correlation lengths, densities of empty sites, Fourier analysis of species densities, interface fluctuations) that allows us to gain a rather complete understanding of the spatial arrangements and the temporal evolution of the system. A possible generalization of the interaction scheme yielding dynamically generated hierarchies is discussed. As cyclic interactions occur spontaneously in systems with competing strategies, the mechanism discussed in this work should contribute to our understanding of various social and biological systems.

Bacterial antagonism in host-associated microbial communities

Antagonistic interactions are abundant in microbial communities and contribute not only to the composition and relative proportions of their members but also to the longer-term stability of a community. This Review will largely focus on bacterial antagonism mediated by ribosomally synthesized peptides and proteins produced by members of host-associated microbial communities. We discuss recent findings on their diversity, functions, and ecological impacts. These systems play key roles in ecosystem defense, pathogen invasion, spatial segregation, and diversity but also confer indirect gains to the aggressor from products released by killed cells. Investigations into antagonistic bacterial interactions are important for our understanding of how the microbiota establish within hosts, influence health and disease, and offer insights into potential translational applications.

Invasion waves in the biochemical warfare between living organisms

Microorganisms and plants very commonly release toxic secondary chemical compounds (allelochemicals) that inhibit or kill sensitive strains or individuals from their own or other species. In this work we study a model that describes two species interacting through allelopathic suppression and competing for resources. Employing linear stability analysis, the conditions for coexistence or extinction of species in spatially homogeneous systems were determined. We found that the borders between the regimes of bistability, coexistence, and the extinction of the weaker by the stronger competitor, are altered by allelopathic interactions. In addition, traveling wave solutions for one species invasion were obtained considering the spatially explicit nature of the model. Our findings indicate that the minimum speed of the invasion wavefronts depends primarily on the competition coefficients and the parameters characterizing the species' functional responses to their allelochemicals. As a general rule, the species provided with the most effective chemical weapons dominates the population dynamics. Finally, we found a tristability at the coexistence region due to the combination of allelopathy and patchy population distributions in space. So, our model provides a distinct mechanism, independent of social behaviors, that produces such unexpected tristability impossible in classical competition models involving one-to-one individual interactions.

Coarsening with nontrivial in-domain dynamics: Correlations and interface fluctuations

Using numerical simulations we investigate the space-time properties of a system in which spirals emerge within coarsening domains, thus giving rise to nontrivial internal dynamics. Initially proposed in the context of population dynamics, the studied six-species model exhibits growing domains composed of three species in a rock-paper-scissors relationship. Through the investigation of different quantities, such as space-time correlations and the derived characteristic length, autocorrelation, density of empty sites, and interface width, we demonstrate that the nontrivial dynamics inside the domains affects the coarsening process as well as the properties of the interfaces separating different domains. Domain growth, aging, and interface fluctuations are shown to be governed by exponents whose values differ from those expected in systems with curvature driven coarsening.

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.

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.

Specialization and bet hedging in heterogeneous populations

Phenotypic heterogeneity is a strategy commonly used by bacteria to rapidly adapt to changing environmental conditions. Here, we study the interplay between phenotypic heterogeneity and genetic diversity in spatially extended populations. By analyzing the spatiotemporal dynamics, we show that the level of mobility and the type of competition qualitatively influence the persistence of phenotypic heterogeneity. While direct competition generally promotes persistence of phenotypic heterogeneity, specialization dominates in models with indirect competition irrespective of the degree of mobility.

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.

Diverging fluctuations in a spatial five-species cyclic dominance game

A five-species predator-prey model is studied on a square lattice where each species has two prey and two predators on the analogy to the rock-paper-scissors-lizard-Spock game. The evolution of the spatial distribution of species is governed by site exchange and invasion between the neighboring predator-prey pairs, where the cyclic symmetry can be characterized by two different invasion rates. The mean-field analysis has indicated periodic oscillations in the species densities with a frequency becoming zero for a specific ratio of invasion rates. When varying the ratio of invasion rates, the appearance of this zero-eigenvalue mode is accompanied by neutrality between the species associations. Monte Carlo simulations of the spatial system reveal diverging fluctuations at a specific invasion rate, which can be related to the vanishing dominance between all pairs of species associations.

Killing as means of promoting biodiversity

Bacteriocins are usually viewed as the effective weapons of bacterial killers. However, killing competitors with bacteriocins may be not only a means of eliminating other strains, but also a crucial unappreciated mechanism promoting bacterial diversity. In the present short review, we summarize recent empirical and theoretical studies examining the role bacteriocins that may play in driving and maintaining diversity among microbes. We conclude by highlighting limitations of current models and suggest directions for future studies.

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.

Defense mechanisms of empathetic players in the spatial ultimatum game

Experiments on the ultimatum game have revealed that humans are remarkably fond of fair play. When asked to share an amount of money, unfair offers are rare and their acceptance rate small. While empathy and spatiality may lead to the evolution of fairness, thus far considered continuous strategies have precluded the observation of solutions that would be driven by pattern formation. Here we introduce a spatial ultimatum game with discrete strategies, and we show that this simple alteration opens the gate to fascinatingly rich dynamical behavior. In addition to mixed stationary states, we report the occurrence of traveling waves and cyclic dominance, where one strategy in the cycle can be an alliance of two strategies. The highly webbed phase diagram, entailing continuous and discontinuous phase transitions, reveals hidden complexity in the pursuit of human fair play.

Conditional strategies and the evolution of cooperation in spatial public goods games

The fact that individuals will most likely behave differently in different situations begets the introduction of conditional strategies. Inspired by this, we study the evolution of cooperation in the spatial public goods game, where, besides unconditional cooperators and defectors, also different types of conditional cooperators compete for space. Conditional cooperators will contribute to the public good only if other players within the group are likely to cooperate as well but will withhold their contribution otherwise. Depending on the number of other cooperators that are required to elicit cooperation of a conditional cooperator, the latter can be classified in as many types as there are players within each group. We find that the most cautious cooperators, who require all other players within a group to be conditional cooperators, are the undisputed victors of the evolutionary process, even at very low synergy factors. We show that the remarkable promotion of cooperation is due primarily to the spontaneous emergence of quarantining of defectors, who become surrounded by conditional cooperators and are forced into isolated convex "bubbles" from which they are unable to exploit the public good. This phenomenon can be observed only in structured populations, thus adding to the relevance of pattern formation for the successful evolution of cooperation.

The evolution of bacteriocin production in bacterial biofilms

Bacteriocin production is a spiteful behavior of bacteria that is central to the competitive dynamics of many human pathogens. Social evolution predicts that bacteriocin production is favored when bacteriocin-producing cells are mixed at intermediate frequency with their competitors and when competitive neighborhoods are localized. Both predictions are supported by biofilm experiments. However, the means by which physical and biological processes interact to produce conditions that favor the evolution of bacteriocin production remain to be investigated. Here we fill this gap using analytical and computational approaches. We identify and collapse key parameters into a single number, the critical bacteriocin range, that measures the threshold distance from a focal bacteriocin-producing cell within which its fitness is higher than that of a sensitive cell. We develop an agent-based model to test our predictions and confirm that bacteriocin production is most favored when relatedness is intermediate and competition is local. We then use invasion analysis to determine evolutionarily stable strategies for bacteriocin production. Finally, we perform long-term evolutionary simulations to analyze how the critical bacteriocin range and genetic lineage segregation affect biodiversity in multistrain biofilms. We find that biodiversity is maintained in highly segregated biofilms for a wide array of critical bacteriocin ranges. However, under conditions of high nutrient penetration leading to well-mixed biofilms, biodiversity rapidly decreases and becomes sensitive to the critical bacteriocin range.

Interaction rules affect species coexistence in intransitive networks

Intransitive communities, those in which species' abilities cannot be ranked in a hierarchy, have been the focus of theoretical and empirical research, as intransitivity could help explain the maintenance of biodiversity. Here we show that models for intransitive competition embedding slightly different interaction rules can produce opposite patterns. In particular, we find that interactions in which an individual can be outcompeted by its neighbors, but cannot outcompete its neighbors, produce negative frequency dependence that, in turn, promotes coexistence. Whenever the interaction rule is modified toward symmetry (the individual and the neighbors can outcompete each other) the negative frequency dependence vanishes, producing different coexistence levels. Macroscopically, we find that asymmetric interactions yield highest biodiversity if species compete globally, while symmetric interactions favor highest biodiversity if competition takes place locally.

Phase diagrams for the spatial public goods game with pool punishment

The efficiency of institutionalized punishment is studied by evaluating the stationary states in the spatial public goods game comprising unconditional defectors, cooperators, and cooperating pool punishers as the three competing strategies. Fines and costs of pool punishment are considered as the two main parameters determining the stationary distributions of strategies on the square lattice. Each player collects a payoff from five five-person public goods games, and the evolution of strategies is subsequently governed by imitation based on pairwise comparisons at a low level of noise. The impact of pool punishment on the evolution of cooperation in structured populations is significantly different from that reported previously for peer punishment. Representative phase diagrams reveal remarkably rich behavior, depending also on the value of the synergy factor that characterizes the efficiency of investments payed into the common pool. Besides traditional single- and two-strategy stationary states, a rock-paper-scissors type of cyclic dominance can emerge in strikingly different ways.

Self-organizing patterns maintained by competing associations in a six-species predator-prey model

Formation and competition of associations are studied in a six-species ecological model where each species has two predators and two prey. Each site of a square lattice is occupied by an individual belonging to one of the six species. The evolution of the spatial distribution of species is governed by iterated invasions between the neighboring predator-prey pairs with species specific rates and by site exchange between the neutral pairs with a probability X . This dynamical rule yields the formation of five associations composed of two or three species with proper spatiotemporal patterns. For large X a cyclic dominance can occur between the three two-species associations whereas one of the two three-species associations prevails in the whole system for low values of X in the final state. Within an intermediate range of X all the five associations coexist due to the fact that cyclic invasions between the two-species associations reduce their resistance temporarily against the invasion of three-species associations.