Abstract
Background
Antagonistic species interactions can lead to coevolutionary genotype or phenotype frequency oscillations, with important implications for ecological and evolutionary processes. However, direct empirical evidence of such oscillations is rare. The rarity of observations is generally attributed to inherent difficulties of ecological and evolutionary longterm studies, to weak or absent interaction between species, or to the absence of negative frequencydependence.
Results
Here, we show that another factor – nongenetic inheritance, mediated for example by epigenetic mechanisms – can completely eliminate oscillations in the presence of such negative frequency dependence, even if only a small fraction of offspring are affected. We analytically derive the threshold value of this fraction at which the dynamics change from oscillatory to stable, and investigate how selection, mutation and generation times differences between the two species affect the threshold value. These results strongly suggest that the lack of phenotype frequency oscillations should not be attributed to the lack of strong interactions between antagonistic species.
Conclusions
Given increasing evidence of nongenetic effects on the outcomes of antagonistic species interactions, we suggest that these effects should be incorporated into ecological and evolutionary models of interacting species.
Keywords:
Host; Parasite; Coevolution; Epigenetic inheritance; Resistance; Red queen dynamicsBackground
The phenotypes of species are generally assumed to be adapted to their environment by natural selection. A change in an environment can therefore lead to an evolutionary change in phenotypes as species adapt to new circumstances. Environments comprise both biotic and abiotic elements, and evolutionary change in one species is often driven by evolutionary change in another species. Indeed, ecology is dominated by species interactions such as predation, parasitism, mutualism and competition. If species interactions are antagonistic (i.e., one species benefits at the expense of another), the resulting patterns of adaptation and counteradaptation can lead to cyclical dynamics typical of predatorprey or hostparasite systems. Understanding the causes and consequences of such fluctuating population dynamics is crucial in a number of biological phenomena, and particularly also in applied fields such as conservation biology and pest management.
The population dynamics of antagonistic species interactions can be captured with wellestablished models such as the LotkaVolterra model [1], the NicholsonBailey model [2], or the Red Queen model [3]. The Red Queen model stands out as a coevolutionary model because it does not primarily focus on fluctuating population densities, but rather on fluctuating genotype and phenotype frequencies of the interacting species. The causes and consequences of fluctuating genotype and phenotype frequencies in hostparasite, hostparasitoid and predatorprey interactions [4] are increasingly well understood at least in twospecies systems, but direct empirical evidence of longterm dynamics is rare [5], not at least because longterm dynamics are inherently difficult to measure [6,7].
Phenotypic adaptations to changing environments need not be driven by natural selection alone. This is because many phenotypes are plastic and can change due to adverse environmental conditions, a property generally referred to as phenotypic plasticity. Interestingly, phenotypic change can be stably transmitted across generations at various levels of specificity. Transgenerational induction of defences has been reported in animals and plants [810]. The maternal transfer of antibodies in vertebrates is a well known phenomenon, and in recent years, it has become clear that both vertebrates and invertebrates exhibit transgenerational immunity (see [11,12] and references therein). In the context of hostparasite coevolution, one of the most remarkable demonstrations has been given by [13] who have provided evidence for strainspecific immunity transmitted from mother to offspring in Daphnia magna infected with the pathogenic bacteria Pasteuria ramosa. Also, transgenerational phenotypic inheritance of virulence traits has been demonstrated in the malaria parasite Plasmodium falciparum[14]. In the microbial world, phenotypic switching has been reported both as a direct response to environmental change [15] and as a stochastic event [16] anticipating environmental change, and phenotypic states are often inherited stably across generations [17,18]. For a recent review of nongenetic inheritance and its evolutionary implications see [19].
Our goal here is to understand the effect of nongenetic inheritance on patterns of antagonistic coevolution. We develop a simple model where two species (e.g., host and parasite) are interacting, and each species has two alternative phenotypes. If their phenotypes match, the outcome of the interaction has negative fitness consequences for one species (host) and positive for the other species (parasite). As a result of this, the phenotypes harmed by the interaction may switch to the alternative phenotype in the offspring. We are purposefully ignorant about the nature of the phenotype (e.g., molecular, developmental, behavioral) and about the underlying form of nongenetic inheritance responsible for the phenotype switch in the offspring. In the absence of nongenetic inheritance, this model reduces to the most basic model of antagonistic coevolution exhibiting negative frequency dependence and resulting in the classical Red Queen dynamics (i.e., oscillations of phenotypes). We find that nongenetic inheritance can strongly affect cycling behavior typical of Red Queen dynamics by dampening the phenotype frequency oscillations. To examine this in detail, we derive analytical expressions of the threshold rate at which this elimination occurs.
Methods
In order to understand how nongenetic inheritance affects the patterns of antagonistic coevolution, we consider a simple, discretegeneration, coevolutionary model of two species X and Y which interact antagonistically, e.g., hosts and parasites or predators and their preys. Each species is represented as a haploid, singlelocus genotype with two possible alleles. The locus can be a genetic factor (gene or genotype) encoding for a given phenotype, or simply a phenotype itself. The two populations, X (host or prey) and Y (parasite or predator) thus carry two alternative phenotypes, 1 and 2, and the model tracks the frequency of each phenotype in every generation. We assume both populations to be infinitely large and initiate their phenotype frequencies at random. To approach its longterm dynamics, coevolution of X and Y proceeds for 11000 generations, and only during the last 1000 generations are the measurements taken. At each generation, both species undergo selection and reproduce; the crucial feature of the latter process is the ability to switch phenotypes due to antagonistic interaction.
Antagonistic interactions induce fitness costs on both species: successful interactions come at a cost for species X while unsuccessful interactions come at a cost for species Y . We assume that successful interaction occurs between the corresponding phenotypes. We denote the frequencies of phenotypes 1 and 2 from species X as x_{1} and x_{2}, and phenotype frequencies from species Y as y_{1}and y_{2}. Only individuals from species Y with phenotype 1 successfully attack individuals from species X with phenotype 1, and only individuals from species Y with phenotype 2 successfully attack individuals from species X with phenotype 2. This results in selection against the matching phenotypes in X and the nonmatching phenotypes in Y , such that only a fraction 1−s_{X }and 1−s_{Y}, respectively, survive in the next generation (see Table 1); the frequencies in the two population are changed accordingly. Both species then undergo reproduction, which may involve induced phenotype switching and stochastic phenotype switching. Induced switching is induced by the antagonistic interaction: individuals which are harmed by the interaction transmit the opposite phenotype to the next generation in the following proportions: α_{X }for species X and α_{Y }for species Y . The phenotype frequencies after selection and induced switching read
where we have assumed that x_{1} + x_{2} = y_{1} + y_{2} = 1, x = x_{1}, and y = y_{1}. Equation (1) can be derived in three steps. First, one calculates the proportion of individuals of species X which are matched by species Y , here equal to xy. Second, this proportion of individuals undergoes selection, which means that the frequencies of the matched individuals are multiplied by , while the frequencies of the unaffected individuals are multiplied by , where is the species X mean population fitness. Third, the matched individuals undergo induced switching, meaning that the proportion α_{X} of them switches to an alternative phenotype. This yields the first equation in (1), and the calculation for species Y is analogous. As those individuals which undergo induced switching are also those which have undergone selection, the selection coefficient s_{X} becomes equivalent to a cost of induced switching, as seen in equation (1) above (see also Discussion).
Table 1. The fitness values resulting from the antagonistic interaction between X and Y
In contrast, stochastic switching occurs independently of antagonistic interactions, and in proportion μ in both species X and Y . Therefore, the frequencies after stochastic switching (and thus after one generation) are given by
This step can be also interpreted as mutation, and we generally assume that μ = 10^{−8} unless mentioned otherwise.
Finally, we allow for asymmetry in the generation time between the two species by defining a parameter g, which denotes the number of generations that species Y undergoes in a single generation of species X. During one generation of Y , a fraction 1/gof the population X is updated according to the equations given above, while the fraction 1−1/gremains unchanged. This process is then repeated g times, and the resulting frequencies x^{′′} and y^{′′ }yield the phenotype frequencies after an entire generation of species X[20]. By default, we assume g = 1 unless mentioned otherwise.
Results
It is generally expected that antagonistic interactions can result in cyclic allele frequency dynamics, reflecting a continuing arms race between the two species [21]. In the absence of induced phenotypic switching (α_{X }= α_{Y }= 0) our model reveals such a pattern (see Figure 1A). In this situation, if the common phenotype of species X (say phenotype 1) is more likely to interact antagonistically with the corresponding matching phenotype of species Y (phenotype 1), then another phenotype 2 of species X has a selective advantage causing a gradual increase of the frequency x_{2} and a decrease of the frequency x_{1}. Such change will in turn drive the frequency change in species Y by selecting for phenotype 2, causing y_{2} to increase, and so on. These oscillations are, in the absence of random genetic drift and mutation, expected to continue indefinitely, otherwise fixation or extinction of one of the two phenotypes occurs [22].
Figure 1. Impact of induced phenotypic switching on the cyclic phenotype frequency dynamics. (A): Coevolutionary dynamics between antagonistic phenotypes are predicted to continue indefinitely when no induced switching occurs due to timelagged, negative frequency dependent selection (α_{X }= 0). (B): When induced switching occurs (here in only in species X, thus α_{Y }= 0) at a low rate ( ), allele frequency cycles persist in time but at an altered amplitude and speed. (C): When the switching rate exceeds a threshold value , cycles begin to dampen and reach a stable equilibrium. (DE): Increased levels of induced switching also decrease the amplitude and increase the speed of the cycles. Note that in panels AC and the one in panels DF is different because of different selection coefficients used; the relation between the strength of selection and the persistence of cycles is examined in detail in the subsequent figure. The following parameter values were used: (AC) s_{X }= s_{Y }= 0.3, (DE) s_{X }= s_{Y }= 0.65; (A) α_{X }= 0, (B) α_{X }= 0.03, (C) α_{X }= 0.1; In all panels we used α_{Y }= 0. Period is defined as a number of generations during which the phenotype frequency cycles around to its original value.
Consider now a situation where induced phenotypic switching is possible in a single species. Figure 1BC shows the impact of such a process on the frequency dynamics between species X and Y . We see that as the switching rate increases in species X (α_{X }> 0, α_{Y }= 0), the cycles become faster and of lower amplitude, eventually leading to a stable state (x^{∗},y^{∗}) = (1/2,1/2); (Figure 1C). This happens when the switching rate α_{X} exceeds a certain threshold value, , such that when the cycles are maintained (even though with altered amplitude and frequency), and when , cycles dampen and reach a stable equilibrium. The changes in speed and amplitude of cycles are directly measured in Figure 1DE, and show that as α_{X }increases the amplitude gradually decreases to zero and the speed increases until the cycles disappear. This already illustrates that induced switching can fundamentally affect the oscillatory dynamics in the system.
In order to examine the persistence of cyclic dynamics in more detail, we derive an analytical expression for the stability of the cyclic behaviour as a function of α_{X}, α_{Y}, s_{X}, and s_{Y}.
In the case of the model considered here, the stability requires that x^{′′ }= x and y^{′′ }= y, where x^{′′ }≡ x(t + 1), x ≡ x(t), y^{′′ }≡ y(t + 1), and y ≡ y(t). It can be shown that the four trivial equilibria of this system are (x^{∗},y^{∗}) = (0,0), (x^{∗},y^{∗}) = (0,1), (x^{∗},y^{∗}) = (1,0), (x^{∗},y^{∗}) = (1,1), and that (x^{∗},y^{∗}) = (1/2,1/2) is a nontrivial equilibrium with the Jacobian of the form
The corresponding eigenvalues are
where
This has been derived under the assumptions of 0 ≤ s_{X},s_{Y }≤ 1, 0 ≤ α_{X},α_{Y }≤ 1, and 0 ≤ μ < 0.5. The condition for the stability at (x^{∗},y^{∗}) = (1/2,1/2) requires that the absolute value of both eigenvalues be smaller than one, or
The inequality (5) yields constraints on the values of α_{X} and α_{Y} for which, given s_{X }and s_{Y}, the equilibrium (x^{∗},y^{∗}) = (1/2,1/2) is unstable, resulting in persisting phenotype frequency oscillations, or for which the equilibrium is stable, resulting in the cessation of the cycles.
The induced switching values for which the stability of the system is lost or regained can be calculated analytically for special cases of the stability condition (5), and otherwise either numerically or estimated from the simulation results. For example, when α_{X }= α_{Y }= α, s_{X }= s_{Y }= s, and μ = 0, the condition (5) is equivalent to
It can be shown that the lefthand side of the inequality (6) is an increasing function of s in the range of α for 0 ≤ s≤ 1, and 0 ≤ α ≤ 1. Therefore, an increasing selection strength will tend to induce cycles rather than destroy them. It can be also shown that the lefthand side of inequality (6) is an decreasing function of α for 0 ≤ α < 0.5, and an increasing function of α for 0.5 < α ≤ 1. Therefore, as induced switching α increases, cycles can be lost at low values of αand regained at high values of α. The inequality (6) can also be expressed as
which allows a precise calculation of the threshold levels of induced switching at which cycles disappear and reappear in this particular example.
To examine the persistence of cycles for a general case of α_{X }≠ α_{Y}, we solve the relation (5) numerically and compare it with the simulation results. Figure 2A shows the combinations of induced switching values α_{X }and α_{Y }for which oscillations dampen, with different combinations of selection coefficients, based on the stability condition (5). Figure 2B shows the analogous results which are extracted from the simulations. The presence or absence of allele frequency cycles in simulations is measured by the amplitude of the cycles after 11000. We define the presence of cycles if such amplitude exceeds the threshold of 5×10^{−2}, and otherwise we consider the cycles to be absent. A comparison between Figure 2A and 2B shows that simulation results and analytical predictions are in good accordance, and we see this in all examined regions of the parameter space. We also see that the results are largely independent of the value of the threshold used to measure the presence of allele cycles, provided the simulations run long enough.
Figure 2. Impact of the strength of selection on the persistence of cyclic phenotype frequency dynamics (A1–A4): Persistence of phenotype oscillations (black regions) for increased selection coefficients based on analytical predictions. When selection is weak (A1) cycles will continue indefinitely only for minimal or maximal levels of induced switching in both species. As selection becomes stronger (A2–A3), cyclic dynamics become more difficult to destroy and persist at higher levels of induced switching. For exceptionally high values of selection coefficients (A4), allele frequency oscillations will become unaffected by induced switching, including their speed and amplitude. As explained in the main text, the nature of cycles for minimal and maximal values of induced switching is different because of being driven by different evolutionary forces. (B1–B4): Numerical calculation of the amplitude of cycles from computer simulations confirms the analytical calculations in panels A1A4. It also reveals that the amplitude of cycles decreases as the cycles reemerge at high values of induced switching in both species X and Y . The following parameter values were used: (A1,B1) s_{X }= s_{Y }= 0.30, (A2,B2) s_{X }= s_{Y }= 0.57, (A3,B3) s_{X }= s_{Y }= 0.60, (A4,B4) s_{X }= s_{Y }= 0.90; panels A1A4 were made in resolution 101×101; in panels B1B4 were made in resolution 51×51, and cycles of amplitude 0.5×10^{−2} where considered nonexistent.
The results in Figure 2 illustrate a few important points. First, when both species undergo induced switching, lower rates of switching are needed to destroy the cyclic behaviour (cf. Figure 2A2). Second, as shown above, when the two species switch phenotypes at the same rate α = α_{X }= α_{Y}, the cycles can reemerge as α → 1. Finally, as also shown above, an increased selection pressure makes the cyclic dynamics more robust to higher levels of induced switching. In fact, as our calculations reveal, this dependence is so strong that for selection coefficient of 0.1 in both species as little as 0.5% of induced switching in species X is enough to eradicate the cyclic dynamics, while for selection coefficient of 0.9 induced switching will never dampen the cycles. The results for asymmetric selection coefficients are qualitatively identical.
Interestingly, the nature of cycles for low and high levels of induced switching is very different. In the case of α_{X }= α_{Y }= 0, the oscillatory behaviour will persist due to timedelayed negative frequencydependent selection (being rare is advantageous, being common is disadvantageous), whereas for α_{X }= α_{Y }= 1 oscillations will occur even in the absence of a selective force. The reason for this is that the latter situation represents the case where the phenotype frequency of one species in the next generation will be fully determined by the frequency of the phenotype of the other species. This will result in one species being constantly adapted to the other species population in the previous generation. However, since the other species does exactly the same, the two species will constantly crossreact even in the absence of any evolutionary force, leading to oscillatory “mirror dynamics” (see Discussion). Interestingly, in the parameter regime where these dynamics are dominant, increasing induced switching increases the amplitude of the allele cycles, while in the parameter regime where selectioninduced cycles are dominant, increasing induced switching decreases the amplitude of the allele cycle, as seen in Figure 2B1B4.
Stochastic switching can also destroy cyclic frequency dynamics. This is illustrated in Figure 3A1A4. As the stochastic switching rate μincreases, the parameter space with persistent phenotypic cycles shrinks down. This result is not surprising. Mathematically speaking, the absolute value of eigenvalues (4) decreases as μincreases, and therefore the stable state can be sometimes reached for high values of μ. This argument has been formulated some time ago in earlier theoretical hostparasite studies when discussing the potential impact of genetic mutation on antagonistic coevolution (see e.g., [23] and references therein). However, as genetic mutation rates are usually thought to be small (with some notable exceptions), they are typically not expected to cause any dampening effects in natural populations. Stochastic switching rates, however, can be orders of magnitude higher [16], and therefore the disappearance of cycles due to stochastic switching is likely to be more pervasive in nature than the disappearance of cycles due to genetic mutation.
Figure 3. Impact of stochastic phenotypic switching and asymmetric generation times on the persistence of phenotype oscillations. (A1–A4): Persistence and amplitude of cycles for different stochastic switching values. Subsequent panels show the results for increased values of μ, illustrating that high levels of stochastic switching can also destroy the persistence of the cyclic phenotype frequency dynamics. (B1B4): Persistence and amplitude of cycles for different asymmetric generation times between species X and Y . Subsequent panels show the results for increased values of g (number of Y generations per X generation). Rapid adaptation of one species also destroys cyclic dynamics, although the impact on cycles at high values of induces switching (especially α_{X}) is much less prominent. The following parameter values were used: (A1A4) s_{X }= s_{Y }= 0.66, (B1B4) s_{X }= s_{Y }= 0.75; (A1) μ = 10^{−8}, (A2) μ = 2×10^{−2}, (A3) μ = 3×10^{−2}, (A4) μ = 4×10^{−2}; (B1) g = 1, (B2) g = 2, (B3) g = 3, (B4) g = 10; all panels were produced in resolution 51×51, and cycles of amplitude 0.5×10^{−2}where considered nonexistent.
Finally, we examined the effect of asymmetric generation times between species X and Y , a situation that is certainly to be expected in hostparasite systems, and not uncommon in predatorprey systems. The results, shown in Figure 3B1B4, illustrate that an increased speed of evolution of species Y again makes the cyclic dynamics more sensitive to increased levels of induced switching. Interestingly, in this case the oscillatory “mirror dynamics” described above do not emerge for very high values of α_{X} and α_{Y}. This is because when one of the species evolves faster (here Y ), the symmetry of these dynamics is violated: species Y will always react more quickly, thereby immediately adapting to the other species. Furthermore, induced switching in the species which adapts more slowly (here X) has now a minor impact on the phenotype frequency dynamics observed in the model.
Discussion
Antagonistic coevolution is pervasive in nature, and oscillatory dynamics are generally thought to be one of its key signatures. The stability of this pattern is of fundamental importance in biology because the dynamics of phenotypes and genotypes are central to evolutionary and ecological processes. Furthermore, the absence of oscillations could be interpreted as the absence of an antagonistic interaction. We have shown here that in a simple model of antagonistic coevolution between two species, phenotypic switching – transmitted to the next generation through nongenetic inheritance – can have a dramatic effect on the patterns of antagonistic coevolution. Minimal levels of induced phenotypic switching can completely eliminate oscillatory dynamics and result in stable frequencies. This therefore suggests that even in the presence of strong links between the two species (i.e., strong selection, high specificity, etc.), antagonistic coevolution need not result in fluctuations of genotypes and phenotypes.
We have identified three parameters that affect the threshold level of induced switching at which cycles disappear. The first is the strength of selection in an antagonistic species interaction. For the threshold level to be high, both species need to suffer large fitness costs, to the extent that when selection is strong enough cycles will never be affected. Parasites may indeed pay such costs because their reproduction often depends on a successful antagonistic interaction with a host (see e.g., [24]). On the other hand, while both hosts and predators suffer fitness costs from being infected, or not being able to predate, their costs are arguably much lower. Second, an increase in g, the number of generations of the faster evolving species (e.g., the parasite) per generation of the other species (e.g., the host), typically reduces that threshold value. This is particularly relevant in the case of microparasites whose generation times can be many orders of magnitude shorter than that of their hosts. Finally, stochastic events affecting phenotypic switching can also reduce the threshold value. As stochastic switching events are increasingly being discovered in the microbial world, this effect might again be most relevant in the case of hostparasite interactions.
What makes the cycles disappear? Fundamentally, cycles depend on timelagged, negative frequencydependent selection (see e.g., [25]). Any factor that acts to reduce the timelag will act to reduce the amplitudes of cycles. In the absence of induced phenotypic switching, the speed at which the rare phenotype with a fitness advantage will increase in frequency depends on the strength of the antagonistic interaction. Lower fitness costs, higher discrepancy in generation times (i.e., higher g) and higher mutation rates all act to reduce the realised strength of interaction. For example, fast evolution in one species can lead to dampened cycles, masking interactions such that even though two species might be tightly linked (i.e., under strong selective pressure), the realized strength of interaction is low [26]. In the presence of induced phenotypic switching, there is limited scope for selection to reduce the frequency of the disadvantaged (common) phenotype; for example when induced switching occurs, counteradaptation occurs instantly at rate α, without the action of natural selection.
A more formal way to describe this phenomenon is to realize that under selection, the change of phenotype frequency in a single generation in one species, say X, depends on the variance of phenotype frequencies, x(1−x). This leads to a certain inertia characteristic of natural selection that produces the timelags in counteradaptation and thus the cyclical dynamics. Imagine a very common phenotype of X that is confronted with its equally common matching phenotype of Y , reducing its fitness. While this matching is obviously detrimental to this phenotype, the change in x in the next generation will be relatively small because individuals with this phenotype are mostly competing against individuals of the same phenotype (see upper right and lower left corners of Figure 4A). Under induced switching, however, this is not the case: under the same scenario, the change in x does not depend on the variance in phenotypes frequencies: a proportion α of the common X phenotype that is matched by Y will immediately switch to the other phenotype in the next generation (which is why the arrows are so much longer in the corners of Figure 4B). As a consequence, the timelag will be reduced and oscillations of phenotype frequencies x will gradually vanish. Similarly, high levels of stochastic switching will destroy oscillations in phenotype frequencies (Figure 4C).
Figure 4. Impact of the three evolutionary forces of our model on the change in phenotype frequency x and y during a single time step (one generation), as given by equations (1), and (2). The beginning and end of each arrow marks the phenotype frequencies before and after the respective step. As can be seen, selection tends to produce oscillations, whereas both induced and stochastic phenotypic switching tend to dampen these oscillation. Parameters used are (A) s_{X }= s_{Y }= 0.5, (B) α_{X }= α_{Y }= 0.2, and (C) μ = 0.1.
Overall, one of the most striking findings of this study is just how little phenotypic switching, especially interactioninduced, is necessary to completely eliminate cycles. One is tempted to speculate that such a process could be one of the reasons why evidence of dynamic polymorphisms is so rare, apart from the fact that longterm observations are difficult [7]. However, there are a number of caveats to consider. First, some antagonistic systems are characterised by strong selection [24], in which case we would expect that cycles would be maintained even in the presence of induced phenotypic switching. Second, evidence for induced phenotype switching as envisioned in this model is still rare, despite the fact that the number of demonstrations of strain or pathogenspecific immunity has been steadily increasing. Third, not every type of induced switching fits the implementation in our model. For example, maternal transfer of antibodies can make the offspring resistant to a pathogen strain encountered by the mother, but it does not come at the cost of becoming susceptible to another strain. However, such a tradeoff assumption is necessary for oscillations to appear in the first place – the model simply argues that if these tradeoffs do exist such that oscillations could be expected all else being equal, then phenotypic switching can dampen the oscillations altogether. Fourth, to what extent phenotype switching is stable across generations is currently largely unknown, and its adaptive value is an open question as well. Fifth, antagonistic fitness interactions are often resulting in fluctuating population densities, which may in turn affect themselves evolutionary dynamics [2630]. In order to understand the nature of the dynamics of phenotype frequency oscillations, we have purposefully ignored such population dynamics. Furthermore, how these results extend to complex communities of multiple species currently remains unknown. Finally, costs of induced switching may further reduce its dampening effects, provided that these costs are paid only by those individuals who are actually switching. Since we assume that only individuals affected by the interaction transmit the opposite phenotype, the impact of such costs can by easily calculated by multiplying relative fitness coefficient 1−s by 1−c, where c is the cost of switching. In the case of species X, the selection coefficient in eq. (1), s_{X}, would by substituted by s_{X} + c−s_{X}c. In contrast, in the presence of a general cost of maintaining a sensory mechanism for an antagonistic interaction, every individual would pay the same cost, and relative fitness would not be affected.
One of the important assumptions of this study is that the model underlying the antagonistic interaction is of a ‘matchingalleles’ type. Such a model is mostly applicable in the case of hosts with a specific immune system, and antigenic parasites, which have to specifically match the host in order to infect it. By contrast, interactions in many plantpathogen systems are usually thought to be of a ‘geneforgene’ type, where a host needs to recognise specific ‘effectors’ of the parasite in order to launch its defence [31]. In spite of this difference, the implications of this study bear similarity to the studies of plantpathogen models, where the conditions for the persistence of oscillatory dynamics and polymorphisms were thoroughly investigated. In particular, it has been previously noted that uncoupling of host and parasite life cycles in time or space can lead to a stabilization of allele cycles [32]. One good example is a high level of polycyclicity in a parasite life cycle, which was shown to induce stable polymorphism over time [33], in analogy to the results of our study (cf., Figure 3B). Analogously, high mutation rates can lead to stable equilibria of allele frequencies in plantpathogen systems [34,35]. Altogether, the analogies between the ‘matching allele’based systems and the ‘geneforgene’based systems point to the importance of empirical studies of nongenetic inheritance in both plantpathogen as well as animalparasite systems.
Conclusion
Environmentally induced phenotypic change that is stable across generations has recently been demonstrated in a number of cases, many of them involving stable epigenetic modifications [36,37]. Given the recent advances in this field, we expect many more demonstrations of these phenomena, and we see no obvious reason why they should not be observed in the realm of antagonistic interactions, especially since all species are likely to suffer severe fitness consequences if they are at the losing end of these interactions.
Competing interests
The authors declare no competing interests.
Authors’ contributions
RM conceived and designed the study, carried out the simulations, analysed the model and the results, and drafted the manuscript. JE analysed the results, and helped to draft the manuscript. MS conceived and designed the study, analysed the model and the results, and drafted the manuscript. All authors read and approved the final manuscript.
Acknowledgements
This work was supported by the Swiss National Science Foundation (RM & JE), and the Branco Weiss Fellowship (MS).
References

May RM: Stability and Complexity in Model Ecosystems. Princeton, NJ: Princeton University Press; 2001.

Hassell MP, Comins HN, May RM: Spatial structure and chaos in insect population dynamics.

Nuismer SL, Thompson JN: Coevolutionary alternation in antagonistic interactions.

SchmidHempel P: Evolutionary Parasitology: The Integrated Study of Infections, Immunology, Ecology, and Genetics. New York: Oxford University Press; 2011.

Abrams PA: The evolution of predatorprey interactions: theory and evidence.

Little TJ: The evolutionary significance of parasitism: do parasitedriven genetic dynamics occur ex silico?

Agrawal AA, Laforsch C, Tollrian R: Transgenerational induction of defences in animals and plants.

Agrawal AA: Phenotypic plasticity in the interactions and evolution of species.

Poulin R, Thomas F: Epigenetic effects of infection on the phenotype of host offspring: parasites reaching across host generations.

Hasselquist D, Nilsson JÅ: Maternal transfer of antibodies in vertebrates: transgenerational effects on offspring immunity.

Schulenburg H, Kurtz J, Moret Y, SivaJothy MT: Introduction. ecological immunology.

Little TJ, O’Connor B, Colegrave N, Watt K, Read AF: Maternal transfer of strainspecific immunity in an invertebrate.

LopezRubio JJ, Gontijo AM, Nunes MC, Issar N, Hernandez Rivas R, Scherf A: 5? flanking region of var genes nucleate histone modification patterns linked to phenotypic inheritance of virulence traits in malaria parasites.

Miller C, Thomsen LE, Gaggero C, Mosseri R, Ingmer H, Cohen SN: SOS response induction by betalactams and bacterial defense against antibiotic lethality.

Acar M, Mettetal JT, van Oudenaarden A: Stochastic switching as a survival strategy in fluctuating environments.

Lim HN, van Oudenaarden A: A multistep epigenetic switch enables the stable inheritance of DNA methylation states.

Rando OJ, Verstrepen KJ: Timescales of genetic and epigenetic inheritance.

Bonduriansky R, Day T: Nongenetic inheritance and its evolutionary implications.

Kouyos RD, Salathé M, Bonhoeffer S: The Red Queen and the persistence of linkagedisequilibrium oscillations in finite and infinite populations.

Woolhouse MEJ, Webster JP, Domingo E, Charlesworth B, Levin BR: Biological and biomedical implications of the coevolution of pathogens and their hosts.

Nee S: Antagonistic coevolution and the evolution of genotypic randomization.

King KC, Jokela J, Lively CM: Trematode parasites infect or die in snail hosts.

Bell G: The Masterpiece of Nature: The Evolution and Genetics of Sexuality. Berkeley: Croom Helm and University of California Press; 1982.

Yoshida T, Ellner SP, Jones LE, Bohannan BJM, Lenski RE, Hairston NG: Cryptic population dynamics: rapid evolution masks trophic interactions.

Abrams PA: The effects of switching behavior on the evolutionary diversification of generalist consumers.

Mougi A, Kishida O: Reciprocal phenotypic plasticity can lead to stable predatorprey interaction.

Mougi A, Iwasa Y: Evolution towards oscillation or stability in a predatorprey system.

Mougi A, Kishida O, Iwasa Y: Coevolution of phenotypic plasticity in predator and prey: why are inducible offenses rarer than inducible defenses?

Thompson JN, Burdon JJ: Geneforgene coevolution between plants and parasites.

Brown JK, Tellier A: Plantparasite coevolution: bridging the gap between genetics and ecology.

Tellier A, Brown JK: Stability of genetic polymorphism in hostparasite interactions.

Kirby GC, Burdon JJ: Effects of mutation and random drift on leonard’s geneforgene coevolution model.

Segarra J: Stable polymorphisms in a twolocus geneforgene system.

Kaufmann BB, Yang Q, Mettetal JT, Van Oudenaarden A: Heritable stochastic switching revealed by singlecell genealogy.

Jablonka E, Raz G: Transgenerational epigenetic inheritance: prevalence, mechanisms, and implications for the study of heredity and evolution.