The Evolution of Parasite Recognition Genes in the Innate Immune System: Purifying Selection on <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> Peptidoglycan Recognition Proteins

Genes involved in the recognition of parasites by the acquired immune system are often subject to intense selection pressures. In some cases, selection to recognize a diverse range of parasites has resulted in high levels of polymorphism, while elsewhere the protein sequence has changed rapidly under directional selection. We tested whether parasite recognition genes in the innate immune system show similar patterns of evolution. We sequenced seven peptidoglycan recognition protein genes (PGRPs) from 12 lines of Drosophila melanogaster and one line of D. simulans and used a variety of tests to determine whether the observed mutations were selectively neutral. We were unable to detect either balancing or directional selection. This suggests that the molecular cues used by insects to detect parasites are highly conserved and probably under strong functional constraints which prevent their evolving to evade the host immune response. Therefore, interactions between these genes are unlikely to be the focus of host–parasite coevolution, at least in Drosophila. We also found evidence of gene conversion occurring between two genes, PGRP-SC1A and PGRP-SC1B.     

Long-Term and Short-Term Evolutionary Impacts of Transposable Elements on Drosophila

Transposable elements (TEs) are considered to be genomic parasites and their interactions with their hosts have been likened to the coevolution between host and other nongenomic, horizontally transferred pathogens. TE families, however, are vertically inherited as integral segments of the nuclear genome. This transmission strategy has been suggested to weaken the selective benefits of host alleles repressing the transposition of specific TE variants. On the other hand, the elevated rates of TE transposition and high incidences of deleterious mutations observed during the rare cases of horizontal transfers of TE families between species could create at least a transient process analogous to the influence of horizontally transmitted pathogens. Here, we formally address this analogy, using empirical and theoretical analysis to specify the mechanism of how host–TE interactions may drive the evolution of host genes. We found that host TE-interacting genes actually have more pervasive evidence of adaptive evolution than immunity genes that interact with nongenomic pathogens in Drosophila. Yet, both our theoretical modeling and empirical observations comparing Drosophila melanogaster populations before and after the horizontal transfer of P elements, which invaded D. melanogaster early last century, demonstrated that horizontally transferred TEs have only a limited influence on host TE-interacting genes. We propose that the more prevalent and constant interaction with multiple vertically transmitted TE families may instead be the main force driving the fast evolution of TE-interacting genes, which is fundamentally different from the gene-for-gene interaction of host–pathogen coevolution.     

Balancing Selection for Pathogen Resistance Reveals an Intercontinental Signature of Red Queen Coevolution

The link between long-term host–parasite coevolution and genetic diversity is key to understanding genetic epidemiology and the evolution of resistance. The model of Red Queen host–parasite coevolution posits that high genetic diversity is maintained when rare host resistance variants have a selective advantage, which is believed to be the mechanistic basis for the extraordinarily high levels of diversity at disease-related genes such as the major histocompatibility complex in jawed vertebrates and R-genes in plants. The parasites that drive long-term coevolution are, however, often elusive. Here we present evidence for long-term balancing selection at the phenotypic (variation in resistance) and genomic (resistance locus) level in a particular host–parasite system: the planktonic crustacean Daphnia magna and the bacterium Pasteuria ramosa. The host shows widespread polymorphisms for pathogen resistance regardless of geographic distance, even though there is a clear genome-wide pattern of isolation by distance at other sites. In the genomic region of a previously identified resistance supergene, we observed consistent molecular signals of balancing selection, including higher genetic diversity, older coalescence times, and lower differentiation between populations, which set this region apart from the rest of the genome. We propose that specific long-term coevolution by negative-frequency-dependent selection drives this elevated diversity at the host''s resistance loci on an intercontinental scale and provide an example of a direct link between the host’s resistance to a virulent pathogen and the large-scale diversity of its underlying genes.     

Major histocompatibility complex (MHC) heterozygote superiority to natural multi-parasite infections in the water vole (Arvicola terrestris)

The fundamental role of the major histocompatibility complex (MHC) in immune recognition has led to a general consensus that the characteristically high levels of functional polymorphism at MHC genes is maintained by balancing selection operating through host–parasite coevolution. However, the actual mechanism by which selection operates is unclear. Two hypotheses have been proposed: overdominance (or heterozygote superiority) and negative frequency-dependent selection. Evidence for these hypotheses was evaluated by examining MHC–parasite relationships in an island population of water voles (Arvicola terrestris). Generalized linear mixed models were used to examine whether individual variation at an MHC class II DRB locus explained variation in the individual burdens of five different parasites. MHC genotype explained a significant amount of variation in the burden of gamasid mites, fleas (Megabothris walkeri) and nymphs of sheep ticks (Ixodes ricinus). Additionally, MHC heterozygotes were simultaneously co-infected by fewer parasite types than homozygotes. In each case where an MHC-dependent effect on parasite burden was resolved, the heterozygote genotype was associated with fewer parasites, and the heterozygote outperformed each homozygote in two of three cases, suggesting an overall superiority against parasitism for MHC heterozygote genotypes. This is the first demonstration of MHC heterozygote superiority against multiple parasites in a natural population, a mechanism that could help maintain high levels of functional MHC genetic diversity in natural populations.     

Molecular Evolution of Daphnia Immunity Genes: Polymorphism in a Gram-Negative Binding Protein Gene and an α-2-Macroglobulin Gene

Studies of DNA polymorphism have shown that some immune system genes of mammals and plants are exceptionally diverse, indicating that coevolution between these taxa and their parasites mediates positive selective sweeps and/or balancing selection. The genes of the arthropod immune system remain comparatively unstudied. We isolated two putative immune system genes from the cladoceran crustacean Daphnia and examined DNA sequence diversity. For one gene, encoding a putative gram-negative binding protein, we found evidence of only purifying selection, indicating that this gene is under strong functional constraint and that selection acts to eliminate amino acid variation. For another gene, encoding a putative -2-macroglobulin, we found evidence of positive selection, indicating the possible involvement of this gene in a host–parasite arms race. We discuss the assumed function of these genes and offer speculation regarding which components of the arthropod immune system might experience diversifying adaptive evolution.[Reviewing Editor: Dr. Martin Kreitman]     

Parasite and host assemblages: embracing the reality will improve our knowledge of parasite transmission and virulence

Interactions involving several parasite species (multi-parasitized hosts) or several host species (multi-host parasites) are the rule in nature. Only a few studies have investigated these realistic, but complex, situations from an evolutionary perspective. Consequently, their impact on the evolution of parasite virulence and transmission remains poorly understood. The mechanisms by which multiple infections may influence virulence and transmission include the dynamics of intrahost competition, mediation by the host immune system and an increase in parasite genetic recombination. Theoretical investigations have yet to be conducted to determine which of these mechanisms are likely to be key factors in the evolution of virulence and transmission. In contrast, the relationship between multi-host parasites and parasite virulence and transmission has seen some theoretical investigation. The key factors in these models are the trade-off between virulence across different host species, variation in host species quality and patterns of transmission. The empirical studies on multi-host parasites suggest that interspecies transmission plays a central role in the evolution of virulence, but as yet no complete picture of the phenomena involved is available. Ultimately, determining how complex host–parasite interactions impact the evolution of host–parasite relationships will require the development of cross-disciplinary studies linking the ecology of quantitative networks with the evolution of virulence.     

Protein-poor diet reduces host-specific immune gene expression in Bombus terrestris

Parasites infect hosts non-randomly as genotypes of hosts vary in susceptibility to the same genotypes of parasites, but this specificity may be modulated by environmental factors such as nutrition. Nutrition plays an important role for any physiological investment. As immune responses are costly, resource limitation should negatively affect immunity through trade-offs with other physiological requirements. Consequently, nutritional limitation should diminish immune capacity in general, but does it also dampen differences among hosts? We investigated the effect of short-term pollen deprivation on the immune responses of our model host Bombus terrestris when infected with the highly prevalent natural parasite Crithidia bombi. Bumblebees deprived of pollen, their protein source, show reduced immune responses to infection. They failed to upregulate a number of genes, including antimicrobial peptides, in response to infection. In particular, they also showed less specific immune expression patterns across individuals and colonies. These findings provide evidence for how immune responses on the individual-level vary with important elements of the environment and illustrate how nutrition can functionally alter not only general resistance, but also alter the pattern of specific host–parasite interactions.     

Genetic architecture of resistance in Daphnia hosts against two species of host-specific parasites

Understanding the genetic architecture of host resistance is key for understanding the evolution of host–parasite interactions. Evolutionary models often assume simple genetics based on few loci and strong epistasis. It is unknown, however, whether these assumptions apply to natural populations. Using a quantitative trait loci (QTL) approach, we explore the genetic architecture of resistance in the crustacean Daphnia magna to two of its natural parasites: the horizontally transmitted bacterium Pasteuria ramosa and the horizontally and vertically transmitted microsporidium Hamiltosporidium tvaerminnensis. These two systems have become models for studies on the evolution of host–parasite interactions. In the QTL panel used here, Daphnia''s resistance to P. ramosa is controlled by a single major QTL (which explains 50% of the observed variation). Resistance to H. tvaerminnensis horizontal infections shows a signature of a quantitative trait based in multiple loci with weak epistatic interactions (together explaining 38% variation). Resistance to H. tvaerminnensis vertical infections, however, shows only one QTL (explaining 13.5% variance) that colocalizes with one of the QTLs for horizontal infections. QTLs for resistance to Pasteuria and Hamiltosporidium do not colocalize. We conclude that the genetics of resistance in D. magna are drastically different for these two parasites. Furthermore, we infer that based on these and earlier results, the mechanisms of coevolution differ strongly for the two host–parasite systems. Only the Pasteuria–Daphnia system is expected to follow the negative frequency-dependent selection (Red Queen) model. How coevolution works in the Hamiltosporidium–Daphnia system remains unclear.     

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.     

Lobesia botrana Larvae Develop Faster in the Presence of Parasitoids

To combat parasitism hosts often rely on their immune system, which is the last line of defense. However, the immune system may not always be effective, and other non-immunological defenses might be favored to reduce the cost of parasite infection. Here we report that larvae of the moth Lobesia botrana can rapidly accelerate their development and reach maturity earlier in response to cues perceived at a distance from parasitoids. Such a phenotypically plastic life history shift, induced by the perception of deadly enemies in the environment, is likely to be an adaptive defensive strategy to prevent parasitoid attack, and has important implications in host–parasite dynamics.     

Noninvasively measured immune responses reflect current parasite infections in a wild carnivore and are linked to longevity

Host immune defenses are important components of host–parasite interactions that affect the outcome of infection and may have fitness consequences for hosts when increased allocation of resources to immune responses undermines other essential life processes. Research on host–parasite interactions in large free‐ranging wild mammals is currently hampered by a lack of verified noninvasive assays. We successfully adapted existing assays to measure innate and adaptive immune responses produced by the gastrointestinal mucosa in spotted hyena (Crocuta crocuta) feces, including enzyme‐linked immunosorbent assays (ELISAs), to quantify fecal immunoglobulins (total IgA, total IgG) and total fecal O‐linked oligosaccharides (mucin). We investigated the effect of infection load by an energetically costly hookworm (Ancylostoma), parasite richness, host age, sex, year of sampling, and clan membership on immune responses and asked whether high investment in immune responses during early life affects longevity in individually known spotted hyenas in the Serengeti National Park, Tanzania. Fecal concentrations of IgA, IgG, and mucin increased with Ancylostoma egg load and were higher in juveniles than in adults. Females had higher mucin concentrations than males. Juvenile females had higher IgG concentrations than juvenile males, whereas adult females had lower IgG concentrations than adult males. High IgA concentrations during the first year of life were linked to reduced longevity after controlling for age at sampling and Ancylostoma egg load. Our study demonstrates that the use of noninvasive methods can increase knowledge on the complex relationship between gastrointestinal parasites and host local immune responses in wild large mammals and reveal fitness‐relevant effects of these responses.     

Genomic Signatures of Selective Pressures and Introgression from Archaic Hominins at Human Innate Immunity Genes

Human genes governing innate immunity provide a valuable tool for the study of the selective pressure imposed by microorganisms on host genomes. A comprehensive, genome-wide study of how selective constraints and adaptations have driven the evolution of innate immunity genes is missing. Using full-genome sequence variation from the 1000 Genomes Project, we first show that innate immunity genes have globally evolved under stronger purifying selection than the remainder of protein-coding genes. We identify a gene set under the strongest selective constraints, mutations in which are likely to predispose individuals to life-threatening disease, as illustrated by STAT1 and TRAF3. We then evaluate the occurrence of local adaptation and detect 57 high-scoring signals of positive selection at innate immunity genes, variation in which has been associated with susceptibility to common infectious or autoimmune diseases. Furthermore, we show that most adaptations targeting coding variation have occurred in the last 6,000–13,000 years, the period at which populations shifted from hunting and gathering to farming. Finally, we show that innate immunity genes present higher Neandertal introgression than the remainder of the coding genome. Notably, among the genes presenting the highest Neandertal ancestry, we find the TLR6-TLR1-TLR10 cluster, which also contains functional adaptive variation in Europeans. This study identifies highly constrained genes that fulfill essential, non-redundant functions in host survival and reveals others that are more permissive to change—containing variation acquired from archaic hominins or adaptive variants in specific populations—improving our understanding of the relative biological importance of innate immunity pathways in natural conditions.     

Immune selection suppresses the emergence of drug resistance in malaria parasites but facilitates its spread

Although drug resistance in Plasmodium falciparum typically evolves in regions of low transmission, resistance spreads readily following introduction to regions with a heavier disease burden. This suggests that the origin and the spread of resistance are governed by different processes, and that high transmission intensity specifically impedes the origin. Factors associated with high transmission, such as highly immune hosts and competition within genetically diverse infections, are associated with suppression of resistant lineages within hosts. However, interactions between these factors have rarely been investigated and the specific relationship between adaptive immunity and selection for resistance has not been explored. Here, we developed a multiscale, agent-based model of Plasmodium parasites, hosts, and vectors to examine how host and parasite dynamics shape the evolution of resistance in populations with different transmission intensities. We found that selection for antigenic novelty (“immune selection”) suppressed the evolution of resistance in high transmission settings. We show that high levels of population immunity increased the strength of immune selection relative to selection for resistance. As a result, immune selection delayed the evolution of resistance in high transmission populations by allowing novel, sensitive lineages to remain in circulation at the expense of the spread of a resistant lineage.In contrast, in low transmission settings, we observed that resistant strains were able to sweep to high population prevalence without interference. Additionally, we found that the relationship between immune selection and resistance changed when resistance was widespread. Once resistance was common enough to be found on many antigenic backgrounds, immune selection stably maintained resistant parasites in the population by allowing them to proliferate, even in untreated hosts, when resistance was linked to a novel epitope. Our results suggest that immune selection plays a role in the global pattern of resistance evolution.     

Faster-X Evolution of Gene Expression in Drosophila

DNA sequences on X chromosomes often have a faster rate of evolution when compared to similar loci on the autosomes, and well articulated models provide reasons why the X-linked mode of inheritance may be responsible for the faster evolution of X-linked genes. We analyzed microarray and RNA–seq data collected from females and males of six Drosophila species and found that the expression levels of X-linked genes also diverge faster than autosomal gene expression, similar to the “faster-X” effect often observed in DNA sequence evolution. Faster-X evolution of gene expression was recently described in mammals, but it was limited to the evolutionary lineages shortly following the creation of the therian X chromosome. In contrast, we detect a faster-X effect along both deep lineages and those on the tips of the Drosophila phylogeny. In Drosophila males, the dosage compensation complex (DCC) binds the X chromosome, creating a unique chromatin environment that promotes the hyper-expression of X-linked genes. We find that DCC binding, chromatin environment, and breadth of expression are all predictive of the rate of gene expression evolution. In addition, estimates of the intraspecific genetic polymorphism underlying gene expression variation suggest that X-linked expression levels are not under relaxed selective constraints. We therefore hypothesize that the faster-X evolution of gene expression is the result of the adaptive fixation of beneficial mutations at X-linked loci that change expression level in cis. This adaptive faster-X evolution of gene expression is limited to genes that are narrowly expressed in a single tissue, suggesting that relaxed pleiotropic constraints permit a faster response to selection. Finally, we present a conceptional framework to explain faster-X expression evolution, and we use this framework to examine differences in the faster-X effect between Drosophila and mammals.     

Cross‐continental comparison of parasite communities in a wide‐ranging carnivore suggests associations with prey diversity and host density

1. Parasites are integral to ecosystem functioning yet often overlooked. Improved understanding of host–parasite associations is important, particularly for wide‐ranging species for which host range shifts and climate change could alter host–parasite interactions and their effects on ecosystem function.
2. Among the most widely distributed mammals with diverse diets, gray wolves (Canis lupus) host parasites that are transmitted among canids and via prey species. Wolf–parasite associations may therefore influence the population dynamics and ecological functions of both wolves and their prey. Our goal was to identify large‐scale processes that shape host–parasite interactions across populations, with the wolf as a model organism.
3. By compiling data from various studies, we examined the fecal prevalence of gastrointestinal parasites in six wolf populations from two continents in relation to wolf density, diet diversity, and other ecological conditions.
4. As expected, we found that the fecal prevalence of parasites transmitted directly to wolves via contact with other canids or their excreta was positively associated with wolf density. Contrary to our expectations, the fecal prevalence of parasites transmitted via prey was negatively associated with prey diversity. We also found that parasite communities reflected landscape characteristics and specific prey items available to wolves.
5. Several parasite taxa identified in this study, including hookworms and coccidian protozoans, can cause morbidity and mortality in canids, especially in pups, or in combination with other stressors. The density–prevalence relationship for parasites with simple life cycles may reflect a regulatory role of gastrointestinal parasites on wolf populations. Our result that fecal prevalence of parasites was lower in wolves with more diverse diets could provide insight into the mechanisms by which biodiversity may regulate disease. A diverse suite of predator–prey interactions could regulate the effects of parasitism on prey populations and mitigate the transmission of infectious agents, including zoonoses, spread via trophic interactions.

     

MHC-I Affects Infection Intensity but Not Infection Status with a Frequent Avian Malaria Parasite in Blue Tits

Host resistance against parasites depends on three aspects: the ability to prevent, control and clear infections. In vertebrates the immune system consists of innate and adaptive immunity. Innate immunity is particularly important for preventing infection and eradicating established infections at an early stage while adaptive immunity is slow, but powerful, and essential for controlling infection intensities and eventually clearing infections. Major Histocompatibility Complex (MHC) molecules are central in adaptive immunity, and studies on parasite resistance and MHC in wild animals have found effects on both infection intensity (parasite load) and infection status (infected or not). It seems MHC can affect both the ability to control infection intensities and the ability to clear infections. However, these two aspects have rarely been considered simultaneously, and their relative importance in natural populations is therefore unclear. Here we investigate if MHC class I genotype affects infection intensity and infection status with a frequent avian malaria infection Haemoproteus majoris in a natural population of blue tits Cyanistes caeruleus. We found a significant negative association between a single MHC allele and infection intensity but no association with infection status. Blue tits that carry a specific MHC allele seem able to suppress H. majoris infection intensity, while we have no evidence that this allele also has an effect on clearance of the H. majoris infection, a result that is in contrast with some previous studies of MHC and avian malaria. A likely explanation could be that the clearance rate of avian malaria parasites differs between avian malaria lineages and/or between avian hosts.     

Within and transgenerational immune priming in an insect to a DNA virus

Invertebrates mount a sophisticated immune response with the potential to exhibit a form of immune memory through ‘priming’. Increased immune protection following early exposure to bacteria has been found both later in life (within generation priming) and in the next generation (transgeneration priming) in a number of invertebrates. However, it is unclear how general immune priming is and whether immune priming occurs in response to different parasites, including viruses. Here, using Plodia interpuctella (Lepidoptera) and its natural DNA virus, Plodia interpunctella granulosis virus, we find evidence for both within generation and transgeneration immune priming. Individuals previously exposed to low doses of virus, as well as the offspring of exposed individuals, are subsequently less susceptible to viral challenge. Relatively little is known about the mechanisms that underpin viral immunity but it is probable that the viral immune response is somewhat different to that of bacteria. We show that immune priming may, however, be a characteristic of both responses, mediated through different mechanisms, suggesting that immune memory may be a general phenomenon of insect immunity. This is important because immune priming may influence both host–parasite population and evolutionary dynamics.     

Lineage-specific positive selection at the merozoite surface protein 1 (msp1) locus of Plasmodium vivax and related simian malaria parasites

Background

The 200 kDa merozoite surface protein 1 (MSP-1) of malaria parasites, a strong vaccine candidate, plays a key role during erythrocyte invasion and is a target of host protective immune response. Plasmodium vivax, the most widespread human malaria parasite, is closely related to parasites that infect Asian Old World monkeys, and has been considered to have become a parasite of man by host switch from a macaque malaria parasite. Several Asian monkey parasites have a range of natural hosts. The same parasite species shows different disease manifestations among host species. This suggests that host immune responses to P. vivax-related malaria parasites greatly differ among host species (albeit other factors). It is thus tempting to invoke that a major immune target parasite protein such as MSP-1 underwent unique evolution, depending on parasite species that exhibit difference in host range and host specificity.

Results

We performed comparative phylogenetic and population genetic analyses of the gene encoding MSP-1 (msp1) from P. vivax and nine P. vivax-related simian malaria parasites. The inferred phylogenetic tree of msp1 significantly differed from that of the mitochondrial genome, with a striking displacement of P. vivax from a position close to P. cynomolgi in the mitochondrial genome tree to an outlier of Asian monkey parasites. Importantly, positive selection was inferred for two ancestral branches, one leading to P. inui and P. hylobati and the other leading to P. vivax, P. fieldi and P. cynomolgi. This ancestral positive selection was estimated to have occurred three to six million years ago, coinciding with the period of radiation of Asian macaques. Comparisons of msp1 polymorphisms between P. vivax, P. inui and P. cynomolgi revealed that while some positively selected amino acid sites or regions are shared by these parasites, amino acid changes greatly differ, suggesting that diversifying selection is acting species-specifically on msp1.

Conclusions

The present results indicate that the msp1 locus of P. vivax and related parasite species has lineage-specific unique evolutionary history with positive selection. P. vivax and related simian malaria parasites offer an interesting system toward understanding host species-dependent adaptive evolution of immune-target surface antigen genes such as msp1.     

Competitive outcome of multiple infections in a behavior‐manipulating virus/wasp interaction

Infections by multiple parasites are common in nature and may impact the evolution of host–parasite interactions. We investigated the existence of multiple infections involving the DNA virus LbFV and the Drosophila parasitoid Leptopilina boulardi. This vertically transmitted virus forces infected females to lay their eggs in already parasitized Drosophila larvae (a behavior called superparasitism), thus favoring its spread through horizontal transmission. Previous theoretical work indicated that the evolution of the level of the manipulation strongly depends on whether infected parasitoids can be re‐infected or not. Here, we describe a strain of LbFV that differs from the reference strain by showing a deletion within the locus used for PCR detection. We used this polymorphism to test for the existence of multiple infections in this system. Viral strains did not differ on their vertical or horizontal transmission rates nor on the way they affect the parasitoid''s phenotype, including their ability to manipulate behavior. Although already infected parasitoids were much less susceptible to new infection than uninfected ones, frequent coinfection was detected. However, following coinfection, competition between viral strains led to the rapid elimination of one strain or the other after a few generations of vertical transmission. We discuss the implications of these results for the evolution of the behavioral manipulation.