Movement can mediate temporal mismatches between resource availability and biological events in host–pathogen interactions

Global change is shifting the timing of biological events, leading to temporal mismatches between biological events and resource availability. These temporal mismatches can threaten species’ populations. Importantly, temporal mismatches not only exert strong pressures on the population dynamics of the focal species, but can also lead to substantial changes in pairwise species interactions such as host–pathogen systems. We adapted an established individual‐based model of host–pathogen dynamics. The model describes a viral agent in a social host, while accounting for the host''s explicit movement decisions. We aimed to investigate how temporal mismatches between seasonal resource availability and host life‐history events affect host–pathogen coexistence, that is, disease persistence. Seasonal resource fluctuations only increased coexistence probability when in synchrony with the hosts’ biological events. However, a temporal mismatch reduced host–pathogen coexistence, but only marginally. In tandem with an increasing temporal mismatch, our model showed a shift in the spatial distribution of infected hosts. It shifted from an even distribution under synchronous conditions toward the formation of disease hotspots, when host life history and resource availability mismatched completely. The spatial restriction of infected hosts to small hotspots in the landscape initially suggested a lower coexistence probability due to the critical loss of susceptible host individuals within those hotspots. However, the surrounding landscape facilitated demographic rescue through habitat‐dependent movement. Our work demonstrates that the negative effects of temporal mismatches between host resource availability and host life history on host–pathogen coexistence can be reduced through the formation of temporary disease hotspots and host movement decisions, with implications for disease management under disturbances and global change.     

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.     

Sexual antagonism for resistance and tolerance to infection in Drosophila melanogaster

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

Population mixing promotes arms race host–parasite coevolution

The consequences of host–parasite coevolution are highly contingent on the qualitative coevolutionary dynamics: whether selection fluctuates (fluctuating selection dynamic; FSD), or is directional towards increasing infectivity/resistance (arms race dynamic; ARD). Both genetics and ecology can play an important role in determining whether coevolution follows FSD or ARD, but the ecological conditions under which FSD shifts to ARD, and vice versa, are not well understood. The degree of population mixing is thought to increase host exposure to parasites, hence selecting for greater resistance and infectivity ranges, and we hypothesize this promotes ARD. We tested this by coevolving bacteria and viruses in soil microcosms and found that population mixing shifted bacteria–virus coevolution from FSD to ARD. A simple theoretical model produced qualitatively similar results, showing that mechanisms that increase host exposure to parasites tend to push dynamics towards ARD. The shift from FSD to ARD with increased population mixing may help to explain variation in coevolutionary dynamics between different host–parasite systems, and more specifically the observed discrepancies between laboratory and field bacteria–virus coevolutionary studies.     

Evolution of parasitoid host preference and performance in response to an invasive host acting as evolutionary trap

The invasion of a novel host species can create a mismatch in host choice and offspring survival (performance) when native parasitoids attempt to exploit the invasive host without being able to circumvent its resistance mechanisms. Invasive hosts can therefore act as evolutionary trap reducing parasitoids'' fitness and this may eventually lead to their extinction. Yet, escape from the trap can occur when parasitoids evolve behavioral avoidance or a physiological strategy compatible with the trap host, resulting in either host‐range expansion or a complete host‐shift. We developed an individual based model to investigate which conditions promote parasitoids to evolve behavioral preference that matches their performance, including host‐trap avoidance, and which conditions lead to adaptations to the unsuitable hosts. The model was inspired by solitary endo‐parasitoids attacking larval host stages. One important aspect of these conditions was reduced host survival during incompatible interaction, where a failed parasitization attempt by a parasitoid resulted not only in death of her offspring but also in host killing. This non‐reproductive host mortality had a strong influence on the likelihood of establishment of novel host–parasitoid relationship, in some cases constraining adaptation to the trap host species. Moreover, our model revealed that host‐search efficiency and genetic variation in host‐preference play a key role in the likelihood that parasitoids will include the suboptimal host in their host range, or will evolve behavioral avoidance resulting in specialization and host‐range conservation, respectively. Hence, invasive species might change the evolutionary trajectory of native parasitoid species, which is important for predicting biocontrol ability of native parasitoids towards novel hosts.     

Adaptive population structure shifts in invasive parasitic mites,Varroa destructor

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

Host phenology regulates parasite–host demographic cycles and eco‐evolutionary feedbacks

Parasite–host interactions can drive periodic population dynamics when parasites overexploit host populations. The timing of host seasonal activity, or host phenology, determines the frequency and demographic impact of parasite–host interactions, which may govern whether parasites sufficiently overexploit hosts to drive population cycles. We describe a mathematical model of a monocyclic, obligate‐killer parasite system with seasonal host activity to investigate the consequences of host phenology on host–parasite dynamics. The results suggest that parasites can reach the densities necessary to destabilize host dynamics and drive cycling as they adapt, but only in some phenological scenarios such as environments with short seasons and synchronous host emergence. Furthermore, only parasite lineages that are sufficiently adapted to phenological scenarios with short seasons and synchronous host emergence can achieve the densities necessary to overexploit hosts and produce population cycles. Host‐parasite cycles also generate an eco‐evolutionary feedback that slows parasite adaptation to the phenological environment as rare advantageous phenotypes can be driven extinct due to a population bottleneck depending on when they are introduced in the cycle. The results demonstrate that seasonal environments can drive population cycling in a restricted set of phenological patterns and provide further evidence that the rate of adaptive evolution depends on underlying ecological dynamics.     

The drivers of avian‐haemosporidian prevalence in tropical lowland forests of New Guinea in three dimensions

Haemosporidians are among the most common parasites of birds and often negatively impact host fitness. A multitude of biotic and abiotic factors influence these associations, but the magnitude of these factors can differ by spatial scales (i.e., local, regional and global). Consequently, to better understand global and regional drivers of avian‐haemosporidian associations, it is key to investigate these associations at smaller (local) spatial scales. Thus, here, we explore the effect of abiotic variables (e.g., temperature, forest structure, and anthropogenic disturbances) on haemosporidian prevalence and host–parasite networks on a horizontal spatial scale, comparing four fragmented forests and five localities within a continuous forest in Papua New Guinea. Additionally, we investigate if prevalence and host–parasite networks differ between the canopy and the understory (vertical stratification) in one forest patch. We found that the majority of Haemosporidian infections were caused by the genus Haemoproteus and that avian‐haemosporidian networks were more specialized in continuous forests. At the community level, only forest greenness was negatively associated with Haemoproteus infections, while the effects of abiotic variables on parasite prevalence differed between bird species. Haemoproteus prevalence levels were significantly higher in the canopy, and an opposite trend was observed for Plasmodium. This implies that birds experience distinct parasite pressures depending on the stratum they inhabit, likely driven by vector community differences. These three‐dimensional spatial analyses of avian‐haemosporidians at horizontal and vertical scales suggest that the effect of abiotic variables on haemosporidian infections are species specific, so that factors influencing community‐level infections are primarily driven by host community composition.     

Indirect effects of parasitism: costs of infection to other individuals can be greater than direct costs borne by the host

Parasitic infection has a direct physiological cost to hosts but may also alter how hosts interact with other individuals in their environment. Such indirect effects may alter both host fitness and the fitness of other individuals in the host''s social network, yet the relative impact of direct and indirect effects of infection are rarely quantified. During reproduction, a host''s social environment includes family members who may be in conflict over resource allocation. In such situations, infection may alter how resources are allocated, thereby redistributing the costs of parasitism between individuals. Here, we experimentally reduce parasite burdens of parent and/or nestling European shags (Phalacrocorax aristotelis) infected with Contracaecum nematodes in a factorial design, then simultaneously measure the impact of an individual''s infection on all family members. We found no direct effect of infection on parent or offspring traits but indirect effects were detected in all group members, with both immediate effects (mass change and survival) and longer-term effects (timing of parents’ subsequent breeding). Our results show that parasite infection can have a major impact on individuals other than the host, suggesting that the effect of parasites on population processes may be greater than previously thought.     

Using Microbial Microcosms to Study Host–parasite Coevolution

Microbial microcosm experiments with bacteria and their viral parasites allow us to observe host–parasite coevolution in action. Laboratory populations of microbes evolve rapidly, thanks to their short generation times and huge population sizes. By taking advantage of a “living fossil record” stored in the laboratory freezer, we can directly compare the fitness of hosts and parasites with their actual evolutionary ancestors. Such experiments demonstrate that host–parasite coevolution is an important evolutionary force and a cause of strong and divergent natural selection.     

Genetic architecture of resistance in Daphnia hosts against two species of host-specific parasites

Understanding the genetic architecture of host resistance is key for understanding the evolution of host–parasite interactions. Evolutionary models often assume simple genetics based on few loci and strong epistasis. It is unknown, however, whether these assumptions apply to natural populations. Using a quantitative trait loci (QTL) approach, we explore the genetic architecture of resistance in the crustacean Daphnia magna to two of its natural parasites: the horizontally transmitted bacterium Pasteuria ramosa and the horizontally and vertically transmitted microsporidium Hamiltosporidium tvaerminnensis. These two systems have become models for studies on the evolution of host–parasite interactions. In the QTL panel used here, Daphnia''s resistance to P. ramosa is controlled by a single major QTL (which explains 50% of the observed variation). Resistance to H. tvaerminnensis horizontal infections shows a signature of a quantitative trait based in multiple loci with weak epistatic interactions (together explaining 38% variation). Resistance to H. tvaerminnensis vertical infections, however, shows only one QTL (explaining 13.5% variance) that colocalizes with one of the QTLs for horizontal infections. QTLs for resistance to Pasteuria and Hamiltosporidium do not colocalize. We conclude that the genetics of resistance in D. magna are drastically different for these two parasites. Furthermore, we infer that based on these and earlier results, the mechanisms of coevolution differ strongly for the two host–parasite systems. Only the Pasteuria–Daphnia system is expected to follow the negative frequency-dependent selection (Red Queen) model. How coevolution works in the Hamiltosporidium–Daphnia system remains unclear.     

Biological control in a disturbed environment

Most ecological and epidemiological models describe systems with continuous uninterrupted interactions between populations. Many systems, though, have ecological disturbances, such as those associated with planting and harvesting of a seasonal crop. In this paper, we introduce host–parasite–hyperparasite systems as models of biological control in a disturbed environment, where the host–parasite interactions are discontinuous. One model is a parasite–hyperparasite system designed to capture the essence of biological control and the other is a host–parasite–hyperparasite system that incorporates many more features of the population dynamics. Two types of discontinuity are included in the models. One corresponds to a pulse of new parasites at harvest and the other reflects the discontinuous presence of the host due to planting and harvesting. Such discontinuities are characteristic of many ecosystems involving parasitism or other interactions with an annual host. The models are tested against data from an experiment investigating the persistent biological control of the fungal plant parasite of lettuce Sclerotinia minor by the fungal hyperparasite Sporidesmium sclerotivorum, over successive crops. Using a combination of mathematical analysis, model fitting and parameter estimation, the factors that contribute the observed persistence of the parasite are examined. Analytical results show that repeated planting and harvesting of the host allows the parasite to persist by maintaining a quantity of host tissue in the system on which the parasite can reproduce. When the host dynamics are not included explicitly in the model, we demonstrate that homogeneous mixing fails to predict the persistence of the parasite population, while incorporating spatial heterogeneity by allowing for heterogeneous mixing prevents fade-out. Including the host''s dynamics lessens the effect of heterogeneous mixing on persistence, though the predicted values for the parasite population are closer to the observed values. An alternative hypothesis for persistence involving a stepped change in rates of infection is also tested and model fitting is used to show that changes in some environmental conditions may contribute to parasite persistence. The importance of disturbances and periodic forcing in models for interacting populations is discussed.     

Evolution towards oscillation or stability in a predator–prey system

We studied a prey–predator system in which both species evolve. We discuss here the conditions that result in coevolution towards a stable equilibrium or towards oscillations. First, we show that a stable equilibrium or population oscillations with small amplitude is likely to occur if the prey''s (host''s) defence is effective when compared with the predator''s (parasite''s) attacking ability at equilibrium, whereas large-amplitude oscillations are likely if the predator''s (parasite''s) attacking ability exceeds the prey''s (host''s) defensive ability. Second, a stable equilibrium is more likely if the prey''s defensive trait evolves faster than the predator''s attack trait, whereas population oscillations are likely if the predator''s trait evolves faster than that of the prey. Third, when the adaptation rates of both species are similar, the amplitude of the fluctuations in their abundances is small when the adaptation rate is either very slow or very fast, but at an intermediate rate of adaptation the fluctuations have a large amplitude. We also show the case in which the prey''s abundance and trait fluctuate greatly, while those of the predator remain almost unchanged. Our results predict that populations and traits in host–parasite systems are more likely than those in prey–predator systems to show large-amplitude oscillations.     

Quantifying Adaptive Evolution in the Drosophila Immune System

It is estimated that a large proportion of amino acid substitutions in Drosophila have been fixed by natural selection, and as organisms are faced with an ever-changing array of pathogens and parasites to which they must adapt, we have investigated the role of parasite-mediated selection as a likely cause. To quantify the effect, and to identify which genes and pathways are most likely to be involved in the host–parasite arms race, we have re-sequenced population samples of 136 immunity and 287 position-matched non-immunity genes in two species of Drosophila. Using these data, and a new extension of the McDonald-Kreitman approach, we estimate that natural selection fixes advantageous amino acid changes in immunity genes at nearly double the rate of other genes. We find the rate of adaptive evolution in immunity genes is also more variable than other genes, with a small subset of immune genes evolving under intense selection. These genes, which are likely to represent hotspots of host–parasite coevolution, tend to share similar functions or belong to the same pathways, such as the antiviral RNAi pathway and the IMD signalling pathway. These patterns appear to be general features of immune system evolution in both species, as rates of adaptive evolution are correlated between the D. melanogaster and D. simulans lineages. In summary, our data provide quantitative estimates of the elevated rate of adaptive evolution in immune system genes relative to the rest of the genome, and they suggest that adaptation to parasites is an important force driving molecular evolution.     

Genotype specificity among hosts,pathogens, and beneficial microbes influences the strength of symbiont‐mediated protection

The microbial symbionts of eukaryotes influence disease resistance in many host‐parasite systems. Symbionts show substantial variation in both genotype and phenotype, but it is unclear how natural selection maintains this variation. It is also unknown whether variable symbiont genotypes show specificity with the genotypes of hosts or parasites in natural populations. Genotype by genotype interactions are a necessary condition for coevolution between interacting species. Uncovering the patterns of genetic specificity among hosts, symbionts, and parasites is therefore critical for determining the role that symbionts play in host‐parasite coevolution. Here, we show that the strength of protection conferred against a fungal pathogen by a vertically transmitted symbiont of an aphid is influenced by both host‐symbiont and symbiont‐pathogen genotype by genotype interactions. Further, we show that certain symbiont phylogenetic clades have evolved to provide stronger protection against particular pathogen genotypes. However, we found no evidence of reciprocal adaptation of co‐occurring host and symbiont lineages. Our results suggest that genetic variation among symbiont strains may be maintained by antagonistic coevolution with their host and/or their host's parasites.     

Developmental timing of extreme temperature events (heat waves) disrupts host–parasitoid interactions

1. When thermal tolerances differ between interacting species, extreme temperature events (heat waves) will alter the ecological outcomes. The parasitoid wasp Cotesia congregata suffers high mortality when reared throughout development at temperatures that are nonstressful for its host, Manduca sexta. However, the effects of short‐term heat stress during parasitoid development are unknown in this host–parasitoid system.
2. Here, we investigate how duration of exposure, daily maximum temperature, and the developmental timing of heat waves impact the performance of C. congregata and its host¸ M. sexta. We find that the developmental timing of short‐term heat waves strongly determines parasitoid and host outcomes.
3. Heat waves during parasitoid embryonic development resulted in complete wasp mortality and the production of giant, long‐lived hosts. Heat waves during the 1st‐instar had little effect on wasp success, whereas heat waves during the parasitoid''s nutritionally and hormonally critical 2nd instar greatly reduced wasp emergence and eclosion. The temperature and duration of heat waves experienced early in development determined what proportion of hosts had complete parasitoid mortality and abnormal phenotypes.
4. Our results suggest that the timing of extreme temperature events will be crucial to determining the ecological impacts on this host–parasitoid system. Discrepancies in thermal tolerance between interacting species and across development will have important ramifications on ecosystem responses to climate change.

     

Cross‐continental comparison of parasite communities in a wide‐ranging carnivore suggests associations with prey diversity and host density

1. Parasites are integral to ecosystem functioning yet often overlooked. Improved understanding of host–parasite associations is important, particularly for wide‐ranging species for which host range shifts and climate change could alter host–parasite interactions and their effects on ecosystem function.
2. Among the most widely distributed mammals with diverse diets, gray wolves (Canis lupus) host parasites that are transmitted among canids and via prey species. Wolf–parasite associations may therefore influence the population dynamics and ecological functions of both wolves and their prey. Our goal was to identify large‐scale processes that shape host–parasite interactions across populations, with the wolf as a model organism.
3. By compiling data from various studies, we examined the fecal prevalence of gastrointestinal parasites in six wolf populations from two continents in relation to wolf density, diet diversity, and other ecological conditions.
4. As expected, we found that the fecal prevalence of parasites transmitted directly to wolves via contact with other canids or their excreta was positively associated with wolf density. Contrary to our expectations, the fecal prevalence of parasites transmitted via prey was negatively associated with prey diversity. We also found that parasite communities reflected landscape characteristics and specific prey items available to wolves.
5. Several parasite taxa identified in this study, including hookworms and coccidian protozoans, can cause morbidity and mortality in canids, especially in pups, or in combination with other stressors. The density–prevalence relationship for parasites with simple life cycles may reflect a regulatory role of gastrointestinal parasites on wolf populations. Our result that fecal prevalence of parasites was lower in wolves with more diverse diets could provide insight into the mechanisms by which biodiversity may regulate disease. A diverse suite of predator–prey interactions could regulate the effects of parasitism on prey populations and mitigate the transmission of infectious agents, including zoonoses, spread via trophic interactions.

     

Network robustness and structure depend on the phenological characteristics of plants and pollinators

Many structural patterns have been found to be important for the stability and robustness of mutualistic plant–pollinator networks. These structural patterns are impacted by a suite of variables, including species traits, species abundances, their spatial configuration, and their phylogenetic history. Here, we consider a specific trait: phenology, or the timing of life history events. We expect that timing and duration of activity of pollinators, or of flowering in plants, could greatly affect the species'' roles within networks in which they are embedded. Using plant–pollinator networks from 33 sites in southern British Columbia, Canada, we asked (a) how phenological species traits, specifically timing of first appearance in the network and duration of activity in a network, were related to species'' roles within a network, and (b) how those traits affected network robustness to phenologically biased species loss. We found that long duration of activity increased connection within modules for both pollinators and plants and among modules for plants. We also found that date of first appearance was positively related to interaction strength asymmetry in plants but negatively related to pollinators. Networks were generally more robust to the loss of pollinators than plants, and robustness increased if the models allow new interactions to form when old ones are lost, constrained by overlapping phenology of plants and pollinators. Robustness declined with the loss of late‐flowering plants, which tended to have higher interaction strength asymmetry. In addition, robustness declined with loss of early‐flying or long‐duration pollinators. These pollinators tended to be among‐module connectors. Our results point to networks being limited by early‐flying pollinators. If plants flower earlier due to climate change, plant fitness may decline as they will depend on early emerging pollinators, unless pollinators also emerge earlier.     

Starvation reveals the cause of infection-induced castration and gigantism

Parasites often induce life-history changes in their hosts. In many cases, these infection-induced life-history changes are driven by changes in the pattern of energy allocation and utilization within the host. Because these processes will affect both host and parasite fitness, it can be challenging to determine who benefits from them. Determining the causes and consequences of infection-induced life-history changes requires the ability to experimentally manipulate life history and a framework for connecting life history to host and parasite fitness. Here, we combine a novel starvation manipulation with energy budget models to provide new insights into castration and gigantism in the Daphnia magna–Pasteuria ramosa host–parasite system. Our results show that starvation primarily affects investment in reproduction, and increasing starvation stress reduces gigantism and parasite fitness without affecting castration. These results are consistent with an energetic structure where the parasite uses growth energy as a resource. This finding gives us new understanding of the role of castration and gigantism in this system, and how life-history variation will affect infection outcome and epidemiological dynamics. The approach of combining targeted life-history manipulations with energy budget models can be adapted to understand life-history changes in other disease systems.     

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.