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