Dynamic genetic differentiation drives the widespread structural and functional convergent evolution of snake venom proteinaceous toxins

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.

The Clot Thickens: Differential Coagulotoxic and Cardiotoxic Activities of Anguimorpha Lizard Venoms

Despite their evolutionary novelty, lizard venoms are much less studied in comparison to the intense research on snake venoms. While the venoms of helodermatid lizards have long been assumed to be for defensive purposes, there is increasing evidence of toxic activities more useful for predation than defence (such as paralytic neurotoxicity). This study aimed to ascertain the effects of Heloderma, Lanthanotus, and Varanus lizard venoms on the coagulation and cardiovascular systems. Anticoagulant toxicity was demonstrated for the Varanus species studied, with the venoms prolonging clotting times in human and bird plasma due to the destructive cleavage of fibrinogen. In contrast, thromboelastographic analyses on human and bird plasmas in this study demonstrated a procoagulant bioactivity for Heloderma venoms. A previous study on Heloderma venom using factor-depleted plasmas as a proxy model suggested a procoagulant factor was present that activated either Factor XI or Factor XII, but could not ascertain the precise target. Our activation studies using purified zymogens confirmed FXII activation. Comparisons of neonate and adult H. exasperatum, revealed the neonates to be more potent in the ability to activate FXII, being more similar to the venom of the smaller species H. suspectum than the adult H. exasperatum. This suggests potent FXII activation a basal trait in the genus, present in the small bodied last common ancestor. This also indicates an ontogenetic difference in prey preferences in the larger Heloderma species paralleing the change in venom biochemistry. In addition, as birds lack Factor XII, the ability to clot avian plasma suggested an additional procoagulant site of action, which was revealed to be the activation of Factor VII, with H. horridum being the most potent. This study also examined the effects upon the cardiovascular system, including the liberation of kinins from kininogen, which contributes to hypotension induction. This form of toxicity was previously described for Heloderma venoms, and was revealed in this study was to also be a pathophysiological effect of Lanthanotus and Varanus venoms. This suggests that this toxic activity was present in the venom of the last common ancestor of the anguimorph lizards, which is consistent with kallikrein enzymes being a shared toxin trait. This study therefore uncovered novel actions of anguimorph lizard venoms, not only contributing to the evolutionary biology body of knowledge but also revealing novel activities to mine for drug design lead compounds.

Highly conserved and extremely variable: The paradoxical pattern of toxin expression revealed by comparative venom-gland transcriptomics of Phalotris (Serpentes: Dipsadidae)

Although non-front fanged snakes account for almost two-thirds of snake diversity, most studies on venom composition and evolution focus exclusively on front-fanged species, which comprise most of the clinically relevant accidents. Comprehensive reports on venom composition of non-front fanged snakes are still scarce for several groups. In this study, we address such shortage of knowledge by providing new insights about the venom composition among species of Phalotris, a poorly studied Neotropical dipsadid genus. Phalotris are known for their specialized venom delivery system and toxic venoms, which can cause life-threatening accidents in humans. We evaluate the venom-gland transcriptome of Phalotris, comparing the following three South American species: P. reticulatus for the Araucaria Pine forests, P. lemniscatus for the Pampa grasslands, and P. mertensi for the Brazilian Cerrado. Our results indicate similar venom profiles, in which they share a high expression level of Kunitz-type inhibitors (KUNZ). On the other hand, comparative analyses revealed substantial differences in the expression levels of C-type lectins (CTL) and snake venom metalloproteinases (SVMP). The diverse set of SVMP and CTL isoforms shows signals of positive selection, and we also identified truncated forms of type III SVMPs, which resemble type II and type I SVMPs of viperids. Additionally, we identified a CNP precursor hosting a proline-rich region containing a BPP motif resembling those commonly detected in viperid venoms with hypotensive activity. Altogether, our results suggest an evolutionary history favoring high expression levels of few KUNZ isoforms in Phalotris venoms, contrasting with a highly diverse set of SVMP and CTL isoforms. Such diversity can be comparable with the venom variability observed in some viperids. Our findings highlight the extreme phenotypic diversity of non-front fanged snakes and the importance to allocate greater effort to study neglected groups of Colubroidea.

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.

De Novo Assembly of Venom Gland Transcriptome of Tropidolaemus wagleri (Temple Pit Viper, Malaysia) and Insights into the Origin of Its Major Toxin, Waglerin

The venom proteome of Temple Pit Viper (Tropidolaemus wagleri) is unique among pit vipers, characterized by a high abundance of a neurotoxic peptide, waglerin. To further explore the genetic diversity of its toxins, the present study de novo assembled the venom gland transcriptome of T. wagleri from west Malaysia. Among the 15 toxin gene families discovered, gene annotation and expression analysis reveal the dominating trend of bradykinin-potentiating peptide/angiotensin-converting enzyme inhibitor-C-type natriuretic peptide (BPP/ACEI-CNP, 76.19% of all-toxin transcription) in the transcriptome, followed by P-III snake venom metalloproteases (13.91%) and other toxins. The transcript TwBNP01 of BPP/ACEI-CNP represents a large precursor gene (209 amino acid residues) containing the coding region for waglerin (24 residues). TwBNP01 shows substantial sequence variations from the corresponding genes of its sister species, Tropidolaemus subannulatus of northern Philippines, and other viperid species which diversely code for proline-rich small peptides such as bradykinin-potentiating peptides (BPPs). The waglerin/waglerin-like peptides, BPPs and azemiopsin are proline-rich, evolving de novo from multiple highly diverged propeptide regions within the orthologous BPP/ACEI-CNP genes. Neofunctionalization of the peptides results in phylogenetic constraints consistent with a phenotypic dichotomy, where Tropidolaemus spp. and Azemiops feae convergently evolve a neurotoxic trait while vasoactive BPPs evolve only in other species.

Evolution of Three-Finger Toxin Genes in Neotropical Colubrine Snakes (Colubridae)

Snake venom research has historically focused on front-fanged species (Viperidae and Elapidae), limiting our knowledge of venom evolution in rear-fanged snakes across their ecologically diverse phylogeny. Three-finger toxins (3FTxs) are a known neurotoxic component in the venoms of some rear-fanged snakes (Colubridae: Colubrinae), but it is unclear how prevalent 3FTxs are both in expression within venom glands and more broadly among colubrine species. Here, we used a transcriptomic approach to characterize the venom expression profiles of four species of colubrine snakes from the Neotropics that were dominated by 3FTx expression (in the genera Chironius, Oxybelis, Rhinobothryum, and Spilotes). By reconstructing the gene trees of 3FTxs, we found evidence of putative novel heterodimers in the sequences of Chironius multiventris and Oxybelis aeneus, revealing an instance of parallel evolution of this structural change in 3FTxs among rear-fanged colubrine snakes. We also found positive selection at sites within structural loops or "fingers" of 3FTxs, indicating these areas may be key binding sites that interact with prey target molecules. Overall, our results highlight the importance of exploring the venoms of understudied species in reconstructing the full evolutionary history of toxins across the tree of life.

Domain loss enabled evolution of novel functions in the snake three-finger toxin gene superfamily

Three-finger toxins (3FTXs) are a functionally diverse family of toxins, apparently unique to venoms of caenophidian snakes. Although the ancestral function of 3FTXs is antagonism of nicotinic acetylcholine receptors, redundancy conferred by the accumulation of duplicate genes has facilitated extensive neofunctionalization, such that derived members of the family interact with a range of targets. 3FTXs are members of the LY6/UPAR family, but their non-toxin ancestor remains unknown. Combining traditional phylogenetic approaches, manual synteny analysis, and machine learning techniques (including AlphaFold2 and ProtT5), we have reconstructed a detailed evolutionary history of 3FTXs. We identify their immediate ancestor as a non-secretory LY6, unique to squamate reptiles, and propose that changes in molecular ecology resulting from loss of a membrane-anchoring domain and changes in gene expression, paved the way for the evolution of one of the most important families of snake toxins.

Medically important snakes and snakebite envenoming in Iran

Snakebite is a common health condition in Iran with a diverse snake fauna, especially in tropical southern and mountainous western areas of the country with plethora of snake species. The list of medically important snakes, circumstances and effects of their bite, and necessary medical care require critical appraisal and should be updated regularly. This study aims to review and map the distributions of medically important snake species of Iran, re-evaluate their taxonomy, review their venomics, describe the clinical effects of envenoming, and discuss medical management and treatment, including the use of antivenom. Nearly 350 published articles and 26 textbooks with information on venomous and mildly venomous snake species and snakebites of Iran, were reviewed, many in Persian (Farsi) language, making them relatively inaccessible to an international readership. This has resulted in a revised updated list of Iran's medically important snake species, with taxonomic revisions of some, compilation of their morphological features, remapping of their geographical distributions, and description of species-specific clinical effects of envenoming. Moreover, the antivenom manufactured in Iran is discussed, together with treatment protocols that have been developed for the hospital management of envenomed patients.

Transcriptome mining extends the host range of the Flaviviridae to non-bilaterians

The flavivirids (family Flaviviridae) are a group of positive-sense RNA viruses that include well-documented agents of human disease. Despite their importance and ubiquity, the timescale of flavivirid evolution is uncertain. An ancient origin, spanning millions of years, is supported by their presence in both vertebrates and invertebrates and by the identification of a flavivirus-derived endogenous viral element in the peach blossom jellyfish genome (Craspedacusta sowerbii, phylum Cnidaria), implying that the flaviviruses arose early in the evolution of the Metazoa. To date, however, no exogenous flavivirid sequences have been identified in these hosts. To help resolve the antiquity of the Flaviviridae, we mined publicly available transcriptome data across the Metazoa. From this, we expanded the diversity within the family through the identification of 32 novel viral sequences and extended the host range of the pestiviruses to include amphibians, reptiles, and ray-finned fish. Through co-phylogenetic analysis we found cross-species transmission to be the predominate macroevolutionary event across the non-vectored flavivirid genera (median, 68 per cent), including a cross-species transmission event between bats and rodents, although long-term virus-host co-divergence was still a regular occurrence (median, 23 per cent). Notably, we discovered flavivirus-like sequences in basal metazoan species, including the first associated with Cnidaria. This sequence formed a basal lineage to the genus Flavivirus and was closer to arthropod and crustacean flaviviruses than those in the tamanavirus group, which includes a variety of invertebrate and vertebrate viruses. Combined, these data attest to an ancient origin of the flaviviruses, likely close to the emergence of the metazoans 750-800 million years ago.

Extreme Procoagulant Potency in Human Plasma of Venoms from the African Viperid Genera Atheris, Cerastes, and Proatheris and the Relative Efficacy of Antivenoms and Synthetic Enzyme-Inhibitors

The African viperid snake genera Atheris, Cerastes, and Proatheris are closely related, similar in size, but occupy extremely divergent ecological niches (arboreal in tropical rainforests, fossorial in deserts, and swamp-dwelling, respectively). Their venoms have not previously been subjected to comparative analyses for their action upon the coagulation of blood, most notably with significant data deficiencies from Atheris and Proatheris. In contrast, the closely related genus Echis is well-documented as capable of producing potent procoagulant effects. In light of this, we set out to compare the coagulotoxic actions of Atheris ceratophora, A. chlorechis, A. desaixi, A. nitschei, A. squamigera, C. cerastes, C. cerastes gasperettii, C. vipera, and Proatheris superciliaris and explore potential pharmacological interventions to reestablish normal blood coagulation. All venoms displayed extremely potent procoagulant effects, over twice as fast as the most potent Echis reported to date. Although Cerastes is used in the immunising mixture of two different regionally available antivenoms (Inoserp-MENA with C. cerastes, C. cerastes gasperettii, C. vipera and Saudi Arabian polyvalent with C. cerastes), none of the other species in this study are included in the immunising mixture of any antivenom. Notably, all the Cerastes species were only neutralised by the Inoserp-MENA antivenom. C. cerastes venom was not neutralised well by the Saudi Arabian antivenom, with the low levels of recognition for any of the Cerastes venoms suggesting a strong regional variation in the venom of this species, as the C. cerastes venom tested was of African (Tunisian) origin versus Saudi locality used in that antivenom's production. The other antivenoms (Micropharm EchiTAbG, ICP EchiTAb-Plus-ICP, Inosan Inoserp Pan-Africa, Premium Serums PANAF Sub-Sahara Africa, South African Vaccine Producers Echis, South African Vaccine Producers Polyvalent) all displayed trivial-to-no ability to neutralise the procoagulant toxicity of any of the Atheris, Cerastes, or Proatheris venoms. Comparative testing of the enzyme inhibitors DMPS, marimastat, and prinomastat, revealed a very potent neutralising capacity of marimastat, with prinomastat showing lower but still significant potency at the same molar concentration, while a 5× molar concentration of DMPS had no apparent effect on procoagulant venom effects normalized by the other inhibitors. These results and methods contribute to the body of knowledge of potential clinical effects and data necessary for evidence-based advancement of clinical management strategies.

Proteomic insight into the venom composition of the largest European rear-fanged snake, Malpolon monspessulanus monspessulanus

Snake envenomations constitute a worldwide neglected tropical disease, with the vast majority of lethal bites inflicted by front-fanged snakes from the viperid and elapid groups. Rear-fanged snakes (colubrids) were often considered harmless and as a result, are much less studied, but several documented deaths have suggested potent venom in this group as well. The largest European snake (Malpolon monspessulanus monspessulanus), known as the “Montpellier snake”, is such a rear-fanged snake that belongs to the Lamprophiidae family. Its venom remains largely unknown but cases of envenomation with neurological symptoms have been reported. Here, we provide the first insights into the composition of its venom using mass spectrometry methods. First, liquid chromatography coupled mass spectrometry analysis of the manually collected venom samples reveals a complex profile, with the majority of masses encompassing the range 500–3000 Da, 4000–8000 Da, and 10 000–30 000 Da. Next, shotgun proteomics allowed the identification of a total of 42 different known families of proteins, including snake venom metalloproteinases, peptidase M1, and cysteine-rich secretory proteins, as the most prominent. Interestingly, three-finger toxins were not detected, suggesting that neurotoxicity may occur via other, yet to be determined, toxin types. Overall, our results provide the basis for a better understanding of the effects of a peculiar snake venom on human symptomatology, but also on the main prey consumed by this species.

•

We investigate the venom composition of the largest European venomous snake.

•

LC-MS analysis of the crude manually collected venom revealed a complex profile.

•

Shotgun proteomic analysis identified 42 different protein families.

•

Major components include SVMPs, consistent with the clinical features.

Keel venom: Rhabdophis subminiatus (red-necked keelback) venom pathophysiologically affects diverse blood clotting pathways

Venoms are evolutionary novelties that have real-world implications due to their impact upon human health. However, relative to the abundant studies of elapid and viperid snake venoms, fewer investigations have been undertaken on those of rear-fanged snakes as they are more problematic for obtaining venom. While most rear-fanged venomous snakes are not considered to be of great medical importance, several species are capable of producing fatalities. Most notable among these are snakes from the genus Rhabdophis, the Asian "keelback" snakes. Prior work have described potent procoagulant toxicity suggesting Factor X and prothrombin activation, but did not investigate the ability to activate other clotting factors. Here we show that in addition to activating both Factor X and prothrombin (with prothrombin twice that of FX), the venom of Rhabdophis subminiatus is able to more potently activate Factor VII (ten times that of prothrombin), while also activating FXII and FIX equipotently to prothrombin, and with FXI also activated but at a much lower level. The ability to activate FVII represents a third convergent evolution of this trait. The Australian elapid clade of [Oxyuranus (taipans) + Pseudonaja (brown snakes)] was the first identified to have evolved this trait. and only recently was it shown to be independently present in another lineage (the Central American viperid species Porthidium volcanicum). In addition, the abilities to activate FXI and FXII are also convergent between R. subminiatus and P. volcanicum, but with R. subminiatus being much more potent. By testing across amphibian, avian, and mammalian plasmas we demonstrate that the venom is potently procoagulant across diverse plasma types. However, consistent with dietary preference, R. subminiatus venom was most potent upon amphibian plasma. While a Rhabdophis antivenom is produced in Japan to treat R. tigrinus envenomings, it is scarce even within Japan and is not exported. As this genus is very wide-ranging in Asia, alternate treatment options are in need of development. Hence we tested the ability of candidate, broad-spectrum enzyme inhibitors to neutralize R. subminiatus venom: marimastat was more effective than prinomastat but both marimastat and prinomastat were significantly more effective than DMPS (2,3-Dimercapto-1-propanesulfonic acid). The findings of this study shed light on the evolution of these fascinating rear-fanged snakes as well as explored their systemic effects upon blood coagulation and point to potential treatment options for the rare, but potentially lethal encounters.

The Target Selects the Toxin: Specific Amino Acids in Snake-Prey Nicotinic Acetylcholine Receptors That Are Selectively Bound by King Cobra Venoms

Snake venom is an adaptive ecological trait that has evolved primarily as a form of prey subjugation. Thus, the selection pressure for toxin diversification is exerted by the prey’s physiological targets, with this pressure being particularly acute for specialist feeders, such as the King Cobra species, all of which are snake-prey specialists. However, while extensive research has been undertaken to elucidate key amino acids that guide toxin structure–activity relationships, reciprocal investigations into the specific sites guiding prey-lineage selective effects have been lacking. This has largely been due to the lack of assay systems amenable to systematic amino acid replacements of targeted proteins in the prey’s physiological pathways. To fill this knowledge gap, we used a recently described approach based upon mimotope peptides corresponding to the orthosteric site of nicotinic acetylcholine receptor alpha-1 subunits, a major binding site for snake venom neurotoxins that cause flaccid paralysis. We investigated the venoms of four different types of King Cobra (Cambodian, Javan, Malaysian, and Thai). This approach allowed for the determination of the key amino acid positions in King Cobra snake prey that are selectively bound by the toxins, whereby replacing these amino acids in the snake-prey orthosteric site with those from lizards or rats resulted in a significantly lower level of binding by the venoms, while conversely replacing the lizard or rat amino acids with those from the snake at that position increased the binding. By doing such, we identified three negatively charged amino acids in the snake orthosteric site that are strongly bound by the positively charged neurotoxic three-finger toxins found in King Cobra venom. This study, thus, sheds light on the selection pressures exerted by a specialist prey item for the evolution of lineage-selective toxins.

Differential Antivenom and Small-Molecule Inhibition of Novel Coagulotoxic Variations in Atropoides, Cerrophidion, Metlapilcoatlus, and Porthidium American Viperid Snake Venoms

Within Neotropical pit-vipers, the Mexican/Central-American clade consisting of Atropoides, Cerrophidion, Metlapilcoatlus, and Porthidium is a wide-ranging, morphologically and ecologically diverse group of snakes. Despite their prevalence, little is known of the functional aspects of their venoms. This study aimed to fill the knowledge gap regarding coagulotoxic effects and to examine the potential of different therapeutic approaches. As a general trait, the venoms were shown to be anticoagulant but were underpinned by diverse biochemical actions. Pseudo-procoagulant activity (i.e., thrombin-like), characterized by the direct cleavage of fibrinogen to form weak fibrin clots, was evident for Atropoides picadoi, Cerrophidiontzotzilorum, Metlapilcoatlus mexicanus, M. nummifer, M. occiduus, M. olmec, and Porthidium porrasi. In contrast, other venoms cleaved fibrinogen in a destructive (non-clotting) manner, with C. godmani and C. wilsoni being the most potent. In addition to actions on fibrinogen, clotting enzymes were also inhibited. FXa was only weakly inhibited by most species, but Cerrophidion godmani and C. wilsoni were extremely strong in their inhibitory action. Other clotting enzymes were more widely inhibited by diverse species spanning the full taxonomical range, but in each case, there were species that had these traits notably amplified relatively to the others. C. godmani and C. wilsoni were the most potent amongst those that inhibited the formation of the prothrombinase complex and were also amongst the most potent inhibitors of Factor XIa. While most species displayed only low levels of thrombin inhibition, Porthidium dunni potently inhibited this clotting factor. The regional polyvalent antivenom produced by Instituto Picado Clodomiro was tested and was shown to be effective against the diverse anticoagulant pathophysiological effects. In contrast to the anticoagulant activities of the other species, Porthidium volcanicum was uniquely procoagulant through the activation of Factor VII and Factor XII. This viperid species is the first snake outside of the Oxyuranus/Pseudonaja elapid snake clade to be shown to activate FVII and the first snake venom of any kind to activate FXII. Interestingly, while small-molecule metalloprotease inhibitors prinomastat and marimastat demonstrated the ability to prevent the procoagulant toxicity of P. volcanicum, neither ICP antivenom nor inhibitor DMPS showed this effect. The extreme variation among the snakes here studied underscores how venom is a dynamic trait and how this can shape clinical outcomes and influence evolving treatment strategies.

Divergent Specialization of Simple Venom Gene Profiles among Rear-Fanged Snake Genera (Helicops and Leptodeira, Dipsadinae, Colubridae)

Many venomous animals express toxins that show extraordinary levels of variation both within and among species. In snakes, most studies of venom variation focus on front-fanged species in the families Viperidae and Elapidae, even though rear-fanged snakes in other families vary along the same ecological axes important to venom evolution. Here we characterized venom gland transcriptomes from 19 snakes across two dipsadine rear-fanged genera (Leptodeira and Helicops, Colubridae) and two front-fanged genera (Bothrops, Viperidae; Micrurus, Elapidae). We compared patterns of composition, variation, and diversity in venom transcripts within and among all four genera. Venom gland transcriptomes of rear-fanged Helicops and Leptodeira and front-fanged Micrurus are each dominated by expression of single toxin families (C-type lectins, snake venom metalloproteinase, and phospholipase A2, respectively), unlike highly diverse front-fanged Bothrops venoms. In addition, expression patterns of congeners are much more similar to each other than they are to species from other genera. These results illustrate the repeatability of simple venom profiles in rear-fanged snakes and the potential for relatively constrained venom composition within genera.