Convergent adaptation to dangerous prey proceeds through the same first‐step mutation in the garter snake Thamnophis sirtalis

Where Does All the Poison Go? Investigating Toxicokinetics of Newt (Taricha) Tetrodotoxin (TTX) in Garter Snakes (Thamnophis)

Animals that consume toxic diets provide models for understanding the molecular and physiological adaptations to ecological challenges. Garter snakes (Thamnophis) in western North America prey on Pacific newts (Taricha), which employ tetrodotoxin (TTX) as an antipredator defense. These snakes possess mutations in voltage-gated sodium channels (Nav), the molecular targets of TTX, that decrease the binding ability of TTX to sodium channels (target-site resistance). However, genetic variation at these loci that cannot explain all the phenotypic variation in TTX resistance in Thamnophis. We explored a separate means of resistance, toxin metabolism, to determine if TTX-resistant snakes either rapidly remove TTX or sequester TTX. We examined the metabolism and distribution of TTX in the body (toxicokinetics), to determine differences between TTX-resistant and TTX-sensitive snakes in the rates at which TTX is eliminated from organs and the whole body (using TTX half-life as our metric). We assayed TTX half-life in snakes from TTX-resistant and TTX-sensitive populations of three garter snake species with a coevolutionary history with newts (T. atratus, T. couchii, T. sirtalis), as well as two non-resistant "outgroup" species (T. elegans, Pituophis catenifer) that seldom (if ever) engage newts. We found TTX half-life varied across species, populations, and tissues. Interestingly, TTX half-life was shortest in T. elegans and P. catenifer compared to all other snakes. Furthermore, TTX-resistant populations of T. couchii and T. sirtalis eliminated TTX faster (shorter TTX half-life) than their TTX-sensitive counterparts, while populations of TTX-resistant and TTX-sensitive T. atratus showed no difference rates of TTX removal (same TTX half-life). The ability to rapidly eliminate TTX may have permitted increased prey consumption, which may have promoted the evolution of additional resistance mechanisms. Finally, snakes still retain substantial amounts of TTX, and we projected that snakes could be dangerous to their own predators days to weeks following the ingestion of a single newt. Thus, aspects of toxin metabolism may have been key in driving predator-prey relationships, and important in determining other ecological interactions.

Genetic architecture, spatial heterogeneity, and the coevolutionary arms race between newts and snakes

Coevolution between two species can lead to exaggerated phenotypes that vary in a correlated manner across space. However, the conditions under which we expect such spatially varying coevolutionary patterns in polygenic traits are not well-understood. We investigate the coevolutionary dynamics between two species undergoing reciprocal adaptation across space and time, using simulations inspired by the Taricha newt - Thamnophis garter snake system. One striking observation from this system is that newts in some areas carry much more tetrodotoxin than in other areas, and garter snakes that live near more toxic newts tend to be more resistant to this toxin, a correlation seen across several broad geographic areas. Furthermore, snakes seem to be "winning" the coevolutionary arms race, i.e., having a high level of resistance compared to local newt toxicity, despite substantial variation in both toxicity and resistance across the range. We explore how possible genetic architectures of the toxin and resistance traits would affect the coevolutionary dynamics by manipulating both mutation rate and effect size of mutations across many simulations. We find that coevolutionary dynamics alone were not sufficient in our simulations to produce the striking mosaic of levels of toxicity and resistance observed in nature, but simulations with ecological heterogeneity (in trait costliness or interaction rate) did produce such patterns. We also find that in simulations, newts tend to "win" across most combinations of genetic architectures, although the species with higher mutational genetic variance tends to have an advantage.

Coevolution with toxic prey produces functional trade-offs in sodium channels of predatory snakes

Seemingly unrelated traits often share the same underlying molecular mechanisms, potentially generating a pleiotropic relationship whereby selection shaping one trait can simultaneously compromise another. While such functional trade-offs are expected to influence evolutionary outcomes, their actual relevance in nature is masked by obscure links between genotype, phenotype, and fitness. Here, we describe functional trade-offs that likely govern a key adaptation and coevolutionary dynamics in a predator-prey system. Several garter snake (Thamnophis spp.) populations have evolved resistance to tetrodotoxin (TTX), a potent chemical defense in their prey, toxic newts (Taricha spp.). Snakes achieve TTX resistance through mutations occurring at toxin-binding sites in the pore of snake skeletal muscle voltage-gated sodium channels (NaV1.4). We hypothesized that these mutations impair basic NaV functions, producing molecular trade-offs that should ultimately scale up to compromised organismal performance. We investigate biophysical costs in two snake species with unique and independently evolved mutations that confer TTX resistance. We show electrophysiological evidence that skeletal muscle sodium channels encoded by toxin-resistant alleles are functionally compromised. Furthermore, skeletal muscles from snakes with resistance genotypes exhibit reduced mechanical performance. Lastly, modeling the molecular stability of these sodium channel variants partially explains the electrophysiological and muscle impairments. Ultimately, adaptive genetic changes favoring toxin resistance appear to negatively impact sodium channel function, skeletal muscle strength, and organismal performance. These functional trade-offs at the cellular and organ levels appear to underpin locomotor deficits observed in resistant snakes and may explain variation in the population-level success of toxin-resistant alleles across the landscape, ultimately shaping the trajectory of snake-newt coevolution.

The ripple effects of clines from coevolutionary hotspots to coldspots

Coevolution has the potential to alter not only the ecological interactions of coevolving partners, but also their interactions with yet other species. The effects of coevolution may ripple throughout networks of interacting species, cascading across trophic levels, swamping competitors, or facilitating survival or reproduction of yet other species linked only indirectly to the coevolving partners. These ripple effects of coevolution may differ among communities, amplifying how the coevolutionary process produces geographic mosaics of traits and outcomes in interactions among species. In a From the Cover article in this issue of Molecular Ecology, Hague et al. (2022) provide a clear example, using the well-studied interactions between Pacific newts (Taricha spp.) and their common garter snake (Thamnophis sirtalis) predators in western North America. Pacific newts harbour tetrodotoxin (TTX), which is highly toxic to vertebrate predators. In coevolutionary hotspots, extreme escalation of toxicity in the newts and resistance to toxicity in the snakes have resulted in snake populations that retain high levels of TTX. In two geographic regions, snakes in these hotspot populations have evolved bright, aposematic colours that may act as warning signals to their own vertebrate predators. The warning signals and toxin-resistance alleles in the snake populations decrease clinally away from the coevolutionary hotpots, shaped by a geographically variable mix of selection imposed by the snakes' prey and by their own predators.

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.

Tissue-specific plant toxins and adaptation in a specialist root herbivore

In coevolution between plants and insects, reciprocal selection often leads to phenotype matching between chemical defense and herbivore offense. Nonetheless, it is not well understood whether distinct plant parts are differentially defended and how herbivores adapted to those parts cope with tissue-specific defense. Milkweed plants produce a diversity of cardenolide toxins and specialist herbivores have substitutions in their target enzyme (Na+/K+-ATPase), each playing a central role in milkweed-insect coevolution. The four-eyed milkweed beetle (Tetraopes tetrophthalmus) is an abundant toxin-sequestering herbivore that feeds exclusively on milkweed roots as larvae and less so on milkweed leaves as adults. Accordingly, we tested the tolerance of this beetle's Na+/K+-ATPase to cardenolide extracts from roots versus leaves of its main host (Asclepias syriaca), along with sequestered cardenolides from beetle tissues. We additionally purified and tested the inhibitory activity of dominant cardenolides from roots (syrioside) and leaves (glycosylated aspecioside). Tetraopes' enzyme was threefold more tolerant of root extracts and syrioside than leaf cardenolides. Nonetheless, beetle-sequestered cardenolides were more potent than those in roots, suggesting selective uptake or dependence on compartmentalization of toxins away from the beetle's enzymatic target. Because Tetraopes has two functionally validated amino acid substitutions in its Na+/K+-ATPase compared to the ancestral form in other insects, we compared its cardenolide tolerance to that of wild-type Drosophila and CRISPR-edited Drosophila with Tetraopes' Na+/K+-ATPase genotype. Those two amino acid substitutions accounted for >50% of Tetraopes' enhanced enzymatic tolerance of cardenolides. Thus, milkweed's tissue-specific expression of root toxins is matched by physiological adaptations in its specialist root herbivore.

Identification of the NA+/K+-ATPase α-Isoforms in Six Species of Poison Dart Frogs and their Sensitivity to Cardiotonic Steroids

Cardiotonic steroids (CTS) are a group of compounds known to be toxic due to their ability to inhibit the Na⁺/K⁺-ATPase (NKA), which is essential to maintain the balance of ions in animal cells. An evolutionary strategy of molecular adaptation to avoid self-intoxication acquired by CTS defended organisms and their predators is the structural modification of their NKA where specific amino acid substitutions confer resistant phenotypes. Several lineages of poison dart frogs (Dendrobatidae) are well known to sequester a wide variety of lipophilic alkaloids from their arthropod diet, however there is no evidence of CTS-sequestration or dietary exposure. Interestingly this study identified the presence of α-NKA isoforms (α₁ and α₂) with amino acid substitutions indicative of CTS-resistant phenotypes in skeletal muscle transcriptomes obtained from six species of dendrobatids: Phyllobates aurotaenia, Oophaga anchicayensis, Epipedobates boulengeri, Andinobates bombetes, Andinobates minutus, and Leucostethus brachistriatus, collected in the Valle del Cauca (Colombia). P. aurotaenia, A. minutus, and E. boulengeri presented two variants for α₁-NKA, with one of them having these substitutions. In contrast, O. anchicayensis and A. bombetes have only one α₁-NKA isoform with an amino acid sequence indicative of CTS susceptibility and an α₂-NKA with one substitution that could confer a reduced affinity for CTS. The α₁ and α₂ isoforms of L. brachistriatus do not contain substitutions imparting CTS resistance. Our findings indicate that poison dart frogs express α-NKA isoforms with different affinities for CTS and the pattern of this expression might be influenced by factors related to evolutionary, physiological, ecological, and geographical burdens.

The road not taken: Evolution of tetrodotoxin resistance in the Sierra garter snake (Thamnophis couchii) by a path less traveled

The repeated evolution of tetrodotoxin (TTX) resistance provides a model for testing hypotheses about the mechanisms of convergent evolution. This poison is broadly employed as a potent antipredator defense, blocking voltage-gated sodium channels (Na_v ) in muscles and nerves, paralyzing and sometimes killing predators. Resistance in taxa bearing this neurotoxin and a few predators appears to come from convergent replacements in specific Na_v residues that interact with TTX. This stereotyped genetic response suggests molecular and phenotypic evolution may be constrained and predictable. Here, we investigate the extent of mechanistic convergence in garter snakes (Thamnophis) that prey on TTX-bearing newts (Taricha) by examining the physiological and genetic basis of TTX resistance in the Sierra garter snake (Th. couchii). We characterize variation in this predatory adaptation across populations at several biological scales: whole-animal TTX resistance; skeletal muscle resistance, functional genetic variation in three Na_v encoding loci; and levels of gene expression for one of these loci. We found Th. couchii possess extensive geographic variation in resistance at the whole-animal and skeletal muscle levels. As in other Thamnophis, resistance at both levels is highly correlated, suggesting convergence across the biological levels linking organism to organ. However, Th. couchii shows no functional variation in Na_v loci among populations or difference in candidate gene expression. Local variation in TTX resistance in Th. couchii cannot be explained by the same relationship between genotype and phenotype seen in other taxa. Thus, historical contingencies may lead different species of Thamnophis down alternative routes to local adaptation.

Every breath you take: assessing metabolic costs of toxin resistance in garter snakes (Thamnophis)

Trait specialization often comes at the expense of original trait function, potentially causing evolutionary tradeoffs that may render specialist populations vulnerable to extinction. However, many specialized adaptations evolve repeatedly, suggesting selection favors specialization in specific environments. Some garter snake (Thamnophis) populations possess specialized mutations in voltage-gated sodium channels that allow them to consume Pacific newts (Taricha) defended by a highly potent neurotoxin (tetrodotoxin). These mutations, however, also decrease protein and muscle function, suggesting garter snakes may suffer evolutionary tradeoffs. We measured a key physiological process, standard metabolic rate (SMR), to investigate whether specialized adaptations in toxin-resistant garter snakes affect baseline energy expenditure. In snakes, skeletal muscles influence metabolism and power ventilation, so inefficiencies of sodium channels in these muscles might impact whole-animal energy expenditure. Further, because sodium channels are membrane-bound proteins, inefficiencies of channel kinetics and performance might be exacerbated at suboptimal temperatures. We measured SMR in 2 species, Thamnophis atratus and Thamnophis sirtalis, that independently evolved tetrodotoxin resistance through unique mutations, providing replicate experiments with distinct underlying genetics and potential physiological costs. Despite our expectations, neither resistance phenotype nor sodium channel genotype affected metabolism and resistant snakes did not perform worse under suboptimal body temperature. Instead, T. atratus and T. sirtalis show nearly identical rates of mass-adjusted energy expenditure at both temperatures, despite differing eco-morphologies, life histories, and distant phylogenetic positions. These findings suggest SMR may be a conserved feature of Thamnophis, and that any organismal tradeoffs may be compensated to retain whole-animal function.

The chemistry and biology of guanidine secondary metabolites

Covering: 2017-2019Guanidine natural products isolated from microorganisms, marine invertebrates and terrestrial plants, amphibians and spiders, represented by non-ribosomal peptides, guanidine-bearing polyketides, alkaloids, terpenoids and shikimic acid derived, are the subject of this review. The topics include the discovery of new metabolites, total synthesis of natural guanidine compounds, biological activity and mechanism-of-action, biosynthesis and ecological functions.

The geographic mosaic of arms race coevolution is closely matched to prey population structure

Reciprocal adaptation is the hallmark of arms race coevolution. Local coadaptation between natural enemies should generate a geographic mosaic pattern where both species have roughly matched abilities across their shared range. However, mosaic variation in ecologically relevant traits can also arise from processes unrelated to reciprocal selection, such as population structure or local environmental conditions. We tested whether these alternative processes can account for trait variation in the geographic mosaic of arms race coevolution between resistant garter snakes (Thamnophis sirtalis) and toxic newts (Taricha granulosa). We found that predator resistance and prey toxin levels are functionally matched in co-occurring populations, suggesting that mosaic variation in the armaments of both species results from the local pressures of reciprocal selection. By the same token, phenotypic and genetic variation in snake resistance deviates from neutral expectations of population genetic differentiation, showing a clear signature of adaptation to local toxin levels in newts. Contrastingly, newt toxin levels are best predicted by genetic differentiation among newt populations, and to a lesser extent, by the local environment and snake resistance. Exaggerated armaments suggest that coevolution occurs in certain hotspots, but prey population structure seems to be of particular influence on local phenotypic variation in both species throughout the geographic mosaic. Our results imply that processes other than reciprocal selection, like historical biogeography and environmental pressures, represent an important source of variation in the geographic mosaic of coevolution. Such a pattern supports the role of "trait remixing" in the geographic mosaic theory, the process by which non-adaptive forces dictate spatial variation in the interactions among species.

The geographic mosaic in parallel: Matching patterns of newt tetrodotoxin levels and snake resistance in multiple predator-prey pairs

The Geographic Mosaic Theory of Coevolution predicts that coevolutionary arms races will vary over time and space because of the diverse ecological settings and population histories of interacting species across the landscape. Thus, understanding coevolution may require investigating broad sets of populations sampled across the range of the interaction. In addition, comparing coevolutionary dynamics between similar systems may reveal the importance of specific factors that structure coevolution. Here, we examine geographic patterns of prey traits and predator traits in the relatively unstudied interaction between the Sierra garter snake (Thamnophis couchii) and sympatric prey, the rough-skinned newt (Taricha granulosa), Sierra newt (Ta. sierrae) and California newt (Ta. torosa). This system parallels, in space and phenotypes, a classic example of coevolution between predatory common garter snakes (Th. sirtalis) and their toxic newt prey exhibiting hotspots of newt tetrodotoxin (TTX) levels and matching snake TTX resistance. We quantified prey and predator traits from hundreds of individuals across their distributions, and functional trait matching at sympatric sites. We show strong regional patterns of trait covariation across the shared ranges of Th. couchii and newt prey. Traits differ significantly among localities, with lower newt TTX levels and snake TTX resistance at the northern latitudes, and higher TTX levels and snake resistance at southern latitudes. Newts and snakes in northern populations show the highest degree of functional trait matching despite possessing the least extreme traits. Conversely, newts and snakes in southern populations show the greatest mismatch despite possessing exaggerated traits, with some snakes so resistant to TTX they would be unaffected by any sympatric newt. Nevertheless, individual variation was substantial, and appears to offer the opportunity for continued reciprocal selection in most populations. Overall, the three species of newts appear to be engaged in a TTX-mediated arms race with Th. couchii. These patterns are congruent with those seen between newts and Th. sirtalis, including the same latitudinal gradient in trait covariation, and the potential 'escape' from the arms race by snake predators. Such concordance in broad scale patterns across two distinct systems suggests common phenomena might structure geographic mosaics in similar ways.

Sex linkage of the skeletal muscle sodium channel gene (SCN4A) explains apparent deviations from Hardy-Weinberg equilibrium of tetrodotoxin-resistance alleles in garter snakes (Thamnophis sirtalis)

The arms race between tetrodotoxin-bearing Pacific newts (Taricha) and their garter snake predators (Thamnophis) in western North America has become a classic example of coevolution, shedding light on predator-prey dynamics, the molecular basis of adaptation, and patterns of convergent evolution. Newts are defended by tetrodotoxin (TTX), a neurotoxin that binds to voltage-gated sodium channels (Na_v proteins), arresting electrical activity in nerves and muscles and paralyzing would-be predators. However, populations of the common garter snake (T. sirtalis) have overcome this defense, largely through polymorphism at the locus SCN4A, which renders the encoded protein (Na_v1.4) less vulnerable to TTX. Previous work suggests that SCN4A commonly shows extreme deviations from Hardy-Weinberg equilibrium (HWE) in these populations, which has been interpreted as the result of intense selection imposed by newts. Here we show that much of this apparent deviation can be attributed to sex linkage of SCN4A. Using genomic data and quantitative PCR, we show that SCN4A is on the Z chromosome in Thamnophis and other advanced snakes. Taking Z-linkage into account, we find that most apparent deviations from HWE can be explained by female hemizygosity rather than low heterozygosity. Sex linkage can affect mutation rates, selection, and drift, and our results suggest that Z-linkage of SCN4A may make significant contributions to the overall dynamics of the coevolutionary arms race between newts and snakes.

Variations in tetrodotoxin levels in populations of Taricha granulosa are expressed in the morphology of their cutaneous glands

Tetrodotoxin (TTX), one of the most toxic substances in nature, is present in bacteria, invertebrates, fishes, and amphibians. Marine organisms seem to bioaccumulate TTX from their food or acquire it from symbiotic bacteria, but its origin in amphibians is unclear. Taricha granulosa can exhibit high TTX levels, presumably concentrated in skin poison glands, acting as an agent of selection upon predatory garter snakes (Thamnophis). This co-evolutionary arms race induces variation in T. granulosa TTX levels, from very high to undetectable. Using morphology and biochemistry, we investigated differences in toxin localization and quality between two populations at the extremes of toxicity. TTX concentration within poison glands is related to the volume of a single cell type in which TTX occurs exclusively in distinctive secretory granules, suggesting a relationship between granule structure and chemical composition. TTX was detected in mucous glands in both populations, contradicting the general understanding that these glands do not secrete defensive chemicals and expanding currently held interpretations of amphibian skin gland functionality. Skin secretions of the two populations differed in low-mass molecules and proteins. Our results demonstrate that interpopulation variation in TTX levels is related to poison gland morphology.

Molecular Adaptations for Sensing and Securing Prey and Insight into Amniote Genome Diversity from the Garter Snake Genome

Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here, we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior, and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.

Large-effect mutations generate trade-off between predatory and locomotor ability during arms race coevolution with deadly prey

Adaptive evolution in response to one selective challenge may disrupt other important aspects of performance. Such evolutionary trade‐offs are predicted to arise in the process of local adaptation, but it is unclear if these phenotypic compromises result from the antagonistic effects of simple amino acid substitutions. We tested for trade‐offs associated with beneficial mutations that confer tetrodotoxin (TTX) resistance in the voltage‐gated sodium channel (Na_V1.4) in skeletal muscle of the common garter snake (Thamnophis sirtalis). Separate lineages in California and the Pacific Northwest independently evolved TTX‐resistant changes to the pore of Na_V1.4 as a result of arms race coevolution with toxic prey, newts of the genus Taricha. Snakes from the California lineage that were homozygous for an allele known to confer large increases in toxin resistance (Na_V1.4^LVNV) had significantly reduced crawl speed compared to individuals with the ancestral TTX‐sensitive channel. Heterologous expression of native snake Na_V1.4 proteins demonstrated that the same Na_V1.4^LVNV allele confers a dramatic increase in TTX resistance and a correlated decrease in overall channel excitability. Our results suggest the same mutations that accumulate during arms race coevolution and beneficially interfere with toxin‐binding also cause changes in electrophysiological function of the channel that may affect organismal performance. This trade‐off was only evident in the predator lineage where coevolution has led to the most extreme resistance phenotype, determined by four critical amino acid substitutions. If these biophysical changes also translate to a fitness cost—for example, through the inability of T. sirtalis to quickly escape predators—then pleiotropy at this single locus could contribute to observed variation in levels of TTX resistance across the mosaic landscape of coevolution.