TETRODOTOXIN RESISTANCE IN GARTER SNAKES: AN EVOLUTIONARY RESPONSE OF PREDATORS TO DANGEROUS PREY

Ecological correlates of cranial evolution in the megaradiation of dipsadine snakes

BACKGROUND

Dipsadine snakes represent one of the most spectacular vertebrate radiations that have occurred in any continental setting, with over 800 species in South and Central America. Their species richness is paralleled by stunning ecological diversity, ranging from arboreal snail-eating and aquatic eel-eating specialists to terrestrial generalists. Despite the ecological importance of this clade, little is known about the extent to which ecological specialization shapes broader patterns of phenotypic diversity within the group. Here, we test how habitat use and diet have influenced morphological diversification in skull shape across 160 dipsadine species using micro-CT and 3-D geometric morphometrics, and we use a phylogenetic comparative approach to test the contributions of habitat use and diet composition to variation in skull shape among species.

RESULTS

We demonstrate that while both habitat use and diet are significant predictors of shape in many regions of the skull, habitat use significantly predicts shape in a greater number of skull regions when compared to diet. We also find that across ecological groupings, fossorial and aquatic behaviors result in the strongest deviations in morphospace for several skull regions. We use simulations to address the robustness of our results and describe statistical anomalies that can arise from the application of phylogenetic generalized least squares to complex shape data.

CONCLUSIONS

Both habitat and dietary ecology are significantly correlated with skull shape in dipsadines; the strongest relationships involved skull shape in snakes with aquatic and fossorial lifestyles. This association between skull morphology and multiple ecological axes is consistent with a classic model of adaptive radiation and suggests that ecological factors were an important component in driving morphological diversification in the dipsadine megaradiation.

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.

Does the Glucocorticoid Stress Response Make Toads More Toxic? An Experimental Study on the Regulation of Bufadienolide Toxin Synthesis

Chemical defense is a crucial component of fitness in many organisms, yet the physiological regulation of defensive toxin synthesis is poorly understood, especially in vertebrates. Bufadienolides, the main defensive compounds of toads, are toxic to many predators and other natural enemies, and their synthesis can be upregulated by stressors, including predation risk, high conspecific density, and pollutants. Thus, higher toxin content may be the consequence of a general endocrine stress response in toads. Therefore, we hypothesized that bufadienolide synthesis may be stimulated by elevated levels of corticosterone (CORT), the main glucocorticoid hormone of amphibians, or by upstream regulators that stimulate CORT production. To test these alternatives, we treated common toad tadpoles with exogenous CORT (exoCORT) or metyrapone (MTP, a CORT-synthesis inhibitor that stimulates upstream regulators of CORT by negative feedback) in the presence or absence of predation cues for 2 or 6 days, and subsequently measured their CORT release rates and bufadienolide content. We found that CORT release rates were elevated by exoCORT, and to a lesser extent also by MTP, regardless of treatment length. Bufadienolide content was significantly decreased by treatment with exoCORT for 6 days but was unaffected by exposure to exoCORT for 2 days or to MTP for either 6 or 2 days. The presence or absence of predation cues affected neither CORT release rate nor bufadienolide content. Our results suggest that changes in bufadienolide synthesis in response to environmental challenges are not driven by CORT but may rather be regulated by upstream hormones of the stress response.

Genomic data reveal local endemism in Southern California Rubber Boas (Serpentes: Boidae, Charina) and the critical need for enhanced conservation actions

The mountains of southern California represent unique, isolated ecosystems that support distinct high-elevation habitats found nowhere else in the area. Analyses of several moisture-dependent species across these sky-islands indicate they exist as locally endemic lineages that occur across these fragmented mountains ranges. The Rubber Boa is a semi-fossorial snake species that is widely distributed in the cooler and more moist ecoregions regions of western North America, including isolated populations across southern California mountain ranges. We developed a genomic and ecological dataset to examine genetic diversity within Rubber Boas and to determine if the endemic Southern Rubber Boa represents a distinct lineage. We quantified current and future habitat suitability under a range of climate change scenarios, and discuss the possible environmental threats facing these unique montane isolates. Our results support four major lineages within Rubber Boas, with genetic breaks that are consistent with biogeographic boundaries observed in other co-distributed, cool-temperature, moisture adapted species. Our data support previous studies that the Southern Rubber Boa is an independent evolutionary unit and now includes multiple locally endemic sky-island populations, restricted to isolated mountain tops and ranges across southern California. Analyses of future habitat suitability indicate that many of these sky-island populations will lose most of their suitable habitat over the next 70 years given predicted increases in drought, rising temperatures, and wildfires. Collectively these data emphasize the critical conservation needs of these montane ecosystems in southern California under current and projected climate change conditions.

Convergent evolution of toxin resistance in animals

Convergence is the phenomenon whereby similar phenotypes evolve independently in different lineages. One example is resistance to toxins in animals. Toxins have evolved many times throughout the tree of life. They disrupt molecular and physiological pathways in target species, thereby incapacitating prey or deterring a predator. In response, molecular resistance has evolved in many species exposed to toxins to counteract their harmful effects. Here, we review current knowledge on the convergence of toxin resistance using examples from a wide range of toxin families. We explore the evolutionary processes and molecular adaptations driving toxin resistance. However, resistance adaptations may carry a fitness cost if they disrupt the normal physiology of the resistant animal. Therefore, there is a trade‐off between maintaining a functional molecular target and reducing toxin susceptibility. There are relatively few solutions that satisfy this trade‐off. As a result, we see a small set of molecular adaptations appearing repeatedly in diverse animal lineages, a phenomenon that is consistent with models of deterministic evolution. Convergence may also explain what has been called ‘autoresistance’. This is often thought to have evolved for self‐protection, but we argue instead that it may be a consequence of poisonous animals feeding on toxic prey. Toxin resistance provides a unique and compelling model system for studying the interplay between trophic interactions, selection pressures and the molecular mechanisms underlying evolutionary novelties.

Venom production and secretion in reptiles

The venom glands of reptiles, particularly those of front-fanged advanced snakes, must satisfy conflicting biological demands: rapid synthesis of potentially labile and highly toxic proteins, storage in the gland lumen for long periods, stabilization of the stored secretions, immediate activation of toxins upon deployment and protection of the animal from the toxic effects of its own venom. This dynamic system could serve as a model for the study of a variety of different phenomena involving exocrine gland activation, protein synthesis, stabilization of protein products and secretory mechanisms. However, these studies have been hampered by a lack of a long-term model that can be propagated in the lab (as opposed to whole-animal studies). Numerous attempts have been made to extend the lifetime of venom gland secretory cells, but only recently has an organoid model been shown to have the requisite qualities of recapitulation of the native system, self-propagation and long-term viability (>1 year). A tractable model is now available for myriad cell- and molecular-level studies of venom glands, protein synthesis and secretion. However, venom glands of reptiles are not identical, and many differ very extensively in overall architecture, microanatomy and protein products produced. This Review summarizes the similarities among and differences between venom glands of helodermatid lizards and of rear-fanged and front-fanged snakes, highlighting those areas that are well understood and identifying areas where future studies can fill in significant gaps in knowledge of these ancient, yet fascinating systems.

Toxic Relationships and Arms-Race Coevolution Revisited

Toxin evolution in animals is one of the most fascinating and complex subjects of scientific inquiry today. Gaining an understanding of toxins poses a multifaceted challenge given the diverse modes of acquisition, evolutionary adaptations, and abiotic components that affect toxin phenotypes. Here, we highlight some of the main genetic and ecological factors that influence toxin evolution and discuss the role of antagonistic interactions and coevolutionary dynamics in shaping the direction and extent of toxicity and resistance in animals. We focus on toxic Pacific newts (family Salamandridae, genus Taricha) as a system to investigate and better evaluate the widely distributed toxin they possess, tetrodotoxin (TTX), and the hypothesized model of arms-race coevolution with snake predators that is used to explain phenotypic patterns of newt toxicity. Finally, we propose an alternative coevolutionary model that incorporates TTX-producing bacteria and draws from an elicitor-receptor concept to explain TTX evolution and ecology.

Chemical Ecology of the North American Newt Genera Taricha and Notophthalmus

The North American newt genera Taricha and Notophthalmus (order Caudata) are well known for the combination of potent toxicity, aposematic coloration, and striking defense postures that protects these animals from predation. This suite of traits is centered around the neurotoxin tetrodotoxin, which causes paralysis and death in metazoans by disrupting the initiation and propagation of electrical signals in the nerves and muscles. Tetrodotoxin defends newts from predation across multiple life history stages and its role in generating arms-race coevolution between Taricha newts and garter snake (genus Thamnophis) predators is well studied. However, understanding the broader picture of chemical defenses in Taricha and Notophthalmus requires an expanded comprehension of the defensive chemical ecology of tetrodotoxin that includes possible coevolutionary interactions with insect egg predators, protection against parasites, as well as mimicry complexes associated with tetrodotoxin and aposematic coloration in both genera. Herein the authors review what is known about the structure, function, and pharmacology of tetrodotoxin to explore its evolution and chemical ecology in the North American newt. Focus is made specifically on the origin and possible biosynthesis of tetrodotoxin in these taxa as well as providing an expanded picture of the web of interactions that contribute to landscape level patterns of toxicity and defense in Taricha and Notophthalmus.

Rapid phenotypic change in a polymorphic salamander over 43 years

Color polymorphic animals offer a unique system for studying intraspecific phenotypic responses to climate change. Discrete color morphs are easy to identify, and correlated trait responses of morphs can indicate how climate warming may facilitate long-term maintenance of polymorphisms. We use a historical dataset spanning 43 years to examine temporal shifts in color morph frequency and body size in response to climate in the Eastern Red-backed Salamander, Plethodon cinereus, which contains a widespread striped/unstriped color polymorphism. We created a pipeline to extract high-throughput trait data from fluid-preserved museum specimens where we batch-photographed salamanders, de-aggregated individual specimens from photographs, and solicited help of community scientists to score color morphs. We used a linear modeling framework that includes information about spatial population structure to demonstrate that color morph frequency and body size vary in response to climate, elevation, and over time, with an overall trend of higher frequency and decreased body size of the striped morph, but increased size of the unstriped morph. These surprising results suggest that morphs may be responding to multiple climate and geographic drivers through co-adapted morphological changes. This work highlights new practices of extracting trait data from museum specimens to demonstrate species phenotypes response to climate change.

Phenotypic and functional variation in venom and venom resistance of two sympatric rattlesnakes and their prey

Predator-prey interactions often lead to the coevolution of adaptations associated with avoiding predation and, for predators, overcoming those defences. Antagonistic coevolutionary relationships are often not simple interactions between a single predator and prey but rather a complex web of interactions between multiple coexisting species. Coevolution between venomous rattlesnakes and small mammals has led to physiological venom resistance in several mammalian taxa. In general, viperid venoms contain large quantities of snake venom metalloproteinase toxins (SVMPs), which are inactivated by SVMP inhibitors expressed in resistant mammals. We explored variation in venom chemistry, SVMP expression, and SVMP resistance across four co-distributed species (California Ground Squirrels, Bryant's Woodrats, Southern Pacific Rattlesnakes, and Red Diamond Rattlesnakes) collected from four different populations in Southern California. Our aim was to understand phenotypic and functional variation in venom and venom resistance in order to compare coevolutionary dynamics of a system involving two sympatric predator-prey pairs to past studies that have focused on single pairs. Proteomic analysis of venoms indicated that these rattlesnakes express different phenotypes when in sympatry, with Red Diamonds expressing more typical viperid venom (with a diversity of SVMPs) and Southern Pacifics expressing a more atypical venom with a broader range of non-enzymatic toxins. We also found that although blood sera from both mammals were generally able to inhibit SVMPs from both rattlesnake species, inhibition depended strongly on the snake population, with snakes from one geographic site expressing SVMPs to which few mammals were resistant. Additionally, we found that Red Diamond venom, rather than woodrat resistance, was locally adapted. Our findings highlight the complexity of coevolutionary relationships between multiple predators and prey that exhibit similar offensive and defensive strategies in sympatry.

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.

Extinction in complex communities as driven by adaptive dynamics

In a complex community, species continuously adapt to each other. On rare occasions, the adaptation of a species can lead to the extinction of others, and even its own. 'Adaptive dynamics' is the standard mathematical framework to describe evolutionary changes in community interactions, and in particular, predict adaptation driven extinction. Unfortunately, most authors implement the equations of adaptive dynamics through computer simulations that require assuming a large number of questionable parameters and fitness functions. In this study, we present analytical solutions to adaptive dynamics equations, thereby clarifying how outcomes depend on any computational input. We develop general formulas that predict equilibrium abundances over evolutionary time scales. Additionally, we predict which species will go extinct next, and when this will happen.

Gene Conversion Facilitates the Adaptive Evolution of Self-Resistance in Highly Toxic Newts

Reconstructing the histories of complex adaptations and identifying the evolutionary mechanisms underlying their origins are two of the primary goals of evolutionary biology. Taricha newts, which contain high concentrations of the deadly toxin tetrodotoxin (TTX) as an antipredator defense, have evolved resistance to self-intoxication, which is a complex adaptation requiring changes in six paralogs of the voltage-gated sodium channel (Na_v) gene family, the physiological target of TTX. Here, we reconstruct the origins of TTX self-resistance by sequencing the entire Na_v gene family in newts and related salamanders. We show that moderate TTX resistance evolved early in the salamander lineage in three of the six Na_v paralogs, preceding the proposed appearance of tetrodotoxic newts by ∼100 My. TTX-bearing newts possess additional unique substitutions across the entire Na_v gene family that provide physiological TTX resistance. These substitutions coincide with signatures of positive selection and relaxed purifying selection, as well as gene conversion events, that together likely facilitated their evolution. We also identify a novel exon duplication within Na_v1.4 encoding an expressed TTX-binding site. Two resistance-conferring changes within newts appear to have spread via nonallelic gene conversion: in one case, one codon was copied between paralogs, and in the second, multiple substitutions were homogenized between the duplicate exons of Na_v1.4. Our results demonstrate that gene conversion can accelerate the coordinated evolution of gene families in response to a common selection pressure.

Stable isotope analysis suggests that tetrodotoxin-resistant Common Gartersnakes (Thamnophis sirtalis) rarely feed on newts in the wild

Toxin-resistant predators may suffer costs from eating chemically defended prey and do not feed exclusively on toxic prey. Common Gartersnakes (Thamnophis sirtalis (Linnaeus, 1758)) have been considered the drivers of an evolutionary arms race with highly toxic newts (genus Taricha Gray, 1850), which they consume with few or no deleterious effects. However, how frequently newts are consumed in nature is less clear. To address this question, we investigated the diets of Th. sirtalis at a site in central Oregon where snakes have high levels of resistance and newts have high levels of tetrodotoxin in the skin. Because snake diets are difficult to quantify using traditional means, we used stable isotopes to estimate the proportion of Th. sirtalis diets made up of newts. Our estimate for the proportion of Th. sirtalis diet made up of Rough-skinned Newts (Taricha granulosa (Skilton, 1849)) at this site is 3.2%. Mole Salamanders (genus Ambystoma Tschudi, 1838) were predicted to be the most important prey, followed by slugs, chorus frogs, and mice, with a very minor role for earthworms. Our results demonstrate that even though Th. sirtalis are physiologically capable of consuming toxic prey, they do not often do so. Generalist predators can be exposed to very strong selection from, and exert reciprocal selection on even rarely eaten, chemically defended prey.

Recovery from discrete wound severities in side-blotched lizards (Uta stansburiana): implications for energy budget, locomotor performance, and oxidative stress

Wounding events (predation attempts, competitive combat) result in injuries and/or infections that induce integrated immune responses for the recovery process. Despite the survival benefits of immunity in this context, the costs incurred may require investment to be diverted from traits contributing to immediate and/or future survival, such as locomotor performance and oxidative status. Yet, whether trait constraints manifest likely depends on wound severity and the implications for energy budget. For this study, food intake, body mass, sprint speed, and oxidative indices (reactive oxygen metabolites, antioxidant capacity) were monitored in male side-blotched lizards (Uta stansburiana) healing from cutaneous wounds of discrete sizes (control, small, large). Results indicate that larger wounds induced faster healing, reduced food consumption, and led to greater oxidative stress over time. Granted wounding did not differentially affect body mass or sprint speed overall, small-wounded lizards with greater wound area healed had faster sprint speeds while large-wounded lizards with greater wound area healed had slower sprint speeds. During recovery from either wound severity, however, healing and sprint performance did not correspond with food consumption, body mass loss, nor oxidative status. These findings provide support that energy budget, locomotor performance, and oxidative status of a reptile are linked to wound recovery to an extent, albeit dependent on wound severity.

Using Integrative Biology to Infer Adaptation from Comparisons of Two (or a Few) Species

AbstractPhylogenetic comparative methods represent a major advance in integrative and comparative biology and have allowed researchers to rigorously test for adaptation in a macroevolutionary framework. However, phylogenetic comparative methods require trait data for many species, which is impractical for certain taxonomic groups and trait types. We propose that the philosophical principle of severity can be implemented in an integrative framework to generate strong inference of adaptation in studies that compare only a few populations or species. This approach requires (1) ensuring that the study system contains species that are relatively closely related; (2) formulating a specific, clear, overarching hypothesis that can be subjected to integrative testing across levels of biological organization (e.g., ecology, behavior, morphology, physiology, and genetics); (3) collecting data that avoid statistical underdetermination and thus allow severe tests of hypotheses; and (4) systematically refining and refuting alternative hypotheses. Although difficult to collect for more than a few species, detailed, integrative data can be used to differentiate among several potential agents of selection. In this way, integrative studies of small numbers of closely related species can complement and even improve on broadscale phylogenetic comparative studies by revealing the specific drivers of adaptation.

Sex-specific growth is mirrored in feeding rate but not moulting frequency in a sexually dimorphic snake

Sexual size dimorphism (SSD), commonly observed in snakes, may arise from a different growth rate between the sexes. This indicates a sex-specific resource intake that is in fact observable in free-living snakes. It is not so well known whether the sexes can express differential feeding rates under conditions unconstrained by spatial accessibility, competition, etc. Here, I studied sex-specific variation in growth, its correlate-moulting frequency, and feeding rate in a captive group of sexually dimorphic banded water snakes (Nerodia fasciata) with access to food unconstrained by predation, competition or space. I showed that the sexes did indeed differ in relative mass growth in that females grew faster than males (p = 0.02), but such differences were not apparent in the moulting rate (p = 0.19). Such differential growth was mirrored in the sex-specific feeding rate, with females ingesting a larger number of meals than males (p = 0.004). Such variation in feeding rate may be governed by an individual's energy expenditure and can be interpreted as a behavioural tendency that contributes to SSD development, independently of other behavioural characteristics. Sex-specific resource demands may drive the differential effects of increasing resource scarcity on both sexes.

Adaptive venom evolution and toxicity in octopods is driven by extensive novel gene formation, expansion, and loss

BACKGROUND

Cephalopods represent a rich system for investigating the genetic basis underlying organismal novelties. This diverse group of specialized predators has evolved many adaptations including proteinaceous venom. Of particular interest is the blue-ringed octopus genus (Hapalochlaena), which are the only octopods known to store large quantities of the potent neurotoxin, tetrodotoxin, within their tissues and venom gland.

FINDINGS

To reveal genomic correlates of organismal novelties, we conducted a comparative study of 3 octopod genomes, including the Southern blue-ringed octopus (Hapalochlaena maculosa). We present the genome of this species and reveal highly dynamic evolutionary patterns at both non-coding and coding organizational levels. Gene family expansions previously reported in Octopus bimaculoides (e.g., zinc finger and cadherins, both associated with neural functions), as well as formation of novel gene families, dominate the genomic landscape in all octopods. Examination of tissue-specific genes in the posterior salivary gland revealed that expression was dominated by serine proteases in non-tetrodotoxin-bearing octopods, while this family was a minor component in H. maculosa. Moreover, voltage-gated sodium channels in H. maculosa contain a resistance mutation found in pufferfish and garter snakes, which is exclusive to the genus. Analysis of the posterior salivary gland microbiome revealed a diverse array of bacterial species, including genera that can produce tetrodotoxin, suggestive of a possible production source.

CONCLUSIONS

We present the first tetrodotoxin-bearing octopod genome H. maculosa, which displays lineage-specific adaptations to tetrodotoxin acquisition. This genome, along with other recently published cephalopod genomes, represents a valuable resource from which future work could advance our understanding of the evolution of genomic novelty in this family.

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.

Predators’ consumption of unpalatable prey does not vary as a function of bitter taste perception

Many prey species contain defensive chemicals that are described as tasting bitter. Bitter taste perception is, therefore, assumed to be important when predators are learning about prey defenses. However, it is not known how individuals differ in their response to bitter taste, and how this influences their foraging decisions. We conducted taste perception assays in which wild-caught great tits (Parus major) were given water with increasing concentrations of bitter-tasting chloroquine diphosphate until they showed an aversive response to bitter taste. This response threshold was found to vary considerably among individuals, ranging from chloroquine concentrations of 0.01 mmol/L to 8 mmol/L. We next investigated whether the response threshold influenced the consumption of defended prey during avoidance learning by presenting birds with novel palatable and defended prey in a random sequence until they refused to attack defended prey. We predicted that individuals with taste response thresholds at lower concentrations would consume fewer defended prey before rejecting them, but found that the response threshold had no effect on the birds’ foraging choices. Instead, willingness to consume defended prey was influenced by the birds’ body condition. This effect was age- and sex-dependent, with adult males attacking more of the defended prey when their body condition was poor, whereas body condition did not have an effect on the foraging choices of juveniles and females. Together, our results suggest that even though taste perception might be important for recognizing prey toxicity, other factors, such as predators’ energetic state, drive the decisions to consume chemically defended prey.

Potential Invasion Risk of Pet Traded Lizards, Snakes, Crocodiles, and Tuatara in the EU on the Basis of a Risk Assessment Model (RAM) and Aquatic Species Invasiveness Screening Kit (AS-ISK)

Because biological invasions can cause many negative impacts, accurate predictions are necessary for implementing effective restrictions aimed at specific high-risk taxa. The pet trade in recent years became the most important pathway for the introduction of non-indigenous species of reptiles worldwide. Therefore, we decided to determine the most common species of lizards, snakes, and crocodiles traded as pets on the basis of market surveys in the Czech Republic, which is an export hub for ornamental animals in the European Union (EU). Subsequently, the establishment and invasion potential for the entire EU was determined for 308 species using proven risk assessment models (RAM, AS-ISK). Species with high establishment potential (determined by RAM) and at the same time with high potential to significantly harm native ecosystems (determined by AS-ISK) included the snakes Thamnophis sirtalis (Colubridae), Morelia spilota (Pythonidae) and also the lizards Tiliqua scincoides (Scincidae) and Intellagama lesueurii (Agamidae).

Serum-based inhibition of pitviper venom by eastern indigo snakes (Drymarchon couperi)

When organisms possess chemical defenses, their predators may eventually evolve resistance to their toxins. Eastern indigo snakes (Drymarchon couperi; EIS) prey on pitvipers and are suspected to possess physiological resistance to their venom. In this study, we formally investigated this hypothesis using microassays that measured the ability of EIS blood sera to inhibit (A) hemolytic and (B) snake venom metalloproteinase (SVMP) activity of copperhead (Agkistrodon contortrix) venom. To serve as controls, we also tested the inhibitory ability of sera from house mice (Mus musculus) and checkered gartersnakes (Thamnophis marcianus), a snake that does not feed on pitvipers. Sera from both EIS and gartersnakes inhibited over 60% of SVMP activity, while only EIS sera also inhibited venom hemolytic activity (78%). Our results demonstrate that EIS serum is indeed capable of inhibiting two of the primary classes of toxins found in copperhead venom, providing the first empirical evidence suggesting that EIS possess physiological resistance to venom upon injection. Because we documented resistance to hemolytic components of pitviper venom within EIS but not gartersnakes, we speculate this resistance may be driven by selection from feeding on pitvipers while resistance to SVMP may be relatively widespread among snakes.

Summary: We confirm eastern indigo snakes possess a physiological resistance to pitviper venom and speculate that venom resistance in snakes may be a product of antagonistic interactions and phylogenetically conserved traits.

Toads phenotypically adjust their chemical defences to anthropogenic habitat change

Despite the well-documented effects of human-induced environmental changes on the morphology, physiology, behaviour and life history of wild animals, next to nothing is known about how anthropogenic habitats influence anti-predatory chemical defence, a crucial fitness component of many species. We investigated the amount and composition of defensive toxins in adult common toads (Bufo bufo) captured in natural, agricultural and urban habitats, and in their offspring raised in a common-garden experiment. We found that, compared to toads captured from natural habitats, adults from both types of anthropogenic habitats had larger toxin glands (parotoids) and their toxin secretion contained higher concentrations of bufagenins, the more potent class of bufadienolide toxins. Furthermore, urban toads had lower concentrations of bufotoxins, the compounds with lower toxicity. None of these differences were present in the captive-raised juveniles; instead, toadlets originating from agricultural habitats had smaller parotoids and lower bufotoxin concentrations. These results suggest that toads' chemical defences respond to the challenges of anthropogenic environments via phenotypic plasticity. These responses may constitute non-adaptive consequences of pollution by endocrine-disrupting chemicals as well as adaptive adjustments to the altered predator assemblages of urban and agricultural habitats.

The evolution and origin of tetrodotoxin acquisition in the blue-ringed octopus (genus Hapalochlaena)

Tetrodotoxin is a potent non-proteinaceous neurotoxin, which is commonly found in the marine environment. Synthesised by bacteria, tetrodotoxin has been isolated from the tissues of several genera including pufferfish, salamanders and octopus. Believed to provide a defensive function, the independent evolution of tetrodotoxin sequestration is poorly understood in most species. Two mechanisms of tetrodotoxin resistance have been identified to date, tetrodotoxin binding proteins in the circulatory system and mutations to voltage gated sodium channels, the binding target of tetrodotoxin with the former potentially succeeding the latter in evolutionary time. This review focuses on the evolution of tetrodotoxin acquisition, in particular how it may have occurred within the blue-ringed octopus genus (Hapalochlaena) and the subsequent impact on venom evolution.

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.

Of poisons and parasites-the defensive role of tetrodotoxin against infections in newts

Classical research on animal toxicity has focused on the role of toxins in protection against predators, but recent studies suggest these same compounds can offer a powerful defense against parasites and infectious diseases. Newts in the genus Taricha are brightly coloured and contain the potent neurotoxin, tetrodotoxin (TTX), which is hypothesized to have evolved as a defense against vertebrate predators such as garter snakes. However, newt populations often vary dramatically in toxicity, which is only partially explained by predation pressure. The primary aim of this study was to evaluate the relationships between TTX concentration and infection by parasites. By systematically assessing micro- and macroparasite infections among 345 adult newts (sympatric populations of Taricha granulosa and T. torosa), we detected 18 unique taxa of helminths, fungi, viruses and protozoans. For both newt species, per-host concentrations of TTX, which varied from undetectable to >60 μg/cm² skin, negatively predicted overall parasite richness as well as the likelihood of infection by the chytrid fungus, Batrachochytrium dendrobatidis, and ranavirus. No such effect was found on infection load among infected hosts. Despite commonly occurring at the same wetlands, T. torosa supported higher parasite richness and average infection load than T. granulosa. Host body size and sex (females > males) tended to positively predict infection levels in both species. For hosts in which we quantified leucocyte profiles, total white blood cell count correlated positively with both parasite richness and total infection load. By coupling data on host toxicity and infection by a broad range of micro- and macroparasites, these results suggest that-alongside its effects on predators-tetrodotoxin may help protect newts against parasitic infections, highlighting the importance of integrative research on animal chemistry, immunological defenses and natural enemy ecology.

Noxious newts and their natural enemies: Experimental effects of tetrodotoxin exposure on trematode parasites and aquatic macroinvertebrates

The dermal glands of many amphibian species secrete toxins or other noxious substances as a defense strategy against natural enemies. Newts in particular possess the potent neurotoxin tetrodotoxin (TTX), for which the highest concentrations are found in species within the genus Taricha. Adult Taricha are hypothesized to use TTX as a chemical defense against vertebrate predators such as garter snakes (Thamnophis spp.). However, less is known about how TTX functions to defend aquatic-developing newt larvae against natural enemies, including trematode parasites and aquatic macroinvertebrates. Here we experimentally investigated the effects of exogenous TTX exposure on survivorship of the infectious stages (cercariae) of five species of trematode parasites that infect larval amphibians. Specifically, we used dose-response curves to test the sensitivity of trematode cercariae to progressively increasing concentrations of TTX (0.0 [control], 0.63, 3.13, 6.26, 31.32, and 62.64 nmol L-1) and how this differed among parasite species. We further compared these results to the effects of TTX exposure (0 and 1000 nmolL-1) over 24 h on seven macroinvertebrate taxa commonly found in aquatic habitats with newt larvae. TTX significantly reduced the survivorship of trematode cercariae for all species, but the magnitude of such effects varied among species. Ribeiroia ondatrae - which causes mortality and limb malformations in amphibians - was the least sensitive to TTX, whereas the kidney-encysting Echinostoma trivolvis was the most sensitive. Among the macroinvertebrate taxa, only mayflies (Ephemeroptera) showed a significant increase in mortality following exogenous TTX exposure, despite the use of a concentration 16x higher than the maximum used for trematodes. Our results suggest that maternal investment of TTX into larval newts may provide protection against certain trematode infections and highlight the importance of future work assessing the effects of newt toxicity on both parasite infection success and the palatability of larval newts to invertebrate predators.

Individual fluctuations in toxin levels affect breeding site fidelity in a chemically defended amphibian

Behaviours that influence habitat selection strongly determine species movement patterns. One component of animal behaviour that largely influences movement patterns and habitat choice is site fidelity. California newts (family Salamandridae) demonstrate remarkable site fidelity, typically homing to the same pool of a stream each breeding season. Individuals often occupy a specific pool throughout the breeding season, but some males shift among breeding pools, altering their set of potential mates, competitors, and predators. In this study, we measured dermal concentrations of the chemical defence compound tetrodotoxin (TTX) in recaptured male California newts (Taricha torosa) over five breeding seasons to evaluate whether relative TTX concentrations are associated with breeding site fidelity in the field. Our five years of field sampling indicates that TTX concentrations of individuals and group means fluctuate tremendously, implying that TTX is not a stable phenotypic trait. Despite such fluctuations, we found that an individual's relative TTX concentration explains fidelity to a breeding pool and suggests that newts may be able to assess both their own concentrations of TTX and that of conspecifics to make decisions about remaining in or abandoning a breeding pool. These results provide us a novel dimension to chemical defence phenotypes in nature and their ecological consequences, potentially requiring a re-evaluation of the coevolutionary dynamics of predation pressure on toxin-laden organisms.

Stressed snakes strike first: Hormone levels and defensive behavior in free ranging cottonmouths (Agkistrodon piscivorus)

Stress is believed to be an important factor mediating animal behavior. Here we explore the relationship between concentrations of a stress hormone and defensive behavior of a snake. The cottonmouth (Agkistrodon piscivorus) is an abundant, large-bodied pitviper that is well known for its intense defensive behaviors. The defensive behaviors and hormonal ecology of cottonmouths have been studied extensively, but the interaction between these is not well understood. We conducted field trials, recording the snake's behavior and obtaining blood samples to quantify plasma CORT concentrations, both upon first encountering a snake and after a 30min standardized confinement stressor. We found that snakes with elevated levels of baseline CORT at first encounter were more likely to strike than exhibit a threat display when approached in the field. However, this behavior was not related to the magnitude of the snake's CORT increase following confinement, suggesting that more stress-prone snakes are not more defensive. Post-stressor antipredator behavior was also not related to any of our CORT measures. This study suggests that baseline CORT levels can be important correlates of defensive behavior. If this is a causative relationship, environmental challenges that increase baseline stress levels of populations may elevate cottonmouth defensive behavior. This would increase costs associated with defensive behavior (energetic, lost opportunity, etc.) and have important consequences for animal-human interactions.

Is there more than one way to skin a newt? Convergent toxin resistance in snakes is not due to a common genetic mechanism

Convergent evolution of tetrodotoxin (TTX) resistance, at both the phenotypic and genetic levels, characterizes coevolutionary arms races between amphibians and their snake predators around the world, and reveals remarkable predictability in the process of adaptation. Here we examine the repeatability of the evolution of TTX resistance in an undescribed predator-prey relationship between TTX-bearing Eastern Newts (Notophthalmus viridescens) and Eastern Hog-nosed Snakes (Heterodon platirhinos). We found that that local newts contain levels of TTX dangerous enough to dissuade most predators, and that Eastern Hog-nosed Snakes within newt range are highly resistant to TTX. In fact, these populations of Eastern Hog-nosed Snakes are so resistant to TTX that the potential for current reciprocal selection might be limited. Unlike all other cases of TTX resistance in vertebrates, H. platirhinos lacks the adaptive amino acid substitutions in the skeletal muscle sodium channel that reduce TTX binding, suggesting that physiological resistance in Eastern Hog-nosed Snakes is conferred by an alternate genetic mechanism. Thus, phenotypic convergence in this case is not due to parallel molecular evolution, indicating that there may be more than one way for this adaptation to arise, even among closely related species.