Snake venom gene expression is coordinated by novel regulatory architecture and the integration of multiple co-opted vertebrate pathways.

Understanding how regulatory mechanisms evolve is critical for understanding the processes that give rise to novel phenotypes. Snake venom systems represent a valuable and tractable model for testing hypotheses related to the evolution of novel regulatory networks, yet the regulatory mechanisms underlying venom production remain poorly understood. Here, we use functional genomics approaches to investigate venom regulatory architecture in the prairie rattlesnake and identify cis-regulatory sequences (enhancers and promoters), trans-regulatory transcription factors, and integrated signaling cascades involved in the regulation of snake venom genes. We find evidence that two conserved vertebrate pathways, the extracellular signal-regulated kinase and unfolded protein response pathways, were co-opted to regulate snake venom. In one large venom gene family (snake venom serine proteases), this co-option was likely facilitated by the activity of transposable elements. Patterns of snake venom gene enhancer conservation, in some cases spanning 50 million yr of lineage divergence, highlight early origins and subsequent lineage-specific adaptations that have accompanied the evolution of venom regulatory architecture. We also identify features of chromatin structure involved in venom regulation, including topologically associated domains and CTCF loops that underscore the potential importance of novel chromatin structure to coevolve when duplicated genes evolve new regulatory control. Our findings provide a model for understanding how novel regulatory systems may evolve through a combination of genomic processes, including tandem duplication of genes and regulatory sequences, cis-regulatory sequence seeding by transposable elements, and diverse transcriptional regulatory proteins controlled by a co-opted regulatory cascade.

Prevalent bee venom genes evolved before the aculeate stinger and eusociality.

BACKGROUND: Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected. RESULTS: Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification. CONCLUSIONS: Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.

Solenodon genome reveals convergent evolution of venom in eulipotyphlan mammals.

Venom systems are key adaptations that have evolved throughout the tree of life and typically facilitate predation or defense. Despite venoms being model systems for studying a variety of evolutionary and physiological processes, many taxonomic groups remain understudied, including venomous mammals. Within the order Eulipotyphla, multiple shrew species and solenodons have oral venom systems. Despite morphological variation of their delivery systems, it remains unclear whether venom represents the ancestral state in this group or is the result of multiple independent origins. We investigated the origin and evolution of venom in eulipotyphlans by characterizing the venom system of the endangered Hispaniolan solenodon (Solenodon paradoxus). We constructed a genome to underpin proteomic identifications of solenodon venom toxins, before undertaking evolutionary analyses of those constituents, and functional assessments of the secreted venom. Our findings show that solenodon venom consists of multiple paralogous kallikrein 1 (KLK1) serine proteases, which cause hypotensive effects in vivo, and seem likely to have evolved to facilitate vertebrate prey capture. Comparative analyses provide convincing evidence that the oral venom systems of solenodons and shrews have evolved convergently, with the 4 independent origins of venom in eulipotyphlans outnumbering all other venom origins in mammals. We find that KLK1s have been independently coopted into the venom of shrews and solenodons following their divergence during the late Cretaceous, suggesting that evolutionary constraints may be acting on these genes. Consequently, our findings represent a striking example of convergent molecular evolution and demonstrate that distinct structural backgrounds can yield equivalent functions.

Co-option of the same ancestral gene family gave rise to mammalian and reptilian toxins.

BACKGROUND: Evolution can occur with surprising predictability when organisms face similar ecological challenges. For most traits, it is difficult to ascertain whether this occurs due to constraints imposed by the number of possible phenotypic solutions or because of parallel responses by shared genetic and regulatory architecture. Exceptionally, oral venoms are a tractable model of trait evolution, being largely composed of proteinaceous toxins that have evolved in many tetrapods, ranging from reptiles to mammals. Given the diversity of venomous lineages, they are believed to have evolved convergently, even though biochemically similar toxins occur in all taxa. RESULTS: Here, we investigate whether ancestral genes harbouring similar biochemical activity may have primed venom evolution, focusing on the origins of kallikrein-like serine proteases that form the core of most vertebrate oral venoms. Using syntenic relationships between genes flanking known toxins, we traced the origin of kallikreins to a single locus containing one or more nearby paralogous kallikrein-like clusters. Additionally, phylogenetic analysis of vertebrate serine proteases revealed that kallikrein-like toxins in mammals and reptiles are genetically distinct from non-toxin ones. CONCLUSIONS: Given the shared regulatory and genetic machinery, these findings suggest that tetrapod venoms evolved by co-option of proteins that were likely already present in saliva. We term such genes 'toxipotent'-in the case of salivary kallikreins they already had potent vasodilatory activity that was weaponized by venomous lineages. Furthermore, the ubiquitous distribution of kallikreins across vertebrates suggests that the evolution of envenomation may be more common than previously recognized, blurring the line between venomous and non-venomous animals.

Differential evolution and neofunctionalization of snake venom metalloprotease domains.

Snake venom metalloproteases (SVMP) are composed of five domains: signal peptide, propeptide, metalloprotease, disintegrin, and cysteine-rich. Secreted toxins are typically combinatorial variations of the latter three domains. The SVMP-encoding genes of Psammophis mossambicus venom are unique in containing only the signal and propeptide domains. We show that the Psammophis SVMP propeptide evolves rapidly and is subject to a high degree of positive selection. Unlike Psammophis, some species of Echis express both the typical multidomain and the unusual monodomain (propeptide only) SVMP, with the result that a lower level of variation is exerted upon the latter. We showed that most mutations in the multidomain Echis SVMP occurred in the protease domain responsible for proteolytic and hemorrhagic activities. The cysteine-rich and disintegrin-like domains, which are putatively responsible for making the P-III SVMPs more potent than the P-I and P-II forms, accumulate the remaining variation. Thus, the binding sites on the molecule's surface are evolving rapidly whereas the core remains relatively conserved. Bioassays conducted on two post-translationally cleaved novel proline-rich peptides from the P. mossambicus propeptide domain showed them to have been neofunctionalized for specific inhibition of mammalian a7 neuronal nicotinic acetylcholine receptors. We show that the proline rich postsynaptic specific neurotoxic peptides from Azemiops feae are the result of convergent evolution within the precursor region of the C-type natriuretic peptide instead of the SVMP. The results of this study reinforce the value of studying obscure venoms for biodiscovery of novel investigational ligands.

Squeezers and leaf-cutters: differential diversification and degeneration of the venom system in toxicoferan reptiles.

Although it has been established that all toxicoferan squamates share a common venomous ancestor, it has remained unclear whether the maxillary and mandibular venom glands are evolving on separate gene expression trajectories or if they remain under shared genetic control. We show that identical transcripts are simultaneously expressed not only in the mandibular and maxillary glands, but also in the enigmatic snake rictal gland. Toxin molecular frameworks recovered in this study were three-finger toxin (3FTx), CRiSP, crotamine (beta-defensin), cobra venom factor, cystatin, epididymal secretory protein, kunitz, L-amino acid oxidase, lectin, renin aspartate protease, veficolin, and vespryn. We also discovered a novel low-molecular weight disulfide bridged peptide class in pythonid snake glands. In the iguanian lizards, the most highly expressed are potentially antimicrobial in nature (crotamine (beta-defensin) and cystatin), with crotamine (beta-defensin) also the most diverse. However, a number of proteins characterized from anguimorph lizards and caenophidian snakes with hemotoxic or neurotoxic activities were recruited in the common toxicoferan ancestor and remain expressed, albeit in low levels, even in the iguanian lizards. In contrast, the henophidian snakes express 3FTx and lectin toxins as the dominant transcripts. Even in the constricting pythonid and boid snakes, where the glands are predominantly mucous-secreting, low-levels of toxin transcripts can be detected. Venom thus appears to play little role in feeding behavior of most iguanian lizards or the powerful constricting snakes, and the low levels of expression argue against a defensive role. However, clearly the incipient or secondarily atrophied venom systems of these taxa may be a source of novel compounds useful in drug design and discovery.

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.

Structural and molecular diversification of the Anguimorpha lizard mandibular venom gland system in the arboreal species Abronia graminea.

In the past, toxinological research on reptiles has focused principally on clinically important species. As a result, our understanding of the evolution of the reptile venom system is limited. Here, for the first time, we describe the structural and molecular evolutionary features of the mandibular toxin-secreting gland of Abronia graminea, a representative of one of the poorly known and entirely arboreal lineages of anguimorph lizards. We show that the mandibular gland is robust and serous, characters consistent with those expected of a toxin-secreting gland in active use. A wide array of transcripts were recovered that were homologous to those encoded by the indisputably venomous helodermatid lizards. We show that some of these toxin transcripts are evolving under active selection and show evidence of rapid diversification. Helokinestatin peptides in particular are revealed to have accumulated residues that have undergone episodic diversifying selections. Conversely, the natriuretic peptides have evolved under tremendous evolutionary constraints despite being encoded in tandem with helokinestatins by the same gene precursor. Of particular note is the sequencing for the first time of kunitz peptides from a lizard toxin-secreting gland. Not only are kunitz peptides shown to be an ancestral toxicoferan toxin, the ancestral state of this peptide is revealed to be a dual domain encoding precursor. This research provides insight into the evolutionary history of the ancient toxicoferan reptile venom system. In addition, it shows that even 'clinically irrelevant' species can be a rich source of novel venom components, worthy of investigation for drug design and biomedical research.

Venom profile of the European carpenter bee Xylocopa violacea: Evolutionary and applied considerations on its toxin components.

Modern venomics is increasing its focus on hymenopterans such as honeybees, bumblebees, parasitoid wasps, ants and true wasps. However solitary bees remain understudied in comparison and the few available venom studies focus on short melittin-like sequences and antimicrobial peptides. Herein we describe the first comprehensive venom profile of a solitary bee, the violet carpenter bee Xylocopa violacea, by using proteo-transcriptomics. We reveal a diverse and complex venom profile with 43 different protein families identified from dissected venom gland extracts of which 32 are also detected in the defensively injected venom. Melittin and apamin are the most highly secreted components, followed by Phospholipase A2, Icarapin, Secapin and three novel components. Other components, including eight novel protein families, are rather lowly expressed. We further identify multiple forms of apamin-like peptides. The melittin-like sequences of solitary bees separate into two clades, one comprised most sequences from solitary bees including xylopin (the variant in Xylocopa), while sequences from Lasioglossa appear closer related to melittin-like peptides from Bombus (Bombolittins). Our study suggests that more proteo-transcriptomic data from other solitary bees should be complemented with corresponding genome data to fully understand the evolution and complexity of bee venom proteins, and is of a particular need to disentangle the ambiguous phylogenetic relations of short peptides.

Modern venomics-Current insights, novel methods, and future perspectives in biological and applied animal venom research.

Venoms have evolved >100 times in all major animal groups, and their components, known as toxins, have been fine-tuned over millions of years into highly effective biochemical weapons. There are many outstanding questions on the evolution of toxin arsenals, such as how venom genes originate, how venom contributes to the fitness of venomous species, and which modifications at the genomic, transcriptomic, and protein level drive their evolution. These questions have received particularly little attention outside of snakes, cone snails, spiders, and scorpions. Venom compounds have further become a source of inspiration for translational research using their diverse bioactivities for various applications. We highlight here recent advances and new strategies in modern venomics and discuss how recent technological innovations and multi-omic methods dramatically improve research on venomous animals. The study of genomes and their modifications through CRISPR and knockdown technologies will increase our understanding of how toxins evolve and which functions they have in the different ontogenetic stages during the development of venomous animals. Mass spectrometry imaging combined with spatial transcriptomics, in situ hybridization techniques, and modern computer tomography gives us further insights into the spatial distribution of toxins in the venom system and the function of the venom apparatus. All these evolutionary and biological insights contribute to more efficiently identify venom compounds, which can then be synthesized or produced in adapted expression systems to test their bioactivity. Finally, we critically discuss recent agrochemical, pharmaceutical, therapeutic, and diagnostic (so-called translational) aspects of venoms from which humans benefit.

Honeybee Cognition as a Tool for Scientific Engagement.

Apis mellifera (honeybees) are a well-established model for the study of learning and cognition. A robust conditioning protocol, the olfactory conditioning of the proboscis extension response (PER), provides a powerful but straightforward method to examine the impact of varying stimuli on learning performance. Herein, we provide a protocol that leverages PER for classroom-based community or student engagement. Specifically, we detail how a class of high school students, as part of the Ryukyu Girls Outreach Program, examined the effects of caffeine and dopamine on learning performance in honeybees. Using a modified version of the PER conditioning protocol, they demonstrated that caffeine, but not dopamine, significantly reduced the number of trials required for a successful conditioning response. In addition to providing an engaging and educational scientific activity, it could be employed, with careful oversight, to garner considerable reliable data examining the effects of varying stimuli on honeybee learning.

How the Toxin got its Toxicity.

Venom systems are functional and ecological traits, typically used by one organism to subdue or deter another. A predominant subset of their constituent molecules-"toxins"-share this ecological function and are therefore molecules that mediate interactions between organisms. Such molecules have been referred to as "exochemicals." There has been debate within the field of toxinology concerning the evolutionary pathways leading to the "recruitment" of a gene product for a toxic role within venom. We review these discussions and the evidence interpreted in support of alternate pathways, along with many of the most popular models describing the origin of novel molecular functions in general. We note that such functions may arise with or without gene duplication occurring and are often the consequence of a gene product encountering a novel "environment," i.e., a range of novel partners for molecular interaction. After stressing the distinction between "activity" and "function," we describe in detail the results of a recent study which reconstructed the evolutionary history of a multigene family that has been recruited as a toxin and argue that these results indicate that a pluralistic approach to understanding the origin of novel functions is advantageous. This leads us to recommend that an expansive approach be taken to the definition of "neofunctionalization"-simply the origins of a novel molecular function by any process-and "recruitment"-the "weaponization" of a molecule via the acquisition of a toxic function in venom, by any process. Recruitment does not occur at the molecular level or even at the level of gene expression, but only when a confluence of factors results in the ecological deployment of a physiologically active molecule as a toxin. Subsequent to recruitment, the evolutionary regime of a gene family may shift into a more dynamic form of "birth-and-death." Thus, recruitment leads to a form of "downwards causation," in which a change at the ecological level at which whole organisms interact leads to a change in patterns of evolution at the genomic level.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins.

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the 'venom-ome' and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 'venom-ome-specific toxins' (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

Snake venom NAD glycohydrolases: primary structures, genomic location, and gene structure.

NAD glycohydrolase (EC 3.2.2.5) (NADase) sequences have been identified in 10 elapid and crotalid venom gland transcriptomes, eight of which are complete. These sequences show very high homology, but elapid and crotalid sequences also display consistent differences. As in Aplysia kurodai ADP-ribosyl cyclase and vertebrate CD38 genes, snake venom NADase genes comprise eight exons; however, in the Protobothrops mucrosquamatus genome, the sixth exon is sometimes not transcribed, yielding a shortened NADase mRNA that encodes all six disulfide bonds, but an active site that lacks the catalytic glutamate residue. The function of this shortened protein, if expressed, is unknown. While many vertebrate CD38s are multifunctional, liberating both ADP-ribose and small quantities of cyclic ADP-ribose (cADPR), snake venom CD38 homologs are dedicated NADases. They possess the invariant TLEDTL sequence (residues 144-149) that bounds the active site and the catalytic residue, Glu228. In addition, they possess a disulfide bond (Cys121-Cys202) that specifically prevents ADP-ribosyl cyclase activity in combination with Ile224, in lieu of phenylalanine, which is requisite for ADPR cyclases. In concert with venom phosphodiesterase and 5'-nucleotidase and their ecto-enzyme homologs in prey tissues, snake venom NADases comprise part of an envenomation strategy to liberate purine nucleosides, and particularly adenosine, in the prey, promoting prey immobilization via hypotension and paralysis.

Proteomic and functional variation within black snake venoms (Elapidae: Pseudechis).

Pseudechis (black snakes) is an Australasian elapid snake genus that inhabits much of mainland Australia, with two representatives confined to Papua New Guinea. The present study is the first to analyse the venom of all 9 described Pseudechis species (plus one undescribed species) to investigate the evolution of venom composition and functional activity. Proteomic results demonstrated that the typical Pseudechis venom profile is dominated by phospholipase A2 toxins. Strong cytotoxicity was the dominant function for most species. P. porphyriacus, the most basal member of the genus, also exhibited the most divergent venom composition, being the only species with appreciable amounts of procoagulant toxins. The relatively high presence of factor Xa recovered in P. porphyriacus venom may be related to a predominantly amphibian diet. Results of this study provide important insights to guide future ecological and toxinological investigations.

The Cardiovascular and Neurotoxic Effects of the Venoms of Six Bony and Cartilaginous Fish Species.

Fish venoms are often poorly studied, in part due to the difficulty in obtaining, extracting, and storing them. In this study, we characterize the cardiovascular and neurotoxic effects of the venoms from the following six species of fish: the cartilaginous stingrays Neotrygon kuhlii and Himantura toshi, and the bony fish Platycephalus fucus, Girella tricuspidata, Mugil cephalus, and Dentex tumifrons. All venoms (10-100 µg/kg, i.v.), except G. tricuspidata and P. fuscus, induced a biphasic response on mean arterial pressure (MAP) in the anesthetised rat. P. fucus venom exhibited a hypotensive response, while venom from G. tricuspidata displayed a single depressor response. All venoms induced cardiovascular collapse at 200 µg/kg, i.v. The in vitro neurotoxic effects of venom were examined using the chick biventer cervicis nerve-muscle (CBCNM) preparation. N. kuhlii, H. toshi, and P. fucus venoms caused concentration-dependent inhibition of indirect twitches in the CBCNM preparation. These three venoms also inhibited responses to exogenous acetylcholine (ACh) and carbachol (CCh), but not potassium chloride (KCl), indicating a post-synaptic mode of action. Venom from G. tricuspidata, M. cephalus, and D. tumifrons had no significant effect on indirect twitches or agonist responses in the CBCNM. Our results demonstrate that envenoming by these species of fish may result in moderate cardiovascular and/or neurotoxic effects. Future studies aimed at identifying the molecules responsible for these effects could uncover potentially novel lead compounds for future pharmaceuticals, in addition to generating new knowledge about the evolutionary relationships between venomous animals.

Correlation between ontogenetic dietary shifts and venom variation in Australian brown snakes (Pseudonaja).

Venom is a key evolutionary trait, as evidenced by its widespread convergent evolution across the animal kingdom. In an escalating prey-predator arms race, venoms evolve rapidly to guarantee predatory or defensive success. Variation in venom composition is ubiquitous among snakes. Here, we tested variation in venom activity on substrates relevant to blood coagulation among Pseudonaja (brown snake) species, Australian elapids responsible for the majority of medically important human envenomations in Australia. A functional approach was employed to elucidate interspecific variation in venom activity in all nine currently recognised species of Pseudonaja. Fluorometric enzymatic activity assays were performed to test variation in whole venom procoagulant activity among species. Analyses confirmed the previously documented ontogenetic shift from non-coagulopathic venom in juveniles to coagulopathic venom as adults, except for the case of P. modesta, which retains non-coagulopathic venom as an adult. These shifts in venom activity correlate with documented ontogenetic shifts in diet among brown snakes from specialisation on reptilian prey as juveniles (and throughout the life cycle of P. modesta), to a more generalised diet in adults that includes mammals. The results of this study bring to light findings relevant to both clinical and evolutionary toxinology.

Enter the Dragon: The Dynamic and Multifunctional Evolution of Anguimorpha Lizard Venoms.

While snake venoms have been the subject of intense study, comparatively little work has been done on lizard venoms. In this study, we have examined the structural and functional diversification of anguimorph lizard venoms and associated toxins, and related these results to dentition and predatory ecology. Venom composition was shown to be highly variable across the 20 species of Heloderma, Lanthanotus, and Varanus included in our study. While kallikrein enzymes were ubiquitous, they were also a particularly multifunctional toxin type, with differential activities on enzyme substrates and also ability to degrade alpha or beta chains of fibrinogen that reflects structural variability. Examination of other toxin types also revealed similar variability in their presence and activity levels. The high level of venom chemistry variation in varanid lizards compared to that of helodermatid lizards suggests that venom may be subject to different selection pressures in these two families. These results not only contribute to our understanding of venom evolution but also reveal anguimorph lizard venoms to be rich sources of novel bioactive molecules with potential as drug design and development lead compounds.

The Evolution of Fangs, Venom, and Mimicry Systems in Blenny Fishes.

Venom systems have evolved on multiple occasions across the animal kingdom, and they can act as key adaptations to protect animals from predators [1]. Consequently, venomous animals serve as models for a rich source of mimicry types, as non-venomous species benefit from reductions in predation risk by mimicking the coloration, body shape, and/or movement of toxic counterparts [2-5]. The frequent evolution of such deceitful imitations provides notable examples of phenotypic convergence and are often invoked as classic exemplars of evolution by natural selection. Here, we investigate the evolution of fangs, venom, and mimetic relationships in reef fishes from the tribe Nemophini (fangblennies). Comparative morphological analyses reveal that enlarged canine teeth (fangs) originated at the base of the Nemophini radiation and have enabled a micropredatory feeding strategy in non-venomous Plagiotremus spp. Subsequently, the evolution of deep anterior grooves and their coupling to venom secretory tissue provide Meiacanthus spp. with toxic venom that they effectively employ for defense. We find that fangblenny venom contains a number of toxic components that have been independently recruited into other animal venoms, some of which cause toxicity via interactions with opioid receptors, and result in a multifunctional biochemical phenotype that exerts potent hypotensive effects. The evolution of fangblenny venom has seemingly led to phenotypic convergence via the formation of a diverse array of mimetic relationships that provide protective (Batesian mimicry) and predatory (aggressive mimicry) benefits to other fishes [2, 6]. Our results further our understanding of how novel morphological and biochemical adaptations stimulate ecological interactions in the natural world.

How the Cobra Got Its Flesh-Eating Venom: Cytotoxicity as a Defensive Innovation and Its Co-Evolution with Hooding, Aposematic Marking, and Spitting.

The cytotoxicity of the venom of 25 species of Old World elapid snake was tested and compared with the morphological and behavioural adaptations of hooding and spitting. We determined that, contrary to previous assumptions, the venoms of spitting species are not consistently more cytotoxic than those of closely related non-spitting species. While this correlation between spitting and non-spitting was found among African cobras, it was not present among Asian cobras. On the other hand, a consistent positive correlation was observed between cytotoxicity and utilisation of the defensive hooding display that cobras are famous for. Hooding and spitting are widely regarded as defensive adaptations, but it has hitherto been uncertain whether cytotoxicity serves a defensive purpose or is somehow useful in prey subjugation. The results of this study suggest that cytotoxicity evolved primarily as a defensive innovation and that it has co-evolved twice alongside hooding behavior: once in the Hemachatus + Naja and again independently in the king cobras (Ophiophagus). There was a significant increase of cytotoxicity in the Asian Naja linked to the evolution of bold aposematic hood markings, reinforcing the link between hooding and the evolution of defensive cytotoxic venoms. In parallel, lineages with increased cytotoxicity but lacking bold hood patterns evolved aposematic markers in the form of high contrast body banding. The results also indicate that, secondary to the evolution of venom rich in cytotoxins, spitting has evolved three times independently: once within the African Naja, once within the Asian Naja, and once in the Hemachatus genus. The evolution of cytotoxic venom thus appears to facilitate the evolution of defensive spitting behaviour. In contrast, a secondary loss of cytotoxicity and reduction of the hood occurred in the water cobra Naja annulata, which possesses streamlined neurotoxic venom similar to that of other aquatic elapid snakes (e.g., hydrophiine sea snakes). The results of this study make an important contribution to our growing understanding of the selection pressures shaping the evolution of snake venom and its constituent toxins. The data also aid in elucidating the relationship between these selection pressures and the medical impact of human snakebite in the developing world, as cytotoxic cobras cause considerable morbidity including loss-of-function injuries that result in economic and social burdens in the tropics of Asia and sub-Saharan Africa.