Rapidly fluctuating environments constrain coevolutionary arms races by impeding selective sweeps

Although pervasive, the impact of temporal environmental heterogeneity on coevolutionary processes is poorly understood. Productivity is a key temporally heterogeneous variable, and increasing productivity has been shown to increase rates of antagonistic arms race coevolution, and lead to the evolution of more broadly resistant hosts and more broadly infectious parasites. We investigated the effects of the grain of environmental heterogeneity, in terms of fluctuations in productivity, on bacteria–phage coevolution. Our findings demonstrate that environmental heterogeneity could constrain antagonistic coevolution, but that its effect was dependent upon the grain of heterogeneity, such that both the rate and extent of coevolution were most strongly limited in fine-grained, rapidly fluctuating heterogeneous environments. We further demonstrate that rapid environmental fluctuations were likely to have impeded selective sweeps of resistance alleles, which occurred over longer durations than the fastest, but not the slowest, frequency of fluctuations used. Taken together our results suggest that fine-grained environmental heterogeneity constrained the coevolutionary arms race by impeding selective sweeps.     

IMPACT OF BACTERIAL MUTATION RATE ON COEVOLUTIONARY DYNAMICS BETWEEN BACTERIA AND PHAGES

Mutator bacteria are frequently found in natural populations of bacteria and although coevolution with parasitic viruses (phages) is thought to be one reason for their persistence, it remains unclear how the presence of mutators affects coevolutionary dynamics. We hypothesized that phages must themselves adapt more rapidly or go extinct, in the face of rapidly evolving mutator bacteria. We compared the coevolutionary dynamics of wild‐type Pseudomonas fluorescens SBW25 with a lytic phage to the dynamics of an isogenic mutator of P. fluorescens SBW25 together with the same phage. At the beginning of the experiment both wild‐type bacteria and mutator bacteria coevolved with phages. However, mutators rapidly evolved higher levels of sympatric resistance to phages. The phages were unable to “keep‐up” with the mutator bacteria, and these rates of coevolution declined to less than the rates of coevolution between the phages and wild‐type bacteria. By the end of the experiment, the sympatric resistance of the mutator bacteria was not significantly different to the sympatric resistance of the wild‐type bacteria. This suggests that the importance of mutators in the coevolutionary interactions with a particular phage population is likely to be short‐lived. More generally, the results demonstrate that coevolving enemies may escape from Red‐Queen dynamics.     

Hot Spots, Cold Spots, and the Geographic Mosaic Theory of Coevolution

Species interactions commonly coevolve as complex geographic mosaics of populations shaped by differences in local selection and gene flow. We use a haploid matching-alleles model for coevolution to evaluate how a pair of species coevolves when fitness interactions are reciprocal in some locations ("hot spots") but not in others ("cold spots"). Our analyses consider mutualistic and antagonistic interspecific interactions and a variety of gene flow patterns between hot and cold spots. We found that hot and cold spots together with gene flow influence coevolutionary dynamics in four important ways. First, hot spots need not be ubiquitous to have a global influence on evolution, although rare hot spots will not have a disproportionate impact unless selection is relatively strong there. Second, asymmetries in gene flow can influence local adaptation, sometimes creating stable equilibria at which species experience minimal fitness in hot spots and maximal fitness in cold spots, or vice versa. Third, asymmetries in gene flow are no more important than asymmetries in population regulation for determining the maintenance of local polymorphisms through coevolution. Fourth, intraspecific allele frequency differences among hot and cold spot populations evolve under some, but not all, conditions. That is, selection mosaics are indeed capable of producing spatially variable coevolutionary outcomes across the landscapes over which species interact. Altogether, our analyses indicate that coevolutionary trajectories can be strongly shaped by the geographic distribution of coevolutionary hot and cold spots, and by the pattern of gene flow among populations.     

Nutrient stoichiometry shapes microbial coevolution

Coevolution is a force contributing to the generation and maintenance of biodiversity. It is influenced by environmental conditions including the scarcity of essential resources, which can drive the evolution of defence and virulence traits. We conducted a long‐term chemostat experiment where the marine cyanobacterium Synechococcus was challenged with a lytic phage under nitrogen (N) or phosphorus (P) limitation. This manipulation of nutrient stoichiometry altered the stability of host–parasite interactions and the underlying mode of coevolution. By assessing the infectivity with > 18 000 pairwise challenges, we documented directional selection for increased phage resistance, consistent with arms‐race dynamics while phage infectivity fluctuated through time, as expected when coevolution is driven by negative frequency‐dependent selection. The resulting infection networks were 50% less modular under N‐ versus P‐limitation reflecting host‐range contraction and asymmetric coevolutionary trajectories. Nutrient stoichiometry affects eco‐evolutionary feedbacks in ways that may alter the dynamics and functioning of environmental and host‐associated microbial communities.     

Antagonistic coevolution across productivity gradients: an experimental test of the effects of dispersal

Coevolution commonly occurs in spatially heterogeneous environments, resulting in variable selection pressures acting on coevolving species. Dispersal across such environments is predicted to have a major impact on local coevolutionary dynamics. Here, we address how co‐dispersal of coevolving populations of host and parasite across an environmental productivity gradient affected coevolution in experimental populations of bacteria and their parasitic viruses (phages). The rate of coevolution between bacteria and phages was greater in high‐productivity environments. High‐productivity immigrants (～2% of the recipient population) caused coevolutionary dynamics (rates of coevolution and degree of generalist evolution) in low‐productivity environments to be largely indistinguishable from high‐productivity environments, whereas immigration from low‐productivity environments (～0.5% of the population) had no discernable impact. These results could not be explained by demography alone, but rather high‐productivity immigrants had a selective advantage in low‐productivity environments, but not vice versa. Coevolutionary interactions in high‐productivity environments are therefore likely to have a disproportionate impact on coevolution across the landscape as a whole.     

Plasmid carriage can limit bacteria–phage coevolution

Coevolution with bacteriophages is a major selective force shaping bacterial populations and communities. A variety of both environmental and genetic factors has been shown to influence the mode and tempo of bacteria–phage coevolution. Here, we test the effects that carriage of a large conjugative plasmid, pQBR103, had on antagonistic coevolution between the bacterium Pseudomonas fluorescens and its phage, SBW25ϕ2. Plasmid carriage limited bacteria–phage coevolution; bacteria evolved lower phage-resistance and phages evolved lower infectivity in plasmid-carrying compared with plasmid-free populations. These differences were not explained by effects of plasmid carriage on the costs of phage resistance mutations. Surprisingly, in the presence of phages, plasmid carriage resulted in the evolution of high frequencies of mucoid bacterial colonies. Mucoidy can provide weak partial resistance against SBW25ϕ2, which may have limited selection for qualitative resistance mutations in our experiments. Taken together, our results suggest that plasmids can have evolutionary consequences for bacteria that go beyond the direct phenotypic effects of their accessory gene cargo.     

Genetic basis of infectivity evolution in a bacteriophage

Antagonistic coevolution between hosts and parasites is probably ubiquitous. However, very little is known of the genetic changes associated with parasite infectivity evolution during adaptation to a coevolving host. We followed the phenotypic and genetic changes in a lytic virus population (bacteriophage; phage Φ2) that coevolved with its bacterial host, Pseudomonas fluorescens SBW25. First, we show the rapid evolution of numerous unique phage infectivity phenotypes, and that both phage host range and bacterial resistance to individual phage increased over coevolutionary time. Second, each of the distinct phage phenotypes in our study had a unique genotype, and molecular evolution did not act uniformly across the phage genome during coevolution. In particular, we detected numerous substitutions on the tail fibre gene, which is involved in the first step of the host-parasite interaction: host adsorption. None of the observed mutations could be directly linked with infection against a particular host, suggesting that the phenotypic effects of infectivity mutations are probably epistatic. However, phage genotypes with the broadest host ranges had the largest number of nonsynonymous amino acid changes on genes implicated in infectivity evolution. An understanding of the molecular genetics of phage infectivity has helped to explain the complex phenotypic coevolutionary dynamics in this system.     

Host-parasite coevolutionary arms races give way to fluctuating selection

Host-parasite coevolution is a key driver of biological diversity and parasite virulence, but its effects depend on the nature of coevolutionary dynamics over time. We used phenotypic data from coevolving populations of the bacterium Pseudomonas fluorescens SBW25 and parasitic phage SBW25Φ2, and genetic data from the phage tail fibre gene (implicated in infectivity evolution) to show that arms race dynamics, typical of short-term studies, decelerate over time. We attribute this effect to increasing costs of generalism for phages and bacteria with increasing infectivity and resistance. By contrast, fluctuating selection on individual host and parasite genotypes was maintained over time, becoming increasingly important for the phenotypic properties of parasite and host populations. Given that costs of generalism are reported for many other systems, arms races may generally give way to fluctuating selection in antagonistically coevolving populations.     

The impact of parasite dispersal on antagonistic host-parasite coevolution

Coevolving populations of hosts and parasites are often subdivided into a set of patches connected by dispersal. Higher relative rates of parasite compared with host dispersal are expected to lead to parasite local adaptation. However, we know of no studies that have considered the implications of higher relative rates of parasite dispersal for other aspects of the coevolutionary process, such as the rate of coevolution and extent of evolutionary escalation of resistance and infectivity traits. We investigated the effect of phage dispersal on coevolution in experimental metapopulations of the bacterium Pseudomonas fluorescens SBW25 and its viral parasite, phage SBW25Phi2. Both the rate of coevolution and the breadth of evolved infectivity and resistance ranges peaked at intermediate rates of parasite dispersal. These results suggest that parasite dispersal can enhance the evolutionary potential of parasites through provision of novel genetic variation, but that high rates of parasite dispersal can impede the evolution of parasites by homogenizing genetic variation between patches, thereby constraining coevolution.     

Antagonistic coevolution limits population persistence of a virus in a thermally deteriorating environment

Understanding the conditions under which rapid evolutionary adaptation can prevent population extinction in deteriorating environments (i.e. evolutionary rescue) is a crucial aim in the face of global climate change. Despite a rapidly growing body of work in this area, little attention has been paid to the importance of interspecific coevolutionary interactions. Antagonistic coevolution commonly observed between hosts and parasites is likely to retard evolutionary rescue because it often reduces population sizes, and results in the evolution of costly host defence and parasite counter-defence. We used experimental populations of a bacterium Pseudomonas fluorescens SBW25 and a bacteriophage virus (SBW25Φ2), to study how host-parasite coevolution impacts viral population persistence in the face of gradually increasing temperature, an environmental stress for the virus but not the bacterium. The virus persisted much longer when it evolved in the presence of an evolutionarily constant host genotype (i.e. in the absence of coevolution) than when the bacterium and virus coevolved. Further experiments suggest that both a reduction in population size and costly infectivity strategies contributed to viral extinction as a result of coevolution. The results highlight the importance of interspecific evolutionary interactions for the evolutionary responses of populations to global climate change.     

Pulsed-resource dynamics increase the asymmetry of antagonistic coevolution between a predatory protist and a prey bacterium

Temporal resource fluctuations could affect the strength of antagonistic coevolution through population dynamics and costs of adaptation. We studied this by coevolving the prey bacterium Serratia marcescens with the predatory protozoa Tetrahymena thermophila in constant and pulsed-resource environments for approximately 1300 prey generations. Consistent with arms race theory, the prey evolved to be more defended, whereas the predator evolved to be more efficient in consuming the bacteria. Coevolutionary adaptations were costly in terms of reduced prey growth in resource-limited conditions and less efficient predator growth on nonliving resource medium. However, no differences in mean coevolutionary changes or adaptive costs were observed between environments, even though resource pulses increased fluctuations and mean densities of coevolving predator populations. Interestingly, a surface-associated prey defence mechanism (bacterial biofilm), to which predators were probably unable to counter-adapt, evolved to be stronger in pulsed-resource environment. These results suggest that temporal resource fluctuations can increase the asymmetry of antagonistic coevolution by imposing stronger selection on one of the interacting species.     

Phenotypic plasticity,but not adaptive tracking,underlies seasonal variation in post‐cold hardening freeze tolerance of Drosophila melanogaster

In temperate regions, an organism's ability to rapidly adapt to seasonally varying environments is essential for its survival. In response to seasonal changes in selection pressure caused by variation in temperature, humidity, and food availability, some organisms exhibit plastic changes in phenotype. In other cases, seasonal variation in selection pressure can rapidly increase the frequency of genotypes that offer survival or reproductive advantages under the current conditions. Little is known about the relative influences of plastic and genetic changes in short‐lived organisms experiencing seasonal environmental fluctuations. Cold hardening is a seasonally relevant plastic response in which exposure to cool, but nonlethal, temperatures significantly increases the organism's ability to later survive at freezing temperatures. In the present study, we demonstrate seasonal variation in cold hardening in Drosophila melanogaster and test the extent to which plasticity and adaptive tracking underlie that seasonal variation. We measured the post‐cold hardening freeze tolerance of flies from outdoor mesocosms over the summer, fall, and winter. We bred outdoor mesocosm‐caught flies for two generations in the laboratory and matched each outdoor cohort to an indoor control cohort of similar genetic background. We cold hardened all flies under controlled laboratory conditions and then measured their post‐cold hardening freeze tolerance. Comparing indoor and field‐caught flies and their laboratory‐reared G1 and G2 progeny allowed us to determine the roles of seasonal environmental plasticity, parental effects, and genetic changes on cold hardening. We also tested the relationship between cold hardening and other factors, including age, developmental density, food substrate, presence of antimicrobials, and supplementation with live yeast. We found strong plastic responses to a variety of field‐ and laboratory‐based environmental effects, but no evidence of seasonally varying parental or genetic effects on cold hardening. We therefore conclude that seasonal variation in post‐cold hardening freeze tolerance results from environmental influences and not genetic changes.     

Effects of predation on real‐time host–parasite coevolutionary dynamics

The impact of community complexity on pairwise coevolutionary dynamics is theoretically dependent on the extent to which species evolve generalised or specialised adaptations to the multiple species they interact with. Here, we show that the bacteria Pseudomonas fluorescens diversifies into defence specialists, when coevolved simultaneously with a virus and a predatory protist, as a result of fitness trade‐offs between defences against the two enemies. Strong bacteria–virus pairwise coevolution persisted, despite strong protist‐imposed selection. However, the arms race dynamic (escalation of host resistance and parasite infectivity ranges) associated with bacteria–virus coevolution broke down to a greater extent in the presence of the protist, presumably through the elevated genetic and demographic costs of increased bacteria resistance ranges. These findings suggest that strong pairwise coevolution can persist even in complex communities, when conflicting selection leads to evolutionary diversification of different defence strategies.     

Host population bottlenecks drive parasite extinction during antagonistic coevolution

Host–parasite interactions are often characterized by large fluctuations in host population size, and we investigated how such host bottlenecks affected coevolution between a bacterium and a virus. Previous theory suggests that host bottlenecks should provide parasites with an evolutionary advantage, but instead we found that phages were rapidly driven to extinction when coevolving with hosts exposed to large genetic bottlenecks. This was caused by the stochastic loss of sensitive bacteria, which are required for phage persistence and infectivity evolution. Our findings emphasize the importance of feedbacks between ecological and coevolutionary dynamics, and how this feedback can qualitatively alter coevolutionary dynamics.     

Host–parasite coevolution favours parasite genetic diversity and horizontal gene transfer

Host–parasite coevolution is predicted to favour genetic diversity and the underlying mechanisms (e.g. sexual reproduction and, more generally, genetic exchange), because diversity enhances the antagonists' potential for rapid adaptation. To date, this prediction has mainly been tested and confirmed for the host. It should similarly apply to the parasite. Indeed, our previous work demonstrated that experimental coevolution between the nematode Caenorhabditis elegans and its microparasite Bacillus thuringiensis selects for genetic diversity in both antagonists. For the parasite, the previous analysis was based on plasmid‐encoded toxin gene markers. Thus, it was restricted to a very small part of the bacterial genome and did not cover the main chromosome, which harbours a large variety of virulence factors. Here, we present new data for chromosomal gene markers of B. thuringiensis and combine this information with the previous results on plasmid‐encoded toxins. Our new results demonstrate that, in comparison with the control treatment, coevolution with a host similarly leads to higher levels of genetic diversity in the bacterial chromosome, thus indicating the relevance of chromosomal genes for coevolution. Furthermore, the frequency of toxin gene gain is significantly elevated during coevolution, highlighting the importance of horizontal gene transfer as a diversity‐generating mechanism. In conclusion, our study emphasizes the strong influence of antagonistic coevolution on parasite genetic diversity and gene exchange.     

Spatial population genetic structure of a bacterial parasite in close coevolution with its host

Knowledge of a species’ population genetic structure can provide insight into fundamental ecological and evolutionary processes including gene flow, genetic drift and adaptive evolution. Such inference is of particular importance for parasites, as an understanding of their population structure can illuminate epidemiological and coevolutionary dynamics. Here, we describe the population genetic structure of the bacterium Pasteuria ramosa, a parasite that infects planktonic crustaceans of the genus Daphnia. This system has become a model for investigations of host–parasite interactions and represents an example of coevolution via negative frequency‐dependent selection (aka “Red Queen” dynamics). To sample P. ramosa, we experimentally infected a panel of Daphnia hosts with natural spore banks from the sediments of 25 ponds throughout much of the species range in Europe and western Asia. Using 12 polymorphic variable number tandem repeat loci (VNTR loci), we identified substantial genetic diversity, both within and among localities, that was structured geographically among ponds. Genetic diversity was also structured among host genotypes within ponds, although this pattern varied by locality, with P. ramosa at some localities partitioned into distinct host‐specific lineages, and other localities where recombination had shuffled genetic variation among different infection phenotypes. Across the sample range, there was a pattern of isolation by distance, and principal components analysis coupled with Procrustes rotation identified congruence between patterns of genetic variation and geography. Our findings support the hypothesis that Pasteuria is an endemic parasite coevolving closely with its host. These results provide important context for previous studies of this model system and inform hypotheses for future research.     

Experimental evolution for generalists and specialists reveals multivariate genetic constraints on thermal reaction norms

Theory predicts the emergence of generalists in variable environments and antagonistic pleiotropy to favour specialists in constant environments, but empirical data seldom support such generalist–specialist trade‐offs. We selected for generalists and specialists in the dung fly Sepsis punctum (Diptera: Sepsidae) under conditions that we predicted would reveal antagonistic pleiotropy and multivariate trade‐offs underlying thermal reaction norms for juvenile development. We performed replicated laboratory evolution using four treatments: adaptation at a hot (31 °C) or a cold (15 °C) temperature, or under regimes fluctuating between these temperatures, either within or between generations. After 20 generations, we assessed parental effects and genetic responses of thermal reaction norms for three correlated life‐history traits: size at maturity, juvenile growth rate and juvenile survival. We find evidence for antagonistic pleiotropy for performance at hot and cold temperatures, and a temperature‐mediated trade‐off between juvenile survival and size at maturity, suggesting that trade‐offs associated with environmental tolerance can arise via intensified evolutionary compromises between genetically correlated traits. However, despite this antagonistic pleiotropy, we found no support for the evolution of increased thermal tolerance breadth at the expense of reduced maximal performance, suggesting low genetic variance in the generalist–specialist dimension.     

Mutualism variation in the nodulation response to nitrate

The evolution of mutualisms under novel selective pressures will play a key role in ecosystem responses to environmental change. Because fixed nitrogen is traded in plant–rhizobium mutualisms, increasing N availability in the soil is predicted to alter coevolution of these interactions. Legumes typically decrease the number of associations (nodules) with rhizobia in response to nitrate, but the evolutionary dynamics of this response remain unknown. We grew plant and rhizobium genotype combinations in three N environments to assess the coevolutionary potential of the nodule nitrate response in natural communities of plants and rhizobia. We found evidence for coevolutionary genetic variation for nodulation in response to nitrate (G × G × E interaction), suggesting that the mutualism response to N deposition will depend on the combination of partner genotypes. Thus, the nitrate response is not a fixed mechanism in plant–rhizobium symbioses, but instead is potentially subject to natural selection and dynamic coevolution.     

INTERLOCUS SEXUALLY ANTAGONISTIC COEVOLUTION CAN CREATE INDIRECT SELECTION FOR INCREASED RECOMBINATION

The ubiquity of recombination in nature is a paradox because it breaks up combinations of alleles favored by natural selection. Theoretical work has shown that antagonistic coevolution between hosts and parasites can result in rapid fluctuations in epistasis that can create a short‐term advantage to recombination. Here, we show that another kind of antagonistic coevolution, interlocus sexually antagonistic coevolution (SAC), can also create indirect selection for modifiers that increase the rate of recombination, and that it can lead to very high levels of recombination at equilibrium. Recombination is favored because interlocus SAC creates heterogeneity in the strength and direction of selection, both within and between generations, which maintains an excess of disadvantageous haplotypes in the population. This result is similar to and consistent with dynamics of fluctuating epistasis produced in models of host–parasite coevolution. However, the conditions under which interlocus SAC provides an advantage to recombination are more permissive.     

No effect of host–parasite co‐evolution on host range expansion

Antagonistic co‐evolution between hosts and parasites (reciprocal selection for resistance and infectivity) is hypothesized to play an important role in host range expansion by selecting for novel infectivity alleles, but tests are lacking. Here, we determine whether experimental co‐evolution between a bacterium (Pseudomonas fluorescens SBW25) and a phage (SBW25Φ2) affects interstrain host range: the ability to infect different strains of P. fluorescens other than SBW25. We identified and tested a genetically and phenotypically diverse suite of co‐evolved phage variants of SBW25Φ2 against both sympatric and allopatric co‐evolving hosts (P. fluorescens SBW25) and a large set of other P. fluorescens strains. Although all co‐evolved phage had a greater host range than the ancestral phage and could differentially infect co‐evolved variants of P. fluorescens SBW25, none could infect any of the alternative P. fluorescens strains. Thus, parasite generalism at one genetic scale does not appear to affect generalism at other scales, suggesting fundamental genetic constraints on parasite adaptation for this virus.