Within-host population dynamics and the evolution of microparasites in a heterogeneous host population

Abstract Why do parasites harm their hosts? The general understanding is that if the transmission rate and virulence of a parasite are linked, then the parasite must harm its host to maximize its transmission. The exact nature of such trade‐offs remains largely unclear, but for vertebrate hosts it probably involves interactions between a microparasite and the host immune system. Previous results have suggested that in a homogeneous host population in the absence of super‐ or coinfection, within‐host dynamics lead to selection of the parasite with an intermediate growth rate that is just being controlled by the immune system before it kills the host (Antia et al. 1994). In this paper, we examine how this result changes when heterogeneity is introduced to the host population. We incorporate the simplest form of heterogeneity–random heterogeneity in the parameters describing the size of the initial parasite inoculum, the immune response of the host, and the lethal density at which the parasite kills the host. We find that the general conclusion of the previous model holds: parasites evolve some intermediate growth rate. However, in contrast with the generally accepted view, we find that virulence (measured by the case mortality or the rate of parasite‐induced host mortality) increases with heterogeneity. Finally, we link the within‐host and between‐host dynamics of parasites. We show how the parameters for epidemiological spread of the disease can be estimated from the within‐host dynamics, and in doing so examine the way in which trade‐offs between these epidemiological parameters arise as a consequence of the interaction of the parasite and the immune response of the host.     

Optimal killing for obligate killers: the evolution of life histories and virulence of semelparous parasites.

Many viral, bacterial and protozoan parasites of invertebrates first propagate inside their host without releasing any transmission stages and then kill their host to release all transmission stages at once. Life history and the evolution of virulence of these obligately killing parasites are modelled, assuming that within-host growth is density dependent. We find that the parasite should kill the host when its per capita growth rate falls to the level of the host mortality rate. The parasite should kill its host later when the carrying capacity, K, is higher, but should kill it earlier when the parasite-independent host mortality increases or when the parasite has a higher birth rate. When K(t), for parasite growth, is not constant over the duration of an infection, but increases with time, the parasite should kill the host around the stage when the growth rate of the carrying capacity decelerates strongly. In case that K(t) relates to host body size, this deceleration in growth is around host maturation.     

A model for the coevolution of immunity and immune evasion in vector-borne diseases with implications for the epidemiology of malaria

We describe a model of host-parasite coevolution, where the interaction depends on the investments by the host in its immune response and by the parasite in its ability to suppress (or evade) its host's immune response. We base our model on the interaction between malaria parasites and their mosquito hosts and thus describe the epidemiological dynamics with the Macdonald-Ross equation of malaria epidemiology. The qualitative predictions of the model are most sensitive to the cost of the immune response and to the intensity of transmission. If transmission is weak or the cost of immunity is low, the system evolves to a coevolutionarily stable equilibrium at intermediate levels of investment (and, generally, at a low frequency of resistance). At a higher cost of immunity and as transmission intensifies, the system is not evolutionarily stable but rather cycles around intermediate levels of investment. At more intense transmission, neither host nor parasite invests any resources in dominating its partner so that no resistance is observed in the population. These results may help to explain the lack of encapsulated malaria parasites generally observed in natural populations of mosquito vectors, despite strong selection pressure for resistance in areas of very intense transmission.     

Reduction of transmission stages concomitant with increased host immune responses to hypervirulent Sarcocystis singaporensis, and natural selection for intermediate virulence.

Parasite virulence (pathogenicity depending on inoculum size) and host immune reactions were examined for the apicomplexan protozoan Sarcocystis singaporensis. This parasite is endemic in southeastern Asia and multiplies as a proliferation (merozoite) and transmission stage (bradyzoite) in rats. Virulence in wild brown rats of parasites freshly isolated in the wild (wild-type) was surprisingly constant within the endemic area and showed an intermediate level. In contrast, serially passaged parasites either became avirulent or virulence increased markedly (hypervirulence). Production of transmission stages was maximal for the wild-type whereas numbers were significantly reduced for hypervirulent and avirulent (shown in a previous study) parasites. Analyses of B and T cell immunity revealed that immune responses of WKY rats to the transmission stage were significantly higher for hypervirulent than for wild-type parasites. These results suggest that it is the immune system of the host that is not only responsible for reduction of transmission stages in individual rats, but also could act as a selective force that maintains intermediate virulence at the population level because reduction of muscle stages challenges transmission of S. singaporensis to the definitive host. Collectively, the presented data support evolutionary theory, which predicts intermediate rates of parasite growth in nature and an ‘arms race’ between host immunity and parasite proliferation.     

Who is in control of the stickleback immune system: interactions between Schistocephalus solidus and its specific vertebrate host

The cestode Schistocephalus solidus is a frequent parasite of three-spined sticklebacks and has a large impact on its host's fitness. Selection pressure should therefore be high on stickleback defence mechanisms, like an efficient immune system, and also on parasite strategies to overcome these. Even though there are indications for manipulation of the immune system of its specific second intermediate host by the cestode, nothing is yet known about the chronology of specific interactions of S. solidus with the stickleback immune system. We here expected sticklebacks to first mount an innate immune response directly post-exposure to the parasite to clear the infection at an early stage and after an initial lag phase to upregulate adaptive immunity. Most interestingly, we did not find any upregulation of the specific lymphocyte-mediated immune response. Also, the pattern of activation of the innate immune system did not match our expectations: the proliferation of monocytes followed fluctuating kinetics suggesting that the parasite repeatedly installs a new surface coat not immunogenic to the host. Furthermore, the respiratory burst activity, which has the potential to clear an early S. solidus infection, was upregulated very late during infection, when the parasite was too big to be cleared but ready for transmission to its final host. We here suggest that the late activation of the innate immune system interferes with the neuroendocrine system, which mediates reduced predation avoidance behaviour and so facilitates the transmission to the final host.     

Reduction of transmission stages concomitant with increased host immune responses to hypervirulent Sarcocystis singaporensis, and natural selection for intermediate virulence

Parasite virulence (pathogenicity depending on inoculum size) and host immune reactions were examined for the apicomplexan protozoan Sarcocystis singaporensis. This parasite is endemic in southeastern Asia and multiplies as a proliferation (merozoite) and transmission stage (bradyzoite) in rats. Virulence in wild brown rats of parasites freshly isolated in the wild (wild-type) was surprisingly constant within the endemic area and showed an intermediate level. In contrast, serially passaged parasites either became avirulent or virulence increased markedly (hypervirulence). Production of transmission stages was maximal for the wild-type whereas numbers were significantly reduced for hypervirulent and avirulent (shown in a previous study) parasites. Analyses of B and T cell immunity revealed that immune responses of WKY rats to the transmission stage were significantly higher for hypervirulent than for wild-type parasites. These results suggest that it is the immune system of the host that is not only responsible for reduction of transmission stages in individual rats, but also could act as a selective force that maintains intermediate virulence at the population level because reduction of muscle stages challenges transmission of S. singaporensis to the definitive host. Collectively, the presented data support evolutionary theory, which predicts intermediate rates of parasite growth in nature and an ‘arms race’ between host immunity and parasite proliferation.     

Predicting optimal transmission investment in malaria parasites

In vertebrate hosts, malaria parasites face a tradeoff between replicating and the production of transmission stages that can be passed onto mosquitoes. This tradeoff is analogous to growth‐reproduction tradeoffs in multicellular organisms. We use a mathematical model tailored to the life cycle and dynamics of malaria parasites to identify allocation strategies that maximize cumulative transmission potential to mosquitoes. We show that plastic strategies can substantially outperform fixed allocation because parasites can achieve greater fitness by investing in proliferation early and delaying the production of transmission stages. Parasites should further benefit from restraining transmission investment later in infection, because such a strategy can help maintain parasite numbers in the face of resource depletion. Early allocation decisions are predicted to have the greatest impact on parasite fitness. If the immune response saturates as parasite numbers increase, parasites should benefit from even longer delays prior to transmission investment. The presence of a competing strain selects for consistently lower levels of transmission investment and dramatically increased exploitation of the red blood cell resource. While we provide a detailed analysis of tradeoffs pertaining to malaria life history, our approach for identifying optimal plastic allocation strategies may be broadly applicable.     

Optimal growth strategies of larval helminths in their intermediate hosts

We consider optimal growth of larval stages in complex parasite life cycles where there is no constraint because of host immune responses. Our model predicts an individual's asymptotic size in its intermediate host, with and without competition from conspecific larvae. We match observed variations in larval growth patterns in pseudophyllid cestodes with theoretical predictions of our model. If survival of the host is vital for transmission, larvae should reduce asymptotic size as intensity increases, to avoid killing the host. The life history strategy (LHS) model predicts a size reduction <1/intensity, thus increasing the parasite burden on the host. We discuss whether body size of competing parasites is an evolved LHS or simply reflects resource constraints (RC) on growth fixed by the host, leading to a constant total burden with intensity. Growth under competition appears comparable with "the tragedy of the commons", much analysed in social sciences. Our LHS prediction suggests that evolution generates a solution that seems cooperative but is actually selfish.     

The role of sex in parasite dynamics: model simulations on transmission of Heligmosomoides polygyrus in populations of yellow-necked mice, Apodemus flavicollis

We investigated possible mechanisms that could cause sex-biased parasite transmission of the helminth Heligmosomoides polygyrus in its rodent host, Apodemus flavicollis, using a modelling approach. Two, not mutually exclusive, hypotheses were examined: that sex-biased parasite transmission is caused by differences in immunity that influence the success of free-living stages and/or is caused by sex differences in host behaviour and the dissemination of infective stages. Model simulations were compared with results from a field manipulation experiment of H. polygyrus in replicated populations of A. flavicollis. Simulations predicted the experimental field results, and both hypotheses explained the pattern observed. Transmission is male-biased if a male immune response increases fertility, hatching or survival of free-living stages. Alternatively, transmission is male-biased if their behavioural characteristics allow them to spread infective larvae in areas more frequently used by females. These results highlight that host sex is not only responsible for differences in parasite susceptibility, but may profoundly influence host-parasite interactions, resulting in a sex bias in parasite transmission.     

Host condition as a constraint for parasite reproduction

Environmental stress has been suggested to increase host susceptibility to infections and reduce host ability to resist parasite growth and reproduction, thus benefiting parasites. This prediction stems from expected costs of immune defence; hosts in poor condition should have less resources to be allocated to immune function. However, the alternative hypothesis for response to environmental stress is that hosts in poor condition provide less resources for parasites and/or suffer higher mortality, leading to reduced parasite growth, reproduction and survival. We contrasted these alternative hypotheses in a trematode–snail ( Diplostomum spathaceum – Lymnaea stagnalis ) system by asking: (1) how host condition affects parasite reproduction (amount and quality of produced transmission stages) and (2) how host condition affects the survival of infected host individuals. We experimentally manipulated host condition by starving the snails, and found that parasites produced fewer and poorer quality transmission stages in stressed hosts. Furthermore, starvation increased snail mortality. These findings indicate that in well-established trematode infections, reduced ability of immune allocation has no effect on host exploitation by parasites. Instead, deteriorating resources for the snail host can directly limit the amount of resources available for the parasite. This, together with increased host mortality, may have negative effects on parasite populations in the wild.     

Sex ratio adjustment in Plasmodium gallinaceum

The sex ratio of the avian malaria parasite, Plasmodium gallinaceum, was examined during the course of infection in its natural host, the chicken. Infections can have two possible outcomes: death of the host resulting from anaemia or self-cure and survival. In lethal infections the sex ratio remained female biased throughout, whereas in self-curing infections, the sex ratio became progressively less female biased. We examined the consequences of altering sex ratio for parasite transmission success using a theoretical fertilisation model and hypothesise that an immune response specifically effective against the male gametes would provide a selective explanation for the observed sex ratio adjustment. Previous studies have demonstrated that there is an efficient anti-gamete antibody response as an infection is cleared by the host. We performed in vitro mosquito infection studies comparing mosquito infection rates with naive serum replacement and showed that decomplemented serum from curing infections is transmission blocking, whereas serum from lethal infections is not. We discuss aspects of the malaria parasite-host interaction which might provide the selective pressure for such observed sex ratio adjustment.     

Modelling the immune response to malaria with ecological concepts: short-term behaviour against long-term equilibrium.

A model for the human immune response to the malaria parasite Plasmodium falciparum is used to analyse the dynamics of an infection within an individual patient. Previous models either looked at competition between two parasite genotypes or at one parasite clone and the immune response to it. This model describes the course of an infection caused by the blood stages of two parasite genotypes differing in reproductive rate and in the immune response they elicit. The interactions between the genotypes can be interpreted as exploitative competition for red blood cells. Interactions between omnipotent immune cells and parasites resemble a predator-prey relation. In analysing these kinds of models, classical theoretical ecology usually deals with long-term behaviours, i.e. looks for equilibria and conditions for coexistence. However, especially in endemic regions with ongoing transmission, an equilibrium state of infections is unlikely. When reinfections with another parasite genotype were considered, the short-term dynamics of the infection changed dramatically, depending on which genotype was first, when the second one appeared, and what kind of immune response was elicited. If the slow development of immunity to malaria really is due to its genotype specificity, the effects of superinfections will be of great importance.     

Host longevity and parasite species richness in mammals

Hosts and parasites co-evolve, with each lineage exerting selective pressures on the other. Thus, parasites may influence host life-history characteristics, such as longevity, and simultaneously host life-history may influence parasite diversity. If parasite burden causes increased mortality, we expect a negative association between host longevity and parasite species richness. Alternatively, if long-lived species represent a more stable environment for parasite establishment, host longevity and parasite species richness may show a positive association. We tested these two opposing predictions in carnivores, primates and terrestrial ungulates using phylogenetic comparative methods and controlling for the potentially confounding effects of sampling effort and body mass. We also tested whether increased host longevity is associated with increased immunity, using white blood cell counts as a proxy for immune investment. Our analyses revealed weak relationships between parasite species richness and longevity. We found a significant negative relationship between longevity and parasite species richness for ungulates, but no significant associations in carnivores or primates. We also found no evidence for a relationship between immune investment and host longevity in any of our three groups. Our results suggest that greater parasite burden is linked to higher host mortality in ungulates. Thus, shorter-lived ungulates may be more vulnerable to disease outbreaks, which has implications for ungulate conservation, and may be applicable to other short-lived mammals.     

Novel epidemiological model of gastrointestinal nematode infection to assess grazing cattle resilience by integrating host growth,parasite, grass and environmental dynamics

Gastrointestinal nematode (GIN) infections are ubiquitous and often cause morbidity and reduced performance in livestock. Emerging anthelmintic resistance and increasing change in climate patterns require evaluation of alternatives to traditional treatment and management practices. Mathematical models of parasite transmission between hosts and the environment have contributed towards the design of appropriate control strategies in ruminants, but have yet to account for relationships between climate, infection pressure, immunity, resources, and growth. Here, we develop a new epidemiological model of GIN transmission in a herd of grazing cattle, including host tolerance (body weight and feed intake), parasite burden and acquisition of immunity, together with weather-dependent development of parasite free-living stages, and the influence of grass availability on parasite transmission. Dynamic host, parasite and environmental factors drive a variable rate of transmission. Using literature sources, the model was parametrised for Ostertagia ostertagi, the prevailing pathogenic GIN in grazing cattle populations in temperate climates. Model outputs were validated on published empirical studies from first season grazing cattle in northern Europe. These results show satisfactory qualitative and quantitative performance of the model; they also indicate the model may approximate the dynamics of grazing systems under co-infection by O. ostertagi and Cooperia oncophora, a second GIN species common in cattle. In addition, model behaviour was explored under illustrative anthelmintic treatment strategies, considering impacts on parasitological and performance variables. The model has potential for extension to explore altered infection dynamics as a result of management and climate change, and to optimise treatment strategies accordingly. As the first known mechanistic model to combine parasitic and free-living stages of GIN with host feed-intake and growth, it is well suited to predict complex system responses under non-stationary conditions. We discuss the implications, limitations and extensions of the model, and its potential to assist in the development of sustainable parasite control strategies.     

The evolution of parasite virulence and transmission rate in a spatially structured population

If the transmission occurs through local contact of the individuals in a spatially structured population, the evolutionarily stable (ESS) traits of parasite might be quite different from what the classical theory with complete mixing predicts. In this paper, we theoretically study the ESS virulence and transmission rate of a parasite in a lattice-structured host population, in which the host can send progeny only to its neighboring vacant site, and the transmission occurs only in between the infected and the susceptible in the nearest-neighbor sites. Infected host is assumed to be infertile. The analysis based on the pair approximation and the Monte Carlo simulation reveal that the ESS transmission rate and virulence in a lattice-structured population are greatly reduced from those in completely mixing population. Unlike completely mixing populations, the spread of parasite can drive the host to extinction, because the local density of the susceptible next to the infected can remain high even when the global density of host becomes very low. This demographic viscosity and group selection between self-organized spatial clusters of host individuals then leads to an intermediate ESS transmission rate even if there is no tradeoff between transmission rate and virulence. The ESS transmission rate is below the region of parasite-driven extinction by a finite amount for moderately large reproductive rate of host; whereas, the evolution of transmission rate leads to the fade out of parasite for small reproductive rate, and the extinction of host for very large reproductive rate.     

Aging in personal and social immunity: do immune traits senesce at the same rate?

How much should an individual invest in immunity as it grows older? Immunity is costly and its value is likely to change across an organism's lifespan. A limited number of studies have focused on how personal immune investment changes with age in insects, but we do not know how social immunity, immune responses that protect kin, changes across lifespan, or how resources are divided between these two arms of the immune response. In this study, both personal and social immune functions are considered in the burying beetle, Nicrophorus vespilloides. We show that personal immune function declines (phenoloxidase levels) or is maintained (defensin expression) across lifespan in nonbreeding beetles but is maintained (phenoloxidase levels) or even upregulated (defensin expression) in breeding individuals. In contrast, social immunity increases in breeding burying beetles up to middle age, before decreasing in old age. Social immunity is not affected by a wounding challenge across lifespan, whereas personal immunity, through PO, is upregulated following wounding to a similar extent across lifespan. Personal immune function may be prioritized in younger individuals in order to ensure survival until reproductive maturity. If not breeding, this may then drop off in later life as state declines. As burying beetles are ephemeral breeders, breeding opportunities in later life may be rare. When allowed to breed, beetles may therefore invest heavily in “staying alive” in order to complete what could potentially be their final reproductive opportunity. As parental care is important for the survival and growth of offspring in this genus, staying alive to provide care behaviors will clearly have fitness payoffs. This study shows that all immune traits do not senesce at the same rate. In fact, the patterns observed depend upon the immune traits measured and the breeding status of the individual.     

Dose effects and density-dependent regulation of two microparasites of<Emphasis Type="Italic"> Daphnia magna</Emphasis>

Individual hosts constitute a limited resource for parasites, suggesting that density-dependent effects may play a role in within-host growth and parasite regulation. This hypothesis has been tested for several helminth parasites, but not for microparasites. We therefore examined dose-response patterns for the microparasitic bacterium Pasteuria ramosa and the fungus Metschnikowiella biscuspidata infecting the planktonic crustacean Daphnia magna. With increasing numbers of transmission stages administered to the host we found that host fecundity and survival and parasite transmission-stage production declined. Using a k-value analysis, a method that quantifies the strength of density dependence, we found for both parasites that density dependence acted at all doses, indicating the absence of a minimum density below which parasite fitness is density- independent. At low doses density was exactly compensated, but it was overcompensated at high doses. Overcompensation at high doses was weak for P. ramosa, but high for M. biscuspidata. At the two highest doses M. biscuspidata killed its hosts before any transmission stages were produced. Our data indicate that density dependence is expressed through retarded spore development in P. ramosa, but through both host mortality and reduced parasite fecundity in M. biscuspidata. A further experiment (P. ramosa only) revealed that in well-fed hosts more parasite transmission stages are produced than in poorly fed hosts, suggesting that competition for host resources retards P. ramosa development. Our data for P. ramosa, but not for M. biscuspidata, are largely consistent with assumptions made in models on microparasite epidemiology. We draw attention to the relevance of dose effects and within-host competition for the evolution of virulence. Received: 15 July 1999 / Accepted: 14 September 1999     

Condition-dependent virulence in a horizontally and vertically transmitted bacterial parasite

Virulence is a key component of parasite fitness. Its expression and selective value may not only depend on the features of the parasite's life cycle, but also on host genotype or environmental conditions. Using the freshwater ciliate Paramecium caudatum and its endonuclear bacterial parasite Holospora undulata , we measured variation in virulence (reduction in host division and survival), parasite load and fidelity of vertical transmission for (i) different stages of infection (associated with different opportunities for vertical and horizontal transmission), (ii) different host clones, and (iii) two food conditions. Later stages of infection dedicated to horizontal transmission were more virulent than earlier stages which rely on vertical transmission only. Besides, investment in horizontal transmission decreased the efficacy of vertical transmission, indicating a tradeoff between the two pathways. This may explain the phenotypic plasticity of transmission mode of this parasite. To some extent, virulence, parasite load and transmission fidelity varied with host clone identity and food treatment (higher virulence at low food). These results suggest that virulence is not a constant property of the parasite, and that a single (and simple) relationship between virulence and transmission does not exist.     

Host lifespan and the evolution of resistance to multiple parasites

Hosts are typically challenged by multiple parasites, but to date theory on the evolution of resistance has mainly focused on single infections. We develop a series of models that examine the impact of multiple parasites on the evolution of resistance under the assumption that parasites coexist at the host population scale as a consequence of superinfection. In this way, we are able to explicitly examine the impact of ecological dynamics on the evolutionary outcome. We use our models to address a key question of how host lifespan affects investment in resistance to multiple parasites. We show that investment in costly resistance depends on the specificity of the immune response and on whether or not the focal parasite leads to more acute infection than the co‐circulating parasite. A key finding is that investment in resistance always increases as the immune response becomes more general independently of whether it is the focal or the co‐circulating parasite that exploits the host most aggressively. Long‐lived hosts always invest more than short‐lived hosts in both general resistance and resistance that is specific to relatively acute focal parasites. However, for specific resistance to parasites that are less acute than co‐circulating parasites it is the short‐lived hosts that are predicted to invest most. We show that these results apply whatever the mode of defence, that is whether it is through avoidance or through increased recovery, with or without acquired immunity, or through acquired immunity itself. As a whole, our results emphasize the importance of considering multiple parasites in determining optimal immune investment in eco‐evolutionary systems.