Lotka–Volterra dynamics kills the Red Queen: population size fluctuations and associated stochasticity dramatically change host-parasite coevolution
© Gokhale et al.; licensee BioMed Central Ltd. 2013
Received: 8 August 2013
Accepted: 13 November 2013
Published: 19 November 2013
Host-parasite coevolution is generally believed to follow Red Queen dynamics consisting of ongoing oscillations in the frequencies of interacting host and parasite alleles. This belief is founded on previous theoretical work, which assumes infinite or constant population size. To what extent are such sustained oscillations realistic?
Here, we use a related mathematical modeling approach to demonstrate that ongoing Red Queen dynamics is unlikely. In fact, they collapse rapidly when two critical pieces of realism are acknowledged: (i) population size fluctuations, caused by the antagonism of the interaction in concordance with the Lotka-Volterra relationship; and (ii) stochasticity, acting in any finite population. Together, these two factors cause fast allele fixation. Fixation is not restricted to common alleles, as expected from drift, but also seen for originally rare alleles under a wide parameter space, potentially facilitating spread of novel variants.
Our results call for a paradigm shift in our understanding of host-parasite coevolution, strongly suggesting that these are driven by recurrent selective sweeps rather than continuous allele oscillations.
KeywordsHost-parasite coevolution Red Queen hypothesis Lotka-Volterra dynamics Genetic drift Population bottleneck
The Red Queen from Lewis Carroll’s tale ‘Through the looking glass’ is commonly used as a metaphor for selection-induced rapid evolution [1–3]. It is based on the observation that persistence in an environment with changing selective constraints requires ongoing adaptations to the encountered challenges . Host-parasite coevolution with antagonistic and inter-dependent interactions represents one of the role models for such rapid evolutionary change [5, 6]. For instance, an increase in host resistance reduces parasite fitness, thus immediately favoring parasite varieties with altered virulence and/or immune-evasion mechanisms. In turn, a novel parasite attack mechanism decreases host fitness, thus favoring host varieties with new counter-defenses. If the interaction persists, then it will lead to continuous parasite adaptations and host counter-adaptations. The rapid evolutionary dynamics associated with these interactions is very well documented in the literature, ranging from field studies on rabbits and their myxoma viruses , snails and their trematode parasites , Daphnia magna waterfleas and their bacterial parasites  to laboratory-based coevolution experiments between bacteria and their phages [10–12], the nematode Caenorhabditis elegans and bacterial parasites [13, 14], or the red flour beetle Tribolium castaneum and its microsporidian parasite [15, 16].
It is thus widely accepted that these interactions evolve fast and continuously. Yet, to date, the exact underlying selection dynamics are not always well understood. These dynamics can generally be influenced by metapopulation structure and environmental variation [17, 18]. Within a particular population and specific environmental context, two main alternatives are thought to be of prime importance: recurrent selective sweeps and negative frequency-dependent selection [5, 6, 19–21]. Both alternatives are consistent with the above original definition of the Red Queen hypothesis by Van Valen , whereas, curiously, only the second alternative is referred to as Red Queen dynamics [5, 6, 20]. The two alternatives are closely related because both assume a selective advantage of a rare genotype, for example a novel host resistance variant. However, they differ fundamentally in the way in which the new variant originates and spreads within the population. The concept of recurrent selective sweeps (often termed arms race dynamics) consists of two steps: the de novo appearance of a beneficial allele (e.g., by mutation or immigration) and its subsequent spread through the population to fixation (i.e., the selective sweep). These sweeps occur repeatedly in host and parasite populations, usually each time with a new beneficial allele. They may only lead to fast changes in absolute time if at least one of the following factors applies: new alleles arise frequently, new alleles become immediately visible and thus selectable at the phenotypic level, the new alleles provide a high selective advantage, and/or the organisms have short generation times. This situation is best met in bacteria-phage interactions, which are usually characterized by large population sizes (i.e., high likelihood of the occurrence of favorable mutations), short generation times, and haploid genomes (i.e., new mutations are immediately expressed phenotypically) [11, 22–24] (but see also ).
In contrast, the dynamics for multicellular host systems are traditionally viewed to be determined by negative frequency-dependent selection leading to sustained oscillations of the same alleles (i.e., Red Queen dynamics [6, 20]), but not to the fixation of single alleles. In this case, standing genetic variation is required, because the population sizes for these hosts are usually comparatively small, their generation times comparatively long, and their genomes diploid. As a consequence, recurrent selective sweeps are commonly thought to be rather slow in these systems. Instead, if standing genetic variation is available, then negative frequency-dependent selection can produce fast and continuous allele frequency changes even in these host systems. Such negative frequency-dependent dynamics seem to be present in some multicellular host systems, including the freshwater snail Potamopyrgus antipodarum[8, 26] and the waterflea Daphnia magna.
Numerous theoretical models have been developed to study the underlying selection dynamics. Interestingly, the current models typically focus on evolutionary change (i.e., the rate of change in host and parasite allele frequencies in response to the type of interaction). These approaches have thus largely neglected ecological dynamics, which can have a huge impact on the evolutionary process. Population size fluctuations deserve particular attention in this context, because they are induced by reciprocal selection among the antagonists and, therefore, represent an inherent property of host-parasite coevolution - irrespective of additional environmental variation [7, 10, 27–30]. Since selection is reciprocal, population size fluctuations should be coupled between the antagonists, and generally follow Lotka-Volterra dynamics [31, 32]. Such demographic variations have the potential to affect the dynamics of host-parasite allele frequency changes by introducing two important effects. Firstly, the rising and falling population sizes produce bottlenecks where selection favours a particular allele. The favored allele may thus reach comparatively high frequencies during the bottleneck, possibly enhancing its spread in the subsequently expanding population. Secondly, the elevated stochasticity during the bottleneck may lead to a further increase and thus spread of the favored allele.
In this manuscript, we aim at understanding in how far Lotka-Volterra population size fluctuations and the associated stochastic effects influence the dynamics of allele frequency changes during host-parasite coevolution. While several previous theoretical models have applied the Lotka-Volterra dynamics to host-parasite coevolution (e.g., [33–37]), their influence on the evolutionary dynamics has not yet been systematically explored by comparison with a model with constant population size. Similarly, stochastic effects during host-parasite coevolution have only been considered in a few theoretical studies (e.g., [38, 39]), yet, to our knowledge, with a single exception  under constant population size and not in combination with Lotka-Volterra dynamics. Hence, while the previous studies have independently utilised stochastic effects or Lotka-Volterra dynamics, a systematic analysis of the consequences of each of these factors, either alone or in combination, is as yet missing - in spite of their potential importance. The novelty of our study lies in bringing together these two aspects and comparing their influence to the traditional model, in which Lotka-Volterra dynamics and stochastic effects are excluded. More specifically, we here use the standard matching-alleles host-parasite interaction model to assess allele frequency dynamics in the presence versus absence of Lotka-Volterra oscillations for a stochastic versus an analogous deterministic model.
where the frequencies of H i and P i are given by h i and p i . The above equations consider interdependence of host and parasite demographies, allowing population sizes to vary in response to the interaction with the antagonist, consistent with the Lotka-Volterra model. The precise nature of the resulting oscillations in population size is determined by the parameters, most importantly by b.
As we are interested in the effects of population size variation induced by the Lotka-Volterra equations, we have to compare this to a scenario in which the population size is constant. Such constant population size models are common, e.g. the Wright-Fisher model or the Moran process. However, microscopically these models are distinct from the Lotka-Volterra equations considered above. Therefore, we used the above approach and enforced constant population size by resetting host and parasite population sizes to their initial values after every generation (Navg transition events, see Appendix), while relative allele frequencies were maintained. The dynamics were subsequently assessed for different average population sizes. To ensure comparability of allele frequency fluctuations across population sizes and evolutionary models, we rescaled the interaction parameters with Navg for the deterministic analogues of the considered stochastic scenarios (Appendix).
The results also highlight that the dynamics are usually determined by fixation of one of the host alleles (red colour is mainly found in middle rather than right panel of the top row of Figure 3). Note that the simulations are stopped as soon as either one of the host or one of the parasite alleles reaches fixation and, thus, the fixation probabilities of both host as well as both parasite alleles sum up to one. In our case, fixation of the host allele is more likely than fixation of the parasite allele because for our parameter combination and initial condition, it is usually the host that first experiences a Lotka-Volterra bottleneck and consequentially a drop in the frequency of one of the alleles (see also Figure 4). Nevertheless, if both the parasite and corresponding host allele are common, then it is the parasite allele that has a high probability of fixation (Figure 3, top right).
The overall pattern looks different in the absence of Lotka-Volterra fluctuations (Figure 3, bottom row). Host allele fixation probability increases with its own high initial frequency and, at the same time, low initial abundance of the corresponding parasite allele (Figure 3 middle panel in bottom row). Parasite allele fixation is enhanced when both parasite and corresponding host alleles are initially common (Figure 3 bottom right). Under these conditions, fixation probabilities of both host and parasite alleles are almost identical at initially intermediate frequencies, most likely due to negative frequency dependent selection, as illustrated in Figure 1.
Population size fluctuations represent an inherent property of host-parasite interactions. Unequivocal evidence for such interaction-dependent demographic variations was obtained from controlled host-parasite mesocosm experiments, for example with E. coli and its phage [10, 30], Hydra hosts and its Hydramoeba parasite , house fly and its parasitic wasps [44, 45], or azuki bean weevil Callosobruchus chinensis and its parasitoid Heteropilus prosopidis. Similar observations were made under field conditions, for example for rabbits and their myxoma viruses , or red grouse and its nematode parasite . Additional examples are summarized by  and . As population size fluctuations produce regular bottlenecks, random genetic drift is likely to influence allele frequencies. Previous theoretical models, developed in a different context, strongly suggest that even large populations are influenced by such stochastic processes [48, 49]. More generally, under natural conditions in a finite population, it is difficult to imagine that changes in population size do not affect evolutionary dynamics. Consequently, an in-depth understanding of the evolution of host-parasite interactions should take account of the associated ecological processes based on Lotka-Volterra fluctuations.
Very few previous theoretical models on host-parasite interactions have considered Lotka-Volterra fluctuations [33–37]. These studies usually used a deterministic approach and thus excluded stochastic effects, which are most prominent during bottlenecks. Similarly, only few theoretical studies considered stochastic effects in this context [38, 39], yet under constant population size, but not in combination with Lotka-Volterra dynamics. We are aware of only one study that looked at host-parasite coevolution in consideration of Lotka-Volterra interactions and stochasticity . However, this study had a different focus, and thus, it did not include a systematic comparison to models without Lotka-Volterra cycles or without stochasticity. Consequently, the interaction of these two aspects for host-parasite coevolution is so far unexplored. At the same time, their relevance was demonstrated for evolution of only one of the antagonists, namely the parasite. For example, the probability of fixation of a beneficial mutation in a bacterial population was shown to be enhanced by periodical bottlenecks [50–52]. Similar results were obtained in a model that explored the effect of bottlenecks during pathogen transmission . Our study explicitly evaluates the influence of both Lotka-Volterra fluctuations and stochastic effects on the dynamics of host-parasite interactions using a comparison to a model with constant population size and/or absence of stochasticity. As the demographic variations are an inherent property of such antagonistic interactions, their influence should apply across a wide range of environmental conditions and thus be of general relevance for our understanding of host-parasite coevolution.
Based on our approach, we obtained evidence that both Lotka-Volterra fluctuations and associated stochastic effects significantly affect the course and pace of coevolutionary adaptations. In particular, both factors facilitate selective sweeps (i.e., the spread and fixation of an allele). Most impressively, this effect appears to be independent of average population size (Figure 2) and occur at a substantially faster rate (Figure 3, left column). Moreover, allele frequency changes are not exclusively due to drift, which should favor fixation of initially frequent alleles and loss of initially rare alleles. In contrast, our results indicate that initially rare host alleles can spread to fixation across a relatively wide range of conditions (Figure 3, top middle panel). Rare parasite alleles may not necessarily go extinct, but have a certain likelihood of spreading contingent on the frequency of the corresponding host allele (Figure 3, top right panel). Based on these results we propose that selective sweeps rather than oscillatory negative frequency-dependent selection may represent the main driving force during host-parasite coevolution.
Recurrent selective sweeps have been repeatedly suggested to determine coevolutionary dynamics for parasite or host systems with large population sizes such as bacterial hosts or microbial parasites, where novel mutations are frequent and often directly exposed to selection because of a haploid genetic system. If these selective dynamics also apply to multicellular host and parasite systems, then two contrasting effects may be expected on the coevolutionary process. On the one hand, these systems usually have much smaller population sizes, facilitating spread of alleles in spite of the often diploid genetic system. On the other hand, continuous coevolution may become difficult because it is usually assumed that small population size results in a reduced likelihood of the occurrence of advantageous novel mutations . However, the latter assumption may not always be true. It is possible that new alleles become available for example by frequent immigration or a high rate of gene duplication. These processes may further favor the formation of novel genotypes if they act in combination with recombination and/or mutation. Interestingly, the possible impact of gene duplications is usually not addressed in theoretical work on host-parasite coevolution, even though such duplications are known to be common in almost all organisms [54–56] and often affect genes of relevance for the interaction such as virulence genes in parasites [57, 58] or immunity genes in animal hosts [54, 59, 60].
Several additional factors may favor ongoing coevolution. One of these is founded on a more complex genetic architecture underlying host-parasite co-adaptation that consists of several interacting loci across the genome (e.g., [61–64]). In this case, allele fixation at one locus may still permit maintenance of variation at other loci, which could then mediate evolutionary responses to the antagonist. Yet another possibility may depend on metapopulation structure, consisting of coevolutionary hot and cold spots and migration among demes, as evidenced for flax and its parasitic rust fungus  or the above mentioned snail-trematode interaction . Such an interconnected network could then maintain allelic diversity across the entire metapopulation, even if alleles become temporarily fixed within single demes. Moreover, environmental gradients or perturbations are known to influence host-parasite coevolutionary dynamics . They could similarly prevent loss of alleles, even if the coevolutionary interaction itself would specifically favour only one of the alleles. Obviously, the above processes act in combination with each other in natural systems. Therefore, it is indeed conceivable that recurrent selective sweeps shape long-lasting coevolutionary dynamics even in multicellular host systems.
In conclusion, decades of empirical efforts have tried to demonstrate the presence of Red Queen dynamics during host-parasite coevolution. This has led to most ingenious experiments which repeatedly and independently confirmed negative frequency dependence as a driving force [8, 9, 26, 67, 68] and such a trend continues to date [21, 69]. These studies yielded impressive evidence that parasite abundance typically increases first and, once the host evolves a defense mechanism, it decreases again. However, sustained allele frequency oscillations of a particular allele, as predicted by numerous theoretical models assuming constant population size in the absence of any stochastic effects, have not been reported. We here propose that Lotka-Volterra population size fluctuations and the associated stochastic effects represent an inherent property of host-parasite interactions that can lead to rapid fixation of alleles, even those initially rare, thus preventing sustained oscillations. Consequently, Lotka-Volterra population size fluctuations have the potential to stop the Red Queen - unless novel variants are introduced into the population and/or additional selective constraints maintain allelic diversity. In retrospect, our findings may not be entirely unexpected. However, to date, they have not yet been directly demonstrated using a systematic analysis approach, as implemented here. More importantly, they are generally neglected in the numerous current empirical studies on the topic, in spite of their potential importance. They clearly deserve specific attention in future theoretical and empirical work aimed at an improved understanding of host-parasite coevolutionary dynamics.
Relating stochastic and deterministic dynamics
Stochastic models are often developed starting from deterministic formulations . Since the same deterministic formulation can be the limiting case of many individual based models, this procedure may be problematic. Instead, beginning from a stochastic, individual based description and then calculating the deterministic analogue will provide only a single direct link between the two approaches and allows for their direct comparison.
The Gillespie algorithm gives an exact numerical solution of the Master equation of the system [41, 70, 71]. Our stochastic simulations are implementations of this algorithm with the transition rates as defined in Eqs. 6. Since the population size is not constrained, this simulation method includes a stochastic analogue of the Lotka-Volterra cycles.
We computationally remove the Lotka-Volterra cycles by culling the population of each species after Navg transitions have taken place. During the Navg transitions the types within a species can evolve to different frequencies. But in the end they are reset to sum up to Navg while maintaining the relative abundances. The Gillespie method is discrete in the number of individuals but continuous in time. The unit of time is the same as in the deterministic system.
CSG and AT acknowledge financial support from the Emmy-Noether program of the Deutsche Forschungsgemeinschaft and the Max Planck Society; AP and HS from the Deutsche Forschungsgemeinschaft within the priority programme SPP1399 on host-parasite coevolution (grant SCHU 1415/8). We thank Christian Hilbe for helpful discussions.
- Carroll L: Through the Looking-glass, and what Alice Found there. 1871, London: Macmillan, Reprinted: Bloomsbury, London, 2001Google Scholar
- Ridley M: Evolution,. 2003, Oxford: Wiley-BlackwellGoogle Scholar
- Barton NH, Briggs DEG, Eisen JA, Goldstein DB, Patel NH: Evolution,. 2007, NY: Cold Spring Harbor Laboratory PressGoogle Scholar
- van Valen L: A new evolutionary law. Evol Theory. 1973, 1: 1-30.Google Scholar
- Woolhouse MEJ, Webster JP, Domingo E, Charlesworth B, Levin BR: Biological and biomedical implications of the co-evolution of pathogens and their hosts. Nat Genet. 2002, 32 (4): 569-577. 10.1038/ng1202-569.PubMedView ArticleGoogle Scholar
- Schmid-Hempel P: Evolutionary Parasitology: The Integrated Study of Infections, Immunology, Ecology, and Genetics. 2011, USA: Oxford Biology. Oxford University PressGoogle Scholar
- Fenner F, Fantini B: Biological Control of Vertebrate pests. The History of Myxomatosis–an Experiment in Evolution. 1999, Oxfordshire: CABI PublishingGoogle Scholar
- Jokela J, Dybdahl MF, Lively CM: The maintenance of sex, clonal dynamics, and host-parasite coevolution in a mixed population of sexual and asexual snails. Am Nat. 2009, 174 (s1): 43-53. 10.1086/599080.View ArticleGoogle Scholar
- Decaestecker E, Gaba S, Raeymaekers JAM, R. Stoks LvK, Ebert D, Meester LD: Host–parasite ‘Red queen’ dynamics archived in pond sediment. Nature. 2007, 450: 870-873. 10.1038/nature06291.PubMedView ArticleGoogle Scholar
- Bohannan BJM, Lenski RE: Linking genetic change to community evolution: insights from studies of bacteria and bacteriophage. Ecol Lett. 2000, 3 (4): 362-377. 10.1046/j.1461-0248.2000.00161.x.View ArticleGoogle Scholar
- Buckling A, Rainey PB: Antagonistic coevolution between a bacterium and a bacteriophage. Proc R Soc B. 2002, 269: 931-936. 10.1098/rspb.2001.1945.PubMed CentralPubMedView ArticleGoogle Scholar
- Brockhurst MA, Morgan AD, Fenton A, Buckling A: Experimental coevolution with bacteria and phage: The pseudomonas fluorescens-phi2 model system. Infect Genet Evol. 2007, 7 (4): 547-552. 10.1016/j.meegid.2007.01.005.PubMedView ArticleGoogle Scholar
- Schulte RD, Makus C, Hasert B, Michiels NK, Schulenburg H: Multiple reciprocal adaptations and rapid genetic change upon experimental coevolution of an animal host and its microbial parasite. Proc Natl Acad Sci USA. 2010, 107 (16): 7359-7364. 10.1073/pnas.1003113107.PubMed CentralPubMedView ArticleGoogle Scholar
- Morran LT, Schmidt OG, Gelarden IA, Parrish-II RC, Lively CM: Running with the red queen: Host-parasite coevolution selects for biparental sex. Science. 2011, 333: 216-218. 10.1126/science.1206360.PubMed CentralPubMedView ArticleGoogle Scholar
- Bérénos C, Schmid-Hempel P, Wegner KM: Experimental coevolution leads to a decrease in parasite-induced host mortality. J Evol Biol. 2011, 24 (8): 1777-1782. 10.1111/j.1420-9101.2011.02306.x.PubMedView ArticleGoogle Scholar
- Fischer O, Schmid-Hempel P: Selection by parasites may increase host recombination frequency. Biol Lett. 2005, 1 (2): 193-195. 10.1098/rsbl.2005.0296.PubMed CentralPubMedView ArticleGoogle Scholar
- King KC, Delph LF, Jokela J, Lively CM: Coevolutionary hotspots and coldspots for host sex and parasite local adaptation in a snail–trematode interaction. Oikos. 2011, 120 (9): 1335-1340. 10.1111/j.1600-0706.2011.19241.x.View ArticleGoogle Scholar
- Koskella B, Thompson JN, Preston GM, Buckling A: Local biotic environment shapes the spatial scale of bacteriophage adaptation to bacteria. Am Nat. 2011, 177 (4): 440-451. 10.1086/658991.PubMedView ArticleGoogle Scholar
- Gandon S, Buckling A, Decaestecker E, Day T: Host-parasite coevolution and patterns of adaptation across time and space. J Evol Biol. 2008, 21: 1861-1866. 10.1111/j.1420-9101.2008.01598.x.PubMedView ArticleGoogle Scholar
- Gaba S, Ebert D: Time-shift experiments as a tool to study antagonistic coevolution. Trends Ecol Evol. 2009, 24 (4): 226-232. 10.1016/j.tree.2008.11.005.PubMedView ArticleGoogle Scholar
- Brockhurst MA, Koskella B: Experimental coevolution of species interactions. Trends Ecol Evol. 2013, 28 (6): 367-375. 10.1016/j.tree.2013.02.009.PubMedView ArticleGoogle Scholar
- Lenski RE, Levin BR: Constraints on the coevolution of bacteria and virulent phage: A model, some experiments and predictions for natural communities. Am Nat. 1985, 125: 585-602. 10.1086/284364.View ArticleGoogle Scholar
- Scanlan PD, Hall AR, Lopez-Pascua LDC, Buckling A: Genetic basis of infectivity evolution in a bacteriophage. Mol Ecol. 2011, 20: 981-989. 10.1111/j.1365-294X.2010.04903.x.PubMedView ArticleGoogle Scholar
- Meyer JR, Dobias DT, Weitz JS, Barrick JE, Quick RT, Lenski RE: Repeatability and contingency in he evolution of a key innovation in phage lambda. Science. 2012, 335: 428-432. 10.1126/science.1214449.PubMed CentralPubMedView ArticleGoogle Scholar
- Hall AR, Scanlan PD, Morgan AD, Buckling A: Host-parasite coevolutionary arms races give way to fluctuating selection. Ecol Lett. 2011, 14: 635-642. 10.1111/j.1461-0248.2011.01624.x.PubMedView ArticleGoogle Scholar
- Lively CM, Dybdahl MF: Parasite adaptation to locally common host genotypes. Nature. 2000, 405: 679-681. 10.1038/35015069.PubMedView ArticleGoogle Scholar
- Anderson RM, May RM: Infectious Diseases of Humans: Dynamics and Control. 1992, USA: Oxford University PressGoogle Scholar
- Dobson AP, Hudson PJ: Microparasites: observed patterns in wild animal populations. Ecology of Infectious Diseases in Natural Populations. Edited by: Grenfell BT, Dobson AP. 1995, Cambridge: Cambridge University Press, 52-89.View ArticleGoogle Scholar
- Hudson PJ, Dobson AP: Macroparasites: observed patterns in naturally fluctuating animal populations. Ecology of Infectious Diseases in Natural Populations. Edited by: Grenfell BT, Dobson AP. 1995, Cambridge: Cambridge University Press, 144-176.View ArticleGoogle Scholar
- Bohannan BJM, Lenski RE: Effect of prey heterogeneity on the response of a model food chain to resource enrichment. Am Nat. 1999, 153: 73-82. 10.1086/303151.View ArticleGoogle Scholar
- Lotka AJ: Analytical note on certain rhythmic relations in organic systems. Proc Nat Acad Sci USA. 1920, 6: 410-415. 10.1073/pnas.6.7.410.PubMed CentralPubMedView ArticleGoogle Scholar
- Volterra V: Variations and fluctuations of the number of individuals in animal species living together. J du conseil international pour l’exploration de la mer. 1928, 3 (1): 3-51.View ArticleGoogle Scholar
- May RM, Anderson RM: Epidemiology and genetics in the coevolution of parasites and hosts. Proc R Soc Lond B Biol Sci. 1983, 219 (1216): 281-313. 10.1098/rspb.1983.0075.PubMedView ArticleGoogle Scholar
- Frank SA: Ecological and genetic models of host-pathogen coevolution. Heredity. 1991, 67: 73-83. 10.1038/hdy.1991.66.PubMedView ArticleGoogle Scholar
- Heesterbeek JAP, Roberts MG: Mathematical models for microparasites of wildlife. Ecology of Infectious Diseases in Natural Populations. Edited by: Grenfell BT, Dobson AP. 1995, Cambridge: Cambridge University Press, 90-122.View ArticleGoogle Scholar
- Roberts MG, Smith G, Grenfell BT: Mathematical models for macroparasites of wildlife. Ecology of Infectious Diseases in Natural Populations. Edited by: Grenfell BT, Dobson AP. 1995, Cambridge: Cambridge University Press, 177-208.View ArticleGoogle Scholar
- Gandon S, Nuismer SL: Interactions between genetic drift, gene flow, and selection mosaics drive parasite local adaptation. Am Nat. 2009, 173 (2): 212-224. 10.1086/593706.PubMedView ArticleGoogle Scholar
- Kirby GC, Burdon JJ: Effects of mutation and random drift on Leonard’s gene-for-gene coevolution model. Phytopathology. 1997, 87 (5): 488-493. 10.1094/PHYTO.19220.127.116.118.PubMedView ArticleGoogle Scholar
- Salathé M, Scherer A, Bonhoeffer S: Neutral drift and polymorphism in gene-for-gene systems. Ecol Lett. 2005, 8: 925-932. 10.1111/j.1461-0248.2005.00794.x.View ArticleGoogle Scholar
- Quigley BJZ, López DG, Buckling A, McKane AJ, Brown SP: The mode of host-parasite interactions shapes coevolutionary dynamics and the fate of host cooperation. Proc R Soc B. 2012, 279 (1743): 3742-3748. 10.1098/rspb.2012.0769.PubMed CentralPubMedView ArticleGoogle Scholar
- Gillespie DT: Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977, 81 (25): 2340-2361. 10.1021/j100540a008.View ArticleGoogle Scholar
- Black AJ, McKane AJ: Stochastic formulation of ecological models and their applications. Trends Ecol Evol. 2012, 27: 337-345. 10.1016/j.tree.2012.01.014.PubMedView ArticleGoogle Scholar
- Stiven AE: Experimental studies on the epidemiology of the host parasite system, hydra and hydramoeba hydroxena (Entz). II. the components of a simple epidemic. Ecol Monogr. 1964, 34 (2): 119-142. 10.2307/1948450.View ArticleGoogle Scholar
- Bach PD, Smith HS: Are population oscillations inherent in the host-parasite relation?. Ecology. 1941, 22 (4): 363-369. 10.2307/1930709.View ArticleGoogle Scholar
- Pimentel D: Population regulation and genetic feedback evolution provides foundation for control of herbivore, parasite, and predator numbers in nature. Science. 1968, 159 (3822): 1432-1437. 10.1126/science.159.3822.1432.PubMedView ArticleGoogle Scholar
- Utida S: Cyclic fluctuations of population density intrinsic to the host-parasite system. Ecology. 1957, 38 (3): 442-449. 10.2307/1929888.View ArticleGoogle Scholar
- Hudson PJ, Dobson AP, Newborn D: Prevention of population cycles by parasite removal. Science. 1998, 282 (5397): 2256-2258.PubMedView ArticleGoogle Scholar
- Gillespie JH: Genetic drift in an infinite population: the pseudohitchhiking model. Genetics. 2000, 155 (2): 909-919.PubMed CentralPubMedGoogle Scholar
- Lenormand T, Roze D, Rousset F: Stochasticity in evolution. Trends Ecol Evol. 2009, 24 (3): 157-165. 10.1016/j.tree.2008.09.014.PubMedView ArticleGoogle Scholar
- Campos PRA, Wahl LM: The effects of population bottlenecks on clonal interference, and the adaptation effective population size. Evolution. 2009, 63 (4): 950-958. 10.1111/j.1558-5646.2008.00595.x.PubMedView ArticleGoogle Scholar
- Wahl LM, Krakauer DC: Models of experimental evolution: the role of genetic chance and selective necessity. Genetics. 2000, 156 (3): 1437-1448.PubMed CentralPubMedGoogle Scholar
- Wahl LM, Gerrish PJ: The probability that beneficial mutations are lost in populations with periodic bottlenecks. Evolution. 2001, 55 (12): 2606-2610.PubMedView ArticleGoogle Scholar
- Handel A, Bennett MR: Surviving the bottleneck: transmission mutants and the evolution of microbial populations. Genetics. 2008, 180 (4): 2193-2200. 10.1534/genetics.108.093013.PubMed CentralPubMedView ArticleGoogle Scholar
- Ponting CP: The functional repertoires of metazoan genomes. Nat Rev Genet. 2008, 9: 689-698. 10.1038/nrg2413.PubMedView ArticleGoogle Scholar
- Levasseur A, Pontarotti P: The role of duplications in the evolution of genomes highlights the need for evolutionary-based approaches in comparative genomics. Biol Direct. 2011, 6 (11): 1-12.Google Scholar
- Andersson DI, Hughes D: Gene amplification and adaptive evolution in bacteria. Ann Rev Genet. 2009, 43: 167-195. 10.1146/annurev-genet-102108-134805.PubMedView ArticleGoogle Scholar
- Moran GP, Coleman DC, Sullivan DJ: Comparative genomics and the evolution of pathogenicity in human pathogenic fungi. Eukaryotic Cell. 2011, 10 (1): 34-42. 10.1128/EC.00242-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Merhej V, Raoult D: Rickettsial evolution in the light of comparative genomics. Biol Rev. 2011, 86: 379-405. 10.1111/j.1469-185X.2010.00151.x.PubMedView ArticleGoogle Scholar
- Iskow RC, Gokcumen O, Lee C: Exploring the role of copy number variants in human adaptation. Trends Genet. 2012, 28 (6): 245-257. 10.1016/j.tig.2012.03.002.PubMed CentralPubMedView ArticleGoogle Scholar
- Korbel JO, Kim PM, Chen X, Urban AE, Weissman S, Snyder M, Gerstein MB: The current excitement about copy-number variation: how it relates to gene duplications and protein families. Curr Opin Struct Biol. 2008, 18: 366-374. 10.1016/j.sbi.2008.02.005.PubMed CentralPubMedView ArticleGoogle Scholar
- Lively CM, Apanius V: Genetic diversity in host-parasite interactions. Ecology of Infectious Diseases in Natural Populations. Edited by: Grenfell BT, Dobson AP. 1995, Cambridge: Cambridge University Press, 421-449.View ArticleGoogle Scholar
- Frank SA: Specific and non-specific defense against parasitic attack. J Theor Biol. 2000, 202 (4): 283-304. 10.1006/jtbi.1999.1054.PubMedView ArticleGoogle Scholar
- Agrawal AF, Lively CM: Modelling infection as a two-step process combining gene-for-gene and matching-allele genetics. Proc R Soc B: Biol Sci. 2003, 270 (1512): 323-334. 10.1098/rspb.2002.2193.View ArticleGoogle Scholar
- Fenton A, Antonovics J, Brockhurst MA: Two-step infection processes can lead to coevolution between functionally independent infection and resistance pathways. Evolution. 2012, 66 (7): 2030-2041. 10.1111/j.1558-5646.2012.01578.x.PubMedView ArticleGoogle Scholar
- Thrall PH, Laine A-L, Ravensdale M, Nemri A, Dodds PN, Barrett LG, Burdon JJ: Rapid genetic change underpins antagonistic coevolution in a natural host-pathogen metapopulation. Ecol Lett. 2012, 15 (5): 425-435. 10.1111/j.1461-0248.2012.01749.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Wolinska J, King KC: Environment can alter selection in host-parasite interactions. Trends Parasitol. 2009, 25 (5): 236-244. 10.1016/j.pt.2009.02.004.PubMedView ArticleGoogle Scholar
- Dybdahl MF, Lively CM: Host-parasite coevolution: evidence for rare advantage and time-lagged selection in a natural population. Evolution. 1998, 52 (4): 1057-1066. 10.2307/2411236.View ArticleGoogle Scholar
- Eizaguirre C, Lenz TL, Kalbe M, Milinski M: Rapid and adaptive evolution of MHC genes under parasite selection in experimental vertebrate populations. Nat Commun. 2012, 3 (621): 1-6.Google Scholar
- Luijckx P, Fienberg H, Duneau D, Ebert D: A matching-allele model explains host resistance to parasites. Curr Biol. 2013, 23 (12): 1085-1088. 10.1016/j.cub.2013.04.064.PubMedView ArticleGoogle Scholar
- Gillespie D: A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. J Comput Phys. 1976, 22: 403-434. 10.1016/0021-9991(76)90041-3.View ArticleGoogle Scholar
- McKane AJ, Newman TJ: Stochastic models in population biology and their deterministic analogs. Phys Rev E. 2004, 70: 19-Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License(http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.