Host-parasite Red Queen dynamics with phase-locked rare genotypes

Abstract

INTRODUCTION

Antagonistic interactions are major factors that drive evolution on different time scales (1–4). The Red Queen hypothesis has been proposed as a model for antagonistic interactions where species (for example, host-parasite, prey-predator, and victim-exploiter) perpetually coevolve in winnerless dynamics (1, 2, 5, 6). In host-parasite interactions, the Red Queen hypothesis suggests that coevolution occurs as a result of time-lagged negative frequency-dependent selection. The fitness of a common host type decreases because of overinfection by a common parasite type, prompting the increase in the fitness of a rare host type, which after its increase becomes the target of other parasite types that will then increase in consequence (1, 3, 5, 7, 8). Winnerless coevolution is widespread in host-parasite interactions because of nearly symmetrical selection, which means that the evolution of hosts is countered by the evolution of parasites (1, 9). Species that fail to catch up with such coevolutionary races may become extinct (4, 10). The Red Queen hypothesis has gained an important position in evolutionary biology, being suggested to explain the evolution of sex (1, 11), the antagonist-mediated diversity of species (12–14), and the emergence of antibiotic resistance (7).

One of the manifestations of the Red Queen hypothesis is oscillatory dynamics, mathematically defined as out-of-phase population/frequency cycles with similar amplitude (1, 5, 15, 16). In canonical Red Queen dynamics (5), all of the host and parasite genotypes undergo negative frequency-dependent selection (represented by the out-of-phase cycles), but their collective fitness remains the same (fig. S1). However, there are cases where only some genotypes are engaged in Red Queen dynamics, whereas others are engaged in phase-locked (synchronized) cycles with minor amplitude (Fig. 1). Here, we investigate the conditions for such cases using mathematical simulations involving multihost and multiparasite types, and we test for the effects of stochastic noise on coevolutionary dynamics. Stochastic noise represents random physical-environmental perturbations introduced in simulation parameters (17–21). Cycles of alternating dominance in the Red Queen dynamics and the effects of noise on coevolution can be associated with the dynamics observed in marine microbial communities (for example, bacteriophage predation), which have been explained by the “killing-the-winner” hypothesis (15, 22, 23). The Red Queen dynamics also explain coevolution in invertebrate-parasite systems (for example, bacterial infection of Daphnia) (8, 24, 25).

Open in a separate window

Fig. 1

Population density time series illustrating Red Queen dynamics with phase-locked rare types (K = 40, σ = 0).

Host 1 (parasite 1) and host 2 (parasite 2) have out-of-phase cycles with dominant amplitude, whereas the rest are inferior and phase-locked. We call the cyclic dominance between the two dominant species types as “Red Queen binary oscillations.” Some inferior phase-locked population plots are not visible because they have almost similar periods and amplitudes as the time series plotted in front of them. The qualitative behaviors of both host and parasite population time series follow Red Queen binary oscillations. (A) n = 3. (B) n = 5. (C) n = 10.

We hypothesize that functional response in host-parasite interactions influences the occurrence of Red Queen dynamics. Functional response represents the rates of infection and host utilization by parasites (26, 27). Type I functional response denotes linear parasitic utilization of hosts. Type II functional response denotes parasitic utilization of hosts following a hyperbolic curve, whereas type III functional response follows a sigmoidal (S-shaped) curve. Type III functional response incorporates host switching and parasite learning, which are characterized by a low parasitism infectivity rate at low host density and by an accelerated parasitism infectivity rate at intermediate host density (27–29). Here, our simulations using type III functional response reveal a novel pattern of coevolution exhibiting the Red Queen dynamics of two host and parasite types, whereas several other genotypes remain subordinate for long periods of time. Discovery of novel coevolutionary dynamics occurring as a result of complex biotic interactions provides an insight into the phenomena seen in epidemiology and evolution, especially the evolution of resistance and virulence (30–32).

RESULTS

To simulate the interaction of n host types and n parasite types (genotypes, strains, varieties, or species), we use an ordinary differential equation representation of an antagonistic system with a type III sigmoidal functional response (see Materials and Methods). We have observed that a limited number of combinations of parameter values lead to Red Queen cycles. A high host basal growth rate (r_i), an intermediate death rate in parasites (d_j), and high parasite specificity (for example, matching-allele interactions) favor the occurrence of Red Queen oscillatory behavior (Fig. 2 and fig. S2). An increased host basal growth rate provides the host the ability to recover from the effects of parasitism, which is a requirement for inducing Red Queen cycles (Fig. 2A). On the other hand, an intermediate degree of parasite mortality is essential for generating Red Queen cycles because a low parasite death rate results in high levels of parasitism that adversely affect the growth of host populations (Fig. 2B). A high parasite death rate limits the antagonistic effect of the parasites, which is not enough to steer negative frequency-dependent selection in hosts, resulting in equilibrium (Fig. 2B). Furthermore, host-parasite interactions with high parasite specificity are favorable for driving cyclic and never-ending negative frequency-dependent selection, compared to interactions with more uniform parasite specificity (fig. S2).

Open in a separate window

Fig. 2

Effects of host basal growth rate (r_i = r) and parasite mortality rate (d_j = d) on the generation of Red Queen binary oscillations (K = 40, σ = 0).

Vertical alternating colors characterize oscillations. (A) A high host basal growth rate r favors the occurrence of Red Queen binary oscillations (for example, r > 0.5; d = 0.5). A lower basal growth rate may result in host extinction or in population dynamics where only one type is dominant. (B) An intermediate parasite mortality rate d favors the occurrence of Red Queen binary oscillations (for example, d = 0.5; r = 1). Relatively low and high death rates may result in host extinction or in population dynamics where only one host type is dominant.

Red Queen dynamics with noise

To test the sensitivity of the deterministic host-parasite interaction model against physical-environmental perturbations, we introduce random noise to the model parameters following the Ornstein-Uhlenbeck process (see Materials and Methods). Under stochastic conditions, it is possible that one of the dominant host/parasite types will be replaced by a subordinate type as a new major player in the Red Queen dynamics (Figs. 3 and ​and4).4). Such a replacement is followed again by Red Queen binary oscillations but with different dominant cyclic types. The occurrence of such transition events becomes more frequent as the degree of physical-environmental perturbations increases (figs. S5 and S6). Moderate degrees of stochastic noise preserve the Red Queen binary oscillations; however, large perturbations can cause extinction of species (fig. S5). The replacement of a dominant type by a subordinate type in the Red Queen dynamics is due to changes in the characteristics of each host type (for example, growth rate) and each parasite type (for example, death rate and specificity).

Open in a separate window

Fig. 3

Illustration of Red Queen dynamics with noise (n = 3, K = 25, σ = 15%).

Here, Red Queen dynamics are represented by winnerless coevolution between two dominant cyclic host/parasite types (binary oscillations). The evolution due to random noise is characterized by a stochastic transition event where a rare cyclic host/parasite type (H₃, P₃) replaces one of the dominant cyclic types (H₂, P₂). The transition event is a result of the accumulated effect of random noise representing uncertain physical-environmental perturbations.

Open in a separate window

Fig. 4

Evolutionary landscape separated into two regions by an entrapment barrier.

The height of the landscape characterizes the population density of species types (balls). The barrier limits the number of dominant cyclic types and constrains subordinate/rare cyclic types to enter the optimal region where Red Queen dynamics occur. In the optimal region, a dominant cyclic type can attain maximum growth potential near its carrying capacity. The population dynamics of the species oscillate (phase-locked subordinate cycles and Red Queen dynamics) as a result of parasitism (the paths of broken arrow lines illustrate oscillatory behavior). (A) Rare cyclic types remain rare because they are trapped at a suboptimal region. Only two species types undergo cyclic dominance of Red Queen dynamics. The population dynamics of rare cyclic types do not exhibit Red Queen dynamics but rather are phase-locked (synchronized). (B) Physical-environmental variability can allow a rare cyclic type to escape the suboptimal region and enter the optimal region by hurdling the evolutionary entrapment barrier. However, the optimal region only allows two dominant cyclic types; hence, one of the types needs to leave the optimal region to accommodate the succeeding rare type (red and green balls).

DISCUSSION

On the basis of our simulations, replacement of a dominant cyclic host type can happen in two ways: through negative frequency-dependent selection following the Red Queen dynamics and through noise-driven replacement of a dominant cyclic host type by a subordinate cyclic type. We can look at these two ways of evolutionary switching in an evolutionary landscape (Fig. 4). In a controlled and deterministic environment, several subordinate/rare cyclic types can be trapped at a suboptimal region in the evolutionary landscape (Fig. 4A). The optimal region in the landscape denotes the situation where a dominant cyclic type can attain maximum growth potential near its carrying capacity. Physical-environmental variability can aid the trapped subordinate cyclic types to escape and move to a new state of dynamics by hurdling the evolutionary entrapment barrier (Fig. 4B). The entrapment barrier is a constraint that only allows two dominant cyclic genotypes. The formation of the barrier is influenced by the characteristics of hosts and parasites and their interactions, such as the growth rate of hosts, the death rate of parasites, parasite specificity, and the competitive ability of host types. The evolutionary switching events can elucidate the mechanisms that drive different ecological and evolutionary phenomena, such as the fluctuating population in marine microbial communities (15, 22, 23, 33), the evolution of water fleas (8, 24, 25), and the abundance of rare species in nature (34, 35).

The structure of parasitism functional response is essential for generating Red Queen dynamics. In type II and type III functional responses, the rates of infection and host utilization by parasites tend to satiate at high host density. However, type III functional response incorporates population size–based parasite learning and host switching (27–29, 36). This is indicated by the case where a host type with intermediate density is more exposed and prone to parasite infection than a host type with lower density. As a result, certain parameter ranges of the model with type III functional response accommodate only a limited number of dominant cyclic host/parasite types, and the rest are subordinate cyclic types (Fig. 1). This finding is in contrast to the dynamics of models with type II functional response, where the Red Queen dynamics arise in every host type and every parasite type given specific parameter values (we call these canonical Red Queen cycles as “Red Queen n-ary oscillations”) (5). Red Queen n-ary oscillations allow all hosts/parasites to undergo never-ending dominance switching, such that each type can be numerically dominant for certain periods of time (fig. S1). In addition, an interaction system involving multihosts and multiparasites with linear functional response (type I) may result in a different scenario where out-of-phase population cycles with identical amplitude do not exist (37).

With a type III functional response, the rate of infectivity is low at low host population densities, but increases at intermediate host population densities and then slows down at very high host population densities as a result of satiation. Various host-parasite interactions in nature conform to this sigmoidal (S-shaped) functional response curve (38–41). The sigmoidal functional response curve could be a logical consequence of mass-action transmission as it is commonly assumed to be valid in host-parasite systems (38–40). Mass-action transmission implies that the number of susceptible hosts that become infected by parasites depends on the host and parasite population densities.

Consequently, both the Red Queen binary oscillations (resulting from type III functional response) and the Red Queen n-ary oscillations (resulting from type II functional response) bring about parasite-regulated temporal biodiversity (12–14). The Red Queen dynamics, thought to be a form of balancing selection (42), reduce the likelihood of extinction. However, binary and n-ary oscillations differ in the number and dynamics of dominant cyclic host/parasite types. Red Queen binary oscillations entail biodiversity where only a limited number of species types can dominate the system and the rest are phase-locked rare types. The formation of binary oscillations, a special mode of Red Queen dynamics, is one of the possible explanations for the “excess of rare species” phenomenon (34, 35) and synchrony in population abundance (43, 44). Moreover, we hypothesize that restructuring the model of interhost competition and parasitism might give rise to m (where 2 < m < n) sequential dominant host/parasite types.

The relationship between interhost competition and parasitism is necessary for the occurrence of Red Queen binary oscillations. Parasitism is not enough to generate out-of-phase oscillations, but rather strong rivalry among host types should also occur (5). In host-parasite coevolution, interhost competition indicates evolutionary selection among host genotypes. In addition to having increased host basal growth rate (Fig. 2), the resources (K) of competing hosts should be large enough to provide each host the ability to attain its full growth potential (figs. S3 and S4). Hosts with high potential for growth are suitable for stimulating effective parasitism that eventually results in negative frequency-dependent selection. Moreover, increasing the host’s carrying capacity generates oscillations with large amplitude, an effect earlier described as the paradox of enrichment (45).

Random noise could represent stochastic physical-environmental perturbations (17, 20). A moderate level of stochastic noise causes one host/parasite type to leave the Red Queen dynamics and another type to enter the winnerless coevolution process (Fig. 3). Red Queen binary oscillations warrant the existence of subordinate host and parasite types, as well as a limited number of dominant cyclic types, whereas random noise provides an opportunity for the subordinate types to replace the dominant ones (Figs. 3 and ​and4).4). A moderate level of stochastic noise affects the model parameters (for example, growth rate and carrying capacity of hosts, mortality rate, specificity, and infectivity of parasites), driving the previously dominant cyclic type to decline and allowing a subordinate cyclic host type to become a new dominant type. The noise-driven changes in parameter values affect the strength of interactions among the hosts and the parasites, resulting in the alteration of players in the Red Queen dynamics. Increasing the degree of physical-environmental perturbations is then expected to increase the occurrence of transition events where a rare cyclic host/parasite type replaces one of the dominant cyclic types (fig. S5).

Many cases of evolutionary switching of dominant populations occur because natural ecosystems are inherently noisy (17–21). Stochastic noise is believed to play a role in regulating the dominant regime in ecological interacting systems (17, 21, 46). Hence, the consequence of the Red Queen binary oscillations with physical-environmental perturbations might be the temporal biodiversity seen in natural ecosystems (for example, see fig. S6) (18, 19, 21). The interplay among biotic and abiotic factors under suitable conditions (genotype-by-genotype-by-environment) shapes the structures and interactions in nature (47–52). However, note that very large noise does not necessarily result in fluctuating densities (fig. S5). Large abiotic perturbations could be outside the survival tolerance range of some species, causing extinction of species (53–55).

In summary, Red Queen binary oscillations illustrate the parasitism-induced evolution of hosts/parasites from one dominant cyclic type to another type (Fig. 1). The population cycles of the two dominant cyclic types have the same amplitude, characterizing a consistent fitness pattern in spite of the continuous evolutionary switching between the two types. The cycle of evolution persists between the two dominant types only, and the remaining types are rare. However, a phase-locked rare type can evolve to become dominant as a result of a considerable degree of physical-environmental noise (Fig. 3). Our mathematical simulation explains the underlying mechanisms of one of the manifestations of Red Queen dynamics acting on host-parasite coevolution. This type of coevolutionary dynamics presents an explanation for the abundance of rare species. Our numerical prediction can guide future empirical studies in finding more evidence for the Red Queen hypothesis and in investigating how winnerless coevolution affects complex ecological interactions in nature (6, 56–58).

MATERIALS AND METHODS

We considered a resource-consumer model of the host-parasite interaction system (5, 29, 37) with up to 10 host types and 10 parasite types (n = 10). The mathematical model with type III functional response is

where i,j = 1,2,…,n (n ≥ 3).

In Eqs. 1 and 2, the population size of a host type (H_i) increases according to basal growth rate (r_i) and decreases as a result of parasitic infections (functional response). The sum of all hosts shares a carrying capacity of size K. On the other hand, the population of a parasite type (P_j) increases with parasitic utilization of hosts and decreases with constant death (d_j) (numerical response). Let the coefficients φ_ik, c_ij, and β_jik (k = 1,2,…,n) represent the competitive ability between hosts i and k, host i to parasite j conversion due to host utilization, and the effect of host type k (k = 1,2,…,n) on the alteration of the functional response associated with the infection of host i by parasite j, respectively (table S1). This model assumes that hosts are immediately eliminated once infected by a parasite. Infected hosts cannot recover from disease, as is often the case in biological systems involving bacteriophages and castrating parasites and in invertebrate-parasite systems (for example, infection of Daphnia magna by Pasteuria ramosa) (8, 37, 59–62).

The complexity of the dynamics of species interaction increases with the number of variables and parameters. We accordingly focused on the conditions under which Red Queen dynamics are likely to occur in restricted parameter space (see the Supplementary Materials for the discussion of simulation parameter values). Various qualitative dynamics and other modes of Red Queen cycles may arise in an interaction system with type III functional response (for example, fig. S7), but we focused on the parameter conditions under which the Red Queen binary oscillations are likely to occur.

The parameter α_ij in the functional response represents the parasitism efficiency of parasite j against host i. Because we were concerned with patterns of oscillations (for example, Red Queen dynamics), we restricted the parasitism efficiency matrix A = [α_ij] to be diagonally dominant such that

for all i.

A diagonally dominant matrix implies that host i is the main host of parasite j = i, whereas all hosts k ≠ i are secondary or nonhosts. This assumption assures that there is neither a unique superior parasite nor an inferior parasite in infectivity, and that there is neither a unique superior host nor an inferior host in parasite resistance. Trade-offs restrict the presence of all-resistant hosts and of superparasites that can infect all hosts (22, 63, 64). The parasitism efficiency matrix can be considered as a representation of the continuum between gene-for-gene and matching-allele models (5, 65–68).

The use of parameter perturbations following an Ornstein-Uhlenbeck process represents physical-environmental stochastic noise. The assumed stochastic differential equation for each variable parameter, say x, is formulated as (20)

The term characterizes drift toward the mean value . The parameter σ_x is the amplitude of the Gaussian white noise dW_x. The variable parameters in our simulations are r_i ≥ 0, d_j ≥ 0, K > 0, and α_ij ≥ 0 for all i and j. Each parameter is perturbed by independent white noise.

Supplementary Material

Acknowledgments

SUPPLEMENTARY MATERIALS

REFERENCES AND NOTES