Random parasite encounters coupled with condition-linked immunity of hosts generate parasite aggregation

Parasite aggregation is viewed as a natural law in parasite-host ecology but is a paradox insofar as parasites should follow the Poisson distribution if hosts are encountered randomly. Much research has focused on whether parasite aggregation in or on hosts is explained by aggregation of infective parasite stages in the environment, or by heterogeneity within host samples in terms of host responses to infection (e.g., through representation of different age classes of hosts). In this paper, we argue that the typically aggregated distributions of parasites may be explained simply. We propose that aggregated distributions can be derived from parasites encountering hosts randomly, but subsequently by parasites being 'lost' from hosts based on condition-linked escape or immunity of hosts. Host condition should be a normally distributed trait even among otherwise homogeneous sets of hosts. Our model shows that mean host condition and variation in host condition have different effects on the different metrics of parasite aggregation. Our model further predicts that as host condition increases, parasites become more aggregated but numbers of attending parasites are reduced overall and this is important for parasite population dynamics. The effects of deviation from random encounter are discussed with respect to the relationship between host condition and final parasite numbers.     

Parasite body size and interspecific variation in levels of aggregation among nematodes

The aggregation of parasites among individual hosts is one of the best documented features of parasite populations; we still do not know, however, why certain parasite species are more highly aggregated than other, related species. Here we search for a general explanation of interspecific variation in aggregation levels, based on the relationship between parasite body size and fecundity, transmission success, and intensity-dependent population regulation. We test the prediction that larger-bodied parasite species are more weakly aggregated than smaller-bodied related species, in a comparative analysis across parasitic nematode species. Across species, the variance-to-mean abundance ratio correlated negatively and significantly with nematode body sizes, as predicted. All other tests, however, including the more robust analyses controlling for phylogenetic influences, failed to support this result. This is mainly because the variance in infection levels is almost completely explained by mean parasite abundance. For this reason, it may prove difficult to identify a general biological explanation for interspecific variability in aggregation levels among parasites.     

THE EFFECTS OF HOST GENOTYPE AND SPATIAL DISTRIBUTION ON TREMATODE PARASITISM IN A BIVALVE POPULATION

A basic assumption underlying models of host-parasite coevolution is the existence of additive genetic variation among hosts for resistance to parasites. However, estimates of additive genetic variation are lacking for natural populations of invertebrates. Testing this assumption is especially important in view of current models that suggest parasites may be responsible for the evolution of sex, such as the Red Queen hypothesis. This hypothesis suggests that the twofold reproductive disadvantage of sex relative to parthenogenesis can be overcome by the more rapid production of rare genotypes resistant to parasites. Here I present evidence of significant levels of additive genetic variance in parasite resistance for an invertebrate host-parasite system in nature. Using families of the bivalve mollusc, Transennella tantilla, cultured in the laboratory, then exposed to parasites in the field, I quantified heritable variation in parasite resistance under natural conditions. The spatial distribution of outplanted hosts was also varied to determine environmental contributions to levels of parasite infection and to estimate potential interactions of host genotype with environment. The results show moderate but significant levels of heritability for resistance to parasites (h² = 0.36). The spatial distribution of hosts also significantly influenced parasite prevalence such that increased host aggregation resulted in decreased levels of parasite infection. Family mean correlations across environments were positive, indicating no genotype-environment interaction. Therefore, these results provide support for important assumptions underlying coevolutionary models of host-parasite systems.     

Partitioning the aggregation of parasites on hosts into intrinsic and extrinsic components via an extended Poisson-gamma mixture model

It is well known that parasites are often highly aggregated on their hosts such that relatively few individuals host the large majority of parasites. When the parasites are vectors of infectious disease, a key consequence of this aggregation can be increased disease transmission rates. The cause of this aggregation, however, is much less clear, especially for parasites such as arthropod vectors, which generally spend only a short time on their hosts. Regression-based analyses of ticks on various hosts have focused almost exclusively on identifying the intrinsic host characteristics associated with large burdens, but these efforts have had mixed results; most host traits examined have some small influence, but none are key. An alternative approach, the Poisson-gamma mixture distribution, has often been used to describe aggregated parasite distributions in a range of host/macroparasite systems, but lacks a clear mechanistic basis. Here, we extend this framework by linking it to a general model of parasite accumulation. Then, focusing on blacklegged ticks (Ixodes scapularis) on mice (Peromyscus leucopus), we fit the extended model to the best currently available larval tick burden datasets via hierarchical Bayesian methods, and use it to explore the relative contributions of intrinsic and extrinsic factors on observed tick burdens. Our results suggest that simple bad luck-inhabiting a home range with high vector density-may play a much larger role in determining parasite burdens than is currently appreciated.     

The effect of spatial heterogenity on the aggregation of ticks on white-footed mice

Parasites are often aggregated on a minority of the individuals in their host populations. Although host characteristics are commonly presumed to explain parasite aggregation on hosts, spatio-temporal aggregation of parasites during their host-seeking stages may have a dominant effect on the aggregation on hosts. We aimed to quantify, using mixed models, repeatability and autocorrelation analyses, the degree to which the aggregation of blacklegged ticks (Ixodes scapularis) on white-footed mice (Peromyscus leucopus) is influenced by spatio-temporal distributions of the host-seeking ticks and by heterogeneity among mice. Host-seeking ticks were spatially aggregated at both the larval and nymphal life-stages. However, this spatial aggregation accounted for little of the variation in larval and nymphal burdens observed on mice (3% and 0%, respectively). Conversely, mouse identity accounted for a substantial proportion of the variance in tick burdens. Mouse identity was a significant explanatory factor as the majority of ticks parasitized a consistent set of mice throughout the activity seasons. Of the characteristics associated with mouse identity investigated, only gender affected larval burdens, and body mass and home range sizes in males were correlated with nymphal burdens. These analyses suggest that aggregation of ticks on a minority of mice does not result from the distribution of host-seeking ticks but from characteristics of the hosts.     

Aggregation patterns of helminth populations in the introduced fish, Liza haematocheilus (Teleostei: Mugilidae): disentangling host–parasite relationships

A number of hypotheses exist to explain aggregated distributions, but they have seldom been used to investigate differences in parasite spatial distribution between native and introduced hosts. We applied two aggregation models, the negative binomial distribution and Taylor’s power law, to study the aggregation patterns of helminth populations from Liza haematocheilus across its native (Sea of Japan) and introduced (Sea of Azov) distribution ranges. In accordance with the enemy release hypothesis, we predicted that parasite populations in the introduced host range would be less aggregated than in the native host area, because aggregation is tightly constrained by abundance. Contrary to our expectation, aggregation of parasite populations was higher in the introduced host range. However, the analyses suggested that the effect of host introduction on parasite aggregation depends on whether parasite species, or higher level taxonomic groups, were acquired in or carried into the new area. The revealed similarity in the aggregation parameters of co-introduced monogeneans can be attributed to the repeatability and identity of the host–parasite systems. In contrast, the degree of aggregation differed markedly between regions for higher level taxa, which are represented by the native parasites in the Sea of Japan versus the acquired species in the Sea of Azov. We propose that the host species plays a crucial role in regulating infra-population sizes of acquired parasites due to the high rate of host-induced mortality. A large part of the introduced host population may remain uninfected due to their resistance to native naïve parasites. The core concept of our study is that the comparative analysis of aggregation patterns of parasites in communities and populations, and macroecological relationships, can provide a useful tool to reveal cryptic relationships in host–parasite systems of invasive hosts and their parasites.     

Effects of resource availability and social aggregation on the species richness of raccoon endoparasite infracommunities

Within populations the contact rate of hosts and infectious parasites is mediated by the interactions of resource availability, host density, and host behavior. Fluctuations in host density can result in the loss or extinction of a parasite population as contact rates between parasites and susceptible individuals drop below thresholds of parasite population persistence. Less understood is how changes in resources and the behavioral ecology of host populations affect parasites. We used food provisioning to experimentally assess the effects of resource availability and of inducing host aggregation on the endoparasite community of free‐ranging raccoons. Twelve independent raccoon populations were subjected to differential resource provisioning for two years: a clumped food distribution to aggregate hosts (n = 5 populations), a dispersed food distribution to add food without aggregating hosts (n = 3), and a no food treatment (n = 4). Remote cameras indicated that aggregation sizes were three to four times greater in aggregated versus non‐aggregated populations. We considered endoparasites with direct and indirect life cycles separately and determined the best‐fit models of parasite species richness in relation to host aggregation, food supplements, and host age and sex. Social aggregation had a negligible impact on the species richness of directly or indirectly transmitted parasites. However, food additions decreased the number of indirectly transmitted parasite species by 35% in the oldest age classes. These results suggest that while resource availability can influence the transmission of indirectly transmitted parasites, an examination of additional factors will be necessary to understand the role of host contact and factors that shape the community structure of endoparasites in natural environments.     

Host patch selection induced by parasitism: basic reproduction ratio r(0) and optimal virulence

The effects of parasites on the behavior of their hosts are well documented. For example, parasites may affect the habitat selection of the host individual. We used variables aggregation methods to investigate the way in which parasites affect the spatial pattern of susceptible hosts. We developed a simple epidemiological model, taking into account both the reproduction processes of hosts (density-dependent birth and death) and infection, considered separately on two different patches, and the migration of susceptible hosts between these two patches. We used the complete model of three equations to generate an aggregated model describing the dynamics of the combined susceptible and infected host populations. We obtained the basic reproduction ratio (R(0)) from the aggregated model, and then studied the effect of the migratory behavior of susceptible hosts on the ability of the parasite to invade the system. We also used the basic reproduction ratio to investigate the evolution of parasite virulence in relation to the migration decisions of susceptible hosts. We found that host investment in avoidance of the infected patch leads to an increase in optimal virulence if host investment is costly.     

From individual heterogeneity to population-level overdispersion: quantifying the relative roles of host exposure and parasite establishment in driving aggregated helminth distributions

In most host-parasite systems, variation in parasite burden among hosts drives transmission dynamics. Heavily infected individuals introduce disproportionate numbers of infective stages into host populations or surrounding environments, causing sharp increases in frequency of infection. Parasite aggregation within host populations may result from variation among hosts in exposure to infective propagules and probability of subsequent establishment of parasites in the host. This is because individual host heterogeneities contribute to a pattern of parasite overdispersion that emerges at the population level. We quantified relative roles of host exposure and parasite establishment in producing variation in parasite burdens, to predict which hosts are more likely to bear heavy burdens, using big brown bats (Eptesicus fuscus) and their helminths as a model system. We captured bats from seven colonies in Michigan and Indiana, USA, assessed their helminth burdens, and collected data on intrinsic and extrinsic variables related to exposure, establishment, or both. Digenetic trematodes had the highest prevalence and mean abundance while cestodes and nematodes had much lower prevalence and mean abundance. Structural equation modeling revealed that best-fitting models to explain variations in parasite burden included genetic heterozygosity and immunocompetence as well as distance to the nearest water source and the year of host capture. Thus, both differential host exposure and differential parasite establishment significantly influence heterogeneous helminth burdens, thus driving population-level patterns of parasite aggregation.     

Immunity,antigenic heterogeneity,and aggregation of helminth parasites

Empirical studies of helminth parasites reveal that the distribution of parasite burdens in their host populations is highly aggregated. This aggregation is fundamental to the ecology and epidemiology of helminth parasites. Results from a stochastic model predict that aggregation of helminth parasites is inversely related to the intensity of host immunity. Aggregation also decreases with antigenic heterogeneity and increases with heterogeneity in transmissibility among parasite strains. It is also found that the degree of aggregation is greater when immunity affects parasite fecundity than when immunity acts on host susceptibility. Potential relevance of this result for assessing the influence of vaccines that target either host susceptibility or parasite fecundity on the level of aggregation and consequent effects on drug resistance and disease prevalence are discussed.     

Modelling variability in lymphatic filariasis: macrofilarial dynamics in the Brugia pahangi--cat model.

A striking feature of lymphatic filariasis is the considerable heterogeneity in infection burden observed between hosts, which greatly complicates the analysis of the population dynamics of the disease. Here, we describe the first application of the moment closure equation approach to model the sources and the impact of this heterogeneity for macrofilarial population dynamics. The analysis is based on the closest laboratory equivalent of the life cycle and immunology of infection in humans--cats chronically infected with the filarial nematode Brugia pahangi. Two sets of long-term experiments are analysed: hosts given either single primary infections or given repeat infections. We begin by quantifying changes in the mean and aggregation of adult parasites (inversely measured by the negative binomial parameter, kappa in cohorts of hosts using generalized linear models. We then apply simple stochastic models to interpret observed patterns. The models and empirical data indicate that parasite aggregation tracks the decline in the mean burden with host age in primary infections. Conversely, in repeat infections, aggregation increases as the worm burden declines with experience of infection. The results show that the primary infection variability is consistent with heterogeneities in parasite survival between hosts. By contrast, the models indicate that the reduction in parasite variability with time in repeat infections is most likely due to the ''filtering'' effect of a strong, acquired immune response, which gradually acts to remove the initial variability generated by heterogeneities in larval mortality. We discuss this result in terms of the homogenizing effect of host immunity-driven density-dependence on macrofilarial burden in older hosts.     

A synthesis of experimental work on parasite local adaptation

The study of parasite local adaptation, whereby parasites perform better on sympatric hosts than on allopatric hosts and/or better on their own host population than do other parasites, is of great importance to both basic and applied biology. Theoretical examination of host-parasite coevolution predicts that parasite migration rate, generation time and virulence all contribute to the pattern of parasite local adaptation, such that parasites with greater dispersal ability, more frequent reproduction and/or high virulence ought to exhibit increased infectivity on local hosts. Here, we present a meta-analysis of experimental work from 57 host-parasite systems across 54 local adaptation studies to directly test theoretical predictions concerning the effect of each attribute on parasite adaptation. As expected, we find that studies of parasites with higher migration rates than their hosts report local adaptation, as measured by infection success, significantly more often than studies of parasites with relatively low migration rates. Furthermore, this synthesis serves to identify biases in the current body of work and highlight areas with the greatest need for further study. We emphasize the importance of unifying the field with regard to experimental methods, local adaptation definitions and reported statistics for cross-infection studies.     

A model of encounters between host and parasite populations

A simulation model of the encounter between host and parasite populations is described. The model is two-dimensional in that it represents hosts and parasites as sums of random numbers. It allows for the manipulation of host and parasite numbers, areas of interaction, congruity of geographic ranges, parasite infectivity, and reproduction, or non-reproduction, of the parasite. The model generates parasite distributions (number of hosts vs. parasite/host classes) and their parameters (prevalence, mean number of parasites/host, variance/mean ratio as a measure of aggregation), and thus reveals the manner in which these parameters vary under different encounter conditions, i.e. their "behavior". Simulation results indicated that the behavior of parasite population mean, prevalence, and degree of aggregation was primarily a function of the rate at which infective stages were supplied to the system. In cases in which infective stages were continuously available, prevalence rose rapidly to nearly 100%, with increasing infectivity and parasite numbers, and the populations were not particularly aggregated. When infective stages were introduced in single large waves, both mean and prevalence remained low and the parasite populations were highly aggregated. Model results were compared with published data sets. The latter were also seen to fall into the two general categories of parameter behavior.     

Strong density-dependent competition and acquired immunity constrain parasite establishment: implications for parasite aggregation

The vast majority of parasites exhibit an aggregated frequency distribution within their host population, such that most hosts have few or no parasites while only a minority of hosts are heavily infected. One exception to this rule is the trophically transmitted parasite Pterygodermatites peromysci of the white-footed mouse (Peromyscus leucopus), which is randomly distributed within its host population. Here, we ask: what are the factors generating the random distribution of parasites in this system when the majority of macroparasites exhibit non-random patterns? We hypothesise that tight density-dependent processes constrain parasite establishment and survival, preventing the build-up of parasites within individual hosts, and preclude aggregation within the host population. We first conducted primary infections in a laboratory experiment using white-footed mice to test for density-dependent parasite establishment and survival of adult worms. Secondary or challenge infection experiments were then conducted to investigate underlying mechanisms, including intra-specific competition and host-mediated restrictions (i.e. acquired immunity). The results of our experimental infections show a dose-dependent constraint on within-host-parasite establishment, such that the proportion of mice infected rose initially with exposure, and then dropped off at the highest dose. Additional evidence of density-dependent competition comes from the decrease in worm length with increasing levels of exposure. In the challenge infection experiment, previous exposure to parasites resulted in a lower prevalence and intensity of infection compared with primary infection of naïve mice; the magnitude of this effect was also density-dependent. Host immune response (IgG levels) increased with the level of exposure, but decreased with the number of worms established. Our results suggest that strong intra-specific competition and acquired host immunity operate in a density-dependent manner to constrain parasite establishment, driving down aggregation and ultimately accounting for the observed random distribution of parasites.     

The host specificity of Gyrodactylus salaris Malmberg (Platyhelminthes, Monogenea): susceptibility of Oncovhynchus mykiss (Walbaum) under experimental conditions

An experimental epidemiological approach was chosen to study the survival and infection dynamics of Gyrodactylus salaris on juvenile rainbow trout, Oncorhynchus mykiss , in the laboratory. A marked heterogeneity in the host stock was apparent. The rainbow trout could be divided into three groups on the basis of parasite survival and infection pattern on individually isolated fish: (1) hosts receptive to initial parasite attachment, but unreceptive to parasite establishment and reproduction; (2) hosts moderately susceptible to parasite establishment and reproduction, but which, after a period of restricted parasite population growth, responded, recovered and eliminated the parasites; and (3) hosts very susceptible to parasite infection and reproduction, but which, after a period of significant parasite population growth, responded, recovered and eliminated the parasites. These different patterns are considered to reflect genetic differences between host individuals. Parasite aggregation was also shown to be an important factor in the outcome of the host-parasite association. The parasites were finally eliminated on the individually isolated hosts, but not on hosts maintained in batches and so host population size and immigration of fresh. previously unexposed, hosts appeared to be important for growth and maintenance of the parasite population. The parasite was not found to cause host mortality. Rainbow trout was a suitable host for G. salaris , capable of transmitting the parasite to new localities as a consequence of stocking programmes or migratory behaviour.     

Evolutionary relationships, cospeciation, and host switching in avian malaria parasites

We used phylogenetic analyses of cytochrome b sequences of malaria parasites and their avian hosts to assess the coevolutionary relationships between host and parasite lineages. Many lineages of avian malaria parasites have broad host distributions, which tend to obscure cospeciation events. The hosts of a single parasite or of closely related parasites were nonetheless most frequently recovered from members of the same host taxonomic family, more so than expected by chance. However, global assessments of the relationship between parasite and host phylogenetic trees, using Component and ParaFit, failed to detect significant cospeciation. The event-based approach employed by TreeFitter revealed significant cospeciation and duplication with certain cost assignments for these events, but host switching was consistently more prominent in matching the parasite tree to the host tree. The absence of a global cospeciation signal despite conservative host distribution most likely reflects relatively frequent acquisition of new hosts by individual parasite lineages. Understanding these processes will require a more refined species concept for malaria parasites and more extensive sampling of parasite distributions across hosts. If parasites can disperse between allopatric host populations through alternative hosts, cospeciation may not have a strong influence on the architecture of host-parasite relationships. Rather, parasite speciation may happen more often in conjunction with the acquisition of new hosts followed by divergent selection between host lineages in sympatry. Detailed studies of the phylogeographic distributions of hosts and parasites are needed to characterize these events.     

Influence of host profitability and microenvironmental conditions on parasite specialization on a main and an alternative hosts

Parasite success depends on both host profitability and the microenvironment provided by the host, which together define host-parasite compatibility and can differ between hosts. We experimentally disentangled the effects of host profitability and microenvironmental conditions provided by nest material on the reproduction of a nest-based ectoparasite when exploiting its main and an alternative avian host species. Parasite reproductive performance was similar on both hosts when breeding in nests of their own species, suggesting no difference in host-parasite compatibility between hosts. The apparent parasite specialization could therefore result from differences in host-parasite encounter processes. However, when hosts were successful, the main host produced more young in infested nests, whereas the alternative host produced less; furthermore, host reproductive performance was higher in nests of the main host species, suggesting that this nest material alleviates parasitism cost. Therefore, our results suggest different evolutionary responses to parasites of the main and alternative hosts, with either higher tolerance or higher resistance, modulated by nest material.     

Minimal selfing, few clones, and no among-host genetic structure in a hermaphroditic parasite with asexual larval propagation

Little is known about actual mating systems in natural populations of parasites or about what constitutes the limits of a parasite deme. These parameters are interesting because they affect levels of genetic diversity, opportunities for local adaptation, and other evolutionary processes. We expect that transmission dynamics and the distribution of parasites among hosts should have a large effect on mating systems and demic structure, but currently we have mostly speculation and very few data. For example, infrapopulations (all the parasites in a single host) should behave as demes if parasite offspring are transmitted as a clump from host to host over several generations. However, if offspring are well mixed, then the parasite component population (all the parasites among a host population) would function as the deme. Similarly, low mean intensities or a high proportion of worms in single infections should increase the selfing rate. For species having an asexual amplification stage, transmission between intermediate and definitive (final) hosts will control the variance in clonal reproductive success, which in turn could have a large influence on effective sizes and rates of inbreeding. We examined demic structure, selfing rates, and the variance in clonal reproductive success in natural populations of Plagioporus shawi, a hermaphroditic trematode that parasitizes salmon. Overall levels of genetic diversity were very high. An a posteriori inference of population structure overwhelmingly supports the component population as the deme, rather than individual infrapopulations. Only a single pair of 597 adult individuals was identified as clones. Thus, the variance in clonal reproductive success was almost zero. Despite being hermaphroditic, P. shawi appears to be almost entirely outcrossing. Genetic estimates of selfing (<5%) were in accordance with the proportion of parasites from single infections. Thus, it appears that individual flukes outcross whenever possible and only resort to selfing when alone. Finally, our data support the hypothesis that aquatic transmission and the use of several intermediate hosts promotes high genetic diversity and well-mixed infrapopulations.     

Host characteristics and environmental factors differentially drive the burden and pathogenicity of an ectoparasite: a multilevel causal analysis

1. Understanding the ecological factors driving the burden and pathogenicity of parasites is challenging. Indeed, the dynamics of host-parasite interactions is driven by factors organized across nested hierarchical levels (e.g. hosts, localities), and indirect effects are expected owing to interactions between levels. 2. In this study, we combined Bayesian multilevel models, path analyses and a model selection procedure to account for these complexities and to decipher the relative effects of host- and environment-related factors on the burden and the pathogenicity of an ectoparasite (Tracheliastes polycolpus) on its fish host (Leuciscus leuciscus). We also tested the year-to-year consistency of the relationships linking these factors to the burden and the pathogenic effects of T. polycolpus. 3. We found significant relationships between the parasite burden and host-related factors: body length and age were positively related to parasite burden and heterozygous hosts displayed a higher parasite burden. In contrast, both host- and environment-related factors were linked to pathogenic effects. Pathogenicity was correlated negatively with host body length and positively with age; this illustrates that some factors (e.g. body length) showed inverse relationships with parasite burden and pathogenicity. Pathogenic effects were stronger in cooler upstream sites and where host density was lower. Path analyses revealed that these relationships between environment-related factors and pathogenic effects were direct and were not indirect relationships mediated by the host characteristics. Finally, we found that the strength and the shape of certain relationships were consistent across years, while they were clearly not for some others. 4. Our study illustrates that considering conjointly causal relationships among factors and the hierarchical structure of host-parasite interactions is appropriate for dissecting the complex links between hosts, parasites and their common environment.     

Relationship between structure and stability of a Halipegus occidualis component population in green frogs: a test of selective treatment

The aggregated nature of helminth parasite populations has led to the suggestion that selectively treating heavily infected hosts can efficiently reduce parasite abundance and morbidity within a host population. Moreover, it has been postulated that a selective treatment protocol might have long-term effects on parasite transmission by disrupting the stability attributed to aggregated parasite population distributions by theoretical models. Long-term investigation has demonstrated year-to-year consistency in the population dynamics of Halipegus occidualis in green frogs from Charlie's Pond, North Carolina. In 1996, removal of all but 1 worm from each frog with > or =15 worms reduced the estimated component worm population by 45%, thereby decreasing mean intensity and aggregation (variance-to-mean ratio) of H. occidualis in the frogs by 85% and 63%, respectively. The following year, mean intensity, aggregation, and host colonization trends returned to pretreatment levels, indicating no effect of worm removal and demonstrating the stability of this host-parasite system. Although this result might be attributable to inefficient treatment or the presence of infection reservoirs, it is suggested that parasite population stability in this system might be governed by prevalence rather than intensity of adult worms. Therefore, repeated selective treatment might effectively modify intensity-dependent morbidity in similar host-parasite systems but should not affect further parasite transmission.