The Red Queen lives: Epistasis between linked resistance loci

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.

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.

Long-read sequencing reveals the evolutionary drivers of intra-host diversity across natural RNA mycovirus infections

Intra-host dynamics are a core component of virus evolution but most intra-host data come from a narrow range of hosts or experimental infections. Gaining broader information on the intra-host diversity and dynamics of naturally occurring virus infections is essential to our understanding of evolution across the virosphere. Here we used PacBio long-read HiFi sequencing to characterize the intra-host populations of natural infections of the RNA mycovirus Cryphonectria hypovirus 1 (CHV1). CHV1 is a biocontrol agent for the chestnut blight fungus (Cryphonectria parasitica), which co-invaded Europe alongside the fungus. We characterized the mutational and haplotypic intra-host virus diversity of thirty-eight natural CHV1 infections spread across four locations in Croatia and Switzerland. Intra-host CHV1 diversity values were shaped by purifying selection and accumulation of mutations over time as well as epistatic interactions within the host genome at defense loci. Geographical landscape features impacted CHV1 inter-host relationships through restricting dispersal and causing founder effects. Interestingly, a small number of intra-host viral haplotypes showed high sequence similarity across large geographical distances unlikely to be linked by dispersal.

Bringing Disciplines and People Together to Characterize the Plastic and Genetic Responses of Molluscs to Environmental Change

Molluscs are remarkably diverse and are found across nearly all ecosystems, meaning that members of this ancient animal phylum provide a powerful means to study genomic-phenotype connections in a climate change framework. Recent advances in genomic sequencing technologies and genome assembly approaches finally allow the relatively cheap and tractable assembly of high-quality mollusc genome resources. After a brief review of these issues and advances, we use a case-study approach to provide some concrete examples of phenotypic plasticity and genomic adaptation in molluscs in response to environmental factors expected to be influenced by climate change. Our goal is to use molluscs as a "common currency" to demonstrate how organismal and evolutionary biologists can use natural systems to make phenotype-genotype connections in the context of changing environments. In parallel, we emphasize the critical need to collaborate and integrate findings across taxa and disciplines in order to use new data and information to advance our understanding of mollusc biology in the context of global environmental change. We end with a brief synthetic summary of the papers inspired by the 2021 SICB Symposium "Genomic Perspectives in Comparative Physiology of Molluscs: Integration across Disciplines".

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.

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.

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.

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.

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.

Rapid change in host specificity in a field population of the biological control organism Pasteuria penetrans

In biological control, populations of both the biological control agent and the pest have the potential to evolve and even to coevolve. This feature marks the most powerful and unpredictable aspect of biological control strategies. In particular, evolutionary change in host specificity of the biological control agent could increase or decrease its efficacy. Here, we tested for change in host specificity in a field population of the biological control organism Pasteuria penetrans. Pasteuria penetrans is an obligate parasite of the plant parasitic nematodes Meloidogyne spp., which are major agricultural pests. From 2013 through 2016, we collected yearly samples of P. penetrans from eight plots in a field infested with M. arenaria. Plots were planted either with peanut (Arachis hypogaea) or with a rotation of peanut and soybean (Glycine max). To detect temporal change in host specificity, we tested P. penetrans samples annually for their ability to attach to (and thereby infect) four clonal lines of M. arenaria. After controlling for temporal variation in parasite abundance, we found that P. penetrans from each of the eight plots showed temporal variation in their attachment specificity to the clonal host lines. The trajectories of change in host specificity were largely unique to each plot. This result suggests that local forces, at the level of individual plots, drive change in specificity. We hypothesize that coevolution with local M. arenaria hosts may be one such force. Lastly, we observed an overall reduction in attachment rate with samples from rotation plots relative to samples from peanut plots. This result may reflect lower abundance of P. penetrans under crop rotation, potentially due to suppressed density of host nematodes. As a whole, the results show local change in specificity on a yearly basis, consistent with evolution of a biological control organism in its ability to infect and suppress its target pest.

Assortment of Flowering Time and Immunity Alleles in Natural Arabidopsis thaliana Populations Suggests Immunity and Vegetative Lifespan Strategies Coevolve

The selective impact of pathogen epidemics on host defenses can be strong but remains transient. By contrast, life-history shifts can durably and continuously modify the balance between costs and benefits of immunity, which arbitrates the evolution of host defenses. Their impact on the evolutionary dynamics of host immunity, however, has seldom been documented. Optimal investment into immunity is expected to decrease with shortening lifespan, because a shorter life decreases the probability to encounter pathogens or enemies. Here, we document that in natural populations of Arabidopsis thaliana, the expression levels of immunity genes correlate positively with flowering time, which in annual species is a proxy for lifespan. Using a novel genetic strategy based on bulk-segregants, we partitioned flowering time-dependent from -independent immunity genes and could demonstrate that this positive covariation can be genetically separated. It is therefore not explained by the pleiotropic action of some major regulatory genes controlling both immunity and lifespan. Moreover, we find that immunity genes containing variants reported to impact fitness in natural field conditions are among the genes whose expression covaries most strongly with flowering time. Taken together, these analyses reveal that natural selection has likely assorted alleles promoting lower expression of immunity genes with alleles that decrease the duration of vegetative lifespan in A. thaliana and vice versa. This is the first study documenting a pattern of variation consistent with the impact that selection on flowering time is predicted to have on diversity in host immunity.

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.

Chaotic provinces in the kingdom of the Red Queen

The interplay between parasites and their hosts is found in all kinds of species and plays an important role in understanding the principles of evolution and coevolution. Usually, the different genotypes of hosts and parasites oscillate in their abundances. The well-established theory of oscillatory Red Queen dynamics proposes an ongoing change in frequencies of the different types within each species. So far, it is unclear under what conditions Red Queen dynamics persists, especially when the number of types per species increases. Some models show that with many types of hosts and parasites or more species chaotic dynamics occur. In our analysis, an arbitrary number of types within two species are examined in a deterministic framework with constant or changing population size and very simple interactions. This general framework allows for analytical solutions for internal fixed points and their stability. The numerical analysis shows that for two species, once more than two types are considered per species, irregular dynamics in their frequencies can be observed in the long run. The nature of the dynamics depends strongly on the initial configuration of the system; the usual regular Red Queen oscillations are only observed when all types initially have similar abundance.

The Evolutionary Consequences of Stepwise Infection Processes

Molecular and cellular studies reveal that the resistance of hosts to parasites and pathogens is a cascade-like process with multiple steps required to be passed for successful infection. By contrast, much of evolutionary reasoning is based on strongly simplified, one- or two-step infection processes with simple genetics or on resistance being a quantitative trait. Here we attempt a conceptual unification of these two perspectives with the aim of cross-fostering research and filling some of the gaps in our concepts of the ecology and evolution of disease. This conceptual unification has a profound impact on the way we understand the genetics and evolution of host resistance, ecological immunity, evolution of virulence, defence portfolios, and host-pathogen coevolution.

Hill-Robertson interference maintained by Red Queen dynamics favours the evolution of sex

Although it is well established theoretically that selective interference among mutations (Hill-Robertson interference) favours meiotic recombination, genomewide mean rates of mutation and strengths of selection appear too low to support this as the mechanism favouring recombination in nature. A possible solution to this discrepancy between theory and observation is that selection is at least intermittently very strong due to the antagonistic coevolution between a host and its parasites. The Red Queen theory posits that such coevolution generates fitness epistasis among loci, which generates negative linkage disequilibrium among beneficial mutations, which in turn favours recombination. This theory has received only limited support. However, Red Queen dynamics without epistasis may provide the ecological conditions that maintain strong and frequent selective interference in finite populations that indirectly selects for recombination. This hypothesis is developed here through the simulation of Red Queen dynamics. This approach required the development of a method to calculate the exact frequencies of multilocus haplotypes after recombination. Simulations show that recombination is favoured by the moderately weak selection of many loci involved in the interaction between a host and its parasites, which results in substitution rates that are compatible with empirical estimates. The model also reproduces the previously reported rapid increase in the rate of outcrossing in Caenorhabditis elegans coevolving with a bacterial pathogen.

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

Negative frequency-dependent selection (NFDS) is an evolutionary mechanism suggested to govern host-parasite coevolution and the maintenance of genetic diversity at host resistance loci, such as the vertebrate MHC and R-genes in plants. Matching-allele interactions of hosts and parasites that prevent the emergence of host and parasite genotypes that are universally resistant and infective are a genetic mechanism predicted to underpin NFDS. The underlying genetics of matching-allele interactions are unknown even in host-parasite systems with empirical support for coevolution by NFDS, as is the case for the planktonic crustacean Daphnia magna and the bacterial pathogen Pasteuria ramosa. We fine-map one locus associated with D. magna resistance to P. ramosa and genetically characterize two haplotypes of the Pasteuria resistance (PR-) locus using de novo genome and transcriptome sequencing. Sequence comparison of PR-locus haplotypes finds dramatic structural polymorphisms between PR-locus haplotypes including a large portion of each haplotype being composed of non-homologous sequences resulting in haplotypes differing in size by 66 kb. The high divergence of PR-locus haplotypes suggest a history of multiple, diverse and repeated instances of structural mutation events and restricted recombination. Annotation of the haplotypes reveals striking differences in gene content. In particular, a group of glycosyltransferase genes that is present in the susceptible but absent in the resistant haplotype. Moreover, in natural populations, we find that the PR-locus polymorphism is associated with variation in resistance to different P. ramosa genotypes, pointing to the PR-locus polymorphism as being responsible for the matching-allele interactions that have been previously described for this system. Our results conclusively identify a genetic basis for the matching-allele interaction observed in a coevolving host-parasite system and provide a first insight into its molecular basis.

Higher rates of sex evolve during adaptation to more complex environments

A leading hypothesis for the evolutionary maintenance of sexual reproduction proposes that sex is advantageous because it facilitates adaptation. Changes in the environment stimulate adaptation but not all changes are equivalent; a change may occur along one or multiple environmental dimensions. In two evolution experiments with the facultatively sexual rotifer Brachionus calyciflorus, we test how environmental complexity affects the evolution of sex by adapting replicate populations to various environments that differ from the original along one, two, or three environmental dimensions. Three different estimates of fitness (growth, lifetime reproduction, and population density) confirmed that populations adapted to their new environment. Growth measures revealed an intriguing cost of complex adaptations: populations that adapted to more complex environments lost greater amounts of fitness in the original environment. Furthermore, both experiments showed that B. calyciflorus became more sexual when adapting to a greater number of environmental dimensions. Common garden experiments confirmed that observed changes in sex were heritable. As environments in nature are inherently complex these findings help explain why sex is maintained in natural populations.