The genetic basis of resistance and matching-allele interactions of a host-parasite system: The Daphnia magna-Pasteuria ramosa model

Long-term balancing selection for pathogen resistance maintains trans-species polymorphisms in a planktonic crustacean

Balancing selection is an evolutionary process that maintains genetic polymorphisms at selected loci and strongly reduces the likelihood of allele fixation. When allelic polymorphisms that predate speciation events are maintained independently in the resulting lineages, a pattern of trans-species polymorphisms may occur. Trans-species polymorphisms have been identified for loci related to mating systems and the MHC, but they are generally rare. Trans-species polymorphisms in disease loci are believed to be a consequence of long-term host-parasite coevolution by balancing selection, the so-called Red Queen dynamics. Here we scan the genomes of three crustaceans with a divergence of over 15 million years and identify 11 genes containing identical-by-descent trans-species polymorphisms with the same polymorphisms in all three species. Four of these genes display molecular footprints of balancing selection and have a function related to immunity. Three of them are located in or close to loci involved in resistance to a virulent bacterial pathogen, Pasteuria, with which the Daphnia host is known to coevolve. This provides rare evidence of trans-species polymorphisms for loci known to be functionally relevant in interactions with a widespread and highly specific parasite. These findings support the theory that specific antagonistic coevolution is able to maintain genetic diversity over millions of years.

The cuticle proteome of a planktonic crustacean

The cuticles of arthropods provide an interface between the organism and its environment. Thus, the cuticle's structure influences how the organism responds to and interacts with its surroundings. Here, we used label-free quantification proteomics to provide a proteome of the moulted cuticle of the aquatic crustacean Daphnia magna, which has long been a prominent subject of studies on ecology, evolution, and developmental biology. We detected a total of 278 high-confidence proteins. Using protein sequence domain and functional enrichment analyses, we identified chitin-binding structural proteins and chitin-modifying enzymes as the most abundant protein groups in the cuticle proteome. Structural cuticular protein families showed a similar distribution to those found in other arthropods and indicated proteins responsible for the soft and flexible structure of the Daphnia cuticle. Finally, cuticle protein genes were also clustered as tandem gene arrays in the D. magna genome. The cuticle proteome presented here will be a valuable resource to the Daphnia research community, informing genome annotations and investigations on diverse topics such as the genetic basis of interactions with predators and parasites.

A common multi-host parasite shows genetic structuring at the host species and population levels

Although individual parasite species commonly infect many populations across physical space as well as multiple host species, the extent to which parasites traverse physical and phylogenetic distances is unclear. Population genetic analyses of parasite populations can reveal how parasites move across space or between host species, including helping assess whether a parasite is more likely to infect a different host species in the same location or the same host species in a different location. Identifying these transmission barriers could be exploited for effective disease control. Here, we analysed population genetic structuring of the parasite Pasteuria ramosa in daphniid host species from different lakes. Outbreaks occurred most often in the common host species Daphnia dentifera and Daphnia retrocurva. The genetic distance between parasite samples tended to be smaller when samples were collected from the same lake, the same host species and closer in time. Within lakes, the parasite showed structure by host species and sampling date; within a host species, the parasite showed structure by lake and sampling date. However, despite this structuring, we found the same parasite genotype infecting closely related host species, and we sometimes found the same genotype in nearby lakes. Thus, P. ramosa experiences challenges infecting different host species and moving between populations, but doing so is possible. In addition, the structuring by sampling date indicates potential adaptation to or coevolution with host populations and supports prior findings that parasite population structure is dynamic during outbreaks.

Ectopical expression of bacterial collagen-like protein supports its role as adhesin in host-parasite coevolution

For a profound understanding of antagonistic coevolution, it is necessary to identify the coevolving genes. The bacterium Pasteuria and its host, the microcrustacean Daphnia, are a well-characterized paradigm for co-evolution, but the underlying genes remain largely unknown. A genome-wide association study suggested a Pasteuria collagen-like protein 7 (Pcl7) as a candidate mediating parasite attachment and driving its coevolution with the host. Since Pasteuria ramosa cannot currently be genetically manipulated, we used Bacillus thuringiensis to express a fusion protein of a Pcl7 carboxy-terminus from P. ramosa and the amino-terminal domain of a B. thuringiensis collagen-like protein (CLP). Mutant B. thuringiensis (Pcl7-Bt) spores but not wild-type B. thuringiensis (WT-Bt) spores attached to the same site of susceptible hosts as P. ramosa. Furthermore, Pcl7-Bt spores attached readily to susceptible host genotypes, but only slightly to resistant host genotypes. These findings indicated that the fusion protein was properly expressed and folded and demonstrated that indeed the C-terminus of Pcl7 mediates attachment in a host genotype-specific manner. These results provide strong evidence for the involvement of a CLP in the coevolution of Daphnia and P. ramosa and open new avenues for genetic epidemiological studies of host-parasite interactions.

Dominance and multi-locus interaction

Dominance is usually considered a constant value that describes the relative difference in fitness or phenotype between heterozygotes and the average of homozygotes at a focal polymorphic locus. However, the observed dominance can vary with the genetic background of the focal locus. Here, alleles at other loci modify the observed phenotype through position effects or dominance modifiers that are sometimes associated with pathogen resistance, lineage, sex, or mating type. Theoretical models have illustrated how variable dominance appears in the context of multi-locus interaction (epistasis). Here, we review empirical evidence for variable dominance and how the observed patterns may be captured by proposed epistatic models. We highlight how integrating epistasis and dominance is crucial for comprehensively understanding adaptation and speciation.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Rapid evolution of a bacterial parasite during outbreaks in two Daphnia populations

Myriad ecological and evolutionary factors can influence whether a particular parasite successfully transmits to a new host during a disease outbreak, with consequences for the structure and diversity of parasite populations. However, even though the diversity and evolution of parasite populations are of clear fundamental and applied importance, we have surprisingly few studies that track how genetic structure of parasites changes during naturally occurring outbreaks in non-human populations. Here, we used population genetic approaches to reveal how genotypes of a bacterial parasite, Pasteuria ramosa, change over time, focusing on how infecting P. ramosa genotypes change during the course of epidemics in Daphnia populations in two lakes. We found evidence for genetic change - and, therefore, evolution - of the parasite during outbreaks. In one lake, P. ramosa genotypes were structured by sampling date; in both lakes, genetic distance between groups of P. ramosa isolates increased with time between sampling. Diversity in parasite populations remained constant over epidemics, although one epidemic (which was large) had low genetic diversity while the other epidemic (which was small) had high genetic diversity. Our findings demonstrate that patterns of parasite evolution differ between outbreaks; future studies exploring the feedbacks among epidemic size, host diversity, and parasite genetic diversity would improve our understanding of parasite dynamics and evolution.

Genetic variation and microbiota in bumble bees cross-infected by different strains of C. bombi

The bumblebee Bombus terrestris is commonly infected by a trypanosomatid gut parasite Crithidia bombi. This system shows a striking degree of genetic specificity where host genotypes are susceptible to different genotypes of parasite. To a degree, variation in host gene expression underlies these differences, however, the effects of standing genetic variation has not yet been explored. Here we report on an extensive experiment where workers of twenty colonies of B. terrestris were each infected by one of twenty strains of C. bombi. To elucidate the host's genetic bases of susceptibility to infection (measured as infection intensity), we used a low-coverage (~2 x) genome-wide association study (GWAS), based on angsd, and a standard high-coverage (~15x) GWAS (with a reduced set from a 8 x 8 interaction matrix, selected from the full set of twenty). The results from the low-coverage approach remained ambiguous. The high-coverage approach suggested potentially relevant genetic variation in cell surface and adhesion processes. In particular, mucin, a surface mucoglycoprotein, potentially affecting parasite binding to the host gut epithelia, emerged as a candidate. Sequencing the gut microbial community of the same bees showed that the abundance of bacterial taxa, such as Gilliamella, Snodgrassella, or Lactobacillus, differed between 'susceptible' and 'resistant' microbiota, in line with earlier studies. Our study suggests that the constitutive microbiota and binding processes at the cell surface are candidates to affect infection intensity after the first response (captured by gene expression) has run its course. We also note that a low-coverage approach may not be powerful enough to analyse such complex traits. Furthermore, testing large interactions matrices (as with the full 20 x 20 combinations) for the effect of interaction terms on infection intensity seems to blur the specific host x parasite interaction effects, likely because the outcome of an infection is a highly non-linear process dominated by variation in individually different pathways of host defence (immune) responses.

Genetic slippage after sex maintains diversity for parasite resistance in a natural host population

Although parasite-mediated selection is a major driver of host evolution, its influence on genetic variation for parasite resistance is not yet well understood. We monitored resistance in a large population of the planktonic crustacean Daphnia magna over 8 years, as it underwent yearly epidemics of the bacterial pathogen Pasteuria ramosa. We observed cyclic dynamics of resistance: Resistance increased throughout the epidemics, but susceptibility was restored each spring when hosts hatched from sexual resting stages. Host resting stages collected across the year showed that largely resistant host populations can produce susceptible sexual offspring. A genetic model of resistance developed for this host-parasite system, based on multiple loci and strong epistasis, is in partial agreement with our findings. Our results reveal that, despite strong selection for resistance in a natural host population, genetic slippage after sexual reproduction can be a strong factor for the maintenance of genetic diversity of host resistance.

MHC class II genotype-by-pathogen genotype interaction for infection prevalence in a natural rodent-Borrelia system

MHC genes are extraordinarily polymorphic in most taxa. Host-pathogen coevolution driven by negative frequency-dependent selection (NFDS) is one of the main hypotheses for the maintenance of such immunogenetic variation. Here, we test a critical but rarely tested assumption of this hypothesis-that MHC alleles affect resistance/susceptibility to a pathogen in a strain-specific way, that is, there is a host genotype-by-pathogen genotype interaction. In a field study of bank voles naturally infected with the tick-transmitted bacterium Borrelia afzelii, we tested for MHC class II (DQB) genotype-by-B. afzelii strain interactions for infection prevalence between 10 DQB alleles and seven strains. One allele (DQB*37) showed an interaction, such that voles carrying DQB*37 had higher prevalence of two strains and lower prevalence of one strain than individuals without the allele. These findings were corroborated by analyses of strain composition of infections, which revealed an effect of DQB*37 in the form of lower β diversity among infections in voles carrying the allele. Taken together, these results provide rare support at the molecular genetic level for a key assumption of models of antagonistic coevolution through NFDS.

Daphnia as a versatile model system in ecology and evolution

Water fleas of the genus Daphnia have been a model system for hundreds of years and is among the best studied ecological model organisms to date. Daphnia are planktonic crustaceans with a cyclic parthenogenetic life-cycle. They have a nearly worldwide distribution, inhabiting standing fresh- and brackish water bodies, from small temporary pools to large lakes. Their predominantly asexual reproduction allows for the study of phenotypes excluding genetic variation, enabling us to separate genetic from non-genetic effects. Daphnia are often used in studies related to ecotoxicology, predator-induced defence, host–parasite interactions, phenotypic plasticity and, increasingly, in evolutionary genomics. The most commonly studied species are Daphnia magna and D. pulex, for which a rapidly increasing number of genetic and genomic tools are available. Here, I review current research topics, where the Daphnia model system plays a critical role.

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.

Symbiosis in Digital Evolution: Past, Present, and Future

Symbiosis, the living together of unlike organisms as symbionts, is ubiquitous in the natural world. Symbioses occur within and across all scales of life, from microbial to macro-faunal systems. Further, the interactions between symbionts are multimodal in both strength and type, can span from parasitic to mutualistic within one partnership, and persist over generations. Studying the ecological and evolutionary dynamics of symbiosis in natural or laboratory systems poses a wide range of challenges, including the long time scales at which symbioses evolve de novo, the limited capacity to experimentally control symbiotic interactions, the weak resolution at which we can quantify interactions, and the idiosyncrasies of current model systems. These issues are especially challenging when seeking to understand the ecological effects and evolutionary pressures on and of a symbiosis, such as how a symbiosis may shift between parasitic and mutualistic modes and how that shift impacts the dynamics of the partner population. In digital evolution, populations of computational organisms compete, mutate, and evolve in a virtual environment. Digital evolution features perfect data tracking and allows for experimental manipulations that are impractical or impossible in natural systems. Furthermore, modern computational power allows experimenters to observe thousands of generations of evolution in minutes (as opposed to several months or years), which greatly expands the range of possible studies. As such, digital evolution is poised to become a keystone technique in our methodological repertoire for studying the ecological and evolutionary dynamics of symbioses. Here, we review how digital evolution has been used to study symbiosis, and we propose a series of open questions that digital evolution is well-positioned to answer.

Infection phenotypes of a coevolving parasite are highly diverse, structured, and specific

Understanding how diversity is maintained in natural populations is a major goal of evolutionary biology. In coevolving hosts and parasites, negative frequency-dependent selection is one mechanism predicted to maintain genetic variation. While much is known about host diversity, parasite diversity remains understudied in coevolutionary research. Here, we survey natural diversity in a bacterial parasite by characterizing infection phenotypes for over 50 isolates in relation to 12 genotypes of their host, Daphnia magna. We find striking phenotypic variation among parasite isolates, and we discover the parasite can infect its host through at least five different attachment sites. Variation in attachment success at each site is explained to varying degrees by host and parasite genotypes. A spatial correlation analysis showed that infectivity of different isolates does not correlate with geographic distance, meaning isolates from widespread populations are equally able to infect the host. Overall, our results reveal that infection phenotypes of this parasite are highly diverse. Our results are consistent with the prediction that under Red Queen coevolutionary dynamics both the host and the parasite should show high genetic diversity for traits of functional importance in their interactions.

Quantitative trait locus analysis of parasitoid counteradaptation to symbiont-conferred resistance

Insect hosts and parasitoids are engaged in an intense struggle of antagonistic coevolution. Infection with heritable bacterial endosymbionts can substantially increase the resistance of aphids to parasitoid wasps, which exerts selection on parasitoids to overcome this symbiont-conferred protection (counteradaptation). Experimental evolution in the laboratory has produced counteradapted populations of the parasitoid wasp Lysiphlebus fabarum. These populations can parasitize black bean aphids (Aphis fabae) protected by the bacterial endosymbiont Hamiltonella defensa, which confers high resistance against L. fabarum. We used two experimentally evolved parasitoid populations to study the genetic architecture of the counteradaptation to symbiont-conferred resistance by QTL analysis. With simple crossing experiments, we showed that the counteradaptation is a recessive trait depending on the maternal genotype. Based on these results, we designed a customized crossing scheme to genotype a mapping population phenotyped for the ability to parasitize Hamiltonella-protected aphids. Using 1835 SNP markers obtained by ddRAD sequencing, we constructed a high-density linkage map consisting of six linkage groups (LGs) with an overall length of 828.3 cM and an average marker spacing of 0.45 cM. We identified a single QTL associated with the counteradaptation to Hamiltonella in L. fabarum on linkage group 2. Out of 120 genes located in this QTL, several genes encoding putative venoms may represent candidates for counteradaptation, as parasitoid wasps inject venoms into their hosts during oviposition.

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.

A Two-Locus System with Strong Epistasis Underlies Rapid Parasite-Mediated Evolution of Host Resistance

Parasites are a major evolutionary force, driving adaptive responses in host populations. Although the link between phenotypic response to parasite-mediated natural selection and the underlying genetic architecture often remains obscure, this link is crucial for understanding the evolution of resistance and predicting associated allele frequency changes in the population. To close this gap, we monitored the response to selection during epidemics of a virulent bacterial pathogen, Pasteuria ramosa, in a natural host population of Daphnia magna. Across two epidemics, we observed a strong increase in the proportion of resistant phenotypes as the epidemics progressed. Field and laboratory experiments confirmed that this increase in resistance was caused by selection from the local parasite. Using a genome-wide association study, we built a genetic model in which two genomic regions with dominance and epistasis control resistance polymorphism in the host. We verified this model by selfing host genotypes with different resistance phenotypes and scoring their F1 for segregation of resistance and associated genetic markers. Such epistatic effects with strong fitness consequences in host–parasite coevolution are believed to be crucial in the Red Queen model for the evolution of genetic recombination.

Novel genomic approaches to study antagonistic coevolution between hosts and parasites

Host-parasite coevolution is ubiquitous, shaping genetic and phenotypic diversity and the evolutionary trajectory of interacting species. With the advances of high throughput sequencing technologies applicable to model and non-model organisms alike, it is now feasible to study in greater detail (a) the genetic underpinnings of coevolution, (b) the speed and type of dynamics at coevolving loci, and (c) the genomic consequences of coevolution. This review focuses on three recently developed approaches that leverage information from host and parasite full genome data simultaneously to pinpoint coevolving loci and draw inference on the coevolutionary history. First, co-genome-wide association study (co-GWAS) methods allow pinpointing the loci underlying host-parasite interactions. These methods focus on detecting associations between genetic variants and the outcome of experimental infection tests or on correlations between genomes of naturally infected hosts and their infecting parasites. Second, extensions to population genomics methods can detect genes under coevolution and infer the coevolutionary history, such as fitness costs. Third, correlations between host and parasite population size in time are indicative of coevolution, and polymorphism levels across independent spatially distributed populations of hosts and parasites can reveal coevolutionary loci and infer coevolutionary history. We describe the principles of these three approaches and discuss their advantages and limitations based on coevolutionary theory. We present recommendations for their application to various host (prokaryotes, fungi, plants, and animals) and parasite (viruses, bacteria, fungi, and macroparasites) species. We conclude by pointing out methodological and theoretical gaps to be filled to extract maximum information from full genome data and thereby to shed light on the molecular underpinnings of coevolution.

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.

Ecology directs host-parasite coevolutionary trajectories across Daphnia-microparasite populations

Host-parasite interactions often fuel coevolutionary change. However, parasitism is one of a myriad of possible ecological interactions in nature. Biotic (for example, predation) and abiotic (for example, temperature) variation can amplify or dilute parasitism as a selective force on hosts and parasites, driving population variation in (co)evolutionary trajectories. We dissected the relationships between wider ecology and coevolutionary trajectory using 16 ecologically complex Daphnia magna-Pasteuria ramosa ponds seeded with an identical starting host (Daphnia) and parasite (Pasteuria) population. We show, using a time-shift experiment and outdoor population data, how multivariate biotic and abiotic ecological differences between ponds caused coevolutionary divergence. Wider ecology drove variation in host evolution of resistance, but not parasite infectivity; parasites subsequently coevolved in response to the changing complement of host genotypes, such that parasites adapted to historically resistant host genotypes. Parasitism was a stronger interaction for the parasite than for its host, probably because the host is the principal environment and selective force, whereas for hosts, parasite-mediated selection is one of many sources of selection. Our findings reveal the mechanisms through which wider ecology creates coevolutionary hotspots and coldspots in biologically realistic arenas of host-parasite interaction, and sheds light on how the ecological theatre can affect the (co)evolutionary play.

Host-parasite co-evolution and its genomic signature

Studies in diverse biological systems have indicated that host-parasite co-evolution is responsible for the extraordinary genetic diversity seen in some genomic regions, such as major histocompatibility (MHC) genes in jawed vertebrates and resistance genes in plants. This diversity is believed to evolve under balancing selection on hosts by parasites. However, the mechanisms that link the genomic signatures in these regions to the underlying co-evolutionary process are only slowly emerging. We still lack a clear picture of the co-evolutionary concepts and of the genetic basis of the co-evolving phenotypic traits in the interacting antagonists. Emerging genomic tools that provide new options for identifying underlying genes will contribute to a fuller understanding of the co-evolutionary process.

Clusters of polymorphic transmembrane genes control resistance to schistosomes in snail vectors

Schistosomiasis is a debilitating parasitic disease infecting hundreds of millions of people. Schistosomes use aquatic snails as intermediate hosts. A promising avenue for disease control involves leveraging innate host mechanisms to reduce snail vectorial capacity. In a genome-wide association study of Biomphalaria glabrata snails, we identify genomic region PTC2 which exhibits the largest known correlation with susceptibility to parasite infection (>15 fold effect). Using new genome assemblies with substantially higher contiguity than the Biomphalaria reference genome, we show that PTC2 haplotypes are exceptionally divergent in structure and sequence. This variation includes multi-kilobase indels containing entire genes, and orthologs for which most amino acid residues are polymorphic. RNA-Seq annotation reveals that most of these genes encode single-pass transmembrane proteins, as seen in another resistance region in the same species. Such groups of hyperdiverse snail proteins may mediate host-parasite interaction at the cell surface, offering promising targets for blocking the transmission of schistosomiasis.

An alternative route of bacterial infection associated with a novel resistance locus in the Daphnia-Pasteuria host-parasite system

To understand the mechanisms of antagonistic coevolution, it is crucial to identify the genetics of parasite resistance. In the Daphnia magna–Pasteuria ramosa host–parasite system, the most important step of the infection process is the one in which P. ramosa spores attach to the host’s foregut. A matching-allele model (MAM) describes the host–parasite genetic interactions underlying attachment success. Here we describe a new P. ramosa genotype, P15, which, unlike previously studied genotypes, attaches to the host’s hindgut, not to its foregut. Host resistance to P15 attachment shows great diversity across natural populations. In contrast to P. ramosa genotypes that use foregut attachment, P15 shows some quantitative variation in attachment success and does not always lead to successful infections, suggesting that hindgut attachment represents a less-efficient infection mechanism than foregut attachment. Using a Quantitative Trait Locus (QTL) approach, we detect two significant QTLs in the host genome: one that co-localizes with the previously described D. magna PR locus of resistance to foregut attachment, and a second, major QTL located in an unlinked genomic region. We find no evidence of epistasis. Fine mapping reveals a genomic region, the D locus, of ~13 kb. The discovery of a second P. ramosa attachment site and of a novel host-resistance locus increases the complexity of this system, with implications for both for the coevolutionary dynamics (e.g., Red Queen and the role of recombination), and for the evolution and epidemiology of the infection process.

Exome-wide association study reveals largely distinct gene sets underlying specific resistance to dengue virus types 1 and 3 in Aedes aegypti

Although specific interactions between host and pathogen genotypes have been well documented in invertebrates, the identification of host genes involved in discriminating pathogen genotypes remains a challenge. In the mosquito Aedes aegypti, the main dengue virus (DENV) vector worldwide, statistical associations between host genetic markers and DENV types or strains were previously detected, but the host genes underlying this genetic specificity have not been identified. In particular, it is unknown whether DENV type- or strain-specific resistance relies on allelic variants of the same genes or on distinct gene sets. Here, we investigated the genetic architecture of DENV resistance in a population of Ae. aegypti from Bakoumba, Gabon, which displays a stronger resistance phenotype to DENV type 1 (DENV-1) than to DENV type 3 (DENV-3) infection. Following experimental exposure to either DENV-1 or DENV-3, we sequenced the exomes of large phenotypic pools of mosquitoes that are either resistant or susceptible to each DENV type. Using variation in single-nucleotide polymorphism (SNP) frequencies among the pools, we computed empirical p values based on average gene scores adjusted for the differences in SNP counts, to identify genes associated with infection in a DENV type-specific manner. Among the top 5% most significant genes, 263 genes were significantly associated with resistance to both DENV-1 and DENV-3, 287 genes were only associated with DENV-1 resistance and 290 were only associated with DENV-3 resistance. The shared significant genes were enriched in genes with ATP binding activity and sulfur compound transmembrane transporter activity, whereas the genes uniquely associated with DENV-3 resistance were enriched in genes with zinc ion binding activity. Together, these results indicate that specific resistance to different DENV types relies on largely non-overlapping sets of genes in this Ae. aegypti population and pave the way for further mechanistic studies.

Compatibility between hosts and pathogens is often genetically specific in invertebrates but host genes underlying this genetic specificity have not been elucidated. We investigated the genetic architecture of dengue virus type-specific resistance in the mosquito vector Aedes aegypti. We used a natural population of Ae. aegypti from Bakoumba, Gabon, which is differentially resistant to dengue virus type 1 and dengue virus type 3. We surveyed genetic variation in protein-coding regions of the mosquito genome and compared the frequency of genetic polymorphisms between groups of mosquitoes that are either resistant or susceptible to each dengue virus type. We found that the Ae. aegypti genes associated with resistance to dengue virus type 1 or dengue virus type 3 were largely non-overlapping. This finding indicates that different sets of host genes, rather than different variants of the same genes, confer pathogen-specific resistance in this population. This study is an important step towards identification of mechanisms underlying the genetic specificity of invertebrate host-pathogen interactions.

Neutral genomic signatures of host-parasite coevolution

Background

Coevolution is a selective process of reciprocal adaptation in hosts and parasites or in mutualistic symbionts. Classic population genetics theory predicts the signatures of selection at the interacting loci of both species, but not the neutral genome-wide polymorphism patterns. To bridge this gap, we build an eco-evolutionary model, where neutral genomic changes over time are driven by a single selected locus in hosts and parasites via a simple biallelic gene-for-gene or matching-allele interaction. This coevolutionary process may lead to cyclic changes in the sizes of the interacting populations.

Results

We investigate if and when these changes can be observed in the site frequency spectrum of neutral polymorphisms from host and parasite full genome data. We show that changes of the host population size are too smooth to be observable in its polymorphism pattern over the course of time. Conversely, the parasite population may undergo a series of strong bottlenecks occurring on a slower relative time scale, which may lead to observable changes in a time series sample. We also extend our results to cases with 1) several parasites per host accelerating relative time, and 2) multiple parasite generations per host generation slowing down rescaled time.

Conclusions

Our results show that time series sampling of host and parasite populations with full genome data are crucial to understand if and how coevolution occurs. This model provides therefore a framework to interpret and draw inference from genome-wide polymorphism data of interacting species.

Independent effects on cellular and humoral immune responses underlie genotype-by-genotype interactions between Drosophila and parasitoids

It is common to find abundant genetic variation in host resistance and parasite infectivity within populations, with the outcome of infection frequently depending on genotype-specific interactions. Underlying these effects are complex immune defenses that are under the control of both host and parasite genes. We have found extensive variation in Drosophila melanogaster's immune response against the parasitoid wasp Leptopilina boulardi. Some aspects of the immune response, such as phenoloxidase activity, are predominantly affected by the host genotype. Some, such as upregulation of the complement-like protein Tep1, are controlled by the parasite genotype. Others, like the differentiation of immune cells called lamellocytes, depend on the specific combination of host and parasite genotypes. These observations illustrate how the outcome of infection depends on independent genetic effects on different aspects of host immunity. As parasite-killing results from the concerted action of different components of the immune response, these observations provide a physiological mechanism to generate phenomena like epistasis and genotype-interactions that underlie models of coevolution.

Selection at behavioural, developmental and metabolic genes is associated with the northward expansion of a successful tropical colonizer

What makes a species able to colonize novel environments? This question is key to understand the dynamics of adaptive radiations and ecological niche shifts, but the mechanisms that underlie expansion into novel habitats remain poorly understood at a genomic scale. Lizards from the genus Anolis are typically tropical, and the green anole (Anolis carolinensis) constitutes an exception since it expanded into temperate North America from subtropical Florida. Thus, we used the green anole as a model to investigate signatures of selection associated with colonization of a new environment, namely temperate North America. To this end, we analysed 29 whole-genome sequences, covering the entire native range of the species. We used a combination of recent methods to quantify both positive and balancing selection in northern populations, including F_ST outlier methods, machine learning and ancestral recombination graphs. We naively scanned for genes of interest and assessed the overlap between multiple tests. Strikingly, we identified many genes involved in behaviour, suggesting that the recent successful colonization of northern environments may have been linked to behavioural shifts as well as physiological adaptation. Using a candidate genes strategy, we determined that genes involved in response to cold or behaviour displayed more frequently signals of selection, while controlling for local recombination rate, gene clustering and gene length. In addition, we found signatures of balancing selection at immune genes in all investigated genetic groups, but also at genes involved in neuronal and anatomical development.

Dissecting the genetic architecture of a stepwise infection process

How a host fights infection depends on an ordered sequence of steps, beginning with attempts to prevent a pathogen from establishing an infection, through to steps that mitigate a pathogen's control of host resources or minimize the damage caused during infection. Yet empirically characterizing the genetic basis of these steps remains challenging. Although each step is likely to have a unique genetic and environmental signature, and may therefore respond to selection in different ways, events that occur earlier in the infection process can mask or overwhelm the contributions of subsequent steps. In this study, we dissect the genetic architecture of a stepwise infection process using a quantitative trait locus (QTL) mapping approach. We control for variation at the first line of defence against a bacterial pathogen and expose downstream genetic variability related to the host's ability to mitigate the damage pathogens cause. In our model, the water-flea Daphnia magna, we found a single major effect QTL, explaining 64% of the variance, that is linked to the host's ability to completely block pathogen entry by preventing their attachment to the host oesophagus; this is consistent with the detection of this locus in previous studies. In susceptible hosts allowing attachment, however, a further 23 QTLs, explaining between 5% and 16% of the variance, were mapped to traits related to the expression of disease. The general lack of pleiotropy and epistasis for traits related to the different stages of the infection process, together with the wide distribution of QTLs across the genome, highlights the modular nature of a host's defence portfolio, and the potential for each different step to evolve independently. We discuss how isolating the genetic basis of individual steps can help to resolve discussion over the genetic architecture of host resistance.

Contrasting patterns of divergence at the regulatory and sequence level in European Daphnia galeata natural populations

Understanding the genetic basis of local adaptation has long been a focus of evolutionary biology. Recently, there has been increased interest in deciphering the evolutionary role of Daphnia's plasticity and the molecular mechanisms of local adaptation. Using transcriptome data, we assessed the differences in gene expression profiles and sequences in four European Daphnia galeata populations. In total, ~33% of 32,903 transcripts were differentially expressed between populations. Among 10,280 differentially expressed transcripts, 5,209 transcripts deviated from neutral expectations and their population-specific expression pattern is likely the result of local adaptation processes. Furthermore, a SNP analysis allowed inferring population structure and distribution of genetic variation. The population divergence at the sequence level was comparatively higher than the gene expression level by several orders of magnitude consistent with strong founder effects and lack of gene flow between populations. Using sequence homology, the candidate transcripts were annotated using a comparative genomics approach. Additionally, we also performed a weighted gene co-expression analysis to identify population-specific regulatory patterns of transcripts in D. galeata. Thus, we identified candidate transcriptomic regions for local adaptation in this key species of aquatic ecosystems in the absence of any laboratory-induced stressor.

A Single Interacting Species Leads to Widespread Parallel Evolution of the Stickleback Genome

Biotic interactions are potent, widespread causes of natural selection and divergent phenotypic evolution and can lead to genetic differentiation with gene flow among wild populations ("isolation by ecology") [1-4]. Biotic selection has been predicted to act on more genes than abiotic selection thereby driving greater adaptation [5]. However, difficulties in isolating the genome-wide effect of single biotic agents of selection have limited our ability to identify and quantify the number and type of genetic regions responding to biotic selection [6-9]. We identified geographically interspersed lakes in which threespine stickleback fish (Gasterosteus aculeatus) have repeatedly adapted to the presence or absence of a single member of the ecological community, prickly sculpin (Cottus asper), a fish that is both a competitor and a predator of the stickleback [10]. Whole-genome sequencing revealed that sculpin presence or absence accounted for the majority of genetic divergence among stickleback populations, more so than geography. The major axis of genomic variation within and between the two lake types was correlated with multiple traits, indicating parallel natural selection across a gradient of biotic environments. A large proportion of the genome-about 1.8%, encompassing more than 600 genes-differentiated stickleback from the two biotic environments. Divergence occurred in 141 discrete genomic clumps located mainly in regions of low recombination, suggesting that genes brought to lakes by the colonizing ancestral population often evolved together in linked blocks. Strong selection and a wealth of standing genetic variation explain how a single member of the biotic community can have such a rapid and profound evolutionary impact.

The genomic basis of Red Queen dynamics during rapid reciprocal host-pathogen coevolution

Pathogens are omnipresent and by definition detrimental to their hosts. Pathogens thus exert high selection on their hosts, which, if adapting, can exert similar levels of selection on the pathogen, resulting in ongoing cycles of reciprocal adaptation between the antagonists. Such coevolutionary interactions have a central influence on the evolution of organisms. Surprisingly, we still know little about the exact selection dynamics and the genome regions involved. Our study uses a controlled experimental approach with an animal host to dissect coevolutionary selection. We find that distinct selective processes underlie rapid coadaptation in the two antagonists, including antagonistic frequency-dependent selection on toxin gene copy number in the pathogen, while the host response is likely influenced by changes in multiple genome regions.

Red Queen dynamics, involving coevolutionary interactions between species, are ubiquitous, shaping the evolution of diverse biological systems. To date, information on the underlying selection dynamics and the involved genome regions is mainly available for bacteria–phage systems or only one of the antagonists of a eukaryotic host–pathogen interaction. We add to our understanding of these important coevolutionary interactions using an experimental host–pathogen model, which includes the nematode Caenorhabditis elegans and its pathogen Bacillus thuringiensis. We combined experimental evolution with time-shift experiments, in which a focal host or pathogen is tested against a coevolved antagonist from the past, present, or future, followed by genomic analysis. We show that (i) coevolution occurs rapidly within few generations, (ii) temporal coadaptation at the phenotypic level is found in parallel across replicate populations, consistent with antagonistic frequency-dependent selection, (iii) genomic changes in the pathogen match the phenotypic pattern and include copy number variations of a toxin-encoding plasmid, and (iv) host genomic changes do not match the phenotypic pattern and likely involve selective responses at more than one locus. By exploring the dynamics of coevolution at the phenotypic and genomic level for both host and pathogen simultaneously, our findings demonstrate a more complex model of the Red Queen, consisting of distinct selective processes acting on the two antagonists during rapid and reciprocal coadaptation.

On the origins and domestication of the olive: a review and perspectives

Background

Unravelling domestication processes is crucial for understanding how species respond to anthropogenic pressures, forecasting crop responses to future global changes and improving breeding programmes. Domestication processes for clonally propagated perennials differ markedly from those for seed-propagated annual crops, mostly due to long generation times, clonal propagation and recurrent admixture with local forms, leading to a limited number of generations of selection from wild ancestors. However, additional case studies are required to document this process more fully.

Scope

The olive is an iconic species in Mediterranean cultural history. Its multiple uses and omnipresence in traditional agrosystems have made this species an economic pillar and cornerstone of Mediterranean agriculture. However, major questions about the domestication history of the olive remain unanswered. New paleobotanical, archeological, historical and molecular data have recently accumulated for olive, making it timely to carry out a critical re-evaluation of the biogeography of wild olives and the history of their cultivation. We review here the chronological history of wild olives and discuss the questions that remain unanswered, or even unasked, about their domestication history in the Mediterranean Basin. We argue that more detailed ecological genomics studies of wild and cultivated olives are crucial to improve our understanding of olive domestication. Multidisciplinary research integrating genomics, metagenomics and community ecology will make it possible to decipher the evolutionary ecology of one of the most iconic domesticated fruit trees worldwide.

Conclusion

The olive is a relevant model for improving our knowledge of domestication processes in clonally propagated perennial crops, particularly those of the Mediterranean Basin. Future studies on the ecological and genomic shifts linked to domestication in olive and its associated community will provide insight into the phenotypic and molecular bases of crop adaptation to human uses.

Parasites driving host diversity: Incidence of disease correlated with Daphnia clonal turnover

According to the Red Queen hypothesis, clonal diversity in asexual populations could be maintained by negative frequency-dependant selection by coevolving parasites. If common clones are selected against and rare clones gain a concomitant advantage, we expect that clonal turnover should be faster during parasite epidemics than between them. We tested this hypothesis exploring field data of the Daphnia-Caullerya host-parasite system. The clonal make-up and turnover of the Daphnia host population was tracked with high temporal resolution from 1998 until 2013, using first allozyme and later microsatellite markers. Significant differences in the clonal composition between random and infected subsamples of Daphnia populations were detected on six of seven tested occasions, confirming genetic specificity of the host-parasite interaction in this system. We used time series analysis to compare the rates of host clonal turnover to the incidence of parasitism, and found that Caullerya prevalence was significantly associated with microsatellite-based clonal turnover. As alternate hypotheses, we further tested whether turnover was related to a variety of biotic, abiotic, and host demographic parameters. Other significant correlates of turnover were cyanobacterial biomass and (weakly) temperature. Overall, parasitism seems to be a strong driver of host clonal turnover, in support of the Red Queen hypothesis.

Back to the future in a petri dish: Origin and impact of resurrected microbes in natural populations

Current natural populations face new interactions because of the re-emergence of ancient microbes and viruses. These risks come from the re-emergence of pathogens kept in laboratories or from pathogens that are retained in the permafrost, which become available upon thawing due to climate change. We here focus on the effects of such re-emergence in natural host populations based on evolutionary theory of virulence and long-term studies, which investigate host-pathogen adaptations. Pathogens tend to be locally and temporally adapted to their co-occurring hosts, but when pathogens from a different environment or different time enter the host community, the degree to which a new host-pathogen interaction is a threat will depend on the specific genotypic associations, the time lag between the host and the pathogen, and the interactions with native or recent host and pathogen species. Some insights can be obtained from long-term studies using a resurrection ecology approach. These long-term studies based on time-shift experiments are essential to obtain insight into the mechanisms underlying host-pathogen coevolution at several ecological and temporal scales. As past pathogens and their corresponding host(s) can differ in infectivity and susceptibility, strong reciprocal selective pressures can be induced by the pathogen. These strong selective pressures often result in an escalating arms race, but do not necessarily result in increased infectivity over time. Human health can also be impacted by these resurrected pathogens as the majority of emerging infectious diseases are zoonoses, which are infectious diseases originating from animal populations naturally transmitted to humans. The sanitary risk associated with pathogen emergence from different environments (spatial or temporal) depends on a combination of socioeconomic, environmental, and ecological factors that affect the virulence or the pathogenic potential of microbes and their ability to infect susceptible host populations.

Keeping Pace with the Red Queen: Identifying the Genetic Basis of Susceptibility to Infectious Disease

The results of genome-wide association studies are known to be affected by epistasis and gene-by-environment interactions. Using a statistical model....

Genome-wide association studies are widely used to identify “disease genes” conferring resistance/susceptibility to infectious diseases. Using a combination of mathematical models and simulations, we demonstrate that genetic interactions between hosts and parasites [genotype-by-genotype (G × G) interactions] can drastically affect the results of these association scans and hamper our ability to detect genetic variation in susceptibility. When hosts and parasites coevolve, these G × G interactions often make genome-wide association studies unrepeatable over time or across host populations. Reanalyzing previously published data on Daphnia magna susceptibility to infection by Pasteuria ramosa, we identify genomic regions consistent with G × G interactions. We conclude by outlining possible avenues for designing more powerful and more repeatable association studies.

Environmental variation causes different (co) evolutionary routes to the same adaptive destination across parasite populations

Epidemics are engines for host-parasite coevolution, where parasite adaptation to hosts drives reciprocal adaptation in host populations. A key challenge is to understand whether parasite adaptation and any underlying evolution and coevolution is repeatable across ecologically realistic populations that experience different environmental conditions, or if each population follows a completely unique evolutionary path. We established twenty replicate pond populations comprising an identical suite of genotypes of crustacean host, Daphnia magna, and inoculum of their parasite, Pasteuria ramosa. Using a time-shift experiment, we compared parasite infection traits before and after epidemics and linked patterns of parasite evolution with shifts in host genotype frequencies. Parasite adaptation to the sympatric suite of host genotypes came at a cost of poorer performance on foreign genotypes across populations and environments. However, this consistent pattern of parasite adaptation was driven by different types of frequency-dependent selection that was contingent on an ecologically relevant environmental treatment (whether or not there was physical mixing of water within ponds). In unmixed ponds, large epidemics drove rapid and strong host-parasite coevolution. In mixed ponds, epidemics were smaller and host evolution was driven mainly by the mixing treatment itself; here, host evolution and parasite evolution were clear, but coevolution was absent. Population mixing breaks an otherwise robust coevolutionary cycle. These findings advance our understanding of the repeatability of (co)evolution across noisy, ecologically realistic populations.

QTL mapping of a natural genetic polymorphism for long-term parasite persistence in Daphnia populations

Knowing the determinants of the geographic ranges of parasites is important for understanding their evolutionary ecology, epidemiology and their potential to expand their range. Here we explore the determinants of geographic range in the peculiar case of a parasite species – the microsporidian Hamiltosporidium tvaerminnensis – that has a limited geographic distribution in a wide-spread host – Daphnia magna. We conducted a quantitative trait loci (QTLs) analysis with monoclonal F2D. magna populations originating from a cross between a susceptible northern European genotype and a resistant central European genotype. Contrary to our expectations, long-term persistence turned out to be a quantitative trait across the F2 offspring. Evidence for two QTLs, one epistatic interaction and for further minor QTL was found. This finding contrasts markedly with the previously described bimodal pattern for long-term parasite persistence in natural host genotypes across Europe and leaves open the question of how a quantitative genetic trait could determine the disjunct geographic distribution of the parasite across Europe.