Convergent evolution of toxin resistance in animals

Sugar-coated survival: N-glycosylation as a unique bearded dragon venom resistance trait within Australian agamid lizards

In the ongoing evolutionary arms race between predators and prey, adaptive innovations often trigger a reciprocal response. For instance, the emergence of α-neurotoxins in snake venom has driven prey species targeted by these snakes to evolve sophisticated defense mechanisms. This study zeroes in on the particular motifs within the orthosteric sites of post-synaptic nicotinic acetylcholine receptors (nAChR) that confer resistance to α-neurotoxins, often through structural alterations of nAChR. This research examined Australian agamid lizards, a primary prey group for Australian elapid snakes, which are subject to predatory selection pressures. We previously showed that Pogona vitticeps (Central bearded dragon) was resistant to α-neurotoxic snake venoms through a steric hindrance form resistance evolving within the nAChR orthosteric, specifically through the 187-189NVT motif resulting in the presence of N-glycosylation, with the branching carbohydrate chains impeding the binding by the neurotoxins. This adaptive trait is thought to be a compensatory mechanism for the lizard's limited escape capabilities. Despite the significance of this novel adaptation, the prevalence and evolutionary roots of such venom resistance in Australian agamids have not been thoroughly investigated. To fill this knowledge gap, we undertook a comprehensive sequencing analysis of the nAChR ligand-binding domain across the full taxonomical diversity of Australian agamid species. Our findings reveal that the N-glycosylation resistance mechanism is a trait unique to the Pogona genus and absent in other Australian agamids. This aligns with Pogona's distinctive morphology, which likely increases vulnerability to neurotoxic elapid snakes, thereby increasing selective pressures for resistance. In contrast, biolayer interferometry experiments with death adder (Acanthophis species) venoms did not indicate any resistance-related binding patterns in other agamids, suggesting a lack of similar resistance adaptations, consistent with these lineages either being fast-moving, covered with large defensive spines, or being arboreal. This research not only uncovers a novel α-neurotoxin resistance mechanism in Australian agamids but also highlights the complex dynamics of the predator-prey chemical arms race. It provides a deeper understanding of how evolutionary pressures shape the interactions between venomous snakes and their prey.

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.

Studying Venom Toxin Variation Using Accurate Masses from Liquid Chromatography-Mass Spectrometry Coupled with Bioinformatic Tools

This study provides a new methodology for the rapid analysis of numerous venom samples in an automated fashion. Here, we use LC-MS (Liquid Chromatography-Mass Spectrometry) for venom separation and toxin analysis at the accurate mass level combined with new in-house written bioinformatic scripts to obtain high-throughput results. This analytical methodology was validated using 31 venoms from all members of a monophyletic clade of Australian elapids: brown snakes (Pseudonaja spp.) and taipans (Oxyuranus spp.). In a previous study, we revealed extensive venom variation within this clade, but the data was manually processed and MS peaks were integrated into a time-consuming and labour-intensive approach. By comparing the manual approach to our new automated approach, we now present a faster and more efficient pipeline for analysing venom variation. Pooled venom separations with post-column toxin fractionations were performed for subsequent high-throughput venomics to obtain toxin IDs correlating to accurate masses for all fractionated toxins. This workflow adds another dimension to the field of venom analysis by providing opportunities to rapidly perform in-depth studies on venom variation. Our pipeline opens new possibilities for studying animal venoms as evolutionary model systems and investigating venom variation to aid in the development of better antivenoms.

High-Voltage Toxin'Roll: Electrostatic Charge Repulsion as a Dynamic Venom Resistance Trait in Pythonid Snakes

The evolutionary interplay between predator and prey has significantly shaped the development of snake venom, a critical adaptation for subduing prey. This arms race has spurred the diversification of the components of venom and the corresponding emergence of resistance mechanisms in the prey and predators of venomous snakes. Our study investigates the molecular basis of venom resistance in pythons, focusing on electrostatic charge repulsion as a defense against α-neurotoxins binding to the alpha-1 subunit of the postsynaptic nicotinic acetylcholine receptor. Through phylogenetic and bioactivity analyses of orthosteric site sequences from various python species, we explore the prevalence and evolution of amino acid substitutions that confer resistance by electrostatic repulsion, which initially evolved in response to predatory pressure by Naja (cobra) species (which occurs across Africa and Asia). The small African species Python regius retains the two resistance-conferring lysines (positions 189 and 191) of the ancestral Python genus, conferring resistance to sympatric Naja venoms. This differed from the giant African species Python sebae, which has secondarily lost one of these lysines, potentially due to its rapid growth out of the prey size range of sympatric Naja species. In contrast, the two Asian species Python brongersmai (small) and Python bivittatus (giant) share an identical orthosteric site, which exhibits the highest degree of resistance, attributed to three lysine residues in the orthosteric sites. One of these lysines (at orthosteric position 195) evolved in the last common ancestor of these two species, which may reflect an adaptive response to increased predation pressures from the sympatric α-neurotoxic snake-eating genus Ophiophagus (King Cobras) in Asia. All these terrestrial Python species, however, were less neurotoxin-susceptible than pythons in other genera which have evolved under different predatory pressure as: the Asian species Malayopython reticulatus which is arboreal as neonates and juveniles before rapidly reaching sizes as terrestrial adults too large for sympatric Ophiophagus species to consider as prey; and the terrestrial Australian species Aspidites melanocephalus which occupies a niche, devoid of selection pressure from α-neurotoxic predatory snakes. Our findings underline the importance of positive selection in the evolution of venom resistance and suggest a complex evolutionary history involving both conserved traits and secondary evolution. This study enhances our understanding of the molecular adaptations that enable pythons to survive in environments laden with venomous threats and offers insights into the ongoing co-evolution between venomous snakes and their prey.

Tissue damaging toxins in snake venoms: mechanisms of action, pathophysiology and treatment strategies

Snakebite envenoming is an important public health issue responsible for mortality and severe morbidity. Where mortality is mainly caused by venom toxins that induce cardiovascular disturbances, neurotoxicity, and acute kidney injury, morbidity is caused by toxins that directly or indirectly destroy cells and degrade the extracellular matrix. These are referred to as 'tissue-damaging toxins' and have previously been classified in various ways, most of which are based on the tissues being affected (e.g., cardiotoxins, myotoxins). This categorisation, however, is primarily phenomenological and not mechanistic. In this review, we propose an alternative way of classifying cytotoxins based on their mechanistic effects rather than using a description that is organ- or tissue-based. The mechanisms of toxin-induced tissue damage and their clinical implications are discussed. This review contributes to our understanding of fundamental biological processes associated with snakebite envenoming, which may pave the way for a knowledge-based search for novel therapeutic options.

A plethora of rodents: Rattlesnake predators generate unanticipated patterns of venom resistance in a grassland ecosystem

Predation has the potential to impart strong selective pressures on organisms within their environments, resulting in adaptive changes in prey that minimize risk of predation. Pressures from venomous snakes present an exceptional challenge to prey, as venom represents a unique chemical arsenal evolutionarily tailored to incapacitate prey. In response, venom resistance has been detected in various snake prey species, and to varying degrees. This study analyzes venom resistance in an eastern Colorado grassland habitat, where the Prairie Rattlesnake (Crotalus viridis) and Desert Massasauga Rattlesnake (Sistrurus tergeminus edwardsii) co-occur with a suite of grassland rodents. We test for venom resistance across rodent and snake pairings using two geographically distant field sites to determine the role of 1) predation pressure and trophic ecology, and 2) sympatric and allopatric patterns of venom resistance. Resistance was measured using serum-based metalloproteinase inhibition assays to determine potential inhibition of proteolytic activity, augmented by median lethal dose (LD50) assays on rodent species to assess toxicity of crude venoms. Resistance is present in several rodent species, with strong resistance present in populations of Eastern Woodrat (Neotoma floridana), Ord's Kangaroo Rat (Dipodomys ordii), and Northern Grasshopper Mouse (Onychomys leucogaster). Resistance is less developed in other species, including the House Mouse (Mus musculus) and Plains Pocket Mouse (Perognathus flavescens). An unexpected differential is present, where Lincoln County Kangaroo Rats are highly resistant to venom of co-occurring Prairie Rattlesnakes yet are sensitive to an allopatric population of Prairie Rattlesnakes in Weld County. Lincoln Co. Northern Grasshopper Mice also demonstrate extremely elevated resistance to Weld Co. Prairie Rattlesnake venoms, and they may possess resistance mechanisms for myotoxin a, an abundant component of Weld Co. C. v viridis venoms. This study illustrates the complexity of venom resistance in biological communities that can exist when incorporating multiple species interactions. Future studies aimed at characterizing resistance mechanisms at the molecular level will provide a more detailed physiological context for understanding mechanisms by which resistance to venoms occurs.

Parasitic fish embryos do a "front-flip" on the yolk to resist expulsion from the host

Embryonic development is often considered shielded from the effects of natural selection, being selected primarily for reliable development. However, embryos sometimes represent virulent parasites, triggering a coevolutionary "arms race" with their host. We have examined embryonic adaptations to a parasitic lifestyle in the bitterling fish. Bitterlings are brood parasites that lay their eggs in the gill chamber of host mussels. Bitterling eggs and embryos have adaptations to resist being flushed out by the mussel. These include a pair of projections from the yolk sac that act as an anchor. Furthermore, bitterling eggs all adopt a head-down position in the mussel gills which further increases their chances of survival. To examine these adaptations in detail, we have studied development in the rosy bitterling (Rhodeus ocellatus) using molecular markers, X-ray tomography, and time-lapse imaging. We describe a suite of developmental adaptations to brood parasitism in this species. We show that the mechanism underlying these adaptions is a modified pattern of blastokinesis-a process unique, among fish, to bitterlings. Tissue movements during blastokinesis cause the embryo to do an extraordinary "front-flip" on the yolk. We suggest that this movement determines the spatial orientation of the other developmental adaptations to parasitism, ensuring that they are optimally positioned to help resist the ejection of the embryo from the mussel. Our study supports the notion that natural selection can drive the evolution of a suite of adaptations, both embryonic and extra-embryonic, via modifications in early development.

A Russian Doll of Resistance: Nested Gains and Losses of Venom Immunity in Varanid Lizards

The interplay between predator and prey has catalyzed the evolution of venom systems, with predators honing their venoms in response to the evolving resistance of prey. A previous study showed that the African varanid species Varanus exanthematicus has heightened resistance to snake venoms compared to the Australian species V. giganteus, V. komodoensis, and V. mertensi, likely due to increased predation by sympatric venomous snakes on V. exanthematicus. To understand venom resistance among varanid lizards, we analyzed the receptor site targeted by venoms in 27 varanid lizards, including 25 Australian varanids. The results indicate an active evolutionary arms race between Australian varanid lizards and sympatric neurotoxic elapid snakes. Large species preying on venomous snakes exhibit inherited neurotoxin resistance, a trait potentially linked to their predatory habits. Consistent with the 'use it or lose it' aspect of venom resistance, this trait was secondarily reduced in two lineages that had convergently evolved gigantism (V. giganteus and the V. komodoensis/V. varius clade), suggestive of increased predatory success accompanying extreme size and also increased mechanical protection against envenomation due to larger scale osteoderms. Resistance was completely lost in the mangrove monitor V. indicus, consistent with venomous snakes not being common in their arboreal and aquatic niche. Conversely, dwarf varanids demonstrate a secondary loss at the base of the clade, with resistance subsequently re-evolving in the burrowing V. acanthurus/V. storri clade, suggesting an ongoing battle with neurotoxic predators. Intriguingly, within the V. acanthurus/V. storri clade, resistance was lost again in V. kingorum, which is morphologically and ecologically distinct from other members of this clade. Resistance was also re-evolved in V. glebopalma which is terrestrial in contrast to the arboreal/cliff dwelling niches occupied by the other members of its clade (V. glebopalma, V. mitchelli, V. scalaris, V. tristis). This 'Russian doll' pattern of venom resistance underscores the dynamic interaction between dwarf varanids and Australian neurotoxic elapid snakes. Our research, which included testing Acanthophis (death adder) venoms against varanid receptors as models for alpha-neurotoxic interactions, uncovered a fascinating instance of the Red Queen Hypothesis: some death adders have developed more potent toxins specifically targeting resistant varanids, a clear sign of the relentless predator-prey arms race. These results offer new insight into the complex dynamics of venom resistance and highlight the intricate ecological interactions that shape the natural world.

Sodium Channel β Subunits-An Additional Element in Animal Tetrodotoxin Resistance?

Tetrodotoxin (TTX) is a neurotoxic molecule used by many animals for defense and/or predation, as well as an important biomedical tool. Its ubiquity as a defensive agent has led to repeated independent evolution of tetrodotoxin resistance in animals. TTX binds to voltage-gated sodium channels (VGSC) consisting of α and β subunits. Virtually all studies investigating the mechanisms behind TTX resistance have focused on the α subunit of voltage-gated sodium channels, where tetrodotoxin binds. However, the possibility of β subunits also contributing to tetrodotoxin resistance was never explored, though these subunits act in concert. In this study, we present preliminary evidence suggesting a potential role of β subunits in the evolution of TTX resistance. We gathered mRNA sequences for all β subunit types found in vertebrates across 12 species (three TTX-resistant and nine TTX-sensitive) and tested for signatures of positive selection with a maximum likelihood approach. Our results revealed several sites experiencing positive selection in TTX-resistant taxa, though none were exclusive to those species in subunit β1, which forms a complex with the main physiological target of TTX (VGSC Nav1.4). While experimental data validating these findings would be necessary, this work suggests that deeper investigation into β subunits as potential players in tetrodotoxin resistance may be worthwhile.

Toxicogenomic Effects of Dissolved Saxitoxin on the Early Life Stages of the Longfin Yellowtail (Seriola rivoliana)

Harmful algal blooms (HABs) can produce a variety of noxious effects and, in some cases, the massive mortality of wild and farmed marine organisms. Some HAB species produce toxins that are released into seawater or transferred via food webs (particulate toxin fraction). The objective of the present study was to identify the toxicological effects of subacute exposure to saxitoxin (STX) during embryonic and early larval stages in Seriola rivoliana. Eggs were exposed to dissolved 19 STX (100 μg L-1). The toxic effects of STX were evaluated via the hatching percentage, the activity of three enzymes (protein and alkaline phosphatases and peroxidase), and the expression of four genes (HSF2, Nav1.4b, PPRC1, and DUSP8). A low hatching percentage (less than 5%) was observed in 44 hpf (hours post fertilization) embryos exposed to STX compared to 71% in the unexposed control. At this STX concentration, no oxidative stress in the embryos was evident. However, STX induced the expression of the NaV1.4 channel α-subunit (NaV1.4b), which is the primary target of this toxin. Our results revealed the overexpression of all four candidate genes in STX-intoxicated lecithotrophic larvae, reflecting the activation of diverse cellular processes involved in stress responses (HSF2), lipid metabolism (PPRC1), and MAP kinase signaling pathways associated with cell proliferation and differentiation (DUSP8). The effects of STX were more pronounced in young larvae than in embryos, indicating a stage-specific sensitivity to the toxin.

Neurotoxin-Derived Optical Probes for Biological and Medical Imaging

The superb specificity and potency of biological toxins targeting various ion channels and receptors are of major interest for the delivery of therapeutics to distinct cell types and subcellular compartments. Fused with reporter proteins or labelled with fluorophores and nanocomposites, animal toxins and their detoxified variants also offer expanding opportunities for visualisation of a range of molecular processes and functions in preclinical models, as well as clinical studies. This article presents state-of-the-art optical probes derived from neurotoxins targeting ion channels, with discussions of their applications in basic and translational biomedical research. It describes the design and production of probes and reviews their applications with advantages and limitations, with prospects for future improvements. Given the advances in imaging tools and expanding research areas benefiting from the use of optical probes, described here resources should assist the discovery process and facilitate high-precision interrogation and therapeutic interventions.

One Size Fits All-Venomics of the Iberian Adder (Vipera seoanei, Lataste 1878) Reveals Low Levels of Venom Variation across Its Distributional Range

European vipers (genus Vipera) are medically important snakes displaying considerable venom variation, occurring at different levels in this group. The presence of intraspecific venom variation, however, remains understudied in several Vipera species. Vipera seoanei is a venomous snake endemic to the northern Iberian Peninsula and south-western France, presenting notable phenotypic variation and inhabiting several diverse habitats across its range. We analysed the venoms of 49 adult specimens of V. seoanei from 20 localities across the species' Iberian distribution. We used a pool of all individual venoms to generate a V. seoanei venom reference proteome, produced SDS-PAGE profiles of all venom samples, and visualised patterns of variation using NMDS. By applying linear regression, we then assessed presence and nature of venom variation between localities, and investigated the effect of 14 predictors (biological, eco-geographic, genetic) on its occurrence. The venom comprised at least 12 different toxin families, of which five (i.e., PLA2, svSP, DI, snaclec, svMP) accounted for about 75% of the whole proteome. The comparative analyses of the SDS-PAGE venom profiles showed them to be remarkably similar across the sampled localities, suggesting low geographic variability. The regression analyses suggested significant effects of biological and habitat predictors on the little variation we detected across the analysed V. seoanei venoms. Other factors were also significantly associated with the presence/absence of individual bands in the SDS-PAGE profiles. The low levels of venom variability we detected within V. seoanei might be the result of a recent population expansion, or of processes other than directional positive selection.

Evolutionary repeatability of emergent properties of ecological communities

Most species belong to ecological communities where their interactions give rise to emergent community-level properties, such as diversity and productivity. Understanding and predicting how these properties change over time has been a major goal in ecology, with important practical implications for sustainability and human health. Less attention has been paid to the fact that community-level properties can also change because member species evolve. Yet, our ability to predict long-term eco-evolutionary dynamics hinges on how repeatably community-level properties change as a result of species evolution. Here, we review studies of evolution of both natural and experimental communities and make the case that community-level properties at least sometimes evolve repeatably. We discuss challenges faced in investigations of evolutionary repeatability. In particular, only a handful of studies enable us to quantify repeatability. We argue that quantifying repeatability at the community level is critical for approaching what we see as three major open questions in the field: (i) Is the observed degree of repeatability surprising? (ii) How is evolutionary repeatability at the community level related to repeatability at the level of traits of member species? (iii) What factors affect repeatability? We outline some theoretical and empirical approaches to addressing these questions. Advances in these directions will not only enrich our basic understanding of evolution and ecology but will also help us predict eco-evolutionary dynamics. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Female-female aggression in the Gila monster (Heloderma suspectum)

Historically, the role of aggression in the social lives of animals overwhelmingly focused on males. In recent years, however, female-female aggression in vertebrates, particularly lizards, has received increasing attention. This growing body of literature shows both similarities and differences to aggressive behaviours between males. Here, we document female-female aggression in captive Gila monsters (Heloderma suspectum). Based on four unique dyadic trials (eight adult female subjects), we developed a qualitative ethogram. Unexpected and most intriguing were the prevalence and intensity of aggressive acts that included brief and sustained biting, envenomation, and lateral rotation (i.e. rolling of body while holding onto opponent with closed jaws). Given specific behavioural acts (i.e. biting) and the results of bite-force experiments, we postulate that osteoderms (bony deposits in the skin) offer some degree of protection and reduce the likelihood of serious injury during female-female fights. Male-male contests in H. suspectum, in contrast, are more ritualized, and biting is rarely reported. Female-female aggression in other lizards has a role in territoriality, courtship tactics, and nest and offspring guarding. Future behavioural research on aggression in female Gila monsters is warranted to test these and other hypotheses in the laboratory and field.

Historical Perspective of the Characterization of Conotoxins Targeting Voltage-Gated Sodium Channels

Marine toxins have potent actions on diverse sodium ion channels regulated by transmembrane voltage (voltage-gated ion channels) or by neurotransmitters (nicotinic acetylcholine receptor channels). Studies of these toxins have focused on varied aspects of venom peptides ranging from evolutionary relationships of predator and prey, biological actions on excitable tissues, potential application as pharmacological intervention in disease therapy, and as part of multiple experimental approaches towards an understanding of the atomistic characterization of ion channel structure. This review examines the historical perspective of the study of conotoxin peptides active on sodium channels gated by transmembrane voltage, which has led to recent advances in ion channel research made possible with the exploitation of the diversity of these marine toxins.

Highly Evolvable: Investigating Interspecific and Intraspecific Venom Variation in Taipans (Oxyuranus spp.) and Brown Snakes (Pseudonaja spp.)

Snake venoms are complex mixtures of toxins that differ on interspecific (between species) and intraspecific (within species) levels. Whether venom variation within a group of closely related species is explained by the presence, absence and/or relative abundances of venom toxins remains largely unknown. Taipans (Oxyuranus spp.) and brown snakes (Pseudonaja spp.) represent medically relevant species of snakes across the Australasian region and provide an excellent model clade for studying interspecific and intraspecific venom variation. Using liquid chromatography with ultraviolet and mass spectrometry detection, we analyzed a total of 31 venoms covering all species of this monophyletic clade, including widespread localities. Our results reveal major interspecific and intraspecific venom variation in Oxyuranus and Pseudonaja species, partially corresponding with their geographical regions and phylogenetic relationships. This extensive venom variability is generated by a combination of the absence/presence and differential abundance of venom toxins. Our study highlights that venom systems can be highly dynamical on the interspecific and intraspecific levels and underscores that the rapid toxin evolvability potentially causes major impacts on neglected tropical snakebites.

Venom biotechnology: casting light on nature's deadliest weapons using synthetic biology

Venoms are complex chemical arsenals that have evolved independently many times in the animal kingdom. Venoms have attracted the interest of researchers because they are an important innovation that has contributed greatly to the evolutionary success of many animals, and their medical relevance offers significant potential for drug discovery. During the last decade, venom research has been revolutionized by the application of systems biology, giving rise to a novel field known as venomics. More recently, biotechnology has also made an increasing impact in this field. Its methods provide the means to disentangle and study venom systems across all levels of biological organization and, given their tremendous impact on the life sciences, these pivotal tools greatly facilitate the coherent understanding of venom system organization, development, biochemistry, and therapeutic activity. Even so, we lack a comprehensive overview of major advances achieved by applying biotechnology to venom systems. This review therefore considers the methods, insights, and potential future developments of biotechnological applications in the field of venom research. We follow the levels of biological organization and structure, starting with the methods used to study the genomic blueprint and genetic machinery of venoms, followed gene products and their functional phenotypes. We argue that biotechnology can answer some of the most urgent questions in venom research, particularly when multiple approaches are combined together, and with other venomics technologies.

A novel broad spectrum venom metalloproteinase autoinhibitor in the rattlesnake Crotalus atrox evolved via a shift in paralog function

The complexity of snake venom composition reflects adaptation to the diversity of prey and may be driven at times by a coevolutionary arms race between snakes and venom-resistant prey. However, many snakes are also resistant to their own venom due to serum-borne inhibitors of venom toxins, which raises the question of how snake autoinhibitors maintain their efficacy as venom proteins evolve. To investigate this potential three-way arms race among venom, prey, and autoinhibitors, we have identified and traced the evolutionary origin of serum inhibitors of snake venom metalloproteinases (SVMPs) in the Western Diamondback rattlesnake Crotalus atrox which possesses the largest known battery of SVMP genes among crotalids examined. We found that C. atrox expresses five members of a Fetuin A-related metalloproteinase inhibitor family but that one family member, FETUA-3, is the major SVMP inhibitor that binds to approximately 20 different C. atrox SVMPs and inhibits activities of all three SVMP classes. We show that the fetua-3 gene arose deep within crotalid evolution before the origin of New World species but, surprisingly, fetua-3 belongs to a different paralog group than previously identified SVMP inhibitors in Asian and South American crotalids. Conversely, the C. atrox FETUA-2 ortholog of previously characterized crotalid SVMP inhibitors shows limited activity against C. atrox SVMPs. These results reveal that there has been a functional evolutionary shift in the major SVMP inhibitor in the C. atrox lineage as the SVMP family expanded and diversified in the Crotalus lineage. This broad-spectrum inhibitor may be of potential therapeutic interest.

In Silico Analysis of Tetrodotoxin Binding in Voltage-Gated Sodium Ion Channels from Toxin-Resistant Animal Lineages

Multiple animal species have evolved resistance to the neurotoxin tetrodotoxin (TTX) through changes in voltage-gated sodium ion channels (VGSCs). Amino acid substitutions in TTX-resistant lineages appear to be positionally convergent with changes in homologous residues associated with reductions in TTX block. We used homology modeling coupled with docking simulations to test whether positionally convergent substitutions generate functional convergence at the level of TTX–channel interactions. We found little evidence that amino acids at convergent positions generated similar patterns among TTX-resistant animal lineages across several metrics, including number of polar contacts, polar contact position, and estimates of binding energy. Though binding energy values calculated for TTX docking were reduced for some TTX-resistant channels, not all TTX-resistant channels and not all of our analyses returned reduced binding energy values for TTX-resistant channels. Our results do not support a simple model of toxin resistance where a reduced number of bonds between TTX and the channel protein prevents blocking. Rather models that incorporate flexibility and movement of the protein overall may better describe how homologous substitutions in the channel cause changes in TTX block.

Widespread convergent evolution of alpha-neurotoxin resistance in African mammals

Convergent evolution is central to the study of adaptation and has been used to understand both the limits of evolution and the diverse patterns and processes which result in adaptive change. Resistance to snake venom alpha-neurotoxins (αNTXs) is a case of widespread convergence having evolved several times in snakes, lizards and mammals. Despite extreme toxicity of αNTXs, substitutions in its target, the nicotinic acetylcholine receptor (nAChR), prevent αNTX binding and render species resistant. Recently, the published meerkat (Herpestidae) genome revealed that meerkats have the same substitutions in nAChR as the venom-resistant Egyptian mongoose (Herpestidae), suggesting that venom-resistant nAChRs may be ancestral to Herpestids. Like the mongoose, many other species of feliform carnivores prey on venomous snakes, though their venom resistance has never been explored. To evaluate the prevalence and ancestry of αNTX resistance in mammals, we generate a dataset of mammalian nAChR using museum specimens and public datasets. We find five instances of convergent evolution within feliform carnivores, and an additional eight instances across all mammals sampled. Tests of selection show that these substitutions are evolving under positive selection. Repeated convergence suggests that this adaptation played an important role in the evolution of mammalian physiology and potentially venom evolution.