In Vivo Microbial Coevolution Favors Host Protection and Plastic Downregulation of Immunity

Microbiota can protect their hosts from infection. The short timescales in which microbes can evolve presents the possibility that “protective microbes” can take-over from the immune system of longer-lived hosts in the coevolutionary race against pathogens. Here, we found that coevolution between a protective bacterium (Enterococcus faecalis) and a virulent pathogen (Staphylococcus aureus) within an animal population (Caenorhabditis elegans) resulted in more disease suppression than when the protective bacterium adapted to uninfected hosts. At the same time, more protective E. faecalis populations became costlier to harbor and altered the expression of 134 host genes. Many of these genes appear to be related to the mechanism of protection, reactive oxygen species production. Crucially, more protective E. faecalis populations downregulated a key immune gene, , known to be effective against S. aureus infection. These results suggest that a microbial line of defense is favored by microbial coevolution and may cause hosts to plastically divest of their own immunity.     

Adaptive population structure shifts in invasive parasitic mites,Varroa destructor

Comparative studies of genetic diversity and population structure can shed light on the ecological and evolutionary factors governing host–parasite interactions. Even though invasive parasites are considered of major biological importance, little is known about their adaptative potential when infesting the new hosts. Here, the genetic diversification of Varroa destructor, a novel parasite of Apis mellifera originating from Asia, was investigated using population genetics to determine how the genetic structure of the parasite changed in distinct European populations of its new host. To do so, mites infesting two categories of hosts in four European regions were compared: (a) adapted hosts surviving through means of natural selection, thereby expected to impose strong selective pressure on the mites, and (b) treated host populations, surviving mite infestations because acaricides are applied, therefore characterized by a relaxed selection imposed by the host on the mites. Significant genetic divergence was found across regions, partially reflecting the invasion pattern of V. destructor throughout Europe and indicating local adaptation of the mite to the host populations. Additionally, varying degrees of genotypic changes were found between mites from adapted and treated colonies. Altogether, these results indicate that V. destructor managed to overcome the genetic bottlenecks following its introduction in Europe and that host‐mediated selection fostered changes in the genetic structure of this mite at diverse geographic scales. These findings highlight the potential of parasites to adapt to their local host populations and confirm that adaptations developed within coevolutionary dynamics are a major determinant of population genetic changes.     

Genome-Wide Association Analysis Identifies a Genetic Basis of Infectivity in a Model Bacterial Pathogen

Knowledge of the genetic architecture of pathogen infectivity and host resistance is essential for a mechanistic understanding of coevolutionary processes, yet the genetic basis of these interacting traits remains unknown for most host–pathogen systems. We used a comparative genomic approach to explore the genetic basis of infectivity in Pasteuria ramosa, a Gram-positive bacterial pathogen of planktonic crustaceans that has been established as a model for studies of Red Queen host–pathogen coevolution. We sequenced the genomes of a geographically, phenotypically, and genetically diverse collection of P. ramosa strains and performed a genome-wide association study to identify genetic correlates of infection phenotype. We found multiple polymorphisms within a single gene, Pcl7, that correlate perfectly with one common and widespread infection phenotype. We then confirmed this perfect association via Sanger sequencing in a large and diverse sample set of P. ramosa clones. Pcl7 codes for a collagen-like protein, a class of adhesion proteins known or suspected to be involved in the infection mechanisms of a number of important bacterial pathogens. Consistent with expectations under Red Queen coevolution, sequence variation of Pcl7 shows evidence of balancing selection, including extraordinarily high diversity and absence of geographic structure. Based on structural homology with a collagen-like protein of Bacillus anthracis, we propose a hypothesis for the structure of Pcl7 and the physical location of the phenotype-associated polymorphisms. Our results offer strong evidence for a gene governing infectivity and provide a molecular basis for further study of Red Queen dynamics in this model host–pathogen system.     

Arabidopsis thaliana as a model for the study of plant–virus co-evolution

Understanding plant–virus coevolution requires wild systems in which there is no human manipulation of either host or virus. To develop such a system, we analysed virus infection in six wild populations of Arabidopsis thaliana in Central Spain. The incidence of five virus species with different life-styles was monitored during four years, and this was analysed in relation to the demography of the host populations. Total virus incidence reached 70 per cent, which suggests a role of virus infection in the population structure and dynamics of the host, under the assumption of a host fitness cost caused by the infection. Maximum incidence occurred at early growth stages, and co-infection with different viruses was frequent, two factors often resulting in increased virulence. Experimental infections under controlled conditions with two isolates of the most prevalent viruses, cauliflower mosaic virus and cucumber mosaic virus, showed that there is genetic variation for virus accumulation, although this depended on the interaction between host and virus genotypes. Comparison of Q_ST-based genetic differentiations between both host populations with F_ST genetic differentiation based on putatively neutral markers suggests different selection dynamics for resistance against different virus species or genotypes. Together, these results are compatible with a hypothesis of plant–virus coevolution.     

Polymorphisms at the innate immune receptor TLR2 are associated with Borrelia infection in a wild rodent population

The discovery of the key role of Toll-like receptors (TLRs) in initiating innate immune responses and modulating adaptive immunity has revolutionized our understanding of vertebrate defence against pathogens. Yet, despite their central role in pathogen recognition and defence initiation, there is little information on how variation in TLRs influences disease susceptibility in natural populations. Here, we assessed the extent of naturally occurring polymorphisms at TLR2 in wild bank voles (Myodes glareolus) and tested for associations between TLR2 variants and infection with Borrelia afzelii, a common tick-transmitted pathogen in rodents and one of the causative agents of human Lyme disease. Bank voles in our population had 15 different TLR2 haplotypes (10 different haplotypes at the amino acid level), which grouped in three well-separated clusters. In a large-scale capture–mark–recapture study, we show that voles carrying TLR2 haplotypes of one particular cluster (TLR2_c2) were almost three times less likely to be Borrelia infected than animals carrying other haplotypes. Moreover, neutrality tests suggested that TLR2 has been under positive selection. This is, to our knowledge, the first demonstration of an association between TLR polymorphism and parasitism in wildlife, and a striking example that genetic variation at innate immune receptors can have a large impact on host resistance.     

Evaluating the within-host fitness effects of mutations fixed during virus adaptation to different ecotypes of a new host

The existence of genetic variation for resistance in host populations is assumed to be essential to the spread of an emerging virus. Models predict that the rate of spread slows down with the increasing frequency and higher diversity of resistance alleles in the host population. We have been using the experimental pathosystem Arabidopsis thaliana—tobacco etch potyvirus (TEV) to explore the interplay between genetic variation in host''s susceptibility and virus diversity. We have recently shown that TEV populations evolving in A. thaliana ecotypes that differ in susceptibility to infection gained within-host fitness, virulence and infectivity in a manner compatible with a gene-for-gene model of host–parasite interactions: hard-to-infect ecotypes were infected by generalist viruses, whereas easy-to-infect ecotypes were infected by every virus. We characterized the genomes of the evolved viruses and found cases of host-driven convergent mutations. To gain further insights in the mechanistic basis of this gene-for-gene model, we have generated all viral mutations individually as well as in specific combinations and tested their within-host fitness effects across ecotypes. Most of these mutations were deleterious or neutral in their local ecotype and only a very reduced number had a host-specific beneficial effect. We conclude that most of the mutations fixed during the evolution experiment were so by drift or by selective sweeps along with the selected driver mutation. In addition, we evaluated the ruggedness of the underlying adaptive fitness landscape and found that mutational effects were mostly multiplicative, with few cases of significant epistasis.     

Sexual antagonism for resistance and tolerance to infection in Drosophila melanogaster

A critical task in evolutionary genetics is to explain the persistence of heritable variation in fitness-related traits such as immunity. Ecological factors can maintain genetic variation in immunity, but less is known about the role of other factors, such as antagonistic pleiotropy, on immunity. Sexually dimorphic immunity—with females often being more immune-competent—may maintain variation in immunity in dioecious populations. Most eco-immunological studies assess host resistance to parasites rather than the host''s ability to maintain fitness during infection (tolerance). Distinguishing between resistance and tolerance is important as they are thought to have markedly different evolutionary and epidemiological outcomes. Few studies have investigated tolerance in animals, and the extent of sexual dimorphism in tolerance is unknown. Using males and females from 50 Drosophila melanogaster genotypes, we investigated possible sources of genetic variation for immunity by assessing both resistance and tolerance to the common bacterial pathogen Pseudomonas aeruginosa. We found evidence of sexual dimorphism and sexual antagonism for resistance and tolerance, and a trade-off between the two traits. Our findings suggest that antagonistic pleiotropy may be a major contributor to variation in immunity, with implications for host–parasite coevolution.     

Effect of spatial connectivity on host resistance in a highly fragmented natural pathosystem

Both theory and experimental evolution studies predict migration to influence the outcome of antagonistic coevolution between hosts and their parasites, with higher migration rates leading to increased diversity and evolutionary potential. Migration rates are expected to vary in spatially structured natural pathosystems, yet how spatial structure generates variation in coevolutionary trajectories across populations occupying the same landscape has not been tested. Here, we studied the effect of spatial connectivity on host evolutionary potential in a natural pathosystem characterized by a stable Plantago lanceolata host network and a highly dynamic Podosphaera plantaginis parasite metapopulation. We designed a large inoculation experiment to test resistance of five isolated and five well‐connected host populations against sympatric and allopatric pathogen strains, over 4 years. Contrary to our expectations, we did not find consistently higher resistance against sympatric pathogen strains in the well‐connected populations. Instead, host local adaptation varied considerably among populations and through time with greater fluctuations observed in the well‐connected populations. Jointly, our results suggest that in populations where pathogens have successfully established, they have the upper hand in the coevolutionary arms race, but hosts may be better able to respond to pathogen‐imposed selection in the well‐connected than in the isolated populations. Hence, the ongoing and extensive fragmentation of natural habitats may increase vulnerability to diseases.     

Local adaptation in the Linum marginale-Melampsora lini host-pathogen interaction

The potential for local adaptation between pathogens and their hosts has generated strong theoretical and empirical interest with evidence both for and against local adaptation reported for a range of systems. We use the Linum marginale-Melampsora lini plant-pathogen system and a hierarchical spatial structure to investigate patterns of local adaptation within a metapopulation characterised by epidemic dynamics and frequent extinction of pathogen populations. Based on large sample sizes and comprehensive cross-inoculation trials, our analyses demonstrate strong local adaptation by Melampsora to its host populations, with this effect being greatest at regional scales, as predicted from the broader spatial scales at which M. lini disperses relative to L. marginale. However, there was no consistent trend for more distant pathogen populations to perform more poorly. Our results further show how the coevolutionary interaction between hosts and pathogens can be influenced by local structure such that resistant hosts select for generally virulent pathogens, while susceptible hosts select for more avirulent pathogens. Empirically, local adaptation has generally been tested in two contrasting ways: (1) pathogen performance on sympatric versus allopatric hosts; and (2) sympatric versus allopatric pathogens on a given host population. In situations where no host population is more resistant or susceptible than others when averaged across pathogen populations (and likewise, no pathogen population is more virulent or avirulent than others), results from these tests should generally be congruent. We argue that this is unlikely to be the case in the metapopulation situations that predominate in natural host-pathogen interactions, thus requiring tests that control simultaneously for variation in plant and pathogen populations.     

Adaptation of tobacco etch potyvirus to a susceptible ecotype of Arabidopsis thaliana capacitates it for systemic infection of resistant ecotypes

Viral pathogens continue to emerge among humans, domesticated animals and cultivated crops. The existence of genetic variance for resistance in the host population is crucial to the spread of an emerging virus. Models predict that rapid spread decreases with the frequency and diversity of resistance alleles in the host population. However, empirical tests of this hypothesis are scarce. Arabiodpsis thaliana—tobacco etch potyvirus (TEV) provides an experimentally suitable pathosystem to explore the interplay between genetic variation in host''s susceptibility and virus diversity. Systemic infection of A. thaliana with TEV is controlled by three dominant loci, with different ecotypes varying in susceptibility depending on the genetic constitution at these three loci. Here, we show that the TEV adaptation to a susceptible ecotype allowed the virus to successfully infect, replicate and induce symptoms in ecotypes that were fully resistant to the ancestral virus. The value of these results is twofold. First, we showed that the existence of partially susceptible individuals allows for the emerging virus to bypass resistance alleles that the virus has never encountered. Second, the concept of resistance genes may only be valid for a well-defined viral genotype but not for polymorphic viral populations.     

The Impact of Recombination Hotspots on Genome Evolution of a Fungal Plant Pathogen

Recombination has an impact on genome evolution by maintaining chromosomal integrity, affecting the efficacy of selection, and increasing genetic variability in populations. Recombination rates are a key determinant of the coevolutionary dynamics between hosts and their pathogens. Historic recombination events created devastating new pathogens, but the impact of ongoing recombination in sexual pathogens is poorly understood. Many fungal pathogens of plants undergo regular sexual cycles, and sex is considered to be a major factor contributing to virulence. We generated a recombination map at kilobase-scale resolution for the haploid plant pathogenic fungus Zymoseptoria tritici. To account for intraspecific variation in recombination rates, we constructed genetic maps from two independent crosses. We localized a total of 10,287 crossover events in 441 progeny and found that recombination rates were highly heterogeneous within and among chromosomes. Recombination rates on large chromosomes were inversely correlated with chromosome length. Short accessory chromosomes often lacked evidence for crossovers between parental chromosomes. Recombination was concentrated in narrow hotspots that were preferentially located close to telomeres. Hotspots were only partially conserved between the two crosses, suggesting that hotspots are short-lived and may vary according to genomic background. Genes located in hotspot regions were enriched in genes encoding secreted proteins. Population resequencing showed that chromosomal regions with high recombination rates were strongly correlated with regions of low linkage disequilibrium. Hence, genes in pathogen recombination hotspots are likely to evolve faster in natural populations and may represent a greater threat to the host.     

Adaptation to abiotic conditions drives local adaptation in bacteria and viruses coevolving in heterogeneous environments

Parasite local adaptation, the greater performance of parasites on their local compared with foreign hosts, has important consequences for the maintenance of diversity and epidemiology. While the abiotic environment may significantly affect local adaptation, most studies to date have failed either to incorporate the effects of the abiotic environment, or to separate them from those of the biotic environment. Here, we tease apart biotic and abiotic components of local adaptation using the bacterium Pseudomonas fluorescens and its viral parasite bacteriophage Φ2. We coevolved replicate populations of bacteria and phages at three different temperatures, and determined their performance against coevolutionary partners from the same and different temperatures. Crucially, we measured performance at different assay temperatures, which allowed us to disentangle adaptation to biotic and abiotic habitat components. Our results show that bacteria and phages are more resistant and infectious, respectively, at the temperature at which they previously coevolved, confirming that local adaptation to abiotic conditions can play a crucial role in determining parasite infectivity and host resistance. Our work underlines the need to assess host–parasite interactions across multiple relevant abiotic environments, and suggests that microbial adaption to local temperatures can create ecological barriers to dispersal across temperature gradients.     

Enemies make you stronger: Coevolution between fruit fly host and bacterial pathogen increases postinfection survivorship in the host

Multiple laboratory studies have evolved hosts against a nonevolving pathogen to address questions about evolution of immune responses. However, an ecologically more relevant scenario is one where hosts and pathogens can coevolve. Such coevolution between the antagonists, depending on the mutual selection pressure and additive variance in the respective populations, can potentially lead to a different pattern of evolution in the hosts compared to a situation where the host evolves against a nonevolving pathogen. In the present study, we used Drosophila melanogaster as the host and Pseudomonas entomophila as the pathogen. We let the host populations either evolve against a nonevolving pathogen or coevolve with the same pathogen. We found that the coevolving hosts on average evolved higher survivorship against the coevolving pathogen and ancestral (nonevolving) pathogen relative to the hosts evolving against a nonevolving pathogen. The coevolving pathogens evolved greater ability to induce host mortality even in nonlocal (novel) hosts compared to infection by an ancestral (nonevolving) pathogen. Thus, our results clearly show that the evolved traits in the host and the pathogen under coevolution can be different from one‐sided adaptation. In addition, our results also show that the coevolving host–pathogen interactions can involve certain general mechanisms in the pathogen, leading to increased mortality induction in nonlocal or novel hosts.     

Partitioning the net effect of host diversity on an emerging amphibian pathogen

The ‘dilution effect’ (DE) hypothesis predicts that diverse host communities will show reduced disease. The underlying causes of pathogen dilution are complex, because they involve non-additive (driven by host interactions and differential habitat use) and additive (controlled by host species composition) mechanisms. Here, we used measures of complementarity and selection traditionally employed in the field of biodiversity–ecosystem function (BEF) to quantify the net effect of host diversity on disease dynamics of the amphibian-killing fungus Batrachochytrium dendrobatidis (Bd). Complementarity occurs when average infection load in diverse host assemblages departs from that of each component species in uniform populations. Selection measures the disproportionate impact of a particular species in diverse assemblages compared with its performance in uniform populations, and therefore has strong additive and non-additive properties. We experimentally infected tropical amphibian species of varying life histories, in single- and multi-host treatments, and measured individual Bd infection loads. Host diversity reduced Bd infection in amphibians through a mechanism analogous to complementarity (sensu BEF), potentially by reducing shared habitat use and transmission among hosts. Additionally, the selection component indicated that one particular terrestrial species showed reduced infection loads in diverse assemblages at the expense of neighbouring aquatic hosts becoming heavily infected. By partitioning components of diversity, our findings underscore the importance of additive and non-additive mechanisms underlying the DE.     

Culture–gene coevolution of individualism–collectivism and the serotonin transporter gene

Culture–gene coevolutionary theory posits that cultural values have evolved, are adaptive and influence the social and physical environments under which genetic selection operates. Here, we examined the association between cultural values of individualism–collectivism and allelic frequency of the serotonin transporter functional polymorphism (5-HTTLPR) as well as the role this culture–gene association may play in explaining global variability in prevalence of pathogens and affective disorders. We found evidence that collectivistic cultures were significantly more likely to comprise individuals carrying the short (S) allele of the 5-HTTLPR across 29 nations. Results further show that historical pathogen prevalence predicts cultural variability in individualism–collectivism owing to genetic selection of the S allele. Additionally, cultural values and frequency of S allele carriers negatively predict global prevalence of anxiety and mood disorder. Finally, mediation analyses further indicate that increased frequency of S allele carriers predicted decreased anxiety and mood disorder prevalence owing to increased collectivistic cultural values. Taken together, our findings suggest culture–gene coevolution between allelic frequency of 5-HTTLPR and cultural values of individualism–collectivism and support the notion that cultural values buffer genetically susceptible populations from increased prevalence of affective disorders. Implications of the current findings for understanding culture–gene coevolution of human brain and behaviour as well as how this coevolutionary process may contribute to global variation in pathogen prevalence and epidemiology of affective disorders, such as anxiety and depression, are discussed.     

Recent range expansion and agricultural landscape heterogeneity have only minimal effect on the spatial genetic structure of the plant pathogenic fungus Mycosphaerella fijiensis

Understanding how geographical and environmental features affect genetic variation at both the population and individual levels is crucial in biology, especially in the case of pathogens. However, distinguishing between these factors and the effects of historical range expansion on spatial genetic structure remains challenging. In the present study, we investigated the case of Mycosphaerella fijiensis—a plant pathogenic fungus that has recently colonized an agricultural landscape characterized by the presence of potential barriers to gene flow, including several commercial plantations in which disease control practises such as the use of fungicides are applied frequently, and low host density areas. We first genotyped 300 isolates sampled at a global scale on untreated plants in two dimensions over a 50 × 80-km area. Using two different clustering algorithms, no genetic structure was detected in the studied area, suggesting expansion of large populations and/or no influence of potential barriers. Second, we investigated the potential effect of disease control practises on M. fijiensis diversity by comparing populations sampled in commercial vs food-crop plantations. At this local scale, we detected significantly higher allelic richness inside commercial plantations compared with the surrounding food-crop plantation populations. Analysis of molecular variance indicated that 99% of the total genetic variance occurred within populations. We discuss the suggestion that high population size and/or high migration rate between populations might be responsible for the absence of any effect of disease control practises on genetic diversity and differentiation.     

Heritable variation in host tolerance and resistance inferred from a wild host–parasite system

Hosts have evolved two distinct defence strategies against parasites: resistance (which prevents infection or limit parasite growth) and tolerance (which alleviates the fitness consequences of infection). However, heritable variation in resistance and tolerance and the genetic correlation between these two traits have rarely been characterized in wild host populations. Here, we estimate these parameters for both traits in Leuciscus burdigalensis, a freshwater fish parasitized by Tracheliastes polycolpus. We used a genetic database to construct a full-sib pedigree in a wild L. burdigalensis population. We then used univariate animal models to estimate inclusive heritability (i.e. all forms of genetic and non-genetic inheritance) in resistance and tolerance. Finally, we assessed the genetic correlation between these two traits using a bivariate animal model. We found significant heritability for resistance (H = 17.6%; 95% CI: 7.2–32.2%) and tolerance (H = 18.8%; 95% CI: 4.4–36.1%), whereas we found no evidence for the existence of a genetic correlation between these traits. Furthermore, we confirm that resistance and tolerance are strongly affected by environmental effects. Our results demonstrate that (i) heritable variation exists for parasite resistance and tolerance in wild host populations, and (ii) these traits can evolve independently in populations.     

Intron sequences of arginine kinase in an intertidal snail suggest an ecotype-specific selective sweep and a gene duplication

Many species with restricted gene flow repeatedly respond similarly to local selection pressures. To fully understand the genetic mechanisms behind this process, the phylogeographic history of the species (inferred from neutral markers) as well as the loci under selection need to be known. Here we sequenced an intron in the arginine kinase gene (Ark), which shows strong clinal variation between two locally adapted ecotypes of the flat periwinkle, Littorina fabalis. The ‘small-sheltered'' ecotype was almost fixed for one haplotype, H1, in populations on both sides of the North Sea, unlike the ‘large-moderately exposed ecotype'', which segregated for ten different haplotypes. This contrasts with neutral markers, where the two ecotypes are equally variable. H1 could have been driven to high frequency in an ancestral population and then repeatedly spread to sheltered habitats due to local selection pressures with the colonization of both sides of the North Sea, after the last glacial maximum (∼18 000 years ago). An alternative explanation is that a positively selected mutation, in or linked to Ark, arose after the range expansion and secondarily spread through sheltered populations throughout the distribution range, causing this ecotype to evolve in a concerted fashion. Also, we were able to sequence up to four haplotypes consistently from some individuals, suggesting a gene duplication in Ark.     

Host behaviour and physiology underpin individual variation in avian influenza virus infection in migratory Bewick's swans

Individual variation in infection modulates both the dynamics of pathogens and their impact on host populations. It is therefore crucial to identify differential patterns of infection and understand the mechanisms responsible. Yet our understanding of infection heterogeneity in wildlife is limited, even for important zoonotic host-pathogen systems, owing to the intractability of host status prior to infection. Using novel applications of stable isotope ecology and eco-immunology, we distinguish antecedent behavioural and physiological traits associated with avian influenza virus (AIV) infection in free-living Bewick's swans (Cygnus columbianus bewickii). Swans infected with AIV exhibited higher serum δ13C (-25.3±0.4) than their non-infected counterparts (-26.3±0.2). Thus, individuals preferentially foraging in aquatic rather than terrestrial habitats experienced a higher risk of infection, suggesting that the abiotic requirements of AIV give rise to heterogeneity in pathogen exposure. Juveniles were more likely to be infected (30.8% compared with 11.3% for adults), shed approximately 15-fold higher quantity of virus and exhibited a lower specific immune response than adults. Together, these results demonstrate the potential for heterogeneity in infection to have a profound influence on the dynamics of pathogens, with concomitant impacts on host habitat selection and fitness.     

Host–Pathogen Coevolution: The Selective Advantage of Bacillus thuringiensis Virulence and Its Cry Toxin Genes

Reciprocal coevolution between host and pathogen is widely seen as a major driver of evolution and biological innovation. Yet, to date, the underlying genetic mechanisms and associated trait functions that are unique to rapid coevolutionary change are generally unknown. We here combined experimental evolution of the bacterial biocontrol agent Bacillus thuringiensis and its nematode host Caenorhabditis elegans with large-scale phenotyping, whole genome analysis, and functional genetics to demonstrate the selective benefit of pathogen virulence and the underlying toxin genes during the adaptation process. We show that: (i) high virulence was specifically favoured during pathogen–host coevolution rather than pathogen one-sided adaptation to a nonchanging host or to an environment without host; (ii) the pathogen genotype BT-679 with known nematocidal toxin genes and high virulence specifically swept to fixation in all of the independent replicate populations under coevolution but only some under one-sided adaptation; (iii) high virulence in the BT-679-dominated populations correlated with elevated copy numbers of the plasmid containing the nematocidal toxin genes; (iv) loss of virulence in a toxin-plasmid lacking BT-679 isolate was reconstituted by genetic reintroduction or external addition of the toxins. We conclude that sustained coevolution is distinct from unidirectional selection in shaping the pathogen''s genome and life history characteristics. To our knowledge, this study is the first to characterize the pathogen genes involved in coevolutionary adaptation in an animal host–pathogen interaction system.