PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(8): e12440.
Published online 2010 August 27. doi:  10.1371/journal.pone.0012440
PMCID: PMC2929189

Parasite Replication and the Evolutionary Epidemiology of Parasite Virulence

Mike B. Gravenor, Editor

Abstract

Parasite virulence evolution is shaped by both within-host and population-level processes yet the link between these differing scales of infection is often neglected. Population structure and heterogeneity in both parasites and hosts will affect how hosts are exploited by pathogens and the intensity of infection. Here, it is shown how the degree of relatedness among parasites together with epidemiological parameters such as pathogen yield and longevity influence the evolution of virulence. Furthermore, the role of kin competition and the degree of cheating within highly structured parasite populations also influences parasite fitness and infectivity patterns. Understanding how the effects of within-host processes scale up to affect the epidemiology has importance for understanding host-pathogen interactions.

Introduction

Virulence, the pathogen-induced reduction in host fitness (through either increased host mortality or reduced host growth rate) is an emergent epidemiological property of a parasite through its interaction with its host. Alternatively, virulence is defined as the ability of a parasite to overcome a given level of host resistance. This concept has given rise to co-evolutionary genetic (e.g. gene-for-gene) mechanisms for the evolution of virulence (particularly in plant-parasite interactions [1]). However, the traditional view that evolutionary change in parasite virulence leads to a reduction in the capacity to induce harm to hosts and prolong host longevity favouring parasite fitness has been challenged on many occasions and in many different studies [2][7]. Parasite virulence is a complex product of many interactions within a host and can depend on a range of different factors that are responsible for the expression of the trait. These include the sublethal expression of protein toxins [8], [9] that affect the life history characteristics of the host [10], the life history of the parasite [6], [11][12], the multiplicity of infections [13], [14] and temporal and/or spatial environmental heterogeneity [15][18].

Although parasite virulence is explicitly integrated to within-host processes, the evolution of this parasite trait may not necessarily be linked to other parasite life history strategies operating at different scales such as pathogen transmission between hosts [19], [20]. Virulence can act in a local, narrow way and lead to high-levels of morbidity (or mortality) within an individual host such that optimal transmission at an epidemiological level might be compromised [21]. The structure of the microenvironment of individual hosts [22], interactions amongst parasites through competition and cooperation [23][25] and the effects of the host immune system [7], [26][27] can all contribute to the differential evolution of virulence amongst hosts. Broadly, heterogeneities will affect virulence evolution [16] as patterns of parasite replication alter host exploitation and affect the intensity of infection.

This heterogeneity and patterns of host exploitation will affect the genetic structure of pathogen populations. Under low levels of pathogen recombination [28], parasite replication within hosts allows population structure to form [29]. This population-level viscosity is known to affect the degree of relatedness and the ecology of local interactions [30][34]. While viscosity affects inclusive fitness by increasing the degree of relatedness [30][31] local competition between kin can decrease any benefits to social interactions [32][34]. More recently, and in the context of host-pathogen interactions, the role of local interactions and competition has been shown to affect patterns of pathogen infectivity [18] and epidemiology (see [35] for a review).

While the epidemiology of host-pathogen interactions is well established [36], [37] the effect of coupling the within-host and epidemiological consequences of parasitic infections is a relatively novel advance in parasite evolutionary ecology [35]. Here, the goal of this study is to address this issue and explore how patterns of parasite replication affect the broader epidemiological consequences of host-pathogen interactions. In explicitly considering this link it becomes essential to separate out the effects of competition from the effects of replication on the evolution of parasite life history traits. For instance, the effects of kin competition and parasite replication can be considered as two mechanisms by which cooperation and collective actions might be manifest. Reducing kin competition can be viewed as a cooperative behaviour and replication within a host allows successful resource exploitation (often through the production of shared public goods).

As outlined, the aim of this study is to explore, theoretically, how within-host pathogen replication affects parasite fitness and epidemiology. After introducing a description of pathogen fitness involving replicating and non-replicating strategies, the optimal replication strategy together with the effects of relatedness on parasite fitness is derived. Following on from this, the epidemiological consequences of within-host replication are explored in a host-pathogen model in which the pathogen is an obligate killer and has free-living infectious stages. These results are discussed with reference to recent developments in host-pathogen epidemiology and virulence evolution.

Within-Host Dynamics

Parasite Replication

Parasite fitness is assumed to be dependent on the parasite's own replication strategy and on the effects of the average replication strategy across the group of parasites within a host. In more detail, the fitness function can be derived in terms of whether a parasite adopts to undergo costly replication or is a strategy that chooses not to replicate to kill the host but utilizes resources from the dead cadaver. The relative success of a parasite that replicates is a function of the costs of replication, the strength of competition from these related (kin) parasites and the proportion of non-kin in the host. This biology is described by:

equation image
(1)

where ρih is the replication fraction of the ith parasite in the hth host, An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e002.jpg is the average replication fraction across the group of parasites in a host, c is the costs associated with replicating, α is the strength of competition. Thus An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e003.jpg is the fitness benefit to parasites that replicate relative to the strength of competition amongst related parasites (An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e004.jpg) and the strength of interactions with replicating non-kin (An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e005.jpg).

Importantly, the strength of competition (α) amongst the focal parasite strain is the key feature of the fitness function and determines the relative strength of the interaction amongst kin and consequently pathogen growth and abundance within a host. If α>1 kin compete more strongly amongst themselves, whereas if α<1 kin compete less intensely amongst themselves and, therefore, may act cooperatively.

The impact of non-replicating parasites An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e006.jpg on fitness can be determined from the proportion of non-replicators relative to the strength of competition amongst non-replicators and the presence of non-kin. This biology is captured with:

equation image
(2)

where An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e008.jpg is the strength of non-replicating parasites (0<An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e009.jpg<1) relative to the replicating strategy. Non replicating parasites do not any pay costs associated with growth but may affect overall host morbidity and mortality (see below). Given all of this, the overall fitness (An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e010.jpg) of the ith parasite phenotype in the hth host can be defined as:

equation image
(3)

Parasite fitness is a non-linear function of replication strategy (Figure 1a). Low replication strategies have low fitness as parasite growth rate is restricted. Similarly high levels of replication have low fitness as the costs of replication act to restrict fitness. Increases in the strength of kin competition (α) also act to decrease fitness (Figure 1a).

Figure 1
The role of parasite replication and relatedness on fitness.

The neighbour-modulated (inclusive fitness) effects of replicating and non-replicating parasites [38] can be determined by considering how an alternative parasite strategy (with a different replication phenotype: An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e013.jpg) alters fitness [39]. In doing so it becomes necessary to evaluate how a parasite alters both direct and indirect fitness under weak selection. This is done through a fitness maximisation method using the multivariate chain rule [39], [40]:

equation image
(4)

where the ratio An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e015.jpg is a statistical measure of relatedness [41] (and is the change in the average phenotype with respect to changes in the focal phenotype – [39]).

The optimal replication fraction is found by setting equation 4 equal to zero and solving for An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e016.jpg, which in the absence of the non-replicator strategy (η = 0), yields

equation image
(5)

where R is the degree of relatedness. The optimal replication strategy declines with increasing relatedness. The presence of non-replicators affects the optimal parasite strategy and as the strength of non-replicators increases (η = 1), the optimal replication fraction (under minimal replication costs, c→0) is then the positive solution from:

equation image
(6)

Under different levels of kin competition the optimal replication fraction declines as relatedness increases (Figure 1b). In unrelated environments (R→0), the optimal replication strategy is high leading to high host exploitation. When parasites are related (R→1), the optimal replication strategy is low favouring prudent host exploitation (Figure 1b). Reduced kin competition (α<1) gives rise to lower optimal strategies as cooperative actions are more prevalent and leads to the prudent exploitation of the host by the parasite. Increasing the effects of the non-replicator effects (increasing η) reduces the prudent exploitation strategy and facilitates higher optimal replication strategies leading to high host exploitation.

Relatedness

Relatedness is not a fixed quantity but a statistical relationship dependent on parasite density [41]. As such, the degree of relatedness (R) amongst (haploid) parasites within a host can be determined from standard population genetic theory [42]:

equation image
(7)

where V0 is the effective parasite population size and x is the mutation rate. At equilibrium (R*),

equation image
(8)

Together with equation (6), this expression (equation 8) allows the optimal parasite replication strategy to be expressed in terms of the degree of relatedness and the effective parasite population size (within a host). Effective population size is driven by the strength of competition amongst parasites within a host (α) and hence understanding the within- and between-host epidemiological processes is critical in order to determine an appropriate measure of pathogen fitness.

Host-Pathogen Dynamics

In order to link the within-host replication dynamics to the epidemiology of host-pathogen interaction a host-pathogen epidemiological framework (e.g., [36]) is developed in which the pathogen (V) is lethal to the host (H) and has free-living infectious stages. Following infection, infected hosts (I) die and cadavers yield a number of free-living pathogen particles. Pathogen replication (ρ) directly affects pathogen virulence (m(ρ)), and indirectly affects pathogen yield (through virulence) (f(m)). The epidemiological model is of the form:

equation image
(9–11)

where a is the host birth rate, ρ is the pathogen transmission rate, dH is the host death rate (independent of both host and parasite density) and dV is the death rate of the free-living pathogen. Virulence, the pathogen induced mortality rate (An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e022.jpg) is assumed to be a decelerating (convex) function of within-host replication fraction such that An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e023.jpg and pathogen yield (f(m)) following death of the host is an increasing linear function of pathogen replication rate such that An external file that holds a picture, illustration, etc.
Object name is pone.0012440.e024.jpg, where ε is a positive scaling constant linking virulence and yield (yield strength).

As outlined, virulence increases as the pathogen replication rate increases but saturates at high level of replication. However, the presence of non-replicating pathogens has a further increased effect on virulence as the optimal replication fraction increases (Figure 2). Given this, it is important to separate the effects of non-replication from cheating and kin competition. Non-replicating parasites may still contribute to pathogen-induced morbidity or mortality within a host (and consequently affect virulence) even if the effects of kin competition may be severe and act to limit overall parasite fitness.

Figure 2
The role of non-replicating parasites on virulence (m(ρ)).

The consequences of non-replicating pathogens on the host-pathogen interaction can be explored further by considering the conditions that would allow a rare novel variant of the disease (In, Vn) to invade and spread. This occurs if the net population growth rate of the novel variant is greater than zero (that is dIn/dt>0 and dVn/dt>0). This invasion analysis approach (which is similar to a standard local stability analysis where the resident strategy is at equilibrium) is an appropriate measure of fitness [43], [44] when density-dependent processes operate. Fitness is evaluated by taking the determinant of:

equation image
(12)

and solving for the dominant eigenvalue, λ (where In and Vn are the infected hosts and free-living pathogen stages associated with the novel invading pathogen, respectively and H* and I* are the equilibrial abundances of susceptible hosts and hosts infected with the ancestral pathogen, respectively).

In the absence of non-replicators, the invasion consequences of novel pathogens can be shown to depend on three main parasite life-history traits (Figure 3): free-living pathogen death rate (dV), yield strength (ε) and the mutation rate (a factor governing the degree of relatedness within a host – equation 8). Long-lived free-living pathogen stages (low dV) favour the invasion of novel disease (Figure 3a) particularly when parasites within a host are of intermediate relatedness or unrelated (~0.6<x<1.0). Similarly, pathogens of intermediate to high relatedness together with high yields further increase the likelihood that a novel disease strategy can invade (Figure 3b).

Figure 3
Fitness contours for rare disease variants in the absence of non-replicators.

In the presence of non-replicators (Figure 4), the previous patterns described hold: high yields and low death rates lead to pathogen invasion and spread. The presence of non-replicating parasites leads to a maximum in fitness when parasites are unrelated (x→1.0) (Figure 4). Alterations in the strength of kin competition (α) can also affect disease invasion potential (Figure 5). Increases in competition can offset the benefits to high yield and favour the invasion of pathogens that produce lower numbers of infectious particles. Furthermore, the degree of competition can also affect the optimal levels of relatedness within hosts. Strong competition within hosts leads to optimal yields at intermediate levels of relatedness. In contrast, reduced competition (cooperation) favours high yields at high levels of relatedness.

Figure 4
Fitness contours for rare disease variants in the presence of non-replicators.
Figure 5
Fitness contours (Yield strength (ε) – Mutation probability (x)) for rare disease variants in the presence of non-replicators under varying degrees of kin competition.

Discussion

Here, the effects of how pathogen dynamics within-host scale up to affect dynamics between hosts have been explored. By considering the effects of parasite relatedness, it has been shown that the epidemiological outcome can depend on the degree of similarity between parasites in a host. In particular, increases in parasite relatedness led to a decline in the optimal replication fraction as the interaction between closely related individuals intensifies. Increases in competition can alleviate these interactions and contribute to parasite fitness.

The consequences of these within-host processes have important epidemiological consequences [5], [45] particularly in terms of parasite life history traits such as virulence. Several studies have highlighted that the effects of virulence can not be understood in isolation but must be considered within the context of host and parasite life histories [5], [10][11]. For instance, following death from infection, hosts release a yield of free-living pathogen, however this yield is constrained by the pathogen's speed of kill. Highly virulent parasites that kill hosts quickly yield fewer pathogens whereas more benign parasites may yield a higher number of transmissible agents. Selection will act differently depending on the interaction among factors such as the mode of transmission, speed of kill (virulence), free-living pathogen persistence time and within-host rate of replication.

Coupled with this is the fact that parasites within a host will be related (through the almost clonal expansion and proliferation on host resources) particularly if levels of recombination are low [28]. This will lead to competition between kin and non-kin and has consequences for virulence evolution and host-pathogen epidemiology. Here, it has been shown that increasing kin competition affects the invasion potential of pathogens such that maximum fitness is achieved under high pathogen yields and differing levels of relatedness. The presence of non-replicating parasites can create conditions where within-host levels of parasite relatedness are low even when kin competition is relatively weak (α<1). Increasing kin competition leads to an optimal invasion strategy similar to that observed when non-replicating parasites are absent.

Cooperation (and the potential for cheating) amongst parasites can take many forms – for instance it could involve participation in replication (and the production of shared toxins) to exploit hosts [46], [47], interference between non-kin through competition or collective action to avoid host defences. Cooperation between parasites increases the utilization of the host resource and can lead to a positive relationship between relatedness and virulence [23], [48]. Multiple infections are expected to reduce relatedness and have consequences for pathogen reproduction and virulence [47], [49]. However, the important issue here is not the relationship between virulence and relatedness but how individual parasites perform in individual hosts and how this translates to affect parasite life-history strategies. The presence of non-replicating (but otherwise equal) pathogens may still have positive fitness contributions to host exploitation and might act (sub)additively to increase virulence. The presence of non-replicators (and potentially non-kin) will alter the competitive environment and appropriate consideration of this alteration in the strength of competition is actually likely to be system-specific. Competition for limited resources even amongst closely related individuals will lead to winners and losers and within hosts this is likely to favour not only parasites that act to rapidly exploit hosts (and thereby the more virulent strains) but also virulence polymorphisms [14]. Strong competition increases parasite virulence as selection acts to promote more aggressive (resource-capturing) genotypes.

Virulence polymorphisms may arise through other manifestation of competition and host exploitation. Competition and exploitation of resources by pathogens will be different amongst hosts (due to host heterogeneity) and may lead to differential population and genetic structure between parasites [50]. As such, different levels of relatedness within hosts will lead to emergent phenomena such as virulence being localised and myopic [19]. Changes in the degree of relatedness will affect the optimal replication fraction and influence the observed level of virulence. As parasite virulence may be influenced by host and parasite life history characteristics [10][12] this also has important implications for host exploitation strategies [51]. High yielding and/or long-lived free-living infectious stages influence the evolutionary optimal strategy. This occurs through inclusive fitness benefits. Such parasite life histories have both direct effects on fitness (virulence-life history trait correlation) and indirect effects through changes in the levels of relatedness within hosts. Levels of selection may operate differently within and between hosts [52] to favour extended bouts of successful infections amongst hosts in a population [4].

Understanding the epidemiology of host-pathogen dynamics necessitates an appreciation of both the within-host [26], [27] and between host [37] dynamics. Understanding how selection operates on these various aspects of parasite fitness will reveal how parasites interact to affect epidemiological patterns [35]. Appreciating the finer implications of population and genetic structures will have important consequences not only for understanding pathogen evolutionary ecology but also for developing public health intervention programmes.

Acknowledgments

I thank Bob May for comments on this work.

Footnotes

Competing Interests: The author has declared that no competing interests exist.

Funding: The work was supported by the Royal Society and the UK Biotechnology and Biological Sciences Research Council (BB/C005172/1). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Thrall PH, Burdon JJ. Evolution of virulence in a plant host-pathogen metapopulaiton. Science. 2003;299:1735–1737. doi: 10.1126/science.1080070. [PubMed]
2. Levin S, Pimentel D. Selection of intermediate rates of increase in parasite-host systems. Am Nat. 1981;117:308–315. doi: 10.1086/283708.
3. Anderson RM, May RM. Coevolution of hosts and parasites. Parasitology. 1982;85:411–426. doi: 10.1017/S0031182000055360. [PubMed]
4. Nesse RM, Williams GC. Why we get sick. The new science of Darwinian medicine. New York: Vintage Books; 1996.
5. Frank SA. Models of parasite virulence. Q Rev Biol. 1996;71:37–78. doi: 10.1086/419267. [PubMed]
6. Day T. Parasite transmission modes and the evolution of virulence. Evolution. 2001;55:2389–2400. doi 10.1111/j.0014-3820.2001.tb00754.x. [PubMed]
7. Andre JB, Godelle B. Within-host evolution and virulence in microparasites. J Theor Biol. 2006;241:402–409. doi: 10.1016/j.jtbi.2005.12.011. [PubMed]
8. Porcelli SA. Tuberculosis: shrewd survival strategy. Nature. 2008;454:702–703. doi: 10.1038/454702a. [PubMed]
9. Raghavan S, Manzanillo P, Chan K, Dovey C, Cox JS. Secreted transcription factor controls Mycobacterium tuberculosis virulence. Nature. 2008;454:717–721. doi: 10.1038/nature07219. [PMC free article] [PubMed]
10. Gandon S, Jansen VAA, van Baalen M. Host life history and the evolution of parasite virulence. Evolution. 2001;55:1056–1062. doi: 10.1111/j.0014-3820.2001.tb00622.x. [PubMed]
11. Bonhoeffer S, Lenski RE, Ebert D. The curse of the pharaoh: the evolution of virulence in pathogens with long living propagules. Proc R Soc Lond B. 1996;263:715–721. doi: 10.1098/rspb.1996.0107. [PubMed]
12. Gandon S. The curse of the pharaoh hypothesis. Proc R Soc Lond B. 1998;265:1545–1552. doi: 10.1098/rspb.1998.0470. [PMC free article] [PubMed]
13. Nowak MA, May RM. Superinfection and the evolution of parasite virulence. Proc R Soc Lond B. 1994;255:81–89. doi: 10.1098/rspb.1994.0012. [PubMed]
14. van Baalen M, Sabelis MW. The dynamics of multiple infection and the evolution of virulence. Am Nat. 1995;146:881–910. doi: 10.1086/285830.
15. van Baalen M, Rand DA. The unit of selection in viscous populations and the evolution of altruism. J Theor Biol. 1998;193:631–648. doi: 10.1006/jtbi.1998.0730. [PubMed]
16. Ganusov VV, Bergstrom CT, Antia R. Within-host population dynamics and the evolution of microparasites in a heterogeneous host population. Evolution. 2002;56:213–223. doi: 10.1111/j.0014-3820.2002.tb01332.x. [PubMed]
17. Kamo M, Boots M. The curse of the pharaoh in space: free-living infectious stages and the evolution of virulence in spatially explicit populations. J Theor Biol. 2004;231:435–441. doi: 10.1016/j.jtbi.2004.07.005. [PubMed]
18. Boots M, Mealor M. Local interactions select for lower pathogen infectivity. Science. 2007;315:1284–1286. doi: 10.1126/science.1137126. [PubMed]
19. Levin BR, Bull JJ. Short-sighted evolution and the virulence of pathogenic microorganisms. Trends Microbiol. 1994;2:76–81. doi: 10.1016/0966-842X(94)90538-X. [PubMed]
20. Ebert D, Bull JJ. Challenging the trade-off model for the evolution of virulence: is virulence management feasible? Trends. Microbiol. 2003;11:15–20. doi 10.1016/S0966-842X(02)00003-3. [PubMed]
21. Levin BR, Lipsitch M, Bonhoeffer S. Evolution and disease - population biology, evolution, and infectious disease: convergence and synthesis. Science. 1999;283:806–809. doi: 10.1126/science.283.5403.806. [PubMed]
22. Raymond B, Lijek RS, Griffiths RI, Bonsall MB. Ecological consequences of ingestion of Bacillus cereus on Bacillus thuringiensis infections and on the gut flora of a lepidopteran host. J Invert Pathol. 2008;99:103–111. doi: 10.1016/j.jip.2008.04.007. [PubMed]
23. West SA, Buckling A. Cooperation, virulence and siderophore production in bacterial parasites. Proc R Soc Lond B. 2003;270:37–44. doi: 10.1098/rspb.2002.2209. [PMC free article] [PubMed]
24. Griffin AS, West SA, Buckling A. Cooperation and competition in pathogenic bacteria. Nature. 2004;430:1024–1027. doi: 10.1038/nature02744. [PubMed]
25. Harrison F, Paul J, Massey RC, Buckling A. Interspecific competition and siderophore-mediated competition in Pseudomonas aeruginosa. ISME Journal. 2008;2:49–55. doi: 10.1038/ismej.2007.96. [PubMed]
26. Antia R, Levin BR, May RM. Within-host population dynamics and the evolution and maintenance of microparasite virulence. Am Nat. 1994;144:457–472. doi: 10.1086/285686.
27. Fenton A, Lello J, Bonsall MB. Pathogen responses to host immunity: the impact of time delays and memory on the evolution of virulence. Proc R Soc Lond B. 2006;273:2083–2090. doi: 10.1098/rspb.2006.3552. [PMC free article] [PubMed]
28. Spratt BG, Maiden MCJ. Bacterial population genetics, evolution and epidemiology. Phil Trans R Soc Lond B. 1999;354:701–710. doi: 10.1098/rstb.1999.0423. [PMC free article] [PubMed]
29. Feil EJ, Spratt BG. Recombination and the population structure of bacterial pathogens. Ann Rev Microbiol. 2001;55:561–590. doi: 10.1146/annurev.micro.55.1.561. [PubMed]
30. Hamilton WD. Genetical evolution of social behaviour. I. J Theor Biol. 1964a;7:1–16. doi: 10.1016/0022-5193(64)90038-4. [PubMed]
31. Hamilton WD. Genetical evolution of social behaviour. II. J Theor Biol. 1964a;7:17–54. doi: 10.1016/0022-5193(64)90039-6. [PubMed]
32. Taylor PD. Inclusive fitness in a homogeneous environment. Proc R Soc Lond B. 1992a;249:299–302. doi: 10.1098/rspb.1992.0118.
33. Taylor PD. Altruism in viscous populations – an inclusive fitness model. Evol Ecol. 1992b;6:352–356. doi: 10.1007/BF02270971.
34. Queller DC. Genetic relatedness in viscous populations. Evol Ecol. 1994;8:70–73. doi: 10.1007/BF01237667.
35. Mideo N, Alizon S, Day T. Linking within- and between-host dynamics in the evolutionary epidemiology of infectious disease. Trends Ecol Evol. 2008;23:511–517. doi: 10.1016/j.tree.2008.05.009. [PubMed]
36. Anderson RM, May RM. The population dynamics of microparasites and their invertebrate hosts. Phil Trans R Soc Lond B. 1981;291:451–524. doi: 10.1098/rstb.1981.0005.
37. Anderson RM, May RM. Infectious diseases of humans. Dynamics and control. Oxford: Oxford University Press; 1991.
38. Taylor PD, Wild G, Gardner A. Direct fitness or inclusive fitness: how shall we model kin selection? J Evol Biol. 2007;20:301–309. doi: 10.1111/j.1420-9101.2006.01196.x. [PubMed]
39. Taylor PD, Frank SA. How to make a kin selection model. J Theor Biol. 1996;180:27–37. doi: 10.1006/jtbi.1996.0075. [PubMed]
40. Binmore K, Davis J. Calculus: concepts and methods. Cambridge: Cambridge University Press; 2002.
41. Michod RE, Hamilton WD. Coefficients of relatedness in sociobiology. Nature. 1980;288:694–697. doi: 10.1038/288694a0.
42. Crow JF, Kimura M. An introduction to population genetic theory. New York: Harper & Row Ltd; 1970.
43. Metz JAJ, Nisbet RM, Geritz SAH. How should we define fitness for general ecological scenarios. Trends Ecol Evol. 1992;7:198–202. doi: 10.1016/0169-5347(92)90073-K. [PubMed]
44. Pasztor L, Meszena G, Kisdi E. R(0) or r: A matter of taste? J Evol Biol. 1996;9:511–516. (doi: 10.1046/j.1420-9101.1996.9040511.x)
45. Alizon S, van Baalen M. Acute or chronic? Within-host models with immune dynamics, infection outcome, and parasite evolution. Am Nat. 2008;172:E244–E256. doi: 10.1086/592404. [PubMed]
46. Brown SP, Hochberg ME, Grenfell BT. Does multiple infection select for raised virulence? Trends Microbiol. 2002;10:401–405. doi: 10.1016/S0966-842X(02)02413-7. [PubMed]
47. Raymond B, Davis D, Bonsall MB. Competition and reproduction in mixed infections of pathogenic and non-pathogenic Bacillus spp. J Invert Pathol. 2007;96:151–155. doi: 10.1016/j.jip.2007.03.003. [PubMed]
48. Chao L, Hanley KA, Burch CL, Dahlberg C, Turner PE. Kin selection and parasite evolution: Higher and lower virulence with hard and soft selection. Q Rev Biol. 2000;75:261–275. doi: 10.1086/393499. [PubMed]
49. Massey RC, Buckling A, Ffrench-Constant R. Interference competition and parasite virulence. Proc R Soc Lond B. 2004;271:785–788. doi: 10.1098/rspb.2004.2676. [PMC free article] [PubMed]
50. Roussett F. Genetic structure and selection in subdivided populations. Princeton: Princeton University Press; 2004.
51. Walther BA, Ewald PW. Pathogen survival in the external environment and the evolution of virulence. Biol Rev. 2004;79:849–869. doi: 10.1017/S1464793104006475. [PubMed]
52. Bremermann HJ, Pickering J. A game-theoretical model of parasite virulence. J Theor Biol. 1983;100:411–426. doi: 10.1016/0022-5193(83)90438-1. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science