Novel snakebite therapeutics must be tested in appropriate rescue models to robustly assess their preclinical efficacy

In the field of antivenom research, development, and manufacture, it is often advised to follow the World Health Organization’s (WHO) guidelines for the production, control, and regulation of snake antivenom immunoglobulins, which recommend the use of preincubation assays to assess the efficacy of snakebite therapeutics. In these assays, venom and antivenom are mixed and incubated prior to in vivo administration to rodents, which allows for a standardizable comparison of antivenoms with similar characteristics. However, these assays are not necessarily sufficient for therapeutics with significantly different pharmacological properties than antibody-based antivenoms, such as small molecule inhibitors, nanoparticles, and other modalities. To ensure that the in vivo therapeutic utility of completely novel toxin-neutralizing molecules with no history of use in envenoming therapy and variable pharmacokinetics is properly evaluated, such molecules must also be tested in preclinical rescue assays, where rodents are first challenged with appropriate doses of venoms or toxins, followed by the administration of neutralizing modalities after an appropriate time delay to better mimic the real-life scenarios faced by human snakebite victims. Such an approach takes the venom (or toxin) toxicokinetics, the drug pharmacokinetics, and the drug pharmacodynamics into consideration. If new modalities are only assessed in preincubation assays and not subjected to evaluation in rescue assays, the publication of neutralization data may unintentionally misrepresent the actual therapeutic efficacy and suitability of the modality being tested, and thus potentially misguide strategic decision making in the research and development of novel therapies for snakebite envenomingVillum Foundation/[Project 00025302]//DinamarcaBioPorto Diagnostics A/S/[]//DinamarcaUCR::Vicerrectoría de Investigación::Unidades de Investigación::Ciencias de la Salud::Instituto Clodomiro Picado (ICP

Widespread vulnerability of Malagasy predators to the toxins of an introduced toad

Invasive species are a key factor contributing to the global decline of biodiversity, and understanding the underlying mechanisms is crucial to mitigate detrimental effects [1]. One such mechanism is the introduction of invasive species with defensive strategies, such as novel toxins, that can disrupt native predator communities [2]. Disruption of such communities can produce trophic cascades, impacting a diverse array of taxa [2]. Madagascar, a globally significant biodiversity hotspot, has recently experienced the introduction of a toxic bufonid amphibian, the Asian common toad (Duttaphrynus melanostictus) [3]. Since its invasion, the toad population has expanded rapidly, making control efforts problematic and eradication extremely difficult [4]. Previous cases of bufonid introductions, such as the ongoing spread of the cane toad (Rhinella marina) in Australia, have resulted in the decimation of many indigenous species [2], prompting fears that Madagascar may be similarly impacted [4]. Here we show that these fears are warranted: we demonstrate that many Malagasy vertebrates are likely to be susceptible to the toxins of this invasive toad

Biogeographical venom variation in the Indian spectacled cobra (Naja naja) underscores the pressing need for pan-India efficacious snakebite therapy

Background: Snake venom composition is dictated by various ecological and environmental factors, and can exhibit dramatic variation across geographically disparate populations of the same species. This molecular diversity can undermine the efficacy of snakebite treatments, as antivenoms produced against venom from one population may fail to neutralise others. India is the world’s snakebite hotspot, with 58,000 fatalities and 140,000 morbidities occurring annually. Spectacled cobra (Naja naja) and Russell’s viper (Daboia russelii) are known to cause the majority of these envenomations, in part due to their near country-wide distributions. However, the impact of differing ecologies and environment on their venom compositions has not been comprehensively studied. Methods: Here, we used a multi-disciplinary approach consisting of venom proteomics, biochemical and pharmacological analyses, and in vivo research to comparatively analyse N. naja venoms across a broad region (>6000 km; seven populations) covering India’s six distinct biogeographical zones. Findings: By generating the most comprehensive pan-Indian proteomic and toxicity profiles to date, we unveil considerable differences in the composition, pharmacological effects and potencies of geographically-distinct venoms from this species and, through the use of immunological assays and preclinical experiments, demonstrate alarming repercussions on antivenom therapy. We find that commercially-available antivenom fails to effectively neutralise envenomations by the pan-Indian populations of N. naja, including a complete lack of neutralisation against the desert Naja population. Conclusion: Our findings highlight the significant influence of ecology and environment on snake venom composition and potency, and stress the pressing need to innovate pan-India effective antivenoms to safeguard the lives, limbs and livelihoods of the country’s 200,000 annual snakebite victims

An interactive database for the investigation of high-density peptide microarray guided interaction patterns and antivenom cross-reactivity

Snakebite envenoming is a major neglected tropical disease that affects millions of people every year. The only effective treatment against snakebite envenoming consists of unspecified cocktails of polyclonal antibodies purified from the plasma of immunized production animals. Currently, little data exists on the molecular interactions between venom-toxin epitopes and antivenom-antibody paratopes. To address this issue, high-density peptide microarray (hdpm) technology has recently been adapted to the field of toxinology. However, analysis of such valuable datasets requires expert understanding and, thus, complicates its broad application within the field. In the present study, we developed a user-friendly, and high-throughput web application named “Snake Toxin and Antivenom Binding Profiles” (STAB Profiles), to allow straight-forward analysis of hdpm datasets. To test our tool and evaluate its performance with a large dataset, we conducted hdpm assays using all African snake toxin protein sequences available in the UniProt database at the time of study design, together with eight commercial antivenoms in clinical use in Africa, thus representing the largest venom-antivenom dataset to date. Furthermore, we introduced a novel method for evaluating raw signals from a peptide microarray experiment and a data normalization protocol enabling intra-microarray and even inter-microarray chip comparisons. Finally, these data, alongside all the data from previous similar studies by Engmark et al., were preprocessed according to our newly developed protocol and made publicly available for download through the STAB Profiles web application (http://tropicalpharmacology.com/tools/stab-profiles/). With these data and our tool, we were able to gain key insights into toxin-antivenom interactions and were able to differentiate the ability of different antivenoms to interact with certain toxins of interest. The data, as well as the web application, we present in this article should be of significant value to the venom-antivenom research community. Knowledge gained from our current and future analyses of this dataset carry the potential to guide the improvement and optimization of current antivenoms for maximum patient benefit, as well as aid the development of next-generation antivenoms.UCR::Vicerrectoría de Investigación::Unidades de Investigación::Ciencias de la Salud::Instituto Clodomiro Picado (ICP

Convergent evolution of toxin resistance in animals

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

Venom complexity in a pitviper produced by facultative parthenogenesis

Facultative parthenogenesis (FP) is asexual reproduction in plant and animal species that would otherwise reproduce sexually. This process in vertebrates typically results from automictic development (likely terminal fusion) and is phylogenetically widespread. In squamate reptiles and chondrichthyan fishes, FP has been reported to occur in nature and can result in the production of reproductively viable offspring; suggesting that it is of ecological and evolutionary significance. However, terminal fusion automixis is believed to result in near genome-wide reductions in heterozygosity; thus, FP seems likely to affect key phenotypic characters, yet this remains almost completely unstudied. Snake venom is a complex phenotypic character primarily used to subjugate prey and is thus tightly linked to individual fitness. Surprisingly, the composition and function of venom produced by a parthenogenetic pitviper exhibits a high degree of similarity to that of its mother and conspecifics from the same population. Therefore, the apparent loss of allelic diversity caused by FP appears unlikely to have a significant impact on the prey-capturing ability of this snake. Accordingly, the pitviper offspring produced by FP retained complex phenotypic characteristics associated with fitness. This result reinforces the potential ecological and evolutionary importance of FP and questions our understanding of the inheritance of venom-associated genes

Physiological constraints dictate toxin spatial heterogeneity in snake venom glands

Background: Venoms are ecological innovations that have evolved numerous times, on each occasion accompanied by the co-evolution of specialised morphological and behavioural characters for venom production and delivery. The close evolutionary interdependence between these characters is exemplified by animals that control the composition of their secreted venom. This ability depends in part on the production of different toxins in different locations of the venom gland, which was recently documented in venomous snakes. Here, we test the hypothesis that the distinct spatial distributions of toxins in snake venom glands are an adaptation that enables the secretion of venoms with distinct ecological functions. Results: We show that the main defensive and predatory peptide toxins are produced in distinct regions of the venom glands of the black-necked spitting cobra (Naja nigricollis), but these distributions likely reflect developmental effects. Indeed, we detected no significant differences in venom collected via defensive ‘spitting’ or predatory ‘biting’ events from the same specimens representing multiple lineages of spitting cobra. We also found the same spatial distribution of toxins in a non-spitting cobra and show that heterogeneous toxin distribution is a feature shared with a viper with primarily predatory venom. Conclusions: Our findings suggest that heterogeneous distributions of toxins are not an adaptation to controlling venom composition in snakes. Instead, it likely reflects physiological constraints on toxin production by the venom glands, opening avenues for future research on the mechanisms of functional differentiation of populations of protein-secreting cells within adaptive contexts

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

Visual pigments, ocular filters and the evolution of snake vision

Much of what is known about the molecular evolution of vertebrate vision comes from studies of mammals, birds and fish. Reptiles (especially snakes) have barely been sampled in previous studies despite their exceptional diversity of retinal photoreceptor complements. Here we analyse opsin gene sequences and ocular media transmission for up to 69 species to investigate snake visual evolution. Most snakes express three visual opsin genes (rh1, sws1, lws). These opsin genes (especially rh1 and sws1) have undergone much evolutionary change, including modifications of amino acid residues at sites of known importance for spectral tuning, with several tuning site combinations unknown elsewhere among vertebrates. These changes are particularly common among dipsadine and colubrine ‘higher’ snakes. All three opsin genes are inferred to be under purifying selection, though dN/dS varies with respect to some lineages, ecologies, and retinal anatomy. Positive selection was inferred at multiple sites in all three opsins, these being concentrated in transmembrane domains and thus likely to have a substantial effect on spectral tuning and other aspects of opsin function. Snake lenses vary substantially in their spectral transmission. Snakes active at night and some of those active by day have very transmissive lenses, while some primarily diurnal species cut out shorter wavelengths (including UVA). In terms of retinal anatomy, lens transmission, visual pigment spectral tuning and opsin gene evolution the visual system of snakes is exceptionally diverse compared to all other extant tetrapod orders

A long-term observational study of paediatric snakebite in Kilifi County, south-east Kenya

INTRODUCTION: Estimates suggest that one-third of snakebite cases in sub-Saharan Africa affect children. Despite children being at a greater risk of disability and death, there are limited published data. This study has determined the: population-incidence and mortality rate of hospital-attended paediatric snakebite; clinical syndromes of snakebite envenoming; and predictors of severe local tissue damage. METHODS: All children presenting to Kilifi County Hospital, Kenya with snakebite were identified through the Kilifi Health and Demographic Surveillance System (KHDSS). Cases were prospectively registered, admitted for at least 24-hours, and managed on a paediatric high dependency unit (HDU). Households within the KHDSS study area have been included in 4-monthly surveillance and verbal autopsy, enabling calculation of population-incidence and mortality. Predictors of severe local tissue damage were identified using a multivariate logistic regression analysis. RESULTS: Between 2003 and 2021, there were 19,606 admissions to the paediatric HDU, of which 584 were due to snakebite. Amongst young children (≤5-years age) the population-incidence of hospital-attended snakebite was 11.3/100,000 person-years; for children aged 6-12 years this was 29.1/100,000 person-years. Incidence remained consistent over the study period despite the population size increasing (98,967 person-years in 2006; and 153,453 person-years in 2021). Most cases had local envenoming alone, but there were five snakebite associated deaths. Low haemoglobin; raised white blood cell count; low serum sodium; high systolic blood pressure; and an upper limb bite-site were independently associated with the development of severe local tissue damage. CONCLUSION: There is a substantial burden of disease due to paediatric snakebite, and the annual number of cases has increased in-line with population growth. The mortality rate was low, which may reflect the species causing snakebite in this region. The identification of independent predictors of severe local tissue damage can help to inform future research to better understand the pathophysiology of this important complication