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.     

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.     

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 dynamics of preferential host switching: Host phylogeny as a key predictor of parasite distribution*

Parasites often jump to and become established in a new host species. There is much evidence that the probability of such host shifts decreases with increasing phylogenetic distance between donor and recipient hosts, but the consequences of such preferential host switching remain little explored. We develop a computational model to investigate the dynamics of parasite host shifts in the presence of this phylogenetic distance effect. In this model, a clade of parasites evolves on an evolving clade of host species where parasites can cospeciate with their hosts, switch to new hosts, speciate within hosts or become extinct. Our model predicts that host phylogenies are major determinants of parasite distributions across trees. In particular, we predict that trees consisting of few large clades of host species and those with fast species turnover should harbor more parasites than trees with many small clades and those that diversify more slowly. Within trees, large clades are predicted to exhibit a higher fraction of infected species than small clades. We discuss our results in the light of recent cophylogenetic studies in a wide range of host–parasite systems.     

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.     

Patterns of parasite aggregation in the wild European rabbit (Oryctolagus cuniculus)

Understanding the factors controlling the distribution of parasites within their host population is fundamental to the wider understanding of parasite epidemiology and ecology. To explore changes in parasite aggregation, Taylor's power law was used to examine the distributions of five gut helminths of the wild rabbit. Aggregation was found to be a dynamic process that varied with year, season, host sex, age class, and myxomatosis. Yearly and seasonal changes are thought, in the main, to be the result of variations in weather conditions acting upon infectious stages (or intermediate hosts). Evidence in support of this was the comparatively low degree of fluctuation in the aggregation of the pinworm, Passalurus ambiguus, as the infectious stage of this parasite is likely to be less susceptible to environmental variation. Host age had a marked effect on the level of aggregation of all parasites, but this effect varied between parasite species. P. ambiguus, Trichostrongylus retortaeformis and Cittotaenia denticulata aggregation were lower in adult than juvenile rabbits whilst Graphidium strigosum and Mosgovoyia pectinata aggregation tended to increase with age. Host immunity is thought to be responsible for these differences. Differences in aggregation for different parasites were also seen when the rabbit population was split into males and females. Myxomatosis had a marked effect on helminth distribution with substantially less aggregation in rabbits showing clinical signs of the disease.     

Spatiotemporal patterns of avian host–parasite interactions in the face of biogeographical range expansions

Exploration of interactions between hosts and parasitic symbionts is important for our understanding of the temporal and spatial distribution of organisms. For example, host colonization of new geographical regions may alter levels of infections and parasite specificity, and even allow hosts to escape from co‐evolved parasites, consequently shaping spatial distributions and community structure of both host and parasite. Here we investigate the effect of host colonization of new regions and the elevational distribution of host–parasite associations between birds and their vector‐transmitted haemosporidian blood parasites in two geological and geographical settings: mountains of New Guinea and the Canary Islands. Our results demonstrate that bird communities in younger regions have significantly lower levels of parasitism compared to those of older regions. Furthermore, host–parasite network analyses demonstrate that blood parasites may respond differently after arriving to a new region, through adaptations that allow for either expanding (Canary Islands) or retaining (New Guinea) their host niches. The spatial prevalence patterns along elevational gradients differed in the two regions, suggesting that region‐specific biotic (e.g., host community) and abiotic factors (e.g., temperature) govern prevalence patterns. Our findings suggest that the spatiotemporal range dynamics in host–parasite systems are driven by multiple factors, but that host and parasite community compositions and colonization histories are of particular importance.     

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.     

Impact of forest size on parasite biodiversity: implications for conservation of hosts and parasites

Studies of biodiversity traditionally focus on charismatic megafauna. By comparison, little is known about parasite biodiversity. Recent studies suggest that co-extinction of host specific parasites with their hosts should be common and that parasites may even go extinct before their hosts. The few studies examining the relationship between parasite diversity and habitat quality have focused on parasites that require intermediate hosts and pathogens that require vectors to complete their life-cycles. Declines in parasite and pathogen richness in these systems could be due to the decline of any of the definitive hosts, intermediate hosts, or vectors. Here we focus on avian ectoparasites, primarily lice, which are host specific parasites with simple, direct, life-cycles. By focusing on these parasites we gain a clearer understanding of how parasites are linked to their hosts and their hosts’ environment. We compare parasite richness on birds from fragmented forests in southern China. We show that parasite richness correlates with forest size, even among birds that are locally common. The absence of some ectoparasite genera in small forests suggests that parasites can go locally extinct even if their hosts persist. Our data suggest that the conservation of parasite biodiversity may require preservation of habitat fragments that are sufficiently large to maintain parasite populations, not just their host populations.     

Ecological determinants of parasite acquisition by exotic fish species

Disease‐mediated threats posed by exotic species to native counterparts are not limited to introduced parasites alone, since exotic hosts frequently acquire native parasites with possible consequences for infection patterns in native hosts. Several biological and geographical factors are thought to explain both the richness of parasites in native hosts, and the invasion success of free‐living exotic species. However, the determinants of native parasite acquisition by exotic hosts remain unknown. Here, we investigated native parasite communities of exotic freshwater fish to determine which traits influence acquisition of native parasites by exotic hosts. Model selection suggested that five factors (total body length, time since introduction, phylogenetic relatedness to the native fish fauna, trophic level and native fish species richness) may be linked to native parasite acquisition by exotic fish, but 95% confidence intervals of coefficient estimates indicated these explained little of the variance in parasite richness. Based on R²‐values, weak positive relationships may exist only between the number of parasites acquired and either host size or time since introduction. Whilst our results suggest that factors influencing parasite richness in native host communities may be less important for exotic species, it seems that analyses of general ecological factors currently fail to adequately incorporate the physiological and immunological complexity of whether a given animal species will become a host for a new parasite.     

The scaling of total parasite biomass with host body mass

The selective pressure exerted by parasites on their hosts will to a large extent be influenced by the abundance or biomass of parasites supported by the hosts. Predicting how much parasite biomass can be supported by host individuals or populations should be straightforward: ultimately, parasite biomass must be controlled by resource supply, which is a direct function of host metabolism. Using comparative data sets on the biomass of metazoan parasites in vertebrate hosts, we determined how parasite biomass scales with host body mass. If the rate at which host resources are converted into parasite biomass is the same as that at which host resources are channelled toward host growth, then on a log-log plot parasite biomass should increase with host mass with a slope of 0.75 when corrected for operating temperature. Average parasite biomass per host scaled with host body mass at a lower rate than expected (across 131 vertebrate species, slope=0.54); this was true independently of phylogenetic influences and also within the major vertebrate groups separately. Since most host individuals in a population harbour a parasite load well below that allowed by their metabolic rate, because of the stochastic nature of infection, it is maximum parasite biomass, and not average biomass, that is predicted to scale with metabolic rate among host species. We found that maximum parasite biomass scaled isometrically (i.e., slope=1) with host body mass. Thus, larger host species can potentially support the same parasite biomass per gram of host tissues as small host species. The relationship found between maximum parasite biomass and host body mass, with its slope greater than 0.75, suggests that parasites are not like host tissues: they are able to appropriate more host resources than expected from metabolically derived host growth rates.     

Does the early frog catch the worm? Disentangling potential drivers of a parasite age–intensity relationship in tadpoles

The manner in which parasite intensity and aggregation varies with host age can provide insights into parasite dynamics and help identify potential means of controlling infections in humans and wildlife. A significant challenge is to distinguish among competing mechanistic hypotheses for the relationship between age and parasite intensity or aggregation. Because different mechanisms can generate similar relationships, testing among competing hypotheses can be difficult, particularly in wildlife hosts, and often requires a combination of experimental and model fitting approaches. We used field data, experiments, and model fitting to distinguish among ten plausible drivers of a curvilinear age–intensity relationship and increasing aggregation with host age for echinostome trematode infections of green frogs. We found little support for most of these proposed drivers but did find that the parsimonious explanation for the observed age–intensity relationship was seasonal exposure to echinostomes. The parsimonious explanation for the aggregated distribution of parasites in this host population was heterogeneity in exposure. A predictive model incorporating seasonal exposure indicated that tadpoles hatching early or late in the breeding season should have lower trematode burdens at metamorphosis, particularly with simulated warmer climates. Application of this multi-pronged approach (field surveys, lab experiments, and modeling) to additional parasite–host systems could lead to discovery of general patterns in the drivers of parasite age–intensity and age–distribution relationships.     

How parasite interaction strategies alter virulence evolution in multi‐parasite communities

The majority of organisms host multiple parasite species, each of which can interact with hosts and competitors through a diverse range of direct and indirect mechanisms. These within‐host interactions can directly alter the mortality rate of coinfected hosts and alter the evolution of virulence (parasite‐induced host mortality). Yet we still know little about how within‐host interactions affect the evolution of parasite virulence in multi‐parasite communities. Here, we modeled the virulence evolution of two coinfecting parasites in a host population in which parasites interacted through cross immunity, immune suppression, immunopathology, or spite. We show (1) that these within‐host interactions have different effects on virulence evolution when all parasites interact with each other in the same way versus when coinfecting parasites have unique interaction strategies, (2) that these interactions cause the evolution of lower virulence in some hosts, and higher virulence in other hosts, depending on the hosts infection status, and (3) that for cross immunity and spite, whether parasites increase or decrease the evolutionarily stable virulence in coinfected hosts depended on interaction strength. These results improve our understanding of virulence evolution in complex parasite communities, and show that virulence evolution must be understood at the community scale.     

Phylogenetic signals in host–parasite associations for Neotropical bats and Nearctic desert rodents

Hosts and their parasites have strong ecological and evolutionary relationships, with hosts representing habitats and resources for parasites. In the present study, we use approaches developed to evaluate the statistical dependence of species trait values on phylogenetic relationships to determine whether host–parasite relationships (i.e. parasite infections) are contingent on host phylogeny. If host–parasite relationships are contingent on the ability of hosts to provide habitat or resources to parasites, and if host phylogeny is an effective surrogate for among‐host variation in habitat and resource quality, host–parasite relationships should evince phylogenetic signals (i.e. be contingent on host phylogeny). Because the strength of ecological relationships between parasites and their hosts may affect the likelihood of phylogenetic signals occurring in host–parasite relationships, we hypothesized that (1) host specificity would be positively correlated with the strength of phylogenetic signals and (2) the strength of phylogenetic signals will be greater for parasites that rely more on their host throughout their life cycle. Analyses were conducted for ectoparasites from tropical bats and for ectoparasites, helminths, and coccidians from desert rodents. Phylogenetic signals were evaluated for parasite presence and for parasite prevalence. The frequency of phylogenetic signal occurrence was similar for parasite presence and prevalence, with a signal detected in 24–27% of cases at the species level and in 67% and 15% of cases at the genus level for parasites of bats and rodents, respectively. No differences in signal strength or the likelihood of detecting a signal existed between groups of parasites. Phylogenetic signal strength was correlated with host specificity, suggesting that mechanisms increasing host specificity also increase the likelihood of a phylogenetic signal in host use by parasites. Differences in the transmission mode did not affect signal strength or the likelihood of detecting a signal, indicating that variation in host switching opportunities associated with the transmission mode does not affect signal strength.     

Whether larval amphibians school does not affect the parasite aggregation rule: testing the effects of host spatial heterogeneity in field and experimental studies

Almost all macroparasites show over‐dispersed infections within natural host populations such that most parasites are distributed among a few heavily‐infected individuals. Despite the importance of parasite aggregation for understanding system stability, the potential for population regulation, and super‐spreading events, many questions persist about its underlying drivers. Theoretically, aggregation results from heterogeneity in host exposure, resistance, and tolerance. However, few studies have examined how host spatial arrangement – which likely affects both parasite encounter and density‐dependent interactions – influences infection and dispersion, representing a critical gap in our current knowledge regarding the possible drivers of parasite aggregation. Using field data from over 165 ponds and 8000 hosts, we evaluated how the spatial clustering of amphibian larvae within ponds 1) varied among different amphibian species, and 2), affected the distribution of parasites within the host population using Taylor's power law. A complementary mesocosm experiment used field‐guided manipulations of the spatial arrangement of larval amphibians to create a gradient in host clustering while controlling host density, thereby testing for spatial effects on both infection success and aggregation by three different trematode species. Our field data indicated that larval amphibians exhibited significant spatial clustering that was well captured by Taylor's power law (R² 0.92 to 0.97 for different host species), but the residual variation only weakly correlated with observed patterns of trematode parasite over‐dispersion. Correspondingly, experimental manipulation of host clustering had no effects on parasite infection success or the degree of parasite aggregation among cages or mesocosms. Given the importance of parasite over‐dispersion for host populations and disease dynamics, we advocate for further investigations of host and parasite spatial aggregation, particularly studies that incorporate and/or control for heterogeneity in exposure and susceptibility.     

Modelling parasite aggregation: disentangling statistical and ecological approaches

The overdispersion in macroparasite infection intensity among host populations is commonly simulated using a constant negative binomial aggregation parameter. We describe an alternative to utilising the negative binomial approach and demonstrate important disparities in intervention efficacy projections that can come about from opting for pattern-fitting models that are not process-explicit. We present model output in the context of the epidemiology and control of soil-transmitted helminths due to the significant public health burden imposed by these parasites, but our methods are applicable to other infections with demonstrable aggregation in parasite numbers among hosts.     

Coinfecting parasites can modify fluctuating selection dynamics in host–parasite coevolution

Genetically specific interactions between hosts and parasites can lead to coevolutionary fluctuations in their genotype frequencies over time. Such fluctuating selection dynamics are, however, expected to occur only under specific circumstances (e.g., high fitness costs of infection to the hosts). The outcomes of host–parasite interactions are typically affected by environmental/ecological factors, which could modify coevolutionary dynamics. For instance, individual hosts are often infected with more than one parasite species and interactions between them can alter host and parasite performance. We examined the potential effects of coinfections by genetically specific (i.e., coevolving) and nonspecific (i.e., generalist) parasite species on fluctuating selection dynamics using numerical simulations. We modeled coevolution (a) when hosts are exposed to a single parasite species that must genetically match the host to infect, (b) when hosts are also exposed to a generalist parasite that increases fitness costs to the hosts, and (c) when coinfecting parasites compete for the shared host resources. Our results show that coinfections can enhance fluctuating selection dynamics when they increase fitness costs to the hosts. Under resource competition, coinfections can either enhance or suppress fluctuating selection dynamics, depending on the characteristics (i.e., fecundity, fitness costs induced to the hosts) of the interacting parasites.