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.     

Genetic variation in a host-parasite association: potential for coevolution and frequency-dependent selection

Abstract.— Models of host‐parasite coevolution assume the presence of genetic variation for host resistance and parasite infectivity, as well as genotype‐specific interactions. We used the freshwater crustacean Daphnia magna and its bacterial microparasite Pasteuria ramosa to study genetic variation for host susceptibility and parasite infectivity within each of two populations. We sought to answer the following questions: Do host clones differ in their susceptibility to parasite isolates? Do parasite isolates differ in their ability to infect different host clones? Are there host clone‐parasite isolate interactions? The analysis revealed considerable variation in both host resistance and parasite infectivity. There were significant host clone‐parasite isolate interactions, such that there was no single host clone that was superior to all other clones in the resistance to every parasite isolate. Likewise, there was no parasite isolate that was superior to all other isolates in infectivity to every host clone. This form of host clone‐parasite isolate interaction indicates the potential for coevolution based on frequency‐dependent selection. Infection success of original host clone‐parasite isolate combinations (i.e., those combinations that were isolated together) was significantly higher than infection success of novel host clone‐parasite isolate combinations (i.e., those combinations that were created in the laboratory). This finding is consistent with the idea that parasites track specific host genotypes under natural conditions. In addition, correspondence analysis revealed that some host clones, although distinguishable with neutral genetic markers, were susceptible to the same set of parasite isolates and thus probably shared resistance genes.     

Pathogen exposure reduces sexual dimorphism in a host’s upper thermal limits

The climate is warming at an unprecedented rate, pushing many species toward and beyond the upper temperatures at which they can survive. Global change is also leading to dramatic shifts in the distribution of pathogens. As a result, upper thermal limits and susceptibility to infection should be key determinants of whether populations continue to persist, or instead go extinct. Within a population, however, individuals vary in both their resistance to both heat stress and infection, and their contributions to vital growth rates. No more so is this true than for males and females. Each sex often varies in their response to pathogen exposure, thermal tolerances, and particularly their influence on population growth, owing to the higher parental investment that females typically make in their offspring. To date, the interplay between host sex, infection, and upper thermal limits has been neglected. Here, we explore the response of male and female Daphnia to bacterial infection and static heat stress. We find that female Daphnia, when uninfected, are much more resistant to static heat stress than males, but that infection negates any advantage that females are afforded. We discuss how the capacity of a population to cope with multiple stressors may be underestimated unless both sexes are considered simultaneously.     

Antagonistic coevolution with parasites maintains host genetic diversity: an experimental test

Genetic variation in natural populations is a prime prerequisite allowing populations to respond to selection, but is under constant threat from forces that tend to reduce it, such as genetic drift and many types of selection. Haldane emphasized the potential importance of parasites as a driving force of genetic diversity. His theory has been taken for granted ever since, but despite numerous studies showing correlations between genetic diversity and parasitism, Haldane''s hypothesis has rarely been tested experimentally for unambiguous support. We experimentally staged antagonistic coevolution between the host Tribolium castaneum and its natural microsporidian parasite, Nosema whitei, to test for the relative importance of two separate evolutionary forces (drift and parasite-induced selection) on the maintenance of genetic variation. Our results demonstrate that coevolution with parasites indeed counteracts drift as coevolving populations had significantly higher levels of heterozygosity and allelic diversity. Genetic drift remained a strong force, strongly reducing genetic variation and increasing genetic differentiation in small populations. To our surprise, differentiation between the evolving populations was smaller when they coevolved with parasites, suggesting parallel balancing selection. Hence, our results experimentally vindicate Haldane''s original hypothesis 60 years after its conception.     

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.     

Population-Genomic Analysis Identifies a Low Rate of Global Adaptive Fixation in the Proteins of the Cyclical Parthenogen Daphnia magna

Daphnia are well-established ecological and evolutionary models, and the interaction between D. magna and its microparasites is widely considered a paragon of the host-parasite coevolutionary process. Like other well-studied arthropods such as Drosophila melanogaster and Anopheles gambiae, D. magna is a small, widespread, and abundant species that is therefore expected to display a large long-term population size and high rates of adaptive protein evolution. However, unlike these other species, D. magna is cyclically asexual and lives in a highly structured environment (ponds and lakes) with moderate levels of dispersal, both of which are predicted to impact upon long-term effective population size and adaptive protein evolution. To investigate patterns of adaptive protein fixation, we produced the complete coding genomes of 36 D. magna clones sampled from across the European range (Western Palaearctic), along with draft sequences for the close relatives D. similis and D. lumholtzi, used as outgroups. We analyzed genome-wide patterns of adaptive fixation, with a particular focus on genes that have an a priori expectation of high rates, such as those likely to mediate immune responses, RNA interference against viruses and transposable elements, and those with a strongly male-biased expression pattern. We find that, as expected, D. magna displays high levels of diversity and that this is highly structured among populations. However, compared with Drosophila, we find that D. magna proteins appear to have a high proportion of weakly deleterious variants and do not show evidence of pervasive adaptive fixation across its entire range. This is true of the genome as a whole, and also of putative ‘arms race’ genes that often show elevated levels of adaptive substitution in other species. In addition to the likely impact of extensive, and previously documented, local adaptation, we speculate that these findings may reflect reduced efficacy of selection associated with cyclical asexual reproduction.     

Genetic variation in the cellular response of Daphnia magna (Crustacea: Cladocera) to its bacterial parasite

Linking measures of immune function with infection, and ultimately, host and parasite fitness is a major goal in the field of ecological immunology. In this study, we tested for the presence and timing of a cellular immune response in the crustacean Daphnia magna following exposure to its sterilizing endoparasite Pasteuria ramosa. We found that D. magna possesses two cell types circulating in the haemolymph: a spherical one, which we call a granulocyte and an irregular-shaped amoeboid cell first described by Metchnikoff over 125 years ago. Daphnia magna mounts a strong cellular response (of the amoeboid cells) just a few hours after parasite exposure. We further tested for, and found, considerable genetic variation for the magnitude of this cellular response. These data fostered a heuristic model of resistance in this naturally coevolving host–parasite interaction. Specifically, the strongest cellular responses were found in the most susceptible hosts, indicating resistance is not always borne from a response that destroys invading parasites, but rather stems from mechanisms that prevent their initial entry. Thus, D. magna may have a two-stage defence—a genetically determined barrier to parasite establishment and a cellular response once establishment has begun.     

Balancing selection in Pattern Recognition Receptor signalling pathways is associated with gene function and pleiotropy in a wild rodent

Pathogen‐mediated balancing selection is commonly considered to play an important role in the maintenance of genetic diversity, in particular in immune genes. However, the factors that may influence which immune genes are the targets of such selection are largely unknown. To address this, here we focus on Pattern Recognition Receptor (PRR) signalling pathways, which play a key role in innate immunity. We used whole‐genome resequencing data from a population of bank voles (Myodes glareolus) to test for associations between balancing selection, pleiotropy and gene function in a set of 123 PRR signalling pathway genes. To investigate the effect of gene function, we compared genes encoding (a) receptors for microbial ligands versus downstream signalling proteins, and (b) receptors recognizing components of microbial cell walls, flagella and capsids versus receptors recognizing features of microbial nucleic acids. Analyses based on the nucleotide diversity of full coding sequences showed that balancing selection primarily targeted receptor genes with a low degree of pleiotropy. Moreover, genes encoding receptors recognizing components of microbial cell walls etc. were more important targets of balancing selection than receptors recognizing nucleic acids. Tests for localized signatures of balancing selection in coding and noncoding sequences showed that such signatures were mostly located in introns, and more evenly distributed among different functional categories of PRR pathway genes. The finding that signatures of balancing selection in full coding sequences primarily occur in receptor genes, in particular those encoding receptors for components of microbial cell walls etc., is consistent with the idea that coevolution between hosts and pathogens is an important cause of balancing selection on immune genes.     

Too much of a good thing: resource provisioning alters infectious disease dynamics in wildlife

Provisioning of abundant food resources in human-altered landscapes can have profound effects on wildlife ecology, with important implications for pathogen transmission. While empirical studies have quantified the effects of provisioning on host behaviour and immunology, the net interactive effect of these components on host–pathogen dynamics is unknown. We use simple compartmental models to investigate how provisioning-induced changes to host demography, contact behaviour and immune defence influence pathogen invasion and persistence. We show that pathogen invasion success and equilibrium prevalence depend critically on how provisioning affects host immune defence and that moderate levels of provisioning can lead to drastically different outcomes of pathogen extinction or maximizing prevalence. These results highlight the need for further empirical studies to fully understand how provisioning affects pathogen transmission in urbanized environments.     

Resistance to RHD virus in wild Australian rabbits: Comparison of susceptible and resistant individuals using a genomewide approach

Deciphering the genes involved in disease resistance is essential if we are to understand host–pathogen coevolutionary processes. The rabbit haemorrhagic disease virus (RHDV) was imported into Australia in 1995 as a biocontrol agent to manage one of the most successful and devastating invasive species, the European rabbit (Oryctolagus cuniculus). During the first outbreaks of the disease, RHDV caused mortality rates of up to 97%. Recently, however, increased genetic resistance to RHDV has been reported. Here, we have aimed to identify genomic differences between rabbits that survived a natural infection with RHDV and those that died in the field using a genomewide next‐generation sequencing (NGS) approach. We detected 72 SNPs corresponding to 133 genes associated with survival of a RHD infection. Most of the identified genes have known functions in virus infections and replication, immune responses or apoptosis, or have previously been found to be regulated during RHD. Some of the genes identified in experimental studies, however, did not seem to play a role under natural selection regimes, highlighting the importance of field studies to complement the genomic background of wildlife diseases. Our study provides a set of candidate markers as a tool for the future scanning of wild rabbits for their resistance to RHDV. This is important both for wild rabbit populations in southern Europe where RHD is regarded as a serious problem decimating the prey of endangered predator species and for assessing the success of currently planned RHDV variant biocontrol releases in Australia.     

Red Queen dancing in the lek: effects of mating skew on host–parasite interactions

The RQH (Red Queen hypothesis), which argues that hosts need to be continuously finding new ways to avoid parasites that are able to infect common host genotypes, has been at the center of discussions on the maintenance of sex. This is because diversity is favored under the host–parasite coevolution based on negative frequency‐dependent selection, and sexual reproduction is a mechanism that generates genetic diversity in the host population. Together with parasite infections, sexual organisms are usually under sexual selection, which leads to mating skew or mating success biased toward males with a particular phenotype. Thus, strong mating skew would affect genetic variance in a population and should affect the benefit of the RQH. However, most models have investigated the RQH under a random mating system and not under mating skew. In this study, I show that sexual selection and the resultant mating skew may increase parasite load in the hosts. An IBM (individual‐based model), which included host–parasite interactions and sexual selection among hosts, demonstrates that mating skew influenced parasite infection in the hosts under various conditions. Moreover, the IBM showed that the mating skew evolves easily in cases of male–male competition and female mate choice, even though it imposes an increased risk of parasite infection on the hosts. These findings indicated that whether the RQH favored sexual reproduction depended on the condition of mating skew. That is, consideration of the host mating system would provide further understanding of conditions in which the RQH favors sexual reproduction in real organisms.     

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa,Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a “core” genome underlying basic cellular processes and a “flexible” genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.     

Thermoregulatory behaviour affects prevalence of chytrid fungal infection in a wild population of Panamanian golden frogs

Predicting how climate change will affect disease dynamics requires an understanding of how the environment affects host–pathogen interactions. For amphibians, global declines and extinctions have been linked to a pathogenic chytrid fungus, Batrachochytrium dendrobatidis. Using a combination of body temperature measurements and disease assays conducted before and after the arrival of B. dendrobatidis, this study tested the hypothesis that body temperature affects the prevalence of infection in a wild population of Panamanian golden frogs (Atelopus zeteki). The timing of first detection of the fungus was consistent with that of a wave of epidemic infections spreading south and eastward through Central America. During the epidemic, many golden frogs modified their thermoregulatory behaviour, raising body temperatures above their normal set point. Odds of infection decreased with increasing body temperature, demonstrating that even slight environmental or behavioural changes have the potential to affect an individual''s vulnerability to infection. The thermal dependency of the relationship between B. dendrobatidis and its amphibian hosts demonstrates how the progression of an epidemic can be influenced by complex interactions between host and pathogen phenotypes and the environments in which they are found.