Avian host composition,local speciation and dispersal drive the regional assembly of avian malaria parasites in South American birds

Identifying the ecological factors that shape parasite distributions remains a central goal in disease ecology. These factors include dispersal capability, environmental filters and geographic distance. Using 520 haemosporidian parasite genetic lineages recovered from 7,534 birds sampled across tropical and temperate South America, we tested (a) the latitudinal diversity gradient hypothesis and (b) the distance–decay relationship (decreasing proportion of shared species between communities with increasing geographic distance) for this host–parasite system. We then inferred the biogeographic processes influencing the diversity and distributions of this cosmopolitan group of parasites across South America. We found support for a latitudinal gradient in diversity for avian haemosporidian parasites, potentially mediated through higher avian host diversity towards the equator. Parasite similarity was correlated with climate similarity, geographic distance and host composition. Local diversification in Amazonian lineages followed by dispersal was the most frequent biogeographic events reconstructed for haemosporidian parasites. Combining macroecological patterns and biogeographic processes, our study reveals that haemosporidian parasites are capable of circumventing geographic barriers and dispersing across biomes, although constrained by environmental filtering. The contemporary diversity and distributions of haemosporidian parasites are mainly driven by historical (speciation) and ecological (dispersal) processes, whereas the parasite community assembly is largely governed by host composition and to a lesser extent by environmental conditions.     

Climate,host phylogeny and the connectivity of host communities govern regional parasite assembly

Aim

Identifying barriers that govern parasite community assembly and parasite invasion risk is critical to understand how shifting host ranges impact disease emergence. We studied regional variation in the phylogenetic compositions of bird species and their blood parasites (Plasmodium and Haemoproteus spp.) to identify barriers that shape parasite community assembly.

Location

Australasia and Oceania.

Methods

We used a data set of parasite infections from >10,000 host individuals sampled across 29 bioregions. Hierarchical models and matrix regressions were used to assess the relative influences of interspecies (host community connectivity and local phylogenetic distinctiveness), climate and geographic barriers on parasite local distinctiveness and composition.

Results

Parasites were more locally distinct (co‐occurred with distantly related parasites) when infecting locally distinct hosts, but less distinct (co‐occurred with closely related parasites) in areas with increased host diversity and community connectivity (a proxy for parasite dispersal potential). Turnover and the phylogenetic symmetry of parasite communities were jointly driven by host turnover, climate similarity and geographic distance.

Main conclusions

Interspecies barriers linked to host phylogeny and dispersal shape parasite assembly, perhaps by limiting parasite establishment or local diversification. Infecting hosts that co‐occur with few related species decreases a parasite's likelihood of encountering related competitors, perhaps increasing invasion potential but decreasing diversification opportunity. While climate partially constrains parasite distributions, future host range expansions that spread distinct parasites and diminish barriers to host shifting will likely be key drivers of parasite invasions.     

Parasite establishment and host extinction in model communities

Studies of host–parasite dynamics usually consider one, or at most two, host species, neglecting the possible effects of other species on the focal hosts and vice versa. To explore the interaction of community structure with host–parasite dynamics, we model the invasion of stable communities of varying size by a parasite. The communities are generated with random interaction coefficients and connectance 0.5. Each community is invaded by parasites with different values of virulence (disease-induced host mortality rate), specificity and transmission rate. The result of each invasion is determined by numerically simulating the dynamics of the community. We classify the outcomes by whether the parasite successfully establishes in the focal host population(s), and, if so, by the proportion of host and non-host species that go extinct as a result of the parasite's introduction. We discuss how the structure of the community and the interaction between hosts and other species affect several important processes of disease ecology: the density threshold for parasite invasion, extinction cascades caused by the parasite, and the frequency of extinctions of hosts and non-hosts. In our simulated communities, non-host species went extinct more frequently than hosts, suggesting the importance of the community context of disease. In some cases, the parasite's invasion induced regular population cycles in the previously stable community.     

Host ecology and life-history traits associated with blood parasite species richness in birds

Identifying host traits associated with the number of different parasite species or strains harboured by a particular host species can have important implications for understanding the impact of parasitism on hosts. We investigated associations between host ecology and life history, and parasite richness and prevalence of the four major avian blood parasite genera. We used an extensive data on blood parasite infections and host ecology in 263 bird species from the Western Palearctic, combining species-specific data with a comparative approach to control for similarity in phenotype among host species due to the effects of common phylogenetic descent. Adult survival rate negatively correlated with the number of parasite species infecting a host species when controlling for similarity due to common descent and body mass. In addition, the prevalence of Haemoproteus, Plasmodium and Leucocytozoon was higher in species harbouring a richer parasite assemblage. These results suggest that the impact on host fitness caused by avian haematozoa may be underestimated in natural populations if the exacerbated virulence associated with exposure to multiple parasites is not taken into account.     

Diversity and host assemblage of avian haemosporidians in different terrestrial ecoregions of Peru

Characterizing the diversity and structure of host–parasite communities is crucial to understanding their eco-evolutionary dynamics. Malaria and related haemosporidian parasites are responsible for fitness loss and mortality in bird species worldwide. However, despite exhibiting the greatest ornithological biodiversity, avian haemosporidians from Neotropical regions are quite unexplored. Here, we analyze the genetic diversity of bird haemosporidian parasites (Plasmodium and Haemoproteus) in 1,336 individuals belonging to 206 bird species to explore for differences in diversity of parasite lineages and bird species across 5 well-differentiated Peruvian ecoregions. We detected 70 different haemosporidian lineages infecting 74 bird species. We showed that 25 out of the 70 haplotypes had not been previously recorded. Moreover, we also identified 81 new host–parasite interactions representing new host records for these haemosporidian parasites. Our outcomes revealed that the effective diversity (as well as the richness, abundance, and Shannon–Weaver index) for both birds and parasite lineages was higher in Amazon basin ecoregions. Furthermore, we also showed that ecoregions with greater diversity of bird species also had high parasite richness, hence suggesting that host community is crucial in explaining parasite richness. Generalist parasites were found in ecoregions with lower bird diversity, implying that the abundance and richness of hosts may shape the exploitation strategy followed by haemosporidian parasites. These outcomes reveal that Neotropical region is a major reservoir of unidentified haemosporidian lineages. Further studies analyzing host distribution and specificity of these parasites in the tropics will provide important knowledge about phylogenetic relationships, phylogeography, and patterns of evolution and distribution of haemosporidian parasites.     

Metabarcoding of eukaryotic parasite communities describes diverse parasite assemblages spanning the primate phylogeny

Despite their ubiquity, in most cases little is known about the impact of eukaryotic parasites on their mammalian hosts. Comparative approaches provide a powerful method to investigate the impact of parasites on host ecology and evolution, though two issues are critical for such efforts: controlling for variation in methods of identifying parasites and incorporating heterogeneity in sampling effort across host species. To address these issues, there is a need for standardized methods to catalogue eukaryotic parasite diversity across broad phylogenetic host ranges. We demonstrate the feasibility of a metabarcoding approach for describing parasite communities by analysing faecal samples from 11 nonhuman primate species representing divergent lineages of the primate phylogeny and the full range of sampling effort (i.e. from no parasites reported in the literature to the best‐studied primates). We detected a number of parasite families and regardless of prior sampling effort, metabarcoding of only ten faecal samples identified parasite families previously undescribed in each host (x? = 8.5 new families per species). We found more overlap between parasite families detected with metabarcoding and published literature when more research effort—measured as the number of publications—had been conducted on the host species' parasites. More closely related primates and those from the same continent had more similar parasite communities, highlighting the biological relevance of sampling even a small number of hosts. Collectively, results demonstrate that metabarcoding methods are sensitive and powerful enough to standardize studies of eukaryotic parasite communities across host species, providing essential new tools for macroecological studies of parasitism.     

Host behaviour–parasite feedback: an essential link between animal behaviour and disease ecology

Animal behaviour and the ecology and evolution of parasites are inextricably linked. For this reason, animal behaviourists and disease ecologists have been interested in the intersection of their respective fields for decades. Despite this interest, most research at the behaviour–disease interface focuses either on how host behaviour affects parasites or how parasites affect behaviour, with little overlap between the two. Yet, the majority of interactions between hosts and parasites are probably reciprocal, such that host behaviour feeds back on parasites and vice versa. Explicitly considering these feedbacks is essential for understanding the complex connections between animal behaviour and parasite ecology and evolution. To illustrate this point, we discuss how host behaviour–parasite feedbacks might operate and explore the consequences of feedback for studies of animal behaviour and parasites. For example, ignoring the feedback of host social structure on parasite dynamics can limit the accuracy of predictions about parasite spread. Likewise, considering feedback in studies of parasites and animal personalities may provide unique insight about the maintenance of variation in personality types. Finally, applying the feedback concept to links between host behaviour and beneficial, rather than pathogenic, microbes may shed new light on transitions between mutualism and parasitism. More generally, accounting for host behaviour–parasite feedbacks can help identify critical gaps in our understanding of how key host behaviours and parasite traits evolve and are maintained.     

Manipulation of host-resource dynamics impacts transmission of trophic parasites

Many complex life cycle parasites rely on predator–prey interactions for transmission, whereby definitive hosts become infected via the consumption of an infected intermediate host. As such, these trophic parasites are embedded in the larger community food web. We postulated that exposure to infection and, hence, parasite transmission are inherently linked to host foraging ecology, and that perturbation of the host-resource dynamic will impact parasite transmission dynamics. We employed a field manipulation experiment in which natural populations of the eastern chipmunk (Tamias striatus) were provisioned with a readily available food resource in clumped or uniform spatial distributions. Using replicated longitudinal capture-mark-recapture techniques, replicated supplemented and unsupplemented control sites were monitored before and after treatment for changes in infection levels with three gastro-intestinal helminth parasites. We predicted that definitive hosts subject to food supplementation would experience lower rates of exposure to infective intermediate hosts, presumably because they shifted their diet away from the intermediate host towards the more readily available resource (sunflower seeds). As predicted, prevalence of infection by the trophically transmitted parasite decreased in response to supplemental food treatment, but no such change in infection prevalence was detected for the two directly transmitted parasites in the system. The fact that food supplementation only had an impact on the transmission of the trophically transmitted parasite, and not the directly transmitted parasites, supports our hypothesis that host foraging ecology directly affects exposure to parasites that rely on the ingestion of intermediate hosts for transmission. We concluded that the relative availability of different food resources has important consequences for the transmission of parasites and, more specifically, parasites that are embedded in the food web. The broader implications of these findings for food web dynamics and disease ecology are discussed.     

The functional importance of parasites in animal communities: many roles at many levels?

Past research on parasites and community ecology has focussed on two distinct levels of the overall community. First, it has been shown that parasites can have a role in structuring host communities. They can have differential effects on the different hosts that they exploit, they can directly debilitate a host that itself is a key structuring force in the community, or they can indirectly alter the phenotype of their host and change the importance of the host for the community. Second, certain parasite species can be important in shaping parasite communities. Dominant parasite species can directly compete with other parasite species inside the host and reduce their abundance to some extent, and parasites that alter host phenotype can indirectly make the host more or less suitable for other parasite species. The possibility that a parasite species simultaneously affects the structure of all levels of the overall community, i.e. the parasite community and the community of free-living animals, is never considered. Given the many direct and indirect ways in which a parasite species can modulate the abundance of other species, it is conceivable that some parasite species have functionally important roles in a community, and that their removal would change the relative composition of the whole community. An example from a soft-sediment intertidal community is used to illustrate how the subtle, indirect effects of a parasite species on non-host species can be very important to the structure of the overall community. Future community studies addressing the many potential influences of parasites will no doubt identify other functionally important parasite species that serve to maintain biodiversity.     

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 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.     

Host traits associated with species roles in parasite sharing networks

The community of host species that a parasite infects is often explained by functional traits and phylogeny, predicting that closely related hosts or those with particular traits share more parasites with other hosts. Previous research has examined parasite community similarity by regressing pairwise parasite community dissimilarity between two host species against host phylogenetic distance. However, pairwise approaches cannot target specific host species responsible for disproportionate levels of parasite sharing. To better identify why some host species contribute differentially to parasite diversity patterns, we represent parasite sharing using ecological networks consisting of host species connected by instances of shared parasitism. These networks can help identify host species and traits associated with high levels of parasite sharing that may subsequently identify important hosts for parasite maintenance and transmission within communities. We used global‐scale parasite sharing networks of ungulates, carnivores, and primates to determine if host importance – encapsulated by the network measures degree, closeness, betweenness, and eigenvector centrality – was predictable based on host traits. Our findings suggest that host centrality in parasite sharing networks is a function of host population density and range size, with range size reflecting both species geographic range and the home range of those species. In the full network, host taxonomic family became an important predictor of centrality, suggesting a role for evolutionary relationships between host and parasite species. More broadly, these findings show that trait data predict key properties of ecological networks, thus highlighting a role for species traits in understanding network assembly, stability, and structure.     

Geographic Variation in Social Parasite Pressure Predicts Intraspecific but not Interspecific Aggressive Responses in Hosts of a Slavemaking Ant

Variation in community composition over a species' geographic range leads to divergent selection pressures, resulting in interpopulation variation in trait expression. One of the most pervasive selective forces stems from antagonists such as parasites. Whereas hosts of microparasites developed sophisticated immune systems, social parasites select for behavioural host defences. Here, we investigated the link between parasite pressure exerted by the socially parasitic slavemaking ant Protomognathus americanus and colony‐level aggression in Temnothorax ants from 17 populations. We studied almost the entire geographic range of two host species, including unparasitized populations. As previous studies have demonstrated that host colonies responding highly aggressively towards conspecifics fare better during slavemaker attacks, we predicted higher aggression levels in severely parasitized populations. Indeed, we demonstrate an increase in aggression towards conspecifics with parasite pressure, a pattern that was consistent over the two host species. In contrast to other studies, aggression against the parasite itself did not shift with parasite pressure. This may be explained by an absence of costs of parasite‐specific aggression in parasite‐free populations. The preferred host species T. longispinosus was generally more aggressive; however, the association between parasite pressure and aggression was found for both species, suggesting convergent co‐adaptation. Two potentially confounding factors, colony density and the co‐occurrence of a competing Temnothorax species in the community, could not explain the level of colony aggression in intra‐ and interspecific interactions. Instead, our study points to social parasite pressure as the determining factor shaping antagonistic interactions within, but not between, host species.     

Host sympatry and body size influence parasite straggling rate in a highly connected multihost,multiparasite system

Parasite lineages commonly diverge when host lineages diverge. However, when large clades of hosts and parasites are analyzed, some cases suggest host switching as another major diversification mechanism. The first step in host switching is the appearance of a parasite on an atypical host, or “straggling.” We analyze the conditions associated with straggling events. We use five species of colonially nesting seabirds from the Galapagos Archipelago and two genera of highly specific ectoparasitic lice to examine host switching. We use both genetic and morphological identification of lice, together with measurements of spatial distribution of hosts in mixed breeding colonies, to test: (1) effects of local host community composition on straggling parasite identity; (2) effects of relative host density within a mixed colony on straggling frequency and parasite species identity; and (3) how straggling rates are influenced by the specifics of louse attachment. Finally, we determine whether there is evidence of breeding in cases where straggling adult lice were found, which may indicate a shift from straggling to the initial stages of host switching. We analyzed more than 5,000 parasite individuals and found that only ~1% of lice could be considered stragglers, with ~5% of 436 host individuals having straggling parasites. We found that the presence of the typical host and recipient host in the same locality influenced straggling. Additionally, parasites most likely to be found on alternate hosts are those that are smaller than the typical parasite of that host, implying that the ability of lice to attach to the host might limit host switching. Given that lice generally follow Harrison's rule, with larger parasites on larger hosts, parasites infecting the larger host species are less likely to successfully colonize smaller host species. Moreover, our study supports the general perception that successful colonization of a novel host is extremely rare, as we found only one nymph of a straggling species, which may indicate successful reproduction.     

Host immunization shapes interspecific associations in trematode parasites

1. Individuals of free-living organisms are commonly infected by multiple parasite species. Under such circumstances, positive or negative associations between the species are possible because of direct or indirect interactions, details in parasite transmission ecology and host-mediated factors. One possible mechanism underlying these processes is host immunity, but its role in shaping these associations has rarely been tackled experimentally.
2. In this study, we tested the effect of host immunization on associations between trematode parasites infecting eyes of fish. We first analysed the associations between three species ( Diplostomum spathaceum , Diplostomum gasterostei and Tylodelphys clavata ) in wild hosts, roach ( Rutilus rutilus ) and perch ( Perca fluviatilis ). Second, using rainbow trout ( Oncorhynchus mykiss ) as a model fish species, we experimentally investigated how sequential immunization of the host (i.e. one parasite species infects and immunizes the host first) could affect the associations between two of the species.
3. The results indicated that most of the associations were positive in wild hosts, which supports between-individual variation in host susceptibility, rather than competitive exclusion between the parasite species. However, positive associations were more common in roach than in perch, possibly reflecting differences in ecological conditions of exposure between the host species. The experimental data showed that positive associations between two of the species were eroded by host immunization against one of the parasite species.
4. We conclude that sequential immunization of hosts has a marked effect on interspecific parasite associations and basically can determine if positive associations are detected or not. This implies that correlative results suggesting non-interactive community structure in general may be obscured by the sequence of previous parasite exposure and corresponding dynamics of host immunization.     

Parasites and invasions: a biogeographic examination of parasites and hosts in native and introduced ranges

Aim To use a comparative approach to understand parasite demographic patterns in native versus introduced populations, evaluating the potential roles of host invasion history and parasite life history. Location North American east and west coasts with a focus on San Francisco Bay (SFB). Methods Species richness and prevalence of trematode parasites were examined in the native and introduced ranges of two gastropod host species, Ilyanassa obsoleta and Littorina saxatilis. We divided the native range into the putative source area for introduction and areas to the north and south; we also sampled the overlapping introduced range in SFB. We dissected 14,781 snails from 103 populations and recorded the prevalence and identity of trematode parasites. We compared trematode species richness and prevalence across the hosts’ introduced and native ranges, and evaluated the influence of host availability on observed patterns. Results Relative to the native range, both I. obsoleta and L. saxatilis have escaped (lost) parasites in SFB, and L. saxatilis demonstrated a greater reduction of trematode diversity and infection prevalence than I. obsoleta. This was not due to sampling inequalities between the hosts. Instead, rarefaction curves suggested complete capture of trematode species in native source and SFB subregions, except for L. saxatilis in SFB, where infection was extremely rare. For I. obsoleta, infection prevalence of trematodes using fish definitive hosts was significantly lower in SFB compared to the native range, unlike those using bird hosts. Host availability partly explained the presence of introduced trematodes in SFB. Main conclusions Differential losses of parasite richness and prevalence for the two gastropod host species in their introduced range is probably the result of several mechanistic factors: time since introduction, propagule pressure, vector of introduction, and host availability. Moreover, the recent occurrence of L. saxatilis’ invasion and its active introduction vector suggest that its parasite diversity and distribution will probably increase over time. Our study suggests that host invasion history and parasite life history play key roles in the extent and diversity of trematodes transferred to introduced populations. Our results also provide vital information for understanding community‐level influences of parasite introductions, as well as for disease ecology in general.     

Parasite predators exhibit a rapid numerical response to increased parasite abundance and reduce transmission to hosts

Predators of parasites have recently gained attention as important parts of food webs and ecosystems. In aquatic systems, many taxa consume free‐living stages of parasites, and can thus reduce parasite transmission to hosts. However, the importance of the functional and numerical responses of parasite predators to disease dynamics is not well understood. We collected host–parasite–predator cooccurrence data from the field, and then experimentally manipulated predator abundance, parasite abundance, and the presence of alternative prey to determine the consequences for parasite transmission. The parasite predator of interest was a ubiquitous symbiotic oligochaete of mollusks, Chaetogaster limnaei limnaei, which inhabits host shells and consumes larval trematode parasites. Predators exhibited a rapid numerical response, where predator populations increased or decreased by as much as 60% in just 5 days, depending on the parasite:predator ratio. Furthermore, snail infection decreased substantially with increasing parasite predator densities, where the highest predator densities reduced infection by up to 89%. Predators of parasites can play an important role in regulating parasite transmission, even when infection risk is high, and especially when predators can rapidly respond numerically to resource pulses. We suggest that these types of interactions might have cascading effects on entire disease systems, and emphasize the importance of considering disease dynamics at the community level.     

Do threatened hosts have fewer parasites? A comparative study in primates

1. Parasites and infectious diseases have become a major concern in conservation biology, in part because they can trigger or accelerate species or population declines. Focusing on primates as a well-studied host clade, we tested whether the species richness and prevalence of parasites differed between threatened and non-threatened host species. 2. We collated data on 386 species of parasites (including viruses, bacteria, protozoa, helminths and arthropods) reported to infect wild populations of 36 threatened and 81 non-threatened primate species. Analyses controlled for uneven sampling effort and host phylogeny. 3. Results showed that total parasite species richness was lower among threatened primates, supporting the prediction that small, isolated host populations harbour fewer parasite species. This trend was consistent across three major parasite groups found in primates (helminths, protozoa and viruses). Counter to our predictions, patterns of parasite species richness were independent of parasite transmission mode and the degree of host specificity. 4. We also examined the prevalence of selected parasite genera among primate sister-taxa that differed in their ranked threat categories, but found no significant differences in prevalence between threatened and non-threatened hosts. 5. This study is the first to demonstrate differences in parasite richness relative to host threat status. Results indicate that human activities and host characteristics that increase the extinction risk of wild animal species may lead simultaneously to the loss of parasites. Lower average parasite richness in threatened host taxa also points to the need for a better understanding of the cascading effects of host biodiversity loss for affiliated parasite species.     

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.