Abstract

Ecological processes, such as reproduction, mobility, and interaction between species, play important roles in the maintenance of biodiversity. Classically, the cyclic dominance of species has been modelled using the nonhierarchical interactions among competing species, represented by the “Rock-Paper-Scissors” (RPS) game. Here we propose a cascaded channel model for analyzing the existence of biodiversity in the RPS game. The transition between successive generations is modelled as communication of information over a noisy communication channel. The rate of transfer of information over successive generations is studied using mutual information and it is found that “greedy” information transfer between successive generations may lead to conditions for extinction. This generalized framework can be used to study biodiversity in any number of interacting species, ecosystems with unequal rates for different species, and also competitive networks.

1. Introduction

Ecological processes, such as reproduction, mobility, and interaction between species, have been shown to play important roles in the maintenance of biodiversity [15]. One of the useful methods for investigating codevelopment of populations is to use evolutionary game theory [68]. In this context, the Rock-Paper-Scissors (RPS) game is a paradigm to describe species diversity [3, 4]. The cyclic dominance of species is modelled using the nonhierarchical interactions among competing species, represented by the “Rock-Paper-Scissors” (RPS) game in which rock crushes scissors, scissors cut paper, and paper wraps rock [916]. Communities of subpopulations or species exhibiting such dynamics have been identified in several ecosystems [17]. For example, cyclic dominance of three male strategies has been reported in lizards [18].

It has been shown that cyclic dominance alone is not sufficient to preserve biodiversity. When the interactions between individuals are local, spatially separated domains are dominated by one species form, leading to stable coexistence [13]. It is shown that beyond a critical value of mobility biodiversity is lost [10]. Researchers have shown that spiral formation takes place without a conservation law for the total density and, in general, fast diffusion can destroy species coexistence [11]. The final survival probabilities characterizing different species in cyclic competition have also been reported [12]. The interplay of evolutionary dynamics, discreteness of the population, and the nature of the interactions also influence the coexistence of strategies [13]. Pattern formations due to cyclic dominance of three species, operating near a bifurcation point, have also been investigated [14]. The spreading of epidemic and its effect on biodiversity has been studied using the RPS game [15]. The role of swarming in the context of biodiversity has also been investigated [16]. A recently published review paper on the cyclic dominance in evolutionary games covers in detail the current advances on RPS and related evolutionary games. In particular, it focuses on pattern formation, the impact of mobility, and the spontaneous emergence of cyclic dominance. The authors also review mean-field and zero-dimensional RPS models and the application of the complex Ginzburg-Landau equation and highlight the importance of statistical physics for the study of large-scale ecological systems [19]. Theoretical explorations of the relevance of spatial structure in evolutionary games have also received a lot of attention in the recent past [2022]. Realistic models, beyond well-mixed populations, have been studied recently [23, 24]. Researchers have also used mean-field theory to show that group interactions at the mesoscopic scale have been taken into account for the coexistence of a large number of species [25]. Recent work on mobility and velocity-enhanced cooperation of moving agents has further linked mobility to biodiversity [26, 27]. The Ginzburg-Landau theory has also been applied successfully to study biodiversity using models of spatial RPS games [28, 29]. Investigation of alliance formation among multiple number of interacting species shows that stable solutions are possible, depending on the parameter values that determine invasion and mixing in the multispecies predator-prey systems [30].

The primary motivation of this work is to model the information transfer between successive generations of evolution and link it to biodiversity. Further, we wish to explore the effect of the rate of information transfer on the existence of biodiversity. Yet another motivation is to study the nonhierarchical cyclic interactions between species using this information-theoretic approach. In this paper, we provide a novel approach using a cascaded communication channel for analyzing biodiversity in the classic Lotka-Volterra model. Specifically, we investigate the transitions from one generation to the next based on the transition probability matrix between two consecutive generations. The transition between successive generations is modelled as communication of information through a noisy communication channel. We derive explicit relation between the probabilities of the species in two consecutive generations in terms of the three basic processes: selection, reproduction, and mobility. This gives us an elegant method to study the evolution of the species for long waiting times and ultimately leads to the condition for the existence of biodiversity. We also show how selection, reproduction, and mobility together influence the existence of biodiversity.

The paper is organized as follows. Section 1 is the Introduction. In Section 2, the system model is discussed in detail. In Section 3, the transition from one generation to the next is modelled as a problem of communicating symbols across a noisy channel. The information transfer across generations is linked to the existence of biodiversity in Section 4. The proposed method is generalized for species in Section 5. Results from lattice simulations are presented in Section 6. Finally, the paper is concluded in Section 7.

2. System Model

We consider the cyclic Lotka-Volterra model where three states , , and cyclically dominate each other [31]. Let the time densities for the three species be , , and , respectively. The deterministic rate equations describing the time evolution of these three species can be written as

As the spatial version of the Lotka-Volterra model, we consider three subpopulations or species (, , and ), the members of which are distributed randomly over an spatial lattice. There are also empty spaces () on the lattice, which are essentially unoccupied. For the possible interactions with the neighbours, we consider a version of the RPS game, namely, a stochastic spatial variant of the model [4]. The Rock, Paper, and Scissors are depicted by , , and , respectively. Figure 1 shows the basic processes of selection, reproduction, and mobility, modelled as the Poisson processes, occurring at rates , , and , respectively. The individuals interact only with their immediate neighbors.

3. Modeling Transition between Generations

Let us model the transition from generation to as a problem of communicating four types of symbols (, and ) across a noisy channel. This channel occasionally makes errors which result in one symbol being received as a different symbol. This is similar to a standard communication channel problem encountered in electrical engineering [32]. This channel model is shown in Figure 2. In our case, the symbols are the lattice locations occupied by the individuals of the three species (or the empty locations). In the case of transition from one generation to the next, a site in the lattice currently occupied by a particular species may become occupied by an individual from another species due to mobility or fall vacant due to selection. Alternatively, an empty site may become occupied by a member from one of the species due to reproduction. It is also possible that there is no change in the occupancy of a lattice site. These transition probabilities are influenced by the selection rate, reproduction rate, and mobility.

The transition probabilities are based on averaged probabilities (densities) ignoring the spatial structure. This is a mean-field approach, assuming a well-mixed population with the limit of population size . The idea can be extended by considering joint densities of pairs of neighbours, and, to approximate higher order densities, one can consider products of lower order densities (so-called pair approximation and moment closure) [33].

The channel transition probability matrix, for a well-mixed population, can be written asAny element of the channel transition probability matrix represents the probability of transition of a particular site occupied by a member of species to the same site occupied by a member of species in the next generation. Here, , and 4 correspond to , , , and , respectively. For example, represents the probability that the site occupied by species will become occupied by species in the subsequent generation. This is possible only by mobility, and hence . Similarly, represents the probability that the site occupied by species will translate to in the next generation. This is possible either by selection () or by mobility (), which are mutually exclusive events, and hence the probabilities add up. It is important to note that the entries in the transition probability matrix () are probabilities consisting of normalized such that each row adds up to unity. In order to differentiate, the rates have been denoted by , , and .

From (2), the eigenvalues [34] of matrix are found out to be , , and . In the case of biodiversity, ; that is, in the process of cyclic dominance, the occupation of the sites in the lattice by the three species is equiprobable. Let be the probability of occurrence of species , , or at a particular lattice site for generation and let be the probability of occurrence of species , , or at a particular lattice site for generation . We can alternately write the transition probability matrix as , where is the matrix of the eigenvectors of and is the diagonal matrix with the diagonal elements as the eigenvalues [34]. Thus, can be written, in terms of the eigenvalues, aswhere We note that the transition probability matrices given in (2) and (3) are equivalent. Using (3) we can express the probabilities of individual species. For example, for species , the next generation may result from , , , or , as dictated by the transition probability matrix given in (3). Thus, for species , the probability of the next generation can be expressed asSimilar expressions can be derived for species and . Solving (5), we obtainEquation (6) gives and explicit relation between the probabilities of species in generation to the probabilities of species, , in the consecutive generation .

So far we have considered two consecutive generations only. We now develop a mathematical framework for large waiting times, that is, for a large number of generations (). As we let the generations evolve using the stochastic spatial variant of the of the RPS game, it is equivalent to communicating through a cascade of noisy channels, as depicted in Figure 3(a). The equivalent channel transition probability matrix for cascaded channels [32] is given by . Since we can express , the eigenvalues of are given simply by , where are the eigenvalues of . Thus, cascaded channels can be represented as one “equivalent channel,” as shown in Figure 3(b). From (6), we can easily write the probabilities of species , , and , after generations, as

Since the matrix is a probability transition matrix, [22]. Thus, . Hence, from (7) we can conclude that a well-mixed population must always result in biodiversity.

We now consider the real-life scenario when the population is not well-mixed. In fact, clustering is routinely observed as the population evolves [10, 14, 31]. In this case, also the assumption of large population sizes exists; that is, size . Let us represent the resultant transition matrix for a non-well-mixed (clustered) scenario aswhere represents the transition matrix corresponding to a well-mixed population and represents the effect of local clustering at generation . We observe that the effect of the matrix is to alter the transition probabilities (the branch labels of Figure 2). Since the resultant matrix, , is also a probability transition matrix with the individual rows adding up to unity, must be a matrix with individual rows adding up to zero. After generations, the effective channel can be represented by where is the resultant transition matrix for generation . The relation between the probabilities for species , , and and the vacant locations, , can be expressed aswhere represents the resulting probabilities for species , , and and the vacant locations, , after generations, while represents the starting probabilities. Extinction implies that one of the species dies out; that is, the left hand side of (9) should be or or for arbitrary starting probabilities . This is possible if one of the first three columns of is zero, which gives us the necessary condition for extinction (uniformity), . Thus, the condition for the existence of biodiversity can be stated asSince , the condition for biodiversity can also be written as for all . This is one of our central results. An important contribution of this paper is to model the “transitions” of cells like the “transitions” of symbols over a noisy channel in communications. This leads to a matrix representation, which eventually has an elegant solution in terms of eigenvalues, when a cascade of channels is considered.

5. Generalization to Species

The given theory can be easily extended to study the nonhierarchical cyclic interactions between species in general. In this case, the two consecutive generations of interacting species can be linked using an transition probability matrix with eigenvalues given by , , and . For a well-mixed population, the expression for the probabilities of species, after generations, is calculated to beWe surmise from (11) that, in case of biodiversity, the lattice area occupied by each species, on an average, will be where is the total number of lattice sites.

Another interesting aspect of cascaded channels is study of mutual information between consecutive generations. A discrete channel is a system with input alphabet , output alphabet , and a probability transition matrix . The amount of information conveyed by the discrete channel is quantified by the mutual information between and and is defined as [32]where is the entropy of and is the entropy of given . The entropy, (in bits), is calculated usingand, in general, , where is a discrete random variable. The notion of mutual information is routinely used to determine the capacity of a channel in the area of communications. Here we wish to study how the rate of transfer of information over successive generations affects the possibility of coexistence among species. In our case, for a    matrix, the expression for mutual information can be written asSolving for mutual information, we obtainwhere , and  are given in (4).

6. Lattice Simulations

Extensive computer simulations were carried out to test the condition for biodiversity, as predicted by (10). Cyclic dominance is modeled using a stochastic lattice [2, 4]. In our stochastic lattice simulations, we have arranged the three species on a two-dimensional square lattice with periodic boundary conditions. By periodic boundary conditions, we imply that the sites on the left edge of the grid are viewed as adjacent to ones on the right edge and those on the bottom are viewed as adjacent to those on the top. Every lattice site is occupied by an individual of species , species , or species or left empty (). At each simulation step, a random individual interacts with one of its eight nearest neighbours (corresponding to the move of the King in a chess game). The choice of the neighbour is also randomly determined (each of the eight neighbours are equally likely). The interaction with the neighbour, reproduction or migration, as well as the corresponding waiting time, is carried out according to a built-in function in MATLAB that generates random values from a Poisson distribution. In our simulations, one generation is counted when every individual has reacted, on an average, once.

Figure 4(a) shows a typical plot of the eigenvalues of : (red), (green), and (blue) when biodiversity occurs. Here the lattice size is 100 × 100. It is clear from the plot that , and hence biodiversity is observed. In Figure 4(b), we show the case when extinction happens. Here, as soon as , one of the species dies out.

Simulations were also carried out for different values of , , and , and the system was tested for the existence of biodiversity (or uniformity) after long waiting periods. These operating points are plotted in Figure 5(a). We show the results for different system sizes: 50 × 50 lattice sites, 70 × 70 lattice sites, and 100 × 100 lattice sites. To obtain the extinction probability, the results are averaged over 100 realizations. Simulation results corroborate the theory. The probability of extinction, , tends to zero for the cases where , while leads to uniformity. We observe that as the system size increases, the transition between biodiversity and uniformity sharpens. In order to determine , the actual values of the normalized rates , and were extracted by carrying out a frequency count of the number of lattice-site transitions from one generation to the next and then averaging these values over 500 generations. Figure 5(b) shows a typical snapshot from lattice simulations after a long waiting time for the case when biodiversity exists.

The proposed mathematical framework can be conveniently used to study systems with unequal reaction rates by appropriately modifying the transition matrix, . Suppose the three species have different normalized rates for selection, reproduction, and mobility, given by , , and 3. Then the corresponding transition probability matrix, for unequal normalized rates, can be written as

The eigenvalues of this generalized transition probability matrix can be used to study the effects of the different parameters on biodiversity. As an illustrative example, suppose we wish to investigate the effect of unequal (normalized) reproduction rates , , and for the three species ( and are the same for all three species). From (12) we can calculate the eigenvalues as , , and .

Another interesting observation from nature is the fact that resource competitors can benefit one another through containment of shared competitors [35, 36]. When several factors determine the outcome of competition, the interacting species can be modeled as a competitive network, as shown in Figure 6(a). Here, species are shown as nodes, and arrows connect the competitive inferior to the superior competitor. Researchers have shown that intransitivity in competitive networks can maintain diversity [3, 37]. The proposed cascaded channel model can be easily adapted to include different numbers of competitive relationships in a competitive network, as illustrated in Figure 6(a). Figure 6(b) gives the generic channel model for species and competitive relationships which can be represented as a transition probability matrix and analyzed for the existence of biodiversity as discussed earlier.

Figure 7 shows the plot of mutual information versus generation number for typical cases where coexistence is observed (cases (a) and (b)) and typical cases where extinction occurs (cases (c) and (d)). These simulations are for lattice size 100 × 100. In these experiments, for different sets of values for , we determine whether the lattice simulations lead to biodiversity or uniformity. It is interesting to observe that between successive generations is fairly stable for the cases when biodiversity occurs, while it shows an increasing trend for the cases when uniformity occurs. The trend is depicted by the dotted line in the figure. It appears that “greedy” information transfer between successive generations may lead to conditions for extinction.

7. Conclusion

We have developed a framework using cascaded channel model for analyzing the existence of biodiversity in nonhierarchical interactions among competing species. This mathematical framework gives a clear insight into how the three parameters, reproduction rate, selection rate, and mobility, influence biodiversity. To give a biological interpretation, we can draw a parallel to colicinogenic strains of E. coli growing on a Petri dish [3]. These strains can produce a toxin and the “antidote” which determine the selection rate (). However, both the toxin and the “antidote” are generated by using resources and hence at the cost of their reproduction rate (). The mobility of E. coli can be increased by using supersoft agar. On the other hand, this mobility can be lowered by increasing the agar concentration. This is modeled by our mobility parameter .

The significant contributions of this paper are as follows.(1)The transition between successive generations is modelled as communication of information through a noisy communication channel. This information theoretic model is a paradigm shift and is not linked to any of the previously reported models. For example, using this model one can understand the interdependence on reproduction rate, selection rate, and mobility together on the existence of biodiversity.(2)The cascaded channel model exploits a very unique characteristic of matrices: the ability to express square matrices in terms of their eigenvalues. This approach permits a cascade of channels to be modelled as a product of several matrices. This approach leads to an elegant solution in terms of the power of the eigenvalues.(3)The versatility of this approach allows it to be extended to study more than three interacting species, species with ecosystems with unequal rates for different species and competitive networks with different competitive relationships.(4)The rate of transfer of information over successive generations is studied using mutual information, and it is found that “greedy” information transfer between successive generations may lead to conditions for extinction.(5)Explicit dependence of the probability of population distribution on the transition probabilities linked to the reproduction rate, selection rate, and mobility is derived in this paper. The condition for biodiversity is derived, which is corroborated by simulations.(6)This generalized mathematical framework can be used to study biodiversity in any number of interacting species, ecosystems with unequal rates for different species, and also competitive networks.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.