Resistance and tolerance in a host plant-holoparasitic plant interaction: genetic variation and costs

Host organisms are believed to evolve defense mechanisms (i.e., resistance and/or tolerance) under selective pressures exerted by natural enemies. A prerequisite for the evolution of resistance and tolerance is the existence of genetic variation in these traits for natural selection to act. However, selection for resistance and/or tolerance may be constrained by negative genetic correlations with other traits that affect host fitness. We studied genetic variation in resistance and tolerance against parasitic infection and the potential fitness costs associated with these traits using a novel study system, namely the interaction between a flowering plant and a parasitic plant. In this system, parasitic infection has significant negative effects on host growth and reproduction and may thus act as a selective agent. We conducted a greenhouse experiment in which we grew host plants, Urtica dioica, that originated from a single natural population and represented 20 maternal families either uninfected or infected with the holoparasitic dodder, Cuscuta europaea. that originated from the same site. We calculated correlations among resistance, tolerance, and host performance to test for costs of resistance and tolerance. We measured resistance as parasite performance (quantitative resistance) and tolerance as the slopes of regressions relating the vegetative and reproductive biomass of host plants to damage level (measured as parasite biomass). We observed significant differences among host families in parasite resistance and in parasite tolerance in terms of reproductive biomass, a result that suggests genetic variation in these traits. Furthermore, we found differences in resistance and tolerance between female and male host plants. In addition, the correlations indicate costs of resistance in terms of host growth and reproduction and costs of tolerance in terms of host reproduction. Our results thus indicate that host tolerance and resistance can evolve as a response to infection by a parasitic plant and that costs of resistance and tolerance may be one factor maintaining genetic variation in these traits.     

Genetic variation in resistance, but not tolerance, to a protozoan parasite in the monarch butterfly

Natural selection should strongly favour hosts that can protect themselves against parasites. Most studies on animals so far have focused on resistance, a series of mechanisms through which hosts prevent infection, reduce parasite growth or clear infection. However, animals may instead evolve tolerance, a defence mechanism by which hosts do not reduce parasite infection or growth, but instead alleviate the negative fitness consequences of such infection and growth. Here, we studied genetic variation in resistance and tolerance in the monarch butterfly (Danaus plexippus) to its naturally occurring protozoan parasite, Ophryocystis elektroscirrha. We exposed 560 monarch larvae of 19 different family lines to one of five different parasite inoculation doses (0, 1, 5, 10 and 100 infective spores) to create a range of parasite loads in infected butterflies. We then used two proxies of host fitness (adult lifespan and body mass) to quantify: (i) qualitative resistance (the ability to prevent infection; also known as avoidance or anti-infection resistance); (ii) quantitative resistance (the ability to limit parasite growth upon infection; also known as control or anti-growth resistance); and (iii) tolerance (the ability to maintain fitness with increasing parasite infection intensity). We found significant differences among host families in qualitative and quantitative resistance, indicating genetic variation in resistance. However, we found no genetic variation in tolerance. This may indicate that all butterflies in our studied population have evolved maximum tolerance, as predicted by some theoretical models.     

Host species,and not environment,predicts variation in blood parasite prevalence,distribution, and diversity along a humidity gradient in northern South America

Environmental factors strongly influence the ecology and evolution of vector‐borne infectious diseases. However, our understanding of the influence of climatic variation on host–parasite interactions in tropical systems is rudimentary. We studied five species of birds and their haemosporidian parasites (Plasmodium and Haemoproteus) at 16 sampling sites to understand how environmental heterogeneity influences patterns of parasite prevalence, distribution, and diversity across a marked gradient in water availability in northern South America. We used molecular methods to screen for parasite infections and to identify parasite lineages. To characterize spatial heterogeneity in water availability, we used weather‐station and remotely sensed climate data. We estimated parasite prevalence while accounting for spatial autocorrelation, and used a model selection approach to determine the effect of variables related to water availability and host species on prevalence. The prevalence, distribution, and lineage diversity of haemosporidian parasites varied among localities and host species, but we found no support for the hypothesis that the prevalence and diversity of parasites increase with increasing water availability. Host species and host × climate interactions had stronger effects on infection prevalence, and parasite lineages were strongly associated with particular host species. Because climatic variables had little effect on the overall prevalence and lineage diversity of haemosporidian parasites across study sites, our results suggest that independent host–parasite dynamics may influence patterns in parasitism in environmentally heterogeneous landscapes.     

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.     

Making the right choice: testing the drivers of asymmetric infections within hosts and their consequences for pathology

Despite the ubiquity of bilateral symmetry among animals, a long‐standing mystery centers on why parasites that infect paired organs often do so non‐randomly. Examples from diverse host and parasite taxa continue to accumulate, yet little is known about their causes or implications for host–parasite fitness. We combined field surveys, experimental infections, and parasite choice assays to evaluate both competing explanations for – and consequences of – asymmetric infections of amphibian kidneys by echinostome trematodes, which are widespread and potentially pathogenic infections of larval amphibians. Samples from 6001 hosts representing 26 species indicated that echinostome infections exhibit a consistent, right‐kidney bias, with ? 62% of parasites in the right kidney. This pattern could not be explained by variation in kidney size or total infection. Experimental infections of three anuran species reproduced this pattern, with 64% of infections in the right kidney, and indicated it was not the result of differential host or parasite mortality. Based on sequential infection experiments and parasite choice assays, we further showed that earlier infections did not affect the distribution of subsequently colonizing parasites and that echinostome cercariae followed host‐derived cues rather than exhibiting congenital ‘sidedness’. We advance the hypothesis that variation in the position of the right kidney along the anterior–posterior axis controls cue strength in the right nephric duct and thus determines parasite encystment. Correspondingly, anatomical measurements from a subset of larval amphibian hosts revealed that the relative position of the right kidney explained 83% of the variation in infection bias, with no additional contributions associated with kidney volume or host size. We also show that the degree of right‐kidney bias associated positively with host growth in experiments. Morphological asymmetries could therefore function as a unique form of tolerance to mitigate the consequences of infection, despite the oft‐cited costs of asymmetry for mate selection and enemy vulnerability.     

Patterns of parasite community dissimilarity: the significant role of land use and lack of distance‐decay in a bat–helminth system

Increasing community dissimilarity across geographic distance has been described for a wide variety of organisms and understanding its underlying causes is key to understanding mechanisms driving patterns of biodiversity. Both niche‐based and neutral processes may produce a distance decay relationship; however, disentangling their relative influence requires simultaneous examination of multiple potential drivers. Parasites represent a unique opportunity in which to study distance decay because community dissimilarity may depend on environmental requirements and dispersal capability of parasites as well also those of their hosts. We used big brown bats Eptesicus fuscus and their intestinal helminths to investigate: 1) independent contributions of geographic and environmental distances on dissimilarity of intestinal helminth component communities between populations of big brown bats; 2) which environmental variables best explained variation in community dissimilarity; and 3) whether similar patterns of decay with geographic or environmental distance were observed for within‐host population and within‐individual host parasite communities. We used compositional measures of community dissimilarity to examine how parasite communities may change with geographic distance and varying environmental conditions. Non‐spatial variables strongly influenced compositional parasite community dissimilarity over multiple community scales, and we observed little evidence for spatial processes such as distance decay. Environment surrounding roost sites better predicted helminth community dissimilarity than any other class of variables and landcover classes representing anthropogenic modification consistently explained variation in community structure. Our results indicate that human disturbance drives significant patterns of parasite community dissimilarity, most likely by changing the presence or abundance of intermediate hosts in an area.     

Heterogeneous selection promotes maintenance of polymorphism in host–parasite interactions

Polymorphism in loci affecting host resistance and parasite virulence is characteristic for nearly all species and this genetic variation is considered to have profound consequences for the patterns of disease incidence, prevalence and evolution. The gene-for-gene (GFG) system is a well-characterized genetic interaction of host recognition and parasite antigenic loci for a wide range of plant-parasite interactions. Long-term maintenance of polymorphism in GFG systems has remained puzzling for both theoreticians and empiricists. Traditionally this diversity has been explained by tradeoffs with other life-history traits closely linked with fitness, yet empirical evidence for such costs has remained mixed. Here we argue that incorporating simple ecological reality – spatial structuring and gradient of environmental conditions – into host–parasite research will help us understand how polymorphism is maintained. While environmental conditions (biotic and abiotic factors) have been studied in depth in plant pathology for their influence on disease severity and plant yield, they have been rarely set into an evolutionary framework. We briefly review recent data on natural plant–parasite metapopulations and theoretical models moving from single population models towards metapopulation theory to reveal in just how many ways spatial structuring may affect the coevolutionary process. We clarify also how spatially heterogeneous selection, through G×E (or G×G×E) interactions, may be particularly important for natural host–parasite interactions and suggest that this provides the unifying ground upon which future theoretical and empirical work should be build on.     

Local adaptation of a parasite to solar radiation impacts disease transmission potential,spore yield,and host fecundity*

Environmentally transmitted parasites spend time in the abiotic environment, where they are subjected to a variety of stressors. Learning how they face this challenge is essential if we are to understand how host–parasite interactions may vary across environmental gradients. We used a zooplankton–bacteria host–parasite system where availability of sunlight (solar radiation) influences disease dynamics to look for evidence of parasite local adaptation to sunlight exposure. We also examined how variation in sunlight tolerance among parasite strains impacted host reproduction. Parasite strains collected from clearer lakes (with greater sunlight penetration) were most tolerant of the negative impacts of sunlight exposure, suggesting local adaptation to sunlight conditions. This adaptation came with both a cost and a benefit for parasites: parasite strains from clearer lakes produced relatively fewer transmission stages (spores) but these strains were more infective. After experimental sunlight exposure, the most sunlight-tolerant parasite strains reduced host fecundity just as much as spores that were never exposed to sunlight. Sunlight availability varies greatly among lakes around the world. Our results suggest that the selective pressure sunlight exposure exerts on parasites may impact both parasite and host fitness, potentially driving variation in disease epidemics and host population dynamics across sunlight availability gradients.     

Disentangling the influence of parasite genotype, host genotype and maternal environment on different stages of bacterial infection in Daphnia magna

Individuals naturally vary in the severity of infectious disease when exposed to a parasite. Dissecting this variation into genetic and environmental components can reveal whether or not this variation depends on the host genotype, parasite genotype or a range of environmental conditions. Complicating this task, however, is that the symptoms of disease result from the combined effect of a series of events, from the initial encounter between a host and parasite, through to the activation of the host immune system and the exploitation of host resources. Here, we use the crustacean Daphnia magna and its parasite Pasteuria ramosa to show how disentangling genetic and environmental factors at different stages of infection improves our understanding of the processes shaping infectious disease. Using compatible host-parasite combinations, we experimentally exclude variation in the ability of a parasite to penetrate the host, from measures of parasite clearance, the reduction in host fecundity and the proliferation of the parasite. We show how parasite resistance consists of two components that vary in environmental sensitivity, how the maternal environment influences all measured aspects of the within-host infection process and how host-parasite interactions following the penetration of the parasite into the host have a distinct temporal component.     

Invasion and persistence of plant parasites in a spatially structured host population

Spatial structure is of central importance in the dynamics of plant-parasite interactions and is imposed by the growth habit and distribution of host plants and by parasite dispersal which is frequently restricted. To investigate the effects of spatial heterogeneity on the dynamics of plant parasites we introduce a simple model for epidemic development within a spatially structured host population. Here the host population is subdivided into a number of patches which are linked to allow for transmission from one patch to another with the connections defining the spatial structure of the host population. Three key parameters are identified that play a critical role in the ability of the parasite to invade and persist within the host population: the within-patch parasite basic reproductive number which characterises the infection dynamics at the local spatial scale; and the neighbourhood of interaction which describes which patches interact with which and the strength of coupling between patches within the neighbourhood which together characterise the spread of the parasite over larger spatial scales. Using both deterministic and stochastic formulations of the model, we investigate how the thresholds and probabilities of invasion and persistence are affected by these parameters, by demographic stochasticity and by differences in the initial level of infection.     

Fecundity compensation and tolerance to a sterilizing pathogen in Daphnia

Hosts are armed with several lines of defence in the battle against parasites: they may prevent the establishment of infection, reduce parasite growth once infected or persevere through mechanisms that reduce the damage caused by infection, called tolerance. Studies on tolerance in animals have focused on mortality, and sterility tolerance has not been investigated experimentally. Here, we tested for genetic variation in the multiple steps of defence when the invertebrate Daphnia magna is infected with the sterilizing bacterial pathogen Pasteuria ramosa: anti-infection resistance, anti-growth resistance and the ability to tolerate sterilization once infected. When exposed to nine doses of a genetically diverse pathogen inoculum, six host genotypes varied in their average susceptibility to infection and in their parasite loads once infected. How host fecundity changed with increasing parasite loads did not vary between genotypes, indicating that there was no genetic variation for this measure of fecundity tolerance. However, genotypes differed in their level of fecundity compensation under infection, and we discuss how, by increasing host fitness without targeting parasite densities, fecundity compensation is consistent with the functional definition of tolerance. Such infection-induced life-history shifts are not traditionally considered to be part of the immune response, but may crucially reduce harm (in terms of fitness loss) caused by disease, and are a distinct source of selection on pathogens.     

The Relationship between Parasite Fitness and Host Condition in an Insect - Virus System

Research in host-parasite evolutionary ecology has demonstrated that environmental variation plays a large role in mediating the outcome of parasite infection. For example, crowding or low food availability can reduce host condition and make them more vulnerable to parasite infection. This observation that poor-condition hosts often suffer more from parasite infection compared to healthy hosts has led to the assumption that parasite productivity is higher in poor-condition hosts. However, the ubiquity of this negative relationship between host condition and parasite fitness is unknown. Moreover, examining the effect of environmental variation on parasite fitness has been largely overlooked in the host-parasite literature. Here we investigate the relationship between parasite fitness and host condition by using a laboratory experiment with the cabbage looper Trichoplusia ni and its viral pathogen, AcMNPV, and by surveying published host-parasite literature. Our experiments demonstrated that virus productivity was positively correlated with host food availability and the literature survey revealed both positive and negative relationships between host condition and parasite fitness. Together these data demonstrate that contrary to previous assumptions, parasite fitness can be positively or negatively correlated with host fitness. We discuss the significance of these findings for host-parasite population biology.     

Ecology and dynamics of the blood parasite, Hepatozoon tuatarae (Apicomplexa), in tuatara (Sphenodon punctatus) on Stephens Island, New Zealand

We explored infection patterns and temporal dynamics of the protozoan blood parasite Hepatozoon tuatarae (Apicomplexa) infecting the tuatara (Sphenodon punctatus), a protected reptile living on Stephens Island, New Zealand. In March 2006, we surveyed tuatara in five study sites to examine spatial variation in infection prevalence, and four times, from May 2005 to November 2006, we recaptured marked individuals within three study sites to examine the temporal dynamics of infection. We also examined how blood-parasite infection patterns were influenced by host sex, body size, and host infestation with ticks (Amblyomma sphenodonti) and mites (Neotrombicula spp.), which are potential vectors of the blood parasite. Infection prevalence (16.9-24% infected) and intensity (<0.01-0.1% blood cells infected) were low in all samples. Infection intensity varied among the five sampled sites in March 2006, but prevalence did not. Neither infection prevalence nor intensity varied with time, and infections were detected in consecutive samples from recaptured individuals for up to 18 mo. Neither survey showed an influence of host sex on infection, but both surveys showed infection intensity declined with increasing host body size, as did infection prevalence in the spatial survey. In the temporal survey, we found a positive relationship between the tick numbers on hosts and blood-parasite infection intensity, which were stronger in two of the sampling periods and among larger hosts. These data suggest that exposure and susceptibility to infection decreases with host size and that ticks, but not mites, are probably the vectors in this ancient host-parasite association of a long-lived (>50 yr) host.     

Cooperation,virulence and siderophore production in bacterial parasites

Kin selection theory predicts that the damage to a host resulting from parasite infection (parasite virulence) will be negatively correlated to the relatedness between parasites within the host. This occurs because a lower relatedness leads to greater competition for host resources, which favours rapid growth to achieve greater relative success within the host, and that higher parasite growth rate leads to higher virulence. We show that a biological feature of bacterial infections can lead to the opposite prediction: a positive correlation between relatedness and virulence. This occurs because a high relatedness can favour greater (cooperative) production of molecules that scavenge iron (siderophores), which results in higher growth rates and virulence. More generally, the same underlying idea can predict a positive relationship between relatedness and virulence in any case where parasites can cooperate to increase their growth rate; other examples include immune suppression and the production of biofilms to aid colonization.     

Spatial determinants of infection risk in a multi‐species avian malaria system

Spatially‐variable processes can be an important element of host–parasite interactions, but their longer term demographic and evolutionary effects depend on the magnitude of variation in space, the scale at which variation occurs and the degree to which such processes are temporally stable. Here, we use multiple years of data from a study of two closely related tit species (Paridae), infected with two congeneric species of avian malaria (Plasmodium), to evaluate the roles of extrinsic and intrinsic factors in driving spatial heterogeneity in infection risk, and to address questions of scale and temporal stability in these vector‐driven host–parasite interactions. We show that the two malaria parasite species exhibit markedly different spatial epidemiology: P. relictum infections are effectively randomly distributed in space, with no temporal consistency, whereas P. circumflexum infections exhibit pronounced spatial structuring that is stable over the six years of this study and similar in both host species. We show that both conspecific and heterospecific host density contribute to elevated infection risk, but that the main determinants of elevated risk of P. circumflexum infection risk are habitat features probably associated with vector distribution and abundance. We discuss the implications of these findings, both for our understanding of the epidemiology of malaria in the wild, but also in terms of the longer‐term evolutionary and demographic consequences that spatially variable parasite‐mediated selection may have on host populations.     

Expression of parasite genetic variation changes over the course of infection: implications of within-host dynamics for the evolution of virulence

How infectious disease agents interact with their host changes during the course of infection and can alter the expression of disease-related traits. Yet by measuring parasite life-history traits at one or few moments during infection, studies have overlooked the impact of variable parasite growth trajectories on disease evolution. Here we show that infection-age-specific estimates of host and parasite fitness components can reveal new insight into the evolution of parasites. We do so by characterizing the within-host dynamics over an entire infection period for five genotypes of the castrating bacterial parasite Pasteuria ramosa infecting the crustacean Daphnia magna. Our results reveal that genetic variation for parasite-induced gigantism, host castration and parasite spore loads increases with the age of infection. Driving these patterns appears to be variation in how well the parasite maintains control of host reproduction late in the infection process. We discuss the evolutionary consequences of this finding with regard to natural selection acting on different ages of infection and the mechanism underlying the maintenance of castration efficiency. Our results highlight how elucidating within-host dynamics can shed light on the selective forces that shape infection strategies and the evolution of virulence.     

DYNAMIC AND GENETIC CONSEQUENCES OF VARIATION IN HORIZONTAL TRANSMISSION FOR A MICROPARASITIC INFECTION

Transmission to a new host is a critical step in the life cycle of a parasite. Variation in the characteristics of the transmission process, for example, due to host demography, is assumed to select for different variants of the parasite. We have experimentally tested how variation in the time to transmission (early or late after infection) and exposure to adverse conditions outside the host (immediate or delayed contact with new host) interact to determine the success of the infection in the next host, using the trypanosome Crithidia bombi infecting its bumblebee host, Bombus terrestris. These two experimentally manageable steps mimic the processes of within- and among-host selection for the parasite. We found that early transmission led to higher infection success in the next host as did immediate contact with the new host. However, there was no interaction between the two parameters as would be expected if early-transmitted variants, resulting from rapid multiplication within the host, would be less adapted to the conditions encountered during the between-host transfer or infection of the next host. Furthermore, typing the genetic variability of the parasites with microsatellites showed that the four different transmission routes of our experiment selected for different degrees of allelic diversity of the infecting parasite populations. The results support the idea that variation in the transmission process selects for different genotypic variants of the parasite. At the same time, the relationship of allelic diversity with infection intensity suggested that the coinfection model of May and Nowak (1995) may be appropriate, where each parasite is able to infect and multiply independent of others within the same host.     

Priority effects within coinfected hosts can drive unexpected population‐scale patterns of parasite prevalence

Organisms are frequently coinfected by multiple parasite strains and species, and interactions between parasites within hosts are known to influence parasite prevalence and diversity, as well as epidemic timing. Importantly, interactions between coinfecting parasites can be affected by the order in which they infect hosts (i.e. within‐host priority effects). In this study, we use a single‐host, two‐pathogen, SI model with environmental transmission to explore how within‐host priority effects scale up to alter host population‐scale infection patterns. Specifically, we ask how parasite prevalence changes in the presence of different types of priority effects. We consider two scenarios without priority effects and four scenarios with priority effects where there is either an advantage or a disadvantage to being the first to infect in a coinfected host. Models without priority effects always predict negative relationships between the prevalences of both parasites. In contrast, models with priority effects can yield unimodal prevalence relationships where the prevalence of a focal parasite is minimized or maximized at intermediate prevalences of a coinfecting parasite. The mechanism behind this pattern is that as the prevalence of the coinfecting parasite increases, most infections of the focal parasite change from occurring as solo infections, to first arrival coinfections, to second arrival coinfections. The corresponding changes in parasite fitness as the focal parasite moves from one infection class to another then map to changes in focal parasite prevalence. Further, we found that even when parasites interact negatively within a host, they still can have positive prevalence relationships at the population scale. These results suggest that within‐host priority effects can change host population‐scale infection patterns in systematic (and initially counterintuitive) ways, and that taking them into account may improve disease forecasting in coinfected populations.     

Importance of infection of haemosporidia blood parasites during different life history stages for long‐term reproductive fitness of collared flycatchers

The interaction between birds and haemosporidia blood parasites is a well‐used system in the study of parasite biology. However, where, when and how parasites are transmitted is often unclear and defining parasite transmission dynamics is essential because of how they influence parasite‐mediated costs to the host. In this study, we used cross‐sectional and longitudinal data taken from a collared flycatcher Ficedula albicollis population to investigate the temporal dynamics of haemosporidia parasite infection and parasite‐mediated costs to host fitness. We investigated host–parasite interactions starting at the nestling stage of the bird's life‐cycle and then followed their progress over three breeding attempts to quantify their fitness – measured as the number of offspring they produced that recruited back into the breeding population. We found that the majority of haemosporidia blood parasite infections occurred within the first year of life and that the most common parasite lineages that infected the breeding population also infected juvenile birds in the natal environment. Moreover, our findings suggest that collared flycatcher nestlings in poorer condition could be at a higher risk of haemosporidia blood parasite infection. In this study, only female and not male bird fitness was adversely affected by parasite infection and the cost of infection on female fitness depended on the timing of transmission. In conclusion, our study indicates that in collared flycatchers, early‐life is potentially important for many of the interactions with haemosporidia parasite lineages, and evidence of parasite‐mediated costs to fitness suggest that these parasites may have influenced the host population dynamics.     

Parasite co-infection and interaction as drivers of host heterogeneity

We examined the hypothesis that the interaction between concomitant infecting parasites modifies host susceptibility, parasite intensity and the pattern of parasite distribution within the host population. We used a 26 year time series of three common parasites in a natural population of rabbits: two gastrointestinal nematodes (Trichostrongylus retortaeformis and Graphidium strigosum) and the immunosuppressive myxoma virus. The frequency distribution of nematodes in the host population and the relationship between host age and nematode intensity were explored in rabbits with either single or dual nematode infections and rabbits infected with the nematodes and myxoma virus. The aggregation of T. retortaeformis and G. strigosum among the rabbits varied with the nature of the co-infection both in male and female hosts. The two nematodes exhibited different age-intensity profiles: G. strigosum intensity increased exponentially with host age while T. retortaeformis intensity exhibited a convex shape. The presence of a secondary infection did not change the age-intensity profile for G. strigosum but for T. retortaeformis co-infection (either both nematodes or myxoma-nematodes) resulted in significantly greater intensities in adult hosts. Results suggest that multi-species infections contributed to aggregation of parasites in the host population and to seasonal variation in intensity, but also enhanced differences in parasitism between sexes. This effect was apparent for T. retortaeformis, which appears to elicit a strong acquired immune response but not for G. strigosum which does not produce any evident immune reaction. We concluded that concomitant infections mediated by host immunity are important in modifying host susceptibility and influencing heterogeneity amongst individual hosts.