EVOLUTIONARY RESPONSE OF PREDATORS TO DANGEROUS PREY: PREADAPTATION AND THE EVOLUTION OF TETRODOTOXIN RESISTANCE IN GARTER SNAKES

Conspicuous coloration of toxin-resistant predators implicates additional trophic interactions in a predator-prey arms race

Antagonistic coevolution between natural enemies can produce highly exaggerated traits, such as prey toxins and predator resistance. This reciprocal process of adaptation and counter-adaptation may also open doors to other evolutionary novelties not directly involved in the phenotypic interface of coevolution. We tested the hypothesis that predator-prey coevolution coincided with the evolution of conspicuous coloration on resistant predators that retain prey toxins. In western North America, common garter snakes (Thamnophis sirtalis) have evolved extreme resistance to tetrodotoxin (TTX) in the coevolutionary arms race with their deadly prey, Pacific newts (Taricha spp.). TTX-resistant snakes can retain large amounts of ingested TTX, which could serve as a deterrent against the snakes' own predators if TTX toxicity and resistance are coupled with a conspicuous warning signal. We evaluated whether arms race escalation covaries with bright red coloration in snake populations across the geographic mosaic of coevolution. Snake colour variation departs from the neutral expectations of population genetic structure and covaries with escalating clines of newt TTX and snake resistance at two coevolutionary hotspots. In the Pacific Northwest, bright red coloration fits an expected pattern of an aposematic warning to avian predators: TTX-resistant snakes that consume highly toxic newts also have relatively large, reddish-orange dorsal blotches. Snake coloration also seems to have evolved with the arms race in California, but overall patterns are less intuitively consistent with aposematism. These results suggest that interactions with additional trophic levels can generate novel traits as a cascading consequence of arms race coevolution across the geographic mosaic.

Evolution in interacting species alters predator life-history traits, behaviour and morphology in experimental microbial communities

Predator-prey interactions heavily influence the dynamics of many ecosystems. An increasing body of evidence suggests that rapid evolution and coevolution can alter these interactions, with important ecological implications, by acting on traits determining fitness, including reproduction, anti-predatory defence and foraging efficiency. However, most studies to date have focused only on evolution in the prey species, and the predator traits in (co)evolving systems remain poorly understood. Here, we investigated changes in predator traits after approximately 600 generations in a predator-prey (ciliate-bacteria) evolutionary experiment. Predators independently evolved on seven different prey species, allowing generalization of the predator's evolutionary response. We used highly resolved automated image analysis to quantify changes in predator life history, morphology and behaviour. Consistent with previous studies, we found that prey evolution impaired growth of the predator, although the effect depended on the prey species. By contrast, predator evolution did not cause a clear increase in predator growth when feeding on ancestral prey. However, predator evolution affected morphology and behaviour, increasing size, speed and directionality of movement, which have all been linked to higher prey search efficiency. These results show that in (co)evolving systems, predator adaptation can occur in traits relevant to foraging efficiency without translating into an increased ability of the predator to grow on the ancestral prey type.

Pick your trade-offs wisely: Predator-prey eco-evo dynamics are qualitatively different under different trade-offs

In recent decades, myriad studies have explored the population dynamics of coevolving populations of predator and prey. A variety of choices are made in these models: exponential or logistic prey growth in the absence of a predator, various forms of predator functional response, and uni- or bi-directional trait axes. In addition, some form of trade-offs are assumed in order to prevent run-away evolution of the prey and predator traits. While there is a considerable amount of theory regarding various forms of prey growth rates and predator functional responses, only a few studies have explored how different types of trade-offs affect predator-prey dynamics. Here, we compared two ditrophic coevolution models incorporating different trade-offs via dual effects of the prey trait on attack rate and either prey carrying capacity or intrinsic growth rate. We employed a standard dynamical systems approach to analyze the equilibrium conditions of each model and find conditions for non-equilibrium oscillatory coexistence. The exact effect of various parameters on the outcome of predator-prey interactions depends on whether the trade-offs affect the intrinsic growth rate or carrying capacity. In particular, coexistence is more likely when prey growth rate is affected by the evolving trait. In addition, in parameter regimes where cycles occur in both models, oscillations typically have larger periods and amplitudes when prey growth rate is affected by the evolving trait.

Toxicity and population structure of the Rough-Skinned Newt (Taricha granulosa) outside the range of an arms race with resistant predators

Species interactions, and their fitness consequences, vary across the geographic range of a coevolutionary relationship. This spatial heterogeneity in reciprocal selection is predicted to generate a geographic mosaic of local adaptation, wherein coevolutionary traits are phenotypically variable from one location to the next. Under this framework, allopatric populations should lack variation in coevolutionary traits due to the absence of reciprocal selection. We examine phenotypic variation in tetrodotoxin (TTX) toxicity of the Rough-Skinned Newt (Taricha granulosa) in regions of allopatry with its TTX-resistant predator, the Common Garter Snake (Thamnophis sirtalis). In sympatry, geographic patterns of phenotypic exaggeration in toxicity and toxin-resistance are closely correlated in prey and predator, implying that reciprocal selection drives phenotypic variation in coevolutionary traits. Therefore, in allopatry with TTX-resistant predators, we expect to find uniformly low levels of newt toxicity. We characterized TTX toxicity in northwestern North America, including the Alaskan panhandle where Ta. granulosa occur in allopatry with Th. sirtalis. First, we used microsatellite markers to estimate population genetic structure and determine if any phenotypic variation in toxicity might be explained by historical divergence. We found northern populations of Ta. granulosa generally lacked population structure in a pattern consistent with northern range expansion after the Pleistocene. Next, we chose a cluster of sites in Alaska, which uniformly lacked genetic divergence, to test for phenotypic divergence in toxicity. As predicted, overall levels of newt toxicity were low; however, we also detected unexpected among- and within-population variation in toxicity. Most notably, a small number of individuals contained large doses of TTX that rival means of toxic populations in sympatry with Th. sirtalis. Phenotypic variation in toxicity, despite limited neutral genetic divergence, suggests that factors other than reciprocal selection with Th. sirtalis likely contribute to geographic patterns of toxicity in Ta. granulosa.

Constraint shapes convergence in tetrodotoxin-resistant sodium channels of snakes

Natural selection often produces convergent changes in unrelated lineages, but the degree to which such adaptations occur via predictable genetic paths is unknown. If only a limited subset of possible mutations is fixed in independent lineages, then it is clear that constraint in the production or function of molecular variants is an important determinant of adaptation. We demonstrate remarkably constrained convergence during the evolution of resistance to the lethal poison, tetrodotoxin, in six snake species representing three distinct lineages from around the globe. Resistance-conferring amino acid substitutions in a voltage-gated sodium channel, Na(v)1.4, are clustered in only two regions of the protein, and a majority of the replacements are confined to the same three positions. The observed changes represent only a small fraction of the experimentally validated mutations known to increase Na(v)1.4 resistance to tetrodotoxin. These results suggest that constraints resulting from functional tradeoffs between ion channel function and toxin resistance led to predictable patterns of evolutionary convergence at the molecular level. Our data are consistent with theoretical predictions and recent microcosm work that suggest a predictable path is followed during an adaptive walk along a mutational landscape, and that natural selection may be frequently constrained to produce similar genetic outcomes even when operating on independent lineages.

Genetic architecture of a feeding adaptation: garter snake (Thamnophis) resistance to tetrodotoxin bearing prey

Detailing the genetic basis of adaptive variation in natural populations is a first step towards understanding the process of adaptive evolution, yet few ecologically relevant traits have been characterized at the genetic level in wild populations. Traits that mediate coevolutionary interactions between species are ideal for studying adaptation because of the intensity of selection and the well-characterized ecological context. We have previously described the ecological context, evolutionary history and partial genetic basis of tetrodotoxin (TTX) resistance in garter snakes (Thamnophis). Derived mutations in a voltage-gated sodium channel gene (Na(v)1.4) in three garter snake species are associated with resistance to TTX, the lethal neurotoxin found in their newt prey (Taricha). Here we evaluate the contribution of Na(v)1.4 alleles to TTX resistance in two of those species from central coastal California. We measured the phenotypes (TTX resistance) and genotypes (Na(v)1.4 and microsatellites) in a local sample of Thamnophis atratus and Thamnophis sirtalis. Allelic variation in Na(v)1.4 explains 23 per cent of the variation in TTX resistance in T. atratus while variation in a haphazard sample of the genome (neutral microsatellite markers) shows no association with the phenotype. Similarly, allelic variation in Na(v)1.4 correlates almost perfectly with TTX resistance in T. sirtalis, but neutral variation does not. These strong correlations suggest that Na(v)1.4 is a major effect locus. The simple genetic architecture of TTX resistance in garter snakes may significantly impact the dynamics of phenotypic coevolution. Fixation of a few alleles of major effect in some garter snake populations may have led to the evolution of extreme phenotypes and an 'escape' from the arms race with newts.

The evolutionary origins of beneficial alleles during the repeated adaptation of garter snakes to deadly prey

Where do the genetic variants underlying adaptive change come from? Are currently adaptive alleles recruited by selection from standing genetic variation within populations, moved through introgression from other populations, or do they arise as novel mutations? Here, we examine the molecular basis of repeated adaptation to the toxin of deadly prey in 3 species of garter snakes (Thamnophis) to determine whether adaptation has evolved through novel mutations, sieving of existing variation, or transmission of beneficial alleles across species. Functional amino acid substitutions in the skeletal muscle sodium channel (Na(v)1.4) are largely responsible for the physiological resistance of garter snakes to tetrodotoxin found in their newt (Taricha) prey. Phylogenetic analyses reject the hypotheses that the unique resistance alleles observed in multiple Thamnophis species were present before the split of these lineages, or that alleles were shared among species through occasional hybridization events. Our results demonstrate that adaptive evolution has occurred independently multiple times in garter snakes via the de novo acquisition of beneficial mutations.

Phenotypic mismatches reveal escape from arms-race coevolution

Because coevolution takes place across a broad scale of time and space, it is virtually impossible to understand its dynamics and trajectories by studying a single pair of interacting populations at one time. Comparing populations across a range of an interaction, especially for long-lived species, can provide insight into these features of coevolution by sampling across a diverse set of conditions and histories. We used measures of prey traits (tetrodotoxin toxicity in newts) and predator traits (tetrodotoxin resistance of snakes) to assess the degree of phenotypic mismatch across the range of their coevolutionary interaction. Geographic patterns of phenotypic exaggeration were similar in prey and predators, with most phenotypically elevated localities occurring along the central Oregon coast and central California. Contrary to expectations, however, these areas of elevated traits did not coincide with the most intense coevolutionary selection. Measures of functional trait mismatch revealed that over one-third of sampled localities were so mismatched that reciprocal selection could not occur given current trait distributions. Estimates of current locality-specific interaction selection gradients confirmed this interpretation. In every case of mismatch, predators were “ahead” of prey in the arms race; the converse escape of prey was never observed. The emergent pattern suggests a dynamic in which interacting species experience reciprocal selection that drives arms-race escalation of both prey and predator phenotypes at a subset of localities across the interaction. This coadaptation proceeds until the evolution of extreme phenotypes by predators, through genes of large effect, allows snakes to, at least temporarily, escape the arms race.

Arms races between natural enemies can lead to the rapid evolution of extreme traits, high degrees of specialization, and the formation of new species. They also serve as the ecological model for the evolution of drug resistance by diseases and for host–pathogen interactions in general. Revealing who wins these arms races and how they do so is critical to our understanding of these processes. Capitalizing on the geographic mosaic of species interactions, we examined the dynamics of the arms race between snakes and their toxic newt prey. Garter snakes in some populations have evolved dramatic resistance to the tetrodotoxin defense of the their local prey. By evaluating the pattern of mismatches between toxicity and resistance, we discovered that predators sometimes escape the arms race through the evolution of extreme resistance, but that prey never come out ahead. The reason for this one-sided outcome appears to depend on the molecular genetic basis of resistance in snakes, wherein changes to a single amino acid residue can confer huge differences in resistance.

Who wins in the arms race between predators and prey? In the interaction between snakes and toxic newts, predators sometimes escape the arms race through the evolution of extreme resistance, but prey never come out ahead.

Parallel arms races between garter snakes and newts involving tetrodotoxin as the phenotypic interface of coevolution

Parallel "arms races" involving the same or similar phenotypic interfaces allow inference about selective forces driving coevolution, as well as the importance of phylogenetic and phenotypic constraints in coevolution. Here, we report the existence of apparent parallel arms races between species pairs of garter snakes and their toxic newt prey that indicate independent evolutionary origins of a key phenotype in the interface. In at least one area of sympatry, the aquatic garter snake, Thamnophis couchii, has evolved elevated resistance to the neurotoxin tetrodotoxin (TTX), present in the newt Taricha torosa. Previous studies have shown that a distantly related garter snake, Thamnophis sirtalis, has coevolved with another newt species that possesses TTX, Taricha granulosa. Patterns of within population variation and phenotypic tradeoffs between TTX resistance and sprint speed suggest that the mechanism of resistance is similar in both species of snake, yet phylogenetic evidence indicates the independent origins of elevated resistance to TTX.

Predator Mixes and the Conspicuousness of Aposematic Signals

Conspicuous warning signals of unprofitable prey are a defense against visually hunting predators. They work because predators learn to associate unprofitability with bright coloration and because strong signals are detectable and memorable. However, many species that can be considered defended are not very conspicuous; they have weak warning signals. This phenomenon has previously been ignored in models and experiments. In addition, there is significant within- and among-species variation among predators in their search behavior, in their visual, cognitive, and learning abilities, and in their resistance to defenses. In this article we explore the effects of variable predators on models that combine positive frequency-dependent, frequency-independent, and negative frequency-dependent predation and show that weak signaling of aposematic species can evolve if predators vary in their tendency to attack defended prey.

The evolutionary response of predators to dangerous prey: hotspots and coldspots in the geographic mosaic of coevolution between garter snakes and newts

The "geographic mosaic" approach to understanding coevolution is predicated on the existence of variable selection across the landscape of an interaction between species. A range of ecological factors, from differences in resource availability to differences in community composition, can generate such a mosaic of selection among populations, and thereby differences in the strength of coevolution. The result is a mixture of hotspots, where reciprocal selection is strong, and coldspots, where reciprocal selection is weak or absent, throughout the ranges of species. Population subdivision further provides the opportunity for nonadaptive forces, including gene flow, drift, and metapopulation dynamics, to influence the coevolutionary interaction between species. Some predicted results of this geographic mosaic of coevolution include maladapted or mismatched phenotypes, maintenance of high levels of polymorphism, and prevention of stable equilibrium trait combinations. To evaluate the potential for the geographic mosaic to influence predator-prey coevolution, we investigated the geographic pattern of genetically determined TTX resistance in the garter snake Thamnophis sirtalis over much of the range of its ecological interaction with toxic newts of genus Taricha. We assayed TTX resistance in over 2900 garter snakes representing 333 families from 40 populations throughout western North America. Our results provide dramatic evidence that geographic structure is an important component in coevolutionary interactions between predators and prey. Resistance levels vary substantially (over three orders of magnitude) among populations and over short distances. The spatial array of variation is consistent with two areas of intense evolutionary response by predators ("hotspots") surrounded by clines of decreasing resistance. Some general predictions of the geographic mosaic process are supported, including clinal variation in phenotypes, polymorphism in some populations, and divergent outcomes of the interaction between predator and prey. Conversely, our data provide little support for one of the major predictions, mismatched values of interacting traits. Two lines of evidence suggest selection is paramount in determining population variation in resistance. First, phylogenetic information indicates that two hotspots of TTX resistance have evolved independently. Second, in the one region that TTX levels in prey have been quantified, resistance and toxicity levels match almost perfectly over a wide phenotypic and geographic range. However, these results do not preclude the role the nonadaptive forces in generating the overall geographic mosaic of TTX resistance. Much work remains to fill in the geographic pattern of variation among prey populations and, just as importantly, to explore the variation in the ecology of the interaction that occurs within populations.