Convergent evolution of pain-inducing defensive venom components in spitting cobras

Ontogeny, not prey availability, underlies allopatric venom variability in insular and mainland populations of Vipera ammodytes

Allopatric populations living under distinct ecological conditions are excellent systems to infer factors underlying intraspecific venom variation. The venom composition of two populations of Vipera ammodytes, insular with a diet based on ectotherms and mainland with a diet based on ectotherms and endotherms, was compared considering the sex and age of individuals. Ten toxin families, dominated by PLA2, svMP, svSP, and DI, were identified through a bottom-up approach. The venom profiles of adult females and males were similar. Results from 58 individual SDS-PAGE profiles and venom pool analysis revealed significant differences between juveniles compared to subadults and adults. Two venom phenotypes were identified: a juvenile svMP-dominated and KUN-lacking phenotype and an adult PLA2/svMP-balanced and KUN-containing phenotype. Despite differences in prey availability (and, therefore, diet) between populations, no significant differences in venom composition were found. As the populations are geographically isolated, the lack of venom diversification could be explained by insufficient time for natural selection and/or genetic drift to act on the venom composition of island vipers. However, substantial differences in proteomes were observed when compared to venoms from geographically distant populations inhabiting different conditions. These findings highlight the need to consider ecological and evolutionary processes when studying venom variability. SIGNIFICANCE: This study provides the first comprehensive analysis of the venom composition of two allopatric populations of Vipera ammodytes, living under similar abiotic (climate) but distinct biotic (prey availability) conditions. The ontogenetic changes in venom composition, coupled with the lack of differences between sex and between populations, shed light on the main determinants of venom evolution in this medically important snake. Seven new proteomes may facilitate future comparative studies of snake venom evolution. This study highlights the importance of considering ecological and evolutionary factors to understand snake venom variation.

Genetic variability in snake venom and its implications for antivenom development in sub-Saharan Africa

Snake venom, a complex mixture of proteins, has attracted human attention for centuries due to its associated mortality, morbidity and other therapeutic properties. In sub-Saharan Africa (SSA), where snakebites pose a significant health risk, understanding the genetic variability of snake venoms is crucial for developing effective antivenoms. The wide geographic distribution of venomous snake species in SSA countries demonstrates the need to develop specific and broad antivenoms. However, the development of broad antivenoms has been hindered by different factors, such as antivenom cross-reactivity and polygenic paratopes. While specific antivenoms have been hindered by the numerous snake species across the SSA region, current antivenoms, such as SAIMR polyvalent and Premium Serums & Vaccines, exhibit varying degrees of cross-reactivity. Such ability to cross-react enables the antivenoms to target multiple components from the different snake species. The advent of biotechnological innovations, including recombinant antibodies, small-molecule drugs, monoclonal antibodies and synthetic antivenoms, presents options for eliminating limitations associated with traditional plasma-derived antivenoms. However, challenges still persist, especially in SSA, in addressing genetic variability, as evidenced by inadequate testing capacity and limited genomic research facilities. This comprehensive review explores the genetic variability of snake venoms in SSA, emphasizing the venom composition of various snake species and their interactions. This information is critical in developing multiple strategies during antivenom development. Finally, it offers information concerning the need for extensive collaborative engagements, technological advancements and comprehensive genomic evaluations to produce targeted and effective antivenoms.

Insight Into Body Size Evolution in Aves: Based on Some Body Size-Related Genes

Birds exhibit remarkable variations in body size, making them an ideal group for the study of adaptive evolution. However, the genetic mechanisms underlying body size evolution in avian species remain inadequately understood. This study investigates the evolutionary patterns of avian body size by analyzing 15 body-size-related genes, including GHSR, IGF2BP1, and IGFBP7 from the growth hormone/insulin-like growth factor axis, EIF2AK3, GALNS, NCAPG, PLOD1, and PLAG1 associated with tall stature, and ACAN, OBSL1, and GRB10 associated with short stature, four genes previously reported in avian species: ATP11A, PLXDC2, TNS3, and TUBGCP3. The results indicate significant adaptive evolution of body size-related genes across different avian lineages. Notably, in the IGF2BP1 gene, a significant positive correlation was observed between the evolutionary rate and body size, suggesting that larger bird species exhibit higher evolutionary rates of the IGF2BP1 gene. Furthermore, the IGFBP7 and PLXDC2 genes demonstrated accelerated evolution in large- and medium-sized birds, respectively, indicating distinct evolutionary patterns for these genes among birds of different sizes. The branch-site model analysis identified numerous positively selected sites, primarily concentrated near functional domains, thereby reinforcing the critical role of these genes in body size evolution. Interestingly, extensive convergent evolution was detected in lineages with larger body sizes. This study elucidates the genetic basis of avian body size evolution for the first time, identifying adaptive evolutionary patterns of body size-related genes across birds of varying sizes and documenting patterns of convergent evolution. These findings provide essential genetic data and novel insights into the adaptive evolution of body size in birds.

Hooding cobras can get ahead of other snakes in the ability to evoke human fear

Fear of snakes is common not only in humans but also in other primates. Consequently, snakes are salient stimuli associated with prioritized attention, early detection and emotional significance. This has been interpreted as an adaptive evolutionary response of the primate brain to a risk of envenoming by a hidden snake. However, the struggle between mammals and snakes is not one-sided. Humans and carnivores regularly kill snakes, and thus snakes develop deterring defensive behaviour that may directly evoke enhanced fear. Here, we show that snakes depicted in threatening posture evoked on average more fear than those in resting posture. Significantly, African (Somali) and European (Czech) respondents considerably agreed on the relative fear elicited by various snakes. Nonetheless, not all defensive postures are equally efficient. Threatening cobras were perceived as top fear-evoking stimuli, even though most of them are not considered very frightening in resting posture. This effect can be attributed to their conspicuous hooding posture which evolved into an efficient warning signal for mammalian predators. Our result demonstrates that cobras are more effective than other snakes in the ability to evoke human fear by a simple behavioural display-hooding. This can be primarily explained by the behavioural evolution of cobras which successfully exploited pre-existing cognitive mechanisms of mammals. Whether human ancestors cohabiting with deadly venomous cobras further improved their fear response to hooding is uncertain, but likely.

Age Is Just a Number: Ontogenetic Conservation in Activation of Blood Clotting Factors VII, X, and XII by Caucasus Blunt-Nosed Viper (Macrovipera lebetina obtusa) Venoms

This study examined the pathophysiological effects of venoms from neonate and adult specimens of the viperid snake Macrovipera lebetina obtusa, focusing on their ability to activate various blood clotting factors in human plasma. All venoms exhibited strong procoagulant properties. In concentration-response tests, the clotting potency of the neonate venoms fell within the range of their parents' maximum clotting velocities and areas under the curve. Intriguingly, females were more potent than males within each age group, but this requires a larger sample size to confirm. Antivenom neutralization efficacy was equipotent across age groups. The venoms potently activated Factor X (FX) robustly, consistent with previous knowledge of this genus. For the first time, the ability to activate Factors VII (FVII) and XII (FXII) was identified in this genus, with FXII exhibiting particularly strong activation. The study found no significant ontogenetic variation in procoagulant venom potency on human plasma, convergent with the Daboia genus, the other large-bodied lineage within the Palearctic viperid clade. However, the activation of FXII and FVII reveals previously undocumented pathways in the procoagulant activity of these venoms, contributing to the broader understanding of venom evolution and its clinical impacts. These findings have implications for venom biodiscovery and the development of antivenoms, highlighting the complexity of clotting factor activation beyond traditional investigations that have myopically focused upon FX and prothrombin pathways, thereby underscoring the importance of exploring additional clotting factors.

Venom variation among the three subspecies of the North African mountain viper Vipera monticola Saint Girons 1953

The North African mountain viper (Vipera monticola) is a medically relevant venomous snake distributed in Morocco, Algeria, and Tunisia. Three subspecies of V. monticola, exhibiting differences in morphotypes and dietary regimes, are currently recognised: V. m. monticola, V. m. atlantica, and V. m. saintgironsi. Through the application of snake venomics, we analysed the venoms of specimens of Moroccan origin belonging to each of the three subspecies. Snake venom metalloproteinases (svMP), snake venom serine proteases (svSP), C-type lectin and C-type lectin-related proteins (CTL), and phospholipases A2 (PLA2) were predominant, with PLA2 being the most abundant toxin family overall. Disintegrins (DI) and cysteine-rich secretory proteins (CRISP) were exclusive to V. m. monticola and V. m. atlantica, while l-amino-acid oxidases (LAAO) were only found in V. m. saintgironsi. The differences detected in the venom profiles, as well as in presence/absence and relative abundances of toxin families, indicate the occurrence of intraspecific venom variation within V. monticola. The identified patterns of venom similarity between subspecies seem to align more with their phylogenetic relationships than with the reported differences in their feeding habits.

Adaptive trade-offs between vertebrate defence and insect predation drive Amazonian ant venom evolution

Stinging ants have diversified into various ecological niches, and selective pressures may have contributed to shape the composition of their venom. To explore the drivers underlying venom variation in ants, we sampled 15 South American rainforest species and recorded a range of traits, including ecology, morphology and venom bioactivities. Principal component analysis of both morphological and venom bioactivity traits reveals that stinging ants display two functional strategies where species have evolved towards either an exclusively offensive venom or a multi-functional venom. Additionally, phylogenetic comparative analysis indicates that venom function (predatory, defensive or both) and mandible morphology correlate with venom bioactivity and volume. Further analysis of the venom biochemistry of the 15 species revealed switches between cytotoxic and neurotoxic venom compositions among species. Our study supports an evolutionary trade-off between the ability of venom to deter vertebrate predators and to paralyse insect prey which are correlated with different venom compositions and life-history strategies among Formicidae.

The Toxin Diversity, Cytotoxicity, and Enzymatic Activity of Cape Cobra (Naja nivea) Venom

"True" cobras (genus Naja) are among the venomous snakes most frequently involved in snakebite accidents in Africa and Asia. The Cape cobra (Naja nivea) is one of the African cobras of highest medical importance, but much remains to be learned about its venom. Here, we used a shotgun proteomics approach to better understand the qualitative composition of N. nivea venom and tested its cytotoxicity and protease activity as well as its effect on intracellular Ca2+ release and NO synthesis. We identified 156 venom components representing 17 protein families, with the dominant ones being three-finger toxins, mostly of the short-chain type. Two-thirds of the three-finger toxin entries identified were assigned as cytotoxins, while the remainder were categorized as neurotoxins, including short-chain, long-chain, and ancestral three-finger toxins. We also identified snake venom metalloproteinases and members of CRISP, l-amino acid oxidase, and other families. Protease activity and its effect on intracellular Ca2+ release and NO synthesis were low. Phospholipase A2 activity was surprisingly high, despite this toxin family being marginally recovered in the analyzed venom. Cytotoxicity was relevant only at higher venom concentrations, with macrophage and neuroblastoma cell lines showing the lowest viability. These results are in line with the predominantly neurotoxic envenomation symptoms caused by Cape cobra bites. The present overview of the qualitatively complex and functionally intriguing venom of N. nivea may provide insights into the pathobiochemistry of this species' venom.

Exploring Top-Down Mass Spectrometric Approaches To Probe Forest Cobra (Naja melanoleuca) Venom Proteoforms

Snake venoms are comprised of bioactive proteins and peptides that facilitate severe snakebite envenomation symptoms. A comprehensive understanding of venom compositions and the subtle heterogeneity therein is important. While bottom-up proteomics has been the well-established approach to catalogue venom compositions, top-down proteomics has emerged as a complementary strategy to characterize venom heterogeneity at the intact protein level. However, top-down proteomics has not been as widely implemented in the snake venom field as bottom-up proteomics, with various emerging top-down methods yet to be developed for venom systems. Here, we have explored three main top-down mass spectrometry methodologies in a proof-of-concept study to characterize selected three-finger toxin and phospholipase A2 proteoforms from the forest cobra (Naja melanoleuca) venom. We demonstrated the utility of a data-independent acquisition mode "MSE" for untargeted fragmentation on a chromatographic time scale and its improvement in protein sequence coverage compared to conventional targeted tandem mass spectrometry analysis. We also showed that protein identification can be further improved using a hybrid fragmentation approach, combining electron-capture dissociation and collision-induced dissociation. Lastly, we reported the promising application of multifunctional cyclic ion mobility separation and post-ion mobility fragmentation on snake venom proteins for the first time.

Intraspecific venom variation in the medically important puff adder (Bitis arietans): Comparative venom gland transcriptomics, in vitro venom activity and immunological recognition by antivenom

BACKGROUND

Variation in snake venoms is well documented, both between and within species, with intraspecific venom variation often correlated with geographically distinct populations. The puff adder, Bitis arietans, is widely distributed across sub-Saharan Africa and into the Arabian Peninsula where it is considered a leading cause of the ~310,000 annual snakebites across the region, with its venom capable of causing substantial morbidity and mortality. Despite its medical importance and wide geographic distribution, there is little known about venom variation between different B. arietans populations and the potential implications of this variation on antivenom efficacy.

METHODOLOGY

We applied a range of analyses, including venom gland transcriptomics, in vitro enzymatic assays and reverse phase chromatography to comparatively analyse B. arietans venoms originating from Nigeria, Tanzania, and South Africa. Immunological assays and in vitro enzymatic neutralisation assays were then applied to investigate the impact of venom variation on the potential efficacy of three antivenom products; SAIMR Polyvalent, EchiTAb-Plus and Fav-Afrique.

FINDINGS

Through the first comparison of venom gland transcriptomes of B. arietans from three geographically distinct regions (Nigeria, Tanzania, and South Africa), we identified substantial variation in toxin expression. Findings of venom variation were further supported by chromatographic venom profiling, and the application of enzymatic assays to quantify the activity of three pathologically relevant toxin families. However, the use of western blotting, ELISA, and in vitro enzymatic inhibition assays revealed that variation within B. arietans venom does not appear to substantially impact upon the efficacy of three African polyvalent antivenoms.

CONCLUSIONS

The large distribution and medical importance of B. arietans makes this species ideal for understanding venom variation and the impact this has on therapeutic efficacy. The findings in this study highlight the likelihood for considerable venom toxin variation across the range of B. arietans, but that this may not dramatically impact upon the utility of treatment available in the region.

Quantification of snake venom proteomes by mass spectrometry-considerations and perspectives

The advent of soft ionization mass spectrometry-based proteomics in the 1990s led to the development of a new dimension in biology that conceptually allows for the integral analysis of whole proteomes. This transition from a reductionist to a global-integrative approach is conditioned to the capability of proteomic platforms to generate and analyze complete qualitative and quantitative proteomics data. Paradoxically, the underlying analytical technique, molecular mass spectrometry, is inherently nonquantitative. The turn of the century witnessed the development of analytical strategies to endow proteomics with the ability to quantify proteomes of model organisms in the sense of "an organism for which comprehensive molecular (genomic and/or transcriptomic) resources are available." This essay presents an overview of the strategies and the lights and shadows of the most popular quantification methods highlighting the common misuse of label-free approaches developed for model species' when applied to quantify the individual components of proteomes of nonmodel species (In this essay we use the term "non-model" organisms for species lacking comprehensive molecular (genomic and/or transcriptomic) resources, a circumstance that, as we detail in this review-essay, conditions the quantification of their proteomes.). We also point out the opportunity of combining elemental and molecular mass spectrometry systems into a hybrid instrumental configuration for the parallel identification and absolute quantification of venom proteomes. The successful application of this novel mass spectrometry configuration in snake venomics represents a proof-of-concept for a broader and more routine application of hybrid elemental/molecular mass spectrometry setups in other areas of the proteomics field, such as phosphoproteomics, metallomics, and in general in any biological process where a heteroatom (i.e., any atom other than C, H, O, N) forms integral part of its mechanism.

Comparative venom analysis between melanistic and normally coloured phenotypes of the common adder (Vipera berus)

Snake venom is an ecologically relevant functional trait directly linked with a snake's fitness and survival, facilitating predation and defence. Snake venom variation occurs at all taxonomic levels, but the study at the intraspecific level is still in its early stages. The common adder (Vipera berus) exhibits considerable variation in colour phenotypes across its distribution range. Melanistic (fully black) individuals are the subject of myths and fairytales, and in German folklore such 'hell adders' are considered more toxic than their normally coloured conspecifics despite any formal investigation. Here, we provide the first comparative analysis of venoms from melanistic and normally coloured common adders. Specifically, we compared the venom profiles by sodium dodecylsulfate polyacrylamide gel electrophoresis and reversed-phase high-performance liquid chromatography and tested the venoms' protease, phospholipase A2 and cytotoxic activities. Phospholipase A2 activity was similar in both phenotypes, whereas general protease activity was higher in the melanistic venom, which was also more cytotoxic at two concentrations (6.25 and 12.5 µg ml-1). These minor differences between the venoms of melanistic and normally coloured adders are unlikely to be of clinical relevance in the context of human envenomation. In light of our results, the claim that melanistic adders produce more toxic venom than their normally coloured conspecifics appears rooted entirely in folklore.

Importance of the Cysteine-Rich Domain of Snake Venom Prothrombin Activators: Insights Gained from Synthetic Neutralizing Antibodies

Snake venoms are cocktails of biologically active molecules that have evolved to immobilize prey, but can also induce a severe pathology in humans that are bitten. While animal-derived polyclonal antivenoms are the primary treatment for snakebites, they often have limitations in efficacy and can cause severe adverse side effects. Building on recent efforts to develop improved antivenoms, notably through monoclonal antibodies, requires a comprehensive understanding of venom toxins. Among these toxins, snake venom metalloproteinases (SVMPs) play a pivotal role, particularly in viper envenomation, causing tissue damage, hemorrhage and coagulation disruption. One of the current challenges in the development of neutralizing monoclonal antibodies against SVMPs is the large size of the protein and the lack of existing knowledge of neutralizing epitopes. Here, we screened a synthetic human antibody library to isolate monoclonal antibodies against an SVMP from saw-scaled viper (genus Echis) venom. Upon characterization, several antibodies were identified that effectively blocked SVMP-mediated prothrombin activation. Cryo-electron microscopy revealed the structural basis of antibody-mediated neutralization, pinpointing the non-catalytic cysteine-rich domain of SVMPs as a crucial target. These findings emphasize the importance of understanding the molecular mechanisms of SVMPs to counter their toxic effects, thus advancing the development of more effective antivenoms.

Molecular dissection of cobra venom highlights heparinoids as an antidote for spitting cobra envenoming

Snakebites affect about 1.8 million people annually. The current standard of care involves antibody-based antivenoms, which can be difficult to access and are generally not effective against local tissue injury, the primary cause of morbidity. Here, we used a pooled whole-genome CRISPR knockout screen to define human genes that, when targeted, modify cell responses to spitting cobra venoms. A large portion of modifying genes that conferred resistance to venom cytotoxicity was found to control proteoglycan biosynthesis, including EXT1, B4GALT7, EXT2, EXTL3, XYLT2, NDST1, and SLC35B2, which we validated independently. This finding suggested heparinoids as possible inhibitors. Heparinoids prevented venom cytotoxicity through binding to three-finger cytotoxins, and the US Food and Drug Administration-approved heparinoid tinzaparin was found to reduce tissue damage in mice when given via a medically relevant route and dose. Overall, our systematic molecular dissection of cobra venom cytotoxicity provides insight into how we can better treat cobra snakebite envenoming.

A comparison of the venom proteomes and potential therapeutics of 3 African naja subgenera

African cobras (Naja species) represent one of the most encountered medically important snakes in Africa. They are classified as African spitting (Afronaja subgenus) and non-spitting cobras (Uraeus and Boulengerina subgenera) with similar and different characteristics. Snake venom toxins including three-finger toxin (3FTx), phospholipase A2 (PLA2), and snake venom metalloproteinase (SVMP) cause snakebite envenomation leading to morbidity and mortality. The profile of the proteome of African cobra venoms will help to develop safer and more effective antivenoms. The approval of Captopril by the US Food and Drug Administration (FDA) for the treatment of cardiovascular diseases, has led to intensified research towards possible use of venom toxins as therapeutics. In this review, we compare the venom proteome profile of 3 African Naja subgenera. In both Afronaja and Boulengerina subgenera, 3FTx (Afronaja-69.79%; Boulengerina-60.56%) followed by PLA2 (Afronaja-21.15%; Boulengerina-20.21%) dominated the venoms compared to the Uraeus subgenus dominated by 3FTx (84.55%) with little to no PLA2 abundance (0.8%). The venom of subgenus Uraeus was distinct from the other two subgenera by the almost total absence of PLA2, thus indicating little or no contribution of PLA2 in the envenomation caused by Uraeus compared to Afronaja and Boulengerina. Furthermore, we report studies on the experimental testing of African cobra venoms and toxins against diseases including anti-cancer properties.

Investigating Snake-Venom-Induced Dermonecrosis and Inflammation Using an Ex Vivo Human Skin Model

Snakebite envenoming is a neglected tropical disease that causes >100,000 deaths and >400,000 cases of morbidity annually. Despite the use of mouse models, severe local envenoming, defined by morbidity-causing local tissue necrosis, remains poorly understood, and human-tissue responses are ill-defined. Here, for the first time, an ex vivo, non-perfused human skin model was used to investigate temporal histopathological and immunological changes following subcutaneous injections of venoms from medically important African vipers (Echis ocellatus and Bitis arietans) and cobras (Naja nigricollis and N. haje). Histological analysis of venom-injected ex vivo human skin biopsies revealed morphological changes in the epidermis (ballooning degeneration, erosion, and ulceration) comparable to clinical signs of local envenoming. Immunostaining of these biopsies confirmed cell apoptosis consistent with the onset of necrosis. RNA sequencing, multiplex bead arrays, and ELISAs demonstrated that venom-injected human skin biopsies exhibited higher rates of transcription and expression of chemokines (CXCL5, MIP1-ALPHA, RANTES, MCP-1, and MIG), cytokines (IL-1β, IL-1RA, G-CSF/CSF-3, and GM-CSF), and growth factors (VEGF-A, FGF, and HGF) in comparison to non-injected biopsies. To investigate the efficacy of antivenom, SAIMR Echis monovalent or SAIMR polyvalent antivenom was injected one hour following E. ocellatus or N. nigricollis venom treatment, respectively, and although antivenom did not prevent venom-induced dermal tissue damage, it did reduce all pro-inflammatory chemokines, cytokines, and growth factors to normal levels after 48 h. This ex vivo skin model could be useful for studies evaluating the progression of local envenoming and the efficacy of snakebite treatments.

Specific Amino Acid Residues in the Three Loops of Snake Cytotoxins Determine Their Membrane Activity and Provide a Rationale for a New Classification of These Toxins

Cytotoxins (CTs) are three-finger membrane-active toxins present mainly in cobra venom. Our analysis of the available CT amino acid sequences, literature data on their membrane activity, and conformational equilibria in aqueous solution and detergent micelles allowed us to identify specific amino acid residues which interfere with CT incorporation into membranes. They include Pro9, Ser28, and Asn/Asp45 within the N-terminal, central, and C-terminal loops, respectively. There is a hierarchy in the effect of these residues on membrane activity: Pro9 > Ser28 > Asn/Asp45. Taking into account all the possible combinations of special residues, we propose to divide CTs into eight groups. Group 1 includes toxins containing all of the above residues. Their representatives demonstrated the lowest membrane activity. Group 8 combines CTs that lack these residues. For the toxins from this group, the greatest membrane activity was observed. We predict that when solely membrane activity determines the cytotoxic effects, the activity of CTs from a group with a higher number should exceed that of CTs from a group with a lower number. This classification is supported by the available data on the cytotoxicity and membranotropic properties of CTs. We hypothesize that the special amino acid residues within the loops of the CT molecule may indicate their involvement in the interaction with non-lipid targets.

Development of a membrane-disruption assay using phospholipid vesicles as a proxy for the detection of cellular membrane degradation

Snakebite envenoming is a global health issue that affects millions of people worldwide, and that causes morbidity rates surpassing 450,000 individuals annually. Patients suffering from snakebite morbidities may experience permanent disabilities such as pain, blindness and amputations. The (local) tissue damage that causes these life-long morbidities is the result of cell- and tissue-damaging toxins present in the venoms. These compounds belong to a variety of toxin classes and may affect cells in various ways, for example, by affecting the cell membrane. In this study, we have developed a high-throughput in vitro assay that can be used to study membrane disruption caused by snake venoms using phospholipid vesicles from egg yolk as a substrate. Resuspended chicken egg yolk was used to form these vesicles, which were fluorescently stained to allow monitoring of the degradation of egg yolk vesicles on a plate reader. The assay proved to be suitable for studying phospholipid vesicle degradation of crude venoms and was also tested for its applicability for neutralisation studies of varespladib, which is a PLA2 inhibitor. We additionally made an effort to identify the responsible toxins using liquid chromatography, followed by post-column bioassaying and protein identification using high-throughput venomics. We successfully identified various toxins in the venoms of C. rhodostoma and N. mossambica, which are likely to be involved in the observed vesicle-degrading effect. This indicates that the assay can be used for screening the membrane degrading activity of both crude and fractionated venoms as well as for neutralisation studies.

De novo designed proteins neutralize lethal snake venom toxins

Snakebite envenoming remains a devastating and neglected tropical disease, claiming over 100,000 lives annually and causing severe complications and long-lasting disabilities for many more1,2. Three-finger toxins (3FTx) are highly toxic components of elapid snake venoms that can cause diverse pathologies, including severe tissue damage3 and inhibition of nicotinic acetylcholine receptors (nAChRs) resulting in life-threatening neurotoxicity4. Currently, the only available treatments for snakebite consist of polyclonal antibodies derived from the plasma of immunized animals, which have high cost and limited efficacy against 3FTxs5,6,7. Here, we use deep learning methods to de novo design proteins to bind short- and long-chain α-neurotoxins and cytotoxins from the 3FTx family. With limited experimental screening, we obtain protein designs with remarkable thermal stability, high binding affinity, and near-atomic level agreement with the computational models. The designed proteins effectively neutralize all three 3FTx sub-families in vitro and protect mice from a lethal neurotoxin challenge. Such potent, stable, and readily manufacturable toxin-neutralizing proteins could provide the basis for safer, cost-effective, and widely accessible next-generation antivenom therapeutics. Beyond snakebite, our computational design methodology should help democratize therapeutic discovery, particularly in resource-limited settings, by substantially reducing costs and resource requirements for development of therapies to neglected tropical diseases.

Dermonecrosis caused by a spitting cobra snakebite results from toxin potentiation and is prevented by the repurposed drug varespladib

Snakebite envenoming is a neglected tropical disease that causes substantial mortality and morbidity globally. The venom of African spitting cobras often causes permanent injury via tissue-destructive dermonecrosis at the bite site, which is ineffectively treated by current antivenoms. To address this therapeutic gap, we identified the etiological venom toxins in Naja nigricollis venom responsible for causing local dermonecrosis. While cytotoxic three-finger toxins were primarily responsible for causing spitting cobra cytotoxicity in cultured keratinocytes, their potentiation by phospholipases A2 toxins was essential to cause dermonecrosis in vivo. This evidence of probable toxin synergism suggests that a single toxin-family inhibiting drug could prevent local envenoming. We show that local injection with the repurposed phospholipase A2-inhibiting drug varespladib significantly prevents local tissue damage caused by several spitting cobra venoms in murine models of envenoming. Our findings therefore provide a therapeutic strategy that may effectively prevent life-changing morbidity caused by snakebite in rural Africa.

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.

Distinct cardiotoxic effects by venoms of a spitting cobra (Naja pallida) and a rattlesnake (Crotalus atrox) revealed using an ex vivo Langendorff heart model

Here we describe the acute myocardial effects of an elapid (red spitting cobra, Naja pallida) and a viper (western diamondback rattlesnake, Crotalus atrox) venom using an ex vivo heart model. Our results reveal two different pathophysiological trajectories that influence heart function and morphology. While cobra venom causes a drop in contractile force, rattlesnake venom causes enhanced contractility and frequency that coincides with differences in myocellular morphology. This highlights the medical complexity of snake venom-induced cardiotoxicity.

Synthetic development of a broadly neutralizing antibody against snake venom long-chain α-neurotoxins

Snakebite envenoming is a major global public health concern for which improved therapies are urgently needed. The antigenic diversity present in snake venom toxins from various species presents a considerable challenge to the development of a universal antivenom. Here, we used a synthetic human antibody library to find and develop an antibody that neutralizes long-chain three-finger α-neurotoxins produced by numerous medically relevant snakes. Our antibody bound diverse toxin variants with high affinity, blocked toxin binding to the nicotinic acetylcholine receptor in vitro, and protected mice from lethal venom challenge. Structural analysis of the antibody-toxin complex revealed a binding mode that mimics the receptor-toxin interaction. The overall workflow presented is generalizable for the development of antibodies that target conserved epitopes among antigenically diverse targets, and it offers a promising framework for the creation of a monoclonal antibody-based universal antivenom to treat snakebite envenoming.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Peptide toxins that target vertebrate voltage-gated sodium channels underly the painful stings of harvester ants

Harvester ants (genus Pogonomyrmex) are renowned for their stings which cause intense, long-lasting pain, and other neurotoxic symptoms in vertebrates. Here, we show that harvester ant venoms are relatively simple and composed largely of peptide toxins. One class of peptides is primarily responsible for the long-lasting local pain of envenomation via activation of peripheral sensory neurons. These hydrophobic, cysteine-free peptides potently modulate mammalian voltage-gated sodium (NaV) channels, reducing the voltage threshold for activation and inhibiting channel inactivation. These toxins appear to have evolved specifically to deter vertebrates.

Recombinant human scFv antibody fragments against phospholipase A2 from Naja naja and Echis carinatus snake venoms: In vivo neutralization and mechanistic insights

Snake envenomation results in a range of clinical sequelae, and widely used animal-based conventional antivenoms exhibit several limitations including the adverse immunological effects in human snake bite victims. Therefore, human monoclonal anti-snake venom antibodies or fragments can be an alternate therapy for overcoming the existing limitations. We developed venom-neutralizing humanized scFv antibodies and analyzed biochemical mechanisms associated with the inhibition of toxicity. Tomlinson I and J human scFv antibody libraries were screened against Naja naja and Echis carinatus venoms, and seven unique scFv antibodies were obtained. Further, specific toxins of snake venom interacting with each of these scFvs were identified, and phospholipase A2 (PLA2) was found to be prominently captured by the phage-anchored scFv antibodies. Our study indicated PLA2 to be one of the abundant toxins in Naja naja and Echis carinatus venom samples. The scFvs binding to PLA2 were used to perform in vivo survival assay using the mouse model and in vitro toxin inhibition assays. scFv N194, which binds to acidic PLA2, protected 50% of mice treated with Naja naja venom. Significant prolongation of survival time and 16% survival were observed in Echis carinatus venom-challenged mice treated with scFv E113 and scFv E10, respectively. However, a combination comprised of an equal amount of two scFvs, E113 and E10, both interacting with basic PLA2, exhibited synergistically enhanced survival of 33% in Echis carinatus venom-challenged mice. No such synergistically enhanced survival was observed in the case of combinatorial treatment with anti-Naja naja scFvs, N194, and N248. These scFvs demonstrated partial inhibition of venom-induced myotoxicity, and E113 also inhibited hemolysis by 50%, which corroborates the enhanced survival during combinatorial treatment in Echis carinatus venom-challenged mice.

Repurposed drugs and their combinations prevent morbidity-inducing dermonecrosis caused by diverse cytotoxic snake venoms

Morbidity from snakebite envenoming affects approximately 400,000 people annually. Tissue damage at the bite-site often leaves victims with catastrophic life-long injuries and is largely untreatable by current antivenoms. Repurposed small molecule drugs that inhibit specific snake venom toxins show considerable promise for tackling this neglected tropical disease. Using human skin cell assays as an initial model for snakebite-induced dermonecrosis, we show that the drugs 2,3-dimercapto-1-propanesulfonic acid (DMPS), marimastat, and varespladib, alone or in combination, inhibit the cytotoxicity of a broad range of medically important snake venoms. Thereafter, using preclinical mouse models of dermonecrosis, we demonstrate that the dual therapeutic combinations of DMPS or marimastat with varespladib significantly inhibit the dermonecrotic activity of geographically distinct and medically important snake venoms, even when the drug combinations are delivered one hour after envenoming. These findings strongly support the future translation of repurposed drug combinations as broad-spectrum therapeutics for preventing morbidity caused by snakebite.

Venomous Noodles: The Evolution of Toxins in Nemertea through Positive Selection and Gene Duplication

Some, probably most and perhaps all, members of the phylum Nemertea are poisonous, documented so far from marine and benthic specimens. Although the toxicity of these animals has been long known, systematic studies on the characterization of toxins, mechanisms of toxicity, and toxin evolution for this group are scarce. Here, we present the first investigation of the molecular evolution of toxins in Nemertea. Using a proteo-transcriptomic approach, we described toxins in the body and poisonous mucus of the pilidiophoran Lineus sanguineus and the hoplonemertean Nemertopsis pamelaroeae. Using these new and publicly available transcriptomes, we investigated the molecular evolution of six selected toxin gene families. In addition, we also characterized in silico the toxin genes found in the interstitial hoplonemertean, Ototyphlonemertes erneba, a meiofaunal taxa. We successfully identified over 200 toxin transcripts in each of these species. Evidence of positive selection and gene duplication was observed in all investigated toxin genes. We hypothesized that the increased rates of gene duplications observed for Pilidiophora could be involved with the expansion of toxin genes. Studies concerning the natural history of Nemertea are still needed to understand the evolution of their toxins. Nevertheless, our results show evolutionary mechanisms similar to other venomous groups.

An in vitro assay to investigate venom neurotoxin activity on muscle-type nicotinic acetylcholine receptor activation and for the discovery of toxin-inhibitory molecules

Snakebite envenoming is a neglected tropical disease that causes over 100,000 deaths annually. Envenomings result in variable pathologies, but systemic neurotoxicity is among the most serious and is currently only treated with difficult to access and variably efficacious commercial antivenoms. Venom-induced neurotoxicity is often caused by α-neurotoxins antagonising the muscle-type nicotinic acetylcholine receptor (nAChR), a ligand-gated ion channel. Discovery of therapeutics targeting α-neurotoxins is hampered by relying on binding assays that do not reveal restoration of receptor activity or more costly and/or lower throughput electrophysiology-based approaches. Here, we report the validation of a screening assay for nAChR activation using immortalised TE671 cells expressing the γ-subunit containing muscle-type nAChR and a fluorescent dye that reports changes in cell membrane potential. Assay validation using traditional nAChR agonists and antagonists, which either activate or block ion fluxes, was consistent with previous studies. We then characterised antagonism of the nAChR by a variety of elapid snake venoms that cause muscle paralysis in snakebite victims, before defining the toxin-inhibiting activities of commercial antivenoms, and new types of snakebite therapeutic candidates, namely monoclonal antibodies, decoy receptors, and small molecules. Our findings show robust evidence of assay uniformity across 96-well plates and highlight the amenability of this approach for the future discovery of new snakebite therapeutics via screening campaigns. The described assay therefore represents a useful first-step approach for identifying α-neurotoxins and their inhibitors in the context of snakebite envenoming, and it should provide wider value for studying modulators of nAChR activity from other sources.

Museum DNA reveals a new, potentially extinct species of rinkhals (Serpentes: Elapidae: Hemachatus) from the Eastern Highlands of Zimbabwe

Genetic information plays a pivotal role in species recognition and delimitation, but rare or extinct animals can be difficult to obtain genetic samples from. While natural history wet collections have proven invaluable in the description of novel species, the use of these historical samples in genetic studies has been greatly impeded by DNA degradation, especially because of formalin-fixation prior to preservation. Here, we use recently developed museum genomics approaches to determine the status of an isolated population of the elapid snake genus Hemachatus from Zimbabwe. We used multiple digestion phases followed by single strand sequencing library construction and hybridisation capture to obtain 12S and 16S rDNA sequences from a poorly preserved tissue sample of this population. Phylogenetic and morphological analyses in an integrated taxonomic framework demonstrate that the Zimbabwean rinkhals population represents an old and highly distinct lineage, which we describe as a new species, Hemachatus nyangensis sp. nov. Our phylogenetic dating analysis is compatible with venom spitting having evolved in response to the threat posed by early hominins, although more data are required for a robust test of this hypothesis. This description demonstrates the power of museum genomics in revealing rare or even extinct species: Hemachatus from Zimbabwe are only known from a small area of the Eastern Highlands known for high endemism. No living specimens have been seen since the 1980s, most likely due to dramatic land-use changes in the Eastern Highlands, suggesting that the species could be extinct. In view of its recognition as a highly distinct lineage, urgent action is required to determine whether any populations survive, and to safeguard remaining habitat.

On the importance of types and the perils of en passant taxonomy: a brief history of the typification of Coluber naja Linnaeus, 1758 (Serpentes: Elapidae) and its implications, with the designation of a lectotype

In response to the recent in passing (en passant) taxonomic decision to split Naja naja (Linnaeus) and recognise the Sri Lankan populations as a separate species, N. polyocellata Deraniyagala, we analyse the evidence underlying the proposal and its nomenclatural implications. The proposed split is weakly supported by the available evidence, so that retaining N. naja as a single species seems appropriate until further analysis. Moreover, the proposal raises several issues concerning types, type locality and nomenclature. Linnaeus description of Coluber naja was based on a single preserved specimen seen by him (now lost) and several illustrations in Sebas Thesaurus. The specimens that were the basis of these illustrations constitute part of the type series. Two of the latter specimens, ZMB 2795 and 2796, have been rediscovered in the collections of the Museum fr Naturkunde, Berlin. Here, we describe them, and determine that both are of Sri Lankan origin. To settle the question of the type and type locality of this iconic taxon, we designate ZMB 2796 as lectotype for the species, thereby implicitly restricting the type locality to Sri Lanka. The name polyocellata thus becomes a subjective junior synonym of Coluber naja, and the name Naja brasiliensis Laurenti, 1768 an objective junior synonym thereof. Any taxonomic recognition of additional diversity within N. naja would thus require the renaming of Indian, not Sri Lankan spectacled cobras, but should await a significant body of convincing evidence. We caution against taxonomic decisions taken in passing, based on limited evidence and without in-depth assessment of their nomenclatural implications.

Malaysian and Chinese King Cobra Venom Cytotoxicity in Melanoma and Neonatal Foreskin Fibroblasts Is Mediated by Age and Geography

Snake venoms constitute a complex, rapidly evolving trait, whose composition varies between and within populations depending on geographical location, age and preys (diets). These factors have determined the adaptive evolution for predatory success and link venom heterogeneity with prey specificity. Moreover, understanding the evolutionary drivers of animal venoms has streamlined the biodiscovery of venom-derived compounds as drug candidates in biomedicine and biotechnology. The king cobra (Ophiophagus hannah; Cantor, 1836) is distributed in diverse habitats, forming independent populations, which confer differing scale markings, including between hatchlings and adults. Furthermore, king cobra venoms possess unique cytotoxic properties that are used as a defensive trait, but their toxins may also have utility as promising anticancer-agent candidates. However, the impact of geographical distribution and age on these potential venom applications has been typically neglected. In this study, we hypothesised that ontogenetic venom variation accompanies the morphological distinction between hatchlings and adults. We used non-transformed neonatal foreskin (NFF) fibroblasts to examine and compare the variability of venom cytotoxicity between adult captive breeding pairs from Malaysian and Chinese lineages, along with that of their progeny upon hatching. In parallel, we assessed the anticancer potential of these venoms in human-melanoma-patient-derived cells (MM96L). We found that in a geographical distribution and gender-independent manner, venoms from hatchlings were significantly less cytotoxic than those from adults (NFF; ~Log EC50: 0.5-0.6 vs. 0.2-0.35 mg/mL). This is consistent with neonates occupying a semifossorial habitat, while adults inhabit more above-ground habitats and are therefore more conspicuous to potential predators. We also observed that Malaysian venoms exhibited a slightly higher cytotoxicity than those from the Chinese cobra cohorts (NFF; Log EC50: 0.1-0.3 vs. 0.3-0.4 mg/mL), which is consistent with Malaysian king cobras being more strongly aposematically marked. These variations are therefore suggestive of differential anti-predator strategies associated with the occupation of distinct niches. However, all cobra venoms were similarly cytotoxic in both melanoma cells and fibroblasts, limiting their potential medical applications in their native forms.

Slowly Making Sense: A Review of the Two-Step Venom System within Slow (Nycticebus spp.) and Pygmy Lorises (Xanthonycticebus spp.)

Since the early 2000s, studies of the evolution of venom within animals have rapidly expanded, offering new revelations on the origins and development of venom within various species. The venomous mammals represent excellent opportunities to study venom evolution due to the varying functional usages, the unusual distribution of venom across unrelated mammals and the diverse variety of delivery systems. A group of mammals that excellently represents a combination of these traits are the slow (Nycticebus spp.) and pygmy lorises (Xanthonycticebus spp.) of south-east Asia, which possess the only confirmed two-step venom system. These taxa also present one of the most intriguing mixes of toxic symptoms (cytotoxicity and immunotoxicity) and functional usages (intraspecific competition and ectoparasitic defence) seen in extant animals. We still lack many pieces of the puzzle in understanding how this venom system works, why it evolved what is involved in the venom system and what triggers the toxic components to work. Here, we review available data building upon a decade of research on this topic, focusing especially on why and how this venom system may have evolved. We discuss that research now suggests that venom in slow lorises has a sophisticated set of multiple uses in both intraspecific competition and the potential to disrupt the immune system of targets; we suggest that an exudate diet reveals several toxic plants consumed by slow and pygmy lorises that could be sequestered into their venom and which may help heal venomous bite wounds; we provide the most up-to-date visual model of the brachial gland exudate secretion protein (BGEsp); and we discuss research on a complement component 1r (C1R) protein in saliva that may solve the mystery of what activates the toxicity of slow and pygmy loris venom. We conclude that the slow and pygmy lorises possess amongst the most complex venom system in extant animals, and while we have still a lot more to understand about their venom system, we are close to a breakthrough, particularly with current technological advances.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila , many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila . We first inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of these species could develop on diet with toxin. Third, we examined fly development on a diet containing a natural mix of toxins extracted from the Death Cap Amanita phalloides mushroom. Both tolerant and susceptible species developed on diet with this mix, though tolerant species survived at significantly higher concentrations. Finally, we asked how cyclopeptide tolerance might have evolved across the immigrans-tripunctata radiation and inferred that toxin tolerance was ancestral and subsequently lost multiple times. Our results suggest the evolutionary history of cyclopeptide tolerance is complex, and simply describing this trait as present or absent does not fully capture the occurrence or impact on this adaptive radiation. More broadly, the evolution of novelty can be more complex than previously thought, and that accurate descriptions of such novelties are critical in studies examining their evolution.

Resistance Is Not Futile: Widespread Convergent Evolution of Resistance to Alpha-Neurotoxic Snake Venoms in Caecilians (Amphibia: Gymnophiona)

Predatory innovations impose reciprocal selection pressures upon prey. The evolution of snake venom alpha-neurotoxins has triggered the corresponding evolution of resistance in the post-synaptic nicotinic acetylcholine receptors of prey in a complex chemical arms race. All other things being equal, animals like caecilians (an Order of legless amphibians) are quite vulnerable to predation by fossorial elapid snakes and their powerful alpha-neurotoxic venoms; thus, they are under strong selective pressure. Here, we sequenced the nicotinic acetylcholine receptor alpha-1 subunit of 37 caecilian species, representing all currently known families of caecilians from across the Americas, Africa, and Asia, including species endemic to the Seychelles. Three types of resistance were identified: (1) steric hindrance from N-glycosylated asparagines; (2) secondary structural changes due to the replacement of proline by another amino acid; and (3) electrostatic charge repulsion of the positively charged neurotoxins, through the introduction of a positively charged amino acid into the toxin-binding site. We demonstrated that resistance to alpha-neurotoxins convergently evolved at least fifteen times across the caecilian tree (three times in Africa, seven times in the Americas, and five times in Asia). Additionally, as several species were shown to possess multiple resistance modifications acting synergistically, caecilians must have undergone at least 20 separate events involving the origin of toxin resistance. On the other hand, resistance in non-caecilian amphibians was found to be limited to five origins. Together, the mutations underlying resistance in caecilians constitute a robust signature of positive selection which strongly correlates with elapid presence through both space (sympatry with caecilian-eating elapids) and time (Cenozoic radiation of elapids). Our study demonstrates the extent of convergent evolution that can be expected when a single widespread predatory adaptation triggers parallel evolutionary arms races at a global scale.

Backward masking reveals coarse-to-fine dynamics in human V1

Natural images exhibit luminance variations aligned across a broad spectrum of spatial frequencies (SFs). It has been proposed that, at early stages of processing, the coarse signals carried by the low SF (LSF) of the visual input are sent rapidly from primary visual cortex (V1) to ventral, dorsal and frontal regions to form a coarse representation of the input, which is later sent back to V1 to guide the processing of fine-grained high SFs (i.e., HSF). We used functional resonance imaging (fMRI) to investigate the role of human V1 in the coarse-to-fine integration of visual input. We disrupted the processing of the coarse and fine content of full-spectrum human face stimuli via backward masking of selective SF ranges (LSFs: <1.75cpd and HSFs: >1.75cpd) at specific times (50, 83, 100 or 150ms). In line with coarse-to-fine proposals, we found that (1) the selective masking of stimulus LSF disrupted V1 activity in the earliest time window, and progressively decreased in influence, while (2) an opposite trend was observed for the masking of stimulus' HSF. This pattern of activity was found in V1, as well as in ventral (i.e. the Fusiform Face area, FFA), dorsal and orbitofrontal regions. We additionally presented subjects with contrast negated stimuli. While contrast negation significantly reduced response amplitudes in the FFA, as well as coupling between FFA and V1, coarse-to-fine dynamics were not affected by this manipulation. The fact that V1 response dynamics to strictly identical stimulus sets differed depending on the masked scale adds to growing evidence that V1 role goes beyond the early and quasi-passive transmission of visual information to the rest of the brain. It instead indicates that V1 may yield a 'spatially registered common forum' or 'blackboard' that integrates top-down inferences with incoming visual signals through its recurrent interaction with high-level regions located in the inferotemporal, dorsal and frontal regions.

Multilevel Comparison of Indian Naja Venoms and Their Cross-Reactivity with Indian Polyvalent Antivenoms

Snake envenoming is caused by many biological species, rather than a single infectious agent, each with a multiplicity of toxins in their venom. Hence, developing effective treatments is challenging, especially in biodiverse and biogeographically complex countries such as India. The present study represents the first genus-wide proteomics analysis of venom composition across Naja species (N. naja, N. oxiana, and N. kaouthia) found in mainland India. Venom proteomes were consistent between individuals from the same localities in terms of the toxin families present, but not in the relative abundance of those in the venom. There appears to be more compositional variation among N. naja from different locations than among N. kaouthia. Immunoblotting and in vitro neutralization assays indicated cross-reactivity with Indian polyvalent antivenom, in which antibodies raised against N. naja are present. However, we observed ineffective neutralization of PLA2 activities of N. naja venoms from locations distant from the source of immunizing venoms. Antivenom immunoprofiling by antivenomics revealed differential antigenicity of venoms from N. kaouthia and N. oxiana, and poor reactivity towards 3FTxs and PLA2s. Moreover, there was considerable variation between antivenoms from different manufacturers. These data indicate that improvements to antivenom manufacturing in India are highly desirable.

Diversity Begets Diversity When Diet Drives Snake Venom Evolution, but Evenness Rather Than Richness Is What Counts

Snake venoms are primarily used to subjugate prey, and consequently, their evolution has been shown to be predominantly driven by diet-related selection pressure. Venoms tend to be more lethal to prey than non-prey species (except in cases of toxin resistance), prey-specific toxins have been identified, and preliminary work has demonstrated an association between the diversity of diet classes and that of toxicological activities of whole venom. However, venoms are complex mixtures of many toxins, and it remains unclear how toxin diversity is driven by diet. Prey-specific toxins do not encompass the molecular diversity of venoms, and whole venom effects could be driven by one, few, or all components, so the link between diet and venom diversity remains minimally understood. Here, we collated a database of venom composition and diet records and used a combination of phylogenetic comparative methods and two quantitative diversity indices to investigate whether and how diet diversity relates to the toxin diversity of snake venoms. We reveal that venom diversity is negatively related to diet diversity using Shannon's index but positively related using Simpson's index. Since Shannon's index predominantly considers the number of prey/toxins, whereas Simpson's index more strongly reflects evenness, we provide insights into how the diet-venom diversity link is driven. Specifically, species with low diet diversity tend to have venoms dominated by a few abundant (possibly specialised) toxin families, whereas species with diverse diets tend to 'hedge their bets' by having venoms with a more even composition of different toxin classes.

An expanded description, natural history, and genetic variation of the recently described cobra species Naja fuxi Shi et al., 2022

The morphological variation, extended distribution, and sequence divergence of a recently described of cobra Naja fuxiShi et al., 2022 captured from mountainous areas in Thailand are evaluated by using molecular and morphological analyses. We investigated the genetic variation and affinities of 72 specimens in the genus Naja by using mitochondrial DNA (cytochrome b and control region) and the nuclear DNA gene, C-mos. Morphological examination was conducted for 33 cobra specimens obtained from the northern, western, and north-eastern regions, and data on their natural history were gathered during field surveys. A high degree of genetic differentiation was shown to exist between the cobras collected from lowlands and those from mountainous areas. N. fuxi occurs in uplands bordering Thailand’s Central Basin, whereas the similar looking N. kaouthia Lesson, 1831 is more or less restricted to the lowlands. All phylogenetic and network analyses supported a distinct clade of N. fuxi from north, west, and, north-east regions. In addition, N. fuxi seems to exhibit a split between the north-eastern population and those from the north and west. The range of N. fuxi probably extends far into the mountainous areas of the neighbouring countries Myanmar, Laos, and Vietnam. Morphologically, N. fuxi in Thailand can be distinguished from all other cobra species in the adjacent Oriental Region. The speciation of cobras in Thailand likely reflects key events in the region’s geographical, climate and environmental history.

Metabolome-Based Classification of Snake Venoms by Bioinformatic Tools

Snakebite is considered a neglected tropical disease, and it is one of the most intricate ones. The variability found in snake venom is what makes it immensely complex to study. These variations are present both in the big and the small molecules found in snake venom. This study focused on examining the variability found in the venom's small molecules (i.e., mass range of 100-1000 Da) between two main families of venomous snakes-Elapidae and Viperidae-managing to create a model able to classify unknown samples by means of specific features, which can be extracted from their LC-MS data and output in a comprehensive list. The developed model also allowed further insight into the composition of snake venom by highlighting the most relevant metabolites of each group by clustering similarly composed venoms. The model was created by means of support vector machines and used 20 features, which were merged into 10 principal components. All samples from the first and second validation data subsets were correctly classified. Biological hypotheses relevant to the variation regarding the metabolites that were identified are also given.

Micro and macroevolution of sea anemone venom phenotype

Venom is a complex trait with substantial inter- and intraspecific variability resulting from strong selective pressures acting on the expression of many toxic proteins. However, understanding the processes underlying toxin expression dynamics that determine the venom phenotype remains unresolved. By interspecific comparisons we reveal that toxin expression in sea anemones evolves rapidly and that in each species different toxin family dictates the venom phenotype by massive gene duplication events. In-depth analysis of the sea anemone, Nematostella vectensis, revealed striking variation of the dominant toxin (Nv1) diploid copy number across populations (1-24 copies) resulting from independent expansion/contraction events, which generate distinct haplotypes. Nv1 copy number correlates with expression at both the transcript and protein levels with one population having a near-complete loss of Nv1 production. Finally, we establish the dominant toxin hypothesis which incorporates observations in other venomous lineages that animals have convergently evolved a similar strategy in shaping their venom.

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.

Intra-colony venom diversity contributes to maintaining eusociality in a cooperatively breeding ant

BACKGROUND

Eusociality is widely considered to evolve through kin selection, where the reproductive success of an individual's close relative is favored at the expense of its own. High genetic relatedness is thus considered a prerequisite for eusociality. While ants are textbook examples of eusocial animals, not all ants form colonies of closely related individuals. One such example is the ectatommine ant Rhytidoponera metallica, which predominantly forms queen-less colonies that have such a low intra-colony relatedness that they have been proposed to represent a transient, unstable form of eusociality. However, R. metallica is among the most abundant and widespread ants on the Australian continent. This apparent contradiction provides an example of how inclusive fitness may not by itself explain the maintenance of eusociality and raises the question of what other selective advantages maintain the eusocial lifestyle of this species.

RESULTS

We provide a comprehensive portrait of the venom of R. metallica and show that the colony-wide venom consists of an exceptionally high diversity of functionally distinct toxins for an ant. These toxins have evolved under strong positive selection, which is normally expected to reduce genetic variance. Yet, R. metallica exhibits remarkable intra-colony variation, with workers sharing only a relatively small proportion of toxins in their venoms. This variation is not due to the presence of chemical castes, but has a genetic foundation that is at least in part explained by toxin allelic diversity.

CONCLUSIONS

Taken together, our results suggest that the toxin diversity contained in R. metallica colonies may be maintained by a form of group selection that selects for colonies that can exploit more resources and defend against a wider range of predators. We propose that increased intra-colony genetic variance resulting from low kinship may itself provide a selective advantage in the form of an expanded pharmacological venom repertoire. These findings provide an example of how group selection on adaptive phenotypes may contribute to maintaining eusociality where a prerequisite for kin selection is diminished.

Spatial Venomics─Cobra Venom System Reveals Spatial Differentiation of Snake Toxins by Mass Spectrometry Imaging

Among venomous animals, toxic secretions have evolved as biochemical weapons associated with various highly specialized delivery systems on many occasions. Despite extensive research, there is still limited knowledge of the functional biology of most animal toxins, including their venom production and storage, as well as the morphological structures within sophisticated venom producing tissues that might underpin venom modulation. Here, we report on the spatial exploration of a snake venom gland system by matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), in combination with standard proteotranscriptomic approaches, to enable in situ toxin mapping in spatial intensity maps across a venom gland sourced from the Egyptian cobra (Naja haje). MALDI-MSI toxin visualization on the elapid venom gland reveals a high spatial heterogeneity of different toxin classes at the proteoform level, which may be the result of physiological constraints on venom production and/or storage that reflects the potential for venom modulation under diverse stimuli.

Venom composition and pain-causing toxins of the Australian great carpenter bee Xylocopa aruana

Most species of bee are capable of delivering a defensive sting which is often painful. A solitary lifestyle is the ancestral state of bees and most extant species are solitary, but information on bee venoms comes predominantly from studies on eusocial species. In this study we investigated the venom composition of the Australian great carpenter bee, Xylocopa aruana Ritsema, 1876. We show that the venom is relatively simple, composed mainly of one small amphipathic peptide (XYTX₁-Xa1a), with lesser amounts of an apamin homologue (XYTX₂-Xa2a) and a venom phospholipase-A₂ (PLA₂). XYTX₁-Xa1a is homologous to, and shares a similar mode-of-action to melittin and the bombilitins, the major components of the venoms of the eusocial Apis mellifera (Western honeybee) and Bombus spp. (bumblebee), respectively. XYTX₁-Xa1a and melittin directly activate mammalian sensory neurons and cause spontaneous pain behaviours in vivo, effects which are potentiated in the presence of venom PLA₂. The apamin-like peptide XYTX₂-Xa2a was a relatively weak blocker of small conductance calcium-activated potassium (K_Ca) channels and, like A. mellifera apamin and mast cell-degranulating peptide, did not contribute to pain behaviours in mice. While the composition and mode-of-action of the venom of X. aruana are similar to that of A. mellifera, the greater potency, on mammalian sensory neurons, of the major pain-causing component in A. mellifera venom may represent an adaptation to the distinct defensive pressures on eusocial Apidae.

AHA: AI-guided tool for the quantification of venom-induced haemorrhage in mice

Venom-induced haemorrhage constitutes a severe pathology in snakebite envenomings, especially those inflicted by viperid species. To both explore venom activity accurately and evaluate the efficacy of viperid antivenoms for the neutralisation of haemorrhagic activity it is essential to have available a precise, quantitative tool for empirically determining venom-induced haemorrhage. Thus, we have built on our prior approach and developed a new AI-guided tool (AHA) for the quantification of venom-induced haemorrhage in mice. Using a smartphone, it takes less than a minute to take a photo, upload the image, and receive accurate information on the magnitude of a venom-induced haemorrhagic lesion in mice. This substantially decreases analysis time, reduces human error, and does not require expert haemorrhage analysis skills. Furthermore, its open access web-based graphical user interface makes it easy to use and implement in laboratories across the globe. Together, this will reduce the resources required to preclinically assess and control the quality of antivenoms, whilst also expediting the profiling of haemorrhagic activity in venoms for the wider toxinology community.

High-throughput proteomics and in vitro functional characterization of the 26 medically most important elapids and vipers from sub-Saharan Africa

Venomous snakes are important parts of the ecosystem, and their behavior and evolution have been shaped by their surrounding environments over the eons. This is reflected in their venoms, which are typically highly adapted for their biological niche, including their diet and defense mechanisms for deterring predators. Sub-Saharan Africa is rich in venomous snake species, of which many are dangerous to humans due to the high toxicity of their venoms and their ability to effectively deliver large amounts of venom into their victims via their bite. In this study, the venoms of 26 of sub-Saharan Africa's medically most relevant elapid and viper species were subjected to parallelized toxicovenomics analysis. The analysis included venom proteomics and in vitro functional characterization of whole venom toxicities, enabling a robust comparison of venom profiles between species. The data presented here corroborate previous studies and provide biochemical details for the clinical manifestations observed in envenomings by the 26 snake species. Moreover, two new venom proteomes (Naja anchietae and Echis leucogaster) are presented here for the first time. Combined, the presented data can help shine light on snake venom evolutionary trends and possibly be used to further improve or develop novel antivenoms.

Description of a New Cobra (Naja Laurenti, 1768; Squamata, Elapidae) from China with Designation of a Neotype for Naja atra

Taxonomic frameworks for medically important species such as cobras (genus Naja Laurenti, 1768; Squamata, Elapidae) are essential for the medical treatment of snake bites and accurate antivenin development. In this paper, we described the former N. kaouthia populations recorded from China as a new species and designated a neotype for N. atra-based morphological and mitochondrial phylogenetic analysis. The new species N. fuxisp. nov. was morphologically diagnosed from N. kaouthia by (1) regular single narrow crossband present on the middle and posterior parts of the dorsum (3-15, 7.9 ± 2.7, n = 32) and the dorsal surface of the tail (1-6, 4.2 ± 1.1, n = 32) of both adults and juveniles, buff-colored with dark fringes on both edges, vs. South Asian populations (n = 39) and Southeast Asian populations (n = 35) without cross bands, with irregular cross bands or multiple light-colored crossbands pairs, or densely woven lines; (2) small scales between the posterior chin shields, usually three (40%) or two (37%), rarely four (13%), or one (10%) (n = 30) vs. mostly one (81%) and rarely two (19%) (n = 28); (3) ventrals 179-205 (195.4 ± 6.7, n = 33) vs. South Asian populations 179-199 (188.7 ± 5.9, n = 12); Southeast Asian populations 168-186 (177.8 ± 4.9, n = 18). Phylogenetically, the new species forms an independent sister clade to the clade including N. atra, N. kaouthia, N. oxiana and N. sagittifera. Furthermore, the subspecies N. naja polyocellata should be resurrected and recognized as a full species, N. polyocellatacomb. nov., and the subspecies N. sumatrana miolepis should be resurrected.